Quantitative firing pattern phenotyping of hippocampal neuron types

Formation and retrieval of cell assemblies in a biologically realistic spiking neural network model of area CA3 in the mouse hippocampus

The hippocampal formation is critical for episodic memory, with area Cornu Ammonis 3 (CA3) a necessary substrate for auto-associative pattern completion. Recent theoretical and experimental evidence suggests that the formation and retrieval of cell assemblies enable these functions. Yet, how cell assemblies are formed and retrieved in a full-scale spiking neural network (SNN) of CA3 that incorporates the observed diversity of neurons and connections within this circuit is not well understood. Here, we demonstrate that a data-driven SNN model quantitatively reflecting the neuron type-specific population sizes, intrinsic electrophysiology, connectivity statistics, synaptic signaling, and long-term plasticity of the mouse CA3 is capable of robust auto-association and pattern completion via cell assemblies. Our results show that a broad range of assembly sizes could successfully and systematically retrieve patterns from heavily incomplete or corrupted cues after a limited number of presentations. Furthermore, performance was robust with respect to partial overlap of assemblies through shared cells, substantially enhancing memory capacity. These novel findings provide computational evidence that the specific biological properties of the CA3 circuit produce an effective neural substrate for associative learning in the mammalian brain.

A Continuous Attractor Model with Realistic Neural and Synaptic Properties Quantitatively Reproduces Grid Cell Physiology

Computational simulations with data-driven physiological detail can foster a deeper understanding of the neural mechanisms involved in cognition. Here, we utilize the wealth of cellular properties from Hippocampome.org to study neural mechanisms of spatial coding with a spiking continuous attractor network model of medial entorhinal cortex circuit activity. The primary goal is to investigate if adding such realistic constraints could produce firing patterns similar to those measured in real neurons. Biological characteristics included in the work are excitability, connectivity, and synaptic signaling of neuron types defined primarily by their axonal and dendritic morphologies. We investigate the spiking dynamics in specific neuron types and the synaptic activities between groups of neurons. Modeling the rodent hippocampal formation keeps the simulations to a computationally reasonable scale while also anchoring the parameters and results to experimental measurements. Our model generates grid cell activity that well matches the spacing, size, and firing rates of grid fields recorded in live behaving animals from both published datasets and new experiments performed for this study. Our simulations also recreate different scales of those properties, e.g., small and large, as found along the dorsoventral axis of the medial entorhinal cortex. Computational exploration of neuronal and synaptic model parameters reveals that a broad range of neural properties produce grid fields in the simulation. These results demonstrate that the continuous attractor network model of grid cells is compatible with a spiking neural network implementation sourcing data-driven biophysical and anatomical parameters from Hippocampome.org. The software (version 1.0) is released as open source to enable broad community reuse and encourage novel applications.

Metabolic regulation of single synaptic vesicle exo- and endocytosis in hippocampal synapses

Glucose has long been considered a primary energy source for synaptic function. However, it remains unclear to what extent alternative fuels, such as lactate/pyruvate, contribute to powering synaptic transmission. By detecting individual release events in hippocampal synapses, we find that mitochondrial ATP production regulates basal vesicle release probability and release location within the active zone (AZ), evoked by single action potentials. Mitochondrial inhibition shifts vesicle release closer to the AZ center and alters the efficiency of vesicle retrieval by increasing the occurrence of ultrafast endocytosis. Furthermore, we uncover that terminals can use oxidative fuels to maintain the vesicle cycle during trains of activity. Mitochondria are sparsely distributed along hippocampal axons, and we find that terminals containing mitochondria display enhanced vesicle release and reuptake during high-frequency trains. Our findings suggest that mitochondria not only regulate several fundamental features of synaptic transmission but may also contribute to modulation of short-term synaptic plasticity.

A Continuous Attractor Model with Realistic Neural and Synaptic Properties Quantitatively Reproduces Grid Cell Physiology

