Dissecting functional connectivity of neuronal microcircuits: experimental and theoretical insights

Disrupting glioblastoma networks with tumor treating fields (TTFields) in in vitro models

PURPOSE

This study investigates the biological effect of Tumor Treating Fields (TTFields) on key drivers of glioblastoma's malignancy-tumor microtube (TM) formation-and on the function and overall integrity of the tumor cell network.

METHOD

Using a two-dimensional monoculture GB cell network model (2DTM) of primary glioblastoma cell (GBC) cultures (S24, BG5 or T269), we evaluated the effects of TTFields on cell density, interconnectivity and structural integrity of the tumor network. We also analyzed calcium (Ca2+) transient dynamics and network morphology, validating findings in patient-derived tumoroids and brain tumor organoids.

RESULTS

In the 2DTM assay, TTFields reduced cell density by 85-88% and disrupted network interconnectivity, particularly in cells with multiple TMs. A "crooked TM" phenotype emerged in 5-6% of treated cells, rarely seen in controls. Ca2+ transients were significantly compromised, with global Ca2+ activity reduced by 51-83%, active and periodic cells by over 50%, and intercellular co-activity by 52% in S24, and almost completely in BG5 GBCs. The effects were more pronounced at 200 kHz compared to a 50 kHz TTFields. Similar reductions in Ca2+ activity were observed in patient-derived tumoroids. In brain tumor organoids, TTFields significantly reduced tumor cell proliferation and infiltration.

CONCLUSION

Our comprehensive study provides new insights into the multiple effects of Inovitro-modeled TTFields on glioma progression, morphology and network dynamics in vitro. Future in vivo studies to verify our in vitro findings may provide the basis for a deeper understanding and optimization of TTFields as a therapeutic modality in the treatment of GB.

Association study of brain structure-function coupling and glymphatic system function in patients with mild cognitive impairment due to Alzheimer's disease

Background

Mild cognitive impairment (MCI) is a critical transitional phase from healthy cognitive aging to dementia, offering a unique opportunity for early intervention. However, few studies focus on the correlation of brain structure and functional activity in patients with MCI due to Alzheimer's disease (AD). Elucidating the complex interactions between structural-functional (SC-FC) brain connectivity and glymphatic system function is crucial for understanding this condition.

Method

The aims of this study were to explore the relationship among SC-FC coupling values, glymphatic system function and cognitive function. 23 MCI patients and 18 healthy controls (HC) underwent diffusion tensor imaging (DTI) and resting-state functional MRI (fMRI). DTI analysis along the perivascular space (DTI-ALPS) index and SC-FC coupling values were calculated using DTI and fMRI. Correlation analysis was conducted to assess the relationship between Mini-Mental State Examination (MMSE) scores, DTI-ALPS index, and coupling values. Receiver operating characteristic (ROC) curves was conducted on the SC-FC coupling between the whole brain and subnetworks. The correlation of coupling values with MMSE scores was also analyzed.

Result

MCI patients (67.74 ± 6.99 years of age) exhibited significantly lower coupling in the whole-brain network and subnetworks, such as the somatomotor network (SMN) and ventral attention network (VAN), than HCs (63.44 ± 6.92 years of age). Whole-brain network coupling was positively correlated with dorsal attention network (DAN), SMN, and visual network (VN) coupling. MMSE scores were significantly positively correlated with whole-brain coupling and SMN coupling. In MCI, whole-brain network demonstrated the highest performance, followed by the SMN and VAN, with the VN, DAN, limbic network (LN), frontoparietal network (FPN), and default mode network (DMN). Compared to HCs, lower DTI-ALPS index was observed in individuals with MCI. Additionally, the left DTI-ALPS index showed a significant positive correlation with MMSE scores and coupling values in the whole-brain network and SMN.

Conclusion

These findings reveal the critical role of SC-FC coupling values and the ALPS index in cognitive function of MCI. The positive correlations observed in the left DTI-ALPS and whole-brain and SMN coupling values provide a new insight for investigating the asymmetrical nature of cognitive impairments.

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.

Exploring Kainic Acid-Induced Alterations in Circular Tripartite Networks with Advanced Analysis Tools

Brain activity implies the orchestrated functioning of interconnected brain regions. Typical in vitro models aim to mimic the brain using single human pluripotent stem cell-derived neuronal networks. However, the field is constantly evolving to model brain functions more accurately through the use of new paradigms, e.g., brain-on-a-chip models with compartmentalized structures and integrated sensors. These methods create novel data requiring more complex analysis approaches. The previously introduced circular tripartite network concept models the connectivity between spatially diverse neuronal structures. The model consists of a microfluidic device allowing axonal connectivity between separated neuronal networks with an embedded microelectrode array to record both local and global electrophysiological activity patterns in the closed circuitry. The existing tools are suboptimal for the analysis of the data produced with this model. Here, we introduce advanced tools for synchronization and functional connectivity assessment. We used our custom-designed analysis to assess the interrelations between the kainic acid (KA)-exposed proximal compartment and its nonexposed distal neighbors before and after KA. Novel multilevel circuitry bursting patterns were detected and analyzed in parallel with the inter- and intracompartmental functional connectivity. The effect of KA on the proximal compartment was captured, and the spread of this effect to the nonexposed distal compartments was revealed. KA induced divergent changes in bursting behaviors, which may be explained by distinct baseline activity and varied intra- and intercompartmental connectivity strengths. The circular tripartite network concept combined with our developed analysis advances importantly both face and construct validity in modeling human epilepsy in vitro.

A clinically applicable connectivity signature for glioblastoma includes the tumor network driver CHI3L1

Tumor microtubes (TMs) connect glioma cells to a network with considerable relevance for tumor progression and therapy resistance. However, the determination of TM-interconnectivity in individual tumors is challenging and the impact on patient survival unresolved. Here, we establish a connectivity signature from single-cell RNA-sequenced (scRNA-Seq) xenografted primary glioblastoma (GB) cells using a dye uptake methodology, and validate it with recording of cellular calcium epochs and clinical correlations. Astrocyte-like and mesenchymal-like GB cells have the highest connectivity signature scores in scRNA-sequenced patient-derived xenografts and patient samples. In large GB cohorts, TM-network connectivity correlates with the mesenchymal subtype and dismal patient survival. CHI3L1 gene expression serves as a robust molecular marker of connectivity and functionally influences TM networks. The connectivity signature allows insights into brain tumor biology, provides a proof-of-principle that tumor cell TM-connectivity is relevant for patients' prognosis, and serves as a robust prognostic biomarker.

Prominent in vivo influence of single interneurons in the developing barrel cortex

