High-Throughput Mapping of Long-Range Neuronal Projection Using In Situ Sequencing

Massive Multiplexing of Spatially Resolved Single Neuron Projections with Axonal BARseq

Neurons in the cortex are heterogenous, sending diverse axonal projections to multiple brain regions. Unraveling the logic of these projections requires single-neuron resolution. Although a growing number of techniques have enabled high-throughput reconstruction, these techniques are typically limited to dozens or at most hundreds of neurons per brain, requiring that statistical analyses combine data from different specimens. Here we present axonal BARseq, a high-throughput approach based on reading out nucleic acid barcodes using in situ RNA sequencing, which enables analysis of even densely labeled neurons. As a proof of principle, we have mapped the long-range projections of >8000 mouse primary auditory cortex neurons from a single brain. We identified major cell types based on projection targets and axonal trajectory. The large sample size enabled us to systematically quantify the projections of intratelencephalic (IT) neurons, and revealed that individual IT neurons project to different layers in an area-dependent fashion. Axonal BARseq is a powerful technique for studying the heterogeneity of single neuronal projections at high throughput within individual brains.

Multimodal mapping of cell types and projections in the central nucleus of the amygdala

The central nucleus of the amygdala (CEA) is a brain region that integrates external and internal sensory information and executes innate and adaptive behaviors through distinct output pathways. Despite its complex functions, the diversity of molecularly defined neuronal types in the CEA and their contributions to major axonal projection targets have not been examined systematically. Here, we performed single-cell RNA-sequencing (scRNA-seq) to classify molecularly defined cell types in the CEA and identified marker genes to map the location of these neuronal types using expansion-assisted iterative fluorescence in situ hybridization (EASI-FISH). We developed new methods to integrate EASI-FISH with 5-plex retrograde axonal labeling to determine the spatial, morphological, and connectivity properties of ~30,000 molecularly defined CEA neurons. Our study revealed spatiomolecular organization of the CEA, with medial and lateral CEA associated with distinct molecularly defined cell families. We also found a long-range axon projection network from the CEA, where target regions receive inputs from multiple molecularly defined cell types. Axon collateralization was found primarily among projections to hindbrain targets, which are distinct from forebrain projections. This resource reports marker gene combinations for molecularly defined cell types and axon-projection types, which will be useful for selective interrogation of these neuronal populations to study their contributions to the diverse functions of the CEA.

Activation of MAP2K signaling by genetic engineering or HF-rTMS promotes corticospinal axon sprouting and functional regeneration

Facilitating axon regeneration in the injured central nervous system remains a challenging task. RAF-MAP2K signaling plays a key role in axon elongation during nervous system development. Here, we show that conditional expression of a constitutively kinase-activated BRAF in mature corticospinal neurons elicited the expression of a set of transcription factors previously implicated in the regeneration of zebrafish retinal ganglion cell axons and promoted regeneration and sprouting of corticospinal tract (CST) axons after spinal cord injury in mice. Newly sprouting axon collaterals formed synaptic connections with spinal interneurons, resulting in improved recovery of motor function. Noninvasive suprathreshold high-frequency repetitive transcranial magnetic stimulation (HF-rTMS) activated the BRAF canonical downstream effectors MAP2K1/2 and modulated the expression of a set of regeneration-related transcription factors in a pattern consistent with that induced by BRAF activation. HF-rTMS enabled CST axon regeneration and sprouting, which was abolished in MAP2K1/2 conditional null mice. These data collectively demonstrate a central role of MAP2K signaling in augmenting the growth capacity of mature corticospinal neurons and suggest that HF-rTMS might have potential for treating spinal cord injury by modulating MAP2K signaling.

Application of spatial transcriptome technologies to neurological diseases

Spatial transcriptome technology acquires gene expression profiles while retaining spatial location information, it displays the gene expression properties of cells in situ. Through the investigation of cell heterogeneity, microenvironment, function, and cellular interactions, spatial transcriptome technology can deeply explore the pathogenic mechanisms of cell-type-specific responses and spatial localization in neurological diseases. The present article overviews spatial transcriptome technologies based on microdissection, in situ hybridization, in situ sequencing, in situ capture, and live cell labeling. Each technology is described along with its methods, detection throughput, spatial resolution, benefits, and drawbacks. Furthermore, their applications in neurodegenerative disease, neuropsychiatric illness, stroke and epilepsy are outlined. This information can be used to understand disease mechanisms, pick therapeutic targets, and establish biomarkers.

Deconvolution algorithms for inference of the cell-type composition of the spatial transcriptome

The spatial transcriptome has enabled researchers to resolve transcriptome expression profiles while preserving information about cell location to better understand the complex biological processes that occur in organisms. Due to technical limitations, the current high-throughput spatial transcriptome sequencing methods (known as next-generation sequencing with spatial barcoding methods or spot-based methods) cannot achieve single-cell resolution. A single measurement site, called a spot, in these technologies frequently contains multiple cells of various types. Computational tools for determining the cellular composition of a spot have emerged as a way to break through these limitations. These tools are known as deconvolution tools. Recently, a couple of deconvolution tools based on different strategies have been developed and have shown promise in different aspects. The resulting single-cell resolution expression profiles and/or single-cell composition of spots will significantly affect downstream data mining; thus, it is crucial to choose a suitable deconvolution tool. In this review, we present a list of currently available tools for spatial transcriptome deconvolution, categorize them based on the strategies they employ, and explain their advantages and limitations in detail in order to guide the selection of these tools in future studies.

The continuing evolution of barcode applications: Functional toxicology to cell lineage

