Emergence of rich-club topology and coordinated dynamics in development of hippocampal functional networks in vitro

Repetitive training enhances the pattern recognition capability of cultured neural networks

Cultured neural networks in vitro have demonstrated the biocomputing capability to recognize patterns. However, the underlying mechanisms behind information processing and pattern recognition remain less understood. Here, we developed an in vitro neural network integrated with microelectrode arrays (MEAs) to explore the network's classification capability and elucidate the mechanisms underlying this classification. After applying different stimulation patterns using MEAs, the network exhibited structural alterations and distinct electrical responses that recognized various stimulation patterns. Alongside the reshaping of network structures, repeated training increased recognition accuracy for each stimulation pattern. Additionally, it was reported for the first time that spontaneous networks after stimulation are more closely related to the structures of evoked networks. This work provides new insights into the structural changes underlying information processing and contributes to our understanding of how cultured neural networks respond to different patterns.

Neural network embedding of functional microconnectome

Our brains operate as a complex network of interconnected neurons. To gain a deeper understanding of this network architecture, it is essential to extract simple rules from its intricate structure. This study aimed to compress and simplify the architecture, with a particular focus on interpreting patterns of functional connectivity in 2.5 hr of electrical activity from a vast number of neurons in acutely sliced mouse brains. Here, we combined two distinct methods together: automatic compression and network analysis. Firstly, for automatic compression, we trained an artificial neural network named NNE (neural network embedding). This allowed us to reduce the connectivity to features, be represented only by 13% of the original neuron count. Secondly, to decipher the topology, we concentrated on the variability among the compressed features and compared them with 15 distinct network metrics. Specifically, we introduced new metrics that had not previously existed, termed as indirect-adjacent degree and neighbor hub ratio. Our results conclusively demonstrated that these new metrics could better explain approximately 40%-45% of the features. This finding highlighted the critical role of NNE in facilitating the development of innovative metrics, because some of the features extracted by NNE were not captured by the currently existed network metrics.

Brain Connectivity, Neural Networks, and Resilience in Aging and Neurodegeneration

The importance of complex systems has become increasingly evident in recent years. The nervous system is one such example, with neural networks sitting at the intersection of complex networks and biology. A particularly exciting feature is the resilience of complex systems. For example, the ability of the nervous system to perform even in the face of challenges that include neuronal loss, neuroinflammation, protein accumulation, axonal disruptions, and metabolic stress is an intriguing and exciting line of investigation. In neurodegenerative diseases, neural network resilience is responsible for the time between the earliest disease-linked changes and clinical symptom onset and disease diagnosis. In this way, connectivity resilience of neurons within the complex network of cells that make up the nervous system has significant implications. This review provides an overview of relevant concepts related to complex systems with a focus on the connectivity of the nervous system. It discusses the development of the neural network and how a delicate balance determines how this complex system responds to injury, with examples illustrating maladaptive plasticity. The review then addresses the implications of these concepts, methods to understand brain connectivity and neural networks, and recent research efforts aimed at understanding neurodegeneration from this perspective. This study aims to provide foundational knowledge and an overview of current research directions in this evolving and exciting area of neuroscience.

Protosequences in brain organoids model intrinsic brain states Authors

Neuronal firing sequences are thought to be the basic building blocks of neural coding and information broadcasting within the brain. However, when sequences emerge during neurodevelopment remains unknown. We demonstrate that structured firing sequences are present in spontaneous activity of human and murine brain organoids and ex vivo neonatal brain slices from the murine somatosensory cortex. We observed a balance between temporally rigid and flexible firing patterns that are emergent phenomena in human and murine brain organoids and early postnatal murine somatosensory cortex, but not in primary dissociated cortical cultures. Our findings suggest that temporal sequences do not arise in an experience-dependent manner, but are rather constrained by an innate preconfigured architecture established during neurogenesis. These findings highlight the potential for brain organoids to further explore how exogenous inputs can be used to refine neuronal circuits and enable new studies into the genetic mechanisms that govern assembly of functional circuitry during early human brain development.

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.

Response inhibition in premotor cortex corresponds to a complex reshuffle of the mesoscopic information network

Recent studies have explored functional and effective neural networks in animal models; however, the dynamics of information propagation among functional modules under cognitive control remain largely unknown. Here, we addressed the issue using transfer entropy and graph theory methods on mesoscopic neural activities recorded in the dorsal premotor cortex of rhesus monkeys. We focused our study on the decision time of a Stop-signal task, looking for patterns in the network configuration that could influence motor plan maturation when the Stop signal is provided. When comparing trials with successful inhibition to those with generated movement, the nodes of the network resulted organized into four clusters, hierarchically arranged, and distinctly involved in information transfer. Interestingly, the hierarchies and the strength of information transmission between clusters varied throughout the task, distinguishing between generated movements and canceled ones and corresponding to measurable levels of network complexity. Our results suggest a putative mechanism for motor inhibition in premotor cortex: a topological reshuffle of the information exchanged among ensembles of neurons.

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.

Ensemble learning and ground-truth validation of synaptic connectivity inferred from spike trains

Probing the architecture of neuronal circuits and the principles that underlie their functional organization remains an important challenge of modern neurosciences. This holds true, in particular, for the inference of neuronal connectivity from large-scale extracellular recordings. Despite the popularity of this approach and a number of elaborate methods to reconstruct networks, the degree to which synaptic connections can be reconstructed from spike-train recordings alone remains controversial. Here, we provide a framework to probe and compare connectivity inference algorithms, using a combination of synthetic ground-truth and in vitro data sets, where the connectivity labels were obtained from simultaneous high-density microelectrode array (HD-MEA) and patch-clamp recordings. We find that reconstruction performance critically depends on the regularity of the recorded spontaneous activity, i.e., their dynamical regime, the type of connectivity, and the amount of available spike-train data. We therefore introduce an ensemble artificial neural network (eANN) to improve connectivity inference. We train the eANN on the validated outputs of six established inference algorithms and show how it improves network reconstruction accuracy and robustness. Overall, the eANN demonstrated strong performance across different dynamical regimes, worked well on smaller datasets, and improved the detection of synaptic connectivity, especially inhibitory connections. Results indicated that the eANN also improved the topological characterization of neuronal networks. The presented methodology contributes to advancing the performance of inference algorithms and facilitates our understanding of how neuronal activity relates to synaptic connectivity.

