Epizootic pneumonia of bighorn sheep is a devastating disease of uncertain etiology. To help clarify the etiology, we used culture and culture-independent methods to compare the prevalence of the bacterial respiratory pathogens Mannheimia haemolytica, Bibersteinia trehalosi, Pasteurella multocida, and Mycoplasma ovipneumoniae in lung tissue from 44 bighorn sheep from herds affected by 8 outbreaks in the western United States. M. ovipneumoniae, the only agent detected at significantly higher prevalence in animals from outbreaks (95%) than in animals from unaffected healthy populations (0%), was the most consistently detected agent and the only agent that exhibited single strain types within each outbreak. The other respiratory pathogens were frequently but inconsistently detected, as were several obligate anaerobic bacterial species, all of which might represent secondary or opportunistic infections that could contribute to disease severity. These data provide evidence that M. ovipneumoniae plays a primary role in the etiology of epizootic pneumonia of bighorn sheep.

The neocortex contains a multitude of cell types that are segregated into layers and functionally distinct areas. To investigate the diversity of cell types across the mouse neocortex, here we analysed 23,822 cells from two areas at distant poles of the mouse neocortex: the primary visual cortex and the anterior lateral motor cortex. We define 133 transcriptomic cell types by deep, single-cell RNA sequencing. Nearly all types of GABA (γ-aminobutyric acid)-containing neurons are shared across both areas, whereas most types of glutamatergic neurons were found in one of the two areas. By combining single-cell RNA sequencing and retrograde labelling, we match transcriptomic types of glutamatergic neurons to their long-range projection specificity. Our study establishes a combined transcriptomic and projectional taxonomy of cortical cell types from functionally distinct areas of the adult mouse cortex.

Exposure to trauma is a risk factor for the development of a number of mood disorders, and may enhance vulnerability to future adverse life events. Recent data demonstrate that ventral tegmental area (VTA) neurons expressing the vesicular glutamate transporter 2 (VGluT2) signal and causally contribute to behaviors that involve aversive or threatening stimuli. However, it is unknown whether VTA VGluT2 neurons regulate transsituational outcomes of stress and whether these neurons are sensitive to stressor controllability. This work adapted an operant mouse paradigm to examine the impact of stressor controllability on VTA VGluT2 neuron function as well as the role of VTA VGluT2 neurons in mediating transsituational stressor outcomes. Uncontrollable (inescapable) stress, but not physically identical controllable (escapable) stress, produced social avoidance and exaggerated fear in male mice. Uncontrollable stress in females led to exploratory avoidance of a novel brightly lit environment. Both controllable and uncontrollable stressors increased VTA VGluT2 neuronal activity, and chemogenetic silencing of VTA VGluT2 neurons prevented the behavioral sequelae of uncontrollable stress in male and female mice. Further, we show that stress activates multiple genetically-distinct subtypes of VTA VGluT2 neurons, especially those that are VGluT2+VGaT+, as well as lateral habenula neurons receiving synaptic input from VTA VGluT2 neurons. Our results provide causal evidence that mice can be used for identifying stressor controllability circuitry and that VTA VGluT2 neurons contribute to transsituational stressor outcomes, such as social avoidance, exaggerated fear, or anxiety-like behavior that are observed within trauma-related disorders.

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.

Dendritic and axonal morphology reflects the input and output of neurons and is a defining feature of neuronal types1,2, yet our knowledge of its diversity remains limited. Here, to systematically examine complete single-neuron morphologies on a brain-wide scale, we established a pipeline encompassing sparse labelling, whole-brain imaging, reconstruction, registration and analysis. We fully reconstructed 1,741 neurons from cortex, claustrum, thalamus, striatum and other brain regions in mice. We identified 11 major projection neuron types with distinct morphological features and corresponding transcriptomic identities. Extensive projectional diversity was found within each of these major types, on the basis of which some types were clustered into more refined subtypes. This diversity follows a set of generalizable principles that govern long-range axonal projections at different levels, including molecular correspondence, divergent or convergent projection, axon termination pattern, regional specificity, topography, and individual cell variability. Although clear concordance with transcriptomic profiles is evident at the level of major projection type, fine-grained morphological diversity often does not readily correlate with transcriptomic subtypes derived from unsupervised clustering, highlighting the need for single-cell cross-modality studies. Overall, our study demonstrates the crucial need for quantitative description of complete single-cell anatomy in cell-type classification, as single-cell morphological diversity reveals a plethora of ways in which different cell types and their individual members may contribute to the configuration and function of their respective circuits.

