Metrics for comparing neuronal tree shapes based on persistent homology

Accelerating the continuous community sharing of digital neuromorphology data

The tree-like morphology of neurons and glia is a key cellular determinant of circuit connectivity and metabolic function in the nervous system of essentially all animals. To elucidate the contribution of specific cell types to both physiological and pathological brain states, it is important to access detailed neuroanatomy data for quantitative analysis and computational modeling. NeuroMorpho.Org is the largest online collection of freely available digital neural reconstructions and related metadata and is continuously updated with new uploads. Earlier in the project, we released multiple datasets together yearly, but this process caused an average delay of several months in making the data public. Moreover, in the past 5 years, >80% of invited authors agreed to share their data with the community via NeuroMorpho.Org, up from <20% in the first 5 years of the project. In the same period, the average number of reconstructions per publication increased 600%, creating the need for automatic processing to release more reconstructions in less time. The progressive automation of our pipeline enabled the transition to agile releases of individual datasets as soon as they are ready. The overall time from data identification to public sharing decreased by 63.7%; 78% of the datasets are now released in less than 3 months with an average workflow duration below 40 days. Furthermore, the mean processing time per reconstruction dropped from 3 h to 2 min. With these continuous improvements, NeuroMorpho.Org strives to forge a positive culture of open data. Most importantly, the new, original research enabled through reuse of datasets across the world has a multiplicative effect on science discovery, benefiting both authors and users.

Accelerating the continuous community sharing of digital neuromorphology data

The tree-like morphology of neurons and glia is a key cellular determinant of circuit connectivity and metabolic function in the nervous system of essentially all animals. To elucidate the contribution of specific cell types to both physiological and pathological brain states, it is important to access detailed neuroanatomy data for quantitative analysis and computational modeling. NeuroMorpho.Org is the largest online collection of freely available digital neural reconstructions and related metadata and is continuously updated with new uploads. Earlier in the project, we released multiple datasets together yearly, but this process caused an average delay of several months in making the data public. Moreover, in the past 5 years, >80% of invited authors agreed to share their data with the community via NeuroMorpho.Org, up from <20% in the first 5 years of the project. In the same period, the average number of reconstructions per publication increased 600%, creating the need for automatic processing to release more reconstructions in less time. The progressive automation of our pipeline enabled the transition to agile releases of individual datasets as soon as they are ready. The overall time from data identification to public sharing decreased by 63.7%; 78% of the datasets are now released in less than 3 months with an average workflow duration below 40 days. Furthermore, the mean processing time per reconstruction dropped from 3 hours to 2 minutes. With these continuous improvements, NeuroMorpho.Org strives to forge a positive culture of open data. Most importantly, the new, original research enabled through reuse of datasets across the world has a multiplicative effect on science discovery, benefiting both authors and users.

A robust approach to 3D neuron shape representation for quantification and classification

We consider the problem of finding an accurate representation of neuron shapes, extracting sub-cellular features, and classifying neurons based on neuron shapes. In neuroscience research, the skeleton representation is often used as a compact and abstract representation of neuron shapes. However, existing methods are limited to getting and analyzing "curve" skeletons which can only be applied for tubular shapes. This paper presents a 3D neuron morphology analysis method for more general and complex neuron shapes. First, we introduce the concept of skeleton mesh to represent general neuron shapes and propose a novel method for computing mesh representations from 3D surface point clouds. A skeleton graph is then obtained from skeleton mesh and is used to extract sub-cellular features. Finally, an unsupervised learning method is used to embed the skeleton graph for neuron classification. Extensive experiment results are provided and demonstrate the robustness of our method to analyze neuron morphology.

BigNeuron: a resource to benchmark and predict performance of algorithms for automated tracing of neurons in light microscopy datasets