Cortical timescales and the modular organization of structural and functional brain networks

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

The neuropeptidergic connectome of C. elegans

Efforts are ongoing to map synaptic wiring diagrams, or connectomes, to understand the neural basis of brain function. However, chemical synapses represent only one type of functionally important neuronal connection; in particular, extrasynaptic, "wireless" signaling by neuropeptides is widespread and plays essential roles in all nervous systems. By integrating single-cell anatomical and gene-expression datasets with biochemical analysis of receptor-ligand interactions, we have generated a draft connectome of neuropeptide signaling in the C. elegans nervous system. This network is characterized by high connection density, extended signaling cascades, autocrine foci, and a decentralized topology, with a large, highly interconnected core containing three constituent communities sharing similar patterns of input connectivity. Intriguingly, several key network hubs are little-studied neurons that appear specialized for peptidergic neuromodulation. We anticipate that the C. elegans neuropeptidergic connectome will serve as a prototype to understand how networks of neuromodulatory signaling are organized.

High-resolution CMOS-based biosensor for assessing hippocampal circuit dynamics in experience-dependent plasticity

Experiential richness creates tissue-level changes and synaptic plasticity as patterns emerge from rhythmic spatiotemporal activity of large interconnected neuronal assemblies. Despite numerous experimental and computational approaches at different scales, the precise impact of experience on network-wide computational dynamics remains inaccessible due to the lack of applicable large-scale recording methodology. We here demonstrate a large-scale multi-site biohybrid brain circuity on-CMOS-based biosensor with an unprecedented spatiotemporal resolution of 4096 microelectrodes, which allows simultaneous electrophysiological assessment across the entire hippocampal-cortical subnetworks from mice living in an enriched environment (ENR) and standard-housed (SD) conditions. Our platform, empowered with various computational analyses, reveals environmental enrichment's impacts on local and global spatiotemporal neural dynamics, firing synchrony, topological network complexity, and large-scale connectome. Our results delineate the distinct role of prior experience in enhancing multiplexed dimensional coding formed by neuronal ensembles and error tolerance and resilience to random failures compared to standard conditions. The scope and depth of these effects highlight the critical role of high-density, large-scale biosensors to provide a new understanding of the computational dynamics and information processing in multimodal physiological and experience-dependent plasticity conditions and their role in higher brain functions. Knowledge of these large-scale dynamics can inspire the development of biologically plausible computational models and computational artificial intelligence networks and expand the reach of neuromorphic brain-inspired computing into new applications.

Cortical timescales and the modular organization of structural and functional brain networks

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

Modeling the three-dimensional connectivity of in vitro cortical ensembles coupled to Micro-Electrode Arrays

Nowadays, in vitro three-dimensional (3D) neuronal networks are becoming a consolidated experimental model to overcome most of the intrinsic limitations of bi-dimensional (2D) assemblies. In the 3D environment, experimental evidence revealed a wider repertoire of activity patterns, characterized by a modulation of the bursting features, than the one observed in 2D cultures. However, it is not totally clear and understood what pushes the neuronal networks towards different dynamical regimes. One possible explanation could be the underlying connectivity, which could involve a larger number of neurons in a 3D rather than a 2D space and could organize following well-defined topological schemes. Driven by experimental findings, achieved by recording 3D cortical networks organized in multi-layered structures coupled to Micro-Electrode Arrays (MEAs), in the present work we developed a large-scale computational network model made up of leaky integrate-and-fire (LIF) neurons to investigate possible structural configurations able to sustain the emerging patterns of electrophysiological activity. In particular, we investigated the role of the number of layers defining a 3D assembly and the spatial distribution of the connections within and among the layers. These configurations give rise to different patterns of activity that could be compared to the ones emerging from real in vitro 3D neuronal populations. Our results suggest that the introduction of three-dimensionality induced a global reduction in both firing and bursting rates with respect to 2D models. In addition, we found that there is a minimum number of layers necessary to obtain a change in the dynamics of the network. However, the effects produced by a 3D organization of the cells is somewhat mitigated if a scale-free connectivity is implemented in either one or all the layers of the network. Finally, the best matching of the experimental data is achieved supposing a 3D connectivity organized in structured bundles of links located in different areas of the 2D network.

Nanomaterial-based microelectrode arrays for in vitro bidirectional brain-computer interfaces: a review

A bidirectional in vitro brain-computer interface (BCI) directly connects isolated brain cells with the surrounding environment, reads neural signals and inputs modulatory instructions. As a noninvasive BCI, it has clear advantages in understanding and exploiting advanced brain function due to the simplified structure and high controllability of ex vivo neural networks. However, the core of ex vivo BCIs, microelectrode arrays (MEAs), urgently need improvements in the strength of signal detection, precision of neural modulation and biocompatibility. Notably, nanomaterial-based MEAs cater to all the requirements by converging the multilevel neural signals and simultaneously applying stimuli at an excellent spatiotemporal resolution, as well as supporting long-term cultivation of neurons. This is enabled by the advantageous electrochemical characteristics of nanomaterials, such as their active atomic reactivity and outstanding charge conduction efficiency, improving the performance of MEAs. Here, we review the fabrication of nanomaterial-based MEAs applied to bidirectional in vitro BCIs from an interdisciplinary perspective. We also consider the decoding and coding of neural activity through the interface and highlight the various usages of MEAs coupled with the dissociated neural cultures to benefit future developments of BCIs.

Selective inhibition of excitatory synaptic transmission alters the emergent bursting dynamics of in vitro neural networks