Spontaneous synchronous activity is a hallmark of developing brain circuits and promotes their formation. Ex vivo, synchronous activity was shown to be orchestrated by a sparse population of highly connected GABAergic 'hub' neurons. The recent development of all-optical methods to record and manipulate neuronal activity in vivo now offers the unprecedented opportunity to probe the existence and function of hub cells in vivo. Using calcium imaging, connectivity analysis and holographic optical stimulation, we show that single GABAergic, but not glutamatergic, neurons influence population dynamics in the barrel cortex of non-anaesthetized mouse pups. Single GABAergic cells mainly exert an inhibitory influence on both spontaneous and sensory-evoked population bursts. Their network influence scales with their functional connectivity, with highly connected hub neurons displaying the strongest impact. We propose that hub neurons function in tailoring intrinsic cortical dynamics to external sensory inputs.

Stimulus-dependent functional network topology in mouse visual cortex

Information is processed by networks of neurons in the brain. On the timescale of sensory processing, those neuronal networks have relatively fixed anatomical connectivity, while functional connectivity, which defines the interactions between neurons, can vary depending on the ongoing activity of the neurons within the network. We thus hypothesized that different types of stimuli, which drive different neuronal activities in the network, could lead those networks to display stimulus-dependent functional connectivity patterns. To test this hypothesis, we analyzed electrophysiological data from the Allen Brain Observatory, which utilized Neuropixels probes to simultaneously record stimulus-evoked activity from hundreds of neurons across 6 different regions of mouse visual cortex. The recordings had single-cell resolution and high temporal fidelity, enabling us to determine fine-scale functional connectivity. Comparing the functional connectivity patterns observed when different stimuli were presented to the mice, we made several nontrivial observations. First, while the frequencies of different connectivity motifs (i.e., the patterns of connectivity between triplets of neurons) were preserved across stimuli, the identities of the neurons within those motifs changed. This means that functional connectivity dynamically changes along with the input stimulus, but does so in a way that preserves the motif frequencies. Secondly, we found that the degree to which functional modules are contained within a single brain region (as opposed to being distributed between regions) increases with increasing stimulus complexity. This suggests a mechanism for how the brain could dynamically alter its computations based on its inputs. Altogether, our work reveals unexpected stimulus-dependence to the way groups of neurons interact to process incoming sensory information.

Incubator-independent perfusion system integrated with microfluidic device for continuous electrophysiology and microscopy readouts

Advances in primary and stem cell derived neuronal cell culture techniques and abundance of available neuronal cell types have enabledin vitroneuroscience as a substantial approach to modelin vivoneuronal networks. Survival of the cultured neurons is inevitably dependent on the cell culture incubators to provide stable temperature and humidity and to supply required CO₂levels for controlling the pH of culture medium. Therefore, imaging and electrophysiology recordings outside of the incubator are often limited to the short-term experimental sessions. This restricts our understanding of physiological events to the short snapshots of recorded data while the major part of temporal data is neglected. Multiple custom-made and commercially available platforms like integrated on-stage incubators have been designed to enable long-term microscopy. Nevertheless, long-term high-spatiotemporal electrophysiology recordings from developing neuronal networks needs to be addressed. In the present work an incubator-independent polydimethylsiloxane-based double-wall perfusion chamber was designed and integrated with multi-electrode arrays (MEAs) electrophysiology and compartmentalized microfluidic device to continuously record from engineered neuronal networks at sub-cellular resolution. Cell culture media underwent iterations of conditioning to the ambient CO₂and adjusting its pH to physiological ranges to retain a stable pH for weeks outside of the incubator. Double-wall perfusion chamber and an integrated air bubble trapper reduced media evaporation and osmolality drifts of the conditioned media for two weeks. Aligned microchannel-microfluidic device on MEA electrodes allowed neurite growth on top of the planar electrodes and amplified their extracellular activity. This enabled continuous sub-cellular resolution imaging and electrophysiology recordings from developing networks and their growing neurites. The on-chip versatile and self-contained system provides long-term, continuous and high spatiotemporal access to the network data and offers a robustin vitroplatform with many potentials to be applied on advanced cell culture systems including organ-on-chip and organoid models.

Changing subplate circuits: Early activity dependent circuit plasticity

Early neural activity in the developing sensory system comprises spontaneous bursts of patterned activity, which is fundamental for sculpting and refinement of immature cortical connections. The crude early connections that are initially refined by spontaneous activity, are further elaborated by sensory-driven activity from the periphery such that orderly and mature connections are established for the proper functioning of the cortices. Subplate neurons (SPNs) are one of the first-born mature neurons that are transiently present during early development, the period of heightened activity-dependent plasticity. SPNs are well integrated within the developing sensory cortices. Their structural and functional properties such as relative mature intrinsic membrane properties, heightened connectivity via chemical and electrical synapses, robust activation by neuromodulatory inputs-place them in an ideal position to serve as crucial elements in monitoring and regulating spontaneous endogenous network activity. Moreover, SPNs are the earliest substrates to receive early sensory-driven activity from the periphery and are involved in its modulation, amplification, and transmission before the maturation of the direct adult-like thalamocortical connectivity. Consequently, SPNs are vulnerable to sensory manipulations in the periphery. A broad range of early sensory deprivations alters SPN circuit organization and functions that might be associated with long term neurodevelopmental and psychiatric disorders. Here we provide a comprehensive overview of SPN function in activity-dependent development during early life and integrate recent findings on the impact of early sensory deprivation on SPNs that could eventually lead to neurodevelopmental disorders.

Predicting in vitro single-neuron firing rates upon pharmacological perturbation using Graph Neural Networks

Modern Graph Neural Networks (GNNs) provide opportunities to study the determinants underlying the complex activity patterns of biological neuronal networks. In this study, we applied GNNs to a large-scale electrophysiological dataset of rodent primary neuronal networks obtained by means of high-density microelectrode arrays (HD-MEAs). HD-MEAs allow for long-term recording of extracellular spiking activity of individual neurons and networks and enable the extraction of physiologically relevant features at the single-neuron and population level. We employed established GNNs to generate a combined representation of single-neuron and connectivity features obtained from HD-MEA data, with the ultimate goal of predicting changes in single-neuron firing rate induced by a pharmacological perturbation. The aim of the main prediction task was to assess whether single-neuron and functional connectivity features, inferred under baseline conditions, were informative for predicting changes in neuronal activity in response to a perturbation with Bicuculline, a GABA A receptor antagonist. Our results suggest that the joint representation of node features and functional connectivity, extracted from a baseline recording, was informative for predicting firing rate changes of individual neurons after the perturbation. Specifically, our implementation of a GNN model with inductive learning capability (GraphSAGE) outperformed other prediction models that relied only on single-neuron features. We tested the generalizability of the results on two additional datasets of HD-MEA recordings-a second dataset with cultures perturbed with Bicuculline and a dataset perturbed with the GABA A receptor antagonist Gabazine. GraphSAGE models showed improved prediction accuracy over other prediction models. Our results demonstrate the added value of taking into account the functional connectivity between neurons and the potential of GNNs to study complex interactions between neurons.

Differential Patterns of Change in Brain Connectivity Resulting from Severe Traumatic Brain Injury

