Development of input connections in neural cultures

Cluster formation due to repulsive spanning trees in attractively coupled networks

Ensembles of coupled nonlinear oscillators are a popular paradigm and an ideal benchmark for analyzing complex collective behaviors. The onset of cluster synchronization is found to be at the core of various technological and biological processes. The current literature has investigated cluster synchronization by focusing mostly on the case of attractive coupling among the oscillators. However, the case of two coexisting competing interactions is of practical interest due to their relevance in diverse natural settings, including neuronal networks consisting of excitatory and inhibitory neurons, the coevolving social model with voters of opposite opinions, and ecological plant communities with both facilitation and competition, to name a few. In the present article, we investigate the impact of repulsive spanning trees on cluster formation within a connected network of attractively coupled limit-cycle oscillators. We successfully predict which nodes belong to each cluster and the emergent frustration of the connected networks independent of the particular local dynamics at the network nodes. We also determine local asymptotic stability of the cluster states using an approach based on the formulation of a master stability function. We additionally validate the emergence of solitary states and antisynchronization for some specific choices of spanning trees and networks.

Preparation and Mechano-Functional Characterization of PEGylated Fibrin Hydrogels: Impact of Thrombin Concentration

Three-dimensional (3D) neuronal cultures grown in hydrogels are promising platforms to design brain-like neuronal networks in vitro. However, the optimal properties of such cultures must be tuned to ensure a hydrogel matrix sufficiently porous to promote healthy development but also sufficiently rigid for structural support. Such an optimization is difficult since it implies the exploration of different hydrogel compositions and, at the same time, a functional analysis to validate neuronal culture viability. To advance in this quest, here we present a combination of a rheological protocol and a network-based functional analysis to investigate PEGylated fibrin hydrogel networks with gradually higher stiffness, achieved by increasing the concentration of thrombin. We observed that moderate thrombin concentrations of 10% and 25% in volume shaped healthy networks, although the functional traits depended on the hydrogel stiffness, which was much higher for the latter concentration. Thrombin concentrations of 65% or higher led to networks that did not survive. Our results illustrate the difficulties and limitations in preparing 3D neuronal networks, and stress the importance of combining a mechano-structural characterization of a biomaterial with a functional one.

Time-dependent homeostatic mechanisms underlie brain-derived neurotrophic factor action on neural circuitry

Plasticity and homeostatic mechanisms allow neural networks to maintain proper function while responding to physiological challenges. Despite previous work investigating morphological and synaptic effects of brain-derived neurotrophic factor (BDNF), the most prevalent growth factor in the central nervous system, how exposure to BDNF manifests at the network level remains unknown. Here we report that BDNF treatment affects rodent hippocampal network dynamics during development and recovery from glutamate-induced excitotoxicity in culture. Importantly, these effects are not obvious when traditional activity metrics are used, so we delve more deeply into network organization, functional analyses, and in silico simulations. We demonstrate that BDNF partially restores homeostasis by promoting recovery of weak and medium connections after injury. Imaging and computational analyses suggest these effects are caused by changes to inhibitory neurons and connections. From our in silico simulations, we find that BDNF remodels the network by indirectly strengthening weak excitatory synapses after injury. Ultimately, our findings may explain the difficulties encountered in preclinical and clinical trials with BDNF and also offer information for future trials to consider.

Investigating the reliability of the evoked response in human iPSCs-derived neuronal networks coupled to micro-electrode arrays

In vitro models of neuronal networks have emerged as a potent instrument for gaining deeper insights into the intricate mechanisms governing the human brain. Notably, the integration of human-induced pluripotent stem cells (hiPSCs) with micro-electrode arrays offers a means to replicate and dissect both the structural and functional elements of the human brain within a controlled in vitro environment. Given that neuronal communication relies on the emission of electrical (and chemical) stimuli, the employment of electrical stimulation stands as a mean to comprehensively interrogate neuronal assemblies, to better understand their inherent electrophysiological dynamics. However, the establishment of standardized stimulation protocols for cultures derived from hiPSCs is still lacking, thereby hindering the precise delineation of efficacious parameters to elicit responses. To fill this gap, the primary objective of this study resides in delineating effective parameters for the electrical stimulation of hiPSCs-derived neuronal networks, encompassing the determination of voltage amplitude and stimulation frequency able to evoke reliable and stable responses. This study represents a stepping-stone in the exploration of efficacious stimulation parameters, thus broadening the electrophysiological activity profiling of neural networks sourced from human-induced pluripotent stem cells.

Modular architecture facilitates noise-driven control of synchrony in neuronal networks

High-level information processing in the mammalian cortex requires both segregated processing in specialized circuits and integration across multiple circuits. One possible way to implement these seemingly opposing demands is by flexibly switching between states with different levels of synchrony. However, the mechanisms behind the control of complex synchronization patterns in neuronal networks remain elusive. Here, we use precision neuroengineering to manipulate and stimulate networks of cortical neurons in vitro, in combination with an in silico model of spiking neurons and a mesoscopic model of stochastically coupled modules to show that (i) a modular architecture enhances the sensitivity of the network to noise delivered as external asynchronous stimulation and that (ii) the persistent depletion of synaptic resources in stimulated neurons is the underlying mechanism for this effect. Together, our results demonstrate that the inherent dynamical state in structured networks of excitable units is determined by both its modular architecture and the properties of the external inputs.

Early path dominance as a principle for neurodevelopment

We perform targeted attack, a systematic computational unlinking of the network, to analyze its effects on global communication across the brain network through its giant cluster. Across diffusion magnetic resonance images from individuals in the UK Biobank, Adolescent Brain Cognitive Development Study and Developing Human Connectome Project, we find that targeted attack procedures on increasing white matter tract lengths and densities are remarkably invariant to aging and disease. Time-reversing the attack computation suggests a mechanism for how brains develop, for which we derive an analytical equation using percolation theory. Based on a close match between theory and experiment, our results demonstrate that tracts are limited to emanate from regions already in the giant cluster and tracts that appear earliest in neurodevelopment are those that become the longest and densest.

Neuronal Cultures: Exploring Biophysics, Complex Systems, and Medicine in a Dish

