Evolving Simple Models of Diverse Intrinsic Dynamics in Hippocampal Neuron Types

Mathematical generation of data-driven hippocampal CA1 pyramidal neurons and interneurons copies via A-GLIF models for large-scale networks covering the experimental variability range

Efficient and accurate large-scale networks are a fundamental tool in modeling brain areas, to advance our understanding of neuronal dynamics. However, their implementation faces two key issues: computational efficiency and heterogeneity. Computational efficiency is achieved using simplified neurons, whereas there are no practical solutions available to solve the problem of reproducing in a large-scale network the experimentally observed heterogeneity of the intrinsic properties of neurons. This is important, because the use of identical nodes in a network can generate artifacts which can hinder an adequate representation of the properties of a real network. To this aim, we introduce a mathematical procedure to generate an arbitrary large number of copies of simplified hippocampal CA1 pyramidal neurons and interneurons models, which exhibit the full range of firing dynamics observed in these cells - including adapting, non-adapting and bursting. For this purpose, we rely on a recently published adaptive generalized leaky integrate-and-fire (A-GLIF) modeling approach, leveraging on its ability to reproduce the rich set of electrophysiological behaviors of these types of neurons under a variety of different stimulation currents. The generation procedure is based on a perturbation of model's parameters related to the initial data, firing block, and internal dynamics, and suitably validated against experimental data to ensure that the firing dynamics of any given cell copy remains within the experimental range. A classification procedure confirmed that the firing behavior of most of the pyramidal/interneuron copies was consistent with the experimental data. This approach allows to obtain heterogeneous copies with mathematically controlled firing properties. A full set of heterogeneous neurons composing the CA1 region of a rat hippocampus (approximately 1.2 million neurons), are provided in a database freely available in the live paper section of the EBRAINS platform. By adapting the underlying A-GLIF framework, it will be possible to extend the numerical approach presented here to create, in a mathematically controlled manner, an arbitrarily large number of non-identical copies of cell populations with firing properties related to other brain areas.

Brain-Region-Specific Genes Form the Major Pathways Featuring Their Basic Functional Role: Their Implication in Animal Chronic Stress Model

The analysis of RNA-Sec data from murine bulk tissue samples taken from five brain regions associated with behavior and stress response was conducted. The focus was on the most contrasting brain region-specific genes (BRSG) sets in terms of their expression rates. These BRSGs are identified as genes with a distinct outlying (high) expression rate in a specific region compared to others used in the study. The analysis suggested that BRSG sets form non-randomly connected compact gene networks, which correspond to the major neuron-mediated functional processes or pathways in each brain region. The number of BRSGs and the connection rate were found to depend on the heterogeneity and coordinated firing rate of neuron types in each brain region. The most connected pathways, along with the highest BRSG number, were observed in the Striatum, referred to as Medium Spiny Neurons (MSNs), which make up 95% of neurons and exhibit synchronous firing upon dopamine influx. However, the Ventral Tegmental Area/Medial Raphe Nucleus (VTA/MRN) regions, although primarily composed of monoaminergic neurons, do not fire synchronously, leading to a smaller BRSG number. The Hippocampus (HPC) region, on the other hand, displays significant neuronal heterogeneity, with glutamatergic neurons being the most numerous and synchronized. Interestingly, the two monoaminergic regions involved in the study displayed a common BRSG subnetwork architecture, emphasizing their proximity in terms of axonal throughput specifics and high-energy metabolism rates. This finding suggests the concerted evolution of monoaminergic neurons, leading to unique adaptations at the genic repertoire scale. With BRSG sets, we were able to highlight the contrasting features of the three groups: control, depressive, and aggressive mice in the animal chronic stress model. Specifically, we observed a decrease in serotonergic turnover in both the depressed and aggressive groups, while dopaminergic emission was high in both groups. There was also a notable absence of dopaminoceptive receptors on the postsynaptic membranes in the striatum in the depressed group. Additionally, we confirmed that neurogenesis BRSGs are specific to HPC, with the aggressive group showing attenuated neurogenesis rates compared to the control/depressive groups. We also confirmed that immune-competent cells like microglia and astrocytes play a crucial role in depressed phenotypes, including mitophagy-related gene Prkcd. Based on this analysis, we propose the use of BRSG sets as a suitable framework for evaluating case-control group-wise assessments of specific brain region gene pathway responses.

Hippocampome.org 2.0 is a knowledge base enabling data-driven spiking neural network simulations of rodent hippocampal circuits