Background: Traumatic brain injury (TBI) damages white matter tracts, disrupting brain network structure and communication. There exists a wide heterogeneity in the pattern of structural damage associated with injury, as well as a large heterogeneity in behavioral outcomes. However, little is known about the relationship between changes in network connectivity and clinical outcomes. Materials and Methods: We utilize the rat lateral fluid-percussion injury model of severe TBI to study differences in brain connectivity in 8 animals that received the insult and 11 animals that received only a craniectomy. Diffusion tensor imaging is performed 5 weeks after the injury and network theory is used to investigate changes in white matter connectivity. Results: We find that (1) global network measures are not able to distinguish between healthy and injured animals; (2) injury induced alterations predominantly exist in a subset of connections (subnetworks) distributed throughout the brain; and (3) injured animals can be divided into subgroups based on changes in network motifs-measures of local structural connectivity. In addition, alterations in predicted functional connectivity indicate that the subgroups have different propensities to synchronize brain activity, which could relate to the heterogeneity of clinical outcomes. Discussion: These results suggest that network measures can be used to quantify progressive changes in brain connectivity due to injury and differentiate among subpopulations with similar injuries, but different pathological trajectories.

Abnormal innate and learned behavior induced by neuron-microglia miscommunication is related to CA3 reconfiguration

Neuron-microglia communication through the Cx3cr1-Cx3cl1 axis is essential for the development and refinement of neural circuits, which determine their function into adulthood. In the present work we set out to extend the behavioral characterization of Cx3cr1^-/- mice evaluating innate behaviors and spatial navigation, both dependent on hippocampal function. Our results show that Cx3cr1-deficient mice, which show some changes in microglial and synaptic terminals morphology and density, exhibit alterations in activities of daily living and in the rapid encoding of novel spatial information that, nonetheless, improves with training. A neural substrate for these cognitive deficiencies was found in the form of synaptic dysfunction in the CA3 region of the hippocampus, with a marked impact on the mossy fiber (MF) pathway. A network analysis of the CA3 microcircuit reveals the effect of these synaptic alterations on the functional connectivity among CA3 neurons with diminished strength and topological reorganization in Cx3cr1-deficient mice. Neonatal population activity of the CA3 region in Cx3cr1-deficient mice shows a marked reorganization around the giant depolarizing potentials, the first form of network-driven activity of the hippocampus, suggesting that alterations found in adult subjects arise early on in postnatal development, a critical period of microglia-dependent neural circuit refinement. Our results show that interruption of the Cx3cr1-Cx3cl1/neuron-microglia axis leads to changes in CA3 configuration that affect innate and learned behaviors.

Small-world connectivity dictates collective endothelial cell signaling

Every blood vessel is lined by a single layer of highly specialized, yet adaptable and multifunctional endothelial cells. These cells, the endothelium, control vascular contractility, hemostasis, and inflammation and regulate the exchange of oxygen, nutrients, and waste products between circulating blood and tissue. To control each function, the endothelium processes endlessly arriving requests from multiple sources using separate clusters of cells specialized to detect specific stimuli. A well-developed but poorly understood communication system operates between cells to integrate multiple lines of information and coordinate endothelial responses. Here, the nature of the communication network has been addressed using single-cell Ca2+ imaging across thousands of endothelial cells in intact blood vessels. Cell activities were cross-correlated and compared to a stochastic model to determine network connections. Highly correlated Ca2+ activities occurred in scattered cell clusters, and network communication links between them exhibited unexpectedly short path lengths. The number of connections between cells (degree distribution) followed a power-law relationship revealing a scale-free network topology. The path length and degree distribution revealed an endothelial network with a “small-world” configuration. The small-world configuration confers particularly dynamic endothelial properties including high signal-propagation speed, stability, and a high degree of synchronizability. Local activation of small clusters of cells revealed that the short path lengths and rapid signal transmission were achieved by shortcuts via connecting extensions to nonlocal cells. These findings reveal that the endothelial network design is effective for local and global efficiency in the interaction of the cells and rapid and robust communication between endothelial cells in order to efficiently control cardiovascular activity.

A 3D in vitro model of the device-tissue interface: Functional and structural symptoms of innate neuroinflammation are mitigated by antioxidant ceria nanoparticles

OBJECTIVE

The recording instability of neural implants due to neuroinflammation at the device-tissue interface is a primary roadblock to broad adoption of brain-machine interfaces. While a multiphasic immune response, marked by glial scaring, oxidative stress (OS), and neurodegeneration, is well-characterized, the independent contributions of systemic and local "innate" immune responses are not well-understood. We aimed to understand and mitigate the isolated the innate neuroinflammatory response to devices.

APPROACH

Three-dimensional primary neural cultures provide a unique environment for studying the drivers of neuroinflammation by decoupling the innate and systemic immune systems, while conserving an endogenous extracellular matrix and structural and functional network complexity. We created a three-dimensional in vitro model of the DTI by seeding primary cortical cells around microwires. Live imaging of both dye and AAV-mediated functional, structural, and lipid peroxidation fluorescence was employed to characterize the neuroinflammatory response.

MAIN RESULTS

Live imaging of microtissues over time revealed independent innate neuroinflammation, marked by increased OS, decreased neuronal density, and increased functional connectivity. We demonstrated the use of this model for therapeutic screening by directly applying drugs to neural tissue, bypassing low bioavailability through the in vivo blood brain barrier. As there is growing interest in long-acting antioxidant therapies, we tested efficacy of "perpetual" antioxidant ceria nanoparticles, which reduced OS, increased neuronal density, and protected functional connectivity.

SIGNIFICANCE

Our 3D in vitro model of the device-tissue interface exhibited symptoms of OS-mediated innate neuroinflammation, indicating a significant local immune response to devices. The dysregulation of functional connectivity of microcircuits surround implants suggests the presence of an observer effect, in which the process of recording neural activity may fundamentally change the neural signal. Finally, the demonstration of antioxidant ceria nanoparticle treatment exhibited substantial promise as a neuroprotective and anti-inflammatory treatment strategy.

In-vitro engineered human cerebral tissues mimic pathological circuit disturbances in 3D

In-vitro modeling of brain network disorders such as epilepsy remains a major challenge. A critical step is to develop an experimental approach that enables recapitulation of in-vivo-like three-dimensional functional complexity while allowing local modulation of the neuronal networks. Here, by promoting matrix-supported active cell reaggregation, we engineered multiregional cerebral tissues with intact 3D neuronal networks and functional interconnectivity characteristic of brain networks. Furthermore, using a multi-chambered tissue-culture chip, we show that our separated but interconnected cerebral tissues can mimic neuropathological signatures such as the propagation of epileptiform discharges.

A method is developed to engineer cerebral tissues with intact 3D neuronal networks, mimicking neuropathological signatures such as the propagation of epileptiform discharges, using a multi-chambered tissue culture chip.

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.