DNA barcoding is a method to identify biological entities, including individual cells, tissues, organs, or species, by unique DNA sequences. With the advent of next generation sequencing (NGS), there has been an exponential increase in data acquisition pertaining to medical diagnosis, genetics, toxicology, ecology, cancer, and developmental biology. While barcoding first gained wide access in identifying species, signature tagged mutagenesis has been useful in elucidating gene function, particularly in microbes. With the advent of CRISPR/CAS9, methodology to profile eukaryotic genes has made a broad impact in toxicology and cancer biology. Designed homing guide RNAs (hgRNAs) that self-target DNA sequences facilitate cell lineage barcoding by introducing stochastic mutations within cell identifiers. While each of these applications has their limitations, the potential of sequence barcoding has yet to be realized. This review will focus on signature-tagged mutagenesis and briefly discuss the history of barcoding, experimental problems, novel detection methods, and future directions.

Ascertaining cells' synaptic connections and RNA expression simultaneously with barcoded rabies virus libraries

Brain function depends on synaptic connections between specific neuron types, yet systematic descriptions of synaptic networks and their molecular properties are not readily available. Here, we introduce SBARRO (Synaptic Barcode Analysis by Retrograde Rabies ReadOut), a method that uses single-cell RNA sequencing to reveal directional, monosynaptic relationships based on the paths of a barcoded rabies virus from its "starter" postsynaptic cell to that cell's presynaptic partners. Thousands of these partner relationships can be ascertained in a single experiment, alongside genome-wide RNAs. We use SBARRO to describe synaptic networks formed by diverse mouse brain cell types in vitro, finding that different cell types have presynaptic networks with differences in average size and cell type composition. Patterns of RNA expression suggest that functioning synapses are critical for rabies virus uptake. By tracking individual rabies clones across cells, SBARRO offers new opportunities to map the synaptic organization of neural circuits.

Pooled genetic screens with image‐based profiling

Spatial structure in biology, spanning molecular, organellular, cellular, tissue, and organismal scales, is encoded through a combination of genetic and epigenetic factors in individual cells. Microscopy remains the most direct approach to exploring the intricate spatial complexity defining biological systems and the structured dynamic responses of these systems to perturbations. Genetic screens with deep single‐cell profiling via image features or gene expression programs have the capacity to show how biological systems work in detail by cataloging many cellular phenotypes with one experimental assay. Microscopy‐based cellular profiling provides information complementary to next‐generation sequencing (NGS) profiling and has only recently become compatible with large‐scale genetic screens. Optical screening now offers the scale needed for systematic characterization and is poised for further scale‐up. We discuss how these methodologies, together with emerging technologies for genetic perturbation and microscopy‐based multiplexed molecular phenotyping, are powering new approaches to reveal genotype–phenotype relationships.

High-throughput sequencing of single neuron projections reveals spatial organization in the olfactory cortex

In most sensory modalities, neuronal connectivity reflects behaviorally relevant stimulus features, such as spatial location, orientation, and sound frequency. By contrast, the prevailing view in the olfactory cortex, based on the reconstruction of dozens of neurons, is that connectivity is random. Here, we used high-throughput sequencing-based neuroanatomical techniques to analyze the projections of 5,309 mouse olfactory bulb and 30,433 piriform cortex output neurons at single-cell resolution. Surprisingly, statistical analysis of this much larger dataset revealed that the olfactory cortex connectivity is spatially structured. Single olfactory bulb neurons targeting a particular location along the anterior-posterior axis of piriform cortex also project to matched, functionally distinct, extra-piriform targets. Moreover, single neurons from the targeted piriform locus also project to the same matched extra-piriform targets, forming triadic circuit motifs. Thus, as in other sensory modalities, olfactory information is routed at early stages of processing to functionally diverse targets in a coordinated manner.

Integrating Micro and Nano Technologies for Cell Engineering and Analysis: Toward the Next Generation of Cell Therapy Workflows

The emerging field of cell therapy offers the potential to treat and even cure a diverse array of diseases for which existing interventions are inadequate. Recent advances in micro and nanotechnology have added a multitude of single cell analysis methods to our research repertoire. At the same time, techniques have been developed for the precise engineering and manipulation of cells. Together, these methods have aided the understanding of disease pathophysiology, helped formulate corrective interventions at the cellular level, and expanded the spectrum of available cell therapeutic options. This review discusses how micro and nanotechnology have catalyzed the development of cell sorting, cellular engineering, and single cell analysis technologies, which have become essential workflow components in developing cell-based therapeutics. The review focuses on the technologies adopted in research studies and explores the opportunities and challenges in combining the various elements of cell engineering and single cell analysis into the next generation of integrated and automated platforms that can accelerate preclinical studies and translational research.

Advances in approaches to study cell-type specific cortical circuits throughout development

Neurons in the neocortex and hippocampus are diverse and form synaptic connections that depend on their type. Recent work has improved our understanding of neuronal cell-types and how to target them for experiments. This is crucial for investigating cortical circuit architecture, as the current catalog of established cell-type specific circuit motifs is small relative to the diversity of neuronal subtypes. Some of these motifs are found throughout the cortex, suggesting they are canonical circuits necessary for basic computations. However, the extent to which circuit organization is stereotyped across the brain or varies by cortical region remains unclear. Cortical circuits are also plastic, and their organization evolves throughout each developmental stage. Thus, experimental access to neuronal subtypes with temporal control is essential for studying cortical structure and function. In this mini review, we highlight several recent advances to target specific neuronal subtypes and study their synaptic connectivity and physiology throughout development. We emphasize approaches that combine multiple techniques, provide examples of successful applications, and describe potential future applications of novel tools.

Spatial transcriptomics technology in cancer research

