Shifting the paradigm: new approaches for characterizing and classifying neurons

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.

Classification of electrophysiological and morphological neuron types in the mouse visual cortex

Understanding the diversity of cell types in the brain has been an enduring challenge and requires detailed characterization of individual neurons in multiple dimensions. To systematically profile morpho-electric properties of mammalian neurons, we established a single-cell characterization pipeline using standardized patch-clamp recordings in brain slices and biocytin-based neuronal reconstructions. We built a publicly accessible online database, the Allen Cell Types Database, to display these datasets. Intrinsic physiological properties were measured from 1,938 neurons from the adult laboratory mouse visual cortex, morphological properties were measured from 461 reconstructed neurons, and 452 neurons had both measurements available. Quantitative features were used to classify neurons into distinct types using unsupervised methods. We established a taxonomy of morphologically and electrophysiologically defined cell types for this region of the cortex, with 17 electrophysiological types, 38 morphological types and 46 morpho-electric types. There was good correspondence with previously defined transcriptomic cell types and subclasses using the same transgenic mouse lines.

Single-nucleus and single-cell transcriptomes compared in matched cortical cell types

Transcriptomic profiling of complex tissues by single-nucleus RNA-sequencing (snRNA-seq) affords some advantages over single-cell RNA-sequencing (scRNA-seq). snRNA-seq provides less biased cellular coverage, does not appear to suffer cell isolation-based transcriptional artifacts, and can be applied to archived frozen specimens. We used well-matched snRNA-seq and scRNA-seq datasets from mouse visual cortex to compare cell type detection. Although more transcripts are detected in individual whole cells (~11,000 genes) than nuclei (~7,000 genes), we demonstrate that closely related neuronal cell types can be similarly discriminated with both methods if intronic sequences are included in snRNA-seq analysis. We estimate that the nuclear proportion of total cellular mRNA varies from 20% to over 50% for large and small pyramidal neurons, respectively. Together, these results illustrate the high information content of nuclear RNA for characterization of cellular diversity in brain tissues.

Comparing the Expression of Genes Related to Serotonin (5-HT) in C57BL/6J Mice and Humans Based on Data Available at the Allen Mouse Brain Atlas and Allen Human Brain Atlas

Brain atlases are tools based on comprehensive studies used to locate biological characteristics (structures, connections, proteins, and gene expression) in different regions of the brain. These atlases have been disseminated to the point where tools have been created to store, manage, and share the information they contain. This study used the data published by the Allen Mouse Brain Atlas (2004) for mice (C57BL/6J) and Allen Human Brain Atlas (2010) for humans (6 donors) to compare the expression of serotonin-related genes. Genes of interest were searched for manually in each case (in situ hybridization for mice and microarrays for humans), normalized expression data (z-scores) were extracted, and the results were graphed. Despite the differences in methodology, quantification, and subjects used in the process, a high degree of similarity was found between expression data. Here we compare expression in a way that allows the use of translational research methods to infer and validate knowledge. This type of study allows part of the relationship between structures and functions to be identified, by examining expression patterns and comparing levels of expression in different states, anatomical correlations, and phenotypes between different species. The study concludes by discussing the importance of knowing, managing, and disseminating comprehensive, open-access studies in neuroscience.

NPY+-, but not PV+- GABAergic neurons mediated long-range inhibition from infra- to prelimbic cortex

Anxiety disorders are thought to reflect deficits in the regulation of fear memories. While the amygdala has long been considered a site of storage of fear memories, newer findings suggest that the prefrontal cortex (PFC) is essential in the regulation of amygdala-dependent memories and fear expression. Here, activation of the prelimbic cortex (PrL) enhances the expression of fear, while an elevated activity in the infralimbic cortex (IL) enhances fear extinction. Despite the presence of these facts, we still know very little about the synaptic interconnectivity within the PFC. The aim of the present study was to investigate the inhibitory circuits between prelimbic and IL using morphological and electrophysiological methods. Our immunohistochemical analysis revealed that the distribution of PV(+)- and NPY(+)-GABAergic neurons was strikingly different within the PFC. In addition, we provided the first experimental evidence that the pyramidal neurons in the PrL received a direct inhibitory input mediated by bipolar NPY(+)-GABAergic projection neurons in the IL. Deletion of the anxiety-related neuroligin 2 gene caused a decrease of this direct synaptic inhibition that originated from the IL. Thus, our data suggested that activation of the IL might not only directly activate the corresponding downstream anxiolytic pathway, but also suppress the PrL-related anxiogenic pathway and thus could differentially bias the regulation of fear expression and extinction.