Hippocampome.org is a mature open-access knowledge base of the rodent hippocampal formation focusing on neuron types and their properties. Previously, Hippocampome.org v1.0 established a foundational classification system identifying 122 hippocampal neuron types based on their axonal and dendritic morphologies, main neurotransmitter, membrane biophysics, and molecular expression (Wheeler et al., 2015). Releases v1.1 through v1.12 furthered the aggregation of literature-mined data, including among others neuron counts, spiking patterns, synaptic physiology, in vivo firing phases, and connection probabilities. Those additional properties increased the online information content of this public resource over 100-fold, enabling numerous independent discoveries by the scientific community. Hippocampome.org v2.0, introduced here, besides incorporating over 50 new neuron types, now recenters its focus on extending the functionality to build real-scale, biologically detailed, data-driven computational simulations. In all cases, the freely downloadable model parameters are directly linked to the specific peer-reviewed empirical evidence from which they were derived. Possible research applications include quantitative, multiscale analyses of circuit connectivity and spiking neural network simulations of activity dynamics. These advances can help generate precise, experimentally testable hypotheses and shed light on the neural mechanisms underlying associative memory and spatial navigation.

Hippocampome.org v2.0: a knowledge base enabling data-driven spiking neural network simulations of rodent hippocampal circuits

Hippocampome.org is a mature open-access knowledge base of the rodent hippocampal formation focusing on neuron types and their properties. Hippocampome.org v1.0 established a foundational classification system identifying 122 hippocampal neuron types based on their axonal and dendritic morphologies, main neurotransmitter, membrane biophysics, and molecular expression. Releases v1.1 through v1.12 furthered the aggregation of literature-mined data, including among others neuron counts, spiking patterns, synaptic physiology, in vivo firing phases, and connection probabilities. Those additional properties increased the online information content of this public resource over 100-fold, enabling numerous independent discoveries by the scientific community. Hippocampome.org v2.0, introduced here, besides incorporating over 50 new neuron types, now recenters its focus on extending the functionality to build real-scale, biologically detailed, data-driven computational simulations. In all cases, the freely downloadable model parameters are directly linked to the specific peer-reviewed empirical evidence from which they were derived. Possible research applications include quantitative, multiscale analyses of circuit connectivity and spiking neural network simulations of activity dynamics. These advances can help generate precise, experimentally testable hypotheses and shed light on the neural mechanisms underlying associative memory and spatial navigation.

Robust Resting-State Dynamics in a Large-Scale Spiking Neural Network Model of Area CA3 in the Mouse Hippocampus

Hippocampal area CA3 performs the critical auto-associative function underlying pattern completion in episodic memory. Without external inputs, the electrical activity of this neural circuit reflects the spontaneous spiking interplay among glutamatergic pyramidal neurons and GABAergic interneurons. However, the network mechanisms underlying these resting-state firing patterns are poorly understood. Leveraging the Hippocampome.org knowledge base, we developed a data-driven, large-scale spiking neural network (SNN) model of mouse CA3 with 8 neuron types, 90,000 neurons, 51 neuron-type specific connections, and 250,000,000 synapses. We instantiated the SNN in the CARLsim4 multi-GPU simulation environment using the Izhikevich and Tsodyks-Markram formalisms for neuronal and synaptic dynamics, respectively. We analyzed the resultant population activity upon transient activation. The SNN settled into stable oscillations with a biologically plausible grand-average firing frequency, which was robust relative to a wide range of transient activation. The diverse firing patterns of individual neuron types were consistent with existing knowledge of cell type-specific activity in vivo. Altered network structures that lacked neuron- or connection-type specificity were neither stable nor robust, highlighting the importance of neuron type circuitry. Additionally, external inputs reflecting dentate mossy fibers shifted the observed rhythms to the gamma band. We freely released the CARLsim4-Hippocampome framework on GitHub to test hippocampal hypotheses. Our SNN may be useful to investigate the circuit mechanisms underlying the computational functions of CA3. Moreover, our approach can be scaled to the whole hippocampal formation, which may contribute to elucidating how the unique neuronal architecture of this system subserves its crucial cognitive roles.

Polynomial, piecewise-Linear, Step (PLS): A Simple, Scalable, and Efficient Framework for Modeling Neurons

