Dynamical networks: Finding, measuring, and tracking neural population activity using network science

MEA-NAP: A flexible network analysis pipeline for neuronal 2D and 3D organoid multielectrode recordings

Microelectrode array (MEA) recordings are commonly used to compare firing and burst rates in neuronal cultures. MEA recordings can also reveal microscale functional connectivity, topology, and network dynamics-patterns seen in brain networks across spatial scales. Network topology is frequently characterized in neuroimaging with graph theoretical metrics. However, few computational tools exist for analyzing microscale functional brain networks from MEA recordings. Here, we present a MATLAB MEA network analysis pipeline (MEA-NAP) for raw voltage time series acquired from single- or multi-well MEAs. Applications to 3D human cerebral organoids or 2D human-derived or murine cultures reveal differences in network development, including topology, node cartography, and dimensionality. MEA-NAP incorporates multi-unit template-based spike detection, probabilistic thresholding for determining significant functional connections, and normalization techniques for comparing networks. MEA-NAP can identify network-level effects of pharmacologic perturbation and/or disease-causing mutations and thus can provide a translational platform for revealing mechanistic insights and screening new therapeutic approaches. VIDEO ABSTRACT.

Stimulus type shapes the topology of cellular functional networks in mouse visual cortex

On the timescale of sensory processing, neuronal networks have relatively fixed anatomical connectivity, while functional interactions between neurons can vary depending on the ongoing activity of the neurons within the network. We thus hypothesized that different types of stimuli could lead those networks to display stimulus-dependent functional connectivity patterns. To test this hypothesis, we analyzed single-cell resolution electrophysiological data from the Allen Institute, with simultaneous recordings of stimulus-evoked activity from neurons across 6 different regions of mouse visual cortex. Comparing the functional connectivity patterns during different stimulus types, we made several nontrivial observations: (1) while the frequencies of different functional motifs were preserved across stimuli, the identities of the neurons within those motifs changed; (2) the degree to which functional modules are contained within a single brain region increases with stimulus complexity. Altogether, our work reveals unexpected stimulus-dependence to the way groups of neurons interact to process incoming sensory information.

MEA-NAP compares microscale functional connectivity, topology, and network dynamics in organoid or monolayer neuronal cultures

Microelectrode array (MEA) recordings are commonly used to compare firing and burst rates in neuronal cultures. MEA recordings can also reveal microscale functional connectivity, topology, and network dynamics-patterns seen in brain networks across spatial scales. Network topology is frequently characterized in neuroimaging with graph theoretical metrics. However, few computational tools exist for analyzing microscale functional brain networks from MEA recordings. Here, we present a MATLAB MEA network analysis pipeline (MEA-NAP) for raw voltage time-series acquired from single- or multi-well MEAs. Applications to 3D human cerebral organoids or 2D human-derived or murine cultures reveal differences in network development, including topology, node cartography, and dimensionality. MEA-NAP incorporates multi-unit template-based spike detection, probabilistic thresholding for determining significant functional connections, and normalization techniques for comparing networks. MEA-NAP can identify network-level effects of pharmacologic perturbation and/or disease-causing mutations and, thus, can provide a translational platform for revealing mechanistic insights and screening new therapeutic approaches.

The Structure of Hippocampal CA1 Interactions Optimizes Spatial Coding across Experience

Although much is known about how single neurons in the hippocampus represent an animal's position, how circuit interactions contribute to spatial coding is less well understood. Using a novel statistical estimator and theoretical modeling, both developed in the framework of maximum entropy models, we reveal highly structured CA1 cell-cell interactions in male rats during open field exploration. The statistics of these interactions depend on whether the animal is in a familiar or novel environment. In both conditions the circuit interactions optimize the encoding of spatial information, but for regimes that differ in the informativeness of their spatial inputs. This structure facilitates linear decodability, making the information easy to read out by downstream circuits. Overall, our findings suggest that the efficient coding hypothesis is not only applicable to individual neuron properties in the sensory periphery, but also to neural interactions in the central brain.SIGNIFICANCE STATEMENT Local circuit interactions play a key role in neural computation and are dynamically shaped by experience. However, measuring and assessing their effects during behavior remains a challenge. Here, we combine techniques from statistical physics and machine learning to develop new tools for determining the effects of local network interactions on neural population activity. This approach reveals highly structured local interactions between hippocampal neurons, which make the neural code more precise and easier to read out by downstream circuits, across different levels of experience. More generally, the novel combination of theory and data analysis in the framework of maximum entropy models enables traditional neural coding questions to be asked in naturalistic settings.

Structure and function in artificial, zebrafish and human neural networks