Classifying neuronal subclasses of the cerebellum through constellation pharmacology

A pressing need in neurobiology is the comprehensive identification and characterization of neuronal subclasses within the mammalian nervous system. To this end, we used constellation pharmacology as a method to interrogate the neuronal and glial subclasses of the mouse cerebellum individually and simultaneously. We then evaluated the data obtained from constellation-pharmacology experiments by cluster analysis to classify cells into neuronal and glial subclasses, based on their functional expression of glutamate, acetylcholine, and GABA receptors, among other ion channels. Conantokin peptides were used to identify N-methyl-d-aspartate (NMDA) receptor subtypes, which revealed that neurons of the young mouse cerebellum expressed NR2A and NR2B NMDA receptor subunits. Additional pharmacological tools disclosed differential expression of α-amino-3-hydroxy-5-methyl-4-isoxazloepropionic, nicotinic acetylcholine, and muscarinic acetylcholine receptors in different neuronal and glial subclasses. Certain cell subclasses correlated with known attributes of granule cells, and we combined constellation pharmacology with genetically labeled neurons to identify and characterize Purkinje cells. This study illustrates the utility of applying constellation pharmacology to classify neuronal and glial subclasses in specific anatomical regions of the brain.

Constellation pharmacology: a new paradigm for drug discovery

Constellation pharmacology is a cell-based high-content phenotypic-screening platform that utilizes subtype-selective pharmacological agents to elucidate the cell-specific combinations (constellations) of key signaling proteins that define specific cell types. Heterogeneous populations of native cells, in which the different individual cell types have been identified and characterized, are the foundation for this screening platform. Constellation pharmacology is useful for screening small molecules or for deconvoluting complex mixtures of biologically active natural products. This platform has been used to purify natural products and discover their molecular mechanisms. In the ongoing development of constellation pharmacology, there is a positive feedback loop between the pharmacological characterization of cell types and screening for new drug candidates. As constellation pharmacology is used to discover compounds with novel targeting-selectivity profiles, those new compounds then further help to elucidate the constellations of specific cell types, thereby increasing the content of this high-content platform.

Towards the automatic classification of neurons

The classification of neurons into types has been much debated since the inception of modern neuroscience. Recent experimental advances are accelerating the pace of data collection. The resulting growth of information about morphological, physiological, and molecular properties encourages efforts to automate neuronal classification by powerful machine learning techniques. We review state-of-the-art analysis approaches and the availability of suitable data and resources, highlighting prominent challenges and opportunities. The effective solution of the neuronal classification problem will require continuous development of computational methods, high-throughput data production, and systematic metadata organization to enable cross-laboratory integration.

Defining modulatory inputs into CNS neuronal subclasses by functional pharmacological profiling

Previously we defined neuronal subclasses within the mouse peripheral nervous system using an experimental strategy called "constellation pharmacology." Here we demonstrate the broad applicability of constellation pharmacology by extending it to the CNS and specifically to the ventral respiratory column (VRC) of mouse brainstem, a region containing the neuronal network controlling respiratory rhythm. Analysis of dissociated cells from this locus revealed three major cell classes, each encompassing multiple subclasses. We broadly analyzed the combinations (constellations) of receptors and ion channels expressed within VRC cell classes and subclasses. These were strikingly different from the constellations of receptors and ion channels found in subclasses of peripheral neurons from mouse dorsal root ganglia. Within the VRC cell population, a subset of dissociated neurons responded to substance P, putatively corresponding to inspiratory pre-Bötzinger complex (preBötC) neurons. Using constellation pharmacology, we found that these substance P-responsive neurons also responded to histamine, and about half responded to bradykinin. Electrophysiological studies conducted in brainstem slices confirmed that preBötC neurons responsive to substance P exhibited similar responsiveness to bradykinin and histamine. The results demonstrate the predictive utility of constellation pharmacology for defining modulatory inputs into specific neuronal subclasses within central neuronal networks.

Using constellation pharmacology to define comprehensively a somatosensory neuronal subclass

