A Novel Method for Clustering Cellular Data to Improve Classification.

Many fields, such as neuroscience, are experiencing the vast proliferation of cellular data, underscoring the need for organizing and interpreting large datasets. A popular approach partitions data into manageable subsets via hierarchical clustering, but objective methods to determine the appropriate classification granularity are missing. We recently introduced a technique to systematically identify when to stop subdividing clusters based on the fundamental principle that cells must differ more between than within clusters. Here we present the corresponding protocol to classify cellular datasets by combining data-driven unsupervised hierarchical clustering with statistical testing. These general-purpose functions are applicable to any cellular dataset that can be organized as two-dimensional matrices of numerical values, including molecular, physiological, and anatomical datasets. We demonstrate the protocol using cellular data from the Janelia MouseLight project to characterize morphological aspects of neurons.

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.

Unsupervised classification of brain-wide axons reveals the presubiculum neuronal projection blueprint.

We present a quantitative strategy to identify all projection neuron types from a given region with statistically different patterns of anatomical targeting. We first validate the technique with mouse primary motor cortex layer 6 data, yielding two clusters consistent with cortico-thalamic and intra-telencephalic neurons. We next analyze the presubiculum, a less-explored region, identifying five classes of projecting neurons with unique patterns of divergence, convergence, and specificity. We report several findings: individual classes target multiple subregions along defined functions; all hypothalamic regions are exclusively targeted by the same class also invading midbrain and agranular retrosplenial cortex; Cornu Ammonis receives input from a single class of presubicular axons also projecting to granular retrosplenial cortex; path distances from the presubiculum to the same targets differ significantly between classes, as do the path distances to distinct targets within most classes; the identified classes have highly non-uniform abundances; and presubicular somata are topographically segregated among classes. This study thus demonstrates that statistically distinct projections shed light on the functional organization of their circuit.

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.

Unsupervised classification of brain-wide axons reveals neuronal projection blueprint.

Long-range axonal projections are quintessential determinants of network connectivity, linking cellular organization and circuit architecture. Here we introduce a quantitative strategy to identify, from a given source region, all "projection neuron types" with statistically different patterns of anatomical targeting. We first validate the proposed technique with well-characterized data from layer 6 of the mouse primary motor cortex. The results yield two clusters, consistent with previously discovered cortico-thalamic and intra-telencephalic neuron classes. We next analyze neurons from the presubiculum, a less-explored region. Extending sparse knowledge from earlier retrograde tracing studies, we identify five classes of presubicular projecting neurons, revealing unique patterns of divergence, convergence, and specificity. We thus report several findings: (1) individual classes target multiple subregions along defined functions, such as spatial representation vs. sensory integration and visual vs. auditory input; (2) all hypothalamic regions are exclusively targeted by the same class also invading midbrain, a sharp subset of thalamic nuclei, and agranular retrosplenial cortex; (3) Cornu Ammonis, in contrast, receives input from the same presubicular axons projecting to granular retrosplenial cortex, also the purview of a single class; (4) path distances from the presubiculum to the same targets differ significantly between classes, as do the path distances to distinct targets within most classes, suggesting fine temporal coordination in activating distant areas; (5) the identified classes have highly non-uniform abundances, with substantially more neurons projecting to midbrain and hypothalamus than to medial and lateral entorhinal cortex; (6) lastly, presubicular soma locations are segregated among classes, indicating topographic organization of projections. This study thus demonstrates that classifying neurons based on statistically distinct axonal projection patterns sheds light on the functional organizational of their circuit.

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.

Cellular anatomy of the mouse primary motor cortex.

An essential step toward understanding brain function is to establish a structural framework with cellular resolution on which multi-scale datasets spanning molecules, cells, circuits and systems can be integrated and interpreted1. Here, as part of the collaborative Brain Initiative Cell Census Network (BICCN), we derive a comprehensive cell type-based anatomical description of one exemplar brain structure, the mouse primary motor cortex, upper limb area (MOp-ul). Using genetic and viral labelling, barcoded anatomy resolved by sequencing, single-neuron reconstruction, whole-brain imaging and cloud-based neuroinformatics tools, we delineated the MOp-ul in 3D and refined its sublaminar organization. We defined around two dozen projection neuron types in the MOp-ul and derived an input-output wiring diagram, which will facilitate future analyses of motor control circuitry across molecular, cellular and system levels. This work provides a roadmap towards a comprehensive cellular-resolution description of mammalian brain architecture.