Neuronal cultures are one of the most important experimental models in modern interdisciplinary neuroscience, allowing to investigate in a control environment the emergence of complex behavior from an ensemble of interconnected neurons. Here, I review the research that we have conducted at the neurophysics laboratory at the University of Barcelona over the last 15 years, describing first the neuronal cultures that we prepare and the associated tools to acquire and analyze data, to next delve into the different research projects in which we actively participated to progress in the understanding of open questions, extend neuroscience research on new paradigms, and advance the treatment of neurological disorders. I finish the review by discussing the drawbacks and limitations of neuronal cultures, particularly in the context of brain-like models and biomedicine.

Interlayer antisynchronization in degree-biased duplex networks

With synchronization being one of nature's most ubiquitous collective behaviors, the field of network synchronization has experienced tremendous growth, leading to significant theoretical developments. However, most previous studies consider uniform connection weights and undirected networks with positive coupling. In the present article, we incorporate the asymmetry in a two-layer multiplex network by assigning the ratio of the adjacent nodes' degrees as the weights to the intralayer edges. Despite the presence of degree-biased weighting mechanism and attractive-repulsive coupling strengths, we are able to find the necessary conditions for intralayer synchronization and interlayer antisynchronization and test whether these two macroscopic states can withstand demultiplexing in a network. During the occurrence of these two states, we analytically calculate the oscillator's amplitude. In addition to deriving the local stability conditions for interlayer antisynchronization via the master stability function approach, we also construct a suitable Lyapunov function to determine a sufficient condition for global stability. We provide numerical evidence to show the necessity of negative interlayer coupling strength for the occurrence of antisynchronization, and such repulsive interlayer coupling coefficients cannot destroy intralayer synchronization.

Rich dynamics and functional organization on topographically designed neuronal networks in vitro

Neuronal cultures are a prominent experimental tool to understand complex functional organization in neuronal assemblies. However, neurons grown on flat surfaces exhibit a strongly coherent bursting behavior with limited functionality. To approach the functional richness of naturally formed neuronal circuits, here we studied neuronal networks grown on polydimethylsiloxane (PDMS) topographical patterns shaped as either parallel tracks or square valleys. We followed the evolution of spontaneous activity in these cultures along 20 days in vitro using fluorescence calcium imaging. The networks were characterized by rich spatiotemporal activity patterns that comprised from small regions of the culture to its whole extent. Effective connectivity analysis revealed the emergence of spatially compact functional modules that were associated with both the underpinned topographical features and predominant spatiotemporal activity fronts. Our results show the capacity of spatial constraints to mold activity and functional organization, bringing new opportunities to comprehend the structure-function relationship in living neuronal circuits.

Spike Timing-Dependent Plasticity with Enhanced Long-Term Depression Leads to an Increase of Statistical Complexity

Synaptic plasticity is characterized by remodeling of existing synapses caused by strengthening and/or weakening of connections. This is represented by long-term potentiation (LTP) and long-term depression (LTD). The occurrence of a presynaptic spike (or action potential) followed by a temporally nearby postsynaptic spike induces LTP; conversely, if the postsynaptic spike precedes the presynaptic spike, it induces LTD. This form of synaptic plasticity induction depends on the order and timing of the pre- and postsynaptic action potential, and has been termed spike time-dependent plasticity (STDP). After an epileptic seizure, LTD plays an important role as a depressor of synapses, which may lead to their complete disappearance together with that of their neighboring connections until days after the event. Added to the fact that after an epileptic seizure the network seeks to regulate the excess activity through two key mechanisms: depressed connections and neuronal death (eliminating excitatory neurons from the network), LTD becomes of great interest in our study. To investigate this phenomenon, we develop a biologically plausible model that privileges LTD at the triplet level while maintaining the pairwise structure in the STPD and study how network dynamics are affected as neuronal damage increases. We find that the statistical complexity is significantly higher for the network where LTD presented both types of interactions. While in the case where the STPD is defined with purely pairwise interactions an increase is observed as damage becomes higher for both Shannon Entropy and Fisher information.

Spontaneous Activity Predicts Survival of Developing Cortical Neurons

Spontaneous activity plays a crucial role in brain development by coordinating the integration of immature neurons into emerging cortical networks. High levels and complex patterns of spontaneous activity are generally associated with low rates of apoptosis in the cortex. However, whether spontaneous activity patterns directly encode for survival of individual cortical neurons during development remains an open question. Here, we longitudinally investigated spontaneous activity and apoptosis in developing cortical cultures, combining extracellular electrophysiology with calcium imaging. These experiments demonstrated that the early occurrence of calcium transients was strongly linked to neuronal survival. Silent neurons exhibited a higher probability of cell death, whereas high frequency spiking and burst behavior were almost exclusively detected in surviving neurons. In local neuronal clusters, activity of neighboring neurons exerted a pro-survival effect, whereas on the functional level, networks with a high modular topology were associated with lower cell death rates. Using machine learning algorithms, cell fate of individual neurons was predictable through the integration of spontaneous activity features. Our results indicate that high frequency spiking activity constrains apoptosis in single neurons through sustained calcium rises and thereby consolidates networks in which a high modular topology is reached during early development.

Self-organization of in vitro neuronal assemblies drives to complex network topology

Activity-dependent self-organization plays an important role in the formation of specific and stereotyped connectivity patterns in neural circuits. By combining neuronal cultures, and tools with approaches from network neuroscience and information theory, we can study how complex network topology emerges from local neuronal interactions. We constructed effective connectivity networks using a transfer entropy analysis of spike trains recorded from rat embryo dissociated hippocampal neuron cultures between 6 and 35 days in vitro to investigate how the topology evolves during maturation. The methodology for constructing the networks considered the synapse delay and addressed the influence of firing rate and population bursts as well as spurious effects on the inference of connections. We found that the number of links in the networks grew over the course of development, shifting from a segregated to a more integrated architecture. As part of this progression, three significant aspects of complex network topology emerged. In agreement with previous in silico and in vitro studies, a small-world architecture was detected, largely due to strong clustering among neurons. Additionally, the networks developed in a modular topology, with most modules comprising nearby neurons. Finally, highly active neurons acquired topological characteristics that made them important nodes to the network and integrators of modules. These findings leverage new insights into how neuronal effective network topology relates to neuronal assembly self-organization mechanisms.

Mechanistic insights into ultrasonic neurostimulation of disconnected neurons using single short pulses