Change is intrinsic to nervous systems; change is required for learning and conditioning and occurs with disease progression, normal development, and aging. To better understand mammalian nervous systems and effectively treat nervous-system disorders, it is essential to track changes in relevant individual neurons. A critical challenge is to identify and characterize the specific cell types involved and the molecular-level changes that occur in each. Using an experimental strategy called constellation pharmacology, we demonstrate that we can define a specific somatosensory neuronal subclass, cold thermosensors, across different species and track changes in these neurons as a function of development. Cold thermosensors are uniformly responsive to menthol and innocuous cool temperature (17 °C), indicating that they express TRPM8 channels. A subset of cold thermosensors expressed α7 nicotinic acetylcholine receptors (nAChRs) but not other nAChR subtypes. Differences in temperature threshold of cold thermosensors correlated with functional expression of voltage-gated K channels Kv1.1/1.2: Relatively higher expression of KV1.1/1.2 channels resulted in a higher threshold response to cold temperature. Other signaling components varied during development and between species. In cold thermosensors of neonatal mice and rats, ATP receptors were functionally expressed, but the expression disappeared with development. This developmental change occurred earlier in low-threshold than high-threshold cold thermosensors. Most rat cold thermosensors expressed TRPA1 channels, whereas mouse cold thermosensors did not. The broad implications of this study are that it is now feasible to track changes in receptor and ion-channel expression in individual neuronal subclasses as a function of development, learning, disease, or aging.

The neuron identity problem: form meets function

A complete understanding of nervous system function cannot be achieved without the identification of its component cell types. In this Perspective, we explore a series of related issues surrounding cell identity and how revolutionary methods for labeling and probing specific neuronal types have clarified this question. Specifically, we ask the following questions: what is the purpose of such diversity, how is it generated, how is it maintained, and, ultimately, how can one unambiguously identity one cell type from another? We suggest that each cell type can be defined by a unique and conserved molecular ground state that determines its capabilities. We believe that gaining an understanding of these molecular barcodes will advance our ability to explore brain function, enhance our understanding of the biochemical basis of CNS disorders, and aid in the development of novel therapeutic strategies.

Morphological and physiological characterization of pyramidal neuron subtypes in rat medial prefrontal cortex

The medial prefrontal cortex (mPFC) has been implicated in cognitive and executive processes including decision making, working memory and behavioral flexibility. Cortical processing depends on the interaction between distinct neuronal cell types in different cortical layers. To better understand cortical processing in the rat mPFC, we studied the diversity of pyramidal neurons using in vitro whole-cell patch clamp recordings and biocytin staining of neurons, followed by morphological analysis. Using unsupervised cluster analysis for the objective grouping of neurons, we identified more than 10 different pyramidal subtypes spread across the different cortical layers. Layer 2 pyramidal neurons possessed a unique morphology with wide apical dendritic field spans and a narrow basal field span. Layer 3 contained the only subtype that showed a burst of action potentials upon current injection. Layer 5 pyramidal neurons showed the largest voltage sags. Finally, pyramidal neurons in layer 6 (L6) showed a great variety in their morphology with 39% of L6 neurons possessing tall apical dendrites that extend into layer 1. Future experiments on the functional role of the mPFC should take into account the great diversity of pyramidal neurons.

Selective knockout of the casein kinase 2 in d1 medium spiny neurons controls dopaminergic function

BACKGROUND

Dopamine, crucial for the regulation of motor function and reward, acts through receptors mainly expressed in striatum as well as cortex. Dysregulation of dopaminergic signaling is associated with various neuropsychiatric disorders. Consequently, dopamine-regulating drugs are effectively used in treating these disorders, such as L-DOPA for Parkinson's disease, methylphenidate for attention-deficit/hyperactivity disorder, or antipsychotics for schizophrenia. As a result, there has been much interest in dissecting signaling networks in the two morphologically indistinguishable D1- and D2-receptor-expressing medium spiny neurons. Our previous results highlighted a role for casein kinase 2 (CK2) in the modulation of dopamine D1 receptor (D1R) signaling in cells.

METHODS

To study the importance of CK2 in vivo, we have selectively knocked out CK2, in either D1- or D2-medium spiny neurons (MSNs) and characterized the mice behaviorally and biochemically (n = 4-18).

RESULTS

The D1-MSN knockout mice exhibited distinct behavioral phenotypes including novelty-induced hyperlocomotion and exploratory behavior, defective motor control, and motor learning. All of these behavioral traits are indicative of dysregulated dopamine signaling and the underlying mechanism appears to be an alteration of D1R signaling. In support of this hypothesis, D1R levels were upregulated in the knockout mice, as well as phosphorylation of DARPP-32 (dopamine- and cyclic adenosine monophosphate [cAMP]-regulated phospho-protein of 32 kDa), most of the behavioral phenotypes were abolished by the D1R antagonist, SCH23390, and the D2-MSN knockout mice displayed no obvious behavioral phenotype.

CONCLUSIONS

A single kinase, CK2, in D1-MSNs significantly alters dopamine signaling, a finding that could have therapeutic implications for disorders characterized by dopamine imbalance such as Parkinson's disease, attention-deficit/hyperactivity disorder, and schizophrenia.