In recent years, spatial transcriptomics (ST) technologies have developed rapidly and have been widely used in constructing spatial tissue atlases and characterizing spatiotemporal heterogeneity of cancers. Currently, ST has been used to profile spatial heterogeneity in multiple cancer types. Besides, ST is a benefit for identifying and comprehensively understanding special spatial areas such as tumor interface and tertiary lymphoid structures (TLSs), which exhibit unique tumor microenvironments (TMEs). Therefore, ST has also shown great potential to improve pathological diagnosis and identify novel prognostic factors in cancer. This review presents recent advances and prospects of applications on cancer research based on ST technologies as well as the challenges.

The landscape of aging

Aging is characterized by a progressive deterioration of physiological integrity, leading to impaired functional ability and ultimately increased susceptibility to death. It is a major risk factor for chronic human diseases, including cardiovascular disease, diabetes, neurological degeneration, and cancer. Therefore, the growing emphasis on “healthy aging” raises a series of important questions in life and social sciences. In recent years, there has been unprecedented progress in aging research, particularly the discovery that the rate of aging is at least partly controlled by evolutionarily conserved genetic pathways and biological processes. In an attempt to bring full-fledged understanding to both the aging process and age-associated diseases, we review the descriptive, conceptual, and interventive aspects of the landscape of aging composed of a number of layers at the cellular, tissue, organ, organ system, and organismal levels.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

What is a cell type and how to define it?

Cell types are the basic functional units of an organism. Cell types exhibit diverse phenotypic properties at multiple levels, making them challenging to define, categorize, and understand. This review provides an overview of the basic principles of cell types rooted in evolution and development and discusses approaches to characterize and classify cell types and investigate how they contribute to the organism's function, using the mammalian brain as a primary example. I propose a roadmap toward a conceptual framework and knowledge base of cell types that will enable a deeper understanding of the dynamic changes of cellular function under healthy and diseased conditions.

The emerging landscape of spatial profiling technologies

Improved scale, multiplexing and resolution are establishing spatial nucleic acid and protein profiling methods as a major pillar for cellular atlas building of complex samples, from tissues to full organisms. Emerging methods yield omics measurements at resolutions covering the nano- to microscale, enabling the charting of cellular heterogeneity, complex tissue architectures and dynamic changes during development and disease. We present an overview of the developing landscape of in situ spatial genome, transcriptome and proteome technologies, exemplify their impact on cell biology and translational research, and discuss current challenges for their community-wide adoption. Among many transformative applications, we envision that spatial methods will map entire organs and enable next-generation pathology.

A transcriptomic axis predicts state modulation of cortical interneurons

Transcriptomics has revealed that cortical inhibitory neurons exhibit a great diversity of fine molecular subtypes^1-6, but it is not known whether these subtypes have correspondingly diverse patterns of activity in the living brain. Here we show that inhibitory subtypes in primary visual cortex (V1) have diverse correlates with brain state, which are organized by a single factor: position along the main axis of transcriptomic variation. We combined in vivo two-photon calcium imaging of mouse V1 with a transcriptomic method to identify mRNA for 72 selected genes in ex vivo slices. We classified inhibitory neurons imaged in layers 1-3 into a three-level hierarchy of 5 subclasses, 11 types and 35 subtypes using previously defined transcriptomic clusters³. Responses to visual stimuli differed significantly only between subclasses, with cells in the Sncg subclass uniformly suppressed, and cells in the other subclasses predominantly excited. Modulation by brain state differed at all hierarchical levels but could be largely predicted from the first transcriptomic principal component, which also predicted correlations with simultaneously recorded cells. Inhibitory subtypes that fired more in resting, oscillatory brain states had a smaller fraction of their axonal projections in layer 1, narrower spikes, lower input resistance and weaker adaptation as determined in vitro⁷, and expressed more inhibitory cholinergic receptors. Subtypes that fired more during arousal had the opposite properties. Thus, a simple principle may largely explain how diverse inhibitory V1 subtypes shape state-dependent cortical processing.

Single cell- and spatial 'Omics revolutionize physiology

Single cell multi- 'Omics and Spatial Transcriptomics are prominent technological highlights of recent years, and both fields still witness a ceaseless firework of novel approaches for high resolution profiling of additional omics layers. As all life processes in organs and organisms are based on the functions of their fundamental building blocks, the individual cells and their interactions, these methods are of utmost worth for the study of physiology in health and disease. Recent discoveries on embryonic development, tumor immunology, the detailed cellular composition and function of complex tissues like for example the kidney or the brain, different roles of the same cell type in different organs, the oncogenic program of individual tumor entities, or the architecture of immunopathology in infected tissue are based on single cell and spatial transcriptomics experiments. In this review, we will give a broad overview of technological concepts for single cell and spatial analysis, showing both advantages and limitations, and illustrate their impact with some particularly impressive case studies.

Advancements in the Quest to Map, Monitor, and Manipulate Neural Circuitry

Neural circuits and the cells that comprise them represent the functional units of the brain. Circuits relay and process sensory information, maintain homeostasis, drive behaviors, and facilitate cognitive functions such as learning and memory. Creating a functionally-precise map of the mammalian brain requires anatomically tracing neural circuits, monitoring their activity patterns, and manipulating their activity to infer function. Advancements in cell-type-specific genetic tools allow interrogation of neural circuits with increased precision. This review provides a broad overview of recombination-based and activity-driven genetic targeting approaches, contemporary viral tracing strategies, electrophysiological recording methods, newly developed calcium, and voltage indicators, and neurotransmitter/neuropeptide biosensors currently being used to investigate circuit architecture and function. Finally, it discusses methods for acute or chronic manipulation of neural activity, including genetically-targeted cellular ablation, optogenetics, chemogenetics, and over-expression of ion channels. With this ever-evolving genetic toolbox, scientists are continuing to probe neural circuits with increasing resolution, elucidating the structure and function of the incredibly complex mammalian brain.