Network science provides a set of tools for the characterization of the structure and functional behavior of complex systems. Yet a major problem is to quantify how the structural domain is related to the dynamical one. In other words, how the diversity of dynamical states of a system can be predicted from the static network structure? Or the reverse problem: starting from a set of signals derived from experimental recordings, how can one discover the network connections or the causal relations behind the observed dynamics? Despite the advances achieved over the last two decades, many challenges remain concerning the study of the structure-dynamics interplay of complex systems. In neuroscience, progress is typically constrained by the low spatio-temporal resolution of experiments and by the lack of a universal inferring framework for empirical systems. To address these issues, applications of network science and artificial intelligence to neural data have been rapidly growing. In this article, we review important recent applications of methods from those fields to the study of the interplay between structure and functional dynamics of human and zebrafish brain. We cover the selection of topological features for the characterization of brain networks, inference of functional connections, dynamical modeling, and close with applications to both the human and zebrafish brain. This review is intended to neuroscientists who want to become acquainted with techniques from network science, as well as to researchers from the latter field who are interested in exploring novel application scenarios in neuroscience.

Identifying Steady State in the Network Dynamics of Spiking Neural Networks

Analysis of the dynamics of complex networks can provide valuable information. For example, the dynamics can be used to characterize and differentiate between different network inputs and configurations. However, without quantitatively delineating the network's dynamic regimes, analysis of the network's dynamics is based on heuristics and qualitative signatures of transient or steady-state regimes. This is not ideal because interesting phenomena can occur during the transient regime, steady-state regime, or at the transition between the two dynamic regimes. Moreover, for simulated and observed systems, precise knowledge of the network's dynamical regime is imperative when considering metrics on minimal mathematical descriptions of the dynamics, otherwise either too much or too little data is analyzed. Here, we develop quantitative methods to ascertain the starting point and period of steady-state network activity. Using the precise knowledge of the network's dynamic regimes, we build minimal representations of the network dynamics that form the basis for future work. We show applications of our techniques on idealized signals and on the dynamics of a biologically inspired spiking neural network.

Scalable and accurate method for neuronal ensemble detection in spiking neural networks

We propose a novel, scalable, and accurate method for detecting neuronal ensembles from a population of spiking neurons. Our approach offers a simple yet powerful tool to study ensemble activity. It relies on clustering synchronous population activity (population vectors), allows the participation of neurons in different ensembles, has few parameters to tune and is computationally efficient. To validate the performance and generality of our method, we generated synthetic data, where we found that our method accurately detects neuronal ensembles for a wide range of simulation parameters. We found that our method outperforms current alternative methodologies. We used spike trains of retinal ganglion cells obtained from multi-electrode array recordings under a simple ON-OFF light stimulus to test our method. We found a consistent stimuli-evoked ensemble activity intermingled with spontaneously active ensembles and irregular activity. Our results suggest that the early visual system activity could be organized in distinguishable functional ensembles. We provide a Graphic User Interface, which facilitates the use of our method by the scientific community.

Integrating new memories into the hippocampal network activity space

By investigating the topology of neuronal co-activity, we found that mnemonic information spans multiple operational axes in the mouse hippocampus network. High-activity principal cells form the core of each memory along a first axis, segregating spatial contexts and novelty. Low-activity cells join co-activity motifs across behavioral events and enable their crosstalk along two other axes. This reveals an organizational principle for continuous integration and interaction of hippocampal memories.

Nonlinear response characteristics of neural networks and single neurons undergoing optogenetic excitation

Optogenetic stimulation has become the method of choice for investigating neural computation in populations of neurons. Optogenetic experiments often aim to elicit a network response by stimulating specific groups of neurons. However, this is complicated by the fact that optogenetic stimulation is nonlinear, more light does not always equal to more spikes, and neurons that are not directly but indirectly stimulated could have a major impact on how networks respond to optogenetic stimulation. To clarify how optogenetic excitation of some neurons alters the network dynamics, we studied the temporal and spatial response of individual neurons and recurrent neural networks. In individual neurons, we find that neurons show a monotonic, saturating rate response to increasing light intensity and a nonmonotonic rate response to increasing pulse frequency. At the network level, we find that Gaussian light beams elicit spatial firing rate responses that are substantially broader than the stimulus profile. In summary, our analysis and our network simulation code allow us to predict the outcome of an optogenetic experiment and to assess whether the observed effects can be attributed to direct or indirect stimulation of neurons.

Optogenetic circuit manipulation has become a popular tool to manipulate the activity of neurons. During optogenetic stimulation, the firing rate of a neuron can rise because of direct light excitation or indirect activation via other neurons. To disentangle these influences and predict the effects of optogenetic excitation, we set up a spiking network model with controlled connectivity and studied its response to light stimulation. We find that the optogenetically evoked activity in a network can spread far beyond the light-stimulated area. We further found a nonmonotonic rate response of single neurons to increasing light intensities and frequencies. Our results help to interpret optogenetic experiments in vivo, and we provide computer code that can be customized to simulate 2D connectivity scenarios and explore their consequences.

Extracting single-trial neural interaction using latent dynamical systems model