Common network effect-patterns after monoamine reuptake inhibition in dissociated hippocampus cultures

The pharmacological treatment of major depressive disorder with currently available antidepressant drugs is still unsatisfying as response to medication is delayed and in some patients even non-existent. To understand complex psychiatric diseases such as major depressive disorder and their treatment, research focus is shifting from investigating single neurons towards a view of the entire functional and effective neuronal network, because alterations on single synapses through antidepressant drugs may translate to alterations in the entire network. Here, we examined the effects of monoamine reuptake inhibitors on in vitro hippocampal network dynamics using calcium fluorescence imaging and analyzing the data with means of graph theoretical parameters. Hypothesizing that monoamine reuptake inhibitors operate through changes of effective connectivity on micro-scale neuronal networks, we measured the effects of the selective monoamine reuptake inhibitors GBR-12783, Sertraline, Venlafaxine, and Amitriptyline on neuronal networks. We identified a common pattern of effects of the different tested monoamine reuptake inhibitors. After treatment with GBR-12783, Sertraline, and Venlafaxine, the connectivity degree, measuring the number of existing connections in the network, was significantly decreased. All tested substances led to networks with more submodules and a reduced global efficiency. No monoamine reuptake inhibitor did affect network-wide firing rate, the characteristic path length, or the network strength. In our study, we found that monoamine reuptake inhibition in neuronal networks in vitro results in a sharpening of the network structure. These alterations could be the basis for the reorganization of a large-scale miswired network in major depressive disorder.

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.

Thresholding Functional Connectivity Matrices to Recover the Topological Properties of Large-Scale Neuronal Networks

The identification of the organization principles on the basis of the brain connectivity can be performed in terms of structural (i.e., morphological), functional (i.e., statistical), or effective (i.e., causal) connectivity. If structural connectivity is based on the detection of the morphological (synaptically mediated) links among neurons, functional and effective relationships derive from the recording of the patterns of electrophysiological activity (e.g., spikes, local field potentials). Correlation or information theory-based algorithms are typical routes pursued to find statistical dependencies and to build a functional connectivity matrix. As long as the matrix collects the possible associations among the network nodes, each interaction between the neuron i and j is different from zero, even though there was no morphological, statistical or causal connection between them. Hence, it becomes essential to find and identify only the significant functional connections that are predictive of the structural ones. For this reason, a robust, fast, and automatized procedure should be implemented to discard the “noisy” connections. In this work, we present a Double Threshold (DDT) algorithm based on the definition of two statistical thresholds. The main goal is not to lose weak but significant links, whose arbitrary exclusion could generate functional networks with a too small number of connections and altered topological properties. The algorithm allows overcoming the limits of the simplest threshold-based methods in terms of precision and guaranteeing excellent computational performances compared to shuffling-based approaches. The presented DDT algorithm was compared with other methods proposed in the literature by using a benchmarking procedure based on synthetic data coming from the simulations of large-scale neuronal networks with different structural topologies.

Data-driven computational modeling predicts "superhubs" play key role in epileptic dynamics

How individual neurons influence epileptic networks remains an open question. In this issue of Neuron, Hadjiabadi et al. (2021) use data-driven, computational models to predict the presence of "superhubs": highly connected neurons that drive network activity through feedforward motifs.

WGEVIA: A Graph Level Embedding Method for Microcircuit Data

Functional microcircuits are useful for studying interactions among neural dynamics of neighboring neurons during cognition and emotion. A functional microcircuit is a group of neurons that are spatially close, and that exhibit synchronized neural activities. For computational analysis, functional microcircuits are represented by graphs, which pose special challenges when applied as input to machine learning algorithms. Graph embedding, which involves the conversion of graph data into low dimensional vector spaces, is a general method for addressing these challenges. In this paper, we discuss limitations of conventional graph embedding methods that make them ill-suited to the study of functional microcircuits. We then develop a novel graph embedding framework, called Weighted Graph Embedding with Vertex Identity Awareness (WGEVIA), that overcomes these limitations. Additionally, we introduce a dataset, called the five vertices dataset, that helps in assessing how well graph embedding methods are suited to functional microcircuit analysis. We demonstrate the utility of WGEVIA through extensive experiments involving real and simulated microcircuit data.

Thermoplasmonic neural chip platform for in situ manipulation of neuronal connections in vitro

Cultured neuronal networks with a controlled structure have been widely studied as an in vitro model system to investigate the relationship between network structure and function. However, most cell culture techniques lack the ability to control network structures during cell cultivation, making it difficult to assess functional changes induced by specific structural changes. In this study, we present an in situ manipulation platform based on gold-nanorod-mediated thermoplasmonics to interrogate an in vitro network model. We find that it is possible to induce new neurite outgrowths, eliminate interconnecting neurites, and estimate functional relationships in matured neuronal networks. This method is expected to be useful for studying functional dynamics of neural networks under controlled structural changes.

Cultured neuron networks provide insight into network structure and function, but the ability to control network topology is a challenge. Here the authors develop a nanorod-mediated thermoplasmonics platform that enables the formation of new connections, the abolishment of existing connections, and the modulation of network activity during cultivation.

Spontaneous Activity of Neuronal Ensembles in Mouse Motor Cortex: Changes after GABAergic Blockade

The mouse motor cortex exhibits spontaneous activity in the form of temporal sequences of neuronal ensembles in vitro without the need of tissue stimulation. These neuronal ensembles are defined as groups of neurons with a strong correlation between its firing patterns, generating what appears to be a predetermined neural conduction mode that needs study. Each ensemble is commonly accompanied by one or more parvalbumin expressing neurons (PV+) or fast spiking interneurons. Many of these interneurons have functional connections between them, helping to form a circuit configuration similar to a small-world network. However, rich club metrics show that most connected neurons are neurons not expressing parvalbumin, mainly pyramidal neurons (PV-) suggesting feed-forward propagation through pyramidal cells. Ensembles with PV+ neurons are connected to these hubs. When ligand-gated fast GABAergic transmission is blocked, temporal sequences of ensembles collapse into a unique synchronous and recurrent ensemble, showing the need of inhibition for coding cortical spontaneous activity. This new ensemble has a duration and electrophysiological characteristics of brief recurrent interictal epileptiform discharges (IEDs) composed by the coactivity of both PV- and PV+ neurons, demonstrating that GABA transmission impedes its occurrence. Synchronous ensembles are clearly divided into two clusters one of them lasting longer and mainly composed by PV+ neurons. Because an ictal-like event was not recorded after several minutes of IEDs recording, it is inferred that an external stimulus and/or fast GABA transmission are necessary for its appearance, making this preparation ideal to study both the neuronal machinery to encode cortical spontaneous activity and its transformation into brief non-ictal epileptiform discharges.

Inference of synaptic connectivity and external variability in neural microcircuits