In the neocortex, subcerebral axonal projections originate largely from layer 5 (L5) extratelencephalic-projecting (ET) neurons. The unique morpho-electric properties of these neurons have been mainly described in rodents, where retrograde tracers or transgenic lines can label them. Similar labeling strategies are infeasible in the human neocortex, rendering the translational relevance of findings in rodents unclear. We leveraged the recent discovery of a transcriptomically defined L5 ET neuron type to study the properties of human L5 ET neurons in neocortical brain slices derived from neurosurgeries. Patch-seq recordings, where transcriptome, physiology, and morphology were assayed from the same cell, revealed many conserved morpho-electric properties of human and rodent L5 ET neurons. Divergent properties were often subtler than differences between L5 cell types within these two species. These data suggest a conserved function of L5 ET neurons in the neocortical hierarchy but also highlight phenotypic divergence possibly related to functional specialization of human neocortex.

Multiplexed fluorescence in situ hybridization techniques have enabled cell-type identification, linking transcriptional heterogeneity with spatial heterogeneity of cells. However, inaccurate cell segmentation reduces the efficacy of cell-type identification and tissue characterization. Here, we present a method called Spot-based Spatial cell-type Analysis by Multidimensional mRNA density estimation (SSAM), a robust cell segmentation-free computational framework for identifying cell-types and tissue domains in 2D and 3D. SSAM is applicable to a variety of in situ transcriptomics techniques and capable of integrating prior knowledge of cell types. We apply SSAM to three mouse brain tissue images: the somatosensory cortex imaged by osmFISH, the hypothalamic preoptic region by MERFISH, and the visual cortex by multiplexed smFISH. Here, we show that SSAM detects regions occupied by known cell types that were previously missed and discovers new cell types.

Mapping the connectivity of diverse neuronal types provides the foundation for understanding the structure and function of neural circuits. High-throughput and low-cost neuroanatomical techniques based on RNA barcode sequencing have the potential to achieve circuit mapping at cellular resolution and a brain-wide scale, but existing Sindbis virus-based techniques can only map long-range projections using anterograde tracing approaches. Rabies virus can complement anterograde tracing approaches by enabling either retrograde labeling of projection neurons or monosynaptic tracing of direct inputs to genetically targeted postsynaptic neurons. However, barcoded rabies virus has so far been only used to map non-neuronal cellular interactions in vivo and synaptic connectivity of cultured neurons. Here we combine barcoded rabies virus with single-cell and in situ sequencing to perform retrograde labeling and transsynaptic labeling in the mouse brain. We sequenced 96 retrogradely labeled cells and 295 transsynaptically labeled cells using single-cell RNAseq, and 4,130 retrogradely labeled cells and 2,914 transsynaptically labeled cells in situ . We determined the transcriptomic identities of rabies virus-infected cells robustly using both single-cell RNA-seq and in situ sequencing. We then distinguished long-range projecting cortical cell types from multiple cortical areas and identified cell types with converging or diverging synaptic connectivity. Combining in situ sequencing with barcoded rabies virus thus complements existing sequencing-based neuroanatomical techniques and provides a potential path for mapping synaptic connectivity of neuronal types at scale.