Neurons in vitro connect to each other and form neural networks that display emergent electrophysiological activity. This activity begins as spontaneous uncorrelated firing in the early phase of development, and as functional excitatory and inhibitory synapses mature, the activity typically emerges as spontaneous network bursts. Network bursts are events of coordinated global activation among many neurons interspersed with periods of silencing and are important for synaptic plasticity, neural information processing, and network computation. While bursting is the consequence of balanced excitatory-inhibitory (E/I) interactions, the functional mechanisms underlying their evolution from physiological to potentially pathophysiological states, such as decreasing or increasing in synchrony, are still poorly understood. Synaptic activity, especially that related to maturity of E/I synaptic transmission, is known to strongly influence these processes. In this study, we used selective chemogenetic inhibition to target and disrupt excitatory synaptic transmission in in vitro neural networks to study functional response and recovery of spontaneous network bursts over time. We found that over time, inhibition resulted in increases in both network burstiness and synchrony. Our results indicate that the disruption in excitatory synaptic transmission during early network development likely affected inhibitory synaptic maturity which resulted in an overall decrease in network inhibition at later stages. These findings lend support to the importance of E/I balance in maintaining physiological bursting dynamics and, conceivably, information processing capacity in neural networks.

Long-term calcium imaging reveals functional development in hiPSC-derived cultures comparable to human but not rat primary cultures

Models for human brain-oriented research are often established on primary cultures from rodents, which fails to recapitulate cellular specificity and molecular cues of the human brain. Here we investigated whether neuronal cultures derived from human induced pluripotent stem cells (hiPSCs) feature key advantages compared with rodent primary cultures. Using calcium fluorescence imaging, we tracked spontaneous neuronal activity in hiPSC-derived, human, and rat primary cultures and compared their dynamic and functional behavior as they matured. We observed that hiPSC-derived cultures progressively changed upon development, exhibiting gradually richer activity patterns and functional traits. By contrast, rat primary cultures were locked in the same dynamic state since activity onset. Human primary cultures exhibited features in between hiPSC-derived and rat primary cultures, although traits from the former predominated. Our study demonstrates that hiPSC-derived cultures are excellent models to investigate development in neuronal assemblies, a hallmark for applications that monitor alterations caused by damage or neurodegeneration.

Dynamic and Functional Alterations of Neuronal Networks In Vitro upon Physical Damage: A Proof of Concept

There is a growing technological interest in combining biological neuronal networks with electronic ones, specifically for biological computation, human-machine interfacing and robotic implants. A major challenge for the development of these technologies is the resilience of the biological networks to physical damage, for instance, when used in harsh environments. To tackle this question, here, we investigated the dynamic and functional alterations of rodent cortical networks grown in vitro that were physically damaged, either by sequentially removing groups of neurons that were central for information flow or by applying an incision that cut the network in half. In both cases, we observed a remarkable capacity of the neuronal cultures to cope with damage, maintaining their activity and even reestablishing lost communication pathways. We also observed-particularly for the cultures cut in half-that a reservoir of healthy neurons surrounding the damaged region could boost resilience by providing stimulation and a communication bridge across disconnected areas. Our results show the remarkable capacity of neuronal cultures to sustain and recover from damage, and may be inspirational for the development of future hybrid biological-electronic systems.

Flexibility of in vitro cortical circuits influences resilience from microtrauma

Background

Small 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.

Methods

Primary 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-13 days in vitro (DIV), and monitored with Ca²⁺ as a non-invasive proxy for electrical activity. We adjusted two structural factors-island size and cellular density-to evaluate their role in guiding spontaneous activity and network formation in neuronal islands.

Results

We 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 "flexibility." The injury severity required for an island to lose activity in 50% of its population was significantly higher in low-density, high flexibility islands.

Conclusion

Together, 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.

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.

Learning populations with hubs govern the initiation and propagation of spontaneous bursts in neuronal networks after learning

Spontaneous bursts in neuronal networks with propagation involving a large number of synchronously firing neurons are considered to be a crucial feature of these networks both in vivo and in vitro. Recently, learning has been shown to improve the association and synchronization of spontaneous events in neuronal networks by promoting the firing of spontaneous bursts. However, little is known about the relationship between the learning phase and spontaneous bursts. By combining high-resolution measurement with a 4,096-channel complementary metal-oxide-semiconductor (CMOS) microelectrode array (MEA) and graph theory, we studied how the learning phase influenced the initiation of spontaneous bursts in cultured networks of rat cortical neurons in vitro. We found that a small number of selected populations carried most of the stimulus information and contributed to learning. Moreover, several new burst propagation patterns appeared in spontaneous firing after learning. Importantly, these "learning populations" had more hubs in the functional network that governed the initiation of spontaneous burst activity. These results suggest that changes in the functional structure of learning populations may be the key mechanism underlying increased bursts after learning. Our findings could increase understanding of the important role that synaptic plasticity plays in the regulation of spontaneous activity.

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.

Stimulus-Evoked Activity Modulation of In Vitro Engineered Cortical and Hippocampal Networks

The delivery of electrical stimuli is crucial to shape the electrophysiological activity of neuronal populations and to appreciate the response of the different brain circuits involved. In the present work, we used dissociated cortical and hippocampal networks coupled to Micro-Electrode Arrays (MEAs) to investigate the features of their evoked response when a low-frequency (0.2 Hz) electrical stimulation protocol is delivered. In particular, cortical and hippocampal neurons were topologically organized to recreate interconnected sub-populations with a polydimethylsiloxane (PDMS) mask, which guaranteed the segregation of the cell bodies and the connections among the sub-regions through microchannels. We found that cortical assemblies were more reactive than hippocampal ones. Despite both configurations exhibiting a fast (<35 ms) response, this did not uniformly distribute over the MEA in the hippocampal networks. Moreover, the propagation of the stimuli-evoked activity within the networks showed a late (35−500 ms) response only in the cortical assemblies. The achieved results suggest the importance of the neuronal target when electrical stimulation experiments are performed. Not all neuronal types display the same response, and in light of transferring stimulation protocols to in vivo applications, it becomes fundamental to design realistic in vitro brain-on-a-chip devices to investigate the dynamical properties of complex neuronal circuits.

Automatic differentiation of human induced pluripotent stem cells toward synchronous neural networks on an arrayed monolayer of nanofiber membrane

