Emergence of assortative mixing between clusters of cultured neurons

The analysis of the activity of neuronal cultures is considered to be a good proxy of the functional connectivity of in vivo neuronal tissues. Thus, the functional complex network inferred from activity patterns is a promising way to unravel the interplay between structure and functionality of neuronal systems. Here, we monitor the spontaneous self-sustained dynamics in neuronal cultures formed by interconnected aggregates of neurons (clusters). Dynamics is characterized by the fast activation of groups of clusters in sequences termed bursts. The analysis of the time delays between clusters' activations within the bursts allows the reconstruction of the directed functional connectivity of the network. We propose a method to statistically infer this connectivity and analyze the resulting properties of the associated complex networks. Surprisingly enough, in contrast to what has been reported for many biological networks, the clustered neuronal cultures present assortative mixing connectivity values, meaning that there is a preference for clusters to link to other clusters that share similar functional connectivity, as well as a rich-club core, which shapes a"connectivity backbone" in the network. These results point out that the grouping of neurons and the assortative connectivity between clusters are intrinsic survival mechanisms of the culture

Emergence of Small-World Anatomical Networks in Self-Organizing Clustered Neuronal Cultures

In vitro primary cultures of dissociated invertebrate neurons from locust ganglia are used to experimentally investigate the morphological evolution of assemblies of living neurons, as they self-organize from collections of separated cells into elaborated, clustered, networks. At all the different stages of the culture's development, identification of neurons' and neurites' location by means of a dedicated software allows to ultimately extract an adjacency matrix from each image of the culture. In turn, a systematic statistical analysis of a group of topological observables grants us the possibility of quantifying and tracking the progression of the main network's characteristics during the self-organization process of the culture. Our results point to the existence of a particular state corresponding to a small-world network configuration, in which several relevant graph's micro- and meso-scale properties emerge. Finally, we identify the main physical processes ruling the culture's morphological transformations, and embed them into a simplified growth model qualitatively reproducing the overall set of experimental observations

Emergence of Modular Structure in a Large-Scale Brain Network with Interactions between Dynamics and Connectivity

A network of 32 or 64 connected neural masses, each representing a large population of interacting excitatory and inhibitory neurons and generating an electroencephalography/magnetoencephalography like output signal, was used to demonstrate how an interaction between dynamics and connectivity might explain the emergence of complex network features, in particular modularity. Network evolution was modeled by two processes: (i) synchronization dependent plasticity (SDP) and (ii) growth dependent plasticity (GDP). In the case of SDP, connections between neural masses were strengthened when they were strongly synchronized, and were weakened when they were not. GDP was modeled as a homeostatic process with random, distance dependent outgrowth of new connections between neural masses. GDP alone resulted in stable networks with distance dependent connection strengths, typical small-world features, but no degree correlations and only weak modularity. SDP applied to random networks induced clustering, but no clear modules. Stronger modularity evolved only through an interaction of SDP and GDP, with the number and size of the modules depending on the relative strength of both processes, as well as on the size of the network. Lesioning part of the network, after a stable state was achieved, resulted in a temporary disruption of the network structure. The model gives a possible scenario to explain how modularity can arise in developing brain networks, and makes predictions about the time course of network changes during development and following acute lesions

Magnetite-Amyloid-β deteriorates activity and functional organization in an in vitro model for Alzheimer's disease

The understanding of the key mechanisms behind human brain deterioration in Alzheimer' disease (AD) is a highly active field of research. The most widespread hypothesis considers a cascade of events initiated by amyloid-β peptide fibrils that ultimately lead to the formation of the lethal amyloid plaques. Recent studies have shown that other agents, in particular magnetite, can also play a pivotal role. To shed light on the action of magnetite and amyloid-β in the deterioration of neuronal circuits, we investigated their capacity to alter spontaneous activity patterns in cultured neuronal networks. Using a versatile experimental platform that allows the parallel monitoring of several cultures, the activity in controls was compared with the one in cultures dosed with magnetite, amyloid-β and magnetite-amyloid-β complex. A prominent degradation in spontaneous activity was observed solely when amyloid-β and magnetite acted together. Our work suggests that magnetite nanoparticles have a more prominent role in AD than previously thought, and may bring new insights in the understanding of the damaging action of magnetite-amyloid-β complex. Our experimental system also offers new interesting perspectives to explore key biochemical players in neurological disorders through a controlled, model system manner

Linear and nonlinear approaches to unravel dynamics and connectivity in neuronal cultures