Computational simulations with data-driven physiological detail can foster a deeper understanding of the neural mechanisms involved in cognition. Here, we utilize the wealth of cellular properties from Hippocampome.org to study neural mechanisms of spatial coding with a spiking continuous attractor network model of medial entorhinal cortex circuit activity. The primary goal was to investigate if adding such realistic constraints could produce firing patterns similar to those measured in real neurons. Biological characteristics included in the work are excitability, connectivity, and synaptic signaling of neuron types defined primarily by their axonal and dendritic morphologies. We investigate the spiking dynamics in specific neuron types and the synaptic activities between groups of neurons. Modeling the rodent hippocampal formation keeps the simulations to a computationally reasonable scale while also anchoring the parameters and results to experimental measurements. Our model generates grid cell activity that well matches the spacing, size, and firing rates of grid fields recorded in live behaving animals from both published datasets and new experiments performed for this study. Our simulations also recreate different scales of those properties, e.g., small and large, as found along the dorsoventral axis of the medial entorhinal cortex. Computational exploration of neuronal and synaptic model parameters reveals that a broad range of neural properties produce grid fields in the simulation. These results demonstrate that the continuous attractor network model of grid cells is compatible with a spiking neural network implementation sourcing data-driven biophysical and anatomical parameters from Hippocampome.org. The software is released as open source to enable broad community reuse and encourage novel applications.

Formation and Retrieval of Cell Assemblies in a Biologically Realistic Spiking Neural Network Model of Area CA3 in the Mouse Hippocampus

The hippocampal formation is critical for episodic memory, with area Cornu Ammonis 3 (CA3) a necessary substrate for auto-associative pattern completion. Recent theoretical and experimental evidence suggests that the formation and retrieval of cell assemblies enable these functions. Yet, how cell assemblies are formed and retrieved in a full-scale spiking neural network (SNN) of CA3 that incorporates the observed diversity of neurons and connections within this circuit is not well understood. Here, we demonstrate that a data-driven SNN model quantitatively reflecting the neuron type-specific population sizes, intrinsic electrophysiology, connectivity statistics, synaptic signaling, and long-term plasticity of the mouse CA3 is capable of robust auto-association and pattern completion via cell assemblies. Our results show that a broad range of assembly sizes could successfully and systematically retrieve patterns from heavily incomplete or corrupted cues after a limited number of presentations. Furthermore, performance was robust with respect to partial overlap of assemblies through shared cells, substantially enhancing memory capacity. These novel findings provide computational evidence that the specific biological properties of the CA3 circuit produce an effective neural substrate for associative learning in the mammalian brain.

Biophysical modulation and robustness of itinerant complexity in neuronal networks

Transient synchronization of bursting activity in neuronal networks, which occurs in patterns of metastable itinerant phase relationships between neurons, is a notable feature of network dynamics observed in vivo. However, the mechanisms that contribute to this dynamical complexity in neuronal circuits are not well understood. Local circuits in cortical regions consist of populations of neurons with diverse intrinsic oscillatory features. In this study, we numerically show that the phenomenon of transient synchronization, also referred to as metastability, can emerge in an inhibitory neuronal population when the neurons' intrinsic fast-spiking dynamics are appropriately modulated by slower inputs from an excitatory neuronal population. Using a compact model of a mesoscopic-scale network consisting of excitatory pyramidal and inhibitory fast-spiking neurons, our work demonstrates a relationship between the frequency of pyramidal population oscillations and the features of emergent metastability in the inhibitory population. In addition, we introduce a method to characterize collective transitions in metastable networks. Finally, we discuss potential applications of this study in mechanistically understanding cortical network dynamics.

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.

Metabolic Regulation of Single Synaptic Vesicle Exo- and Endocytosis in Hippocampal Synapses

