Slow oscillating population activity in developing cortical networks: models and experimental results

Network Bursts in 3D Neuron Clusters Cultured on Microcontact-Printed Substrates

Microcontact printing (CP) is widely used to guide neurons to form 2D networks for neuroscience research. However, it is still difficult to establish 3D neuronal cultures on the CP substrate even though 3D neuronal structures are able to recapitulate critical aspects of native tissue. Here, we demonstrate that the reduced cell-substrate adhesion caused by the CP substrate could conveniently facilitate the aggregate formation of large-scale 3D neuron cluster networks. Furthermore, based on the quantitative analysis of the calcium activity of the resulting cluster networks, the effect of cell seeding density and local restriction of the CP substrate on network dynamics was investigated in detail. The results revealed that cell aggregation degree, rather than cell number, could take on the main role of the generation of synchronized network-wide calcium oscillation (network bursts) in the 3D neuron cluster networks. This finding may provide new insights for easy and cell-saving construction of in vitro 3D pathological models of epilepsy, and into deciphering the onset and evolution of network bursts in developmental nerve systems.

The Structural E/I Balance Constrains the Early Development of Cortical Network Activity

Neocortical networks have a characteristic constant ratio in the number of glutamatergic projection neurons (PN) and GABAergic interneurons (IN), and deviations in this ratio are often associated with developmental neuropathologies. Cultured networks with defined cellular content allowed us to ask if initial PN/IN ratios change the developmental population dynamics, and how different ratios impact the physiological excitatory/inhibitory (E/I) balance and the network activity development. During the first week in vitro, the IN content modulated PN numbers, increasing their proliferation in networks with higher IN proportions. The proportion of INs in each network set remained similar to the initial plating ratio during the 4 weeks cultivation period. Results from additional networks generated with more diverse cellular composition, including early-born GABA neurons, suggest that a GABA-dependent mechanism may decrease the survival of additional INs. A large variation of the PN/IN ratio did not change the balance between isolated spontaneous glutamatergic and GABAergic postsynaptic currents charge transfer (E/I balance) measured in PNs or INs. In contrast, the E/I balance of multisynaptic bursts reflected differences in IN content. Additionally, the spontaneous activity recorded by calcium imaging showed that higher IN ratios were associated with increased frequency of network bursts combined with a decrease of participating neurons per event. In the 4th week in vitro, bursting activity was stereotypically synchronized in networks with very few INs but was more desynchronized in networks with higher IN proportions. These results suggest that the E/I balance of isolated postsynaptic currents in single cells may be regulated independently of PN/IN proportions, but the network bursts E/I balance and the maturation of spontaneous network activity critically depends upon the structural PN/IN ratio.

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.

Spontaneous activity emerging from an inferred network model captures complex spatio-temporal dynamics of spike data

Inference methods are widely used to recover effective models from observed data. However, few studies attempted to investigate the dynamics of inferred models in neuroscience, and none, to our knowledge, at the network level. We introduce a principled modification of a widely used generalized linear model (GLM), and learn its structural and dynamic parameters from in-vitro spike data. The spontaneous activity of the new model captures prominent features of the non-stationary and non-linear dynamics displayed by the biological network, where the reference GLM largely fails, and also reflects fine-grained spatio-temporal dynamical features. Two ingredients were key for success. The first is a saturating transfer function: beyond its biological plausibility, it limits the neuron's information transfer, improving robustness against endogenous and external noise. The second is a super-Poisson spikes generative mechanism; it accounts for the undersampling of the network, and allows the model neuron to flexibly incorporate the observed activity fluctuations.

Challenges in Reproducibility, Replicability, and Comparability of Computational Models and Tools for Neuronal and Glial Networks, Cells, and Subcellular Structures

The possibility to replicate and reproduce published research results is one of the biggest challenges in all areas of science. In computational neuroscience, there are thousands of models available. However, it is rarely possible to reimplement the models based on the information in the original publication, let alone rerun the models just because the model implementations have not been made publicly available. We evaluate and discuss the comparability of a versatile choice of simulation tools: tools for biochemical reactions and spiking neuronal networks, and relatively new tools for growth in cell cultures. The replicability and reproducibility issues are considered for computational models that are equally diverse, including the models for intracellular signal transduction of neurons and glial cells, in addition to single glial cells, neuron-glia interactions, and selected examples of spiking neuronal networks. We also address the comparability of the simulation results with one another to comprehend if the studied models can be used to answer similar research questions. In addition to presenting the challenges in reproducibility and replicability of published results in computational neuroscience, we highlight the need for developing recommendations and good practices for publishing simulation tools and computational models. Model validation and flexible model description must be an integral part of the tool used to simulate and develop computational models. Constant improvement on experimental techniques and recording protocols leads to increasing knowledge about the biophysical mechanisms in neural systems. This poses new challenges for computational neuroscience: extended or completely new computational methods and models may be required. Careful evaluation and categorization of the existing models and tools provide a foundation for these future needs, for constructing multiscale models or extending the models to incorporate additional or more detailed biophysical mechanisms. Improving the quality of publications in computational neuroscience, enabling progressive building of advanced computational models and tools, can be achieved only through adopting publishing standards which underline replicability and reproducibility of research results.