BigNeuron is an open community bench-testing platform with the goal of setting open standards for accurate and fast automatic neuron tracing. We gathered a diverse set of image volumes across several species that is representative of the data obtained in many neuroscience laboratories interested in neuron tracing. Here, we report generated gold standard manual annotations for a subset of the available imaging datasets and quantified tracing quality for 35 automatic tracing algorithms. The goal of generating such a hand-curated diverse dataset is to advance the development of tracing algorithms and enable generalizable benchmarking. Together with image quality features, we pooled the data in an interactive web application that enables users and developers to perform principal component analysis, t-distributed stochastic neighbor embedding, correlation and clustering, visualization of imaging and tracing data, and benchmarking of automatic tracing algorithms in user-defined data subsets. The image quality metrics explain most of the variance in the data, followed by neuromorphological features related to neuron size. We observed that diverse algorithms can provide complementary information to obtain accurate results and developed a method to iteratively combine methods and generate consensus reconstructions. The consensus trees obtained provide estimates of the neuron structure ground truth that typically outperform single algorithms in noisy datasets. However, specific algorithms may outperform the consensus tree strategy in specific imaging conditions. Finally, to aid users in predicting the most accurate automatic tracing results without manual annotations for comparison, we used support vector machine regression to predict reconstruction quality given an image volume and a set of automatic tracings.

Local Microtubule and F-Actin Distributions Fully Constrain the Spatial Geometry of Drosophila Sensory Dendritic Arbors

Dendritic morphology underlies the source and processing of neuronal signal inputs. Morphology can be broadly described by two types of geometric characteristics. The first is dendrogram topology, defined by the length and frequency of the arbor branches; the second is spatial embedding, mainly determined by branch angles and straightness. We have previously demonstrated that microtubules and actin filaments are associated with arbor elongation and branching, fully constraining dendrogram topology. Here, we relate the local distribution of these two primary cytoskeletal components with dendritic spatial embedding. We first reconstruct and analyze 167 sensory neurons from the Drosophila larva encompassing multiple cell classes and genotypes. We observe that branches with a higher microtubule concentration tend to deviate less from the direction of their parent branch across all neuron types. Higher microtubule branches are also overall straighter. F-actin displays a similar effect on angular deviation and branch straightness, but not as consistently across all neuron types as microtubule. These observations raise the question as to whether the associations between cytoskeletal distributions and arbor geometry are sufficient constraints to reproduce type-specific dendritic architecture. Therefore, we create a computational model of dendritic morphology purely constrained by the cytoskeletal composition measured from real neurons. The model quantitatively captures both spatial embedding and dendrogram topology across all tested neuron groups. These results suggest a common developmental mechanism regulating diverse morphologies, where the local cytoskeletal distribution can fully specify the overall emergent geometry of dendritic arbors.

Local microtubule and F-actin distributions fully determine the spatial geometry of Drosophila sensory dendritic arbors

Dendritic morphology underlies the source and processing of neuronal signal inputs. Morphology can be broadly described by two types of geometric characteristics. The first is dendrogram topology, defined by the length and frequency of the arbor branches; the second is spatial embedding, mainly determined by branch angles and tortuosity. We have previously demonstrated that microtubules and actin filaments are associated with arbor elongation and branching, fully constraining dendrogram topology. Here we relate the local distribution of these two primary cytoskeletal components with dendritic spatial embedding. We first reconstruct and analyze 167 sensory neurons from the Drosophila larva encompassing multiple cell classes and genotypes. We observe that branches with higher microtubule concentration are overall straighter and tend to deviate less from the direction of their parent branch. F-actin displays a similar effect on the angular deviation from the parent branch direction, but its influence on branch tortuosity varies by class and genotype. We then create a computational model of dendritic morphology purely constrained by the cytoskeletal composition imaged from real neurons. The model quantitatively captures both spatial embedding and dendrogram topology across all tested neuron groups. These results suggest a common developmental mechanism regulating diverse morphologies, where the local cytoskeletal distribution can fully specify the overall emergent geometry of dendritic arbors.

Cell morphologies in the nervous system: Glia steal the limelight