Biological neurons can be modeled with different levels of biophysical/biochemical details. The accuracy with which a model reflects the actual physiological processes and ultimately the information function of a neuron, can range from very detailed to a schematic phenomenological representation. This range exists due to the common problem: one needs to find an optimal trade-off between the level of details needed to capture the necessary information processing in a neuron and the computational load needed to compute 1 s of model time. An increase in modeled network size or model-time, for which the solution should be obtained, makes this trade-off pivotal in model development. Numerical simulations become incredibly challenging when an extensive network with a detailed representation of each neuron needs to be modeled over a long time interval to study slow evolving processes, e.g., development of the thalamocortical circuits. Here we suggest a simple, powerful and flexible approach in which we approximate the right-hand sides of differential equations by combinations of functions from three families: Polynomial, piecewise-Linear, Step (PLS). To obtain a single coherent framework, we provide four core principles in which PLS functions should be combined. We show the rationale behind each of the core principles. Two examples illustrate how to build a conductance-based or phenomenological model using the PLS-framework. We use the first example as a benchmark on three different computational platforms: CPU, GPU, and mobile system-on-chip devices. We show that the PLS-framework speeds up computations without increasing the memory footprint and maintains high model fidelity comparable to the fully-computed model or with lookup-table approximation. We are convinced that the full range of neuron models: from biophysical to phenomenological and even to abstract models, may benefit from using the PLS-framework.

On the Modeling of the Three Types of Non-spiking Neurons of the Caenorhabditis elegans

The nematode Caenorhabditis elegans (C. elegans) is a well-known model organism in neuroscience. The relative simplicity of its nervous system, made up of few hundred neurons, shares some essential features with more sophisticated nervous systems, including the human one. If we are able to fully characterize the nervous system of this organism, we will be one step closer to understanding the mechanisms underlying the behavior of living things. Following a recently conducted electrophysiological survey on different C. elegans neurons, this paper aims at modeling the three non-spiking RIM, AIY and AFD neurons (arbitrarily named with three upper case letters by convention). To date, they represent the three possible forms of non-spiking neuronal responses of the C. elegans. To achieve this objective, we propose a conductance-based neuron model adapted to the electrophysiological features of each neuron. These features are based on current biological research and a series of in-silico experiments which use differential evolution to fit the model to experimental data. From the obtained results, we formulate a series of biological hypotheses regarding currents involved in the neuron dynamics. These models reproduce experimental data with a high degree of accuracy while being biologically consistent with state-of-the-art research.

Optimization of Efficient Neuron Models With Realistic Firing Dynamics. The Case of the Cerebellar Granule Cell

Biologically relevant large-scale computational models currently represent one of the main methods in neuroscience for studying information processing primitives of brain areas. However, biologically realistic neuron models tend to be computationally heavy and thus prevent these models from being part of brain-area models including thousands or even millions of neurons. The cerebellar input layer represents a canonical example of large scale networks. In particular, the cerebellar granule cells, the most numerous cells in the whole mammalian brain, have been proposed as playing a pivotal role in the creation of somato-sensorial information representations. Enhanced burst frequency (spiking resonance) in the granule cells has been proposed as facilitating the input signal transmission at the theta-frequency band (4–12 Hz), but the functional role of this cell feature in the operation of the granular layer remains largely unclear. This study aims to develop a methodological pipeline for creating neuron models that maintain biological realism and computational efficiency whilst capturing essential aspects of single-neuron processing. Therefore, we selected a light computational neuron model template (the adaptive-exponential integrate-and-fire model), whose parameters were progressively refined using an automatic parameter tuning with evolutionary algorithms (EAs). The resulting point-neuron models are suitable for reproducing the main firing properties of a realistic granule cell from electrophysiological measurements, including the spiking resonance at the theta-frequency band, repetitive firing according to a specified intensity-frequency (I-F) curve and delayed firing under current-pulse stimulation. Interestingly, the proposed model also reproduced some other emergent properties (namely, silent at rest, rheobase and negligible adaptation under depolarizing currents) even though these properties were not set in the EA as a target in the fitness function (FF), proving that these features are compatible even in computationally simple models. The proposed methodology represents a valuable tool for adjusting AdEx models according to a FF defined in the spiking regime and based on biological data. These models are appropriate for future research of the functional implication of bursting resonance at the theta band in large-scale granular layer network models.

Itinerant complexity in networks of intrinsically bursting neurons

Active neurons can be broadly classified by their intrinsic oscillation patterns into two classes characterized by spiking or bursting. Here, we show that networks of identical bursting neurons with inhibitory pulsatory coupling exhibit itinerant dynamics. Using the relative phases of bursts between neurons, we numerically demonstrate that the network exhibits endogenous transitions between multiple modes of transient synchrony. This is true even for bursts consisting of two spikes. In contrast, our simulations reveal that networks of identical singlet-spiking neurons do not exhibit such complexity. These results suggest a role for bursting dynamics in realizing itinerant complexity in neural circuits.