Glucose has long been considered a primary source of energy for synaptic function. However, it remains unclear under what conditions alternative fuels, such as lactate/pyruvate, contribute to powering synaptic transmission. By detecting individual release events in cultured hippocampal synapses, we found that mitochondrial ATP production from oxidation of lactate/pyruvate regulates basal vesicle release probability and release location within the active zone (AZ) evoked by single action potentials (APs). Mitochondrial inhibition shifted vesicle release closer to the AZ center, suggesting that the energetic barrier for vesicle release is lower in the AZ center that the periphery. Mitochondrial inhibition also altered the efficiency of single AP evoked vesicle retrieval by increasing occurrence of ultrafast endocytosis, while inhibition of glycolysis had no effect. Mitochondria are sparsely distributed along hippocampal axons and we found that nerve terminals containing mitochondria displayed enhanced vesicle release and reuptake during high-frequency trains, irrespective of whether neurons were supplied with glucose or lactate. Thus, synaptic terminals can entirely bypass glycolysis to robustly maintain the vesicle cycle using oxidative fuels in the absence of glucose. These observations further suggest that mitochondrial metabolic function not only regulates several fundamental features of synaptic transmission but may also contribute to modulation of short-term synaptic plasticity.

Phase relations of interneuronal activity relative to theta rhythm

The theta rhythm plays a crucial role in synchronizing neural activity during attention and memory processes. However, the mechanisms behind the formation of neural activity during theta rhythm generation remain unknown. To address this, we propose a mathematical model that explains the distribution of interneurons in the CA1 field during the theta rhythm phase. Our model consists of a network of seven types of interneurons in the CA1 field that receive inputs from the CA3 field, entorhinal cortex, and local pyramidal neurons in the CA1 field. By adjusting the parameters of the connections in the model. We demonstrate that it is possible to replicate the experimentally observed phase relations between interneurons and the theta rhythm. Our model predicts that populations of interneurons receive unimodal excitation and inhibition with coinciding peaks, and that excitation dominates to determine the firing dynamics of interneurons.

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.

Multivesicular release favors short term synaptic depression in hippocampal autapses

Presynaptic terminals of the central nervous system can support univesicular and multivesicular synchronous release of neurotransmitters, however, the functional implications of the prevalence of one mechanism over the other are yet unresolved. Here, we took advantage of the expression of SF-iGluSnFR.S72A in the astrocytic feeder layer of autaptic hippocampal neuronal cultures to associate the liberation of glutamate to excitatory postsynaptic currents. The presence of the glutamate sensor in glial cells avoided any interference with the function of endogenous postsynaptic receptors. It was possible to optically detect changes in neurotransmitter release probability, which was heterogeneous among synaptic boutons studied. For each neuron investigated, the liberation of neurotransmitters occurred through a predominant mechanism. The prevalence of multivesicular over univesicular release increased synaptic strength and enhanced short-term synaptic depression. These results show that the preference of hippocampal boutons to synchronously release one or more vesicles determines the strength and low pass filtering properties of the synapses established.

ProtTrans: Toward Understanding the Language of Life Through Self-Supervised Learning

Computational biology and bioinformatics provide vast data gold-mines from protein sequences, ideal for Language Models (LMs) taken from Natural Language Processing (NLP). These LMs reach for new prediction frontiers at low inference costs. Here, we trained two auto-regressive models (Transformer-XL, XLNet) and four auto-encoder models (BERT, Albert, Electra, T5) on data from UniRef and BFD containing up to 393 billion amino acids. The protein LMs (pLMs) were trained on the Summit supercomputer using 5616 GPUs and TPU Pod up-to 1024 cores. Dimensionality reduction revealed that the raw pLM-embeddings from unlabeled data captured some biophysical features of protein sequences. We validated the advantage of using the embeddings as exclusive input for several subsequent tasks: (1) a per-residue (per-token) prediction of protein secondary structure (3-state accuracy Q3=81%-87%); (2) per-protein (pooling) predictions of protein sub-cellular location (ten-state accuracy: Q10=81%) and membrane versus water-soluble (2-state accuracy Q2=91%). For secondary structure, the most informative embeddings (ProtT5) for the first time outperformed the state-of-the-art without multiple sequence alignments (MSAs) or evolutionary information thereby bypassing expensive database searches. Taken together, the results implied that pLMs learned some of the grammar of the language of life. All our models are available through https://github.com/agemagician/ProtTrans.

