CloudReg: automatic terabyte-scale cross-modal brain volume registration

Modular strategies for spatial mapping of diverse cell type data of the mouse brain

Large-scale, international collaborative efforts by members of the BRAIN Initiative Cell Census Network (BICCN) consortium are aggregating the most comprehensive reference database to date for diverse cell type profiling of the mouse brain, which encompasses over 40 different multi-modal profiling techniques from more than 30 research groups. One central challenge for this integrative effort has been the need to map these unique datasets into common reference spaces such that the spatial, structural, and functional information from different cell types can be jointly analyzed. However, significant variation in the acquisition, tissue processing, and imaging techniques across data types makes mapping such diverse data a multifarious problem. Different data types exhibit unique tissue distortion and signal characteristics that precludes a single mapping strategy from being generally applicable across all cell type data. Tailored mapping approaches are often needed to address the unique barriers present in each modality. This work highlights modular atlas mapping strategies developed across separate BICCN studies using the Advanced Normalization Tools Ecosystem (ANTsX) to map spatial transcriptomic (MERFISH) and high-resolution morphology (fMOST) mouse brain data into the Allen Common Coordinate Framework (AllenCCFv3), and developmental (MRI and LSFM) data into the Developmental Common Coordinate Framework (DevCCF). We discuss common mapping strategies that can be shared across modalities and driven by specific challenges from each data type. These mapping strategies include novel open-source contributions that are made publicly available through ANTSX. These include 1) a velocity flow-based approach for continuously mapping developmental trajectories such as that characterizing the DevCCF and 2) an automated framework for determining structural morphology solely through the leveraging of publicly resources. Finally, we provide general guidance to aid investigators to tailor these strategies to address unique data challenges without the need to develop additional specialized software.

RegBoost: Enhancing mouse brain image registration using geometric priors and Laplacian interpolation

We show in this work that incorporating geometric features and geometry processing algorithms for mouse brain image registration broadens the applicability of registration algorithms and improves the registration accuracy of existing methods. We introduce the preprocessing and postprocessing steps in our proposed framework as RegBoost. We develop a method to align the axis of 3D image stacks by detecting the central planes that pass symmetrically through the image volumes. We then find geometric contours by defining external and internal structures to facilitate image correspondences. We establish Dirichlet boundary conditions at these correspondences and find the displacement map throughout the volume using Laplacian interpolation. We discuss the challenges in our standalone framework and demonstrate how our new approaches can improve the results of existing image registration methods. We expect our new approach and algorithms will have critical applications in brain mapping projects.

Preserving Derivative Information while Transforming Neuronal Curves

The international neuroscience community is building the first comprehensive atlases of brain cell types to understand how the brain functions from a higher resolution, and more integrated perspective than ever before. In order to build these atlases, subsets of neurons (e.g. serotonergic neurons, prefrontal cortical neurons etc.) are traced in individual brain samples by placing points along dendrites and axons. Then, the traces are mapped to common coordinate systems by transforming the positions of their points, which neglects how the transformation bends the line segments in between. In this work, we apply the theory of jets to describe how to preserve derivatives of neuron traces up to any order. We provide a framework to compute possible error introduced by standard mapping methods, which involves the Jacobian of the mapping transformation. We show how our first order method improves mapping accuracy in both simulated and real neuron traces under random diffeomorphisms. Our method is freely available in our open-source Python package brainlit.

A protocol for visualization of murine in situ neurovascular interfaces

Mapping cranial vasculature and adjacent neurovascular interfaces in their entirety will enhance our understanding of central nervous system function in any physiologic state. We present a workflow to visualize in situ murine vasculature and surrounding cranial structures using terminal polymer casting of vessels, iterative sample processing and image acquisition, and automated image registration and processing. While this method does not obtain dynamic imaging due to mouse sacrifice, these studies can be performed before sacrifice and processed with other acquired images. For complete details on the use and execution of this protocol, please refer to Rosenblum et al.1.

Preserving Derivative Information while Transforming Neuronal Curves