Computational genetic neuroanatomy of the developing mouse brain: dimensionality reduction, visualization, and clustering

Background

The structured organization of cells in the brain plays a key role in its functional efficiency. This delicate organization is the consequence of unique molecular identity of each cell gradually established by precise spatiotemporal gene expression control during development. Currently, studies on the molecular-structural association are beginning to reveal how the spatiotemporal gene expression patterns are related to cellular differentiation and structural development.

Results

In this article, we aim at a global, data-driven study of the relationship between gene expressions and neuroanatomy in the developing mouse brain. To enable visual explorations of the high-dimensional data, we map the in situ hybridization gene expression data to a two-dimensional space by preserving both the global and the local structures. Our results show that the developing brain anatomy is largely preserved in the reduced gene expression space. To provide a quantitative analysis, we cluster the reduced data into groups and measure the consistency with neuroanatomy at multiple levels. Our results show that the clusters in the low-dimensional space are more consistent with neuroanatomy than those in the original space.

Conclusions

Gene expression patterns and developing brain anatomy are closely related. Dimensionality reduction and visual exploration facilitate the study of this relationship.

A supplementary circuit rule-set for the neuronal wiring

Limitations of known anatomical circuit rules necessitate the identification of supplementary rules. This is essential for explaining how associative sensory stimuli induce nervous system changes that generate internal sensations of memory, concurrent with triggering specific motor activities in response to specific cue stimuli. A candidate mechanism is rapidly reversible, yet stabilizable membrane hemi-fusion formed between the closely apposed postsynaptic membranes of different neurons at locations of convergence of sensory inputs during associative learning. The lateral entry of activity from the cue stimulus-activated postsynapse re-activates the opposite postsynapse through the hemi-fused area and induces the basic units of internal sensation (namely, semblions) as a systems property. Working, short-term and long-term memories can be viewed as functions of the number of re-activatible hemi-fusions present at the time of memory retrieval. Blocking membrane hemi-fusion either by the insertion of the herpes simplex virus (HSV) glycoproteins or by the deposition of insoluble intermediates of amyloid protein in the inter-postsynaptic extracellular matrix (ECM) space leads to cognitive impairments, supporting this mechanism. The introduction of membrane fusion blockers into the postsynaptic cell cytoplasm that attenuates long-term potentiation (LTP), a correlate of behavioral motor activities in response to memory retrieval, provides further support. The lateral spread of activity through the inter-postsynaptic membrane is capable of contributing to oscillating neuronal activity at certain neuronal orders. At the resting state these oscillations provide sub-threshold activation to many neurons at higher orders, including motor neurons maintaining them at a low initiation threshold for motor activity.

Development and plasticity of the primary visual cortex

Hubel and Wiesel began the modern study of development and plasticity of primary visual cortex (V1), discovering response properties of cortical neurons that distinguished them from their inputs and that were arranged in a functional architecture. Their findings revealed an early innate period of development and a later critical period of dramatic experience-dependent plasticity. Recent studies have used rodents to benefit from biochemistry and genetics. The roles of spontaneous neural activity and molecular signaling in innate, experience-independent development have been clarified, as have the later roles of visual experience. Plasticity produced by monocular visual deprivation (MD) has been dissected into stages governed by distinct signaling mechanisms, some of whose molecular players are known. Many crucial questions remain, but new tools for perturbing cortical cells and measuring plasticity at the level of changes in connections among identified neurons now exist. The future for the study of V1 to illuminate cortical development and plasticity is bright.

An ontological approach to describing neurons and their relationships

The advancement of neuroscience, perhaps one of the most information rich disciplines of all the life sciences, requires basic frameworks for organizing the vast amounts of data generated by the research community to promote novel insights and integrated understanding. Since Cajal, the neuron remains a fundamental unit of the nervous system, yet even with the explosion of information technology, we still have few comprehensive or systematic strategies for aggregating cell-level knowledge. Progress toward this goal is hampered by the multiplicity of names for cells and by lack of a consensus on the criteria for defining neuron types. However, through umbrella projects like the Neuroscience Information Framework (NIF) and the International Neuroinformatics Coordinating Facility (INCF), we have the opportunity to propose and implement an informatics infrastructure for establishing common tools and approaches to describe neurons through a standard terminology for nerve cells and a database (a Neuron Registry) where these descriptions can be deposited and compared. This article provides an overview of the problem and outlines a solution approach utilizing ontological characterizations. Based on illustrative implementation examples, we also discuss the need for consensus criteria to be adopted by the research community, and considerations on future developments. A scalable repository of neuron types will provide researchers with a resource that materially contributes to the advancement of neuroscience.