Ultrasonic neurostimulation is a potentially potent noninvasive therapy, whose mechanism has yet to be elucidated. We designed a system capable of applying ultrasound with minimal reflections to neuronal cultures. Synaptic transmission was pharmacologically controlled, eliminating network effects, enabling examination of single-cell processes. Short single pulses of low-intensity ultrasound were applied, and time-locked responses were examined using calcium imaging. Low-pressure (0.35MPa) ultrasound directly stimulated ∼20% of pharmacologically disconnected neurons, regardless of membrane poration. Stimulation was resistant to the blockade of several purinergic receptor and mechanosensitive ion channel types. Stimulation was blocked, however, by suppression of action potentials. Surprisingly, even extremely short (4μs) pulses were effective, stimulating ∼8% of the neurons. Lower-pressure pulses (0.35MPa) were less effective than higher-pressure ones (0.65MPa). Attrition effects dominated, with no indication of compromised viability. Our results detract from theories implicating cavitation, heating, non-transient membrane pores >1.5nm, pre-synaptic release, or gradual effects. They implicate a post-synaptic mechanism upstream of the action potential, and narrow down the list of possible targets involved.

Explosive Synchronization-Based Brain Modulation Reduces Hypersensitivity in the Brain Network: A Computational Model Study

Fibromyalgia (FM) is a chronic pain condition that is characterized by hypersensitivity to multimodal sensory stimuli, widespread pain, and fatigue. We have previously proposed explosive synchronization (ES), a phenomenon wherein a small perturbation to a network can lead to an abrupt state transition, as a potential mechanism of the hypersensitive FM brain. Therefore, we hypothesized that converting a brain network from ES to general synchronization (GS) may reduce the hypersensitivity of FM brain. To find an effective brain network modulation to convert ES into GS, we constructed a large-scale brain network model near criticality (i.e., an optimally balanced state between order and disorders), which reflects brain dynamics in conscious wakefulness, and adjusted two parameters: local structural connectivity and signal randomness of target brain regions. The network sensitivity to global stimuli was compared between the brain networks before and after the modulation. We found that only increasing the local connectivity of hubs (nodes with intense connections) changes ES to GS, reducing the sensitivity, whereas other types of modulation such as decreasing local connectivity, increasing and decreasing signal randomness are not effective. This study would help to develop a network mechanism-based brain modulation method to reduce the hypersensitivity in FM.

Discrete field potentials produced by coherent activation of spinal dorsal horn neurons

In addition to the action potentials generated by the ongoing activation of single dorsal horn neurons in the anesthetized cat, we often recorded small negative field potentials with a fast-rising phase and a slow decay (dIFPs). These potentials could be separated in different classes, each with a specific and rather constant shape and amplitude. They were largest in spinal laminae III-V and gradually faded at deeper locations, without showing the polarity reversal displayed at these depths by the focal potentials produced by stimulation of muscle and cutaneous afferents. We propose that the dIFPs are postsynaptic field potentials generated by strongly coupled sets of dorsal horn neurons displaying a spatial orientation that generates closed field potentials in response to stimulation of high-threshold cutaneous and muscle afferents. These neuronal sets could form part of the spinal inhibitory circuitry that mediates presynaptic inhibition and Ib non-reciprocal postsynaptic inhibition and could be involved in the sensory-motor transformations activated by stimulation of high-threshold cutaneous afferents.

Percolation in networks with local homeostatic plasticity

Percolation is a process that impairs network connectedness by deactivating links or nodes. This process features a phase transition that resembles paradigmatic critical transitions in epidemic spreading, biological networks, traffic and transportation systems. Some biological systems, such as networks of neural cells, actively respond to percolation-like damage, which enables these structures to maintain their function after degradation and aging. Here we study percolation in networks that actively respond to link damage by adopting a mechanism resembling synaptic scaling in neurons. We explain critical transitions in such active networks and show that these structures are more resilient to damage as they are able to maintain a stronger connectedness and ability to spread information. Moreover, we uncover the role of local rescaling strategies in biological networks and indicate a possibility of designing smart infrastructures with improved robustness to perturbations.

Physical and chemical properties of carbon nanotubes in view of mechanistic neuroscience investigations. Some outlook from condensed matter, materials science and physical chemistry

The open border between non-living and living matter, suggested by increasingly emerging fields of nanoscience interfaced to biological systems, requires a detailed knowledge of nanomaterials properties. An account of the wide spectrum of phenomena, belonging to physical chemistry of interfaces, materials science, solid state physics at the nanoscale and bioelectrochemistry, thus is acquainted for a comprehensive application of carbon nanotubes interphased with neuron cells. This review points out a number of conceptual tools to further address the ongoing advances in coupling neuronal networks with (carbon) nanotube meshworks, and to deepen the basic issues that govern a biological cell or tissue interacting with a nanomaterial. Emphasis is given here to the properties and roles of carbon nanotube systems at relevant spatiotemporal scales of individual molecules, junctions and molecular layers, as well as to the point of view of a condensed matter or materials scientist. Carbon nanotube interactions with blood-brain barrier, drug delivery, biocompatibility and functionalization issues are also regarded.

Inferring Excitatory and Inhibitory Connections in Neuronal Networks

The comprehension of neuronal network functioning, from most basic mechanisms of signal transmission to complex patterns of memory and decision making, is at the basis of the modern research in experimental and computational neurophysiology. While mechanistic knowledge of neurons and synapses structure increased, the study of functional and effective networks is more complex, involving emergent phenomena, nonlinear responses, collective waves, correlation and causal interactions. Refined data analysis may help in inferring functional/effective interactions and connectivity from neuronal activity. The Transfer Entropy (TE) technique is, among other things, well suited to predict structural interactions between neurons, and to infer both effective and structural connectivity in small- and large-scale networks. To efficiently disentangle the excitatory and inhibitory neural activities, in the article we present a revised version of TE, split in two contributions and characterized by a suited delay time. The method is tested on in silico small neuronal networks, built to simulate the calcium activity as measured via calcium imaging in two-dimensional neuronal cultures. The inhibitory connections are well characterized, still preserving a high accuracy for excitatory connections prediction. The method could be applied to study effective and structural interactions in systems of excitable cells, both in physiological and in pathological conditions.

Action potential-coupled Rho GTPase signaling drives presynaptic plasticity