The international neuroscience community is building the first comprehensive atlases of brain cell types to understand how the brain functions from a higher resolution, and more integrated perspective than ever before. In order to build these atlases, subsets of neurons (e.g. serotonergic neurons, prefrontal cortical neurons etc.) are traced in individual brain samples by placing points along dendrites and axons. Then, the traces are mapped to common coordinate systems by transforming the positions of their points, which neglects how the transformation bends the line segments in between. In this work, we apply the theory of jets to describe how to preserve derivatives of neuron traces up to any order. We provide a framework to compute possible error introduced by standard mapping methods, which involves the Jacobian of the mapping transformation. We show how our first order method improves mapping accuracy in both simulated and real neuron traces under random diffeomorphisms. Our method is freely available in our open-source Python package brainlit.

A guide to the BRAIN Initiative Cell Census Network data ecosystem

Characterizing cellular diversity at different levels of biological organization and across data modalities is a prerequisite to understanding the function of cell types in the brain. Classification of neurons is also essential to manipulate cell types in controlled ways and to understand their variation and vulnerability in brain disorders. The BRAIN Initiative Cell Census Network (BICCN) is an integrated network of data-generating centers, data archives, and data standards developers, with the goal of systematic multimodal brain cell type profiling and characterization. Emphasis of the BICCN is on the whole mouse brain with demonstration of prototype feasibility for human and nonhuman primate (NHP) brains. Here, we provide a guide to the cellular and spatial approaches employed by the BICCN, and to accessing and using these data and extensive resources, including the BRAIN Cell Data Center (BCDC), which serves to manage and integrate data across the ecosystem. We illustrate the power of the BICCN data ecosystem through vignettes highlighting several BICCN analysis and visualization tools. Finally, we present emerging standards that have been developed or adopted toward Findable, Accessible, Interoperable, and Reusable (FAIR) neuroscience. The combined BICCN ecosystem provides a comprehensive resource for the exploration and analysis of cell types in the brain.

Merged magnetic resonance and light sheet microscopy of the whole mouse brain

We have developed workflows to align 3D magnetic resonance histology (MRH) of the mouse brain with light sheet microscopy (LSM) and 3D delineations of the same specimen. We start with MRH of the brain in the skull with gradient echo and diffusion tensor imaging (DTI) at 15 μm isotropic resolution which is ~ 1,000 times higher than that of most preclinical MRI. Connectomes are generated with superresolution tract density images of ~5 μm. Brains are cleared, stained for selected proteins, and imaged by LSM at 1.8 μm/pixel. LSM data are registered into the reference MRH space with labels derived from the ABA common coordinate framework. The result is a high-dimensional integrated volume with registration (HiDiver) with alignment precision better than 50 µm. Throughput is sufficiently high that HiDiver is being used in quantitative studies of the impact of gene variants and aging on mouse brain cytoarchitecture and connectomics.

Preserving Derivative Information while Transforming Neuronal Curves

The international neuroscience community is building the first comprehensive atlases of brain cell types to understand how the brain functions from a higher resolution, and more integrated perspective than ever before. In order to build these atlases, subsets of neurons (e.g. serotonergic neurons, prefrontal cortical neurons etc.) are traced in individual brain samples by placing points along dendrites and axons. Then, the traces are mapped to common coordinate systems by transforming the positions of their points, which neglects how the transformation bends the line segments in between. In this work, we apply the theory of jets to describe how to preserve derivatives of neuron traces up to any order. We provide a framework to compute possible error introduced by standard mapping methods, which involves the Jacobian of the mapping transformation. We show how our first order method improves mapping accuracy in both simulated and real neuron traces, though zeroth order mapping is generally adequate in our real data setting. Our method is freely available in our open-source Python package brainlit.

BrainLine: An Open Pipeline for Connectivity Analysis of Heterogeneous Whole-Brain Fluorescence Volumes

