Abstract:Abstract This paper presents a variational framework for dense diffeomorphic atlas-mapping onto high-throughput histology stacks at the 20 μ m meso-scale. The observed sections are modelled as Gaussian random fields conditioned on a sequence of unknown section by section rigid motions and unknown diffeomorphic transformation of a three-dimensional atlas. To regularize over the high-dimensionality of our parameter space (which is a product space of the rigid motion dimensions and the diffeomorphism dimensions), the histology stacks are modelled as arising from a first order Sobolev space smoothness prior. We show that the joint maximum a-posteriori, penalized-likelihood estimator of our high dimensional parameter space emerges as a joint optimization interleaving rigid motion estimation for histology restacking and large deformation diffeomorphic metric mapping to atlas coordinates. We show that joint optimization in this parameter space solves the classical curvature non-identifiability of the histology stacking problem. The algorithms are demonstrated on a collection of whole-brain histological image stacks from the Mouse Brain Architecture Project. Author Summary New developments in neural tracing techniques have motivated the widespread use of histology as a modality for exploring the circuitry of the brain. Automated mapping of pre-labeled atlases onto modern large datasets of histological imagery is a critical step for elucidating the brain’s neural circuitry and shape. This task is challenging as histological sections are imaged independently and the reconstruction of the unsectioned volume is nontrivial. Typically, neuroanatomists use reference volumes of the same subject (e.g. MRI) to guide reconstruction. However, obtaining reference imagery is often non-standard, as in high-throughput animal models like mouse histology. Others have proposed using anatomical atlases as guides, but have not accounted for the intrinsic nonlinear shape difference from atlas to subject. Our method addresses these limitations by jointly optimizing reconstruction informed by an atlas simultaneously with the nonlinear change of coordinates that encapsulates anatomical variation. This accounts for intrinsic shape differences and enables rigorous, direct comparisons of atlas and subject coordinates. Using simulations, we demonstrate that our method recovers the reconstruction parameters more accurately than atlas-free models and innately produces accurate segmentations from simultaneous atlas mapping. We also demonstrate our method on the Mouse Brain Architecture dataset, successfully mapping and reconstructing over 500 brains.

Image Varifolds on Meshes for Mapping Spatial Transcriptomics

Space-feature measures on meshes for mapping spatial transcriptomics

Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections

xIV-LDDMM Toolkit: A Suite of Image-Varifold Based Technologies for Representing and Mapping 3D Imaging and Spatial-omics Data Simultaneously Across Scales

Solving the problem and quantifying geometric variation in neuroanatomy using generative diffeomorphic mapping

Spatial components of molecular tissue biology

On variational solutions for whole brain serial-section histology using the computational anatomy random orbit model

Probing Brain Micro-architecture by Orientation Distribution Invariant Identification of Diffusion Compartments.

Human Brain Mapping With Conformal Geometry And Multivariate Tensor-Based Morphometry

Subdivision meshes for organizing spatial biomedical data

White Matter Fiber Segmentation Using Functional Varifolds

Machine learning of dissection photographs and surface scanning for quantitative 3D neuropathology

A Diffusion MRI-based Spatiotemporal Continuum of the Embryonic Mouse Brain for Probing Gene-Neuroanatomy Connections

Geometry Processing of Conventionally Produced Mouse Brain Slice Images

Parallelized analysis of spatial gene expression patterns by database integration

Geometry of the cumulant series in neuroimaging

Latent Space Inference For Spatial Transcriptomics

Reconstruction of Micron Resolution Mouse Brain Surface from Large-Scale Imaging Dataset Using Resampling-Based Variational Model

SpatialDE2: Fast and localized variance component analysis of spatial transcriptomics

Insights from spatially mapped gene expression in the mouse brain