In systems neuroscience, advances in simultaneous recording technology have helped reveal the population dynamics that underlie the complex neural correlates of animal behavior and cognitive processes. To investigate these correlates, neural interactions are typically abstracted from spike trains of pairs of neurons accumulated over the course of many trials. However, the resultant averaged values do not lead to understanding of neural computation in which the responses of populations are highly variable even under identical external conditions. Accordingly, neural interactions within the population also show strong fluctuations. In the present study, we introduce an analysis method reflecting the temporal variation of neural interactions, in which cross-correlograms on rate estimates are applied via a latent dynamical systems model. Using this method, we were able to predict time-varying neural interactions within a single trial. In addition, the pairwise connections estimated in our analysis increased along behavioral epochs among neurons categorized within similar functional groups. Thus, our analysis method revealed that neurons in the same groups communicate more as the population gets involved in the assigned task. We also showed that the characteristics of neural interaction from our model differ from the results of a typical model employing cross-correlation coefficients. This suggests that our model can extract nonoverlapping information about network topology, unlike the typical model.

Organizing principles of whole-brain functional connectivity in zebrafish larvae

Network science has begun to reveal the fundamental principles by which large-scale brain networks are organized, including geometric constraints, a balance between segregative and integrative features, and functionally flexible brain areas. However, it remains unknown whether whole-brain networks imaged at the cellular level are organized according to similar principles. Here, we analyze whole-brain functional networks reconstructed from calcium imaging data recorded in larval zebrafish. Our analyses reveal that functional connections are distance-dependent and that networks exhibit hierarchical modular structure and hubs that span module boundaries. We go on to show that spontaneous network structure places constraints on stimulus-evoked reconfigurations of connections and that networks are highly consistent across individuals. Our analyses reveal basic organizing principles of whole-brain functional brain networks at the mesoscale. Our overarching methodological framework provides a blueprint for studying correlated activity at the cellular level using a low-dimensional network representation. Our work forms a conceptual bridge between macro- and mesoscale network neuroscience and opens myriad paths for future studies to investigate network structure of nervous systems at the cellular level.

Stability of spontaneous, correlated activity in mouse auditory cortex

Neural systems can be modeled as complex networks in which neural elements are represented as nodes linked to one another through structural or functional connections. The resulting network can be analyzed using mathematical tools from network science and graph theory to quantify the system’s topological organization and to better understand its function. Here, we used two-photon calcium imaging to record spontaneous activity from the same set of cells in mouse auditory cortex over the course of several weeks. We reconstruct functional networks in which cells are linked to one another by edges weighted according to the correlation of their fluorescence traces. We show that the networks exhibit modular structure across multiple topological scales and that these multi-scale modules unfold as part of a hierarchy. We also show that, on average, network architecture becomes increasingly dissimilar over time, with similarity decaying monotonically with the distance (in time) between sessions. Finally, we show that a small fraction of cells maintain strongly-correlated activity over multiple days, forming a stable temporal core surrounded by a fluctuating and variable periphery. Our work indicates a framework for studying spontaneous activity measured by two-photon calcium imaging using computational methods and graphical models from network science. The methods are flexible and easily extended to additional datasets, opening the possibility of studying cellular level network organization of neural systems and how that organization is modulated by stimuli or altered in models of disease.

Neurons coordinate their activity with one another, forming networks that help support adaptive, flexible behavior. Still, little is known about the organization of these networks at the cellular scale and their stability over time. Here, we reconstruct networks from calcium imaging data recorded in mouse primary auditory cortex. We show that these networks exhibit spatially constrained, hierarchical modular structure, which may facilitate specialized information processing. However, we show that connection weights and modular structure are also variable over time, changing on a timescale of days and adopting novel network configurations. Despite this, a small subset of neurons maintain their connections to one another and preserve their modular organization across time, forming a stable temporal core surrounded by a flexible periphery. These findings represent a conceptual bridge linking network analyses of macroscale and cellular-level neuroimaging data. They also represent a complementary approach to existing circuits- and systems-based interrogation of nervous system function, opening the door for deeper and more targeted analysis in the future.

Brain network segregation and integration during an epoch-related working memory fMRI experiment

The characterization of brain subnetwork segregation and integration has previously focused on changes that are detectable at the level of entire sessions or epochs of imaging data. In this study, we applied time-varying functional connectivity analysis together with temporal network theory to calculate point-by-point estimates in subnetwork segregation and integration during an epoch-based (2-back, 0-back, baseline) working memory fMRI experiment as well as during resting-state. This approach allowed us to follow task-related changes in subnetwork segregation and integration at a high temporal resolution. At a global level, the cognitively more taxing 2-back epochs elicited an overall stronger response of integration between subnetworks compared to the 0-back epochs. Moreover, the visual, sensorimotor and fronto-parietal subnetworks displayed characteristic and distinct temporal profiles of segregation and integration during the 0- and 2-back epochs. During the interspersed epochs of baseline, several subnetworks, including the visual, fronto-parietal, cingulo-opercular and dorsal attention subnetworks showed pronounced increases in segregation. Using a drift diffusion model we show that the response time for the 2-back trials are correlated with integration for the fronto-parietal subnetwork and correlated with segregation for the visual subnetwork. Our results elucidate the fast-evolving events with regard to subnetwork integration and segregation that occur in an epoch-related task fMRI experiment. Our findings suggest that minute changes in subnetwork integration are of importance for task performance.