The structure of axonal arbors controls how signals from individual neurons are routed within the mammalian brain. However, the arbors of very few long-range projection neurons have been reconstructed in their entirety, as axons with diameters as small as 100 nm arborize in target regions dispersed over many millimeters of tissue. We introduce a platform for high-resolution, three-dimensional fluorescence imaging of complete tissue volumes that enables the visualization and reconstruction of long-range axonal arbors. This platform relies on a high-speed two-photon microscope integrated with a tissue vibratome and a suite of computational tools for large-scale image data. We demonstrate the power of this approach by reconstructing the axonal arbors of multiple neurons in the motor cortex across a single mouse brain.

The thalamus is the central communication hub of the forebrain and provides the cerebral cortex with inputs from sensory organs, subcortical systems and the cortex itself. Multiple thalamic regions send convergent information to each cortical region, but the organizational logic of thalamic projections has remained elusive. Through comprehensive transcriptional analyses of retrogradely labeled thalamic neurons in adult mice, we identify three major profiles of thalamic pathways. These profiles exist along a continuum that is repeated across all major projection systems, such as those for vision, motor control and cognition. The largest component of gene expression variation in the mouse thalamus is topographically organized, with features conserved in humans. Transcriptional differences between these thalamic neuronal identities are tied to cellular features that are critical for function, such as axonal morphology and membrane properties. Molecular profiling therefore reveals covariation in the properties of thalamic pathways serving all major input modalities and output targets, thus establishing a molecular framework for understanding the thalamus.

The microvasculature underlies the supply networks that support neuronal activity within heterogeneous brain regions. What are common versus heterogeneous aspects of the connectivity, density, and orientation of capillary networks? To address this, we imaged, reconstructed, and analyzed the microvasculature connectome in whole adult mice brains with sub-micrometer resolution. Graph analysis revealed common network topology across the brain that leads to a shared structural robustness against the rarefaction of vessels. Geometrical analysis, based on anatomically accurate reconstructions, uncovered a scaling law that links length density, i.e., the length of vessel per volume, with tissue-to-vessel distances. We then derive a formula that connects regional differences in metabolism to differences in length density and, further, predicts a common value of maximum tissue oxygen tension across the brain. Last, the orientation of capillaries is weakly anisotropic with the exception of a few strongly anisotropic regions; this variation can impact the interpretation of fMRI data.

Behavior requires neural activity across the brain, but most experiments probe neurons in a single area at a time. Here we used multiple Neuropixels probes to record neural activity simultaneously in brain-wide circuits, in mice performing a memory-guided directional licking task. We targeted brain areas that form multi-regional loops with anterior lateral motor cortex (ALM), a key circuit node mediating the behavior. Neurons encoding sensory stimuli, choice, and actions were distributed across the brain. However, in addition to ALM, coding of choice was concentrated in subcortical areas receiving input from ALM, in an ALM-dependent manner. Choice signals were first detected in ALM and the midbrain, followed by the thalamus, and other brain areas. At the time of movement initiation, choice-selective activity collapsed across the brain, followed by new activity patterns driving specific actions. Our experiments provide the foundation for neural circuit models of decision-making and movement initiation.

The striatum shows general topographic organization and regional differences in behavioral functions. How corticostriatal topography differs across cortical areas and cell types to support these distinct functions is unclear. This study contrasted corticostriatal projections from two layer 5 cell types, intratelencephalic (IT-type) and pyramidal tract (PT-type) neurons, using viral vectors expressing fluorescent reporters in Cre-driver mice. Long-range corticostriatal projections from sensory and motor cortex are somatotopic, with a decreasing somatotopic specificity as injections move from sensory to motor and frontal areas. Somatotopic organization differs between IT-type and PT-type neurons, including injections in the same site, with IT-type neurons having higher somatotopic stereotypy than PT-type neurons. Furthermore, IT-type projections from interconnected cortical areas have stronger correlations in corticostriatal targeting than PT-type projections do. Thus, as predicted by a long-standing basal ganglia model, corticostriatal projections of interconnected cortical areas form parallel circuits in basal ganglia-thalamus-cortex loops.

The mammalian hippocampus, comprised of serially connected subfields, participates in diverse behavioral and cognitive functions. It has been postulated that parallel circuitry embedded within hippocampal subfields may underlie such functional diversity. We sought to identify, delineate, and manipulate this putatively parallel architecture in the dorsal subiculum, the primary output subfield of the dorsal hippocampus. Population and single-cell RNA-seq revealed that the subiculum can be divided into two spatially adjacent subregions associated with prominent differences in pyramidal cell gene expression. Pyramidal cells occupying these two regions differed in their long-range inputs, local wiring, projection targets, and electrophysiological properties. Leveraging gene-expression differences across these regions, we use genetically restricted neuronal silencing to show that these regions differentially contribute to spatial working memory. This work provides a coherent molecular-, cellular-, circuit-, and behavioral-level demonstration that the hippocampus embeds structurally and functionally dissociable streams within its serial architecture.

Activity in the motor cortex predicts movements, seconds before they are initiated. This preparatory activity has been observed across cortical layers, including in descending pyramidal tract neurons in layer 5. A key question is how preparatory activity is maintained without causing movement, and is ultimately converted to a motor command to trigger appropriate movements. Here, using single-cell transcriptional profiling and axonal reconstructions, we identify two types of pyramidal tract neuron. Both types project to several targets in the basal ganglia and brainstem. One type projects to thalamic regions that connect back to motor cortex; populations of these neurons produced early preparatory activity that persisted until the movement was initiated. The second type projects to motor centres in the medulla and mainly produced late preparatory activity and motor commands. These results indicate that two types of motor cortex output neurons have specialized roles in motor control.

The neuronal circuits defined by the axonal projections of pyramidal neurons in the cerebral cortex are responsible for processing sensory and other information to plan and execute behavior. Subtypes of cortical pyramidal neurons are organized across layers, with those in different layers distinguished by their patterns of axonal projections and connectivity. For example, those in layers 2 and 3 project between cortical areas to integrate sensory and other information with motor areas; while those in layers 5 and 6 also integrate information between cortical areas, but also project to subcortical structures involved in the generation of behavior. Recent advances in neuroanatomical techniques allow one to target specific subtypes of cortical pyramidal neurons and label both their inputs and projections. Combining these methods with neurophysiological recording techniques and newly introduced atlases of the mouse brain provide the opportunity to achieve a detailed view of the organization of cerebral cortical circuits.

Understanding the principles governing neuronal diversity is a fundamental goal for neuroscience. Here we provide an anatomical and transcriptomic database of nearly 200 genetically identified cell populations. By separately analyzing the robustness and pattern of expression differences across these cell populations, we identify two gene classes contributing distinctly to neuronal diversity. Short homeobox transcription factors distinguish neuronal populations combinatorially, and exhibit extremely low transcriptional noise, enabling highly robust expression differences. Long neuronal effector genes, such as channels and cell adhesion molecules, contribute disproportionately to neuronal diversity, based on their patterns rather than robustness of expression differences. By linking transcriptional identity to genetic strains and anatomical atlases we provide an extensive resource for further investigation of mouse neuronal cell types.