In contrast to their postsynaptic counterparts, the contributions of activity-dependent cytoskeletal signaling to presynaptic plasticity remain controversial and poorly understood. To identify and evaluate these signaling pathways, we conducted a proteomic analysis of the presynaptic cytomatrix using in vivo biotin identification (iBioID). The resultant proteome was heavily enriched for actin cytoskeleton regulators, including Rac1, a Rho GTPase that activates the Arp2/3 complex to nucleate branched actin filaments. Strikingly, we find Rac1 and Arp2/3 are closely associated with synaptic vesicle membranes in adult mice. Using three independent approaches to alter presynaptic Rac1 activity (genetic knockout, spatially restricted inhibition, and temporal optogenetic manipulation), we discover that this pathway negatively regulates synaptic vesicle replenishment at both excitatory and inhibitory synapses, bidirectionally sculpting short-term synaptic depression. Finally, we use two-photon fluorescence lifetime imaging to show that presynaptic Rac1 activation is coupled to action potentials by voltage-gated calcium influx. Thus, this study uncovers a previously unrecognized mechanism of actin-regulated short-term presynaptic plasticity that is conserved across excitatory and inhibitory terminals. It also provides a new proteomic framework for better understanding presynaptic physiology, along with a blueprint of experimental strategies to isolate the presynaptic effects of ubiquitously expressed proteins.

Parkinson's disease patient-specific neuronal networks carrying the LRRK2 G2019S mutation unveil early functional alterations that predate neurodegeneration

A deeper understanding of early disease mechanisms occurring in Parkinson’s disease (PD) is needed to reveal restorative targets. Here we report that human induced pluripotent stem cell (iPSC)-derived dopaminergic neurons (DAn) obtained from healthy individuals or patients harboring LRRK2 PD-causing mutation can create highly complex networks with evident signs of functional maturation over time. Compared to control neuronal networks, LRRK2 PD patients’ networks displayed an elevated bursting behavior, in the absence of neurodegeneration. By combining functional calcium imaging, biophysical modeling, and DAn-lineage tracing, we found a decrease in DAn neurite density that triggered overall functional alterations in PD neuronal networks. Our data implicate early dysfunction as a prime focus that may contribute to the initiation of downstream degenerative pathways preceding DAn loss in PD, highlighting a potential window of opportunity for pre-symptomatic assessment of chronic degenerative diseases.

Neuronal circuits on a chip for biological network monitoring

Cultured neuronal networks (CNNs) are a robust model to closely investigate neuronal circuits' formation and monitor their structural properties evolution. Typically, neurons are cultured in plastic plates or, more recently, in microfluidic platforms with potentially a wide variety of neuroscience applications. As a biological protocol, cell culture integration with a microfluidic system provides benefits such as accurate control of cell seeding area, culture medium renewal, or lower exposure to contamination. The objective of this report is to present a novel neuronal network on a chip device, including a chamber, fabricated from PDMS, vinyl and glass connected to a microfluidic platform to perfuse the continuous flow of culture medium. Network growth is compared in chips and traditional Petri dishes to validate the microfluidic chip performance. The network assessment is performed by computing relevant topological measures like the number of connected neurons, the clustering coefficient, and the shortest path between any pair of neurons throughout the culture's life. The results demonstrate that neuronal circuits on a chip have a more stable network structure and lifespan than developing in conventional settings, and therefore this setup is an advantageous alternative to current culture methods. This technology could lead to challenging applications such as batch drug testing of in vitro cell culture models. From the engineering perspective, a device's advantage is the chance to develop custom designs more efficiently than other microfluidic systems.

Neuronal circuits overcome imbalance in excitation and inhibition by adjusting connection numbers

The interplay between excitation and inhibition is crucial for neuronal circuitry in the brain. Inhibitory cell fractions in the neocortex and hippocampus are typically maintained at 15 to 30%, which is assumed to be important for stable dynamics. We have studied systematically the role of precisely controlled excitatory/inhibitory (E/I) cellular ratios on network activity using mice hippocampal cultures. Surprisingly, networks with varying E/I ratios maintain stable bursting dynamics. Interburst intervals remain constant for most ratios, except in the extremes of 0 to 10% and 90 to 100% inhibitory cells. Single-cell recordings and modeling suggest that networks adapt to chronic alterations of E/I compositions by balancing E/I connectivity. Gradual blockade of inhibition substantiates the agreement between the model and experiment and defines its limits. Combining measurements of population and single-cell activity with theoretical modeling, we provide a clearer picture of how E/I balance is preserved and where it fails in living neuronal networks.

Altered Fast Synaptic Transmission in a Mouse Model of DNM1-Associated Developmental Epileptic Encephalopathy

Developmental epileptic encephalopathies (DEEs) are severe seizure disorders that occur in infants and young children, characterized by developmental delay, cognitive decline, and early mortality. Recent efforts have identified a wide variety of genetic variants that cause DEEs. Among these, variants in the DNM1 gene have emerged as definitive causes of DEEs, including infantile spasms and Lennox–Gastaut syndrome. A mouse model of Dnm1-associated DEE, known as “Fitful” (Dnm1^Ftfl), recapitulates key features of the disease, including spontaneous seizures, early lethality, and neuronal degeneration. Previous work showed that DNM1 is a key regulator of synaptic vesicle (SV) endocytosis and synaptic transmission and suggested that inhibitory neurotransmission may be more reliant on DNM1 function than excitatory transmission. The Dnm1^Ftfl variant is thought to encode a dominant negative DNM1 protein; however, the effects of the Dnm1^Ftfl variant on synaptic transmission are largely unknown. To understand these synaptic effects, we recorded from pairs of cultured mouse cortical neurons and characterized all four major connection types [excitation of excitation (E-E), inhibition of inhibition (I-I), E-I, I-E]. Miniature and spontaneous EPSCs and IPSCs were larger, but less frequent, at all Dnm1^Ftfl synaptic types, and Dnm1^Ftfl neurons had reduced expression of excitatory and inhibitory SV markers. Baseline evoked transmission, however, was reduced only at inhibitory synapses onto excitatory neurons, because of a smaller pool of releasable SVs. In addition to these synaptic alterations, Dnm1^Ftfl neurons degenerated later in development, although their activity levels were reduced, suggesting that Dnm1^Ftfl may impair synaptic transmission and neuronal health through distinct mechanisms.

Antiphase synchronization in multiplex networks with attractive and repulsive interactions