Normalized unitary synaptic signaling of the hippocampus and entorhinal cortex predicted by deep learning of experimental recordings

Biologically realistic computer simulations of neuronal circuits require systematic data-driven modeling of neuron type-specific synaptic activity. However, limited experimental yield, heterogeneous recordings conditions, and ambiguous neuronal identification have so far prevented the consistent characterization of synaptic signals for all connections of any neural system. We introduce a strategy to overcome these challenges and report a comprehensive synaptic quantification among all known neuron types of the hippocampal-entorhinal network. First, we reconstructed >2600 synaptic traces from ∼1200 publications into a unified computational representation of synaptic dynamics. We then trained a deep learning architecture with the resulting parameters, each annotated with detailed metadata such as recording method, solutions, and temperature. The model learned to predict the synaptic properties of all 3,120 circuit connections in arbitrary conditions with accuracy approaching the intrinsic experimental variability. Analysis of data normalized and completed with the deep learning model revealed that synaptic signals are controlled by few latent variables associated with specific molecular markers and interrelating conductance, decay time constant, and short-term plasticity. We freely release the tools and full dataset of unitary synaptic values in 32 covariate settings. Normalized synaptic data can be used in brain simulations, and to predict and test experimental hypothesis.

A deep learning model trained on roughly 2,600 synaptic traces from hippocampal electrophysiology datasets demonstrates how specific covariates influence synaptic signals.

Quantification of neuron types in the rodent hippocampal formation by data mining and numerical optimization

Quantifying the population sizes of distinct neuron types in different anatomical regions is an essential step towards establishing a brain cell census. Although estimates exist for the total neuronal populations in different species, the number and definition of each specific neuron type are still intensively investigated. Hippocampome.org is an open-source knowledge base with morphological, physiological and molecular information for 122 neuron types in the rodent hippocampal formation. While such framework identifies all known neuron types in this system, their relative abundances remain largely unknown. This work quantitatively estimates the counts of all Hippocampome.org neuron types by literature mining and numerical optimization. We report the number of neurons in each type identified by main neurotransmitter (glutamate or GABA) and axonal-dendritic patterns throughout 26 subregions and layers of the dentate gyrus, Ammon's horn, subiculum and entorhinal cortex. We produce by sensitivity analysis reliable numerical ranges for each type and summarize the amounts across broad neuronal families defined by biomarkers expression and firing dynamics. Study of density distributions indicates that the number of dendritic-targeting interneurons, but not of other neuronal classes, is independent of anatomical volumes. All extracted values, experimental evidence and related software code are released on Hippocampome.org.

Determining clinically relevant features in cytometry data using persistent homology

Cytometry experiments yield high-dimensional point cloud data that is difficult to interpret manually. Boolean gating techniques coupled with comparisons of relative abundances of cellular subsets is the current standard for cytometry data analysis. However, this approach is unable to capture more subtle topological features hidden in data, especially if those features are further masked by data transforms or significant batch effects or donor-to-donor variations in clinical data. We present that persistent homology, a mathematical structure that summarizes the topological features, can distinguish different sources of data, such as from groups of healthy donors or patients, effectively. Analysis of publicly available cytometry data describing non-naïve CD8+ T cells in COVID-19 patients and healthy controls shows that systematic structural differences exist between single cell protein expressions in COVID-19 patients and healthy controls. We identify proteins of interest by a decision-tree based classifier, sample points randomly and compute persistence diagrams from these sampled points. The resulting persistence diagrams identify regions in cytometry datasets of varying density and identify protruded structures such as ‘elbows’. We compute Wasserstein distances between these persistence diagrams for random pairs of healthy controls and COVID-19 patients and find that systematic structural differences exist between COVID-19 patients and healthy controls in the expression data for T-bet, Eomes, and Ki-67. Further analysis shows that expression of T-bet and Eomes are significantly downregulated in COVID-19 patient non-naïve CD8+ T cells compared to healthy controls. This counter-intuitive finding may indicate that canonical effector CD8+ T cells are less prevalent in COVID-19 patients than healthy controls. This method is applicable to any cytometry dataset for discovering novel insights through topological data analysis which may be difficult to ascertain otherwise with a standard gating strategy or existing bioinformatic tools.

