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.

Hemodynamic signals are widely used to infer neural activity in the brain. We tested the hypothesis that hemodynamic signals faithfully report neural activity during voluntary behaviors by measuring cerebral blood volume (CBV) and neural activity in the somatosensory cortex and frontal cortex of head-fixed mice during locomotion. Locomotion induced a large and robust increase in firing rate and gamma-band (40 -100 Hz) power in the local field potential in the limb representations in somatosensory cortex, and was accompanied by increases in CBV, demonstrating that hemodynamic signals are coupled with neural activity in this region. However, in the frontal cortex, CBV did not change during locomotion, but firing rate and gamma-band power both increased, indicating a decoupling of neural activity from the hemodynamic signal. These results show that hemodynamic signals are not faithful indicators of the mean neural activity in the frontal cortex during locomotion; thus, the results from fMRI and other hemodynamic imaging methodologies for studying neural processes must be interpreted with caution.

Voluntary locomotion is accompanied by large increases in cortical activity and localized increases in cerebral blood volume (CBV). We sought to quantitatively determine the spatial and temporal dynamics of voluntary locomotion-evoked cerebral hemodynamic changes. We measured single vessel dilations using two-photon microscopy and cortex-wide changes in CBV-related signal using intrinsic optical signal (IOS) imaging in head-fixed mice freely locomoting on a spherical treadmill. During bouts of locomotion, arteries dilated rapidly, while veins distended slightly and recovered slowly. The dynamics of diameter changes of both vessel types could be captured using a simple linear convolution model. Using these single vessel measurements, we developed a novel analysis approach to separate out spatially and temporally distinct arterial and venous components of the location-specific hemodynamic response functions (HRF) for IOS. The HRF of each pixel of was well fit by a sum of a fast arterial and a slow venous component. The HRFs of pixels in the limb representations of somatosensory cortex had a large arterial contribution, while in the frontal cortex the arterial contribution to the HRF was negligible. The venous contribution was much less localized, and was substantial in the frontal cortex. The spatial pattern and amplitude of these HRFs in response to locomotion in the cortex was robust across imaging sessions. Separating the more localized, arterial component from the diffuse venous signals will be useful for dealing with the dynamic signals generated by naturalistic stimuli.

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.

Understanding how changes in the cardiovascular system contribute to cerebral blood flow (CBF) and volume (CBV) increases is critical for interpreting hemodynamic signals. Here we investigated how systemic cardiovascular changes affect the cortical hemodynamic response during voluntary locomotion. In the mouse, voluntary locomotion drives an increase in cortical CBF and arterial CBV that is localized to the forelimb/hindlimb representation in the somatosensory cortex, as well as a diffuse venous CBV increase. To determine if the heart rate increases that accompany locomotion contribute to locomotion-induced CBV and CBF increases, we occluded heart rate increases with the muscarinic cholinergic receptor antagonist glycopyrrolate, and reduced heart rate with the β1-adrenergic receptor antagonist atenolol. We quantified the effects of these cardiovascular manipulations on CBV and CBF dynamics by comparing the hemodynamic response functions (HRF) to locomotion across these conditions. Neither the CBF HRF nor the arterial component of the CBV HRF was significantly affected by pharmacological disruption of the heart rate. In contrast, the amplitude and spatial extent of the venous component of the CBV HRF was decreased by atenolol. These results suggest that the increase in venous CBV during locomotion was partially driven by peripheral cardiovascular changes, whereas CBF and arterial CBV increases associated with locomotion reflect central processes.