Towards the automatic classification of neurons

Automated Analysis of Neuronal Morphology in 2D Fluorescence Micrographs through an Unsupervised Semantic Segmentation of Neurons

Brain function emerges from a highly complex network of specialized cells that are interlinked by billions of synapses. The synaptic connectivity between neurons is established between the elongated processes of their axons and dendrites or, together, neurites. To establish these connections, cellular neurites have to grow in highly specialized, cell-type dependent patterns covering extensive distances and connecting with thousands of other neurons. The outgrowth and branching of neurites are tightly controlled during development and are a commonly used functional readout of imaging in the neurosciences. Manual analysis of neuronal morphology from microscopy images, however, is very time intensive and prone to bias. Most automated analyses of neurons rely on reconstruction of the neuron as a whole without a semantic analysis of each neurite. A fully-automated classification of all neurites still remains unavailable in open-source software. Here we present a standalone, GUI-based software for batch-quantification of neuronal morphology in two-dimensional fluorescence micrographs of cultured neurons with minimal requirements for user interaction. Single neurons are first reconstructed into binarized images using a Hessian-based segmentation algorithm to detect thin neurite structures combined with intensity- and shape-based reconstruction of the cell body. Neurites are then classified into axon, dendrites and their branches of increasing order using a geodesic distance transform of the cell skeleton. The software was benchmarked against a published dataset and reproduced the phenotype observed after manual annotation. Our tool promises accelerated and improved morphometric studies of neuronal morphology by allowing for consistent and automated analysis of large datasets.

Automatical Spike Sorting With Low-Rank and Sparse Representation

Spikesorting is crucial in studying neural individually and synergistically encoding and decoding behaviors. However, existent spike sorting algorithms perform unsatisfactorily in real scenarios where heavy noises and overlapping samples are commonly in the spikes, and the spikes from different neurons are similar. To address such challenging scenarios, we propose an automatic spike sporting method in this paper, which integrally combines low-rank and sparse representation (LRSR) into a unified model. In particular, LRSR models spikes through low-rank optimization, uncovering global data structure for handling similar and overlapped samples. To eliminate the influence of the embedded noises, LRSR uses a sparse constraint, effectively separating spikes from noise. The optimization is solved using alternate augmented Lagrange multipliers methods. Moreover, we conclude with an automatic spike-sorting framework that employs the spectral clustering theorem to estimate the number of neurons. Extensive experiments over various simulated and real-world datasets demonstrate that our proposed method, LRSR, can handle spike sorting effectively and efficiently.

Graph Representation Learning for Large-Scale Neuronal Morphological Analysis

The analysis of neuronal morphological data is essential to investigate the neuronal properties and brain mechanisms. The complex morphologies, absence of annotations, and sheer volume of these data pose significant challenges in neuronal morphological analysis, such as identifying neuron types and large-scale neuron retrieval, all of which require accurate measuring and efficient matching algorithms. Recently, many studies have been conducted to describe neuronal morphologies quantitatively using predefined measurements. However, hand-crafted features are usually inadequate for distinguishing fine-grained differences among massive neurons. In this article, we propose a novel morphology-aware contrastive graph neural network (MACGNN) for unsupervised neuronal morphological representation learning. To improve the retrieval efficiency in large-scale neuronal morphological datasets, we further propose Hash-MACGNN by introducing an improved deep hash algorithm to train the network end-to-end to learn binary hash representations of neurons. We conduct extensive experiments on the largest dataset, NeuroMorpho, which contains more than 100 000 neurons. The experimental results demonstrate the effectiveness and superiority of our MACGNN and Hash-MACGNN for large-scale neuronal morphological analysis.

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.

DSM: Deep sequential model for complete neuronal morphology representation and feature extraction

The full morphology of single neurons is indispensable for understanding cell types, the basic building blocks in brains. Projecting trajectories are critical to extracting biologically relevant information from neuron morphologies, as they provide valuable information for both connectivity and cell identity. We developed an artificial intelligence method, deep sequential model (DSM), to extract concise, cell-type-defining features from projections across brain regions. DSM achieves more than 90% accuracy in classifying 12 major neuron projection types without compromising performance when spatial noise is present. Such remarkable robustness enabled us to efficiently manage and analyze several major full-morphology data sources, showcasing how characteristic long projections can define cell identities. We also succeeded in applying our model to both discovering previously unknown neuron subtypes and analyzing exceptional co-expressed genes involved in neuron projection circuits.

Betz cells of the primary motor cortex

Betz cells, named in honor of Volodymyr Betz (1834-1894), who described them as "giant pyramids" in the primary motor cortex of primates and other mammalian species, are layer V extratelencephalic projection (ETP) neurons that directly innervate α-motoneurons of the brainstem and spinal cord. Despite their large volume and circumferential dendritic architecture, to date, no single molecular criterion has been established that unequivocally distinguishes adult Betz cells from other layer V ETP neurons. In primates, transcriptional signatures suggest the presence of at least two ETP neuron clusters that contain mature Betz cells; these are characterized by an abundance of axon guidance and oxidative phosphorylation transcripts. How neurodevelopmental programs drive the distinct positional and morphological features of Betz cells in humans remains unknown. Betz cells display a distinct biphasic firing pattern involving early cessation of firing followed by delayed sustained acceleration in spike frequency and magnitude. Few cell type-specific transcripts and electrophysiological characteristics are conserved between rodent layer V ETP neurons of the motor cortex and primate Betz cells. This has implications for the modeling of disorders that affect the motor cortex in humans, such as amyotrophic lateral sclerosis (ALS). Perhaps vulnerability to ALS is linked to the evolution of neural networks for fine motor control reflected in the distinct morphomolecular architecture of the human motor cortex, including Betz cells. Here, we discuss histological, molecular, and functional data concerning the position of Betz cells in the emerging taxonomy of neurons across diverse species and their role in neurological disorders.

Sensitivity of CNN image analysis to multifaceted measurements of neurite growth