Automatic differentiation of human-induced pluripotent stem cells (hiPSCs) facilitates the generation of cortical neural networks and studies of brain functions. Here, we present a method of directed differentiation of hiPSCs with a substrate made of a honeycomb microframe and a monolayer of crosslinked gelatin nanofibers in the form of an array of nanofiber membranes. Neural precursor cells (NPCs) were firstly derived from hiPSCs and then placed on the nanofiber membranes for automatically controlled neural differentiation over a long period. Due to the strong modulation of the substrate stiffness and permeability, most cells were found in the center area of the honeycomb compartments, giving rise to regular and inter-connected cortical neural clusters. More importantly, the neural activities of the clusters were synchronized proving the reliability of the method. Our results showed that the self-organization, as well as the neural activities of differentiating neural cells, were more efficient in the nanofiber membrane compared to the types of the substrate such as glass and nanofiber-covered glass. In addition to the inherent advantages such as manpower saving and fewer risks of contamination and human error, automatic differentiation avoided undesired shaking which might have critical effects on the formation of synchronous neural clusters. STATEMENT OF SIGNIFICANCE: : Synchronization of cortical neural activities is essential for information processing and human cognition. By automated differentiation of human induced pluripotent stem cells on arrayed monolayer of nanofiber membrane, synchronous neural clusters could be formed. Such an approach would allow creating a variety of neural networks with regular and interconnected clusters for systematic studies of human cortical functions.

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.

Nanoscale Patterning of In Vitro Neuronal Circuits

Methods for patterning neurons in vitro have gradually improved and are used to investigate questions that are difficult to address in or ex vivo. Though these techniques guide axons between groups of neurons, multiscale control of neuronal connectivity, from circuits to synapses, is yet to be achieved in vitro. As studying neuronal circuits with synaptic resolution in vivo poses significant challenges, we present an in vitro alternative to validate biophysical and computational models. In this work we use a combination of electron beam lithography and photolithography to create polydimethylsiloxane (PDMS) structures with features ranging from 150 nm to a few millimeters. Leveraging the difference between average axon and dendritic spine diameters, we restrict axon growth while allowing spines to pass through nanochannels to guide synapse formation between small groups of neurons (i.e., nodes). We show this technique can be used to generate large numbers of isolated feed-forward circuits where connections between nodes are restricted to regions connected by nanochannels. Using a genetically encoded calcium indicator in combination with fluorescently tagged postsynaptic protein, PSD-95, we demonstrate functional synapses can form in this region.

Early lock-in of structured and specialised information flows during neural development

The brains of many organisms are capable of complicated distributed computation underpinned by a highly advanced information processing capacity. Although substantial progress has been made towards characterising the information flow component of this capacity in mature brains, there is a distinct lack of work characterising its emergence during neural development. This lack of progress has been largely driven by the lack of effective estimators of information processing operations for spiking data. Here, we leverage recent advances in this estimation task in order to quantify the changes in transfer entropy during development. We do so by studying the changes in the intrinsic dynamics of the spontaneous activity of developing dissociated neural cell cultures. We find that the quantity of information flowing across these networks undergoes a dramatic increase across development. Moreover, the spatial structure of these flows exhibits a tendency to lock-in at the point when they arise. We also characterise the flow of information during the crucial periods of population bursts. We find that, during these bursts, nodes tend to undertake specialised computational roles as either transmitters, mediators, or receivers of information, with these roles tending to align with their average spike ordering. Further, we find that these roles are regularly locked-in when the information flows are established. Finally, we compare these results to information flows in a model network developing according to a spike-timing-dependent plasticity learning rule. Similar temporal patterns in the development of information flows were observed in these networks, hinting at the broader generality of these phenomena.

Neurons-on-a-Chip: In Vitro NeuroTools

Neurons-on-a-Chip technology has been developed to provide diverse in vitro neuro-tools to study neuritogenesis, synaptogensis, axon guidance, and network dynamics. The two core enabling technologies are soft-lithography and microelectrode array technology. Soft lithography technology made it possible to fabricate microstamps and microfluidic channel devices with a simple replica molding method in a biological laboratory and innovatively reduced the turn-around time from assay design to chip fabrication, facilitating various experimental designs. To control nerve cell behaviors at the single cell level via chemical cues, surface biofunctionalization methods and micropatterning techniques were developed. Microelectrode chip technology, which provides a functional readout by measuring the electrophysiological signals from individual neurons, has become a popular platform to investigate neural information processing in networks. Due to these key advances, it is possible to study the relationship between the network structure and functions, and they have opened a new era of neurobiology and will become standard tools in the near future.

A multilayer network model of neuron-astrocyte populations in vitro reveals mGluR5 inhibition is protective following traumatic injury

Astrocytes communicate bidirectionally with neurons, enhancing synaptic plasticity and promoting the synchronization of neuronal microcircuits. Despite recent advances in understanding neuron-astrocyte signaling, little is known about astrocytic modulation of neuronal activity at the population level, particularly in disease or following injury. We used high-speed calcium imaging of mixed cortical cultures in vitro to determine how population activity changes after disruption of glutamatergic signaling and mechanical injury. We constructed a multilayer network model of neuron-astrocyte connectivity, which captured distinct topology and response behavior from single-cell-type networks. mGluR₅ inhibition decreased neuronal activity, but did not on its own disrupt functional connectivity or network topology. In contrast, injury increased the strength, clustering, and efficiency of neuronal but not astrocytic networks, an effect that was not observed in networks pretreated with mGluR₅ inhibition. Comparison of spatial and functional connectivity revealed that functional connectivity is largely independent of spatial proximity at the microscale, but mechanical injury increased the spatial-functional correlation. Finally, we found that astrocyte segments of the same cell often belong to separate functional communities based on neuronal connectivity, suggesting that astrocyte segments function as independent entities. Our findings demonstrate the utility of multilayer network models for characterizing the multiscale connectivity of two distinct but functionally dependent cell populations.

Astrocytes communicate bidirectionally with neurons, enhancing synaptic plasticity and promoting the synchronization of neuronal microcircuits. We constructed a multilayer network model of neuron-astrocyte connectivity based on calcium activity in mixed cortical cultures, and used this model to evaluate the effect of glutamatergic inhibition and mechanical injury on network topology. We found that injury increased the strength, clustering, and efficiency of neuronal but not astrocytic networks, an effect that was not observed in injured networks pretreated with a glutamate receptor antagonist. Our findings demonstrate the utility of multilayer network models for characterizing the multiscale connectivity of two distinct but functionally dependent cell populations.