A Method for Estimating the Potential Synaptic Connections Between Axons and Dendrites From 2D Neuronal Images.

Computational neuroscience aims to model, reproduce, and predict network dynamics for different neuronal ensembles by distilling knowledge derived from electrophysiological and morphological evidence. However, analyses and simulations often remain critically limited by the sparsity of direct experimental constraints on essential parameters, such as electron microscopy and electrophysiology pair/multiple recording evidence of connectivity statistics. Notably, available data are particularly scarce regarding quantitative information on synaptic connections among identified neuronal types. Here, we present a user-friendly data-driven pipeline to estimate connection probabilities, number of contacts per connected pair, and distances from the pre- and postsynaptic somas along the axonal and dendritic paths from commonly available two-dimensional tracings and other broadly accessible measurements. The described procedure does not require any computational background and is accessible to all neuroscientists. This protocol therefore fills the important gap from neuronal morphology to circuit organization and can be applied to many different neural systems, brain regions, animal species, and data sources. Graphic abstract: The processing protocol from 2D reconstructions to quantitated synaptic connections.

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.

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.

Molecular expression profiles of morphologically defined hippocampal neuron types: Empirical evidence and relational inferences.

Gene and protein expressions are key determinants of cellular function. Neurons are the building blocks of brain circuits, yet the relationship between their molecular identity and the spatial distribution of their dendritic inputs and axonal outputs remains incompletely understood. The open-source knowledge base Hippocampome.org amasses such transcriptomic data from the scientific literature for morphologically defined neuron types in the rodent hippocampal formation: dentate gyrus, CA3, CA2, CA1, subiculum, and entorhinal cortex. Positive, negative, or mixed expression reports were initially obtained from published articles directly connecting molecular evidence to neurons with known axonal and dendritic patterns across hippocampal layers. Here, we supplement this information by collating, formalizing, and leveraging relational expression inferences that link a gene or protein expression or lack thereof to that of another molecule or to an anatomical location. With these additional interpretations, we freely release online a comprehensive human- and machine-readable molecular profile for more than 100 neuron types in Hippocampome.org. Analysis of these data ascertains the ability to distinguish unequivocally most neuron types in each of the major subdivisions of the hippocampus based on currently known biochemical markers. Moreover, grouping neuron types by expression similarity reveals eight superfamilies characterized by a few defining molecules.

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.

Graph Theoretic and Motif Analyses of the Hippocampal Neuron Type Potential Connectome.

We computed the potential connectivity map of all known neuron types in the rodent hippocampal formation by supplementing scantly available synaptic data with spatial distributions of axons and dendrites from the open-access knowledge base Hippocampome.org. The network that results from this endeavor, the broadest and most complete for a mammalian cortical region at the neuron-type level to date, contains more than 3200 connections among 122 neuron types across six subregions. Analyses of these data using graph theory metrics unveil the fundamental architectural principles of the hippocampal circuit. Globally, we identify a highly specialized topology minimizing communication cost; a modular structure underscoring the prominence of the trisynaptic loop; a core set of neuron types serving as information-processing hubs as well as a distinct group of particular antihub neurons; a nested, two-tier rich club managing much of the network traffic; and an innate resilience to random perturbations. At the local level, we uncover the basic building blocks, or connectivity patterns, that combine to produce complex global functionality, and we benchmark their utilization in the circuit relative to random networks. Taken together, these results provide a comprehensive connectivity profile of the hippocampus, yielding novel insights on its functional operations at the computationally crucial level of neuron types.

In search of a periodic table of the neurons: Axonal-dendritic circuitry as the organizing principle: Patterns of axons and dendrites within distinct anatomical parcels provide the blueprint for circuit-based neuronal classification.