Quantitative analysis of neurite growth and morphology is essential for understanding the determinants of neural development and regeneration, however, it is complicated by the labor-intensive process of measuring diverse parameters of neurite outgrowth. Consequently, automated approaches have been developed to study neurite morphology in a high-throughput and comprehensive manner. These approaches include computer-automated algorithms known as 'convolutional neural networks' (CNNs)-powerful models capable of learning complex tasks without the biases of hand-crafted models. Nevertheless, their complexity often relegates them to functioning as 'black boxes.' Therefore, research in the field of explainable AI is imperative to comprehend the relationship between CNN image analysis output and predefined morphological parameters of neurite growth in order to assess the applicability of these machine learning approaches. In this study, drawing inspiration from the field of automated feature selection, we investigate the correlation between quantified metrics of neurite morphology and the image analysis results from NeuriteNet-a CNN developed to analyze neurite growth. NeuriteNet accurately distinguishes images of neurite growth based on different treatment groups within two separate experimental systems. These systems differentiate between neurons cultured on different substrate conditions and neurons subjected to drug treatment inhibiting neurite outgrowth. By examining the model's function and patterns of activation underlying its classification decisions, we discover that NeuriteNet focuses on aspects of neuron morphology that represent quantifiable metrics distinguishing these groups. Additionally, it incorporates factors that are not encompassed by neuron morphology tracing analyses. NeuriteNet presents a novel tool ideally suited for screening morphological differences in heterogeneous neuron groups while also providing impetus for targeted follow-up studies.

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.

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.

Correlation Analysis of Molecularly-Defined Cortical Interneuron Populations with Morpho-Electric Properties in Layer V of Mouse Neocortex

Cortical interneurons can be categorized into distinct populations based on multiple modalities, including molecular signatures and morpho-electrical (M/E) properties. Recently, many transcriptomic signatures based on single-cell RNA-seq have been identified in cortical interneurons. However, whether different interneuron populations defined by transcriptomic signature expressions correspond to distinct M/E subtypes is still unknown. Here, we applied the Patch-PCR approach to simultaneously obtain the M/E properties and messenger RNA (mRNA) expression of >600 interneurons in layer V of the mouse somatosensory cortex (S1). Subsequently, we identified 11 M/E subtypes, 9 neurochemical cell populations (NCs), and 20 transcriptomic cell populations (TCs) in this cortical lamina. Further analysis revealed that cells in many NCs and TCs comprised several M/E types and were difficult to clearly distinguish morpho-electrically. A similar analysis of layer V interneurons of mouse primary visual cortex (V1) and motor cortex (M1) gave results largely comparable to S1. Comparison between S1, V1, and M1 suggested that, compared to V1, S1 interneurons were morpho-electrically more similar to M1. Our study reveals the presence of substantial M/E variations in cortical interneuron populations defined by molecular expression.

Machine learning classification reveals robust morphometric biomarker of glial and neuronal arbors

Neurons and glia are the two main cell classes in the nervous systems of most animals. Although functionally distinct, neurons and glia are both characterized by multiple branching arbors stemming from the cell bodies. Glial processes are generally known to form smaller trees than neuronal dendrites. However, the full extent of morphological differences between neurons and glia in multiple species and brain regions has not yet been characterized, nor is it known whether these cells can be reliably distinguished based on geometric features alone. Here, we show that multiple supervised learning algorithms deployed on a large database of morphological reconstructions can systematically classify neuronal and glial arbors with nearly perfect accuracy and precision. Moreover, we report multiple morphometric properties, both size related and size independent, that differ substantially between these cell types. In particular, we newly identify an individual morphometric measurement, Average Branch Euclidean Length that can robustly separate neurons from glia across multiple animal models, a broad diversity of experimental conditions, and anatomical areas, with the notable exception of the cerebellum. We discuss the practical utility and physiological interpretation of this discovery.

Cell types as species: Exploring a metaphor

The concept of "cell type," though fundamental to cell biology, is controversial. Cells have historically been classified into types based on morphology, physiology, or location. More recently, single cell transcriptomic studies have revealed fine-scale differences among cells with similar gross phenotypes. Transcriptomic snapshots of cells at various stages of differentiation, and of cells under different physiological conditions, have shown that in many cases variation is more continuous than discrete, raising questions about the relationship between cell type and cell state. Some researchers have rejected the notion of fixed types altogether. Throughout the history of discussions on cell type, cell biologists have compared the problem of defining cell type with the interminable and often contentious debate over the definition of arguably the most important concept in systematics and evolutionary biology, "species." In the last decades, systematics, like cell biology, has been transformed by the increasing availability of molecular data, and the fine-grained resolution of genetic relationships have generated new ideas about how that variation should be classified. There are numerous parallels between the two fields that make exploration of the "cell types as species" metaphor timely. These parallels begin with philosophy, with discussion of both cell types and species as being either individuals, groups, or something in between (e.g., homeostatic property clusters). In each field there are various different types of lineages that form trees or networks that can (and in some cases do) provide criteria for grouping. Developing and refining models for evolutionary divergence of species and for cell type differentiation are parallel goals of the two fields. The goal of this essay is to highlight such parallels with the hope of inspiring biologists in both fields to look for new solutions to similar problems outside of their own field.

Identification of adult spinal Shox2 neuronal subpopulations based on unbiased computational clustering of electrophysiological properties

Spinal cord neurons integrate sensory and descending information to produce motor output. The expression of transcription factors has been used to dissect out the neuronal components of circuits underlying behaviors. However, most of the canonical populations of interneurons are heterogeneous and require additional criteria to determine functional subpopulations. Neurons expressing the transcription factor Shox2 can be subclassified based on the co-expression of the transcription factor Chx10 and each subpopulation is proposed to have a distinct connectivity and different role in locomotion. Adult Shox2 neurons have recently been shown to be diverse based on their firing properties. Here, in order to subclassify adult mouse Shox2 neurons, we performed multiple analyses of data collected from whole-cell patch clamp recordings of visually-identified Shox2 neurons from lumbar spinal slices. A smaller set of Chx10 neurons was included in the analyses for validation. We performed k-means and hierarchical unbiased clustering approaches, considering electrophysiological variables. Unlike the categorizations by firing type, the clusters displayed electrophysiological properties that could differentiate between clusters of Shox2 neurons. The presence of clusters consisting exclusively of Shox2 neurons in both clustering techniques suggests that it is possible to distinguish Shox2+Chx10- neurons from Shox2+Chx10+ neurons by electrophysiological properties alone. Computational clusters were further validated by immunohistochemistry with accuracy in a small subset of neurons. Thus, unbiased cluster analysis using electrophysiological properties is a tool that can enhance current interneuronal subclassifications and can complement groupings based on transcription factor and molecular expression.