[eng] In the present thesis, we propose to explore neuronal circuits at the mesoscale, an approach in which one monitors small populations of few thousand neurons and concentrates in the emergence of collective behavior. In our case, we carried out such an exploration both experimentally and numerically, and by adopting an analysis perspective centered on time series analysis and dynamical systems. Experimentally, we used neuronal cultures and prepared more than 200 of them, which were monitored using fluorescence calcium imaging. By adjusting the experimental conditions, we could set two basic arrangements of neurons, namely homogeneous and aggregated. In the experiments, we carried out two major explorations, namely development and disintegration. In the former we investigated changes in network behavior as it matured; in the latter we applied a drug that reduced neuronal interconnectivity. All the subsequent analyses and modeling along the thesis are based on these experimental data. Numerically, the thesis comprised two aspects. The first one was oriented towards a simulation of neuronal connectivity and dynamics. The second one was oriented towards the development of linear and nonlinear analysis tools to unravel dynamic and connectivity aspects of the measured experimental networks. For the first aspect, we developed a sophisticated software package to simulate single neuronal dynamics using a quadratic integrate–and–fire model with adaptation and depression. This model was plug into a synthetic graph in which the nodes of the network are neurons, and the edges connections. The graph was created using spatial embedding and realistic biology. We carried out hundreds of simulations in which we tuned the density of neurons, their spatial arrangement and the characteristics of the fluorescence signal. As a key result, we observed that homogeneous networks required a substantial number of neurons to fire and exhibit collective dynamics, and that the presence of aggregation significantly reduced the number of required neurons. For the second aspect, data analysis, we analyzed experiments and simulations to tackle three major aspects: network dynamics reconstruction using linear descriptions, dynamics reconstruction using nonlinear descriptors, and the assessment of neuronal connectivity from solely activity data. For the linear study, we analyzed all experiments using the power spectrum density (PSD), and observed that it was sufficiently good to describe the development of the network or its disintegration. PSD also allowed us to distinguish between healthy and unhealthy networks, and revealed dynamical heterogeneities across the network. For the nonlinear study, we used techniques in the context of recurrence plots. We first characterized the embedding dimension m and the time delay δ for each experiment, built the respective recurrence plots, and extracted key information of the dynamics of the system through different descriptors. Experimental results were contrasted with numerical simulations. After analyzing about 400 time series, we concluded that the degree of dynamical complexity in neuronal cultures changes both during development and disintegration. We also observed that the healthier the culture, the higher its dynamic complexity. Finally, for the reconstruction study, we first used numerical simulations to determine the best measure of ‘statistical interdependence’ among any two neurons, and took Generalized Transfer Entropy. We then analyzed the experimental data. We concluded that young cultures have a weak connectivity that increases along maturation. Aggregation increases average connectivity, and more interesting, also the assortativity, i.e. the tendency of highly connected nodes to connect with other highly connected node. In turn, this assortativity may delineates important aspects of the dynamics of the network. Overall, the results show that spatial arrangement and neuronal dynamics are able to shape a very rich repertoire of dynamical states of varying complexity.[cat] L’habilitat dels teixits neuronals de processar i transmetre informació de forma eficient depèn de les propietats dinàmiques intrínseques de les neurones i de la connectivitat entre elles. La present tesi proposa explorar diferents tècniques experimentals i de simulació per analitzar la dinàmica i connectivitat de xarxes neuronals corticals de rata embrionària. Experimentalment, la gravació de l’activitat espontània d’una població de neurones en cultiu, mitjançant una càmera ràpida i tècniques de fluorescència, possibilita el seguiment de forma controlada de l’activitat individual de cada neurona, així com la modificació de la seva connectivitat. En conjunt, aquestes eines permeten estudiar el comportament col.lectiu emergent de la població neuronal. Amb l’objectiu de simular els patrons observats en el laboratori, hem implementat un model mètric aleatori de creixement neuronal per simular la xarxa física de connexions entre neurones, i un model quadràtic d’integració i dispar amb adaptació i depressió per modelar l’ampli espectre de dinàmiques neuronals amb un cost computacional reduït. Hem caracteritzat la dinàmica global i individual de les neurones i l’hem correlacionat amb la seva estructura subjacent mitjançant tècniques lineals i no–lineals de series temporals. L’anàlisi espectral ens ha possibilitat la descripció del desenvolupament i els canvis en connectivitat en els cultius, així com la diferenciació entre cultius sans dels patològics. La reconstrucció de la dinàmica subjacent mitjançant mètodes d’incrustació i l’ús de gràfics de recurrència ens ha permès detectar diferents transicions dinàmiques amb el corresponent guany o pèrdua de la complexitat i riquesa dinàmica del cultiu durant els diferents estudis experimentals. Finalment, a fi de reconstruir la connectivitat interna hem testejat, mitjançant simulacions, diferents quantificadors per mesurar la dependència estadística entre neurona i neurona, seleccionant finalment el mètode de transferència d’entropia gereralitzada. Seguidament, hem procedit a caracteritzar les xarxes amb diferents paràmetres. Malgrat presentar certs tres de xarxes tipus ‘petit món’, els nostres cultius mostren una distribució de grau ‘exponencial’ o ‘esbiaixada’ per, respectivament, cultius joves i madurs. Addicionalment, hem observat que les xarxes homogènies presenten la propietat de disassortativitat, mentre que xarxes amb un creixent nivell d’agregació espaial presenten assortativitat. Aquesta propietat impacta fortament en la transmissió, resistència i sincronització de la xarxa

Dynamical effects of degree correlations in networks of type I model neurons : a dissertation presented in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Mathematics at Massey University, Auckland, New Zealand

