Connecting single-cell transcriptomes to projectomes in mouse visual cortex.

The mammalian brain is composed of diverse neuron types that play different functional roles. Recent single-cell RNA sequencing approaches have led to a whole brain taxonomy of transcriptomically-defined cell types, yet cell type definitions that include multiple cellular properties can offer additional insights into a neuron's role in brain circuits. While the Patch-seq method can investigate how transcriptomic properties relate to the local morphological and electrophysiological properties of cell types, linking transcriptomic identities to long-range projections is a major unresolved challenge. To address this, we collected coordinated Patch-seq and whole brain morphology data sets of excitatory neurons in mouse visual cortex. From the Patch-seq data, we defined 16 integrated morpho-electric-transcriptomic (MET)-types; in parallel, we reconstructed the complete morphologies of 300 neurons. We unified the two data sets with a multi-step classifier, to integrate cell type assignments and interrogate cross-modality relationships. We find that transcriptomic variations within and across MET-types correspond with morphological and electrophysiological phenotypes. In addition, this variation, along with the anatomical location of the cell, can be used to predict the projection targets of individual neurons. We also shed new light on infragranular cell types and circuits, including cell-type-specific, interhemispheric projections. With this approach, we establish a comprehensive, integrated taxonomy of excitatory neuron types in mouse visual cortex and create a system for integrated, high-dimensional cell type classification that can be extended to the whole brain and potentially across species.

Human neocortical expansion involves glutamatergic neuron diversification.

The neocortex is disproportionately expanded in human compared with mouse1,2, both in its total volume relative to subcortical structures and in the proportion occupied by supragranular layers composed of neurons that selectively make connections within the neocortex and with other telencephalic structures. Single-cell transcriptomic analyses of human and mouse neocortex show an increased diversity of glutamatergic neuron types in supragranular layers in human neocortex and pronounced gradients as a function of cortical depth3. Here, to probe the functional and anatomical correlates of this transcriptomic diversity, we developed a robust platform combining patch clamp recording, biocytin staining and single-cell RNA-sequencing (Patch-seq) to examine neurosurgically resected human tissues. We demonstrate a strong correspondence between morphological, physiological and transcriptomic phenotypes of five human glutamatergic supragranular neuron types. These were enriched in but not restricted to layers, with one type varying continuously in all phenotypes across layers 2 and 3. The deep portion of layer 3 contained highly distinctive cell types, two of which express a neurofilament protein that labels long-range projection neurons in primates that are selectively depleted in Alzheimer's disease4,5. Together, these results demonstrate the explanatory power of transcriptomic cell-type classification, provide a structural underpinning for increased complexity of cortical function in humans, and implicate discrete transcriptomic neuron types as selectively vulnerable in disease.

Scaled, high fidelity electrophysiological, morphological, and transcriptomic cell characterization.

The Patch-seq approach is a powerful variation of the patch-clamp technique that allows for the combined electrophysiological, morphological, and transcriptomic characterization of individual neurons. To generate Patch-seq datasets at scale, we identified and refined key factors that contribute to the efficient collection of high-quality data. We developed patch-clamp electrophysiology software with analysis functions specifically designed to automate acquisition with online quality control. We recognized the importance of extracting the nucleus for transcriptomic success and maximizing membrane integrity during nucleus extraction for morphology success. The protocol is generalizable to different species and brain regions, as demonstrated by capturing multimodal data from human and macaque brain slices. The protocol, analysis and acquisition software are compiled at https://githubcom/AllenInstitute/patchseqtools. This resource can be used by individual labs to generate data across diverse mammalian species and that is compatible with large publicly available Patch-seq datasets.

A taxonomy of transcriptomic cell types across the isocortex and hippocampal formation.