A series of recent publications, within the framework of network science, have focused on the coexistence of mixed attractive and repulsive (excitatory and inhibitory) interactions among the units within the same system, motivated by the analogies with spin glasses as well as to neural networks, or ecological systems. However, most of these investigations have been restricted to single layer networks, requiring further analysis of the complex dynamics and particular equilibrium states that emerge in multilayer configurations. This article investigates the synchronization properties of dynamical systems connected through multiplex architectures in the presence of attractive intralayer and repulsive interlayer connections. This setting enables the emergence of antisynchronization, i.e., intralayer synchronization coexisting with antiphase dynamics between coupled systems of different layers. We demonstrate the existence of a transition from interlayer antisynchronization to antiphase synchrony in any connected bipartite multiplex architecture when the repulsive coupling is introduced through any spanning tree of a single layer. We identify, analytically, the required graph topologies for interlayer antisynchronization and its interplay with intralayer and antiphase synchronization. Next, we analytically derive the invariance of intralayer synchronization manifold and calculate the attractor size of each oscillator exhibiting interlayer antisynchronization together with intralayer synchronization. The necessary conditions for the existence of interlayer antisynchronization along with intralayer synchronization are given and numerically validated by considering Stuart-Landau oscillators. Finally, we also analytically derive the local stability condition of the interlayer antisynchronization state using the master stability function approach.

Functional strengthening through synaptic scaling upon connectivity disruption in neuronal cultures

An elusive phenomenon in network neuroscience is the extent of neuronal activity remodeling upon damage. Here, we investigate the action of gradual synaptic blockade on the effective connectivity in cortical networks in vitro. We use two neuronal cultures configurations-one formed by about 130 neuronal aggregates and another one formed by about 600 individual neurons-and monitor their spontaneous activity upon progressive weakening of excitatory connectivity. We report that the effective connectivity in all cultures exhibits a first phase of transient strengthening followed by a second phase of steady deterioration. We quantify these phases by measuring GEFF, the global efficiency in processing network information. We term hyperefficiency the sudden strengthening of GEFF upon network deterioration, which increases by 20-50% depending on culture type. Relying on numerical simulations we reveal the role of synaptic scaling, an activity-dependent mechanism for synaptic plasticity, in counteracting the perturbative action, neatly reproducing the observed hyperefficiency. Our results demonstrate the importance of synaptic scaling as resilience mechanism.

Turing instability in the reaction-diffusion network

It is an established fact that a positive wave number plays an essential role in Turing instability. However, the impact of a negative wave number on Turing instability remains unclear. Here, we investigate the effect of the weights and nodes on Turing instability in the FitzHugh-Nagumo model, and theoretical results reveal genesis of Turing instability due to a negative wave number through the stability analysis and mean-field method. We obtain the Turing instability region in the continuous media system and provide the relationship between degree and eigenvalue of the network matrix by the Gershgorin circle theorem. Furthermore, the Turing instability condition about nodes and the weights is provided in the network-organized system. Additionally, we found chaotic behavior because of interactions between I Turing instability and II Turing instability. Besides, we apply this above analysis to explaining the mechanism of the signal conduction of the inhibitory neuron. We find a moderate coupling strength and corresponding number of links are necessary to the signal conduction.

Impact of Physical Obstacles on the Structural and Effective Connectivity of in silico Neuronal Circuits

Scaffolds and patterned substrates are among the most successful strategies to dictate the connectivity between neurons in culture. Here, we used numerical simulations to investigate the capacity of physical obstacles placed on a flat substrate to shape structural connectivity, and in turn collective dynamics and effective connectivity, in biologically-realistic neuronal networks. We considered μ-sized obstacles placed in mm-sized networks. Three main obstacle shapes were explored, namely crosses, circles and triangles of isosceles profile. They occupied either a small area fraction of the substrate or populated it entirely in a periodic manner. From the point of view of structure, all obstacles promoted short length-scale connections, shifted the in- and out-degree distributions toward lower values, and increased the modularity of the networks. The capacity of obstacles to shape distinct structural traits depended on their density and the ratio between axonal length and substrate diameter. For high densities, different features were triggered depending on obstacle shape, with crosses trapping axons in their vicinity and triangles funneling axons along the reverse direction of their tip. From the point of view of dynamics, obstacles reduced the capacity of networks to spontaneously activate, with triangles in turn strongly dictating the direction of activity propagation. Effective connectivity networks, inferred using transfer entropy, exhibited distinct modular traits, indicating that the presence of obstacles facilitated the formation of local effective microcircuits. Our study illustrates the potential of physical constraints to shape structural blueprints and remodel collective activity, and may guide investigations aimed at mimicking organizational traits of biological neuronal circuits.

A 5' UTR GGN repeat controls localisation and translation of a potassium leak channel mRNA through G-quadruplex formation

RNA G-quadruplexes (G4s) are secondary structures proposed to function as regulators of post-transcriptional mRNA localisation and translation. G4s within some neuronal mRNAs are known to control distal localisation and local translation, contributing to distinct local proteomes that facilitate the synaptic remodelling attributed to normal cellular function. In this study, we characterise the G4 formation of a (GGN)13 repeat found within the 5' UTR of the potassium 2-pore domain leak channel Task3 mRNA. Biophysical analyses show that this (GGN)13 repeat forms a parallel G4 in vitro exhibiting the stereotypical potassium specificity of G4s, remaining thermostable under physiological ionic conditions. Through mouse brain tissue G4-RNA immunoprecipitation, we further confirm that Task3 mRNA forms a G4 structure in vivo. The G4 is inhibitory to translation of Task3 in vitro and is overcome through activity of a G4-specific helicase DHX36, increasing K+ leak currents and membrane hyperpolarisation in HEK293 cells. Further, we observe that this G4 is fundamental to ensuring delivery of Task3 mRNA to distal primary cortical neurites. It has been shown that aberrant Task3 expression correlates with neuronal dysfunction, we therefore posit that this G4 is important in regulated local expression of Task3 leak channels that maintain K+ leak within neurons.

Neurodegeneration exposes firing rate dependent effects on oscillation dynamics in computational neural networks

