A systematic assessment of global neural network connectivity through direct electrophysiological assays has remained technically infeasible, even in simpler systems like dissociated neuronal cultures. We introduce an improved algorithmic approach based on Transfer Entropy to reconstruct structural connectivity from network activity monitored through calcium imaging. We focus in this study on the inference of excitatory synaptic links. Based on information theory, our method requires no prior assumptions on the statistics of neuronal firing and neuronal connections. The performance of our algorithm is benchmarked on surrogate time series of calcium fluorescence generated by the simulated dynamics of a network with known ground-truth topology. We find that the functional network topology revealed by Transfer Entropy depends qualitatively on the time-dependent dynamic state of the network (bursting or non-bursting). Thus by conditioning with respect to the global mean activity, we improve the performance of our method. This allows us to focus the analysis to specific dynamical regimes of the network in which the inferred functional connectivity is shaped by monosynaptic excitatory connections, rather than by collective synchrony. Our method can discriminate between actual causal influences between neurons and spurious non-causal correlations due to light scattering artifacts, which inherently affect the quality of fluorescence imaging. Compared to other reconstruction strategies such as cross-correlation or Granger Causality methods, our method based on improved Transfer Entropy is remarkably more accurate. In particular, it provides a good estimation of the excitatory network clustering coefficient, allowing for discrimination between weakly and strongly clustered topologies. Finally, we demonstrate the applicability of our method to analyses of real recordings of in vitro disinhibited cortical cultures where we suggest that excitatory connections are characterized by an elevated level of clustering compared to a random graph (although not extreme) and can be markedly non-local. Unraveling the general organizing principles of connectivity in neural circuits is a crucial step towards understanding brain function. However, even the simpler task of assessing the global excitatory connectivity of a culture in vitro, where neurons form self-organized networks in absence of external stimuli, remains challenging. Neuronal cultures undergo spontaneous switching between episodes of synchronous bursting and quieter inter-burst periods. We introduce here a novel algorithm which aims at inferring the connectivity of neuronal cultures from calcium fluorescence recordings of their network dynamics. To achieve this goal, we develop a suitable generalization of Transfer Entropy, an information-theoretic measure of causal influences between time series. Unlike previous algorithmic approaches to reconstruction, Transfer Entropy is data-driven and does not rely on specific assumptions about neuronal firing statistics or network topology. We generate simulated calcium signals from networks with controlled ground-truth topology and purely excitatory interactions and show that, by restricting the analysis to inter-bursts periods, Transfer Entropy robustly achieves a good reconstruction performance for disparate network connectivities. Finally, we apply our method to real data and find evidence of non-random features in cultured networks, such as the existence of highly connected hub excitatory neurons and of an elevated (but not extreme) level of clustering.
We introduce an approach for the quantitative assessment of the connectivity in neuronal cultures, based on the statistical mechanics of percolation on a graph. This allows us to monitor the development of the culture and to see the emergence of connectivity in the network. The culture becomes fully connected at a time equivalent to the expected time of birth. The spontaneous bursting activity that characterizes cultures develops in parallel with the connectivity. The average number of inputs per neuron can be quantitatively determined in units of m0 , the number of activated inputs needed to excite the neuron. For m0 ≃ 15 we find that hippocampal neurons have on average ≈60–120 inputs, whereas cortical neurons have ≈75–150, depending on neuronal density. The ratio of excitatory to inhibitory neurons is determined by using the GABA A antagonist bicuculine. This ratio changes during development and reaches the final value at day 7–8, coinciding with the expected time of the GABA switch. For hippocampal cultures the inhibitory cells comprise ≈30% of the neurons in the culture whereas for cortical cultures they are ≈20%. Such detailed global information on the connectivity of networks in neuronal cultures is at present inaccessible by any electrophysiological or other technique.
The epigenomic landscape of Parkinson's disease (PD) remains unknown. We performed a genomewide DNA methylation and a transcriptome studies in induced pluripotent stem cell (iPSC)‐derived dopaminergic neurons (DAn) generated by cell reprogramming of somatic skin cells from patients with monogenic LRRK2‐associated PD (L2PD) or sporadic PD (sPD), and healthy subjects. We observed extensive DNA methylation changes in PD DAn, and of RNA expression, which were common in L2PD and sPD. No significant methylation differences were present in parental skin cells, undifferentiated iPSCs nor iPSC‐derived neural cultures not‐enriched‐in‐DAn. These findings suggest the presence of molecular defects in PD somatic cells which manifest only upon differentiation into the DAn cells targeted in PD. The methylation profile from PD DAn, but not from controls, resembled that of neural cultures not‐enriched‐in‐DAn indicating a failure to fully acquire the epigenetic identity own to healthy DAn in PD. The PD‐associated hypermethylation was prominent in gene regulatory regions such as enhancers and was related to the RNA and/or protein downregulation of a network of transcription factors relevant to PD (FOXA1, NR3C1, HNF4A, and FOSL2). Using a patient‐specific iPSC‐based DAn model, our study provides the first evidence that epigenetic deregulation is associated with monogenic and sporadic PD.
We study living neural networks by measuring the neurons’ response to a global electrical stimulation. Neural connectivity is lowered by reducing the synaptic strength, chemically blocking neurotransmitter receptors. We use a graph-theoretic approach to show that the connectivity undergoes a percolation transition. This occurs as the giant component disintegrates, characterized by a power law with an exponent . is independent of the balance between excitatory and inhibitory neurons and indicates that the degree distribution is Gaussian rather than scale free.