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.

A novel methodology to describe neuronal networks activity reveals spatiotemporal recruitment dynamics of synchronous bursting states.

We propose a novel phase based analysis with the purpose of quantifying the periodic bursts of activity observed in various neuronal systems. The way bursts are intiated and propagate in a spatial network is still insufficiently characterized. In particular, we investigate here how these spatiotemporal dynamics depend on the mean connection length. We use a simplified description of a neuron's state as a time varying phase between firings. This leads to a definition of network bursts, that does not depend on the practitioner's individual judgment as the usage of subjective thresholds and time scales. This allows both an easy and objective characterization of the bursting dynamics, only depending on system's proper scales. Our approach thus ensures more reliable and reproducible measurements. We here use it to describe the spatiotemporal processes in networks of intrinsically oscillating neurons. The analysis rigorously reveals the role of the mean connectivity length in spatially embedded networks in determining the existence of "leader" neurons during burst initiation, a feature incompletely understood observed in several neuronal cultures experiments. The precise definition of a burst with our method allowed us to rigorously characterize the initiation dynamics of bursts and show how it depends on the mean connectivity length. Although presented with simulations, the methodology can be applied to other forms of neuronal spatiotemporal data. As shown in a preliminary study with MEA recordings, it is not limited to in silico modeling.

Understanding the Generation of Network Bursts by Adaptive Oscillatory Neurons.

Experimental and numerical studies have revealed that isolated populations of oscillatory neurons can spontaneously synchronize and generate periodic bursts involving the whole network. Such a behavior has notably been observed for cultured neurons in rodent's cortex or hippocampus. We show here that a sufficient condition for this network bursting is the presence of an excitatory population of oscillatory neurons which displays spike-driven adaptation. We provide an analytic model to analyze network bursts generated by coupled adaptive exponential integrate-and-fire neurons. We show that, for strong synaptic coupling, intrinsically tonic spiking neurons evolve to reach a synchronized intermittent bursting state. The presence of inhibitory neurons or plastic synapses can then modulate this dynamics in many ways but is not necessary for its appearance. Thanks to a simple self-consistent equation, our model gives an intuitive and semi-quantitative tool to understand the bursting behavior. Furthermore, it suggests that after-hyperpolarization currents are sufficient to explain bursting termination. Through a thorough mapping between the theoretical parameters and ion-channel properties, we discuss the biological mechanisms that could be involved and the relevance of the explored parameter-space. Such an insight enables us to propose experimentally-testable predictions regarding how blocking fast, medium or slow after-hyperpolarization channels would affect the firing rate and burst duration, as well as the interburst interval.

Effect of threshold disorder on the quorum percolation model.

We study the modifications induced in the behavior of the quorum percolation model on neural networks with Gaussian in-degree by taking into account an uncorrelated Gaussian thresholds variability. We derive a mean-field approach and show its relevance by carrying out explicit Monte Carlo simulations. It turns out that such a disorder shifts the position of the percolation transition, impacts the size of the giant cluster, and can even destroy the transition. Moreover, we highlight the occurrence of disorder independent fixed points above the quorum critical value. The mean-field approach enables us to interpret these effects in terms of activation probability. A finite-size analysis enables us to show that the order parameter is weakly self-averaging with an exponent independent on the thresholds disorder. Last, we show that the effects of the thresholds and connectivity disorders cannot be easily discriminated from the measured averaged physical quantities.

Combining microfluidics, optogenetics and calcium imaging to study neuronal communication in vitro.

In this paper we report the combination of microfluidics, optogenetics and calcium imaging as a cheap and convenient platform to study synaptic communication between neuronal populations in vitro. We first show that Calcium Orange indicator is compatible in vitro with a commonly used Channelrhodopsine-2 (ChR2) variant, as standard calcium imaging conditions did not alter significantly the activity of transduced cultures of rodent primary neurons. A fast, robust and scalable process for micro-chip fabrication was developed in parallel to build micro-compartmented cultures. Coupling optical fibers to each micro-compartment allowed for the independent control of ChR2 activation in the different populations without crosstalk. By analyzing the post-stimuli activity across the different populations, we finally show how this platform can be used to evaluate quantitatively the effective connectivity between connected neuronal populations.

Memory decay and loss of criticality in quorum percolation.

In this paper, we present the effects of memory decay on a bootstrap percolation model applied to random directed graphs (quorum percolation). The addition of decay was motivated by its natural occurrence in physical systems previously described by percolation theory, such as cultured neuronal networks, where decay originates from ionic leakage through the membrane of neurons and/or synaptic depression. Surprisingly, this feature alone appears to change the critical behavior of the percolation transition, where discontinuities are replaced by steep but finite slopes. Using different numerical approaches, we show evidence for this qualitative change even for very small decay values. In experiments where the steepest slopes can not be resolved and still appear as discontinuities, decay produces nonetheless a quantitative difference on the location of the apparent critical point. We discuss how this shift impacts network connectivity previously estimated without considering decay. In addition to this particular example, we believe that other percolation models are worth reinvestigating, taking into account similar sorts of memory decay.

Effects of random and deterministic discrete scale invariance on the critical behavior of the Potts model.

The effects of disorder on the critical behavior of the q-state Potts model in noninteger dimensions are studied by comparison of deterministic and random fractals sharing the same dimensions in the framework of a discrete scale invariance. We carried out intensive Monte Carlo simulations. In the case of a fractal dimension slightly smaller than two d(f) ~/= 1.974636, we give evidence that the disorder structured by discrete scale invariance does not change the first order transition associated with the deterministic case when q = 7. Furthermore the study of the high value q = 14 shows that the transition is a second order one both for deterministic and random scale invariance, but that their behavior belongs to different university classes.

Critical behavior of the Ising model on random fractals.

We study the critical behavior of the Ising model in the case of quenched disorder constrained by fractality on random Sierpinski fractals with a Hausdorff dimension d(f) is approximately equal to 1.8928. This is a first attempt to study a situation between the borderline cases of deterministic self-similarity and quenched randomness. Intensive Monte Carlo simulations were carried out. Scaling corrections are much weaker than in the deterministic cases, so that our results enable us to ensure that finite-size scaling holds, and that the critical behavior is described by a new universality class. The hyperscaling relation is compatible with an effective dimension equal to the Hausdorff one; moreover the two eigenvalues exponents of the renormalization flows are shown to be different from the ones calculated from ε expansions, and from the ones obtained for fourfold symmetric deterministic fractals. Although the space dimensionality is not integer, lack of self-averaging properties exhibits some features very close to the ones of a random fixed point associated with a relevant disorder.