Neurons and glia have distinct yet interactive functions but are both characterized by branching morphology. Dendritic trees have been digitally traced for over 40 years in many animal species, anatomical regions, and neuron types. Recently, long-range axons also are being reconstructed throughout the brain of many organisms from invertebrates to primates. In contrast, less attention has been paid until lately to glial morphology. Thus, although glia and neurons are similarly abundant in the nervous systems of humans and most animal models, glia have traditionally been much less represented than neurons in morphological reconstruction repositories such as NeuroMorpho.Org. This is rapidly changing with the advent of high-throughput glia tracing. NeuroMorpho.Org introduced glial cells in 2017 and today they constitute nearly a third of the database content. It took NeuroMorpho.Org 10 years to collect the first 40,000 neurons and now that amount of data can be produced in a single publication. This not only demonstrates the spectacular technological progress in data production, but also demands a corresponding advancement in informatics processing. At the same time, these publicly available data also open new opportunities for quantitative analysis and computational modeling to identify universal or cell-type-specific design principles in the cellular architecture of nervous systems. As a first application, we demonstrated that supervised machine learning of tree geometry classifies neurons and glia with practically perfect accuracy. Furthermore, we discovered a new morphometric biomarker capable of robustly separating these cell classes across multiple species, brain regions, and experimental preparations, with only sparse sampling of branch measurements.

Neuron tracing from light microscopy images: automation, deep learning and bench testing

MOTIVATION

Large-scale neuronal morphologies are essential to neuronal typing, connectivity characterization and brain modeling. It is widely accepted that automation is critical to the production of neuronal morphology. Despite previous survey papers about neuron tracing from light microscopy data in the last decade, thanks to the rapid development of the field, there is a need to update recent progress in a review focusing on new methods and remarkable applications.

RESULTS

This review outlines neuron tracing in various scenarios with the goal to help the community understand and navigate tools and resources. We describe the status, examples and accessibility of automatic neuron tracing. We survey recent advances of the increasingly popular deep-learning enhanced methods. We highlight the semi-automatic methods for single neuron tracing of mammalian whole brains as well as the resulting datasets, each containing thousands of full neuron morphologies. Finally, we exemplify the commonly used datasets and metrics for neuron tracing bench testing.

A tool for mapping microglial morphology, morphOMICs, reveals brain-region and sex-dependent phenotypes

Environmental cues influence the highly dynamic morphology of microglia. Strategies to characterize these changes usually involve user-selected morphometric features, which preclude the identification of a spectrum of context-dependent morphological phenotypes. Here we develop MorphOMICs, a topological data analysis approach, which enables semiautomatic mapping of microglial morphology into an atlas of cue-dependent phenotypes and overcomes feature-selection biases and biological variability. We extract spatially heterogeneous and sexually dimorphic morphological phenotypes for seven adult mouse brain regions. This sex-specific phenotype declines with maturation but increases over the disease trajectories in two neurodegeneration mouse models, with females showing a faster morphological shift in affected brain regions. Remarkably, microglia morphologies reflect an adaptation upon repeated exposure to ketamine anesthesia and do not recover to control morphologies. Finally, we demonstrate that both long primary processes and short terminal processes provide distinct insights to morphological phenotypes. MorphOMICs opens a new perspective to characterize microglial morphology.

Colombo et al. build a morphological spectrum of over 40,000 microglia across development and disease with a topological data analysis approach that allows mapping of new conditions along these sex-region-specific and brain-region-specific atlases.

Large scale similarity search across digital reconstructions of neural morphology

Most functions of the nervous system depend on neuronal and glial morphology. Continuous advances in microscopic imaging and tracing software have provided an increasingly abundant availability of 3D reconstructions of arborizing dendrites, axons, and processes, allowing their detailed study. However, efficient, large-scale methods to rank neural morphologies by similarity to an archetype are still lacking. Using the NeuroMorpho.Org database, we present a similarity search software enabling fast morphological comparison of hundreds of thousands of neural reconstructions from any species, brain regions, cell types, and preparation protocols. We compared the performance of different morphological measurements: 1) summary morphometrics calculated by L-Measure, 2) persistence vectors, a vectorized descriptor of branching structure, 3) the combination of the two. In all cases, we also investigated the impact of applying dimensionality reduction using principal component analysis (PCA). We assessed qualitative performance by gauging the ability to rank neurons in order of visual similarity. Moreover, we quantified information content by examining explained variance and benchmarked the ability to identify occasional duplicate reconstructions of the same specimen. We also compared two different methods for selecting the number of principal components using this benchmark. The results indicate that combining summary morphometrics and persistence vectors with applied PCA using maximum likelihood based automatic dimensionality selection provides an information rich characterization that enables efficient and precise comparison of neural morphology. We have deployed the similarity search as open-source online software both through a user-friendly graphical interface and as an API for programmatic access.