No one knows yet how to organize, in a simple yet predictive form, the knowledge concerning the anatomical, biophysical, and molecular properties of neurons that are accumulating in thousands of publications every year. The situation is not dissimilar to the state of Chemistry prior to Mendeleev's tabulation of the elements. We propose that the patterns of presence or absence of axons and dendrites within known anatomical parcels may serve as the key principle to define neuron types. Just as the positions of the elements in the periodic table indicate their potential to combine into molecules, axonal and dendritic distributions provide the blueprint for network connectivity. Furthermore, among the features commonly employed to describe neurons, morphology is considerably robust to experimental conditions. At the same time, this core classification scheme is suitable for aggregating biochemical, physiological, and synaptic information.

Hippocampome.org: a knowledge base of neuron types in the rodent hippocampus.

Hippocampome.org is a comprehensive knowledge base of neuron types in the rodent hippocampal formation (dentate gyrus, CA3, CA2, CA1, subiculum, and entorhinal cortex). Although the hippocampal literature is remarkably information-rich, neuron properties are often reported with incompletely defined and notoriously inconsistent terminology, creating a formidable challenge for data integration. Our extensive literature mining and data reconciliation identified 122 neuron types based on neurotransmitter, axonal and dendritic patterns, synaptic specificity, electrophysiology, and molecular biomarkers. All ∼3700 annotated properties are individually supported by specific evidence (∼14,000 pieces) in peer-reviewed publications. Systematic analysis of this unprecedented amount of machine-readable information reveals novel correlations among neuron types and properties, the potential connectivity of the full hippocampal circuitry, and outstanding knowledge gaps. User-friendly browsing and online querying of Hippocampome.org may aid design and interpretation of both experiments and simulations. This powerful, simple, and extensible neuron classification endeavor is unique in its detail, utility, and completeness.

Bivariate and Multivariate NeuroXidence: A Robust and Reliable Method to Detect Modulations of Spike-Spike Synchronization Across Experimental Conditions.

Synchronous neuronal firing has been proposed as a potential neuronal code. To determine whether synchronous firing is really involved in different forms of information processing, one needs to directly compare the amount of synchronous firing due to various factors, such as different experimental or behavioral conditions. In order to address this issue, we present an extended version of the previously published method, NeuroXidence. The improved method incorporates bi- and multivariate testing to determine whether different factors result in synchronous firing occurring above the chance level. We demonstrate through the use of simulated data sets that bi- and multivariate NeuroXidence reliably and robustly detects joint-spike-events across different factors.

NeuroXidence: reliable and efficient analysis of an excess or deficiency of joint-spike events.

We present a non-parametric and computationally efficient method named NeuroXidence that detects coordinated firing of two or more neurons and tests whether the observed level of coordinated firing is significantly different from that expected by chance. The method considers the full auto-structure of the data, including the changes in the rate responses and the history dependencies in the spiking activity. Also, the method accounts for trial-by-trial variability in the dataset, such as the variability of the rate responses and their latencies. NeuroXidence can be applied to short data windows lasting only tens of milliseconds, which enables the tracking of transient neuronal states correlated to information processing. We demonstrate, on both simulated data and single-unit activity recorded in cat visual cortex, that NeuroXidence discriminates reliably between significant and spurious events that occur by chance.

Behavioral performance modulates spike field coherence in monkey prefrontal cortex.

To investigate neuronal processing during short-term memory, we analyzed behavior-related modulations of coupling between signals on two spatial scales: first, very local multiunit activity and second, local field potentials. Coupling was assessed by spike field coherence using a new approach to overcome limitations in cases of low firing rates. We demonstrate the reliability of our approach with simulated data. Application to recordings in prefrontal cortex of two monkeys revealed that locking of spikes was differentially modulated with two different frequency bands depending on behavioral performance.

Multiparameter estimation using only a chaotic time series and its applications.

An important extension to the techniques of synchronization-based parameter estimation is presented. Based on adaptive chaos synchronization, several methods are proposed to dynamically estimate multiple parameters using only a scalar chaotic time series. In comparison with previous schemes, the presented methods decrease the cost of parameter estimation and are more applicable in practice. Numerical examples are used to demonstrate the effectiveness and robustness of the presented methods. As an example application, an implementation of multichannel digital communication is proposed, where multiparameter modulation is used to simultaneously transmit more than one digital message. From a theoretical perspective, such an encoding increases the difficulty to directly read out the message from the transmitted signal and decreases the implementation cost.