Spatially resolved transcriptomics provide a new method for cancer research

A major feature of cancer is the heterogeneity, both intratumoral and intertumoral. Traditional single-cell techniques have given us a comprehensive understanding of the biological characteristics of individual tumor cells, but the lack of spatial context of the transcriptome has limited the study of cell-to-cell interaction patterns and hindered further exploration of tumor heterogeneity. In recent years, the advent of spatially resolved transcriptomics (SRT) technology has made possible the multidimensional analysis of the tumor microenvironment in the context of intact tissues. Different SRT methods are applicable to different working ranges due to different working principles. In this paper, we review the advantages and disadvantages of various current SRT methods and the overall idea of applying these techniques to oncology studies, hoping to help researchers find breakthroughs. Finally, we discussed the future direction of SRT technology, and deeper investigation into the complex mechanisms of tumor development from different perspectives through multi-omics fusion, paving the way for precisely targeted tumor therapy.

Blind demixing methods for recovering dense neuronal morphology from barcode imaging data

Cellular barcoding methods offer the exciting possibility of 'infinite-pseudocolor' anatomical reconstruction-i.e., assigning each neuron its own random unique barcoded 'pseudocolor,' and then using these pseudocolors to trace the microanatomy of each neuron. Here we use simulations, based on densely-reconstructed electron microscopy microanatomy, with signal structure matched to real barcoding data, to quantify the feasibility of this procedure. We develop a new blind demixing approach to recover the barcodes that label each neuron, and validate this method on real data with known barcodes. We also develop a neural network which uses the recovered barcodes to reconstruct the neuronal morphology from the observed fluorescence imaging data, 'connecting the dots' between discontiguous barcode amplicon signals. We find that accurate recovery should be feasible, provided that the barcode signal density is sufficiently high. This study suggests the possibility of mapping the morphology and projection pattern of many individual neurons simultaneously, at high resolution and at large scale, via conventional light microscopy.

Brain-wide projection reconstruction of single functionally defined neurons

Reconstructing axonal projections of single neurons at the whole-brain level is currently a converging goal of the neuroscience community that is fundamental for understanding the logic of information flow in the brain. Thousands of single neurons from different brain regions have recently been morphologically reconstructed, but the corresponding physiological functional features of these reconstructed neurons are unclear. By combining two-photon Ca²⁺ imaging with targeted single-cell plasmid electroporation, we reconstruct the brain-wide morphologies of single neurons that are defined by a sound-evoked response map in the auditory cortices (AUDs) of awake mice. Long-range interhemispheric projections can be reliably labelled via co-injection with an adeno-associated virus, which enables enhanced expression of indicator protein in the targeted neurons. Here we show that this method avoids the randomness and ambiguity of conventional methods of neuronal morphological reconstruction, offering an avenue for developing a precise one-to-one map of neuronal projection patterns and physiological functional features.

Brain-wide axonal projections of single neurons have been extensively reconstructed without any functional characterization. The authors present a method that allows for developing a precise one-to-one map of both projection patterns and functional features of single neurons in mice.

Emerging strategies for the genetic dissection of gene functions, cell types, and neural circuits in the mammalian brain

The mammalian brain is composed of a large number of highly diverse cell types with different molecular, anatomical, and functional features. Distinct cellular identities are generated during development under the regulation of intricate genetic programs and manifested through unique combinations of gene expression. Recent advancements in our understanding of the molecular and cellular mechanisms underlying the assembly, function, and pathology of the brain circuitry depend on the invention and application of genetic strategies that engage intrinsic gene regulatory mechanisms. Here we review the strategies for gene regulation on DNA, RNA, and protein levels and their applications in cell type targeting and neural circuit dissection. We highlight newly emerged strategies and emphasize the importance of combinatorial approaches. We also discuss the potential caveats and pitfalls in current methods and suggest future prospects to improve their comprehensiveness and versatility.

Cell-type specific transcriptomic signatures of neocortical circuit organization and their relevance to autism

A prevailing challenge in neuroscience is understanding how diverse neuronal cell types select their synaptic partners to form circuits. In the neocortex, major classes of excitatory projection neurons and inhibitory interneurons are conserved across functionally distinct regions. There is evidence these classes form canonical circuit motifs that depend primarily on their identity; however, regional cues likely also influence their choice of synaptic partners. We mined the Allen Institute's single-cell RNA-sequencing database of mouse cortical neurons to study the expression of genes necessary for synaptic connectivity and physiology in two regions: the anterior lateral motor cortex (ALM) and the primary visual cortex (VISp). We used the Allen's metadata to parse cells by clusters representing major excitatory and inhibitory classes that are common to both ALM and VISp. We then performed two types of pairwise differential gene expression analysis: (1) between different neuronal classes within the same brain region (ALM or VISp), and (2) between the same neuronal class in ALM and VISp. We filtered our results for differentially expressed genes related to circuit connectivity and developed a novel bioinformatic approach to determine the sets uniquely enriched in each neuronal class in ALM, VISp, or both. This analysis provides an organized set of genes that may regulate synaptic connectivity and physiology in a cell-type-specific manner. Furthermore, it identifies candidate mechanisms for circuit organization that are conserved across functionally distinct cortical regions or that are region dependent. Finally, we used the SFARI Human Gene Module to identify genes from this analysis that are related to risk for autism spectrum disorder (ASD). Our analysis provides clear molecular targets for future studies to understand neocortical circuit organization and abnormalities that underlie autistic phenotypes.

A genetically defined insula-brainstem circuit selectively controls motivational vigor