Smart imaging to empower brain-wide neuroscience at single-cell levels

A deep understanding of the neuronal connectivity and networks with detailed cell typing across brain regions is necessary to unravel the mechanisms behind the emotional and memorial functions as well as to find the treatment of brain impairment. Brain-wide imaging with single-cell resolution provides unique advantages to access morphological features of a neuron and to investigate the connectivity of neuron networks, which has led to exciting discoveries over the past years based on animal models, such as rodents. Nonetheless, high-throughput systems are in urgent demand to support studies of neural morphologies at larger scale and more detailed level, as well as to enable research on non-human primates (NHP) and human brains. The advances in artificial intelligence (AI) and computational resources bring great opportunity to 'smart' imaging systems, i.e., to automate, speed up, optimize and upgrade the imaging systems with AI and computational strategies. In this light, we review the important computational techniques that can support smart systems in brain-wide imaging at single-cell resolution.

The branching code: A model of actin-driven dendrite arborization

The cytoskeleton is crucial for defining neuronal-type-specific dendrite morphologies. To explore how the complex interplay of actin-modulatory proteins (AMPs) can define neuronal types in vivo, we focused on the class III dendritic arborization (c3da) neuron of Drosophila larvae. Using computational modeling, we reveal that the main branches (MBs) of c3da neurons follow general models based on optimal wiring principles, while the actin-enriched short terminal branches (STBs) require an additional growth program. To clarify the cellular mechanisms that define this second step, we thus concentrated on STBs for an in-depth quantitative description of dendrite morphology and dynamics. Applying these methods systematically to mutants of six known and novel AMPs, we revealed the complementary roles of these individual AMPs in defining STB properties. Our data suggest that diverse dendrite arbors result from a combination of optimal-wiring-related growth and individualized growth programs that are neuron-type specific.

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.

Quantitative neuronal morphometry by supervised and unsupervised learning

We present a protocol to characterize the morphological properties of individual neurons reconstructed from microscopic imaging. We first describe a simple procedure to extract relevant morphological features from digital tracings of neural arbors. Then, we provide detailed steps on classification, clustering, and statistical analysis of the traced cells based on morphological features. We illustrate the pipeline design using specific examples from zebrafish anatomy. Our approach can be readily applied and generalized to the characterization of axonal, dendritic, or glial geometry. For complete context and scientific motivation for the studies and datasets used here, refer to Valera et al. (2021).

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.

Organization of the inputs and outputs of the mouse superior colliculus

The superior colliculus (SC) receives diverse and robust cortical inputs to drive a range of cognitive and sensorimotor behaviors. However, it remains unclear how descending cortical input arising from higher-order associative areas coordinate with SC sensorimotor networks to influence its outputs. Here, we construct a comprehensive map of all cortico-tectal projections and identify four collicular zones with differential cortical inputs: medial (SC.m), centromedial (SC.cm), centrolateral (SC.cl) and lateral (SC.l). Further, we delineate the distinctive brain-wide input/output organization of each collicular zone, assemble multiple parallel cortico-tecto-thalamic subnetworks, and identify the somatotopic map in the SC that displays distinguishable spatial properties from the somatotopic maps in the neocortex and basal ganglia. Finally, we characterize interactions between those cortico-tecto-thalamic and cortico-basal ganglia-thalamic subnetworks. This study provides a structural basis for understanding how SC is involved in integrating different sensory modalities, translating sensory information to motor command, and coordinating different actions in goal-directed behaviors.

Scalar Field Comparison with Topological Descriptors: Properties and Applications for Scientific Visualization