The complex behaviour of human brains arises from the complex interconnection of the well-known building blocks -- neurons. With novel imaging techniques it is possible to monitor firing patterns and link them to brain function or dysfunction. How the network structure affects neuronal activity is, however, poorly understood. In this thesis we study the effects of degree correlations in recurrent neuronal networks on self-sustained activity patterns. Firstly, we focus on correlations between the in- and out-degrees of individual neurons. By using Theta Neurons and Ott/Antonsen theory, we can derive a set of coupled differential equations for the expected dynamics of neurons with equal in-degree. A Gaussian copula is used to introduce correlations between a neuron’s in- and out-degree, and numerical bifurcation analysis is used determine the effects of these correlations on the network's dynamics. We find that positive correlations increase the mean firing rate, while negative correlations have the opposite effect. Secondly, we turn to degree correlations between neurons -- often referred to as degree assortativity -- which describes the increased or decreased probability of connecting two neurons based on their in-or out-degrees, relative to what would be expected by chance. We present an alternative derivation of coarse-grained degree mean field equations utilising Theta Neurons and the Ott/Antonsen ansatz as well, but incorporate actual adjacency matrices. Families of degree connectivity matrices are parametrised by assortativity coefficients and subsequently reduced by singular value decomposition. Thus, we efficiently perform numerical bifurcation analysis on a set of coarse-grained equations. To our best knowledge, this is the first time a study examines the four possible types of degree assortativity separately, showing that two have no effect on the networks' dynamics, while the other two can have a significant effect

Degree correlations in directed scale-free networks

Scale-free networks, in which the distribution of the degrees obeys a power-law, are ubiquitous in the study of complex systems. One basic network property that relates to the structure of the links found is the degree assortativity, which is a measure of the correlation between the degrees of the nodes at the end of the links. Degree correlations are known to affect both the structure of a network and the dynamics of the processes supported thereon, including the resilience to damage, the spread of information and epidemics, and the efficiency of defence mechanisms. Nonetheless, while many studies focus on undirected scale-free networks, the interactions in real-world systems often have a directionality. Here, we investigate the dependence of the degree correlations on the power-law exponents in directed scale-free networks. To perform our study, we consider the problem of building directed networks with a prescribed degree distribution, providing a method for proper generation of power-law-distributed directed degree sequences. Applying this new method, we perform extensive numerical simulations, generating ensembles of directed scale-free networks with exponents between~2 and~3, and measuring ensemble averages of the Pearson correlation coefficients. Our results show that scale-free networks are on average uncorrelated across directed links for three of the four possible degree-degree correlations, namely in-degree to in-degree, in-degree to out-degree, and out-degree to out-degree. However, they exhibit anticorrelation between the number of outgoing connections and the number of incoming ones. The findings are consistent with an entropic origin for the observed disassortativity in biological and technological networks.Comment: 10 pages, 5 figure

Front Cell Neurosci

BackgroundSmall clusters comprising hundreds to thousands of neurons are an important level of brain architecture that correlates single neuronal properties to fulfill brain function, but the specific mechanisms through which this scaling occurs are not well understood. In this study, we developed an in vitro experimental platform of small neuronal circuits (islands) to probe the importance of structural properties for their development, physiology, and response to microtrauma.MethodsPrimary cortical neurons were plated on a substrate patterned to promote attachment in clusters of hundreds of cells (islands), transduced with GCaMP6f, allowed to mature until 10\u201313 days in vitro (DIV), and monitored with Ca2+ as a non-invasive proxy for electrical activity. We adjusted two structural factors\u2013island size and cellular density\u2013to evaluate their role in guiding spontaneous activity and network formation in neuronal islands.ResultsWe found cellular density, but not island size, regulates of circuit activity and network function in this system. Low cellular density islands can achieve many states of activity, while high cellular density biases islands towards a limited regime characterized by low rates of activity and high synchronization, a property we summarized as \u201cflexibility.\u201d The injury severity required for an island to lose activity in 50% of its population was significantly higher in low-density, high flexibility islands.ConclusionTogether, these studies demonstrate flexible living cortical circuits are more resilient to microtrauma, providing the first evidence that initial circuit state may be a key factor to consider when evaluating the consequences of trauma to the cortex

Structure-function relation in a stochastic whole-brain model at criticality

Understanding the relation between brain architecture and function is one of the central issues in neuroscience nowadays. In the last few years, important efforts have been devoted to map the large-scale structure of the human cortex, the so-called "connectome". An example is the neuroanatomical connectivity matrix of the entire human brain obtained through MR diffusion tractography. Recent studies proposed a stochastic model built on top of this connectivity matrix that displays a phase-transition and is able to reproduce several aspects of brain functioning when tuned to its critical point. This master thesis is aimed to review recent results on this subject and to get a deeper insight into the model by studying the distribution of the avalanches, the dynamical range and to investigate how the use of simulated connectivity matrices affects the dynamics. Furthermore, a theoretical description of the dynamics is proposed by introducing a master equation in order to understand the nature of the phase transition and the role of stochasticity.ope