Functional profiling of neurons through cellular neuropharmacology

We describe a functional profiling strategy to identify and characterize subtypes of neurons present in a peripheral ganglion, which should be extendable to neurons in the CNS. In this study, dissociated dorsal-root ganglion neurons from mice were exposed to various pharmacological agents (challenge compounds), while at the same time the individual responses of >100 neurons were simultaneously monitored by calcium imaging. Each challenge compound elicited responses in only a subset of dorsal-root ganglion neurons. Two general types of challenge compounds were used: agonists of receptors (ionotropic and metabotropic) that alter cytoplasmic calcium concentration (receptor-agonist challenges) and compounds that affect voltage-gated ion channels (membrane-potential challenges). Notably, among the latter are K-channel antagonists, which elicited unexpectedly diverse types of calcium responses in different cells (i.e., phenotypes). We used various challenge compounds to identify several putative neuronal subtypes on the basis of their shared and/or divergent functional, phenotypic profiles. Our results indicate that multiple receptor-agonist and membrane-potential challenges may be applied to a neuronal population to identify, characterize, and discriminate among neuronal subtypes. This experimental approach can uncover constellations of plasma membrane macromolecules that are functionally coupled to confer a specific phenotypic profile on each neuronal subtype. This experimental platform has the potential to bridge a gap between systems and molecular neuroscience with a cellular-focused neuropharmacology, ultimately leading to the identification and functional characterization of all neuronal subtypes at a given locus in the nervous system.

Post hoc immunostaining of GABAergic neuronal subtypes following in vivo two-photon calcium imaging in mouse neocortex

GABAergic neurons in the neocortex are diverse with regard to morphology, physiology, and axonal targeting pattern, indicating functional specializations within the cortical microcircuitry. Little information is available, however, about functional properties of distinct subtypes of GABAergic neurons in the intact brain. Here, we combined in vivo two-photon calcium imaging in supragranular layers of the mouse neocortex with post hoc immunohistochemistry against the three calcium-binding proteins parvalbumin, calretinin, and calbindin in order to assign subtype marker profiles to neuronal activity. Following coronal sectioning of fixed brains, we matched cells in corresponding volumes of image stacks acquired in vivo and in fixed brain slices. In GAD67-GFP mice, more than 95% of the GABAergic cells could be unambiguously matched, even in large volumes comprising more than a thousand interneurons. Triple immunostaining revealed a depth-dependent distribution of interneuron subtypes with increasing abundance of PV-positive neurons with depth. Most importantly, the triple-labeling approach was compatible with previous in vivo calcium imaging following bulk loading of Oregon Green 488 BAPTA-1, which allowed us to classify spontaneous calcium transients recorded in vivo according to the neurochemically defined GABAergic subtypes. Moreover, we demonstrate that post hoc immunostaining can also be applied to wild-type mice expressing the genetically encoded calcium indicator Yellow Cameleon 3.60 in cortical neurons. Our approach is a general and flexible method to distinguish GABAergic subtypes in cell populations previously imaged in the living animal. It should thus facilitate dissecting the functional roles of these subtypes in neural circuitry.

A guide to delineate the logic of neurovascular signaling in the brain

The neurovascular system may be viewed as a distributed nervous system within the brain. It transforms local neuronal activity into a change in the tone of smooth muscle that lines the walls of arterioles and microvessels. We review the current state of neurovascular coupling, with an emphasis on signaling molecules that convey information from neurons to neighboring vessels. At the level of neocortex, this coupling is mediated by: (i) a likely direct interaction with inhibitory neurons, (ii) indirect interaction, via astrocytes, with excitatory neurons, and (iii) fiber tracts from subcortical layers. Substantial evidence shows that control involves competition between signals that promote vasoconstriction versus vasodilation. Consistent with this picture is evidence that, under certain circumstances, increased neuronal activity can lead to vasoconstriction rather than vasodilation. This confounds naïve interpretations of functional brain images. We discuss experimental approaches to detect signaling molecules in vivo with the goal of formulating an empirical basis for the observed logic of neurovascular control.

Bridging the gaps between synapses, circuits, and behavior

The decade of the brain may have come and gone, but the final frontier, cracking the neuronal code, still lies ahead. Today, new technologies that allow precise spatiotemporal remote control over the activity of genetically defined populations of neurons within intact neural circuits are providing a means of obtaining a functional wiring diagram of the mammalian brain, bringing us one step closer to understanding precisely how neuronal activity codes for perception, thought, emotion, and action. These technologies and the design principles underlying them are reviewed herein.