In topological data analysis and visualization, topological descriptors such as persistence diagrams, merge trees, contour trees, Reeb graphs, and Morse–Smale complexes play an essential role in capturing the shape of scalar field data. We present a state‐of‐the‐art report on scalar field comparison using topological descriptors. We provide a taxonomy of existing approaches based on visualization tasks associated with three categories of data: single fields, time‐varying fields, and ensembles. These tasks include symmetry detection, periodicity detection, key event/feature detection, feature tracking, clustering, and structure statistics. Our main contributions include the formulation of a set of desirable mathematical and computational properties of comparative measures, and the classification of visualization tasks and applications that are enabled by these measures.

Connectivity characterization of the mouse basolateral amygdalar complex

The basolateral amygdalar complex (BLA) is implicated in behaviors ranging from fear acquisition to addiction. Optogenetic methods have enabled the association of circuit-specific functions to uniquely connected BLA cell types. Thus, a systematic and detailed connectivity profile of BLA projection neurons to inform granular, cell type-specific interrogations is warranted. Here, we apply machine-learning based computational and informatics analysis techniques to the results of circuit-tracing experiments to create a foundational, comprehensive BLA connectivity map. The analyses identify three distinct domains within the anterior BLA (BLAa) that house target-specific projection neurons with distinguishable morphological features. We identify brain-wide targets of projection neurons in the three BLAa domains, as well as in the posterior BLA, ventral BLA, posterior basomedial, and lateral amygdalar nuclei. Inputs to each nucleus also are identified via retrograde tracing. The data suggests that connectionally unique, domain-specific BLAa neurons are associated with distinct behavior networks.

A neuronal blueprint for directional mechanosensation in larval zebrafish

Animals have a remarkable ability to use local cues to orient in space in the absence of a panoramic fixed reference frame. Here we use the mechanosensory lateral line in larval zebrafish to understand rheotaxis, an innate oriented swimming evoked by water currents. We generated a comprehensive light-microscopy cell-resolution projectome of lateralis afferent neurons (LANs) and used clustering techniques for morphological classification. We find surprising structural constancy among LANs. Laser-mediated microlesions indicate that precise topographic mapping of lateral-line receptors is not essential for rheotaxis. Recording neuronal-activity during controlled mechanical stimulation of neuromasts reveals unequal representation of water-flow direction in the hindbrain. We explored potential circuit architectures constrained by anatomical and functional data to suggest a parsimonious model under which the integration of lateralized signals transmitted by direction-selective LANs underlies the encoding of water-flow direction in the brain. These data provide a new framework to understand how animals use local mechanical cues to orient in space.

A Systematic Evaluation of Interneuron Morphology Representations for Cell Type Discrimination

Quantitative analysis of neuronal morphologies usually begins with choosing a particular feature representation in order to make individual morphologies amenable to standard statistics tools and machine learning algorithms. Many different feature representations have been suggested in the literature, ranging from density maps to intersection profiles, but they have never been compared side by side. Here we performed a systematic comparison of various representations, measuring how well they were able to capture the difference between known morphological cell types. For our benchmarking effort, we used several curated data sets consisting of mouse retinal bipolar cells and cortical inhibitory neurons. We found that the best performing feature representations were two-dimensional density maps, two-dimensional persistence images and morphometric statistics, which continued to perform well even when neurons were only partially traced. Combining these feature representations together led to further performance increases suggesting that they captured non-redundant information. The same representations performed well in an unsupervised setting, implying that they can be suitable for dimensionality reduction or clustering.

Event history and topological data analysis

Persistent homology is used to track the appearance and disappearance of features as we move through a nested sequence of topological spaces. Equating the nested sequence to a filtration and the appearance and disappearance of features to events, we show that simple event history methods can be used for the analysis of topological data. We propose a version of the well-known Nelson–Aalen cumulative hazard estimator for the comparison of topological features of random fields and for testing parametric assumptions. We suggest a Cox proportional hazards approach for the analysis of embedded metric trees. The Nelson–Aalen method is illustrated on globally distributed climate data and on neutral hydrogen distribution in the Milky Way. The Cox method is used to compare vascular patterns in fundus images of the eyes of healthy and diabetic retinopathy patients.

Density Visualization Pipeline: A Tool for Cellular and Network Density Visualization and Analysis