RealNeuralNetworks.jl: An Integrated Julia Package for Skeletonization, Morphological Analysis, and Synaptic Connectivity Analysis of Terabyte-Scale 3D Neural Segmentations

Benefiting from the rapid development of electron microscopy imaging and deep learning technologies, an increasing number of brain image datasets with segmentation and synapse detection are published. Most of the automated segmentation methods label voxels rather than producing neuron skeletons directly. A further skeletonization step is necessary for quantitative morphological analysis. Currently, several tools are published for skeletonization as well as morphological and synaptic connectivity analysis using different computer languages and environments. Recently the Julia programming language, notable for elegant syntax and high performance, has gained rapid adoption in the scientific computing community. Here, we present a Julia package, called RealNeuralNetworks.jl, for efficient sparse skeletonization, morphological analysis, and synaptic connectivity analysis. Based on a large-scale Zebrafish segmentation dataset, we illustrate the software features by performing distributed skeletonization in Google Cloud, clustering the neurons using the NBLAST algorithm, combining morphological similarity and synaptic connectivity to study their relationship. We demonstrate that RealNeuralNetworks.jl is suitable for use in terabyte-scale electron microscopy image segmentation datasets.

Efficient metadata mining of web-accessible neural morphologies

Advancements in neuroscience research have led to steadily accelerating data production and sharing. The online community repository of neural reconstructions NeuroMorpho.Org grew from fewer than 1000 digitally traced neurons in 2006 to more than 140,000 cells today, including glia that now constitute 10.1% of the content. Every reconstruction consists of a detailed 3D representation of branch geometry and connectivity in a standardized format, from which a collection of morphometric features is extracted and stored. Moreover, each entry in the database is accompanied by rich metadata annotation describing the animal subject, anatomy, and experimental details. The rapid expansion of this resource in the past decade was accompanied by a parallel rise in the complexity of the available information, creating both opportunities and challenges for knowledge mining. Here, we introduce a new summary reporting functionality, allowing NeuroMorpho.Org users to efficiently download digests of metadata and morphometrics from multiple groups of similar cells for further analysis. We demonstrate the capabilities of the tool for both glia and neurons and present an illustrative statistical analysis of the resulting data.

Neuronal classification from network connectivity via adjacency spectral embedding

This work presents a novel strategy for classifying neurons, represented by nodes of a directed graph, based on their circuitry (edge connectivity). We assume a stochastic block model (SBM) in which neurons belong together if they connect to neurons of other groups according to the same probability distributions. Following adjacency spectral embedding of the SBM graph, we derive the number of classes and assign each neuron to a class with a Gaussian mixture model-based expectation maximization (EM) clustering algorithm. To improve accuracy, we introduce a simple variation using random hierarchical agglomerative clustering to initialize the EM algorithm and picking the best solution over multiple EM restarts. We test this procedure on a large (≈212-215 neurons), sparse, biologically inspired connectome with eight neuron classes. The simulation results demonstrate that the proposed approach is broadly stable to the choice of embedding dimension, and scales extremely well as the number of neurons in the network increases. Clustering accuracy is robust to variations in model parameters and highly tolerant to simulated experimental noise, achieving perfect classifications with up to 40% of swapped edges. Thus, this approach may be useful to analyze and interpret large-scale brain connectomics data in terms of underlying cellular components.

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.

Neuron type classification in rat brain based on integrative convolutional and tree-based recurrent neural networks

The study of cellular complexity in the nervous system based on anatomy has shown more practical and objective advantages in morphology than other perspectives on molecular, physiological, and evolutionary aspects. However, morphology-based neuron type classification in the whole rat brain is challenging, given the significant number of neuron types, limited reconstructed neuron samples, and diverse data formats. Here, we report that different types of deep neural network modules may well process different kinds of features and that the integration of these submodules will show power on the representation and classification of neuron types. For SWC-format data, which are compressed but unstructured, we construct a tree-based recurrent neural network (Tree-RNN) module. For 2D or 3D slice-format data, which are structured but with large volumes of pixels, we construct a convolutional neural network (CNN) module. We also generate a virtually simulated dataset with two classes, reconstruct a CASIA rat-neuron dataset with 2.6 million neurons without labels, and select the NeuroMorpho-rat dataset with 35,000 neurons containing hierarchical labels. In the twelve-class classification task, the proposed model achieves state-of-the-art performance compared with other models, e.g., the CNN, RNN, and support vector machine based on hand-designed features.

A community-based transcriptomics classification and nomenclature of neocortical cell types

To understand the function of cortical circuits, it is necessary to catalog their cellular diversity. Past attempts to do so using anatomical, physiological or molecular features of cortical cells have not resulted in a unified taxonomy of neuronal or glial cell types, partly due to limited data. Single-cell transcriptomics is enabling, for the first time, systematic high-throughput measurements of cortical cells and generation of datasets that hold the promise of being complete, accurate and permanent. Statistical analyses of these data reveal clusters that often correspond to cell types previously defined by morphological or physiological criteria and that appear conserved across cortical areas and species. To capitalize on these new methods, we propose the adoption of a transcriptome-based taxonomy of cell types for mammalian neocortex. This classification should be hierarchical and use a standardized nomenclature. It should be based on a probabilistic definition of a cell type and incorporate data from different approaches, developmental stages and species. A community-based classification and data aggregation model, such as a knowledge graph, could provide a common foundation for the study of cortical circuits. This community-based classification, nomenclature and data aggregation could serve as an example for cell type atlases in other parts of the body.