A major goal in neuroscience is to estimate neural connectivity from large scale extracellular recordings of neural activity in vivo. This is challenging in part because any such activity is modulated by the unmeasured external synaptic input to the network, known as the common input problem. Many different measures of functional connectivity have been proposed in the literature, but their direct relationship to synaptic connectivity is often assumed or ignored. For in vivo data, measurements of this relationship would require a knowledge of ground truth connectivity, which is nearly always unavailable. Instead, many studies use in silico simulations as benchmarks for investigation, but such approaches necessarily rely upon a variety of simplifying assumptions about the simulated network and can depend on numerous simulation parameters. We combine neuronal network simulations, mathematical analysis, and calcium imaging data to address the question of when and how functional connectivity, synaptic connectivity, and latent external input variability can be untangled. We show numerically and analytically that, even though the precision matrix of recorded spiking activity does not uniquely determine synaptic connectivity, it is in practice often closely related to synaptic connectivity. This relation becomes more pronounced when the spatial structure of neuronal variability is jointly considered.

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.

Impact of targeted attack on the spontaneous activity in spatial and biologically-inspired neuronal networks

We study the structural and dynamical consequences of damage in spatial neuronal networks. Inspired by real in vitro networks, we construct directed networks embedded in a two-dimensional space and follow biological rules for designing the wiring of the system. As a result, synthetic cultures display strong metric correlations similar to those observed in real experiments. In its turn, neuronal dynamics is incorporated through the Izhikevich model adopting the parameters derived from observation in real cultures. We consider two scenarios for damage, targeted attacks on those neurons with the highest out-degree and random failures. By analyzing the evolution of both the giant connected component and the dynamical patterns of the neurons as nodes are removed, we observe that network activity halts for a removal of 50% of the nodes in targeted attacks, much lower than the 70% node removal required in the case of random failures. Notably, the decrease of neuronal activity is not gradual. Both damage scenarios portray "boosts" of activity just before full silencing that are not present in equivalent random (Erdös-Rényi) graphs. These boosts correspond to small, spatially compact subnetworks that are able to maintain high levels of activity. Since these subnetworks are absent in the equivalent random graphs, we hypothesize that metric correlations facilitate the existence of local circuits sufficiently integrated to maintain activity, shaping an intrinsic mechanism for resilience.

Single-neuron perturbations reveal feature-specific competition in V1

The computations performed by local neural populations, such as a cortical layer, are typically inferred from anatomical connectivity and observations of neural activity. Here we describe a method-influence mapping-that uses single-neuron perturbations to directly measure how cortical neurons reshape sensory representations. In layer 2/3 of the primary visual cortex (V1), we use two-photon optogenetics to trigger action potentials in a targeted neuron and calcium imaging to measure the effect on spiking in neighbouring neurons in awake mice viewing visual stimuli. Excitatory neurons on average suppressed other neurons and had a centre-surround influence profile over anatomical space. A neuron's influence on its neighbour depended on their similarity in activity. Notably, neurons suppressed activity in similarly tuned neurons more than in dissimilarly tuned neurons. In addition, photostimulation reduced the population response, specifically to the targeted neuron's preferred stimulus, by around 2%. Therefore, V1 layer 2/3 performed feature competition, in which a like-suppresses-like motif reduces redundancy in population activity and may assist with inference of the features that underlie sensory input. We anticipate that influence mapping can be extended to investigate computations in other neural populations.

A Proposed Mechanism for Spontaneous Transitions between Interictal and Ictal Activity

Epileptic networks are characterized by two outputs: brief interictal spikes and rarer, more prolonged seizures. Although either output state is readily modeled in silico and induced experimentally, the transition mechanisms are unknown, in part because no models exhibit both output states spontaneously. In silico small-world neural networks were built using single-compartment neurons whose physiological parameters were derived from dual whole-cell recordings of pyramidal cells in organotypic hippocampal slice cultures that were generating spontaneous seizure-like activity. In silico, neurons were connected by abundant local synapses and rare long-distance synapses. Activity-dependent synaptic depression and gradual recovery delimited synchronous activity. Full synaptic recovery engendered interictal population spikes that spread via long-distance synapses. When synaptic recovery was incomplete, postsynaptic neurons required coincident activation of multiple presynaptic terminals to reach firing threshold. Only local connections were sufficiently dense to spread activity under these conditions. This coalesced network activity into traveling waves whose velocity varied with synaptic recovery. Seizures were comprised of sustained traveling waves that were similar to those recorded during experimental and human neocortical seizures. Sustained traveling waves occurred only when wave velocity, network dimensions, and the rate of synaptic recovery enabled wave reentry into previously depressed areas at precisely ictogenic levels of synaptic recovery. Wide-field, cellular-resolution GCamP7b calcium imaging demonstrated similar initial patterns of activation in the hippocampus, although the anatomical distribution of traveling waves of synaptic activation was altered by the pattern of synaptic connectivity in the organotypic hippocampal cultures.SIGNIFICANCE STATEMENT When computerized distributed neural network models are required to generate both features of epileptic networks (i.e., spontaneous interictal population spikes and seizures), the network structure is substantially constrained. These constraints provide important new hypotheses regarding the nature of epileptic networks and mechanisms of seizure onset.

Reconstruction of Functional Connectivity from Multielectrode Recordings and Calcium Imaging

In the last two decades, increasing research efforts in neuroscience have been focused on determining both structural and functional connectivity of brain circuits, with the main goal of relating the wiring diagram of neuronal systems to their emerging properties, from the microscale to the macroscale. While combining multisite parallel recordings with structural circuits' reconstruction in vivo is still very challenging, the reductionist in vitro approach based on neuronal cultures offers lower technical difficulties and is much more stable under control conditions. In this chapter, we present different approaches to infer the connectivity of cultured neuronal networks using multielectrode array or calcium imaging recordings. We first formally introduce the used methods, and then we will describe into details how those methods were applied in case studies. Since multielectrode array and calcium imaging recordings provide distinct and complementary spatiotemporal features of neuronal activity, in this chapter we present the strategies implemented with the two different methodologies in distinct sections.

Data-driven brain network models differentiate variability across language tasks

The relationship between brain structure and function has been probed using a variety of approaches, but how the underlying structural connectivity of the human brain drives behavior is far from understood. To investigate the effect of anatomical brain organization on human task performance, we use a data-driven computational modeling approach and explore the functional effects of naturally occurring structural differences in brain networks. We construct personalized brain network models by combining anatomical connectivity estimated from diffusion spectrum imaging of individual subjects with a nonlinear model of brain dynamics. By performing computational experiments in which we measure the excitability of the global brain network and spread of synchronization following a targeted computational stimulation, we quantify how individual variation in the underlying connectivity impacts both local and global brain dynamics. We further relate the computational results to individual variability in the subjects' performance of three language-demanding tasks both before and after transcranial magnetic stimulation to the left-inferior frontal gyrus. Our results show that task performance correlates with either local or global measures of functional activity, depending on the complexity of the task. By emphasizing differences in the underlying structural connectivity, our model serves as a powerful tool to assess individual differences in task performances, to dissociate the effect of targeted stimulation in tasks that differ in cognitive demand, and to pave the way for the development of personalized therapeutics.