Cytosine DNA methylation is essential in brain development and has been implicated in various neurological disorders. A comprehensive understanding of DNA methylation diversity across the entire brain in the context of the brain's 3D spatial organization is essential for building a complete molecular atlas of brain cell types and understanding their gene regulatory landscapes. To this end, we employed optimized single-nucleus methylome (snmC-seq3) and multi-omic (snm3C-seq1) sequencing technologies to generate 301,626 methylomes and 176,003 chromatin conformation/methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell type taxonomy that contains 4,673 cell groups and 261 cross-modality-annotated subclasses. We identified millions of differentially methylated regions (DMRs) across the genome, representing potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide multiplexed error-robust fluorescence in situ hybridization (MERFISH2) data validated the association of this spatial epigenetic diversity with transcription and allowed the mapping of the DNA methylation and topology information into anatomical structures more precisely than our dissections. Furthermore, multi-scale chromatin conformation diversities occur in important neuronal genes, highly associated with DNA methylation and transcription changes. Brain-wide cell type comparison allowed us to build a regulatory model for each gene, linking transcription factors, DMRs, chromatin contacts, and downstream genes to establish regulatory networks. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a companion whole-brain SMART-seq3 dataset. Our study establishes the first brain-wide, single-cell resolution DNA methylome and 3D multi-omic atlas, providing an unparalleled resource for comprehending the mouse brain's cellular-spatial and regulatory genome diversity.

The mammalian brain is composed of millions to billions of cells that are organized into numerous cell types with specific spatial distribution patterns and structural and functional properties. An essential step towards understanding brain function is to obtain a parts list, i.e., a catalog of cell types, of the brain. Here, we report a comprehensive and high-resolution transcriptomic and spatial cell type atlas for the whole adult mouse brain. The cell type atlas was created based on the combination of two single-cell-level, whole-brain-scale datasets: a single- cell RNA-sequencing (scRNA-seq) dataset of ∼7 million cells profiled, and a spatially resolved transcriptomic dataset of ∼4.3 million cells using MERFISH. The atlas is hierarchically organized into five nested levels of classification: 7 divisions, 32 classes, 306 subclasses, 1,045 supertypes and 5,200 clusters. We systematically analyzed the neuronal, non-neuronal, and immature neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell type organization in different brain regions, in particular, a dichotomy between the dorsal and ventral parts of the brain: the dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. We also systematically characterized cell-type specific expression of neurotransmitters, neuropeptides, and transcription factors. The study uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types across the brain, suggesting they mediate a myriad of modes of intercellular communications. Finally, we found that transcription factors are major determinants of cell type classification in the adult mouse brain and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole-mouse-brain transcriptomic and spatial cell type atlas establishes a benchmark reference atlas and a foundational resource for deep and integrative investigations of cell type and circuit function, development, and evolution of the mammalian brain.

Biological aging can be defined as a gradual loss of homeostasis across various aspects of molecular and cellular function. Aging is a complex and dynamic process which influences distinct cell types in a myriad of ways. The cellular architecture of the mammalian brain is heterogeneous and diverse, making it challenging to identify precise areas and cell types of the brain that are more susceptible to aging than others. Here, we present a high-resolution single-cell RNA sequencing dataset containing ~1.2 million high-quality single-cell transcriptomic profiles of brain cells from young adult and aged mice across both sexes, including areas spanning the forebrain, midbrain, and hindbrain. We find age-associated gene expression signatures across nearly all 130+ neuronal and non-neuronal cell subclasses we identified. We detect the greatest gene expression changes in non-neuronal cell types, suggesting that different cell types in the brain vary in their susceptibility to aging. We identify specific, age-enriched clusters within specific glial, vascular, and immune cell types from both cortical and subcortical regions of the brain, and specific gene expression changes associated with cell senescence, inflammation, decrease in new myelination, and decreased vasculature integrity. We also identify genes with expression changes across multiple cell subclasses, pointing to certain mechanisms of aging that may occur across wide regions or broad cell types of the brain. Finally, we discover the greatest gene expression changes in cell types localized to the third ventricle of the hypothalamus, including tanycytes, ependymal cells, and Tbx3+ neurons found in the arcuate nucleus that are part of the neuronal circuits regulating food intake and energy homeostasis. These findings suggest that the area surrounding the third ventricle in the hypothalamus may be a hub for aging in the mouse brain. Overall, we reveal a dynamic landscape of cell-type-specific transcriptomic changes in the brain associated with normal aging that will serve as a foundation for the investigation of functional changes in the aging process and the interaction of aging and diseases.

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.

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