The isocortex and hippocampal formation (HPF) in the mammalian brain play critical roles in perception, cognition, emotion, and learning. We profiled ∼1.3 million cells covering the entire adult mouse isocortex and HPF and derived a transcriptomic cell-type taxonomy revealing a comprehensive repertoire of glutamatergic and GABAergic neuron types. Contrary to the traditional view of HPF as having a simpler cellular organization, we discover a complete set of glutamatergic types in HPF homologous to all major subclasses found in the six-layered isocortex, suggesting that HPF and the isocortex share a common circuit organization. We also identify large-scale continuous and graded variations of cell types along isocortical depth, across the isocortical sheet, and in multiple dimensions in hippocampus and subiculum. Overall, our study establishes a molecular architecture of the mammalian isocortex and hippocampal formation and begins to shed light on its underlying relationship with the development, evolution, connectivity, and function of these two brain structures.

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.

Integrated Morphoelectric and Transcriptomic Classification of Cortical GABAergic Cells.

Neurons are frequently classified into distinct types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 4,200 mouse visual cortical GABAergic interneurons and reconstructed the local morphologies of 517 of those neurons. We find that most transcriptomic types (t-types) occupy specific laminar positions within visual cortex, and, for most types, the cells mapping to a t-type exhibit consistent electrophysiological and morphological properties. These properties display both discrete and continuous variation among t-types. Through multimodal integrated analysis, we define 28 met-types that have congruent morphological, electrophysiological, and transcriptomic properties and robust mutual predictability. We identify layer-specific axon innervation pattern as a defining feature distinguishing different met-types. These met-types represent a unified definition of cortical GABAergic interneuron types, providing a systematic framework to capture existing knowledge and bridge future analyses across different modalities.

Systematic Integration of Structural and Functional Data into Multi-scale Models of Mouse Primary Visual Cortex.

Structural rules underlying functional properties of cortical circuits are poorly understood. To explore these rules systematically, we integrated information from extensive literature curation and large-scale experimental surveys into a data-driven, biologically realistic simulation of the awake mouse primary visual cortex. The model was constructed at two levels of granularity, using either biophysically detailed or point neurons. Both variants have identical network connectivity and were compared to each other and to experimental recordings of visual-driven neural activity. While tuning these networks to recapitulate experimental data, we identified rules governing cell-class-specific connectivity and synaptic strengths. These structural constraints constitute hypotheses that can be tested experimentally. Despite their distinct single-cell abstraction, both spatially extended and point models perform similarly at the level of firing rate distributions for the questions we investigated. All data and models are freely available as a resource for the community.

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.

Visual physiology of the layer 4 cortical circuit in silico.

Despite advances in experimental techniques and accumulation of large datasets concerning the composition and properties of the cortex, quantitative modeling of cortical circuits under in-vivo-like conditions remains challenging. Here we report and publicly release a biophysically detailed circuit model of layer 4 in the mouse primary visual cortex, receiving thalamo-cortical visual inputs. The 45,000-neuron model was subjected to a battery of visual stimuli, and results were compared to published work and new in vivo experiments. Simulations reproduced a variety of observations, including effects of optogenetic perturbations. Critical to the agreement between responses in silico and in vivo were the rules of functional synaptic connectivity between neurons. Interestingly, after extreme simplification the model still performed satisfactorily on many measurements, although quantitative agreement with experiments suffered. These results emphasize the importance of functional rules of cortical wiring and enable a next generation of data-driven models of in vivo neural activity and computations.

BioNet: A Python interface to NEURON for modeling large-scale networks.

There is a significant interest in the neuroscience community in the development of large-scale network models that would integrate diverse sets of experimental data to help elucidate mechanisms underlying neuronal activity and computations. Although powerful numerical simulators (e.g., NEURON, NEST) exist, data-driven large-scale modeling remains challenging due to difficulties involved in setting up and running network simulations. We developed a high-level application programming interface (API) in Python that facilitates building large-scale biophysically detailed networks and simulating them with NEURON on parallel computer architecture. This tool, termed "BioNet", is designed to support a modular workflow whereby the description of a constructed model is saved as files that could be subsequently loaded for further refinement and/or simulation. The API supports both NEURON's built-in as well as user-defined models of cells and synapses. It is capable of simulating a variety of observables directly supported by NEURON (e.g., spikes, membrane voltage, intracellular [Ca++]), as well as plugging in modules for computing additional observables (e.g. extracellular potential). The high-level API platform obviates the time-consuming development of custom code for implementing individual models, and enables easy model sharing via standardized files. This tool will help refocus neuroscientists on addressing outstanding scientific questions rather than developing narrow-purpose modeling code.