Quantitative firing pattern phenotyping of hippocampal neuron types

Systematically organizing the anatomical, molecular, and physiological properties of cortical neurons is important for understanding their computational functions. Hippocampome.org defines 122 neuron types in the rodent hippocampal formation based on their somatic, axonal, and dendritic locations, putative excitatory/inhibitory outputs, molecular marker expression, and biophysical properties. We augmented the electrophysiological data of this knowledge base by collecting, quantifying, and analyzing the firing responses to depolarizing current injections for every hippocampal neuron type from published experiments. We designed and implemented objective protocols to classify firing patterns based on 5 transients (delay, adapting spiking, rapidly adapting spiking, transient stuttering, and transient slow-wave bursting) and 4 steady states (non-adapting spiking, persistent stuttering, persistent slow-wave bursting, and silence). This automated approach revealed 9 unique (plus one spurious) families of firing pattern phenotypes while distinguishing potential new neuronal subtypes. Novel statistical associations emerged between firing responses and other electrophysiological properties, morphological features, and molecular marker expression. The firing pattern parameters, experimental conditions, spike times, references to the original empirical evidences, and analysis scripts are released open-source through Hippocampome.org for all neuron types, greatly enhancing the existing search and browse capabilities. This information, collated online in human- and machine-accessible form, will help design and interpret both experiments and model simulations.

Simple models of quantitative firing phenotypes in hippocampal neurons: Comprehensive coverage of intrinsic diversity

Patterns of periodic voltage spikes elicited by a neuron help define its dynamical identity. Experimentally recorded spike trains from various neurons show qualitatively distinguishable features such as delayed spiking, spiking with or without frequency adaptation, and intrinsic bursting. Moreover, the input-dependent responses of a neuron not only show different quantitative features, such as higher spike frequency for a stronger input current injection, but can also exhibit qualitatively different responses, such as spiking and bursting under different input conditions, thus forming a complex phenotype of responses. In previous work, the comprehensive knowledge base of hippocampal neuron types Hippocampome.org systematically characterized various spike pattern phenotypes experimentally identified from 120 neuron types/subtypes. In this paper, we present a complete set of simple phenomenological models that quantitatively reproduce the diverse and complex phenotypes of hippocampal neurons. In addition to point-neuron models, we created compact multi-compartment models with up to four compartments, which will allow spatial segregation of synaptic integration in network simulations. Electrotonic compartmentalization observed in our compact multi-compartment models is qualitatively consistent with experimental observations. The models were created using an automated pipeline based on evolutionary algorithms. This work maps 120 neuron types/subtypes in the rodent hippocampus to a low-dimensional model space and adds another dimension to the knowledge accumulated in Hippocampome.org. Computationally efficient representations of intrinsic dynamics, along with other pieces of knowledge available in Hippocampome.org, provide a biologically realistic platform to explore the large-scale interactions of various neuron types at the mesoscopic level.

Complex Dynamics in Simplified Neuronal Models: Reproducing Golgi Cell Electroresponsiveness

Brain neurons exhibit complex electroresponsive properties – including intrinsic subthreshold oscillations and pacemaking, resonance and phase-reset – which are thought to play a critical role in controlling neural network dynamics. Although these properties emerge from detailed representations of molecular-level mechanisms in “realistic” models, they cannot usually be generated by simplified neuronal models (although these may show spike-frequency adaptation and bursting). We report here that this whole set of properties can be generated by the extended generalized leaky integrate-and-fire (E-GLIF) neuron model. E-GLIF derives from the GLIF model family and is therefore mono-compartmental, keeps the limited computational load typical of a linear low-dimensional system, admits analytical solutions and can be tuned through gradient-descent algorithms. Importantly, E-GLIF is designed to maintain a correspondence between model parameters and neuronal membrane mechanisms through a minimum set of equations. In order to test its potential, E-GLIF was used to model a specific neuron showing rich and complex electroresponsiveness, the cerebellar Golgi cell, and was validated against experimental electrophysiological data recorded from Golgi cells in acute cerebellar slices. During simulations, E-GLIF was activated by stimulus patterns, including current steps and synaptic inputs, identical to those used for the experiments. The results demonstrate that E-GLIF can reproduce the whole set of complex neuronal dynamics typical of these neurons – including intensity-frequency curves, spike-frequency adaptation, post-inhibitory rebound bursting, spontaneous subthreshold oscillations, resonance, and phase-reset – providing a new effective tool to investigate brain dynamics in large-scale simulations.