Whole-brain fluorescence images require several stages of computational processing to fully reveal the neuron morphology and connectivity information they contain. However, these computational tools are rarely part of an integrated pipeline. Here we present BrainLine, an open-source pipeline that interfaces with existing software to provide registration, axon segmentation, soma detection, visualization and analysis of results. By implementing a feedback based training paradigm with BrainLine, we were able to use a single learning algorithm to accurately process a diverse set of whole-brain images generated by light-sheet microscopy. BrainLine is available as part of our Python package brainlit: http://brainlit.neurodata.io/ .

Fluorescent transgenic mouse models for whole-brain imaging in health and disease

A paradigm shift is occurring in neuroscience and in general in life sciences converting biomedical research from a descriptive discipline into a quantitative, predictive, actionable science. Living systems are becoming amenable to quantitative description, with profound consequences for our ability to predict biological phenomena. New experimental tools such as tissue clearing, whole-brain imaging, and genetic engineering technologies have opened the opportunity to embrace this new paradigm, allowing to extract anatomical features such as cell number, their full morphology, and even their structural connectivity. These tools will also allow the exploration of new features such as their geometrical arrangement, within and across brain regions. This would be especially important to better characterize brain function and pathological alterations in neurological, neurodevelopmental, and neurodegenerative disorders. New animal models for mapping fluorescent protein-expressing neurons and axon pathways in adult mice are key to this aim. As a result of both developments, relevant cell populations with endogenous fluorescence signals can be comprehensively and quantitatively mapped to whole-brain images acquired at submicron resolution. However, they present intrinsic limitations: weak fluorescent signals, unequal signal strength across the same cell type, lack of specificity of fluorescent labels, overlapping signals in cell types with dense labeling, or undetectable signal at distal parts of the neurons, among others. In this review, we discuss the recent advances in the development of fluorescent transgenic mouse models that overcome to some extent the technical and conceptual limitations and tradeoffs between different strategies. We also discuss the potential use of these strains for understanding disease.

Whole-Brain Single-Cell Imaging and Analysis of Intact Neonatal Mouse Brains Using MRI, Tissue Clearing, and Light-Sheet Microscopy

Tissue clearing followed by light-sheet microscopy (LSFM) enables cellular-resolution imaging of intact brain structure, allowing quantitative analysis of structural changes caused by genetic or environmental perturbations. Whole-brain imaging results in more accurate quantification of cells and the study of region-specific differences that may be missed with commonly used microscopy of physically sectioned tissue. Using light-sheet microscopy to image cleared brains greatly increases acquisition speed as compared to confocal microscopy. Although these images produce very large amounts of brain structural data, most computational tools that perform feature quantification in images of cleared tissue are limited to counting sparse cell populations, rather than all nuclei. Here, we demonstrate NuMorph (Nuclear-Based Morphometry), a group of analysis tools, to quantify all nuclei and nuclear markers within annotated regions of a postnatal day 4 (P4) mouse brain after clearing and imaging on a light-sheet microscope. We describe magnetic resonance imaging (MRI) to measure brain volume prior to shrinkage caused by tissue clearing dehydration steps, tissue clearing using the iDISCO+ method, including immunolabeling, followed by light-sheet microscopy using a commercially available platform to image mouse brains at cellular resolution. We then demonstrate this image analysis pipeline using NuMorph, which is used to correct intensity differences, stitch image tiles, align multiple channels, count nuclei, and annotate brain regions through registration to publicly available atlases. We designed this approach using publicly available protocols and software, allowing any researcher with the necessary microscope and computational resources to perform these techniques. These tissue clearing, imaging, and computational tools allow measurement and quantification of the three-dimensional (3D) organization of cell-types in the cortex and should be widely applicable to any wild-type/knockout mouse study design.

Non-invasive in situ Visualization of the Murine Cranial Vasculature

Understanding physiologic and pathologic central nervous system function depends on our ability to map the entire in situ cranial vasculature and neurovascular interfaces. To accomplish this, we developed a non-invasive workflow to visualize murine cranial vasculature via polymer casting of vessels, iterative sample processing and micro-computed tomography, and automatic deformable image registration, feature extraction, and visualization. This methodology is applicable to any tissue and allows rapid exploration of normal and altered pathologic states.