Lipopolysaccharide-induced neuroinflammation disrupts functional connectivity and community structure in primary cortical microtissues

Three-dimensional (3D) neural microtissues are a powerful in vitro paradigm for studying brain development and disease under controlled conditions, while maintaining many key attributes of the in vivo environment. Here, we used primary cortical microtissues to study the effects of neuroinflammation on neural microcircuits. We demonstrated the use of a genetically encoded calcium indicator combined with a novel live-imaging platform to record spontaneous calcium transients in microtissues from day 14-34 in vitro. We implemented graph theory analysis of calcium activity to characterize underlying functional connectivity and community structure of microcircuits, which are capable of capturing subtle changes in network dynamics during early disease states. We found that microtissues cultured for 34 days displayed functional remodeling of microcircuits and that community structure strengthened over time. Lipopolysaccharide, a neuroinflammatory agent, significantly increased functional connectivity and disrupted community structure 5-9 days after exposure. These microcircuit-level changes have broad implications for the role of neuroinflammation in functional dysregulation of neural networks.

Early prediction of developing spontaneous activity in cultured neuronal networks

Synchronization and bursting activity are intrinsic electrophysiological properties of in vivo and in vitro neural networks. During early development, cortical cultures exhibit a wide repertoire of synchronous bursting dynamics whose characterization may help to understand the parameters governing the transition from immature to mature networks. Here we used machine learning techniques to characterize and predict the developing spontaneous activity in mouse cortical neurons on microelectrode arrays (MEAs) during the first three weeks in vitro. Network activity at three stages of early development was defined by 18 electrophysiological features of spikes, bursts, synchrony, and connectivity. The variability of neuronal network activity during early development was investigated by applying k-means and self-organizing map (SOM) clustering analysis to features of bursts and synchrony. These electrophysiological features were predicted at the third week in vitro with high accuracy from those at earlier times using three machine learning models: Multivariate Adaptive Regression Splines, Support Vector Machines, and Random Forest. Our results indicate that initial patterns of electrical activity during the first week in vitro may already predetermine the final development of the neuronal network activity. The methodological approach used here may be applied to explore the biological mechanisms underlying the complex dynamics of spontaneous activity in developing neuronal cultures.

Functional Inhibitory Connections Modulate the Electrophysiological Activity Patterns of Cortical-Hippocampal Ensembles

The brain is a complex organ composed of billions of neurons connected through excitatory and inhibitory synapses. Its structure reveals a modular topological organization, where neurons are arranged in interconnected assemblies. The generated patterns of electrophysiological activity are shaped by two main factors: network heterogeneity and the topological properties of the underlying connectivity that strongly push the dynamics toward different brain-states. In this work, we exploited an innovative polymeric structure coupled to Micro-Electrode Arrays (MEAs) to recreate in vitro heterogeneous interconnected (modular) neuronal networks made up of cortical and hippocampal neurons. We investigated the propagation of spike sequences between the two interconnected subpopulations during the networks' development, correlating functional and structural connectivity to dynamics. The simultaneous presence of two neuronal types shaped the features of the functional connections (excitation vs. inhibition), orchestrating the emerging patterns of electrophysiological activity. In particular, we found that hippocampal neurons mostly project inhibitory connections toward the cortical counterpart modulating the temporal scale of the population events (network bursts). In contrast, cortical neurons establish a larger amount of intrapopulation connections. Moreover, we proved topological properties such as small-worldness, degree distribution, and modularity of neuronal assemblies were favored by the physical environment where networks developed and matured.

The Flow of Axonal Information Among Hippocampal Subregions: 1. Feed-Forward and Feedback Network Spatial Dynamics Underpinning Emergent Information Processing

The tri-synaptic pathway in the mammalian hippocampus enables cognitive learning and memory. Despite decades of reports on anatomy and physiology, the functional architecture of the hippocampal network remains poorly understood in terms of the dynamics of axonal information transfer between subregions. Information inputs largely flow from the entorhinal cortex (EC) to the dentate gyrus (DG), and then are processed further in the CA3 and CA1 before returning to the EC. Here, we reconstructed elements of the rat hippocampus in a novel device over an electrode array that allowed for monitoring the directionality of individual axons between the subregions. The direction of spike propagation was determined by the transmission delay of the axons recorded between two electrodes in microfluidic tunnels. The majority of axons from the EC to the DG operated in the feed-forward direction, with other regions developing unexpectedly large proportions of feedback axons to balance excitation. Spike timing in axons between each region followed single exponential log-log distributions over two orders of magnitude from 0.01 to 1 s, indicating that conventional descriptors of mean firing rates are misleading assumptions. Most of the spiking occurred in bursts that required two exponentials to fit the distribution of inter-burst intervals. This suggested the presence of up-states and down-states in every region, with the least up-states in the DG to CA3 feed-forward axons and the CA3 subregion. The peaks of the log-normal distributions of intra-burst spike rates were similar in axons between regions with modes around 95 Hz distributed over an order of magnitude. Burst durations were also log-normally distributed around a peak of 88 ms over two orders of magnitude. Despite the diversity of these spike distributions, spike rates from individual axons were often linearly correlated to subregions. These linear relationships enabled the generation of structural connectivity graphs, not possible previously without the directional flow of axonal information. The rich axonal spike dynamics between subregions of the hippocampus reveal both constraints and broad emergent dynamics of hippocampal architecture. Knowledge of this network architecture may enable more efficient computational artificial intelligence (AI) networks, neuromorphic hardware, and stimulation and decoding from cognitive implants.

Improved Method for Efficient Generation of Functional Neurons from Murine Neural Progenitor Cells