Single-cell analyses parse the brain's billions of neurons into thousands of 'cell-type' clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.

We describe a new repressible binary expression system based on the regulatory genes from the Neurospora qa gene cluster. This "Q system" offers attractive features for transgene expression in Drosophila and mammalian cells: low basal expression in the absence of the transcriptional activator QF, high QF-induced expression, and QF repression by its repressor QS. Additionally, feeding flies quinic acid can relieve QS repression. The Q system offers many applications, including (1) intersectional "logic gates" with the GAL4 system for manipulating transgene expression patterns, (2) GAL4-independent MARCM analysis, and (3) coupled MARCM analysis to independently visualize and genetically manipulate siblings from any cell division. We demonstrate the utility of the Q system in determining cell division patterns of a neuronal lineage and gene function in cell growth and proliferation, and in dissecting neurons responsible for olfactory attraction. The Q system can be expanded to other uses in Drosophila and to any organism conducive to transgenesis.

The mammalian Protocadherin (Pcdh) alpha, beta, and gamma gene clusters encode a large family of cadherin-like transmembrane proteins that are differentially expressed in individual neurons. The 22 isoforms of the Pcdhg gene cluster are diversified into A-, B-, and C-types, and the C-type isoforms differ from all other clustered Pcdhs in sequence and expression. Here, we show that mice lacking the three C-type isoforms are phenotypically indistinguishable from the Pcdhg null mutants, displaying virtually identical cellular and synaptic alterations resulting from neuronal apoptosis. By contrast, mice lacking three A-type isoforms exhibit no detectable phenotypes. Remarkably, however, genetically blocking apoptosis rescues the neonatal lethality of the C-type isoform knockouts, but not that of the Pcdhg null mutants. We conclude that the role of the Pcdhg gene cluster in neuronal survival is primarily, if not specifically, mediated by its C-type isoforms, whereas a separate role essential for postnatal development, likely in neuronal wiring, requires isoform diversity.

Proper function of the mammalian brain relies on the establishment of highly specific synaptic connections among billions of neurons. To understand how complex neural circuits function, it is crucial to precisely describe neuronal connectivity and the distributions of synapses to and from individual neurons.

Genetically encoded voltage indicators (GEVIs) enable monitoring of neuronal activity at high spatial and temporal resolution. However, the utility of existing GEVIs has been limited by the brightness and photostability of fluorescent proteins and rhodopsins. We engineered a GEVI, called Voltron, that uses bright and photostable synthetic dyes instead of protein-based fluorophores, thereby extending the number of neurons imaged simultaneously in vivo by a factor of 10 and enabling imaging for significantly longer durations relative to existing GEVIs. We used Voltron for in vivo voltage imaging in mice, zebrafish, and fruit flies. In the mouse cortex, Voltron allowed single-trial recording of spikes and subthreshold voltage signals from dozens of neurons simultaneously over a 15-minute period of continuous imaging. In larval zebrafish, Voltron enabled the precise correlation of spike timing with behavior.

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.

Brain circuits comprise vast numbers of interconnected neurons with diverse molecular, anatomical and physiological properties. To allow targeting of individual neurons for structural and functional studies, we created light-inducible site-specific DNA recombinases based on Cre, Dre and Flp (RecVs). RecVs can induce genomic modifications by one-photon or two-photon light induction in vivo. They can produce targeted, sparse and strong labeling of individual neurons by modifying multiple loci within mouse and zebrafish genomes. In combination with other genetic strategies, they allow intersectional targeting of different neuronal classes. In the mouse cortex they enable sparse labeling and whole-brain morphological reconstructions of individual neurons. Furthermore, these enzymes allow single-cell two-photon targeted genetic modifications and can be used in combination with functional optical indicators with minimal interference. In summary, RecVs enable spatiotemporally precise optogenomic modifications that can facilitate detailed single-cell analysis of neural circuits by linking genetic identity, morphology, connectivity and function.