The anterior insular cortex (aIC) plays a critical role in cognitive and motivational control of behavior, but the underlying neural mechanism remains elusive. Here, we show that aIC neurons expressing Fezf2 (aICFezf2), which are the pyramidal tract neurons, signal motivational vigor and invigorate need-seeking behavior through projections to the brainstem nucleus tractus solitarii (NTS). aICFezf2 neurons and their postsynaptic NTS neurons acquire anticipatory activity through learning, which encodes the perceived value and the vigor of actions to pursue homeostatic needs. Correspondingly, aIC → NTS circuit activity controls vigor, effort, and striatal dopamine release but only if the action is learned and the outcome is needed. Notably, aICFezf2 neurons do not represent taste or valence. Moreover, aIC → NTS activity neither drives reinforcement nor influences total consumption. These results pinpoint specific functions of aIC → NTS circuit for selectively controlling motivational vigor and suggest that motivation is subserved, in part, by aIC's top-down regulation of dopamine signaling.

Virally encoded connectivity transgenic overlay RNA sequencing (VECTORseq) defines projection neurons involved in sensorimotor integration

Behavior arises from concerted activity throughout the brain. Consequently, a major focus of modern neuroscience is defining the physiology and behavioral roles of projection neurons linking different brain areas. Single-cell RNA sequencing has facilitated these efforts by revealing molecular determinants of cellular physiology and markers that enable genetically targeted perturbations such as optogenetics, but existing methods for sequencing defined projection populations are low throughput, painstaking, and costly. We developed a straightforward, multiplexed approach, virally encoded connectivity transgenic overlay RNA sequencing (VECTORseq). VECTORseq repurposes commercial retrogradely infecting viruses typically used to express functional transgenes (e.g., recombinases and fluorescent proteins) by treating viral transgene mRNA as barcodes within single-cell datasets. VECTORseq is compatible with different viral families, resolves multiple populations with different projection targets in one sequencing run, and identifies cortical and subcortical excitatory and inhibitory projection populations. Our study provides a roadmap for high-throughput identification of neuronal subtypes based on connectivity.

Widening spinal injury research to consider all supraspinal cell types: Why we must and how we can

The supraspinal connectome consists of dozens of neuronal populations that project axons from the brain to the spinal cord to influence a wide range of motor, autonomic, and sensory functions. The complexity and wide distribution of supraspinal neurons present significant technical challenges, leading most spinal cord injury research to focus on a handful of major pathways such as the corticospinal, rubrospinal, and raphespinal. Much less is known about many additional populations that carry information to modulate or compensate for these main pathways, or which carry pre-autonomic and other information of high value to individuals with spinal injury. A confluence of technical developments, however, now enables a whole-connectome study of spinal cord injury. Improved viral labeling, tissue clearing, and automated registration to 3D atlases can quantify supraspinal neurons throughout the murine brain, offering a practical means to track responses to injury and treatment on an unprecedented scale. Here we discuss the need for expanded connectome-wide analyses in spinal injury research, illustrate the potential by discussing a new web-based resource for brain-wide study of supraspinal neurons, and highlight future prospects for connectome analyses.

Temporal controls over inter-areal cortical projection neuron fate diversity

Interconnectivity between neocortical areas is critical for sensory integration and sensorimotor transformations^1-6. These functions are mediated by heterogeneous inter-areal cortical projection neurons (ICPN), which send axon branches across cortical areas as well as to subcortical targets^7-9. Although ICPN are anatomically diverse^10-14, they are molecularly homogeneous¹⁵, and how the diversity of their anatomical and functional features emerge during development remains largely unknown. Here we address this question by linking the connectome and transcriptome in developing single ICPN of the mouse neocortex using a combination of multiplexed analysis of projections by sequencing^16,17 (MAPseq, to identify single-neuron axonal projections) and single-cell RNA sequencing (to identify corresponding gene expression). Focusing on neurons of the primary somatosensory cortex (S1), we reveal a protracted unfolding of the molecular and functional differentiation of motor cortex-projecting ([Formula: see text]) ICPN compared with secondary somatosensory cortex-projecting ([Formula: see text]) ICPN. We identify SOX11 as a temporally differentially expressed transcription factor in [Formula: see text] versus [Formula: see text] ICPN. Postnatal manipulation of SOX11 expression in S1 impaired sensorimotor connectivity and disrupted selective exploratory behaviours in mice. Together, our results reveal that within a single cortical area, different subtypes of ICPN have distinct postnatal paces of molecular differentiation, which are subsequently reflected in distinct circuit connectivities and functions. Dynamic differences in the expression levels of a largely generic set of genes, rather than fundamental differences in the identity of developmental genetic programs, may thus account for the emergence of intra-type diversity in cortical neurons.

Networks of habenula-projecting cortical neurons regulate cocaine seeking

How neurons in the medial prefrontal cortex broadcast stress-relevant information to subcortical brain sites to regulate cocaine relapse remains unclear. The lateral habenula (LHb) serves as a “hub” to filter and propagate stress- and aversion-relevant information in the brain. Here, we show that chemogenetic inhibition of cortical inputs to LHb attenuates relapse-like reinstatement of extinguished cocaine seeking in mice. Using an RNA sequencing–based brain mapping procedure with single-cell resolution, we identify networks of cortical neurons that project to LHb and then preferentially innervate different downstream brain sites, including the ventral tegmental area, median raphe nucleus, and locus coeruleus (LC). By using an intersectional chemogenetics approach, we show that inhibition of cortico-habenular neurons that project to LC, but not to other sites, blocks reinstatement of cocaine seeking. These findings highlight the remarkable complexity of descending cortical inputs to the habenula and identify a cortico-habenulo-hindbrain circuit that regulates cocaine seeking.

Dissecting intercellular and intracellular signaling networks with barcoded genetic tools