Use of Radiation Therapy Within the Last Year of Life Among Cancer Patients.

PURPOSE: We examined radiation therapy (RT) use within the last year of life (LYOL). As palliative RT (PRT) has been well studied in patients with ≥6-week life expectancies, we hypothesized that PRT use would be constant over the LYOL, except for the last 30 days, when use would decline given lack of prospective data supporting it. MATERIALS AND METHODS: At a single institution, 870 cancer patients died between October 2, 2014, and September 30, 2015, and had ≥3 evaluation and management visits within the LYOL. Claims and RT data were extracted and linked. Over the LYOL, we evaluated RT use by intent (curative vs palliative) and indications. RESULTS: Within the LYOL, one-third of patients underwent RT in the last 365 days of life to 444 sites, which decreased to 24.3% and 8.5% in the last 180 and 30 days of life, respectively. Patients who received any RT in the last 365 days of life were younger at death and had a higher proportion of lung, sarcoma, and transplant disease groups. One-quarter of sites were irradiated with curative intent, which remained constant over the LYOL. In contrast, PRT was used at a supralinear rate, in which treatment of bone metastases and use of single-fraction PRT increased closer to death. CONCLUSIONS: PRT appears to be disproportionately used closer to death, with an increasing proportion of irradiated sites being bone metastases. This may be secondary to increased symptoms from advanced cancer toward the end of life. As patients with very poor prognoses (eg, within 30 days of death) are generally not included in RT clinical trials, further studies are warranted to assess whether PRT for bone metastases at the end of life is efficacious.

Systematic generation of biophysically detailed models for diverse cortical neuron types.

The cellular components of mammalian neocortical circuits are diverse, and capturing this diversity in computational models is challenging. Here we report an approach for generating biophysically detailed models of 170 individual neurons in the Allen Cell Types Database to link the systematic experimental characterization of cell types to the construction of cortical models. We build models from 3D morphologies and somatic electrophysiological responses measured in the same cells. Densities of active somatic conductances and additional parameters are optimized with a genetic algorithm to match electrophysiological features. We evaluate the models by applying additional stimuli and comparing model responses to experimental data. Applying this technique across a diverse set of neurons from adult mouse primary visual cortex, we verify that models preserve the distinctiveness of intrinsic properties between subsets of cells observed in experiments. The optimized models are accessible online alongside the experimental data. Code for optimization and simulation is also openly distributed.

Inferring cortical function in the mouse visual system through large-scale systems neuroscience.

The scientific mission of the Project MindScope is to understand neocortex, the part of the mammalian brain that gives rise to perception, memory, intelligence, and consciousness. We seek to quantitatively evaluate the hypothesis that neocortex is a relatively homogeneous tissue, with smaller functional modules that perform a common computational function replicated across regions. We here focus on the mouse as a mammalian model organism with genetics, physiology, and behavior that can be readily studied and manipulated in the laboratory. We seek to describe the operation of cortical circuitry at the computational level by comprehensively cataloging and characterizing its cellular building blocks along with their dynamics and their cell type-specific connectivities. The project is also building large-scale experimental platforms (i.e., brain observatories) to record the activity of large populations of cortical neurons in behaving mice subject to visual stimuli. A primary goal is to understand the series of operations from visual input in the retina to behavior by observing and modeling the physical transformations of signals in the corticothalamic system. We here focus on the contribution that computer modeling and theory make to this long-term effort.

Synchronization of firing in cortical fast-spiking interneurons at gamma frequencies: a phase-resetting analysis.

Fast-spiking (FS) cells in the neocortex are interconnected both by inhibitory chemical synapses and by electrical synapses, or gap-junctions. Synchronized firing of FS neurons is important in the generation of gamma oscillations, at frequencies between 30 and 80 Hz. To understand how these synaptic interactions control synchronization, artificial synaptic conductances were injected in FS cells, and the synaptic phase-resetting function (SPRF), describing how the compound synaptic input perturbs the phase of gamma-frequency spiking as a function of the phase at which it is applied, was measured. GABAergic and gap junctional conductances made distinct contributions to the SPRF, which had a surprisingly simple piecewise linear form, with a sharp midcycle break between phase delay and advance. Analysis of the SPRF showed how the intrinsic biophysical properties of FS neurons and their interconnections allow entrainment of firing over a wide gamma frequency band, whose upper and lower frequency limits are controlled by electrical synapses and GABAergic inhibition respectively.