The landscape of regulatory genes in brain-wide neuronal phenotypes of a vertebrate brain

Multidimensional landscapes of regulatory genes in neuronal phenotypes at whole-brain levels in the vertebrate remain elusive. We generated single-cell transcriptomes of ~67,000 region- and neurotransmitter/neuromodulator-identifiable cells from larval zebrafish brains. Hierarchical clustering based on effector gene profiles ('terminal features') distinguished major brain cell types. Sister clusters at hierarchical termini displayed similar terminal features. It was further verified by a population-level statistical method. Intriguingly, glutamatergic/GABAergic sister clusters mostly expressed distinct transcription factor (TF) profiles ('convergent pattern'), whereas neuromodulator-type sister clusters predominantly expressed the same TF profiles ('matched pattern'). Interestingly, glutamatergic/GABAergic clusters with similar TF profiles could also display different terminal features ('divergent pattern'). It led us to identify a library of RNA-binding proteins that differentially marked divergent pair clusters, suggesting the post-transcriptional regulation of neuron diversification. Thus, our findings reveal multidimensional landscapes of transcriptional and post-transcriptional regulators in whole-brain neuronal phenotypes in the zebrafish brain.

Common cell type nomenclature for the mammalian brain

The advancement of single-cell RNA-sequencing technologies has led to an explosion of cell type definitions across multiple organs and organisms. While standards for data and metadata intake are arising, organization of cell types has largely been left to individual investigators, resulting in widely varying nomenclature and limited alignment between taxonomies. To facilitate cross-dataset comparison, the Allen Institute created the common cell type nomenclature (CCN) for matching and tracking cell types across studies that is qualitatively similar to gene transcript management across different genome builds. The CCN can be readily applied to new or established taxonomies and was applied herein to diverse cell type datasets derived from multiple quantifiable modalities. The CCN facilitates assigning accurate yet flexible cell type names in the mammalian cortex as a step toward community-wide efforts to organize multi-source, data-driven information related to cell type taxonomies from any organism.

Functional genomics, genetic risk profiling and cell phenotypes in neurodegenerative disease

Human genetics provides unbiased insights into the causes of human disease, which can be used to create a foundation for effective ways to more accurately diagnose patients, stratify patients for more successful clinical trials, discover and develop new therapies, and ultimately help patients choose the safest and most promising therapeutic option based on their risk profile. But the process for translating basic observations from human genetics studies into pathogenic disease mechanisms and treatments is laborious and complex, and this challenge has particularly slowed the development of interventions for neurodegenerative disease. In this review, we discuss the many steps in the process, the important considerations at each stage, and some of the latest tools and technologies that are available to help investigators translate insights from human genetics into diagnostic and therapeutic strategies that will lead to the sort of advances in clinical care that make a difference for patients.

Objective Morphological Classification of Neocortical Pyramidal Cells

A consensus on the number of morphologically different types of pyramidal cells (PCs) in the neocortex has not yet been reached, despite over a century of anatomical studies, due to the lack of agreement on the subjective classifications of neuron types, which is based on expert analyses of neuronal morphologies. Even for neurons that are visually distinguishable, there is no common ground to consistently define morphological types. The objective classification of PCs can be achieved with methods from algebraic topology, and the dendritic arborization is sufficient for the reliable identification of distinct types of cortical PCs. Therefore, we objectively identify 17 types of PCs in the rat somatosensory cortex. In addition, we provide a solution to the challenging problem of whether 2 similar neurons belong to different types or to a continuum of the same type. Our topological classification does not require expert input, is stable, and helps settle the long-standing debate on whether cell-types are discrete or continuous morphological variations of each other.

Operations Research Methods for Estimating the Population Size of Neuron Types

Understanding brain computation requires assembling a complete catalog of its architectural components. Although the brain is organized into several anatomical and functional regions, it is ultimately the neurons in every region that are responsible for cognition and behavior. Thus, classifying neuron types throughout the brain and quantifying the population sizes of distinct classes in different regions is a key subject of research in the neuroscience community. The total number of neurons in the brain has been estimated for multiple species, but the definition and population size of each neuron type are still open questions even in common model organisms: the so called "cell census" problem. We propose a methodology that uses operations research principles to estimate the number of neurons in each type based on available information on their distinguishing properties. Thus, assuming a set of neuron type definitions, we provide a solution to the issue of assessing their relative proportions. Specifically, we present a three-step approach that includes literature search, equation generation, and numerical optimization. Solving computationally the set of equations generated by literature mining yields best estimates or most likely ranges for the number of neurons in each type. While this strategy can be applied towards any neural system, we illustrate its usage on the rodent hippocampus.

Axonal Tree Morphology and Signal Propagation Dynamics Improve Interneuron Classification

Neurons are diverse and can be differentiated by their morphological, electrophysiological, and molecular properties. Current morphology-based classification approaches largely rely on the dendritic tree structure or on the overall axonal projection layout. Here, we use data from public databases of neuronal reconstructions and membrane properties to study the characteristics of the axonal and dendritic trees for interneuron classification. We show that combining signal propagation patterns observed by biophysical simulations of the activity along ramified axonal trees with morphological parameters of the axonal and dendritic trees, significantly improve classification results compared to previous approaches. The classification schemes introduced here can be utilized for robust neuronal classification. Our work paves the way for understanding and utilizing form-function principles in realistic neuronal reconstructions.

Cluster tendency assessment in neuronal spike data