Traumatic brain injury (TBI) can lead to neurodegeneration in the injured circuitry, either through primary structural damage to the neuron or secondary effects that disrupt key cellular processes. Moreover, traumatic injuries can preferentially impact subpopulations of neurons, but the functional network effects of these targeted degeneration profiles remain unclear. Although isolating the consequences of complex injury dynamics and long-term recovery of the circuit can be difficult to control experimentally, computational networks can be a powerful tool to analyze the consequences of injury. Here, we use the Izhikevich spiking neuron model to create networks representative of cortical tissue. After an initial settling period with spike-timing-dependent plasticity (STDP), networks developed rhythmic oscillations similar to those seen in vivo. As neurons were sequentially removed from the network, population activity rate and oscillation dynamics were significantly reduced. In a successive period of network restructuring with STDP, network activity levels returned to baseline for some injury levels and oscillation dynamics significantly improved. We next explored the role that specific neurons have in the creation and termination of oscillation dynamics. We determined that oscillations initiate from activation of low firing rate neurons with limited structural inputs. To terminate oscillations, high activity excitatory neurons with strong input connectivity activate downstream inhibitory circuitry. Finally, we confirm the excitatory neuron population role through targeted neurodegeneration. These results suggest targeted neurodegeneration can play a key role in the oscillation dynamics after injury.

The Epilepsy-Related Protein PCDH19 Regulates Tonic Inhibition, GABAAR Kinetics, and the Intrinsic Excitability of Hippocampal Neurons

PCDH19 encodes for protocadherin-19 (PCDH19), a cell-adhesion molecule of the cadherin superfamily preferentially expressed in the brain. PCDH19 mutations cause a neurodevelopmental syndrome named epileptic encephalopathy, early infantile, 9 (EIEE9) characterized by seizures associated with cognitive and behavioral deficits. We recently reported that PCDH19 binds the alpha subunits of GABA_A receptors (GABA_ARs), modulating their surface availability and miniature inhibitory postsynaptic currents (mIPSCs). Here, we investigated whether PCDH19 regulatory function on GABA_ARs extends to the extrasynaptic receptor pool that mediates tonic current. In fact, the latter shapes neuronal excitability and network properties at the base of information processing. By combining patch-clamp recordings in whole-cell and cell-attached configurations, we provided a functional characterization of primary hippocampal neurons from embryonic rats of either sex expressing a specific PCDH19 short hairpin (sh)RNA. We first demonstrated that PCDH19 downregulation reduces GABA_AR-mediated tonic current, evaluated by current shift and baseline noise analysis. Next, by single-channel recordings, we showed that PCDH19 regulates GABA_ARs kinetics without altering their conductance. In particular, GABA_ARs of shRNA-expressing neurons preferentially exhibit brief openings at the expense of long ones, thus displaying a flickering behavior. Finally, we showed that PCDH19 downregulation reduces the rheobase and increases the frequency of action potential firing, thus indicating neuronal hyperexcitability. These findings establish PCDH19 as a critical determinant of GABA_AR-mediated tonic transmission and GABA_ARs gating, and provide the first mechanistic insights into PCDH19-related hyperexcitability and comorbidities.

Astrocytic Regulation of Synchronous Bursting in Cortical Cultures: From Local to Global

Synchronous bursting (SB) is ubiquitous in neuronal networks and independent of network structure. Although it is known to be driven by glutamatergic neurotransmissions, its underlying mechanism remains unclear. Recent studies show that local glutamate recycle by astrocytes affects nearby neuronal activities, which indicate that the local dynamics might also be the origin of SBs in networks. We investigated the effects of local glutamate dynamics on SBs in both cultures developed on multielectrode array (MEA) systems and a tripartite synapse simulation. Local glutamate uptake by astrocytes was altered by pharmacological targeting of GLT-1 glutamate transporters, whereas neuronal firing activities and synaptic glutamate level was simultaneously monitored with MEA and astrocyte-specific glutamate sensors (intensity-based glutamate-sensing fluorescent reporter), respectively. Global SB properties were significantly altered on targeting GLT-1. Detailed simulation of a network with astrocytic glutamate uptake and recycle mechanisms, conforming with the experimental observations, shows that astrocytes function as a slow negative feedback to neuronal activities in the network. SB in the network can be realized as an alternation between positive and negative feedback in the neurons and astrocytes, respectively. An understanding of glutamate trafficking dynamics is of general application to explain how astrocyte malfunction can result in pathological seizure-like phenomena in neuronal systems.

Deficits in coordinated neuronal activity and network topology are striatal hallmarks in Huntington's disease

Background

Network alterations underlying neurodegenerative diseases often precede symptoms and functional deficits. Thus, their early identification is central for improved prognosis. In Huntington’s disease (HD), the cortico-striatal networks, involved in motor function processing, are the most compromised neural substrate. However, whether the network alterations are intrinsic of the striatum or the cortex is not fully understood.

Results

In order to identify early HD neural deficits, we characterized neuronal ensemble calcium activity and network topology of HD striatal and cortical cultures. We used large-scale calcium imaging combined with activity-based network inference analysis. We extracted collective activity events and inferred the topology of the neuronal network in cortical and striatal primary cultures from wild-type and R6/1 mouse model of HD. Striatal, but not cortical, HD networks displayed lower activity and a lessened ability to integrate information. GABA_A receptor blockade in healthy and HD striatal cultures generated similar coordinated ensemble activity and network topology, highlighting that the excitatory component of striatal system is spared in HD. Conversely, NMDA receptor activation increased individual neuronal activity while coordinated activity became highly variable and undefined. Interestingly, by boosting NMDA activity, we rectified striatal HD network alterations.

Conclusions

Overall, our integrative approach highlights striatal defective network integration capacity as a major contributor of basal ganglia dysfunction in HD and suggests that increased excitatory drive may serve as a potential intervention. In addition, our work provides a valuable tool to evaluate in vitro network recovery after treatment intervention in basal ganglia disorders.

Suppression of hypersynchronous network activity in cultured cortical neurons using an ultrasoft silicone scaffold

The spontaneous activity pattern of cortical neurons in dissociated culture is characterized by burst firing that is highly synchronized among a wide population of cells. The degree of synchrony, however, is excessively higher than that in cortical tissues. Here, we employed polydimethylsiloxane (PDMS) elastomers to establish a novel system for culturing neurons on a scaffold with an elastic modulus resembling brain tissue, and investigated the effect of the scaffold's elasticity on network activity patterns in cultured rat cortical neurons. Using whole-cell patch clamp to assess the scaffold effect on the development of synaptic connections, we found that the amplitude of excitatory postsynaptic current, as well as the frequency of spontaneous transmissions, was reduced in neuronal networks grown on an ultrasoft PDMS with an elastic modulus of 0.5 kPa. Furthermore, the ultrasoft scaffold was found to suppress neural correlations in the spontaneous activity of the cultured neuronal network. The dose of GsMTx-4, an antagonist of stretch-activated cation channels (SACs), required to reduce the generation of the events below 1.0 event per min on PDMS substrates was lower than that for neurons on a glass substrate. This suggests that the difference in the baseline level of SAC activation is a molecular mechanism underlying the alteration in neuronal network activity depending on scaffold stiffness. Our results demonstrate the potential application of PDMS with biomimetic elasticity as cell-culture scaffold for bridging the in vivo-in vitro gap in neuronal systems.