Identifying differences between cytometry data seen as a point cloud can be complicated by random variations in data collection and data sources. We apply persistent homology used in topological data analysis to describe the shape and structure of the data representing immune cells in healthy donors and COVID-19 patients. By looking at how the shape and structure differ between healthy donors and COVID-19 patients, we are able to definitively conclude how these groups differ despite random variations in the data. Furthermore, these results are novel in their ability to capture shape and structure of cytometry data, something not described by other analyses.

Reconstruction of the Hippocampus

The hippocampus is a widely studied brain region thought to play an important role in higher cognitive functions such as learning, memory, and navigation. The amount of data on this region increases every day and delineates a complex and fragmented picture, but an integrated understanding of hippocampal function remains elusive. Computational methods can help to move the research forward, and reconstructing a full-scale model of the hippocampus is a challenging yet feasible task that the research community should undertake.

In this chapter, we present strategies for reconstructing a large-scale model of the hippocampus. Based on a previously published approach to reconstruct and simulate brain tissue, which is also explained in Chap. 10, we discuss the characteristics of the hippocampus in the light of its special anatomical and physiological features, data availability, and existing large-scale hippocampus models. A large-scale model of the hippocampus is a compound model of several building blocks: ion channels, morphologies, single cell models, connections, synapses. We discuss each of those building blocks separately and discuss how to merge them back and simulate the resulting network model.

Integrative Models of Brain Structure and Dynamics: Concepts, Challenges, and Methods

The anatomical architecture of the brain constrains the dynamics of interactions between various regions. On a microscopic scale, neural plasticity regulates the connections between individual neurons. This microstructural adaptation facilitates coordinated dynamics of populations of neurons (mesoscopic scale) and brain regions (macroscopic scale). However, the mechanisms acting on multiple timescales that govern the reciprocal relationship between neural network structure and its intrinsic dynamics are not well understood. Studies empirically investigating such relationships on the whole-brain level rely on macroscopic measurements of structural and functional connectivity estimated from various neuroimaging modalities such as Diffusion-weighted Magnetic Resonance Imaging (dMRI), Electroencephalography (EEG), Magnetoencephalography (MEG), and functional Magnetic Resonance Imaging (fMRI). dMRI measures the anisotropy of water diffusion along axonal fibers, from which structural connections are estimated. EEG and MEG signals measure electrical activity and magnetic fields induced by the electrical activity, respectively, from various brain regions with a high temporal resolution (but limited spatial coverage), whereas fMRI measures regional activations indirectly via blood oxygen level-dependent (BOLD) signals with a high spatial resolution (but limited temporal resolution). There are several studies in the neuroimaging literature reporting statistical associations between macroscopic structural and functional connectivity. On the other hand, models of large-scale oscillatory dynamics conditioned on network structure (such as the one estimated from dMRI connectivity) provide a platform to probe into the structure-dynamics relationship at the mesoscopic level. Such investigations promise to uncover the theoretical underpinnings of the interplay between network structure and dynamics and could be complementary to the macroscopic level inquiries. In this article, we review theoretical and empirical studies that attempt to elucidate the coupling between brain structure and dynamics. Special attention is given to various clinically relevant dimensions of brain connectivity such as the topological features and neural synchronization, and their applicability for a given modality, spatial or temporal scale of analysis is discussed. Our review provides a summary of the progress made along this line of research and identifies challenges and promising future directions for multi-modal neuroimaging analyses.

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

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