Assessment of the Effects of Endocrine Disrupting Compounds on the Development of Vertebrate Neural Network Function Using Multi-electrode Arrays

Bis-phenols, such as bis-phenol A (BPA) and bis-phenol-S (BPS), are polymerizing agents widely used in the production of plastics and numerous everyday products. They are classified as endocrine disrupting compounds (EDC) with estradiol-like properties. Long-term exposure to EDCs, even at low doses, has been linked with various health defects including cancer, behavioral disorders, and infertility, with greater vulnerability during early developmental periods. To study the effects of BPA on the development of neuronal function, we used an in vitro neuronal network derived from the early chick embryonic brain as a model. We found that exposure to BPA affected the development of network activity, specifically spiking activity and synchronization. A change in network activity is the crucial link between the molecular target of a drug or compound and its effect on behavioral outcome. Multi-electrode arrays are increasingly becoming useful tools to study the effects of drugs on network activity in vitro. There are several systems available in the market and, although there are variations in the number of electrodes, the type and quality of the electrode array and the analysis software, the basic underlying principles, and the data obtained is the same across the different systems. Although currently limited to analysis of two-dimensional in vitro cultures, these MEA systems are being improved to enable in vivo network activity in brain slices. Here, we provide a detailed protocol for embryonic exposure and recording neuronal network activity and synchrony, along with representative results.

Recurrently connected and localized neuronal communities initiate coordinated spontaneous activity in neuronal networks

Developing neuronal systems intrinsically generate coordinated spontaneous activity that propagates by involving a large number of synchronously firing neurons. In vivo, waves of spikes transiently characterize the activity of developing brain circuits and are fundamental for activity-dependent circuit formation. In vitro, coordinated spontaneous spiking activity, or network bursts (NBs), interleaved within periods of asynchronous spikes emerge during the development of 2D and 3D neuronal cultures. Several studies have investigated this type of activity and its dynamics, but how a neuronal system generates these coordinated events remains unclear. Here, we investigate at a cellular level the generation of network bursts in spontaneously active neuronal cultures by exploiting high-resolution multielectrode array recordings and computational network modelling. Our analysis reveals that NBs are generated in specialized regions of the network (functional neuronal communities) that feature neuronal links with high cross-correlation peak values, sub-millisecond lags and that share very similar structural connectivity motifs providing recurrent interactions. We show that the particular properties of these local structures enable locally amplifying spontaneous asynchronous spikes and that this mechanism can lead to the initiation of NBs. Through the analysis of simulated and experimental data, we also show that AMPA currents drive the coordinated activity, while NMDA and GABA currents are only involved in shaping the dynamics of NBs. Overall, our results suggest that the presence of functional neuronal communities with recurrent local connections allows a neuronal system to generate spontaneous coordinated spiking activity events. As suggested by the rules used for implementing our computational model, such functional communities might naturally emerge during network development by following simple constraints on distance-based connectivity.

Coordinated spontaneous spiking activity is fundamental for the normal formation of brain circuits during development. However, how ensembles of neurons generate these events remains unclear. To address this question, in the present study, we investigated the network properties that might be required to a neuronal system for the generation of these spontaneous waves of spikes. We performed our study on spontaneously active neuronal cell cultures using high-resolution electrical recordings and a computational network model developed to reproduce our experimental data both quantitatively and qualitatively. Through the analysis of both experimental and simulated data, we found that network bursts are initiated in regions of the network, or “functional communities”, characterized by particular local connectivity properties. We also found that these regions can amplify the background asynchronous spiking activity preceding a network burst and, in this way, can give rise to coordinated spiking events.

As a whole, our results suggest the presence of functional communities of neurons in a developing neuronal system that might naturally emerge by following simple constraints on distance-based connectivity. These regions are most likely required for the generation of the spontaneous coordinated activity that can drive activity-dependent circuit formation.

Size-dependent regulation of synchronized activity in living neuronal networks