Neuronal culture was used to investigate neuronal function in physiological and pathological conditions. Despite its inevitability, primary neuronal culture remained a gold standard method that requires laborious preparation, intensive training, and animal resources. To circumvent the shortfalls of primary neuronal preparations and efficiently give rise to functional neurons, we combine a neural stem cell culture method with a direct cell type-conversion approach. The lucidity of this method enables the efficient preparation of functional neurons from mouse neural progenitor cells on demand. We demonstrate that induced neurons (NPC-iNs) by this method make synaptic connections, elicit neuronal activity-dependent cellular responses, and develop functional neuronal networks. This method will provide a concise platform for functional neuronal assessments. This indeed offers a perspective for using these characterized neuronal networks for investigating plasticity mechanisms, drug screening assays, and probing the molecular and biophysical basis of neurodevelopmental and neurodegenerative diseases.

A generative network model of neurodevelopmental diversity in structural brain organization

The formation of large-scale brain networks, and their continual refinement, represent crucial developmental processes that can drive individual differences in cognition and which are associated with multiple neurodevelopmental conditions. But how does this organization arise, and what mechanisms drive diversity in organization? We use generative network modeling to provide a computational framework for understanding neurodevelopmental diversity. Within this framework macroscopic brain organization, complete with spatial embedding of its organization, is an emergent property of a generative wiring equation that optimizes its connectivity by renegotiating its biological costs and topological values continuously over time. The rules that govern these iterative wiring properties are controlled by a set of tightly framed parameters, with subtle differences in these parameters steering network growth towards different neurodiverse outcomes. Regional expression of genes associated with the simulations converge on biological processes and cellular components predominantly involved in synaptic signaling, neuronal projection, catabolic intracellular processes and protein transport. Together, this provides a unifying computational framework for conceptualizing the mechanisms and diversity in neurodevelopment, capable of integrating different levels of analysis-from genes to cognition.

Neural anatomy and optical microscopy (NAOMi) simulation for evaluating calcium imaging methods

BACKGROUND

The past decade has seen a multitude of new in vivo functional imaging methodologies. However, the lack of ground-truth comparisons or evaluation metrics makes the large-scale, systematic validation vital to the continued development and use of optical microscopy impossible.

NEW-METHOD

We provide a new framework for evaluating two-photon microscopy methods via in silico Neural Anatomy and Optical Microscopy (NAOMi) simulation. Our computationally efficient model generates large anatomical volumes of mouse cortex, simulates neural activity, and incorporates optical propagation and scanning to create realistic calcium imaging datasets.

RESULTS

We verify NAOMi simulations against in vivo two-photon recordings from mouse cortex. We leverage this in silico ground truth to directly compare different segmentation algorithms and optical designs. We find modern segmentation algorithms extract strong neural time-courses comparable to estimation using oracle spatial information, but with an increase in the false positive rate. Comparison between optical setups demonstrate improved resilience to motion artifacts in sparsely labeled samples using Bessel beams, increased signal-to-noise ratio and cell-count using low numerical aperture Gaussian beams and nuclear GCaMP, and more uniform spatial sampling with temporal focusing versus multi-plane imaging.

COMPARISON WITH EXISTING METHODS

NAOMi is a first-of-its kind framework for assessing optical imaging modalities. Existing methods are either anatomical simulations or do not address functional imaging. Thus there is no competing method for simulating realistic functional optical microscopy data.

CONCLUSIONS

By leveraging the rich accumulated knowledge of neural anatomy and optical physics, we provide a powerful new tool to assess and develop important methods in neural imaging.

Dopamine Transporter Is a Master Regulator of Dopaminergic Neural Network Connectivity

Dopaminergic neurons of the substantia nigra pars compacta (SNC) and ventral tegmental area (VTA) exhibit spontaneous firing activity. The dopaminergic neurons in these regions have been shown to exhibit differential sensitivity to neuronal loss and psychostimulants targeting dopamine transporter. However, it remains unclear whether these regional differences scale beyond individual neuronal activity to regional neuronal networks. Here, we used live-cell calcium imaging to show that network connectivity greatly differs between SNC and VTA regions with higher incidence of hub-like neurons in the VTA. Specifically, the frequency of hub-like neurons was significantly lower in SNC than in the adjacent VTA, consistent with the interpretation of a lower network resilience to SNC neuronal loss. We tested this hypothesis, in DAT-cre/loxP-GCaMP6f mice of either sex, when activity of an individual dopaminergic neuron is suppressed, through whole-cell patch clamp electrophysiology, in either SNC or VTA networks. Neuronal loss in the SNC increased network clustering, whereas the larger number of hub-neurons in the VTA overcompensated by decreasing network clustering in the VTA. We further show that network properties are regulatable via a dopamine transporter but not a D2 receptor dependent mechanism. Our results demonstrate novel regulatory mechanisms of functional network topology in dopaminergic brain regions.SIGNIFICANCE STATEMENT In this work, we begin to untangle the differences in complex network properties between the substantia nigra pars compacta (SNC) and VTA, that may underlie differential sensitivity between regions. The methods and analysis employed provide a springboard for investigations of network topology in multiple deep brain structures and disorders.

Identifying Emergent Mesoscopic-Macroscopic Functional Brain Network Dynamics in Infants at Term-Equivalent Age with Electric Source Neuroimaging

Aim: To identify and characterize the functional brain networks at the time when the brain is yet to develop higher order functions in term-born and preterm infants at term-equivalent age. Introduction: Although functional magnetic resonance imaging (fMRI) data have revealed the existence of spatially structured resting-state brain activity in infants, the temporal information of fMRI data limits the characterization of fast timescale brain oscillations. In this study, we use infants' high-density electroencephalography (EEG) to characterize spatiotemporal and spectral functional organizations of brain network dynamics. Methods: We used source-reconstructed EEG and graph theoretical analyses in 100 infants (84 preterm, 16 term born) to identify the rich-club topological organization, temporal dynamics, and spectral fingerprints of dynamic functional brain networks. Results: Five dynamic functional brain networks are identified, which have rich-club topological organizations, distinctive spectral fingerprints (in the delta and low-alpha frequency), and scale-invariant temporal dynamics (<0.1 Hz): The default mode, primary sensory-limbic system, thalamo-frontal, thalamo-sensorimotor, and visual-limbic system. The temporal dynamics of these networks are correlated in a hierarchically leading-following organization, showing that infant brain networks arise from long-range synchronization of band-limited cortical oscillation based on interacting fast- and slow-coherent cortical oscillations. Conclusion: Dynamic functional brain networks do not solely depend on the maturation of cognitive networks; instead, the brain network dynamics exist in infants at term age well before the childhood and adulthood, and hence, it offers a quantitative measurement of neurotypical development in infants. Clinical Trial Registration Number: ACTRN12615000591550. Impact statement Our work offers novel functional insights into the brain network characterization in infants, providing a new functional basis for future deployable prognostication approaches.

