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.

DeePhys: A machine learning-assisted platform for electrophysiological phenotyping of human neuronal networks

Reproducible functional assays to study in vitro neuronal networks represent an important cornerstone in the quest to develop physiologically relevant cellular models of human diseases. Here, we introduce DeePhys, a MATLAB-based analysis tool for data-driven functional phenotyping of in vitro neuronal cultures recorded by high-density microelectrode arrays. DeePhys is a modular workflow that offers a range of techniques to extract features from spike-sorted data, allowing for the examination of functional phenotypes both at the individual cell and network levels, as well as across development. In addition, DeePhys incorporates the capability to integrate novel features and to use machine-learning-assisted approaches, which facilitates a comprehensive evaluation of pharmacological interventions. To illustrate its practical application, we apply DeePhys to human induced pluripotent stem cell-derived dopaminergic neurons obtained from both patients and healthy individuals and showcase how DeePhys enables phenotypic screenings.

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.

Functional imaging of brain organoids using high-density microelectrode arrays

ABSTRACT

Studies have provided evidence that human cerebral organoids (hCOs) recapitulate fundamental milestones of early brain development, but many important questions regarding their functionality and electrophysiological properties persist. High-density microelectrode arrays (HD-MEAs) represent an attractive analysis platform to perform functional studies of neuronal networks at the cellular and network scale. Here, we use HD-MEAs to derive large-scale electrophysiological recordings from sliced hCOs. We record the activity of hCO slices over several weeks and probe observed neuronal dynamics pharmacologically. Moreover, we present results on how the obtained recordings can be spike-sorted and subsequently studied across scales. For example, we show how to track single neurons across several days on the HD-MEA and how to infer axonal action potential velocities. We also infer putative functional connectivity from hCO recordings. The introduced methodology will contribute to a better understanding of developing neuronal networks in brain organoids and provide new means for their functional characterization.

IMPACT STATEMENT

Human cerebral organoids (hCOs) represent an attractive in vitro model system to study key physiological mechanisms underlying early neuronal network formation in tissue with healthy or disease-related genetic backgrounds. Despite remarkable advances in the generation of brain organoids, knowledge on the functionality of their neuronal circuits is still scarce. Here, we used complementary metal-oxide-semiconductor (CMOS)-based high-density microelectrode arrays (HD-MEAs) to perform large-scale recordings from sliced hCOs over several weeks and quantified their activity across scales. Using single-cell and network metrics, we were able to probe aspects of hCO neurophysiology that are more difficult to obtain with other techniques, such as patch clamping (lower yield) and calcium imaging (lower temporal resolution). These metrics included, for example, extracellular action potential (AP) waveform features and axonal AP velocity at the cellular level, as well as functional connectivity at the network level. Analysis was enabled by the large sensing area and the high spatiotemporal resolution provided by HD-MEAs, which allowed recordings from hundreds of neurons and spike sorting of their activity. Our results demonstrate that HD-MEAs provide a multi-purpose platform for the functional characterization of hCOs, which will be key in improving our understanding of this model system and assessing its relevance for translational research.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1557/s43577-022-00282-w.