Feed-Forward Propagation of Temporal and Rate Information between Cortical Populations during Coherent Activation in Engineered In Vitro Networks.

Transient propagation of information across neuronal assembles is thought to underlie many cognitive processes. However, the nature of the neural code that is embedded within these transmissions remains uncertain. Much of our understanding of how information is transmitted among these assemblies has been derived from computational models. While these models have been instrumental in understanding these processes they often make simplifying assumptions about the biophysical properties of neurons that may influence the nature and properties expressed. To address this issue we created an in vitro analog of a feed-forward network composed of two small populations (also referred to as assemblies or layers) of living dissociated rat cortical neurons. The populations were separated by, and communicated through, a microelectromechanical systems (MEMS) device containing a strip of microscale tunnels. Delayed culturing of one population in the first layer followed by the second a few days later induced the unidirectional growth of axons through the microtunnels resulting in a primarily feed-forward communication between these two small neural populations. In this study we systematically manipulated the number of tunnels that connected each layer and hence, the number of axons providing communication between those populations. We then assess the effect of reducing the number of tunnels has upon the properties of between-layer communication capacity and fidelity of neural transmission among spike trains transmitted across and within layers. We show evidence based on Victor-Purpura's and van Rossum's spike train similarity metrics supporting the presence of both rate and temporal information embedded within these transmissions whose fidelity increased during communication both between and within layers when the number of tunnels are increased. We also provide evidence reinforcing the role of synchronized activity upon transmission fidelity during the spontaneous synchronized network burst events that propagated between layers and highlight the potential applications of these MEMs devices as a tool for further investigation of structure and functional dynamics among neural populations

Neuronal dynamics and connectivity analysis of neuronal cultures on multi electrode arrays

Despite a number of attempts over the past two decades, research into reliable, controlled induction of long term evoked responses, mimicking low level learning and memory in dissociated cell cultures remains challenging. In addition, a full understanding of the stimulus-response relationships that underlie synaptic plasticity has not yet been achieved, and many of the underlying principles remain largely unknown. Plasticity studies have been predominantly limited to low density Multi/Micro Electrode Arrays (MEAs). With the advent of complementary metal-oxide-semiconductor (CMOS) based High-Density (HD) MEAs, unprecedented spatial and temporal resolution is now possible. In this thesis, an attempt to bridge the gap between studies of neural plasticity and the use of CMOS based HD-MEAs with thousands of electrodes, is reported. Additionally, since such HD-MEAs generate a large volume of data and require advanced analytics to efficiently process and analyse recordings, computational tools and novel algorithms to infer connectivity during plasticity have been developed. The study showed that the responsiveness, stability and initial firing rate of neuronal cultures are the deciding factors to reliably induce evoked responses. With multi-site stimulation, sustained long term potentiation was achieved, which was validated both by evoked response plots and overall firing rates measured at five different time points - before and after repeated stimulation, and at a three day time points. In contrast, while depression responses were observed, it was found that the effects were not sustained over many days. The findings of the study suggest that appropriate selection of neuronal cultures is crucial for inducing desired evoked responses and criteria for this have been developed. Furthermore, it is concluded that the initial responses to test stimuli can be used to determine whether potentiated or depressed responses are to be expected. To analyse the recordings, pipeline of computational tools was developed. Firstly, neuronal synchrony metrics were adapted for the first time for large HD-MEA recordings and shown to correspond effectively to the firing dynamics. To analyse functional connectivity, an information theoretic approach, Transfer Entropy(TE), was utilised. The method showed accurate estimation of functional connectivity with mid 80th percentile accuracy on simulated data. A superimposition method was proposed to enhance confidence in the connectivity estimation. To statistically evaluate connectivity estimation, a new surrogate method, based on ISI distribution approach, was proposed and validated with a simulated Izhikevich network. The method achieved improved accuracy, compared to the existing ISI shuffling method. This newly developed method was later utilised to infer connectivity and refine connections during the learning process of real neuronal cultures over many days of stimulation. The connectivity inference corresponded accurately to both the spontaneous and stimulated networks during evoked responses and the proposed method permitted observation of the evolution of connections for the potentiated network

Spatio-temporal dependencies in functional connectivity in rodent cortical cultures

Models of functional connectivity in cortical cultures on multi-electrodes arrays may aid in understanding how cognitive pathways form and improve techniques that aim to interface with neuronal systems. To enable research on such models, this study uses both data- and model-driven approaches to determine what dependencies are present in and between functional connectivity networks derived from bursts of extracellularly recorded activity. Properties of excitation in bursts were analysed using correlative techniques to assess the degree of linear dependence and then two parallel techniques were used to assess functional connectivity. Three models presenting increasing levels of spatio-temporal dependency were used to capture the dynamics of individual functional connections and their consistencies were verified using surrogate data. By comparing network-wide properties between model generated networks and functional networks from data, complex interdependencies were revealed. This indicates the persistent co-activation of neuronal pathways in spontaneous bursts, as can be found in whole brain structures