Sorting spikes from extracellular recording into clusters associated with distinct single units (putative neurons) is a fundamental step in analyzing neuronal populations. Such spike sorting is intrinsically unsupervised, as the number of neurons are not known a priori. Therefor, any spike sorting is an unsupervised learning problem that requires either of the two approaches: specification of a fixed value k for the number of clusters to seek, or generation of candidate partitions for several possible values of c, followed by selection of a best candidate based on various post-clustering validation criteria. In this paper, we investigate the first approach and evaluate the utility of several methods for providing lower dimensional visualization of the cluster structure and on subsequent spike clustering. We also introduce a visualization technique called improved visual assessment of cluster tendency (iVAT) to estimate possible cluster structures in data without the need for dimensionality reduction. Experimental results are conducted on two datasets with ground truth labels. In data with a relatively small number of clusters, iVAT is beneficial in estimating the number of clusters to inform the initialization of clustering algorithms. With larger numbers of clusters, iVAT gives a useful estimate of the coarse cluster structure but sometimes fails to indicate the presumptive number of clusters. We show that noise associated with recording extracellular neuronal potentials can disrupt computational clustering schemes, highlighting the benefit of probabilistic clustering models. Our results show that t-Distributed Stochastic Neighbor Embedding (t-SNE) provides representations of the data that yield more accurate visualization of potential cluster structure to inform the clustering stage. Moreover, The clusters obtained using t-SNE features were more reliable than the clusters obtained using the other methods, which indicates that t-SNE can potentially be used for both visualization and to extract features to be used by any clustering algorithm.

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.

Morphological Neuron Classification Based on Dendritic Tree Hierarchy

The shape of a neuron can reveal interesting properties about its function. Therefore, morphological neuron characterization can contribute to a better understanding of how the brain works. However, one of the great challenges of neuroanatomy is the definition of morphological properties that can be used for categorizing neurons. This paper proposes a new methodology for neuron morphological analysis by considering different hierarchies of the dendritic tree for characterizing and categorizing neuronal cells. The methodology consists in using different strategies for decomposing the dendritic tree along its hierarchies, allowing the identification of relevant parts (possibly related to specific neuronal functions) for classification tasks. A set of more than 5000 neurons corresponding to 10 classes were examined with supervised classification algorithms based on this strategy. It was found that classification accuracies similar to those obtained by using whole neurons can be achieved by considering only parts of the neurons. Branches close to the soma were found to be particularly relevant for classification.

A Brief History of Simulation Neuroscience

Our knowledge of the brain has evolved over millennia in philosophical, experimental and theoretical phases. We suggest that the next phase is simulation neuroscience. The main drivers of simulation neuroscience are big data generated at multiple levels of brain organization and the need to integrate these data to trace the causal chain of interactions within and across all these levels. Simulation neuroscience is currently the only methodology for systematically approaching the multiscale brain. In this review, we attempt to reconstruct the deep historical paths leading to simulation neuroscience, from the first observations of the nerve cell to modern efforts to digitally reconstruct and simulate the brain. Neuroscience began with the identification of the neuron as the fundamental unit of brain structure and function and has evolved towards understanding the role of each cell type in the brain, how brain cells are connected to each other, and how the seemingly infinite networks they form give rise to the vast diversity of brain functions. Neuronal mapping is evolving from subjective descriptions of cell types towards objective classes, subclasses and types. Connectivity mapping is evolving from loose topographic maps between brain regions towards dense anatomical and physiological maps of connections between individual genetically distinct neurons. Functional mapping is evolving from psychological and behavioral stereotypes towards a map of behaviors emerging from structural and functional connectomes. We show how industrialization of neuroscience and the resulting large disconnected datasets are generating demand for integrative neuroscience, how the scale of neuronal and connectivity maps is driving digital atlasing and digital reconstruction to piece together the multiple levels of brain organization, and how the complexity of the interactions between molecules, neurons, microcircuits and brain regions is driving brain simulation to understand the interactions in the multiscale brain.

Generative models and abstractions for large-scale neuroanatomy datasets

Neural datasets are increasing rapidly in both resolution and volume. In neuroanatomy, this trend has been accelerated by innovations in imaging technology. As full datasets are impractical and unnecessary for many applications, it is important to identify abstractions that distill useful features of neural structure, organization, and anatomy. In this review article, we discuss several such abstractions and highlight recent algorithmic advances in working with these models. In particular, we discuss the use of generative models in neuroanatomy; such models may be considered 'meta-abstractions' that capture distributions over other abstractions.

The roles of supervised machine learning in systems neuroscience

Over the last several years, the use of machine learning (ML) in neuroscience has been rapidly increasing. Here, we review ML's contributions, both realized and potential, across several areas of systems neuroscience. We describe four primary roles of ML within neuroscience: (1) creating solutions to engineering problems, (2) identifying predictive variables, (3) setting benchmarks for simple models of the brain, and (4) serving itself as a model for the brain. The breadth and ease of its applicability suggests that machine learning should be in the toolbox of most systems neuroscientists.

Neuronal cell types in the fly: single-cell anatomy meets single-cell genomics

At around 150 000 neurons, the adult Drosophila melanogaster central nervous system is one of the largest species, for which a complete cellular catalogue is imminent. While numerically much simpler than mammalian brains, its complexity is still difficult to parse without grouping neurons into consistent types, which can number 1-1000 cells per hemisphere. We review how neuroanatomical and gene expression data are being used to discover neuronal types at scale. The correlation among multiple co-varying neuronal properties, including lineage, gene expression, morphology, connectivity, response properties and shared behavioral significance is essential to the definition of neuronal cell type. Initial studies comparing morphological and transcriptomic definitions of neuronal type suggest that these are highly consistent, but there is much to do to match these approaches brain-wide. Matched single-cell transcriptomic and morphological data provide an effective reference point to integrate other data types, including connectomics data. This will significantly enhance our ability to make functional predictions from brain wiring diagrams as well facilitating molecular genetic manipulation of neuronal types.

Towards a supervised classification of neocortical interneuron morphologies

BACKGROUND

The challenge of classifying cortical interneurons is yet to be solved. Data-driven classification into established morphological types may provide insight and practical value.

RESULTS

We trained models using 217 high-quality morphologies of rat somatosensory neocortex interneurons reconstructed by a single laboratory and pre-classified into eight types. We quantified 103 axonal and dendritic morphometrics, including novel ones that capture features such as arbor orientation, extent in layer one, and dendritic polarity. We trained a one-versus-rest classifier for each type, combining well-known supervised classification algorithms with feature selection and over- and under-sampling. We accurately classified the nest basket, Martinotti, and basket cell types with the Martinotti model outperforming 39 out of 42 leading neuroscientists. We had moderate accuracy for the double bouquet, small and large basket types, and limited accuracy for the chandelier and bitufted types. We characterized the types with interpretable models or with up to ten morphometrics.