Signal propagation in Drosophila central neurons.

Drosophila is an important model organism for investigating neural development, neural morphology, neurophysiology, and neural correlates of behaviors. However, almost nothing is known about how electrical signals propagate in Drosophila neurons. Here, we address these issues in antennal lobe projection neurons, one of the most well studied classes of Drosophila neurons. We use morphological and electrophysiological data to deduce the passive membrane properties of these neurons and to build a compartmental model of their electrotonic structure. We find that these neurons are electrotonically extensive and that a somatic recording electrode can only imperfectly control the voltage in the rest of the cell. Simulations predict that action potentials initiate at a location distant from the soma, in the proximal portion of the axon. Simulated synaptic input to a single dendritic branch propagates poorly to the rest of the cell and cannot match the size of real unitary synaptic events, but we can obtain a good fit to data when we model unitary input synapses as dozens of release sites distributed across many dendritic branches. We also show that the true resting potential of these neurons is more hyperpolarized than previously thought, attributable to the experimental error introduced by the electrode seal conductance. A leak sodium conductance also contributes to the resting potential. Together, these findings have fundamental implications for how these neurons integrate their synaptic inputs. Our results also have important consequences for the design and interpretation of experiments aimed at understanding Drosophila neurons and neural circuits.

Sensory processing in the Drosophila antennal lobe increases reliability and separability of ensemble odor representations.

Here we describe several fundamental principles of olfactory processing in the Drosophila melanogaster antennal lobe (the analog of the vertebrate olfactory bulb), through the systematic analysis of input and output spike trains of seven identified glomeruli. Repeated presentations of the same odor elicit more reproducible responses in second-order projection neurons (PNs) than in their presynaptic olfactory receptor neurons (ORNs). PN responses rise and accommodate rapidly, emphasizing odor onset. Furthermore, weak ORN inputs are amplified in the PN layer but strong inputs are not. This nonlinear transformation broadens PN tuning and produces more uniform distances between odor representations in PN coding space. In addition, portions of the odor response profile of a PN are not systematically related to their direct ORN inputs, which probably indicates the presence of lateral connections between glomeruli. Finally, we show that a linear discriminator classifies odors more accurately using PN spike trains than using an equivalent number of ORN spike trains.

The contribution of resurgent sodium current to high-frequency firing in Purkinje neurons: an experimental and modeling study.

Purkinje neurons generate high-frequency action potentials and express voltage-gated, tetrodotoxin-sensitive sodium channels with distinctive kinetics. Their sodium currents activate and inactivate during depolarization, as well as reactivate during repolarization from positive potentials, producing a "resurgent" current. This reopening of channels not only generates inward current after each action potential, but also permits rapid recovery from inactivation, leading to the hypothesis that resurgent current may facilitate high-frequency firing. Mutant med mice are ataxic and lack expression of the Scn8a gene, which encodes the NaV1.6 protein. In med Purkinje cells, transient sodium current inactivates more rapidly than in wild-type cells, and resurgent current is nearly abolished. To investigate how NaV1.6-specific kinetics influence firing patterns, we recorded action potentials of Purkinje neurons isolated from wild-type and med mice. We also recorded non-sodium currents from Purkinje cells of both genotypes to test whether the Scn8a mutation induced changes in other ion channels. Last, we modeled action potential firing by simulating eight currents directly recorded from Purkinje cells in both wild-type and med mice. Regular, high-frequency firing was slowed in med Purkinje neurons. In addition to disrupted sodium currents, med neurons had small but significant changes in potassium and leak currents. Simulations indicated that these modified non-sodium currents could not account for the reduced excitability of med cells but instead slightly facilitated spiking. The loss of NaV1.6-specific kinetics, however, slowed simulated spontaneous activity. Together, the data suggest that across a range of conditions, sodium currents with a resurgent component promote and accelerate firing.