We study the effect of network size on synchronized activity in living neuronal networks. Dissociated cortical neurons form synaptic connections in culture and generate synchronized spontaneous activity within 10 days in vitro. Using micropatterned surfaces to extrinsically control the size of neuronal networks, we show that synchronized activity can emerge in a network as small as 12 cells. Furthermore, a detailed comparison of small (∼20 cells), medium (∼100 cells), and large (∼400 cells) networks reveal that synchronized activity becomes destabilized in the small networks. A computational modeling of neural activity is then employed to explore the underlying mechanism responsible for the size effect. We find that the generation and maintenance of the synchronized activity can be minimally described by: (1) the stochastic firing of each neuron in the network, (2) enhancement in the network activity in a positive feedback loop of excitatory synapses, and (3) Ca-dependent suppression of bursting activity. The model further shows that the decrease in total synaptic input to a neuron that drives the positive feedback amplification of correlated activity is a key factor underlying the destabilization of synchrony in smaller networks. Spontaneous neural activity plays a critical role in cortical information processing, and our work constructively clarifies an aspect of the structural basis behind this.

Network Events on Multiple Space and Time Scales in Cultured Neural Networks and in a Stochastic Rate Model

Cortical networks, in-vitro as well as in-vivo, can spontaneously generate a variety of collective dynamical events such as network spikes, UP and DOWN states, global oscillations, and avalanches. Though each of them has been variously recognized in previous works as expression of the excitability of the cortical tissue and the associated nonlinear dynamics, a unified picture of the determinant factors (dynamical and architectural) is desirable and not yet available. Progress has also been partially hindered by the use of a variety of statistical measures to define the network events of interest. We propose here a common probabilistic definition of network events that, applied to the firing activity of cultured neural networks, highlights the co-occurrence of network spikes, power-law distributed avalanches, and exponentially distributed 'quasi-orbits', which offer a third type of collective behavior. A rate model, including synaptic excitation and inhibition with no imposed topology, synaptic short-term depression, and finite-size noise, accounts for all these different, coexisting phenomena. We find that their emergence is largely regulated by the proximity to an oscillatory instability of the dynamics, where the non-linear excitable behavior leads to a self-amplification of activity fluctuations over a wide range of scales in space and time. In this sense, the cultured network dynamics is compatible with an excitation-inhibition balance corresponding to a slightly sub-critical regime. Finally, we propose and test a method to infer the characteristic time of the fatigue process, from the observed time course of the network's firing rate. Unlike the model, possessing a single fatigue mechanism, the cultured network appears to show multiple time scales, signalling the possible coexistence of different fatigue mechanisms.

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.

Relationship between individual neuron and network spontaneous activity in developing mouse cortex

Spontaneous synchronous activity (SSA) that propagates as electrical waves is found in numerous central nervous system structures and is critical for normal development, but the mechanisms of generation of such activity are not clear. In previous work, we showed that the ventrolateral piriform cortex is uniquely able to initiate SSA in contrast to the dorsal neocortex, which participates in, but does not initiate, SSA (Lischalk JW, Easton CR, Moody WJ. Dev Neurobiol 69: 407-414, 2009). In this study, we used Ca(2+) imaging of cultured embryonic day 18 to postnatal day 2 coronal slices (embryonic day 17 + 1-4 days in culture) of the mouse cortex to investigate the different activity patterns of individual neurons in these regions. In the piriform cortex where SSA is initiated, a higher proportion of neurons was active asynchronously between waves, and a larger number of groups of coactive cells was present compared with the dorsal cortex. When we applied GABA and glutamate synaptic antagonists, asynchronous activity and cellular clusters remained, while synchronous activity was eliminated, indicating that asynchronous activity is a result of cell-intrinsic properties that differ between these regions. To test the hypothesis that higher levels of cell-autonomous activity in the piriform cortex underlie its ability to initiate waves, we constructed a conductance-based network model in which three layers differed only in the proportion of neurons able to intrinsically generate bursting behavior. Simulations using this model demonstrated that a gradient of intrinsic excitability was sufficient to produce directionally propagating waves that replicated key experimental features, indicating that the higher level of cell-intrinsic activity in the piriform cortex may provide a substrate for SSA generation.

Network bursting dynamics in excitatory cortical neuron cultures results from the combination of different adaptive mechanisms