Neuronal Degeneration Impairs Rhythms Between Connected Microcircuits

Synchronization of neural activity across brain regions is critical to processes that include perception, learning, and memory. After traumatic brain injury (TBI), neuronal degeneration is one possible effect and can alter communication between neural circuits. Consequently, synchronization between neurons may change and can contribute to both lasting changes in functional brain networks and cognitive impairment in patients. However, fundamental principles relating exactly how TBI at the cellular scale affects synchronization of mesoscale circuits are not well understood. In this work, we use computational networks of Izhikevich integrate-and-fire neurons to study synchronized, oscillatory activity between clusters of neurons, which also adapt according to spike-timing-dependent plasticity (STDP). We study how the connections within and between these neuronal clusters change as unidirectional connections form between the two neuronal populations. In turn, we examine how neuronal deletion, intended to mimic the temporary or permanent loss of neurons in the mesoscale circuit, affects these dynamics. We determine synchronization of two neuronal circuits requires very modest connectivity between these populations; approximately 10% of neurons projecting from one circuit to another circuit will result in high synchronization. In addition, we find that synchronization level inversely affects the strength of connection between neuronal microcircuits - moderately synchronized microcircuits develop stronger intercluster connections than do highly synchronized circuits. Finally, we find that highly synchronized circuits are largely protected against the effects of neuronal deletion but may display changes in frequency properties across circuits with targeted neuronal loss. Together, our results suggest that strongly and weakly connected regions differ in their inherent resilience to damage and may serve different roles in a larger network.

Jensen's force and the statistical mechanics of cortical asynchronous states

Cortical networks are shaped by the combined action of excitatory and inhibitory interactions. Among other important functions, inhibition solves the problem of the all-or-none type of response that comes about in purely excitatory networks, allowing the network to operate in regimes of moderate or low activity, between quiescent and saturated regimes. Here, we elucidate a noise-induced effect that we call “Jensen’s force” –stemming from the combined effect of excitation/inhibition balance and network sparsity– which is responsible for generating a phase of self-sustained low activity in excitation-inhibition networks. The uncovered phase reproduces the main empirically-observed features of cortical networks in the so-called asynchronous state, characterized by low, un-correlated and highly-irregular activity. The parsimonious model analyzed here allows us to resolve a number of long-standing issues, such as proving that activity can be self-sustained even in the complete absence of external stimuli or driving. The simplicity of our approach allows for a deep understanding of asynchronous states and of the phase transitions to other standard phases it exhibits, opening the door to reconcile, asynchronous-state and critical-state hypotheses, putting them within a unified framework. We argue that Jensen’s forces are measurable experimentally and might be relevant in contexts beyond neuroscience.

Self-organization of modular network architecture by activity-dependent neuronal migration and outgrowth

The spatial distribution of neurons and activity-dependent neurite outgrowth shape long-range interaction, recurrent local connectivity and the modularity in neuronal networks. We investigated how this mesoscale architecture develops by interaction of neurite outgrowth, cell migration and activity in cultured networks of rat cortical neurons and show that simple rules can explain variations of network modularity. In contrast to theoretical studies on activity-dependent outgrowth but consistent with predictions for modular networks, spontaneous activity and the rate of synchronized bursts increased with clustering, whereas peak firing rates in bursts increased in highly interconnected homogeneous networks. As Ca²⁺ influx increased exponentially with increasing network recruitment during bursts, its modulation was highly correlated to peak firing rates. During network maturation, long-term estimates of Ca²⁺ influx showed convergence, even for highly different mesoscale architectures, neurite extent, connectivity, modularity and average activity levels, indicating homeostatic regulation towards a common set-point of Ca²⁺ influx.

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.

Chimera states in coupled logistic maps with additional weak nonlocal topology

We demonstrate the occurrence of coexisting domains of partially coherent and incoherent patterns or simply known as chimera states in a network of globally coupled logistic maps upon addition of weak nonlocal topology. We find that the chimera states survive even after we disconnect nonlocal connections of some of the nodes in the network. Also, we show that the chimera states exist when we introduce symmetric gaps in the nonlocal coupling between predetermined nodes. We ascertain our results, for the existence of chimera states, by carrying out the recurrence quantification analysis and by computing the strength of incoherence. We extend our analysis for the case of small-world networks of coupled logistic maps and found the emergence of chimeralike states under the influence of weak nonlocal topology.

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.

Statistical inference approach to structural reconstruction of complex networks from binary time series

Complex networks hosting binary-state dynamics arise in a variety of contexts. In spite of previous works, to fully reconstruct the network structure from observed binary data remains challenging. We articulate a statistical inference based approach to this problem. In particular, exploiting the expectation-maximization (EM) algorithm, we develop a method to ascertain the neighbors of any node in the network based solely on binary data, thereby recovering the full topology of the network. A key ingredient of our method is the maximum-likelihood estimation of the probabilities associated with actual or nonexistent links, and we show that the EM algorithm can distinguish the two kinds of probability values without any ambiguity, insofar as the length of the available binary time series is reasonably long. Our method does not require any a priori knowledge of the detailed dynamical processes, is parameter-free, and is capable of accurate reconstruction even in the presence of noise. We demonstrate the method using combinations of distinct types of binary dynamical processes and network topologies, and provide a physical understanding of the underlying reconstruction mechanism. Our statistical inference based reconstruction method contributes an additional piece to the rapidly expanding "toolbox" of data based reverse engineering of complex networked systems.

Dependence and Homeostasis of Membrane Impedance on Cell Morphology in Cultured Hippocampal Neurons

The electrical impedance of cell membranes is important for excitable cells, such as neurons, because it strongly influences the amount of membrane potential change upon a flow of ionic current across the membrane. Here, we report on an investigation of how neuronal morphology affects membrane impedance of cultured hippocampal neurons. Microfabricated substrates with patterned scaffolding molecules were used to restrict the neurite growth of hippocampal neurons, and the impedance was measured via whole-cell patch-clamp recording under the inhibition of voltage-dependent ion channels. Membrane impedance was found to depend inversely on the dendrite length and soma area, as would be expected from the fact that its electrical property is equivalent to a parallel RC circuit. Moreover, we found that in biological neurons, the membrane impedance is homeostatically regulated to impede changes in the membrane area. The findings provide direct evidence on cell-autonomous regulation of neuronal impedance and pave the way towards elucidating the mechanism responsible for the resilience of biological neuronal networks.