The power of next-generation sequencing has stimulated the development of many analysis techniques for transcriptomics and genomics. More recently, the concept of 'molecular barcoding' has broadened the spectrum of sequencing-based applications to dissect different aspects of intracellular and intercellular signaling. In these assay formats, barcode reporters replace standard reporter genes. The virtually infinitive number of expressed barcode sequences allows high levels of multiplexing, hence accelerating experimental progress. Furthermore, reporter barcodes are used to quantitatively monitor a variety of biological events in living cells which has already provided much insight into complex cellular signaling and will further increase our knowledge in the future.

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.

Morphological diversity of single neurons in molecularly defined cell types

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.

A multimodal cell census and atlas of the mammalian primary motor cortex

Here we report the generation of a multimodal cell census and atlas of the mammalian primary motor cortex as the initial product of the BRAIN Initiative Cell Census Network (BICCN). This was achieved by coordinated large-scale analyses of single-cell transcriptomes, chromatin accessibility, DNA methylomes, spatially resolved single-cell transcriptomes, morphological and electrophysiological properties and cellular resolution input-output mapping, integrated through cross-modal computational analysis. Our results advance the collective knowledge and understanding of brain cell-type organization^1-5. First, our study reveals a unified molecular genetic landscape of cortical cell types that integrates their transcriptome, open chromatin and DNA methylation maps. Second, cross-species analysis achieves a consensus taxonomy of transcriptomic types and their hierarchical organization that is conserved from mouse to marmoset and human. Third, in situ single-cell transcriptomics provides a spatially resolved cell-type atlas of the motor cortex. Fourth, cross-modal analysis provides compelling evidence for the transcriptomic, epigenomic and gene regulatory basis of neuronal phenotypes such as their physiological and anatomical properties, demonstrating the biological validity and genomic underpinning of neuron types. We further present an extensive genetic toolset for targeting glutamatergic neuron types towards linking their molecular and developmental identity to their circuit function. Together, our results establish a unifying and mechanistic framework of neuronal cell-type organization that integrates multi-layered molecular genetic and spatial information with multi-faceted phenotypic properties.

A multimodal cell census and atlas of the mammalian primary motor cortex

Here we report the generation of a multimodal cell census and atlas of the mammalian primary motor cortex as the initial product of the BRAIN Initiative Cell Census Network (BICCN). This was achieved by coordinated large-scale analyses of single-cell transcriptomes, chromatin accessibility, DNA methylomes, spatially resolved single-cell transcriptomes, morphological and electrophysiological properties and cellular resolution input-output mapping, integrated through cross-modal computational analysis. Our results advance the collective knowledge and understanding of brain cell-type organization1-5. First, our study reveals a unified molecular genetic landscape of cortical cell types that integrates their transcriptome, open chromatin and DNA methylation maps. Second, cross-species analysis achieves a consensus taxonomy of transcriptomic types and their hierarchical organization that is conserved from mouse to marmoset and human. Third, in situ single-cell transcriptomics provides a spatially resolved cell-type atlas of the motor cortex. Fourth, cross-modal analysis provides compelling evidence for the transcriptomic, epigenomic and gene regulatory basis of neuronal phenotypes such as their physiological and anatomical properties, demonstrating the biological validity and genomic underpinning of neuron types. We further present an extensive genetic toolset for targeting glutamatergic neuron types towards linking their molecular and developmental identity to their circuit function. Together, our results establish a unifying and mechanistic framework of neuronal cell-type organization that integrates multi-layered molecular genetic and spatial information with multi-faceted phenotypic properties.

Epigenomic diversity of cortical projection neurons in the mouse brain

Neuronal cell types are classically defined by their molecular properties, anatomy and functions. Although recent advances in single-cell genomics have led to high-resolution molecular characterization of cell type diversity in the brain1, neuronal cell types are often studied out of the context of their anatomical properties. To improve our understanding of the relationship between molecular and anatomical features that define cortical neurons, here we combined retrograde labelling with single-nucleus DNA methylation sequencing to link neural epigenomic properties to projections. We examined 11,827 single neocortical neurons from 63 cortico-cortical and cortico-subcortical long-distance projections. Our results showed unique epigenetic signatures of projection neurons that correspond to their laminar and regional location and projection patterns. On the basis of their epigenomes, intra-telencephalic cells that project to different cortical targets could be further distinguished, and some layer 5 neurons that project to extra-telencephalic targets (L5 ET) formed separate clusters that aligned with their axonal projections. Such separation varied between cortical areas, which suggests that there are area-specific differences in L5 ET subtypes, which were further validated by anatomical studies. Notably, a population of cortico-cortical projection neurons clustered with L5 ET rather than intra-telencephalic neurons, which suggests that a population of L5 ET cortical neurons projects to both targets. We verified the existence of these neurons by dual retrograde labelling and anterograde tracing of cortico-cortical projection neurons, which revealed axon terminals in extra-telencephalic targets including the thalamus, superior colliculus and pons. These findings highlight the power of single-cell epigenomic approaches to connect the molecular properties of neurons with their anatomical and projection properties.

Multiplexing Methods for Simultaneous Large-Scale Transcriptomic Profiling of Samples at Single-Cell Resolution

Barcoding technology has greatly improved the throughput of cells and genes detected in single-cell RNA sequencing (scRNA-seq) studies. Recently, increasing studies have paid more attention to the use of this technology to increase the throughput of samples, as it has greatly reduced the processing time, technical batch effects, and library preparation costs, and lowered the per-sample cost. In this review, the various DNA-based barcoding methods for sample multiplexing are focused on, specifically, on the four major barcoding strategies. A detailed comparison of the barcoding methods is also presented, focusing on aspects such as sample/cell throughput and gene detection, and guidelines for choosing the most appropriate barcoding technique according to the personalized requirements are developed. Finally, the critical applications of sample multiplexing and technical challenges in combinatorial labeling, barcoding in vivo, and multimodal tagging at the spatially resolved resolution, as well as, the future prospects of multiplexed scRNA-seq, for example, prioritizing and predicting the severity of coronavirus disease 2019 (COVID-19) in patients of different gender and age are highlighted.

