Publications

An improved understanding of the human lung necessitates advanced systems models informed by an ever-increasing repertoire of molecular omics, cellular, imaging, and pathological datasets. To centralize and standardize information across broad lung research efforts we expanded the LungMAP.net website into a new gateway portal. This portal connects a broad spectrum of research networks, bulk and single-cell multi-omics data and a diverse collection of image data that span mammalian lung development, and disease. The data are standardized across species and technologies using harmonized data and metadata models that leverage recent advances including those from the Human Cell Atlas, diverse ontologies, and the LungMAP CellCards initiative. To cultivate future discoveries, we have aggregated a diverse collection of single-cell atlases for multiple species (human, rhesus, mouse), to enable consistent queries across technologies, cohorts, age, disease, and drug treatment. These atlases are provided as independent and integrated queryable datasets, with an emphasis on dynamic visualization, figure generation, re-analysis, cell-type curation, and automated reference-based classification of user-provided single-cell genomics datasets (Azimuth). As this resource grows, we intend to increase the breadth of available interactive interfaces, supported data types, data portals and datasets from LungMAP and external research efforts.

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.

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.

Understanding of neuronal circuitry at cellular resolution within the brain has relied on neuron tracing methods which involve careful observation and interpretation by experienced neuroscientists. With recent developments in imaging and digitization, this approach is no longer feasible with the large scale (terabyte to petabyte range) images. Machine learning based techniques, using deep networks, provide an efficient alternative to the problem. However, these methods rely on very large volumes of annotated images for training and have error rates that are too high for scientific data analysis, and thus requires a significant volume of human-in-the-loop proofreading. Here we introduce a hybrid architecture combining prior structure in the form of topological data analysis methods, based on discrete Morse theory, with the best-in-class deep-net architectures for the neuronal connectivity analysis. We show significant performance gains using our hybrid architecture on detection of topological structure (e.g. connectivity of neuronal processes and local intensity maxima on axons corresponding to synaptic swellings) with precision/recall close to 90% compared with human observers. We have adapted our architecture to a high performance pipeline capable of semantic segmentation of light microscopic whole-brain image data into a hierarchy of neuronal compartments. We expect that the hybrid architecture incorporating discrete Morse techniques into deep nets will generalize to other data domains.

Until the late twentieth century, it was believed that different sensory modalities were processed by largely independent pathways in the primate cortex, with cross-modal integration only occurring in specialized polysensory areas. This model was challenged by the finding that the peripheral representation of the primary visual cortex (V1) receives monosynaptic connections from areas of the auditory cortex in the macaque. However, auditory projections to V1 have not been reported in other primates. We investigated the existence of direct interconnections between V1 and auditory areas in the marmoset, a New World monkey. Labelled neurons in auditory cortex were observed following 4 out of 10 retrograde tracer injections involving V1. These projections to V1 originated in the caudal subdivisions of auditory cortex (primary auditory cortex, caudal belt and parabelt areas), and targeted parts of V1 that represent parafoveal and peripheral vision. Injections near the representation of the vertical meridian of the visual field labelled few or no cells in auditory cortex. We also placed 8 retrograde tracer injections involving core, belt and parabelt auditory areas, none of which revealed direct projections from V1. These results confirm the existence of a direct, nonreciprocal projection from auditory areas to V1 in a different primate species, which has evolved separately from the macaque for over 30 million years. The essential similarity of these observations between marmoset and macaque indicate that early-stage audiovisual integration is a shared characteristic of primate sensory processing.

Traditionally, the dorsal lateral geniculate nucleus (LGN) and the inferior pulvinar (IPul) nucleus are considered as anatomically and functionally distinct thalamic nuclei. However, in several primate species it has also been established that the koniocellular (K) layers of LGN and parts of the IPul have a shared pattern of immunoreactivity for the calcium-binding protein calbindin. These calbindin-rich cells constitute a thalamic matrix system which is implicated in thalamocortical synchronisation. Further, the K layers and IPul are both involved in visual processing and have similar connections with retina and superior colliculus. Here, we confirmed the continuity between calbindin-rich cells in LGN K layers and the central lateral division of IPul (IPulCL) in marmoset monkeys. By employing a high-throughput neuronal tracing method, we found that both the K layers and IPulCL form comparable patterns of connections with striate and extrastriate cortices; these connections are largely different to those of the parvocellular and magnocellular laminae of LGN. Retrograde tracer-labelled cells and anterograde tracer-labelled axon terminals merged seamlessly from IPulCL into LGN K layers. These results support continuity between LGN K layers and IPulCL, providing an anatomical basis for functional congruity of this region of the dorsal thalamic matrix and calling into question the traditional segregation between LGN and the inferior pulvinar nucleus.

Understanding the connectivity architecture of entire vertebrate brains is a fundamental but difficult task. Here we present an integrated neuro-histological pipeline as well as a grid-based tracer injection strategy for systematic mesoscale connectivity mapping in the common marmoset (Callithrix jacchus). Individual brains are sectioned into  1700 20 µm sections using the tape transfer technique, permitting high quality 3D reconstruction of a series of histochemical stains (Nissl, myelin) interleaved with tracer labeled sections. Systematic in-vivo MRI of the individual animals facilitates injection placement into reference-atlas defined anatomical compartments. Further, by combining the resulting 3D volumes, containing informative cytoarchitectonic markers, with in-vivo and ex-vivo MRI, and using an integrated computational pipeline, we are able to accurately map individual brains into a common reference atlas despite the significant individual variation. This approach will facilitate the systematic assembly of a mesoscale connectivity matrix together with unprecedented 3D reconstructions of brain-wide projection patterns in a primate brain.

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 were 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.

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 were 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.