Neuron classification is an important component in analyzing network structure and quantifying the effect of neuron topology on signal processing. Current quantification and classification approaches rely on morphology projection onto lower-dimensional spaces. In this paper a 3D visualization and quantification tool is presented. The Density Visualization Pipeline (DVP) computes, visualizes and quantifies the density distribution, i.e., the "mass" of interneurons. We use the DVP to characterize and classify a set of GABAergic interneurons. Classification of GABAergic interneurons is of crucial importance to understand on the one hand their various functions and on the other hand their ubiquitous appearance in the neocortex. 3D density map visualization and projection to the one-dimensional x, y, z subspaces show a clear distinction between the studied cells, based on these metrics. The DVP can be coupled to computational studies of the behavior of neurons and networks, in which network topology information is derived from DVP information. The DVP reads common neuromorphological file formats, e.g., Neurolucida XML files, NeuroMorpho.org SWC files and plain ASCII files. Full 3D visualization and projections of the density to 1D and 2D manifolds are supported by the DVP. All routines are embedded within the visual programming IDE VRL-Studio for Java which allows the definition and rapid modification of analysis workflows.

The natverse, a versatile toolbox for combining and analysing neuroanatomical data

To analyse neuron data at scale, neuroscientists expend substantial effort reading documentation, installing dependencies and moving between analysis and visualisation environments. To facilitate this, we have developed a suite of interoperable open-source R packages called the natverse. The natverse allows users to read local and remote data, perform popular analyses including visualisation and clustering and graph-theoretic analysis of neuronal branching. Unlike most tools, the natverse enables comparison across many neurons of morphology and connectivity after imaging or co-registration within a common template space. The natverse also enables transformations between different template spaces and imaging modalities. We demonstrate tools that integrate the vast majority of Drosophila neuroanatomical light microscopy and electron microscopy connectomic datasets. The natverse is an easy-to-use environment for neuroscientists to solve complex, large-scale analysis challenges as well as an open platform to create new code and packages to share with the community.

A Structural Average of Labeled Merge Trees for Uncertainty Visualization

Physical phenomena in science and engineering are frequently modeled using scalar fields. In scalar field topology, graph-based topological descriptors such as merge trees, contour trees, and Reeb graphs are commonly used to characterize topological changes in the (sub)level sets of scalar fields. One of the biggest challenges and opportunities to advance topology-based visualization is to understand and incorporate uncertainty into such topological descriptors to effectively reason about their underlying data. In this paper, we study a structural average of a set of labeled merge trees and use it to encode uncertainty in data. Specifically, we compute a 1-center tree that minimizes its maximum distance to any other tree in the set under a well-defined metric called the interleaving distance. We provide heuristic strategies that compute structural averages of merge trees whose labels do not fully agree. We further provide an interactive visualization system that resembles a numerical calculator that takes as input a set of merge trees and outputs a tree as their structural average. We also highlight structural similarities between the input and the average and incorporate uncertainty information for visual exploration. We develop a novel measure of uncertainty, referred to as consistency, via a metric-space view of the input trees. Finally, we demonstrate an application of our framework through merge trees that arise from ensembles of scalar fields. Our work is the first to employ interleaving distances and consistency to study a global, mathematically rigorous, structural average of merge trees in the context of uncertainty visualization.

Genetic Single Neuron Anatomy Reveals Fine Granularity of Cortical Axo-Axonic Cells

Parsing diverse nerve cells into biological types is necessary for understanding neural circuit organization. Morphology is an intuitive criterion for neuronal classification and a proxy of connectivity, but morphological diversity and variability often preclude resolving the granularity of neuron types. Combining genetic labeling with high-resolution, large-volume light microscopy, we established a single neuron anatomy platform that resolves, registers, and quantifies complete neuron morphologies in the mouse brain. We discovered that cortical axo-axonic cells (AACs), a cardinal GABAergic interneuron type that controls pyramidal neuron (PyN) spiking at axon initial segments, consist of multiple subtypes distinguished by highly laminar-specific soma position and dendritic and axonal arborization patterns. Whereas the laminar arrangements of AAC dendrites reflect differential recruitment by input streams, the laminar distribution and local geometry of AAC axons enable differential innervation of PyN ensembles. This platform will facilitate genetically targeted, high-resolution, and scalable single neuron anatomy in the mouse brain.