Fractals in the Nervous System: conceptual Implications for Theoretical Neuroscience

This essay is presented with two principal objectives in mind: first, to document the prevalence of fractals at all levels of the nervous system, giving credence to the notion of their functional relevance; and second, to draw attention to the as yet still unresolved issues of the detailed relationships among power law scaling, self-similarity, and self-organized criticality. As regards criticality, I will document that it has become a pivotal reference point in Neurodynamics. Furthermore, I will emphasize the not yet fully appreciated significance of allometric control processes. For dynamic fractals, I will assemble reasons for attributing to them the capacity to adapt task execution to contextual changes across a range of scales. The final Section consists of general reflections on the implications of the reviewed data, and identifies what appear to be issues of fundamental importance for future research in the rapidly evolving topic of this review

Neural Avalanches at the Critical Point between Replay and Non-Replay of Spatiotemporal Patterns

We model spontaneous cortical activity with a network of coupled spiking units, in which multiple spatio-temporal patterns are stored as dynamical attractors. We introduce an order parameter, which measures the overlap (similarity) between the activity of the network and the stored patterns. We find that, depending on the excitability of the network, different working regimes are possible. For high excitability, the dynamical attractors are stable, and a collective activity that replays one of the stored patterns emerges spontaneously, while for low excitability, no replay is induced. Between these two regimes, there is a critical region in which the dynamical attractors are unstable, and intermittent short replays are induced by noise. At the critical spiking threshold, the order parameter goes from zero to one, and its fluctuations are maximized, as expected for a phase transition (and as observed in recent experimental results in the brain). Notably, in this critical region, the avalanche size and duration distributions follow power laws. Critical exponents are consistent with a scaling relationship observed recently in neural avalanches measurements. In conclusion, our simple model suggests that avalanche power laws in cortical spontaneous activity may be the effect of a network at the critical point between the replay and non-replay of spatio-temporal patterns

Experimental analysis and computational modeling of interburst intervals in spontaneous activity of cortical neuronal culture

Rhythmic bursting is the most striking behavior of cultured cortical networks and may start in the second week after plating. In this study, we focus on the intervals between spontaneously occurring bursts, and compare experimentally recorded values with model simulations. In the models, we use standard neurons and synapses, with physiologically plausible parameters taken from literature. All networks had a random recurrent architecture with sparsely connected neurons. The number of neurons varied between 500 and 5,000. We find that network models with homogeneous synaptic strengths produce asynchronous spiking or stable regular bursts. The latter, however, are in a range not seen in recordings. By increasing the synaptic strength in a (randomly chosen) subset of neurons, our simulations show interburst intervals (IBIs) that agree better with in vitro experiments. In this regime, called weakly synchronized, the models produce irregular network bursts, which are initiated by neurons with relatively stronger synapses. In some noise-driven networks, a subthreshold, deterministic, input is applied to neurons with strong synapses, to mimic pacemaker network drive. We show that models with such “intrinsically active neurons” (pacemaker-driven models) tend to generate IBIs that are determined by the frequency of the fastest pacemaker and do not resemble experimental data. Alternatively, noise-driven models yield realistic IBIs. Generally, we found that large-scale noise-driven neuronal network models required synaptic strengths with a bimodal distribution to reproduce the experimentally observed IBI range. Our results imply that the results obtained from small network models cannot simply be extrapolated to models of more realistic size. Synaptic strengths in large-scale neuronal network simulations need readjustment to a bimodal distribution, whereas small networks do not require such change

CELL-BASED SENSOR CHIP FOR NEUROTOXICITY MEASUREMENTS IN DRINKING WATER

Our drinking water contains residues of pharmaceuticals. A sub-group of these contaminants are neuro-active substances, the antiepileptic carbamazepine being one of the most relevant. For assessment of the neurotoxicity of this drug at a sub-therapeutic level, a cell-based sensor chip platform has been realized and characterized. For this purpose, a microelectrode array chip was designed and processed in a clean room and optimized in terms of low processing costs and good recording properties. For characterization of the system neuronal cells were plated on microelectrode array chips and electrical activity was measured as a function of applied carbamazepine concentration. We found that the relative spike rate decreased with increasing drug concentration resulted in IC50 values of around 36 μM. This value is five orders of magnitude higher than the maximal dose found in drinking water. IC50 values for burst rate, burst duration and synchrony were slightly higher, suggesting spike rate being a more sensitive parameter to carbamazepine