On the nature and use of models in network neuroscience

Network theory provides an intuitively appealing framework for studying relationships among interconnected brain mechanisms and their relevance to behaviour. As the space of its applications grows, so does the diversity of meanings of the term network model. This diversity can cause confusion, complicate efforts to assess model validity and efficacy, and hamper interdisciplinary collaboration. In this Review, we examine the field of network neuroscience, focusing on organizing principles that can help overcome these challenges. First, we describe the fundamental goals in constructing network models. Second, we review the most common forms of network models, which can be described parsimoniously along the following three primary dimensions: from data representations to first-principles theory; from biophysical realism to functional phenomenology; and from elementary descriptions to coarse-grained approximations. Third, we draw on biology, philosophy and other disciplines to establish validation principles for these models. We close with a discussion of opportunities to bridge model types and point to exciting frontiers for future pursuits.

Identification of excitatory-inhibitory links and network topology in large-scale neuronal assemblies from multi-electrode recordings

Functional-effective connectivity and network topology are nowadays key issues for studying brain physiological functions and pathologies. Inferring neuronal connectivity from electrophysiological recordings presents open challenges and unsolved problems. In this work, we present a cross-correlation based method for reliably estimating not only excitatory but also inhibitory links, by analyzing multi-unit spike activity from large-scale neuronal networks. The method is validated by means of realistic simulations of large-scale neuronal populations. New results related to functional connectivity estimation and network topology identification obtained by experimental electrophysiological recordings from high-density and large-scale (i.e., 4096 electrodes) microtransducer arrays coupled to in vitro neural populations are presented. Specifically, we show that: (i) functional inhibitory connections are accurately identified in in vitro cortical networks, providing that a reasonable firing rate and recording length are achieved; (ii) small-world topology, with scale-free and rich-club features are reliably obtained, on condition that a minimum number of active recording sites are available. The method and procedure can be directly extended and applied to in vivo multi-units brain activity recordings.

The balance between excitation and inhibition is fundamental for proper brain functions and for this reason is precisely regulated in adult cortices. Impaired excitation/inhibition balance is often associated with several neurological disorders, such as epilepsy, autism and schizophrenia. However, estimating functional inhibitory connections is not an easy task and few methods are available to identify such connections from electrophysiological data. Here we present a cross-correlation based method to identify both excitatory and inhibitory functional connections in large-scale neuronal networks. The method is applicable to both in vitro and in vivo spike data recordings. Once a connectivity map (i.e. a graph) is obtained, we characterized the associated topology by means of classical graph theory metrics to unveil functional architecture. In this work, we analyze in vitro cortical networks probed by means of large-scale microelectrode arrays (i.e., 4096 sensors) and we derive network topologies from spike data. The functional organization found is called “small-world and scale-free” and is the same organization found in cortical in vivo brain regions by means of different experimental methods. We also show that to obtain reliable information about network architecture at least a network with a hundred of nodes-neurons is needed.

Link Prediction Investigation of Dynamic Information Flow in Epilepsy

As a brain disorder, epilepsy is characterized with abnormal hypersynchronous neural firings. It is known that seizures initiate and propagate in different brain regions. Long-term intracranial multichannel electroencephalography (EEG) reflects broadband ictal activity under seizure occurrence. Network-based techniques are efficient in discovering brain dynamics and offering finger-print features for specific individuals. In this study, we adopt link prediction for proposing a novel workflow aiming to quantify seizure dynamics and uncover pathological mechanisms of epilepsy. A dataset of EEG signals was enrolled that recorded from 8 patients with 3 different types of pharmocoresistant focal epilepsy. Weighted networks are obtained from phase locking value (PLV) in subband EEG oscillations. Common neighbor (CN), resource allocation (RA), Adamic-Adar (AA), and Sorenson algorithms are brought in for link prediction performance comparison. Results demonstrate that RA outperforms its rivals. Similarity, matrix was produced from the RA technique performing on EEG networks later. Nodes are gathered to form sequences by selecting the ones with the highest similarity. It is demonstrated that variations are in accordance with seizure attack in node sequences of gamma band EEG oscillations. What is more, variations in node sequences monitor the total seizure journey including its initiation and termination.

Structure-function multi-scale connectomics reveals a major role of the fronto-striato-thalamic circuit in brain aging

Physiological aging affects brain structure and function impacting morphology, connectivity, and performance. However, whether some brain connectivity metrics might reflect the age of an individual is still unclear. Here, we collected brain images from healthy participants (N = 155) ranging from 10 to 80 years to build functional (resting state) and structural (tractography) connectivity matrices, both data sets combined to obtain different connectivity features. We then calculated the brain connectome age-an age estimator resulting from a multi-scale methodology applied to the structure-function connectome, and compared it to the chronological age (ChA). Our results were twofold. First, we found that aging widely affects the connectivity of multiple structures, such as anterior cingulate and medial prefrontal cortices, basal ganglia, thalamus, insula, cingulum, hippocampus, parahippocampus, occipital cortex, fusiform, precuneus, and temporal pole. Second, we found that the connectivity between basal ganglia and thalamus to frontal areas, also known as the fronto-striato-thalamic (FST) circuit, makes the major contribution to age estimation. In conclusion, our results highlight the key role played by the FST circuit in the process of healthy aging. Notably, the same methodology can be generally applied to identify the structural-functional connectivity patterns correlating to other biomarkers than ChA.

Personalized brain network models for assessing structure-function relationships

Many recent efforts in computational modeling of macro-scale brain dynamics have begun to take a data-driven approach by incorporating structural and/or functional information derived from subject data. Here, we discuss recent work using personalized brain network models to study structure-function relationships in human brains. We describe the steps necessary to build such models and show how this computational approach can provide previously unobtainable information through the ability to perform virtual experiments. Finally, we present examples of how personalized brain network models can be used to gain insight into the effects of local stimulation and improve surgical outcomes in epilepsy.

Effective Subnetwork Topology for Synchronizing Interconnected Networks of Coupled Phase Oscillators

A system consisting of interconnected networks, or a network of networks (NoN), appears diversely in many real-world systems, including the brain. In this study, we consider NoNs consisting of heterogeneous phase oscillators and investigate how the topology of subnetworks affects the global synchrony of the network. The degree of synchrony and the effect of subnetwork topology are evaluated based on the Kuramoto order parameter and the minimum coupling strength necessary for the order parameter to exceed a threshold value, respectively. In contrast to an isolated network in which random connectivity is favorable for achieving synchrony, NoNs synchronize with weaker interconnections when the degree distribution of subnetworks is heterogeneous, suggesting the major role of the high-degree nodes. We also investigate a case in which subnetworks with different average natural frequencies are coupled to show that direct coupling of subnetworks with the largest variation is effective for synchronizing the whole system. In real-world NoNs like the brain, the balance of synchrony and asynchrony is critical for its function at various spatial resolutions. Our work provides novel insights into the topological basis of coordinated dynamics in such networks.