Spiking Neural Networks and Hippocampal Function: A Web-Accessible Survey of Simulations, Modeling Methods, and Underlying Theories

Computational modeling has contributed to hippocampal research in a wide variety of ways and through a large diversity of approaches, reflecting the many advanced cognitive roles of this brain region. The intensively studied neuron type circuitry of the hippocampus is a particularly conducive substrate for spiking neural models. Here we present an online knowledge base of spiking neural network simulations of hippocampal functions. First, we overview theories involving the hippocampal formation in subjects such as spatial representation, learning, and memory. Then we describe an original literature mining process to organize published reports in various key aspects, including: (i) subject area (e.g., navigation, pattern completion, epilepsy); (ii) level of modeling detail (Hodgkin-Huxley, integrate-and-fire, etc.); and (iii) theoretical framework (attractor dynamics, oscillatory interference, self-organizing maps, and others). Moreover, every peer-reviewed publication is also annotated to indicate the specific neuron types represented in the network simulation, establishing a direct link with the Hippocampome.org portal. The web interface of the knowledge base enables dynamic content browsing and advanced searches, and consistently presents evidence supporting every annotation. Moreover, users are given access to several types of statistical reports about the collection, a selection of which is summarized in this paper. This open access resource thus provides an interactive platform to survey spiking neural network models of hippocampal functions, compare available computational methods, and foster ideas for suitable new directions of research.

Development, Diversity, and Death of MGE-Derived Cortical Interneurons

In the mammalian brain, cortical interneurons (INs) are a highly diverse group of cells. A key neurophysiological question concerns how each class of INs contributes to cortical circuit function and whether specific roles can be attributed to a selective cell type. To address this question, researchers are integrating knowledge derived from transcriptomic, histological, electrophysiological, developmental, and functional experiments to extensively characterise the different classes of INs. Our hope is that such knowledge permits the selective targeting of cell types for therapeutic endeavours. This review will focus on two of the main types of INs, namely the parvalbumin (PV⁺) or somatostatin (SOM⁺)-containing cells, and summarise the research to date on these classes.

Membrane Resonance in Pyramidal and GABAergic Neurons of the Mouse Perirhinal Cortex

The perirhinal cortex (PRC) is a polymodal associative region of the temporal lobe that works as a gateway between cortical areas and hippocampus. In recent years, an increasing interest arose in the role played by the PRC in learning and memory processes, such as object recognition memory, in contrast with certain forms of hippocampus-dependent spatial and episodic memory. The integrative properties of the PRC should provide all necessary resources to select and enhance the information to be propagated to and from the hippocampus. Among these properties, we explore in this paper the ability of the PRC neurons to amplify the output voltage to current input at selected frequencies, known as membrane resonance. Within cerebral circuits the resonance of a neuron operates as a filter toward inputs signals at certain frequencies to coordinate network activity in the brain by affecting the rate of neuronal firing and the precision of spike timing. Furthermore, the ability of the PRC neurons to resonate could have a fundamental role in generating subthreshold oscillations and in the selection of cortical inputs directed to the hippocampus. Here, performing whole-cell patch-clamp recordings from perirhinal pyramidal neurons and GABAergic interneurons of GAD67-GFP⁺ mice, we found, for the first time, that the majority of PRC neurons are resonant at their resting potential, with a resonance frequency of 0.5–1.5 Hz at 23°C and of 1.5–2.8 Hz at 36°C. In the presence of ZD7288 (blocker of HCN channels) resonance was abolished in both pyramidal neurons and interneurons, suggesting that I_h current is critically involved in resonance generation. Otherwise, application of TTx (voltage-dependent Na⁺ channel blocker) attenuates the resonance in pyramidal neurons but not in interneurons, suggesting that only in pyramidal neurons the persistent sodium current has an amplifying effect. These experimental results have also been confirmed by a computational model. From a functional point of view, the resonance in the PRC would affect the reverberating activity between neocortex and hippocampus, especially during slow wave sleep, and could be involved in the redistribution and strengthening of memory representation in cortical regions.