CONCLUSION

Except for large basket, 50 high-quality reconstructions sufficed to learn an accurate model of a type. Improving these models may require quantifying complex arborization patterns and finding correlates of bouton-related features. Our study brings attention to practical aspects important for neuron classification and is readily reproducible, with all code and data available online.

Characterization of Retinal Functionality at Different Eccentricities in a Diurnal Rodent

Although the properties of the neurons of the visual system that process central and peripheral regions of the visual field have been widely researched in the visual cortex and the LGN, they have scarcely been documented for the retina. The retina is the first step in integrating optical signals, and despite considerable efforts to functionally characterize the different types of retinal ganglion cells (RGCs), a clear account of the particular functionality of cells with central vs. peripheral fields is still wanting. Here, we use electrophysiological recordings, gathered from retinas of the diurnal rodent Octodon degus, to show that RGCs with peripheral receptive fields (RF) are larger, faster, and have shorter transient responses. This translates into higher sensitivity at high temporal frequencies and a full frequency bandwidth when compared to RGCs with more central RF. We also observed that imbalances between ON and OFF cell populations are preserved with eccentricity. Finally, the high diversity of functional types of RGCs highlights the complexity of the computational strategies implemented in the early stages of visual processing, which could inspire the development of bio-inspired artificial systems.

Electrophysiological Profiling of Neocortical Neural Subtypes: A Semi-Supervised Method Applied to in vivo Whole-Cell Patch-Clamp Data

A lot of efforts have been made to understand the structure and function of neocortical circuits. In fact, a promising way to understand the functions of cortical circuits is the classification of the neural types, based on their different properties. Recent studies focused on applying modern computational methods to classify neurons based on molecular, morphological, physiological, or mixed of these criteria. Although there are studies in the literature on in vitro/vivo extracellular or in vitro intracellular recordings, a study on the classification of neuronal types using in vivo whole-cell patch-clamp recordings is still lacking. We thus proposed a novel semi-supervised classification method based on waveform shape of neurons' spikes using in vivo whole-cell patch-clamp recordings. We, first, detected spike candidates. Then discriminative features were extracted from the time samples of the spikes using discrete cosine transform. We then extracted the center of clusters using fuzzy c-mean clustering and finally, the neurons were classified using the minimum distance classifier. We distinguished three types of neurons: excitatory pyramidal cells (Pyr) and two types of inhibitory neurons: GABAergic- parvalbumin positive (PV), and somatostatin positive (SST) non-pyramidal cells in layer II/III of the mice primary visual cortex. We used 10-fold cross validation in our study. The classification accuracy for PV, Pyr, and SST was 91.59 ± 1.69, 97.47 ± 0.67, and 89.06 ± 1.99, respectively. Overall, the algorithm correctly classified 92.67 ± 0.54% of the cells, confirming the relative robustness of the discriminant functions. The performance of the method was further assessed on in vitro recordings by using a pool of 50 neurons from Allen institute Cell Types Database (5 major subtypes of neurons: Pyr, PV, SST, 5HT3a, and vasoactive intestinal peptide (VIP) cells). Its overall accuracy was 84.13 ± 0.81% on this data set using cross validation framework. The proposed algorithm is thus a promising new tool in recognizing cell's type with high accuracy in laboratories using in vivo/vitro whole-cell patch-clamp recording technique. The developed programs and the entire dataset are available online to interested readers.

Flexible Learning-Free Segmentation and Reconstruction of Neural Volumes

Imaging is a dominant strategy for data collection in neuroscience, yielding stacks of images that often scale to gigabytes of data for a single experiment. Machine learning algorithms from computer vision can serve as a pair of virtual eyes that tirelessly processes these images, automatically detecting and identifying microstructures. Unlike learning methods, our Flexible Learning-free Reconstruction of Imaged Neural volumes (FLoRIN) pipeline exploits structure-specific contextual clues and requires no training. This approach generalizes across different modalities, including serially-sectioned scanning electron microscopy (sSEM) of genetically labeled and contrast enhanced processes, spectral confocal reflectance (SCoRe) microscopy, and high-energy synchrotron X-ray microtomography (μCT) of large tissue volumes. We deploy the FLoRIN pipeline on newly published and novel mouse datasets, demonstrating the high biological fidelity of the pipeline’s reconstructions. FLoRIN reconstructions are of sufficient quality for preliminary biological study, for example examining the distribution and morphology of cells or extracting single axons from functional data. Compared to existing supervised learning methods, FLoRIN is one to two orders of magnitude faster and produces high-quality reconstructions that are tolerant to noise and artifacts, as is shown qualitatively and quantitatively.

Classification of GABAergic neuron subtypes from the globus pallidus using wild-type and transgenic mice

KEY POINTS

Classifying different subtypes of neurons in deep brain structures is a challenge and is crucial to better understand brain function. Understanding the diversity of neurons in the globus pallidus (GP), a brain region positioned to influence afferent and efferent information processing within basal ganglia, could help to explain a variety of brain functions. We present a classification of neurons from the GP using electrophysiological data from wild-type mice and confirmation using transgenic mice. This work will help researchers to identify specific neuronal subsets in the GP of wild-type mice when transgenic mice with labelled neurons are lacking.

ABSTRACT