Functional network stability and average minimal distance - A framework to rapidly assess dynamics of functional network representations

BACKGROUND

Recent advances in neurophysiological recording techniques have increased both the spatial and temporal resolution of data. New methodologies are required that can handle large data sets in an efficient manner as well as to make quantifiable, and realistic, predictions about the global modality of the brain from under-sampled recordings.

NEW METHOD

To rectify both problems, we first propose an analytical modification to an existing functional connectivity algorithm, Average Minimal Distance (AMD), to rapidly capture functional network connectivity. We then complement this algorithm by introducing Functional Network Stability (FuNS), a metric that can be used to quickly assess the global network dynamic changes over time, without being constrained by the activities of a specific set of neurons.

RESULTS

We systematically test the performance of AMD and FuNS (1) on artificial spiking data with different statistical characteristics, (2) from spiking data generated using a neural network model, and (3) using in vivo data recorded from mouse hippocampus during fear learning. Our results show that AMD and FuNS are able to monitor the change in network dynamics during memory consolidation.

COMPARISON WITH OTHER METHODS

AMD outperforms traditional bootstrapping and cross-correlation (CC) methods in both significance and computation time. Simultaneously, FuNS provides a reliable way to establish a link between local structural network changes, global dynamics of network-wide representations activity, and behavior.

CONCLUSIONS

The AMD-FuNS framework should be universally useful in linking long time-scale, global network dynamics and cognitive behavior.

Generative models for network neuroscience: prospects and promise

Network neuroscience is the emerging discipline concerned with investigating the complex patterns of interconnections found in neural systems, and identifying principles with which to understand them. Within this discipline, one particularly powerful approach is network generative modelling, in which wiring rules are algorithmically implemented to produce synthetic network architectures with the same properties as observed in empirical network data. Successful models can highlight the principles by which a network is organized and potentially uncover the mechanisms by which it grows and develops. Here, we review the prospects and promise of generative models for network neuroscience. We begin with a primer on network generative models, with a discussion of compressibility and predictability, and utility in intuiting mechanisms, followed by a short history on their use in network science, broadly. We then discuss generative models in practice and application, paying particular attention to the critical need for cross-validation. Next, we review generative models of biological neural networks, both at the cellular and large-scale level, and across a variety of species including Caenorhabditis elegans, Drosophila, mouse, rat, cat, macaque and human. We offer a careful treatment of a few relevant distinctions, including differences between generative models and null models, sufficiency and redundancy, inferring and claiming mechanism, and functional and structural connectivity. We close with a discussion of future directions, outlining exciting frontiers both in empirical data collection efforts as well as in method and theory development that, together, further the utility of the generative network modelling approach for network neuroscience.

Inferring structural connectivity using Ising couplings in models of neuronal networks

Functional connectivity metrics have been widely used to infer the underlying structural connectivity in neuronal networks. Maximum entropy based Ising models have been suggested to discount the effect of indirect interactions and give good results in inferring the true anatomical connections. However, no benchmarking is currently available to assess the performance of Ising couplings against other functional connectivity metrics in the microscopic scale of neuronal networks through a wide set of network conditions and network structures. In this paper, we study the performance of the Ising model couplings to infer the synaptic connectivity in in silico networks of neurons and compare its performance against partial and cross-correlations for different correlation levels, firing rates, network sizes, network densities, and topologies. Our results show that the relative performance amongst the three functional connectivity metrics depends primarily on the network correlation levels. Ising couplings detected the most structural links at very weak network correlation levels, and partial correlations outperformed Ising couplings and cross-correlations at strong correlation levels. The result was consistent across varying firing rates, network sizes, and topologies. The findings of this paper serve as a guide in choosing the right functional connectivity tool to reconstruct the structural connectivity.

Recurrently connected and localized neuronal communities initiate coordinated spontaneous activity in neuronal networks

Developing neuronal systems intrinsically generate coordinated spontaneous activity that propagates by involving a large number of synchronously firing neurons. In vivo, waves of spikes transiently characterize the activity of developing brain circuits and are fundamental for activity-dependent circuit formation. In vitro, coordinated spontaneous spiking activity, or network bursts (NBs), interleaved within periods of asynchronous spikes emerge during the development of 2D and 3D neuronal cultures. Several studies have investigated this type of activity and its dynamics, but how a neuronal system generates these coordinated events remains unclear. Here, we investigate at a cellular level the generation of network bursts in spontaneously active neuronal cultures by exploiting high-resolution multielectrode array recordings and computational network modelling. Our analysis reveals that NBs are generated in specialized regions of the network (functional neuronal communities) that feature neuronal links with high cross-correlation peak values, sub-millisecond lags and that share very similar structural connectivity motifs providing recurrent interactions. We show that the particular properties of these local structures enable locally amplifying spontaneous asynchronous spikes and that this mechanism can lead to the initiation of NBs. Through the analysis of simulated and experimental data, we also show that AMPA currents drive the coordinated activity, while NMDA and GABA currents are only involved in shaping the dynamics of NBs. Overall, our results suggest that the presence of functional neuronal communities with recurrent local connections allows a neuronal system to generate spontaneous coordinated spiking activity events. As suggested by the rules used for implementing our computational model, such functional communities might naturally emerge during network development by following simple constraints on distance-based connectivity.

Coordinated spontaneous spiking activity is fundamental for the normal formation of brain circuits during development. However, how ensembles of neurons generate these events remains unclear. To address this question, in the present study, we investigated the network properties that might be required to a neuronal system for the generation of these spontaneous waves of spikes. We performed our study on spontaneously active neuronal cell cultures using high-resolution electrical recordings and a computational network model developed to reproduce our experimental data both quantitatively and qualitatively. Through the analysis of both experimental and simulated data, we found that network bursts are initiated in regions of the network, or “functional communities”, characterized by particular local connectivity properties. We also found that these regions can amplify the background asynchronous spiking activity preceding a network burst and, in this way, can give rise to coordinated spiking events.

As a whole, our results suggest the presence of functional communities of neurons in a developing neuronal system that might naturally emerge by following simple constraints on distance-based connectivity. These regions are most likely required for the generation of the spontaneous coordinated activity that can drive activity-dependent circuit formation.

Cortical Dynamics in Presence of Assemblies of Densely Connected Weight-Hub Neurons