In the brain, synchronization among cells of an assembly is a common phenomenon, and thought to be functionally relevant. Here we used an in vitro experimental model of cell assemblies, cortical cultures, combined with numerical simulations of a spiking neural network (SNN) to investigate how and why spontaneous synchronization occurs. In order to deal with excitation only, we pharmacologically blocked GABA_Aergic transmission using bicuculline. Synchronous events in cortical cultures tend to involve almost every cell and to display relatively constant durations. We have thus named these “network spikes” (NS). The inter-NS-intervals (INSIs) proved to be a more interesting phenomenon. In most cortical cultures NSs typically come in series or bursts (“bursts of NSs”, BNS), with short (∼1 s) INSIs and separated by long silent intervals (tens of s), which leads to bimodal INSI distributions. This suggests that a facilitating mechanism is at work, presumably short-term synaptic facilitation, as well as two fatigue mechanisms: one with a short timescale, presumably short-term synaptic depression, and another one with a longer timescale, presumably cellular adaptation. We thus incorporated these three mechanisms into the SNN, which, indeed, produced realistic BNSs. Next, we systematically varied the recurrent excitation for various adaptation timescales. Strong excitability led to frequent, quasi-periodic BNSs (CV∼0), and weak excitability led to rare BNSs, approaching a Poisson process (CV∼1). Experimental cultures appear to operate within an intermediate weakly-synchronized regime (CV∼0.5), with an adaptation timescale in the 2–8 s range, and well described by a Poisson-with-refractory-period model. Taken together, our results demonstrate that the INSI statistics are indeed informative: they allowed us to infer the mechanisms at work, and many parameters that we cannot access experimentally.

Thyroid hormone-dependent development of early cortical networks: temporal specificity and the contribution of trkB and mTOR pathways

Early in neocortical network development, triiodothyronine (T3) promotes GABAergic neurons' population increase, their somatic growth and the formation of GABAergic synapses. In the presence of T3, GABAergic interneurons form longer axons and conspicuous axonal arborizations, with an increased number of putative synaptic boutons. Here we show that the increased GABAergic axonal growth is positively correlated with the proximity to non-GABAergic neurons (non-GABA). A differential innervation emerges from a T3-dependent decrease of axonal length in fields with low density of neuronal cell bodies, combined with an increased bouton formation in fields with high density of neuronal somata. T3 addition to deprived networks after the first 2 weeks of development did not rescue deficits in the GABAergic synaptic bouton distribution, or in the frequency and duration of spontaneous bursts. During the critical 2-week-period, GABAergic signaling is depolarizing as revealed by calcium imaging experiments. Interestingly, T3 enhanced the expression of the potassium-chloride cotransporter 2 (KCC2), and accelerated the developmental shift from depolarizing to hyperpolarizing GABAergic signaling in non-GABA. The T3-related increase of spontaneous network activity was remarkably reduced after blockade of either tropomyosin-receptor kinase B (trkB) or mammalian target of rapamycin (mTOR) pathways. T3-dependent increase in GABAergic neurons' soma size was mediated mainly by mTOR signaling. Conversely, the T3-dependent selective increase of GABAergic boutons near non-GABAergic cell bodies is mediated by trkB signaling only. Both trkB and mTOR signaling mediate T3-dependent reduction of the GABAergic axon extension. The circuitry context is relevant for the interaction between T3 and trkB signaling, but not for the interactions between T3 and mTOR signaling.

Structure-dynamics relationships in bursting neuronal networks revealed using a prediction framework

The question of how the structure of a neuronal network affects its functionality has gained a lot of attention in neuroscience. However, the vast majority of the studies on structure-dynamics relationships consider few types of network structures and assess limited numbers of structural measures. In this in silico study, we employ a wide diversity of network topologies and search among many possibilities the aspects of structure that have the greatest effect on the network excitability. The network activity is simulated using two point-neuron models, where the neurons are activated by noisy fluctuation of the membrane potential and their connections are described by chemical synapse models, and statistics on the number and quality of the emergent network bursts are collected for each network type. We apply a prediction framework to the obtained data in order to find out the most relevant aspects of network structure. In this framework, predictors that use different sets of graph-theoretic measures are trained to estimate the activity properties, such as burst count or burst length, of the networks. The performances of these predictors are compared with each other. We show that the best performance in prediction of activity properties for networks with sharp in-degree distribution is obtained when the prediction is based on clustering coefficient. By contrast, for networks with broad in-degree distribution, the maximum eigenvalue of the connectivity graph gives the most accurate prediction. The results shown for small () networks hold with few exceptions when different neuron models, different choices of neuron population and different average degrees are applied. We confirm our conclusions using larger () networks as well. Our findings reveal the relevance of different aspects of network structure from the viewpoint of network excitability, and our integrative method could serve as a general framework for structure-dynamics studies in biosciences.