Classification of the extensive neuronal diversity in the brain is fundamental for neuroscience. The globus pallidus external segment (GPe), also referred to as the globus pallidus in rodents, is a large nucleus located in the core of the basal ganglia whose circuitry is implicated in action control, decision-making and reward. Although considerable progress has been made in characterizing different GPe neuronal subtypes, no work has directly attempted to characterize these neurons in non-transgenic mice. Here, we provide data showing the degree of overlap in expression of neuronal PAS domain protein (Npas1), LIM homeobox 6 (Lhx6), parvalbumin (PV) and transcription factor FoxP2 biomarkers in mouse GPe neurons. We used an unbiased statistical method to classify neurons based on electrophysiological properties from nearly 200 neurons from C57BL/6J mice. In addition, we examined the subregion distribution of the neuronal subtypes. Cluster analysis using firing rate and hyperpolarization-induced membrane potential sag variables revealed three distinct neuronal clusters: type 1, characterized by low firing rate and small sag potential; type 2, with low firing rate and larger sag potential; and type 3, with high firing rate and small sag potential. We used other electrophysiological variables and data from marker-expressing neurons to evaluate the clusters. We propose that the GPe GABAergic neurons should be classified into three subgroups: arkypallidal, low-firing prototypical and high-firing prototypical neurons. This work will help researchers identify GPe neuron subtypes when transgenic mice with labelled neurons cannot be used.

Mining Big Neuron Morphological Data

The advent of automatic tracing and reconstruction technology has led to a surge in the number of neurons 3D reconstruction data and consequently the neuromorphology research. However, the lack of machine-driven annotation schema to automatically detect the types of the neurons based on their morphology still hinders the development of this branch of science. Neuromorphology is important because of the interplay between the shape and functionality of neurons and the far-reaching impact on the diagnostics and therapeutics in neurological disorders. This survey paper provides a comprehensive research in the field of automatic neurons classification and presents the existing challenges, methods, tools, and future directions for automatic neuromorphology analytics. We summarize the major automatic techniques applicable in the field and propose a systematic data processing pipeline for automatic neuron classification, covering data capturing, preprocessing, analyzing, classification, and retrieval. Various techniques and algorithms in machine learning are illustrated and compared to the same dataset to facilitate ongoing research in the field.

New Features for Neuron Classification

This paper addresses the problem of obtaining new neuron features capable of improving results of neuron classification. Most studies on neuron classification using morphological features have been based on Euclidean geometry. Here three one-dimensional (1D) time series are derived from the three-dimensional (3D) structure of neuron instead, and afterwards a spatial time series is finally constructed from which the features are calculated. Digitally reconstructed neurons were separated into control and pathological sets, which are related to three categories of alterations caused by epilepsy, Alzheimer's disease (long and local projections), and ischemia. These neuron sets were then subjected to supervised classification and the results were compared considering three sets of features: morphological, features obtained from the time series and a combination of both. The best results were obtained using features from the time series, which outperformed the classification using only morphological features, showing higher correct classification rates with differences of 5.15, 3.75, 5.33% for epilepsy and Alzheimer's disease (long and local projections) respectively. The morphological features were better for the ischemia set with a difference of 3.05%. Features like variance, Spearman auto-correlation, partial auto-correlation, mutual information, local minima and maxima, all related to the time series, exhibited the best performance. Also we compared different evaluators, among which ReliefF was the best ranked.

Distinct Functional Groups Emerge from the Intrinsic Properties of Molecularly Identified Entorhinal Interneurons and Principal Cells

Inhibitory interneurons are an important source of synaptic inputs that may contribute to network mechanisms for coding of spatial location by entorhinal cortex (EC). The intrinsic properties of inhibitory interneurons in the EC of the mouse are mostly undescribed. Intrinsic properties were recorded from known cell types, such as, stellate and pyramidal cells and 6 classes of molecularly identified interneurons (regulator of calcineurin 2, somatostatin, serotonin receptor 3a, neuropeptide Y neurogliaform (NGF), neuropeptide Y non-NGF, and vasoactive intestinal protein) in acute brain slices. We report a broad physiological diversity between and within cell classes. We also found differences in the ability to produce postinhibitory rebound spikes and in the frequency and amplitude of incoming EPSPs. To understand the source of this intrinsic variability we applied hierarchical cluster analysis to functionally classify neurons. These analyses revealed physiologically derived cell types in EC that mostly corresponded to the lines identified by biomarkers with a few unexpected and important differences. Finally, we reduced the complex multidimensional space of intrinsic properties to the most salient five that predicted the cellular biomolecular identity with 81.4% accuracy. These results provide a framework for the classification of functional subtypes of cortical neurons by their intrinsic membrane properties.

An open repository for single-cell reconstructions of the brain forest

NeuroMorpho.Org was launched in 2006 to provide unhindered access to any and all digital tracings of neuronal morphology that researchers were willing to share freely upon request. Today this database is the largest public inventory of cellular reconstructions in neuroscience with a content of over 80,000 neurons and glia from a representative diversity of animal species, anatomical regions, and experimental methods. Datasets continuously contributed by hundreds of laboratories worldwide are centrally curated, converted into a common non-proprietary format, morphometrically quantified, and annotated with comprehensive metadata. Users download digital reconstructions for a variety of scientific applications including visualization, classification, analysis, and simulations. With more than 1,000 peer-reviewed publications describing data stored in or utilizing data retrieved from NeuroMorpho.Org, this ever-growing repository can already be considered a mature resource for neuroscience.

Learning from connectomics on the fly

Parallels between invertebrates and vertebrates in nervous system development, organisation and circuits are powerful reasons to use insects to study the mechanistic basis of behaviour. The last few years have seen the generation in Drosophila melanogaster of very large light microscopy data sets, genetic driver lines and tools to report or manipulate neural activity. These resources in conjunction with computational tools are enabling large scale characterisation of neuronal types and their functional properties. These are complemented by 3D electron microscopy, providing synaptic resolution data. A whole brain connectome of the fly larva is approaching completion based on manual reconstruction of electron-microscopy data. An adult whole brain dataset is already publicly available and focussed reconstruction is under way, but its 40× greater volume would require ∼500-5000 person-years of manual labour. Nevertheless rapid technical improvements in imaging and especially automated segmentation will likely deliver a complete adult connectome in the next 5 years. To enhance our understanding of the circuit basis of behaviour, light and electron microscopy outputs must be integrated with functional and physiological information into comprehensive databases. We review presently available data, tools and opportunities in Drosophila. We then consider the limits and potential of future progress and how this may impact neuroscience in rich model systems provided by larger insects and vertebrates.