GABAA receptor-mediated inhibition of Dahlgren cells electrical activity in the olive flounder, Paralichthys olivaceus

γ-Aminobutyric acid (GABA) is a major inhibitory neurotransmitter in the central nervous system. We investigated its potential role as a neurotransmitter in the neuroendocrine Dahlgren cell population of the caudal neurosecretory system (CNSS) of the flounder Paralichthys olivaceus. The application of GABA in vitro resulted in a decrease in electrical activity of Dahlgren cells, followed by an increase of the number of silent cells, together with a decreased firing frequency of all three activity patterns (tonic, phasic, bursting). GABA_A receptor agonist etomidate decreased Dahlgren cell firing activity, in a similar way to GABA. The response to GABA was blocked by the GABA_A receptor antagonist bicuculline. GABA_A receptor gamma2 subunit (Gabrg2) and chloride channel (Clcn2) mRNA expression were significantly upregulated in the CNSS after GABA superfusion. These data suggest that GABA may modulate CNSS activity in vivo mediated by GABA_A receptors.

An update to Hippocampome.org by integrating single-cell phenotypes with circuit function in vivo

Understanding brain operation demands linking basic behavioral traits to cell-type specific dynamics of different brain-wide subcircuits. This requires a system to classify the basic operational modes of neurons and circuits. Single-cell phenotyping of firing behavior during ongoing oscillations in vivo has provided a large body of evidence on entorhinal-hippocampal function, but data are dispersed and diverse. Here, we mined literature to search for information regarding the phase-timing dynamics of over 100 hippocampal/entorhinal neuron types defined in Hippocampome.org. We identified missing and unresolved pieces of knowledge (e.g., the preferred theta phase for a specific neuron type) and complemented the dataset with our own new data. By confronting the effect of brain state and recording methods, we highlight the equivalences and differences across conditions and offer a number of novel observations. We show how a heuristic approach based on oscillatory features of morphologically identified neurons can aid in classifying extracellular recordings of single cells and discuss future opportunities and challenges towards integrating single-cell phenotypes with circuit function.

Determining clinically relevant features in cytometry data using persistent homology

Cytometry experiments yield high-dimensional point cloud data that is difficult to interpret manually. Boolean gating techniques coupled with comparisons of relative abundances of cellular subsets is the current standard for cytometry data analysis. However, this approach is unable to capture more subtle topological features hidden in data, especially if those features are further masked by data transforms or significant batch effects or donor-to-donor variations in clinical data. We present that persistent homology, a mathematical structure that summarizes the topological features, can distinguish different sources of data, such as from groups of healthy donors or patients, effectively. Analysis of publicly available cytometry data describing non-naïve CD8+ T cells in COVID-19 patients and healthy controls shows that systematic structural differences exist between single cell protein expressions in COVID-19 patients and healthy controls. Our method identifies proteins of interest by a decision-tree based classifier and passes them to a kernel-density estimator (KDE) for sampling points from the density distribution. We then compute persistence diagrams from these sampled points. The resulting persistence diagrams identify regions in cytometry datasets of varying density and identify protruded structures such as 'elbows'. We compute Wasserstein distances between these persistence diagrams for random pairs of healthy controls and COVID-19 patients and find that systematic structural differences exist between COVID-19 patients and healthy controls in the expression data for T-bet, Eomes, and Ki-67. Further analysis shows that expression of T-bet and Eomes are significantly downregulated in COVID-19 patient non-naïve CD8+ T cells compared to healthy controls. This counter-intuitive finding may indicate that canonical effector CD8+ T cells are less prevalent in COVID-19 patients than healthy controls. This method is applicable to any cytometry dataset for discovering novel insights through topological data analysis which may be difficult to ascertain otherwise with a standard gating strategy or in the presence of large batch effects.

