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.

Effective classification of neuronal cell types requires both molecular and morphological descriptors to be collected in situ at single cell resolution. However, current spatial transcriptomics techniques are not compatible with imaging workflows that successfully reconstruct the morphology of complete axonal projections. Here, we introduce a new methodology that combines tissue clearing, submicron whole-brain two photon imaging, and Expansion-Assisted Iterative Fluorescence In Situ Hybridization (EASI-FISH) to assign molecular identities to fully reconstructed neurons in the mouse brain, which we call morphoFISH. We used morphoFISH to molecularly identify a previously unknown population of cingulate neurons projecting ipsilaterally to the dorsal striatum and contralaterally to higher-order thalamus. By pairing whole-brain morphometry, improved techniques for nucleic acid preservation and spatial gene expression, morphoFISH offers a quantitative solution for discovery of multimodal cell types and complements existing techniques for characterization of increasingly fine-grained cellular heterogeneity in brain circuits.Competing Interest StatementThe authors have declared no competing interest.

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.

HortaCloud is a cloud-based, open-source platform designed to facilitate the collaborative reconstruction of long-range projection neurons from whole-brain light microscopy data. By providing virtual environments directly within the cloud, it eliminates the need for costly and time-consuming data downloads, allowing researchers to work efficiently with terabyte- scale volumetric datasets. Standardization of computational resources in the cloud make deployment easier and more predictable. The pay-as-you-go cloud model reduces adoption barriers by eliminating upfront investments in expensive hardware. Finally, HortaCloud’s decentralized architecture enables global collaboration between researchers and between institutions.