Inhibitory neurons exhibit high controlling ability in the cortical microconnectome

The brain is a network system in which excitatory and inhibitory neurons keep activity balanced in the highly non-random connectivity pattern of the microconnectome. It is well known that the relative percentage of inhibitory neurons is much smaller than excitatory neurons in the cortex. So, in general, how inhibitory neurons can keep the balance with the surrounding excitatory neurons is an important question. There is much accumulated knowledge about this fundamental question. This study quantitatively evaluated the relatively higher functional contribution of inhibitory neurons in terms of not only properties of individual neurons, such as firing rate, but also in terms of topological mechanisms and controlling ability on other excitatory neurons. We combined simultaneous electrical recording (~2.5 hours) of ~1000 neurons in vitro, and quantitative evaluation of neuronal interactions including excitatory-inhibitory categorization. This study accurately defined recording brain anatomical targets, such as brain regions and cortical layers, by inter-referring MRI and immunostaining recordings. The interaction networks enabled us to quantify topological influence of individual neurons, in terms of controlling ability to other neurons. Especially, the result indicated that highly influential inhibitory neurons show higher controlling ability of other neurons than excitatory neurons, and are relatively often distributed in deeper layers of the cortex. Furthermore, the neurons having high controlling ability are more effectively limited in number than central nodes of k-cores, and these neurons also participate in more clustered motifs. In summary, this study suggested that the high controlling ability of inhibitory neurons is a key mechanism to keep balance with a large number of other excitatory neurons beyond simple higher firing rate. Application of the selection method of limited important neurons would be also applicable for the ability to effectively and selectively stimulate E/I imbalanced disease states.

How small numbers of inhibitory neurons functionally keep balance with large numbers of excitatory neurons in the brain by controlling each other is a fundamental question. Especially, this study quantitatively evaluated a topological mechanism of interaction networks in terms of controlling abilities of individual cortical neurons to other neurons. Combination of simultaneous electrical recording of ~1000 neurons and a quantitative evaluation method of neuronal interactions including excitatory-inhibitory categories, enabled us to evaluate the influence of individual neurons not only about firing rate but also about their relative positions in the networks and controllable ability of other neurons. Especially, the result showed that inhibitory neurons have more controlling ability than excitatory neurons, and such neurons were more often observed in deep layers. Because the limited number of neurons in terms controlling ability were much smaller than neurons based on centrality measure and, of course, more directly selected neurons based on their ability to control other neurons, the selection method of important neurons will help not only to produce realistic computational models but also will help to stimulate brain to effectively treat imbalanced disease states.

Configuration and dynamics of dominant inspiratory multineuronal activity patterns during eupnea and gasping generation in vitro

The pre-Bötzinger complex (preBötC), located within the ventral respiratory column, produces inspiratory bursts in varying degrees of synchronization/amplitude. This wide range of population burst patterns reflects the flexibility of the preBötC neurons, which is expressed in variations in the onset/offset times of their activations and their activity during the population bursts, with respiratory neurons exhibiting a large cycle-to-cycle timing jitter both at the population activity onset and at the population activity peak, suggesting that respiratory neurons are stochastically activated before and during the inspiratory bursts. However, it is still unknown whether this stochasticity is maintained while evaluating the coactivity of respiratory neuronal ensembles. Moreover, the preBötC topology also remains unknown. In this study, by simultaneously recording tens of preBötC neurons and using coactivation analysis during the inspiratory periods, we found that the preBötC has a scale-free configuration (mixture of not many highly connected nodes, hubs, with abundant poorly connected elements) exhibiting the rich-club phenomenon (hubs more likely interconnected with each other). PreBötC neurons also produce multineuronal activity patterns (MAPs) that are highly stable and change during the hypoxia-induced reconfiguration. Moreover, preBötC contains a coactivating core network shared by all its MAPs. Finally, we found a distinctive pattern of sequential coactivation of core network neurons at the beginning of the inspiratory periods, indicating that, when evaluated at the multicellular level, the coactivation of respiratory neurons seems not to be stochastic.NEW & NOTEWORTHY By means of multielectrode recordings of preBötC neurons, we evaluated their configuration in normoxia and hypoxia, finding that the preBötC exhibits a scale-free configuration with a rich-club phenomenon. preBötC neurons produce multineuronal activity patterns that are highly stable but change during hypoxia. The preBötC contains a coactivating core network that exhibit a distinctive pattern of coactivation at the beginning of inspirations. These results reveal some network basis of inspiratory rhythm generation and its reconfiguration during hypoxia.

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.

Functional Characterization of Three-Dimensional Cortical Cultures for In Vitro Modeling of Brain Networks

Three-dimensional (3D) in vitro cultures recapitulate key features of the brain including morphology, cell-cell and cell-extracellular matrix interactions, gradients of factors, and mechanical properties. However, there remains a need for experimental and computational tools to investigate network functions in these 3D models. To address this need, we present an experimental system based on 3D scaffold-based cortical neuron cultures in which we expressed the genetically encoded calcium indicator GCaMP6f to record neuronal activity at the millimeter-scale. Functional neural network descriptors were computed with graph-theory-based network analysis methods, showing the formation of functional networks at 3 weeks of culture. Changes to the functional network properties upon perturbations to glutamatergic neurotransmission or GABAergic neurotransmission were quantitatively characterized. The results illustrate the applicability of our 3D experimental system for the study of brain network development, function, and disruption in a biomimetic microenvironment.

Network neuroscience of apathy in cerebrovascular disease