Author summary

Identifying differences between cytometry data seen as a point cloud can be complicated by random variations in data collection and data sources. We apply persistent homology used in topological data analysis to describe the shape and structure of the data representing immune cells in healthy donors and COVID-19 patients. By looking at how the shape and structure differ between healthy donors and COVID-19 patients, we are able to definitively conclude how these groups differ despite random variations in the data. Furthermore, these results are novel in their ability to capture shape and structure of cytometry data, something not described by other analyses.

Comprehensive Estimates of Potential Synaptic Connections in Local Circuits of the Rodent Hippocampal Formation by Axonal-Dendritic Overlap

A quantitative description of the hippocampal formation synaptic architecture is essential for understanding the neural mechanisms of episodic memory. Yet the existing knowledge of connectivity statistics between different neuron types in the rodent hippocampus only captures a mere 5% of this circuitry. We present a systematic pipeline to produce first-approximation estimates for most of the missing information. Leveraging the www.Hippocampome.org knowledge base, we derive local connection parameters between distinct pairs of morphologically identified neuron types based on their axonal-dendritic overlap within every layer and subregion of the hippocampal formation. Specifically, we adapt modern image analysis technology to determine the parcel-specific neurite lengths of every neuron type from representative morphologic reconstructions obtained from either sex. We then compute the average number of synapses per neuron pair using relevant anatomic volumes from the mouse brain atlas and ultrastructurally established interaction distances. Hence, we estimate connection probabilities and number of contacts for >1900 neuron type pairs, increasing the available quantitative assessments more than 11-fold. Connectivity statistics thus remain unknown for only a minority of potential synapses in the hippocampal formation, including those involving long-range (23%) or perisomatic (6%) connections and neuron types without morphologic tracings (7%). The described approach also yields approximate measurements of synaptic distances from the soma along the dendritic and axonal paths, which may affect signal attenuation and delay. Overall, this dataset fills a substantial gap in quantitatively describing hippocampal circuits and provides useful model specifications for biologically realistic neural network simulations, until further direct experimental data become available.SIGNIFICANCE STATEMENT The hippocampal formation is a crucial functional substrate for episodic memory and spatial representation. Characterizing the complex neuron type circuit of this brain region is thus important to understand the cellular mechanisms of learning and navigation. Here we present the first numerical estimates of connection probabilities, numbers of contacts per connected pair, and synaptic distances from the soma along the axonal and dendritic paths, for more than 1900 distinct neuron type pairs throughout the dentate gyrus, CA3, CA2, CA1, subiculum, and entorhinal cortex. This comprehensive dataset, publicly released online at www.Hippocampome.org, constitutes an unprecedented quantification of the majority of the local synaptic circuit for a prominent mammalian neural system and provides an essential foundation for data-driven, anatomically realistic neural 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.

Diabetic encephalopathy causes the imbalance of neural activities between hippocampal glutamatergic neurons and GABAergic neurons in mice

Diabetic encephalopathy is a severe diabetes-related complication in the central nervous system (CNS) that is characterized by the impairment of neurochemical and structural changes leading to cognitive dysfunction. Its cellular and molecular mechanisms are still unclear and clinical approaches are still lacking of promising therapies. In this study, we have investigated the changes of different hippocampal neurons during diabetic encephalopathy in mouse models of diabetes by simultaneously analyzing the activities and synaptic transmission of glutamatergic neurons and GABAergic neurons in brain slices. Compared with the data from a group of control, diabetic encephalopathy permanently impairs the excitability of GABAergic neurons and synaptic transmission mediated by γ-aminobutyric acid (GABA). However, glutamatergic neurons appear to be more excited. Our findings highlight the critical role of the dysfunction of GABAergic neurons and glutamatergic neurons during diabetic encephalopathy in hippocampus to neural impairment as well as a strategy to prevent the function of progress of diabetic encephalopathy by protecting central neurons.