Experimental measurements of pairwise connection probability of pyramidal neurons together with the distribution of synaptic weights have been used to construct randomly connected model networks. However, several experimental studies suggest that both wiring and synaptic weight structure between neurons show statistics that differ from random networks. Here we study a network containing a subset of neurons which we call weight-hub neurons, that are characterized by strong inward synapses. We propose a connectivity structure for excitatory neurons that contain assemblies of densely connected weight-hub neurons, while the pairwise connection probability and synaptic weight distribution remain consistent with experimental data. Simulations of such a network with generalized integrate-and-fire neurons display regular and irregular slow oscillations akin to experimentally observed up/down state transitions in the activity of cortical neurons with a broad distribution of pairwise spike correlations. Moreover, stimulation of a model network in the presence or absence of assembly structure exhibits responses similar to light-evoked responses of cortical layers in optogenetically modified animals. We conclude that a high connection probability into and within assemblies of excitatory weight-hub neurons, as it likely is present in some but not all cortical layers, changes the dynamics of a layer of cortical microcircuitry significantly.

Spontaneous dynamics of neural networks in deep layers of prefrontal cortex

Cortical systems maintain and process information through the sustained activation of recurrent local networks of neurons. Layer 5 is known to have a major role in generating the recurrent activation associated with these functions, but relatively little is known about its intrinsic dynamics at the mesoscopic level of large numbers of neighboring neurons. Using calcium imaging, we measured the spontaneous activity of networks of deep-layer medial prefrontal cortical neurons in an acute slice model. Inferring the simultaneous activity of tens of neighboring neurons, we found that while the majority showed only sporadic activity, a subset of neurons engaged in sustained delta frequency rhythmic activity. Spontaneous activity under baseline conditions was weakly correlated between pairs of neurons, and rhythmic neurons showed little coherence in their oscillations. However, we consistently observed brief bouts of highly synchronous activity that must be attributed to network activity. NMDA-mediated stimulation enhanced rhythmicity, synchrony, and correlation within these local networks. These results characterize spontaneous prefrontal activity at a previously unexplored spatiotemporal scale and suggest that medial prefrontal cortex can act as an intrinsic generator of delta oscillations.NEW & NOTEWORTHY Using calcium imaging and a novel analytic framework, we characterized the spontaneous and NMDA-evoked activity of layer 5 prefrontal cortex at a largely unexplored spatiotemporal scale. Our results suggest that the mPFC microcircuitry is capable of intrinsically generating delta oscillations and sustaining synchronized network activity that is potentially relevant for understanding its contribution to cognitive processes.

Two's company, three (or more) is a simplex : Algebraic-topological tools for understanding higher-order structure in neural data

The language of graph theory, or network science, has proven to be an exceptional tool for addressing myriad problems in neuroscience. Yet, the use of networks is predicated on a critical simplifying assumption: that the quintessential unit of interest in a brain is a dyad - two nodes (neurons or brain regions) connected by an edge. While rarely mentioned, this fundamental assumption inherently limits the types of neural structure and function that graphs can be used to model. Here, we describe a generalization of graphs that overcomes these limitations, thereby offering a broad range of new possibilities in terms of modeling and measuring neural phenomena. Specifically, we explore the use of simplicial complexes: a structure developed in the field of mathematics known as algebraic topology, of increasing applicability to real data due to a rapidly growing computational toolset. We review the underlying mathematical formalism as well as the budding literature applying simplicial complexes to neural data, from electrophysiological recordings in animal models to hemodynamic fluctuations in humans. Based on the exceptional flexibility of the tools and recent ground-breaking insights into neural function, we posit that this framework has the potential to eclipse graph theory in unraveling the fundamental mysteries of cognition.

Signal Propagation between Neuronal Populations Controlled by Micropatterning

The central nervous system consists of an unfathomable number of functional networks enabling highly sophisticated information processing. Guided neuronal growth with a well-defined connectivity and accompanying polarity is essential for the formation of these networks. To investigate how two-dimensional protein patterns influence neuronal outgrowth with respect to connectivity and functional polarity between adjacent populations of neurons, a microstructured model system was established. Exclusive cell growth on patterned substrates was achieved by transferring a mixture of poly-l-lysine and laminin to a cell-repellent glass surface by microcontact printing. Triangular structures with different opening angle, height, and width were chosen as a pattern to achieve network formation with defined behavior at the junction of adjacent structures. These patterns were populated with dissociated primary cortical embryonic rat neurons and investigated with respect to their impact on neuronal outgrowth by immunofluorescence analysis, as well as their functional connectivity by calcium imaging. Here, we present a highly reproducible technique to devise neuronal networks in vitro with a predefined connectivity induced by the design of the gateway. Daisy-chained neuronal networks with predefined connectivity and functional polarity were produced using the presented micropatterning method. Controlling the direction of signal propagation among populations of neurons provides insights to network communication and offers the chance to investigate more about learning processes in networks by external manipulation of cells and signal cascades.

High-Degree Neurons Feed Cortical Computations

Recent work has shown that functional connectivity among cortical neurons is highly varied, with a small percentage of neurons having many more connections than others. Also, recent theoretical developments now make it possible to quantify how neurons modify information from the connections they receive. Therefore, it is now possible to investigate how information modification, or computation, depends on the number of connections a neuron receives (in-degree) or sends out (out-degree). To do this, we recorded the simultaneous spiking activity of hundreds of neurons in cortico-hippocampal slice cultures using a high-density 512-electrode array. This preparation and recording method combination produced large numbers of neurons recorded at temporal and spatial resolutions that are not currently available in any in vivo recording system. We utilized transfer entropy (a well-established method for detecting linear and nonlinear interactions in time series) and the partial information decomposition (a powerful, recently developed tool for dissecting multivariate information processing into distinct parts) to quantify computation between neurons where information flows converged. We found that computations did not occur equally in all neurons throughout the networks. Surprisingly, neurons that computed large amounts of information tended to receive connections from high out-degree neurons. However, the in-degree of a neuron was not related to the amount of information it computed. To gain insight into these findings, we developed a simple feedforward network model. We found that a degree-modified Hebbian wiring rule best reproduced the pattern of computation and degree correlation results seen in the real data. Interestingly, this rule also maximized signal propagation in the presence of network-wide correlations, suggesting a mechanism by which cortex could deal with common random background input. These are the first results to show that the extent to which a neuron modifies incoming information streams depends on its topological location in the surrounding functional network.

We recorded the electrical activity of hundreds of neurons simultaneously in brain tissue from mice and we analyzed these signals using state-of-the-art tools from information theory. These tools allowed us to ascertain which neurons were transmitting information to other neurons and to characterize the computations performed by neurons using the inputs they received from two or more other neurons. We found that computations did not occur equally in all neurons throughout the networks. Surprisingly, neurons that computed large amounts of information tended to be recipients of information from neurons with a large number of outgoing connections. Interestingly, the number of incoming connections to a neuron was not related to the amount of information that neuron computed. To better understand these results, we built a network model to match the data. Unexpectedly, the model also maximized information transfer in the presence of network-wide correlations. This suggested a way that networks of cortical neurons could deal with common random background input. These results are the first to show that the amount of information computed by a neuron depends on where it is located in the surrounding network.