Linking neuronal structure to function in rodent hippocampus: a methodological prospective

Since the discovery of place cells, hippocampus-dependent spatial navigation has proven to be an ideal model system for resolving the relationship between neural coding and behavior. Electrical recordings from the hippocampal formation in freely moving animals have revealed a rich repertoire of spatial firing patterns and have enormously advanced our understanding of the neural principles of spatial representation. However, limited progress has been achieved in resolving the underlying cellular mechanisms. This is partially attributable to the inability of standard recording techniques to link neuronal structure to function directly. In this review, we summarize recent efforts aimed at filling this gap. We also highlight the development of methodologies that allow functional measurements from identified neuronal elements in behaving rodents. Recent progress in the dentate gyrus serves as a showcase to reveal the potential of such methodologies and the necessity of resolving structure-function relationships in order to access the cellular mechanisms of hippocampal circuit computations.

Transcriptional Architecture of Synaptic Communication Delineates GABAergic Neuron Identity

Understanding the organizational logic of neural circuits requires deciphering the biological basis of neuronal diversity and identity, but there is no consensus on how neuron types should be defined. We analyzed single-cell transcriptomes of a set of anatomically and physiologically characterized cortical GABAergic neurons and conducted a computational genomic screen for transcriptional profiles that distinguish them from one another. We discovered that cardinal GABAergic neuron types are delineated by a transcriptional architecture that encodes their synaptic communication patterns. This architecture comprises 6 categories of ∼40 gene families, including cell-adhesion molecules, transmitter-modulator receptors, ion channels, signaling proteins, neuropeptides and vesicular release components, and transcription factors. Combinatorial expression of select members across families shapes a multi-layered molecular scaffold along the cell membrane that may customize synaptic connectivity patterns and input-output signaling properties. This molecular genetic framework of neuronal identity integrates cell phenotypes along multiple axes and provides a foundation for discovering and classifying neuron types.

CA1 pyramidal cells have diverse biophysical properties, affected by development, experience, and aging

Neuron types (e.g., pyramidal cells) within one area of the brain are often considered homogeneous, despite variability in their biophysical properties. Here we review literature demonstrating variability in the electrical activity of CA1 hippocampal pyramidal cells (PCs), including responses to somatic current injection, synaptic stimulation, and spontaneous network-related activity. In addition, we describe how responses of CA1 PCs vary with development, experience, and aging, and some of the underlying ionic currents responsible. Finally, we suggest directions that may be the most impactful in expanding this knowledge, including the use of text and data mining to systematically study cellular heterogeneity in more depth; dynamical systems theory to understand and potentially classify neuron firing patterns; and mathematical modeling to study the interaction between cellular properties and network output. Our goals are to provide a synthesis of the literature for experimentalists studying CA1 PCs, to give theorists an idea of the rich diversity of behaviors models may need to reproduce to accurately represent these cells, and to provide suggestions for future research.

Metrics for comparing neuronal tree shapes based on persistent homology

As more and more neuroanatomical data are made available through efforts such as NeuroMorpho.Org and FlyCircuit.org, the need to develop computational tools to facilitate automatic knowledge discovery from such large datasets becomes more urgent. One fundamental question is how best to compare neuron structures, for instance to organize and classify large collection of neurons. We aim to develop a flexible yet powerful framework to support comparison and classification of large collection of neuron structures efficiently. Specifically we propose to use a topological persistence-based feature vectorization framework. Existing methods to vectorize a neuron (i.e, convert a neuron to a feature vector so as to support efficient comparison and/or searching) typically rely on statistics or summaries of morphometric information, such as the average or maximum local torque angle or partition asymmetry. These simple summaries have limited power in encoding global tree structures. Based on the concept of topological persistence recently developed in the field of computational topology, we vectorize each neuron structure into a simple yet informative summary. In particular, each type of information of interest can be represented as a descriptor function defined on the neuron tree, which is then mapped to a simple persistence-signature. Our framework can encode both local and global tree structure, as well as other information of interest (electrophysiological or dynamical measures), by considering multiple descriptor functions on the neuron. The resulting persistence-based signature is potentially more informative than simple statistical summaries (such as average/mean/max) of morphometric quantities-Indeed, we show that using a certain descriptor function will give a persistence-based signature containing strictly more information than the classical Sholl analysis. At the same time, our framework retains the efficiency associated with treating neurons as points in a simple Euclidean feature space, which would be important for constructing efficient searching or indexing structures over them. We present preliminary experimental results to demonstrate the effectiveness of our persistence-based neuronal feature vectorization framework.

An ontology-based search engine for digital reconstructions of neuronal morphology

Neuronal morphology is extremely diverse across and within animal species, developmental stages, brain regions, and cell types. This diversity is functionally important because neuronal structure strongly affects synaptic integration, spiking dynamics, and network connectivity. Digital reconstructions of axonal and dendritic arbors are thus essential to quantify and model information processing in the nervous system. NeuroMorpho.Org is an established repository containing tens of thousands of digitally reconstructed neurons shared by several hundred laboratories worldwide. Each neuron is annotated with specific metadata based on the published references and additional details provided by data owners. The number of represented metadata concepts has grown over the years in parallel with the increase of available data. Until now, however, the lack of standardized terminologies and of an adequately structured metadata schema limited the effectiveness of user searches. Here we present a new organization of NeuroMorpho.Org metadata grounded on a set of interconnected hierarchies focusing on the main dimensions of animal species, anatomical regions, and cell types. We have comprehensively mapped each metadata term in NeuroMorpho.Org to this formal ontology, explicitly resolving all ambiguities caused by synonymy and homonymy. Leveraging this consistent framework, we introduce OntoSearch, a powerful functionality that seamlessly enables retrieval of morphological data based on expert knowledge and logical inferences through an intuitive string-based user interface with auto-complete capability. In addition to returning the data directly matching the search criteria, OntoSearch also identifies a pool of possible hits by taking into consideration incomplete metadata annotation.