epDevAtlas: Mapping GABAergic cells and microglia in postnatal mouse brains.

During development, brain regions follow encoded growth trajectories. Compared to classical brain growth charts, high-definition growth charts could quantify regional volumetric growth and constituent cell types, improving our understanding of typical and pathological brain development. Here, we create high-resolution 3D atlases of the early postnatal mouse brain, using Allen CCFv3 anatomical labels, at postnatal days (P) 4, 6, 8, 10, 12, and 14, and determine the volumetric growth of different brain regions. We utilize 11 different cell type-specific transgenic animals to validate and refine anatomical labels. Moreover, we reveal region-specific density changes in γ-aminobutyric acid-producing (GABAergic), cortical layer-specific cell types, and microglia as key players in shaping early postnatal brain development. We find contrasting changes in GABAergic neuronal densities between cortical and striatal areas, stabilizing at P12. Moreover, somatostatin-expressing cortical interneurons undergo regionally distinct density reductions, while vasoactive intestinal peptide-expressing interneurons show no significant changes. Remarkably, microglia transition from high density in white matter tracks to gray matter at P10, and show selective density increases in sensory processing areas that correlate with the emergence of individual sensory modalities. Lastly, we create an open-access web-visualization (https://kimlab.io/brain-map/epDevAtlas) for cell-type growth charts and developmental atlases for all postnatal time points.

Sinogram domain angular upsampling of sparse-view micro-CT with dense residual hierarchical transformer and attention-weighted loss.

Reduced angular sampling is a key strategy for increasing scanning efficiency of micron-scale computed tomography (micro-CT). Despite boosting throughput, this strategy introduces noise and extrapolation artifacts due to undersampling. In this work, we present a solution to this issue, by proposing a novel Dense Residual Hierarchical Transformer (DRHT) network to recover high-quality sinograms from 2×, 4× and 8× undersampled scans. DRHT is trained to utilize limited information available from sparsely angular sampled scans and once trained, it can be applied to recover higher-resolution sinograms from shorter scan sessions. Our proposed DRHT model aggregates the benefits of a hierarchical- multi-scale structure along with the combination of local and global feature extraction through dense residual convolutional blocks and non-overlapping window transformer blocks respectively. We also propose a novel noise-aware loss function named KL-L1 to improve sinogram restoration to full resolution. KL-L1, a weighted combination of pixel-level and distribution-level cost functions, leverages inconsistencies in noise distribution and uses learnable spatial weight maps to improve the training of the DRHT model. We present ablation studies and evaluations of our method against other state-of-the-art (SOTA) models over multiple datasets. Our proposed DRHT network achieves an average increase in peak signal to noise ratio (PSNR) of 17.73 dB and a structural similarity index (SSIM) of 0.161, for 8× upsampling, across the three diverse datasets, compared to their respective Bicubic interpolated versions. This novel approach can be utilized to decrease radiation exposure to patients and reduce imaging time for large-scale CT imaging projects.

Developmental Mouse Brain Common Coordinate Framework.

3D standard reference brains serve as key resources to understand the spatial organization of the brain and promote interoperability across different studies. However, unlike the adult mouse brain, the lack of standard 3D reference atlases for developing mouse brains has hindered advancement of our understanding of brain development. Here, we present a multimodal 3D developmental common coordinate framework (DevCCF) spanning mouse embryonic day (E) 11.5, E13.5, E15.5, E18.5, and postnatal day (P) 4, P14, and P56 with anatomical segmentations defined by a developmental ontology. At each age, the DevCCF features undistorted morphologically averaged atlas templates created from Magnetic Resonance Imaging and co-registered high-resolution templates from light sheet fluorescence microscopy. Expert-curated 3D anatomical segmentations at each age adhere to an updated prosomeric model and can be explored via an interactive 3D web-visualizer. As a use case, we employed the DevCCF to unveil the emergence of GABAergic neurons in embryonic brains. Moreover, we integrated the Allen CCFv3 into the P56 template with stereotaxic coordinates and mapped spatial transcriptome cell-type data with the developmental ontology. In summary, the DevCCF is an openly accessible resource that can be used for large-scale data integration to gain a comprehensive understanding of brain development.

Quantitative Geometric Modeling of Blood Cells from X-ray Histotomograms of Whole Zebrafish Larvae.

Tissue phenotyping is foundational to understanding and assessing the cellular aspects of disease in organismal context and an important adjunct to molecular studies in the dissection of gene function, chemical effects, and disease. As a first step toward computational tissue phenotyping, we explore the potential of cellular phenotyping from 3-Dimensional (3D), 0.74 µm isotropic voxel resolution, whole zebrafish larval images derived from X-ray histotomography, a form of micro-CT customized for histopathology. As proof of principle towards computational tissue phenotyping of cells, we created a semi-automated mechanism for the segmentation of blood cells in the vascular spaces of zebrafish larvae, followed by modeling and extraction of quantitative geometric parameters. Manually segmented cells were used to train a random forest classifier for blood cells, enabling the use of a generalized cellular segmentation algorithm for the accurate segmentation of blood cells. These models were used to create an automated data segmentation and analysis pipeline to guide the steps in a 3D workflow including blood cell region prediction, cell boundary extraction, and statistical characterization of 3D geometric and cytological features. We were able to distinguish blood cells at two stages in development (4- and 5-days-post-fertilization) and wild-type vs. polA2 huli hutu ( hht ) mutants. The application of geometric modeling across cell types to and across organisms and sample types may comprise a valuable foundation for computational phenotyping that is more open, informative, rapid, objective, and reproducible.

Quantitative relationship between cerebrovascular network and neuronal cell types in mice.

The cerebrovasculature and its mural cells must meet brain regional energy demands, but how their spatial relationship with different neuronal cell types varies across the brain remains largely unknown. Here we apply brain-wide mapping methods to comprehensively define the quantitative relationships between the cerebrovasculature, capillary pericytes, and glutamatergic and GABAergic neurons, including neuronal nitric oxide synthase-positive (nNOS+) neurons and their subtypes in adult mice. Our results show high densities of vasculature with high fluid conductance and capillary pericytes in primary motor sensory cortices compared with association cortices that show significant positive and negative correlations with energy-demanding parvalbumin+ and vasomotor nNOS+ neurons, respectively. Thalamo-striatal areas that are connected to primary motor sensory cortices also show high densities of vasculature and pericytes, suggesting dense energy support for motor sensory processing areas. Our cellular-resolution resource offers opportunities to examine spatial relationships between the cerebrovascular network and neuronal cell composition in largely understudied subcortical areas.

Whole-Brain Wiring Diagram of Oxytocin System in Adult Mice.

Oxytocin (Oxt) neurons regulate diverse physiological responses via direct connections with different neural circuits. However, the lack of comprehensive input-output wiring diagrams of Oxt neurons and their quantitative relationship with Oxt receptor (Oxtr) expression presents challenges to understanding circuit-specific Oxt functions. Here, we establish a whole-brain distribution and anatomic connectivity map of Oxt neurons, and their relationship with Oxtr expression using high-resolution 3D mapping methods in adult male and female mice. We use a flatmap to describe Oxt neuronal expression in four hypothalamic domains including under-characterized Oxt neurons in the tuberal nucleus (TU). Oxt neurons in the paraventricular hypothalamus (PVH) broadly project to nine functional circuits that control cognition, brain state, and somatic visceral response. In contrast, Oxt neurons in the supraoptic (SO) and accessory (AN) nuclei have limited central projection to a small subset of the nine circuits. Surprisingly, quantitative comparison between Oxt output and Oxtr expression showed no significant correlation across the whole brain, suggesting abundant indirect Oxt signaling in Oxtr-expressing areas. Unlike output, Oxt neurons in both the PVH and SO receive similar monosynaptic inputs from a subset of the nine circuits mainly in the thalamic, hypothalamic, and cerebral nuclei areas. Our results suggest that PVH-Oxt neurons serve as a central modulator to integrate external and internal information via largely reciprocal connection with the nine circuits while the SO-Oxt neurons act mainly as unidirectional Oxt hormonal output. In summary, our Oxt wiring diagram provides anatomic insights about distinct behavioral functions of Oxt signaling in the brain.SIGNIFICANCE STATEMENT Oxytocin (Oxt) neurons regulate diverse physiological functions from prosocial behavior to pain sensation via central projection in the brain. Thus, understanding detailed anatomic connectivity of Oxt neurons can provide insight on circuit-specific roles of Oxt signaling in regulating different physiological functions. Here, we use high-resolution mapping methods to describe the 3D distribution, monosynaptic input and long-range output of Oxt neurons, and their relationship with Oxt receptor (Oxtr) expression across the entire mouse brain. We found Oxt connections with nine functional circuits controlling cognition, brain state, and somatic visceral response. Furthermore, we identified a quantitatively unmatched Oxt-Oxtr relationship, suggesting broad indirect Oxt signaling. Together, our comprehensive Oxt wiring diagram advances our understanding of circuit-specific roles of Oxt neurons.

A wide-field micro-computed tomography detector: micron resolution at half-centimetre scale.

Ideal three-dimensional imaging of complex samples made up of micron-scale structures extending over mm to cm, such as biological tissues, requires both wide field of view and high resolution. For existing optics and detectors used for micro-CT (computed tomography) imaging, sub-micron pixel resolution can only be achieved for fields of view of <2 mm. This article presents a unique detector system with a 6 mm field-of-view image circle and 0.5 µm pixel size that can be used in micro-CT units utilizing both synchrotron and commercial X-ray sources. A resolution-test pattern with linear microstructures and whole adult Daphnia magna were imaged at beamline 8.3.2 of the Berkeley Advanced Light Source. Volumes of 10000 × 10000 × 7096 isotropic 0.5 µm voxels were reconstructed over a 5.0 mm × 3.5 mm field of view. Measurements in the projection domain confirmed a 0.90 µm measured spatial resolution that is largely Nyquist-limited. This unprecedented combination of field of view and resolution dramatically reduces the need for sectional scans and computational stitching for large samples, ultimately offering the means to elucidate changes in tissue and cellular morphology in the context of larger, whole, intact model organisms and specimens. This system is also anticipated to benefit micro-CT imaging in materials science, microelectronics, agricultural science and biomedical engineering.

Whole-organism 3D quantitative characterization of zebrafish melanin by silver deposition micro-CT.

We previously described X-ray histotomography, a high-resolution, non-destructive form of X-ray microtomography (micro-CT) imaging customized for three-dimensional (3D), digital histology, allowing quantitative, volumetric tissue and organismal phenotyping (Ding et al., 2019). Here, we have combined micro-CT with a novel application of ionic silver staining to characterize melanin distribution in whole zebrafish larvae. The resulting images enabled whole-body, computational analyses of regional melanin content and morphology. Normalized micro-CT reconstructions of silver-stained fish consistently reproduced pigment patterns seen by light microscopy, and further allowed direct quantitative comparisons of melanin content across wild-type and mutant samples, including subtle phenotypes not previously noticed. Silver staining of melanin for micro-CT provides proof-of-principle for whole-body, 3D computational phenomic analysis of a specific cell type at cellular resolution, with potential applications in other model organisms and melanocytic neoplasms. Advances such as this in whole-organism, high-resolution phenotyping provide superior context for studying the phenotypic effects of genetic, disease, and environmental variables.

Quantitative cellular-resolution map of the oxytocin receptor in postnatally developing mouse brains.

The oxytocin receptor (OTR) plays critical roles in social behavior development. Despite its significance, brain-wide quantitative understanding of OTR expression remains limited in postnatally developing brains. Here, we develop postnatal 3D template brains to register whole brain images with cellular resolution to systematically quantify OTR cell densities. We utilize fluorescent reporter mice (Otrvenus/+) and find that cortical regions show temporally and spatially heterogeneous patterns with transient postnatal OTR expression without cell death. Cortical OTR cells are largely glutamatergic neurons with the exception of cells in layer 6b. Subcortical regions show similar temporal regulation except the hypothalamus and two hypothalamic nuclei display sexually dimorphic OTR expression. Lack of OTR expression correlates with reduced dendritic spine densities in selected cortical regions of developing brains. Lastly, we create a website to visualize our high-resolution imaging data. In summary, our research provides a comprehensive resource for postnatal OTR expression in the mouse brain.

Enhanced and unified anatomical labeling for a common mouse brain atlas.

Anatomical atlases in standard coordinates are necessary for the interpretation and integration of research findings in a common spatial context. However, the two most-used mouse brain atlases, the Franklin-Paxinos (FP) and the common coordinate framework (CCF) from the Allen Institute for Brain Science, have accumulated inconsistencies in anatomical delineations and nomenclature, creating confusion among neuroscientists. To overcome these issues, we adopt here the FP labels into the CCF to merge the labels in the single atlas framework. We use cell type-specific transgenic mice and an MRI atlas to adjust and further segment our labels. Moreover, detailed segmentations are added to the dorsal striatum using cortico-striatal connectivity data. Lastly, we digitize our anatomical labels based on the Allen ontology, create a web-interface for visualization, and provide tools for comprehensive comparisons between the CCF and FP labels. Our open-source labels signify a key step towards a unified mouse brain atlas.

Computational 3D histological phenotyping of whole zebrafish by X-ray histotomography.

Organismal phenotypes frequently involve multiple organ systems. Histology is a powerful way to detect cellular and tissue phenotypes, but is largely descriptive and subjective. To determine how synchrotron-based X-ray micro-tomography (micro-CT) can yield 3-dimensional whole-organism images suitable for quantitative histological phenotyping, we scanned whole zebrafish, a small vertebrate model with diverse tissues, at ~1 micron voxel resolutions. Micro-CT optimized for cellular characterization (histotomography) allows brain nuclei to be computationally segmented and assigned to brain regions, and cell shapes and volumes to be computed for motor neurons and red blood cells. Striking individual phenotypic variation was apparent from color maps of computed densities of brain nuclei. Unlike histology, the histotomography also allows the study of 3-dimensional structures of millimeter scale that cross multiple tissue planes. We expect the computational and visual insights into 3D cell and tissue architecture provided by histotomography to be useful for reference atlases, hypothesis generation, comprehensive organismal screens, and diagnostics.

Rigid Embedding of Fixed and Stained, Whole, Millimeter-Scale Specimens for Section-free 3D Histology by Micro-Computed Tomography.

For over a hundred years, the histological study of tissues has been the gold standard for medical diagnosis because histology allows all cell types in every tissue to be identified and characterized. Our laboratory is actively working to make technological advances in X-ray micro-computed tomography (micro-CT) that will bring the diagnostic power of histology to the study of full tissue volumes at cellular resolution (i.e., an X-ray Histo-tomography modality). Toward this end, we have made targeted improvements to the sample preparation pipeline. One key optimization, and the focus of the present work, is a straightforward method for rigid embedding of fixed and stained millimeter-scale samples. Many of the published methods for sample immobilization and correlative micro-CT imaging rely on placing the samples in paraffin wax, agarose, or liquids such as alcohol. Our approach extends this work with custom procedures and the design of a 3-dimensional printable apparatus to embed the samples in an acrylic resin directly into polyimide tubing, which is relatively transparent to X-rays. Herein, sample preparation procedures are described for the samples from 0.5 to 10 mm in diameter, which would be suitable for whole zebrafish larvae and juveniles, or other animals and tissue samples of similar dimensions. As proof of concept, we have embedded the specimens from Danio, Drosophila, Daphnia, and a mouse embryo; representative images from 3-dimensional scans for three of these samples are shown. Importantly, our methodology leads to multiple benefits including rigid immobilization, long-term preservation of laboriously-created resources, and the ability to re-interrogate samples.

Comparative analysis of fixation and embedding techniques for optimized histological preparation of zebrafish.

In recognition of the importance of zebrafish as a model organism for studying human disease, we have created zebrafish content for a web-based reference atlas of microanatomy for comparing histology and histopathology between model systems and with humans (http://bio-atlas.psu.edu). Fixation, decalcification, embedding, and sectioning of zebrafish were optimized to maximize section quality. A comparison of protocols involving six fixatives showed that 10% Neutral Buffered Formalin at 21°C for 24h yielded excellent results. Sectioning of juveniles and adults requires bone decalcification; EDTA at 0.35M produced effective decalcification in 21-day-old juveniles through adults (≥~3Months). To improve section plane consistency in sets of larvae, we have developed new array casting molds based on the outside contours of larvae derived from 3D microCT images. Tissue discontinuity in sections, a common barrier to creating quality sections of zebrafish, was minimized by processing and embedding the formalin-fixed zebrafish tissues in plasticized forms of paraffin wax, and by periodic hydration of the block surface in ice water between sets of sections. Optimal H&E (Hematoxylin and Eosin) staining was achieved through refinement of standard protocols. High quality slide scans produced from glass histology slides were digitally processed to maximize image quality, and experimental replicates posted as full slides as part of this publication. Modifications to tissue processing are still needed to eliminate the need for block surface hydration. The further addition of slide collections from other model systems and 3D tools for visualizing tissue architecture would greatly increase the utility of the digital atlas.