Single-Cell Mapping of GLP-1 and GIP Receptor Expression in the Dorsal Vagal Complex

The dorsal vagal complex (DVC) in the hindbrain, composed of the area postrema, nucleus of the solitary tract, and dorsal motor nucleus of the vagus, plays a critical role in modulating satiety. The incretins glucagon-like peptide 1 (GLP-1) and glucose-dependent insulinotropic polypeptide (GIP) act directly in the brain to modulate feeding, and receptors for both are expressed in the DVC. Given the impressive clinical responses to pharmacologic manipulation of incretin signaling, understanding the central mechanisms by which incretins alter metabolism and energy balance is of critical importance. Here, we review recent single-cell approaches used to detect molecular signatures of GLP-1 and GIP receptor-expressing cells in the DVC. In addition, we discuss how current advancements in single-cell transcriptomics, epigenetics, spatial transcriptomics, and circuit mapping techniques have the potential to further characterize incretin receptor circuits in the hindbrain.

Whole-body integration of gene expression and single-cell morphology

Animal bodies are composed of cell types with unique expression programs that implement their distinct locations, shapes, structures, and functions. Based on these properties, cell types assemble into specific tissues and organs. To systematically explore the link between cell-type-specific gene expression and morphology, we registered an expression atlas to a whole-body electron microscopy volume of the nereid Platynereis dumerilii. Automated segmentation of cells and nuclei identifies major cell classes and establishes a link between gene activation, chromatin topography, and nuclear size. Clustering of segmented cells according to gene expression reveals spatially coherent tissues. In the brain, genetically defined groups of neurons match ganglionic nuclei with coherent projections. Besides interneurons, we uncover sensory-neurosecretory cells in the nereid mushroom bodies, which thus qualify as sensory organs. They furthermore resemble the vertebrate telencephalon by molecular anatomy. We provide an integrated browser as a Fiji plugin for remote exploration of all available multimodal datasets.

•

A cellular atlas integrates gene expression and ultrastructure for an entire annelid

•

Morphometry of all segmented cells, nuclei, and chromatin categorizes cell classes

•

Molecular anatomy and projectome of head ganglionic nuclei and mushroom bodies

•

An open-source browser for multimodal big image data exploration and analysis

Exploring tissue architecture using spatial transcriptomics

Deciphering the principles and mechanisms by which gene activity orchestrates complex cellular arrangements in multicellular organisms has far-reaching implications for research in the life sciences. Recent technological advances in next-generation sequencing- and imaging-based approaches have established the power of spatial transcriptomics to measure expression levels of all or most genes systematically throughout tissue space, and have been adopted to generate biological insights in neuroscience, development and plant biology as well as to investigate a range of disease contexts, including cancer. Similar to datasets made possible by genomic sequencing and population health surveys, the large-scale atlases generated by this technology lend themselves to exploratory data analysis for hypothesis generation. Here we review spatial transcriptomic technologies and describe the repertoire of operations available for paths of analysis of the resulting data. Spatial transcriptomics can also be deployed for hypothesis testing using experimental designs that compare time points or conditions-including genetic or environmental perturbations. Finally, spatial transcriptomic data are naturally amenable to integration with other data modalities, providing an expandable framework for insight into tissue organization.

Eavesdropping wires: Recording activity in axons using genetically encoded calcium indicators

Neurons broadcast electrical signals to distal brain regions through extensive axonal arbors. Genetically encoded calcium sensors permit the direct observation of action potential activity at axonal terminals, providing unique insights on the organization and function of neural projections. Here, we consider what information can be gleaned from axonal recordings made at scales ranging from the summed activity extracted from multi-cell axon projections to single boutons. In particular, we discuss the application of different recently developed multi photon and fiber photometry methods for recording neural activity in axons of rodents. We define experimental difficulties associated with imaging approaches in the axonal compartment and highlight the latest methodological advances for addressing these issues. Finally, we reflect on ways in which new technologies can be used in conjunction with axon calcium imaging to address current questions in neurobiology.

Over-representation of fundamental decision variables in the prefrontal cortex underlies decision bias

The brain is organized into anatomically distinct structures consisting of a variety of projection neurons. While such evolutionarily conserved neural circuit organization underlies the innate ability of animals to swiftly adapt to environments, they can cause biased cognition and behavior. Although recent studies have begun to address the causal importance of projection-neuron types as distinct computational units, it remains unclear how projection types are functionally organized in encoding variables during cognitive tasks. This review focuses on the neural computation of decision making in the prefrontal cortex and discusses what decision variables are encoded by single neurons, neuronal populations, and projection type, alongside how specific projection types constrain decision making. We focus particularly on "over-representations" of distinct decision variables in the prefrontal cortex that reflect the biological and subjective significance of the variables for the decision makers. We suggest that task-specific over-representation in the prefrontal cortex involves the refinement of the given decision making, while generalized over-representation of fundamental decision variables is associated with suboptimal decision biases, including pathological ones such as those in patients with psychiatric disorders. Such over-representation of the fundamental decision variables in the prefrontal cortex appear to be tightly constrained by afferent and efferent connections that can be optogenetically intervened on. These ideas may provide critical insights into potential therapeutic targets for psychiatric disorders, including addiction and depression.

Assessing the replicability of spatial gene expression using atlas data from the adult mouse brain