Can a time varying external drive give rise to apparent criticality in neural systems?

The finding of power law scaling in neural recordings lends support to the hypothesis of critical brain dynamics. However, power laws are not unique to critical systems and can arise from alternative mechanisms. Here, we investigate whether a common time-varying external drive to a set of Poisson units can give rise to neuronal avalanches and exhibit apparent criticality. To this end, we analytically derive the avalanche size and duration distributions, as well as additional measures, first for homogeneous Poisson activity, and then for slowly varying inhomogeneous Poisson activity. We show that homogeneous Poisson activity cannot give rise to power law distributions. Inhomogeneous activity can also not generate perfect power laws, but it can exhibit approximate power laws with cutoffs that are comparable to those typically observed in experiments. The mechanism of generating apparent criticality by time-varying external fields, forces or input may generalize to many other systems like dynamics of swarms, diseases or extinction cascades. Here, we illustrate the analytically derived effects for spike recordings in vivo and discuss approaches to distinguish true from apparent criticality. Ultimately, this requires causal interventions, which allow separating internal system properties from externally imposed ones.

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.

Hierarchy measurement for modeling network dynamics under directed attacks

A fundamental issue in the dynamics of complex systems is the resilience of the network in response to targeted attacks. This paper explores the local dynamics of the network attack process by investigating the order of removal of the nodes that have maximal degree, and shows that this dynamic network response can be predicted from the graph's initial connectivity. We demonstrate numerically that the maximal degree M(τ) of the network at time step τ decays exponentially with τ via a topology-dependent exponent. Moreover, the order in which sites are removed can be approximated by considering the network's "hierarchy" function h, which measures for each node V_{i} how many of its initial nearest neighbors have lower degree versus those that have a higher one. Finally, we show that the exponents we identified for the attack dynamics are related to the exponential behavior of spreading activation dynamics. The results suggest that the function h, which has both local and global properties, is a novel nodal measurement for network dynamics and structure.

Mapping functional connectivity of bursting neuronal networks

Using single-cell laser scanning photostimulation (LSPS) combined with broad-field calcium imaging, we measured the functional connectivity of neuronal cultures before and after the developmental appearance of network bursting. From these data, network properties were determined for these relatively large neuronal networks. Based on these properties, we found that although 'small-world' network behavior existed throughout this time period, only average node degree and global efficiency correlate with the development of network bursting while clustering and local efficiency remained relatively constant.

Noise focusing in neuronal tissues: Symmetry breaking and localization in excitable networks with quenched disorder

We introduce a coarse-grained stochastic model for the spontaneous activity of neuronal cultures to explain the phenomenon of noise focusing, which entails localization of the noise activity in excitable networks with metric correlations. The system is modeled as a continuum excitable medium with a state-dependent spatial coupling that accounts for the dynamics of synaptic connections. The most salient feature is the emergence at the mesoscale of a vector field V(r), which acts as an advective carrier of the noise. This entails an explicit symmetry breaking of isotropy and homogeneity that stems from the amplification of the quenched fluctuations of the network by the activity avalanches, concomitant with the excitable dynamics. We discuss the microscopic interpretation of V(r) and propose an explicit construction of it. The coarse-grained model shows excellent agreement with simulations at the network level. The generic nature of the observed phenomena is discussed.

Dominance of Metric Correlations in Two-Dimensional Neuronal Cultures Described through a Random Field Ising Model

We introduce a novel random field Ising model, grounded on experimental observations, to assess the importance of metric correlations in cortical circuits in vitro. Metric correlations arise from both the finite axonal length and the heterogeneity in the spatial arrangement of neurons. The experiments consider the response of neuronal cultures to an external electric stimulation for a gradually weaker connectivity strength between neurons, and in cultures with different spatial configurations. The model can be analytically solved in the metric-free, mean-field scenario. The presence of metric correlations precipitates a strong deviation from the mean field. Null models of the same networks that preserve the distribution of connections recover the mean field. Our results show that metric-inherited correlations in spatial networks dominate the connectivity blueprint, mask the actual distribution of connections, and may emerge as the asset that shapes network dynamics.

External Excitation of Neurons Using Electric and Magnetic Fields in One- and Two-dimensional Cultures

A neuron will fire an action potential when its membrane potential exceeds a certain threshold. In typical activity of the brain, this occurs as a result of chemical inputs to its synapses. However, neurons can also be excited by an imposed electric field. In particular, recent clinical applications activate neurons by creating an electric field externally. It is therefore of interest to investigate how the neuron responds to the external field and what causes the action potential. Fortunately, precise and controlled application of an external electric field is possible for embryonic neuronal cells that are excised, dissociated and grown in cultures. This allows the investigation of these questions in a highly reproducible system. In this paper some of the techniques used for controlled application of external electric field on neuronal cultures are reviewed. The networks can be either one dimensional, i.e. patterned in linear forms or allowed to grow on the whole plane of the substrate, and thus two dimensional. Furthermore, the excitation can be created by the direct application of electric field via electrodes immersed in the fluid (bath electrodes) or by inducing the electric field using the remote creation of magnetic pulses.

A spatially resolved network spike in model neuronal cultures reveals nucleation centers, circular traveling waves and drifting spiral waves

We show that in model neuronal cultures, where the probability of interneuronal connection formation decreases exponentially with increasing distance between the neurons, there exists a small number of spatial nucleation centers of a network spike, from where the synchronous spiking activity starts propagating in the network typically in the form of circular traveling waves. The number of nucleation centers and their spatial locations are unique and unchanged for a given realization of neuronal network but are different for different networks. In contrast, if the probability of interneuronal connection formation is independent of the distance between neurons, then the nucleation centers do not arise and the synchronization of spiking activity during a network spike occurs spatially uniform throughout the network. Therefore one can conclude that spatial proximity of connections between neurons is important for the formation of nucleation centers. It is also shown that fluctuations of the spatial density of neurons at their random homogeneous distribution typical for the experiments in vitro do not determine the locations of the nucleation centers. The simulation results are qualitatively consistent with the experimental observations.