Apathy is a reduction in motivated goal-directed behavior (GDB) that is prevalent in cerebrovascular disease, providing an important opportunity to study the mechanistic underpinnings of motivation in humans. Focal lesions, such as those seen in stroke, have been crucial in developing models of brain regions underlying motivated behavior, while studies of cerebral small vessel disease (SVD) have helped define the connections between brain regions supporting such behavior. However, current lesion-based models cannot fully explain the neurobiology of apathy in stroke and SVD. To address this, we propose a network-based model which conceptualizes apathy as the result of damage to GDB-related networks. A review of the current evidence suggests that cerebrovascular disease-related pathology can lead to network changes outside of initially damaged territories, which may propagate to regions that share structural or functional connections. The presentation and longitudinal trajectory of apathy in stroke and SVD may be the result of these network changes. Distinct subnetworks might support cognitive components of GDB, the disruption of which results in specific symptoms of apathy. This network-based model of apathy may open new approaches for investigating its underlying neurobiology, and presents novel opportunities for its diagnosis and treatment.

Modeling the temporal network dynamics of neuronal cultures

Neurons form complex networks that evolve over multiple time scales. In order to thoroughly characterize these networks, time dependencies must be explicitly modeled. Here, we present a statistical model that captures both the underlying structural and temporal dynamics of neuronal networks. Our model combines the class of Stochastic Block Models for community formation with Gaussian processes to model changes in the community structure as a smooth function of time. We validate our model on synthetic data and demonstrate its utility on three different studies using in vitro cultures of dissociated neurons.

Multiscale Neuroscience of Psychiatric Disorders

The human brain comprises a multiscale network with multiple levels of organization. Neurons with dendritic and axonal connections form the microscale fabric of brain circuitry, and macroscale brain regions and white matter connections form the infrastructure for system-level brain communication and information integration. In this review, we discuss the emerging trend of multiscale neuroscience, the multidisciplinary field that brings together data from these different levels of nervous system organization to form a better understanding of between-scale relationships of brain structure, function, and behavior in health and disease. We provide a broad overview of this developing field and discuss recent findings of exemplary multiscale neuroscience studies that illustrate the importance of studying cross-scale interactions among the genetic, molecular, cellular, and macroscale levels of brain circuitry and connectivity and behavior. We particularly consider a central, overarching goal of these multiscale neuroscience studies of human brain connectivity: to obtain insight into how disease-related alterations at one level of organization may underlie alterations observed at other scales of brain network organization in mental disorders. We conclude by discussing the current limitations, challenges, and future directions of the field.

Methamphetamine regulation of activity and topology of ventral midbrain networks

The ventral midbrain supports a variety of functions through the heterogeneity of neurons. Dopaminergic and GABA neurons within this region are particularly susceptible targets of amphetamine-class psychostimulants such as methamphetamine. While this has been evidenced through single-neuron methods, it remains unclear whether and to what extent the local neuronal network is affected and if so, by which mechanisms. Both GABAergic and dopaminergic neurons were heavily featured within the primary ventral midbrain network model system. Using spontaneous calcium activity, our data suggest methamphetamine decreased total network output via a D2 receptor-dependent manner. Over culture duration, functional connectivity between neurons decreased significantly but was unaffected by methamphetamine. However, across culture duration, exposure to methamphetamine significantly altered changes in network assortativity. Here we have established primary ventral midbrain networks culture as a viable model system that reveals specific changes in network activity, connectivity, and topology modulation by methamphetamine. This network culture system enables control over the type and number of neurons that comprise a network and facilitates detection of emergent properties that arise from the specific organization. Thus, the multidimensional properties of methamphetamine can be unraveled, leading to a better understanding of its impact on the local network structure and function.

Active High-Density Electrode Arrays: Technology and Applications in Neuronal Cell Cultures

Active high-density electrode arrays realized with complementary metal-oxide-semiconductor (CMOS) technology provide electrophysiological recordings from several thousands of closely spaced microelectrodes. This has drastically advanced the spatiotemporal recording resolution of conventional multielectrode arrays (MEAs). Thus, today's electrophysiology in neuronal cultures can exploit label-free electrical readouts from a large number of single neurons within the same network. This provides advanced capabilities to investigate the properties of self-assembling neuronal networks, to advance studies on neurotoxicity and neurodevelopmental alterations associated with human brain diseases, and to develop cell culture models for testing drug- or cell-based strategies for therapies.Here, after introducing the reader to this neurotechnology, we summarize the results of different recent studies demonstrating the potential of active high-density electrode arrays for experimental applications. We also discuss ongoing and possible future research directions that might allow for moving these platforms forward for screening applications.

The effect of acute glutamate treatment on the functional connectivity and network topology of cortical cultures

Microelectrode arrays (MEAs) allow the investigation of the pharmacological and toxicological effects of chemicals on cultured neuronal networks. Understanding the functional connections between neurons and the resulting neuronal networks is important for evaluating drugs that affect synaptic transmission. Therefore, we acutely treated a mature cultured neuronal network on MEAs with accumulating amounts of glutamate and recorded their altered electrophysiology. Subsequently, a cross-covariance analysis was applied to process the spiking activity in the network and to evaluate the connections between neurons. Finally, graph theory was used to assess the functional network properties under acute glutamate treatment. Our data demonstrated that glutamate increased the similarity, connectivity weight, density, and largest-component size of the functional network. In addition, the small-world network topology was altered after glutamate treatment. Our results indicate that the graph theory can advance our understanding of the pharmacological significance of neurotransmitters on neuronal networks.

Multiplane Calcium Imaging Reveals Disrupted Development of Network Topology in Zebrafish pcdh19 Mutants

Functional brain networks self-assemble during development, although the molecular basis of network assembly is poorly understood. Protocadherin-19 (pcdh19) is a homophilic cell adhesion molecule that is linked to neurodevelopmental disorders, and influences multiple cellular and developmental events in zebrafish. Although loss of PCDH19 in humans and model organisms leads to functional deficits, the underlying network defects remain unknown. Here, we employ multiplane, resonant-scanning in vivo two-photon calcium imaging of developing zebrafish, and use graph theory to characterize the development of resting state functional networks in both wild-type and pcdh19 mutant larvae. We find that the brain networks of pcdh19 mutants display enhanced clustering and an altered developmental trajectory of network assembly. Our results show that functional imaging and network analysis in zebrafish larvae is an effective approach for characterizing the developmental impact of lesions in genes of clinical interest.