High-throughput, spatially resolved gene expression techniques are poised to be transformative across biology by overcoming a central limitation in single-cell biology: the lack of information on relationships that organize the cells into the functional groupings characteristic of tissues in complex multicellular organisms. Spatial expression is particularly interesting in the mammalian brain, which has a highly defined structure, strong spatial constraint in its organization, and detailed multimodal phenotypes for cells and ensembles of cells that can be linked to mesoscale properties such as projection patterns, and from there, to circuits generating behavior. However, as with any type of expression data, cross-dataset benchmarking of spatial data is a crucial first step. Here, we assess the replicability, with reference to canonical brain subdivisions, between the Allen Institute's in situ hybridization data from the adult mouse brain (Allen Brain Atlas (ABA)) and a similar dataset collected using spatial transcriptomics (ST). With the advent of tractable spatial techniques, for the first time, we are able to benchmark the Allen Institute's whole-brain, whole-transcriptome spatial expression dataset with a second independent dataset that similarly spans the whole brain and transcriptome. We use regularized linear regression (LASSO), linear regression, and correlation-based feature selection in a supervised learning framework to classify expression samples relative to their assayed location. We show that Allen Reference Atlas labels are classifiable using transcription in both data sets, but that performance is higher in the ABA than in ST. Furthermore, models trained in one dataset and tested in the opposite dataset do not reproduce classification performance bidirectionally. While an identifying expression profile can be found for a given brain area, it does not generalize to the opposite dataset. In general, we found that canonical brain area labels are classifiable in gene expression space within dataset and that our observed performance is not merely reflecting physical distance in the brain. However, we also show that cross-platform classification is not robust. Emerging spatial datasets from the mouse brain will allow further characterization of cross-dataset replicability ultimately providing a valuable reference set for understanding the cell biology of the brain.

Progress towards a cellularly resolved mouse mesoconnectome is empowered by data fusion and new neuroanatomy techniques

Over the past decade there has been a rapid improvement in techniques for obtaining large-scale cellular level data related to the mouse brain connectome. However, a detailed mapping of cell-type-specific projection patterns is lacking, which would, for instance, allow us to study the role of circuit motifs in cognitive processes. In this work, we review advanced neuroanatomical and data fusion techniques within the context of a proposed Multimodal Connectomic Integration Framework for augmenting the cellularly resolved mouse mesoconnectome. First, we emphasize the importance of registering data modalities to a common reference atlas. We then review a number of novel experimental techniques that can provide data for characterizing cell-types in the mouse brain. Furthermore, we examine a number of data integration strategies, which involve fine-grained cell-type classification, spatial inference of cell densities, latent variable models for the mesoconnectome and multi-modal factorisation. Finally, we discuss a number of use cases which depend on connectome augmentation techniques, such as model simulations of functional connectivity and generating mechanistic hypotheses for animal disease models.

Integrating barcoded neuroanatomy with spatial transcriptional profiling enables identification of gene correlates of projections

Functional circuits consist of neurons with diverse axonal projections and gene expression. Understanding the molecular signature of projections requires high-throughput interrogation of both gene expression and projections to multiple targets in the same cells at cellular resolution, which is difficult to achieve using current technology. Here, we introduce BARseq2, a technique that simultaneously maps projections and detects multiplexed gene expression by in situ sequencing. We determined the expression of cadherins and cell-type markers in 29,933 cells and the projections of 3,164 cells in both the mouse motor cortex and auditory cortex. Associating gene expression and projections in 1,349 neurons revealed shared cadherin signatures of homologous projections across the two cortical areas. These cadherins were enriched across multiple branches of the transcriptomic taxonomy. By correlating multigene expression and projections to many targets in single neurons with high throughput, BARseq2 provides a potential path to uncovering the molecular logic underlying neuronal circuits.

Visual intracortical and transthalamic pathways carry distinct information to cortical areas

Sensory processing involves information flow between neocortical areas, assumed to rely on direct intracortical projections. However, cortical areas may also communicate indirectly via higher-order nuclei in the thalamus, such as the pulvinar or lateral posterior nucleus (LP) in the visual system of rodents. The fine-scale organization and function of these cortico-thalamo-cortical pathways remains unclear. We find that responses of mouse LP neurons projecting to higher visual areas likely derive from feedforward input from primary visual cortex (V1) combined with information from many cortical and subcortical areas, including superior colliculus. Signals from LP projections to different higher visual areas are tuned to specific features of visual stimuli and their locomotor context, distinct from the signals carried by direct intracortical projections from V1. Thus, visual transthalamic pathways are functionally specific to their cortical target, different from feedforward cortical pathways, and combine information from multiple brain regions, linking sensory signals with behavioral context.

•

Transthalamic pathway through pulvinar indirectly connects lower to higher cortical areas

•

This pathway combines input from V1 with that of many cortical and subcortical areas

•

Pulvinar conveys distinct visual and motor information to different higher visual areas

•

Direct intracortical and transthalamic pathways convey different information

BARcode DEmixing through Non-negative Spatial Regression (BarDensr)

Modern spatial transcriptomics methods can target thousands of different types of RNA transcripts in a single slice of tissue. Many biological applications demand a high spatial density of transcripts relative to the imaging resolution, leading to partial mixing of transcript rolonies in many voxels; unfortunately, current analysis methods do not perform robustly in this highly-mixed setting. Here we develop a new analysis approach, BARcode DEmixing through Non-negative Spatial Regression (BarDensr): we start with a generative model of the physical process that leads to the observed image data and then apply sparse convex optimization methods to estimate the underlying (demixed) rolony densities. We apply BarDensr to simulated and real data and find that it achieves state of the art signal recovery, particularly in densely-labeled regions or data with low spatial resolution. Finally, BarDensr is fast and parallelizable. We provide open-source code as well as an implementation for the 'NeuroCAAS' cloud platform.