The molecular cytoarchitecture of the adult mouse brain.

The function of the mammalian brain relies upon the specification and spatial positioning of diversely specialized cell types. Yet, the molecular identities of the cell types and their positions within individual anatomical structures remain incompletely known. To construct a comprehensive atlas of cell types in each brain structure, we paired high-throughput single-nucleus RNA sequencing with Slide-seq1,2-a recently developed spatial transcriptomics method with near-cellular resolution-across the entire mouse brain. Integration of these datasets revealed the cell type composition of each neuroanatomical structure. Cell type diversity was found to be remarkably high in the midbrain, hindbrain and hypothalamus, with most clusters requiring a combination of at least three discrete gene expression markers to uniquely define them. Using these data, we developed a framework for genetically accessing each cell type, comprehensively characterized neuropeptide and neurotransmitter signalling, elucidated region-specific specializations in activity-regulated gene expression and ascertained the heritability enrichment of neurological and psychiatric phenotypes. These data, available as an online resource ( www.BrainCellData.org ), should find diverse applications across neuroscience, including the construction of new genetic tools and the prioritization of specific cell types and circuits in the study of brain diseases.

Toward Community-Driven Big Open Brain Science: Open Big Data and Tools for Structure, Function, and Genetics.

As acquiring bigger data becomes easier in experimental brain science, computational and statistical brain science must achieve similar advances to fully capitalize on these data. Tackling these problems will benefit from a more explicit and concerted effort to work together. Specifically, brain science can be further democratized by harnessing the power of community-driven tools, which both are built by and benefit from many different people with different backgrounds and expertise. This perspective can be applied across modalities and scales and enables collaborations across previously siloed communities.

Machine learning approach for obstructive sleep apnea screening using brain diffusion tensor imaging.

Patients with obstructive sleep apnea (OSA) show autonomic, mood, cognitive, and breathing dysfunctions that are linked to increased morbidity and mortality, which can be improved with early screening and intervention. The gold standard and other available methods for OSA diagnosis are complex, require whole-night data, and have significant wait periods that potentially delay intervention. Our aim was to examine whether using faster and less complicated machine learning models, including support vector machine (SVM) and random forest (RF), with brain diffusion tensor imaging (DTI) data can classify OSA from healthy controls. We collected two DTI series from 59 patients with OSA [age: 50.2 ± 9.9 years; body mass index (BMI): 31.5 ± 5.6 kg/m2 ; apnea-hypopnea index (AHI): 34.1 ± 21.2 events/h 23 female] and 96 controls (age: 51.8 ± 9.7 years; BMI: 26.2 ± 4.1 kg/m2 ; 51 female) using a 3.0-T magnetic resonance imaging scanner. Using DTI data, mean diffusivity maps were calculated from each series, realigned and averaged, normalised to a common space, and used to conduct cross-validation for model training and selection and to predict OSA. The RF model showed 0.73 OSA and controls classification accuracy and 0.85 area under the curve (AUC) value on the receiver-operator curve. Cross-validation showed the RF model with comparable fitting over SVM for OSA and control data (SVM; accuracy, 0.77; AUC, 0.84). The RF ML model performs similar to SVM, indicating the comparable statistical fitness to DTI data. The findings indicate that RF model has similar AUC and accuracy over SVM, and either model can be used as a faster OSA screening tool for subjects having brain DTI data.

On variational solutions for whole brain serial-section histology using a Sobolev prior in the computational anatomy random orbit model.

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.

ESTIMATING DIFFEOMORPHIC MAPPINGS BETWEEN TEMPLATES AND NOISY DATA: VARIANCE BOUNDS ON THE ESTIMATED CANONICAL VOLUME FORM.

Anatomy is undergoing a renaissance driven by the availability of large digital data sets generated by light microscopy. A central computational task is to map individual data volumes to standardized templates. This is accomplished by regularized estimation of a diffeomorphic transformation between the coordinate systems of the individual data and the template, building the transformation incrementally by integrating a smooth flow field. The canonical volume form of this transformation is used to quantify local growth, atrophy, or cell density. While multiple implementations exist for this estimation, less attention has been paid to the variance of the estimated diffeomorphism for noisy data. Notably, there is an infinite dimensional unobservable space defined by those diffeomorphisms which leave the template invariant. These form the stabilizer subgroup of the diffeomorphic group acting on the template. The corresponding flat directions in the energy landscape are expected to lead to increased estimation variance. Here we show that a least-action principle used to generate geodesics in the space of diffeomor-phisms connecting the subject brain to the template removes the stabilizer. This provides reduced-variance estimates of the volume form. Using simulations we demonstrate that the asymmetric large deformation diffeomorphic mapping methods (LDDMM), which explicitly incorporate the asymmetry between idealized template images and noisy empirical images, provide lower variance estimators than their symmetrized counterparts (cf. ANTs). We derive Cramer-Rao bounds for the variances in the limit of small deformations. Analytical results are shown for the Jacobian in terms of perturbations of the vector fields and divergence of the vector field.

A Majorization-Minimization Algorithm for Neuroimage Registration.

Intensity-based image registration is critical for neuroimaging tasks, such as 3D reconstruction, times-series alignment, and common coordinate mapping. The gradient-based optimization methods commonly used to solve this problem require a careful selection of step-length. This limitation imposes substantial time and computational costs. Here we propose a gradient-independent rigid-motion registration algorithm based on the majorization-minimization (MM) principle. Each iteration of our intensity-based MM algorithm reduces to a simple point-set rigid registration problem with a closed form solution that avoids the step-length issue altogether. The details of the algorithm are presented, and an error bound for its more practical truncated form is derived. The performance of the MM algorithm is shown to be more effective than gradient descent on simulated images and Nissl stained coronal slices of mouse brain. We also compare and contrast the similarities and differences between the MM algorithm and another gradient-free registration algorithm called the block-matching method. Finally, extensions of this algorithm to more complex problems are discussed.

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.

Applying joint graph embedding to study Alzheimer's neurodegeneration patterns in volumetric data.

Neurodegeneration measured through volumetry in MRI is recognized as a potential Alzheimer's Disease (AD) biomarker, but its utility is limited by lack of specificity. Quantifying spatial patterns of neurodegeneration on a whole brain scale rather than locally may help address this. In this work, we turn to network based analyses and extend a graph embedding algorithm to study morphometric connectivity from volume-change correlations measured with structural MRI on the timescale of years. We model our data with the multiple random eigengraphs framework, as well as modify and implement a multigraph embedding algorithm proposed earlier to estimate a low dimensional embedding of the networks. Our version of the algorithm guarantees meaningful finite-sample results and estimates maximum likelihood edge probabilities from population-specific network modes and subject-specific loadings. Furthermore, we propose and implement a novel statistical testing procedure to analyze group differences after accounting for confounders and locate significant structures during AD neurodegeneration. Family-wise error rate is controlled at 5% using permutation testing on the maximum statistic. We show that results from our analysis reveal networks dominated by known structures associated to AD neurodegeneration, indicating the framework has promise for studying AD. Furthermore, we find network-structure tuples that are not found with traditional methods in the field.

Applying Joint Graph Embedding to Study Alzheimer's Neurodegeneration Patterns in Volumetric Data.

Neurodegeneration measured through volumetry in MRI is recognized as a potential Alzheimer's Disease (AD) biomarker, but its utility is limited by lack of specificity. Quantifying spatial patterns of neurodegeneration on a whole brain scale rather than locally may help address this. In this work, we turn to network based analyses and extend a graph embedding algorithm to study morphometric connectivity from volume-change correlations measured with structural MRI on the timescale of years. We model our data with the multiple random eigengraphs framework, as well as modify and implement a multigraph embedding algorithm proposed earlier to estimate a low dimensional embedding of the networks. Our version of the algorithm guarantees meaningful finite-sample results and estimates maximum likelihood edge probabilities from population-specific network modes and subject-specific loadings. Furthermore, we propose and implement a novel statistical testing procedure to analyze group differences after accounting for confounders and locate significant structures during AD neurodegeneration. Family-wise error rate is controlled at 5% using permutation testing on the maximum statistic. We show that results from our analysis reveal networks dominated by known structures associated to AD neurodegeneration, indicating the framework has promise for studying AD. Furthermore, we find network-structure tuples that are not found with traditional methods in the field.

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.

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.

Alignment of spatial transcriptomics data using diffeomorphic metric mapping.

Spatial transcriptomics (ST) technologies enable high throughput gene expression characterization within thin tissue sections. However, comparing spatial observations across sections, samples, and technologies remains challenging. To address this challenge, we developed STalign to align ST datasets in a manner that accounts for partially matched tissue sections and other local non-linear distortions using diffeomorphic metric mapping. We apply STalign to align ST datasets within and across technologies as well as to align ST datasets to a 3D common coordinate framework. We show that STalign achieves high gene expression and cell-type correspondence across matched spatial locations that is significantly improved over landmark-based affine alignments. Applying STalign to align ST datasets of the mouse brain to the 3D common coordinate framework from the Allen Brain Atlas, we highlight how STalign can be used to lift over brain region annotations and enable the interrogation of compositional heterogeneity across anatomical structures. STalign is available as an open-source Python toolkit at https://github.com/JEFworks-Lab/STalign and as supplementary software with additional documentation and tutorials available at https://jef.works/STalign.

STalign: Alignment of spatial transcriptomics data using diffeomorphic metric mapping.

Spatial transcriptomics (ST) technologies enable high throughput gene expression characterization within thin tissue sections. However, comparing spatial observations across sections, samples, and technologies remains challenging. To address this challenge, we develop STalign to align ST datasets in a manner that accounts for partially matched tissue sections and other local non-linear distortions using diffeomorphic metric mapping. We apply STalign to align ST datasets within and across technologies as well as to align ST datasets to a 3D common coordinate framework. We show that STalign achieves high gene expression and cell-type correspondence across matched spatial locations that is significantly improved over landmark-based affine alignments. Applying STalign to align ST datasets of the mouse brain to the 3D common coordinate framework from the Allen Brain Atlas, we highlight how STalign can be used to lift over brain region annotations and enable the interrogation of compositional heterogeneity across anatomical structures. STalign is available as an open-source Python toolkit at https://github.com/JEFworks-Lab/STalign and as Supplementary Software with additional documentation and tutorials available at https://jef.works/STalign .

Early amygdala and ERC atrophy linked to 3D reconstruction of rostral neurofibrillary tau tangle pathology in Alzheimer's disease.

Previous research has emphasized the unique impact of Alzheimer's Disease (AD) pathology on the medial temporal lobe (MTL), a reflection that tau pathology is particularly striking in the entorhinal and transentorhinal cortex (ERC, TEC) early in the course of disease. However, other brain regions are affected by AD pathology during its early phases. Here, we use longitudinal diffeomorphometry to measure the atrophy rate from MRI of the amygdala compared with that in the ERC and TEC in cognitively unimpaired (CU) controls, CU individuals who progressed to mild cognitive impairment (MCI), and individuals with MCI who progressed to dementia of the AD type (DAT), using a dataset from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Our results show significantly higher atrophy rates of the amygdala in both groups of 'converters' (CU→MCI, MCI→DAT) compared to controls, with rates of volume loss comparable to rates of thickness loss in the ERC and TEC. We localize atrophy within the amygdala within each of these groups using fixed effects modeling. Controlling for the familywise error rate highlights the medial regions of the amygdala as those with significantly higher atrophy in both groups of converters than in controls. Using our recently developed method, referred to as Projective LDDMM, we map measures of neurofibrillary tau tangles (NFTs) from digital pathology to MRI atlases and reconstruct dense 3D spatial distributions of NFT density within regions of the MTL. The distribution of NFTs is consistent with the spatial distribution of MR measured atrophy rates, revealing high densities (and atrophy) in the amygdala (particularly medial), ERC, and rostral third of the MTL. The similarity of the location of NFTs in AD and shape changes in a well-defined clinical population suggests that amygdalar atrophy rate, as measured through MRI may be a viable biomarker for AD.

Disease risk and healthcare utilization among ancestrally diverse groups in the Los Angeles region.

An individual's disease risk is affected by the populations that they belong to, due to shared genetics and environmental factors. The study of fine-scale populations in clinical care is important for identifying and reducing health disparities and for developing personalized interventions. To assess patterns of clinical diagnoses and healthcare utilization by fine-scale populations, we leveraged genetic data and electronic medical records from 35,968 patients as part of the UCLA ATLAS Community Health Initiative. We defined clusters of individuals using identity by descent, a form of genetic relatedness that utilizes shared genomic segments arising due to a common ancestor. In total, we identified 376 clusters, including clusters with patients of Afro-Caribbean, Puerto Rican, Lebanese Christian, Iranian Jewish and Gujarati ancestry. Our analysis uncovered 1,218 significant associations between disease diagnoses and clusters and 124 significant associations with specialty visits. We also examined the distribution of pathogenic alleles and found 189 significant alleles at elevated frequency in particular clusters, including many that are not regularly included in population screening efforts. Overall, this work progresses the understanding of health in understudied communities and can provide the foundation for further study into health inequities.

The cell type composition of the adult mouse brain revealed by single cell and spatial genomics.

The function of the mammalian brain relies upon the specification and spatial positioning of diversely specialized cell types. Yet, the molecular identities of the cell types, and their positions within individual anatomical structures, remain incompletely known. To construct a comprehensive atlas of cell types in each brain structure, we paired high-throughput single-nucleus RNA-seq with Slide-seq-a recently developed spatial transcriptomics method with near-cellular resolution-across the entire mouse brain. Integration of these datasets revealed the cell type composition of each neuroanatomical structure. Cell type diversity was found to be remarkably high in the midbrain, hindbrain, and hypothalamus, with most clusters requiring a combination of at least three discrete gene expression markers to uniquely define them. Using these data, we developed a framework for genetically accessing each cell type, comprehensively characterized neuropeptide and neurotransmitter signaling, elucidated region-specific specializations in activity-regulated gene expression, and ascertained the heritability enrichment of neurological and psychiatric phenotypes. These data, available as an online resource (BrainCellData.org) should find diverse applications across neuroscience, including the construction of new genetic tools, and the prioritization of specific cell types and circuits in the study of brain diseases.

Beyond noise power in 3D computed tomography: the local NPS and off-diagonal elements of the Fourier domain covariance matrix.

PURPOSE: To investigate the correlation and stationarity of noise in volumetric computed tomography (CT) using the local discrete noise-power spectrum (NPS) and off-diagonal elements of the covariance matrix of the discrete Fourier transform of noise-only images (denoted Σ(DFT)). Experimental conditions were varied to affect noise correlation and stationarity, the effects were quantified in terms of the NPS and Σ(DFT), and practical considerations in CT performance characterization were identified. METHODS: Cone-beam CT (CBCT) images were acquired using a benchtop system comprising an x-ray tube and flat-panel detector for a range of acquisition techniques (e.g., dose and x-ray scatter) and three phantom configurations hypothesized to impart distinct effects on the NPS and Σ(DFT): (A) air, (B) a 20-cm-diameter water cylinder with a bowtie filter, and (C) the cylinder without a bowtie filter. The NPS and off-diagonal elements of the Σ(DFT) were analyzed as a function of position within the reconstructions. RESULTS: The local NPS varied systematically throughout the axial plane in a manner consistent with changes in fluence transmitted to the detector and view sampling effects. Variability in fluence was manifest in the NPS magnitude-e.g., a factor of ~2 variation in NPS magnitude within the axial plane for case C (cylinder without bowtie), compared to nearly constant NPS magnitude for case B (bowtie filter matched to the cylinder). View sampling effects were most prominent in case A (air) where the variance increased at greater distance from the center of reconstruction and in case C (cylinder) where the NPS exhibited correlations in the radial direction. The effects of detector lag were observed as azimuthal correlation. The cylinder (without bowtie) had the strongest nonstationarity because of the larger variability in fluence transmitted to the detector. The diagonal elements of the Σ(DFT) were equivalent to the NPS estimated from the periodogram, and the average off-diagonal elements of the Σ(DFT) exhibited amplitude of ~1% of the NPS for the experimental conditions investigated. Furthermore, the off-diagonal elements demonstrated fairly long tails of nearly constant amplitude, with magnitude somewhat reduced for experimental conditions associated with greater stationarity (viz., lower Σ(DFT) tails for cases A and B in comparison to case C). CONCLUSIONS: Volumetric CT exhibits nonstationarity in the NPS as hypothesized in relation to fluence uniformity and view sampling. Measurement of the NPS should seek to minimize such changes in noise correlations and include careful reporting of experimental conditions (e.g., phantom design and use of a bowtie filter) and spatial dependence (e.g., analysis at fixed radius within a phantom). Off-diagonal elements of the Σ(DFT) similarly depend on experimental conditions and can be readily computed from the same data as the NPS. This work begins to check assumptions in NPS analysis examine the extent to which NPS is an appropriate descriptor of noise correlations, and investigate the magnitude of off-diagonal elements of the Σ(DFT). While the magnitude of such off-diagonal elements appears to be low, their cumulative effect on space-variant detectability remains to be investigated-e.g., using task-specific figures of merit.

Effects of protocol and obesity on dose conversion factors in adult body CT.

PURPOSE: In computed tomography (CT), organ dose, effective dose, and risk index can be estimated from volume-weighted CT dose index (CTDI(vol)) or dose-length product (DLP) using conversion coefficients. Studies have investigated how these coefficients vary across scanner models, scan parameters, and patient size. However, their variability across CT protocols has not been systematically studied. Furthermore, earlier studies of the effect of patient size have not included obese individuals, which currently represent more than one-third of U.S. adults. The purpose of this study was to assess the effects of protocol and obesity on dose and risk conversion coefficients in adult body CT. METHODS: Whole-body computational phantoms were created from clinical CT images of six adult patients (three males, three females), representing normal-weight patients and patients of three obesity classes. Body CT protocols at our institution were selected and categorized into ten examination categories based on anatomical region examined. A validated Monte Carlo program was used to estimate organ dose. Organ dose estimates were normalized by CTDI(vol) and size-specific dose estimate (SSDE) to obtain organ dose conversion coefficients (denoted as h and h(ss) factors, respectively). Assuming each phantom to be 20, 40, and 60 years old, effective dose and risk index were calculated and normalized by DLP to obtain effective dose and risk index conversion coefficients (denoted as k and q factors, respectively). Coefficient of variation was used to quantify the variability of each conversion coefficient across examination categories. The effect of obesity was assessed by comparing each obese phantom with the normal-weight phantom of the same gender. RESULTS: For a given organ, the variability of h factor across examination categories that encompassed the entire organ volume was generally within 15%. However, k factor varied more across examination categories (15%-27%). For all three ages, the variability of q factor was small for male (<10%), but large for female phantoms (21%-43%). Relative to the normal-weight phantoms, the reduction in h factor (an average across fully encompassed organs) was 17%-42%, 17%-40%, and 51%-63% for obese-class-I, obese-class-II, and obese-class-III phantoms, respectively. h(ss) factor was not independent of patient diameter and generally decreased with increasing obesity. Relative to the normal-weight phantoms, the reduction in k factor was 12%-40%, 14%-46%, and 44%-59% for obese-class-I, obese-class-II, and obese-class-III phantoms, respectively. The respective reduction in q factor was 11%-36%, 17%-42%, and 48%-59% at 20 years of age and similar at other ages. CONCLUSIONS: In adult body CT, dose to an organ fully encompassed by the primary radiation beam can be estimated from CTDI(vol) using a protocol-independent conversion coefficient. However, fully encompassed organs only account for 50% ± 19% of k factor and 46% ± 24% of q factor. Dose received by partially encompassed organs is also substantial. To estimate effective dose and risk index from DLP, it is necessary to use conversion coefficients specific to the anatomical region examined. Obesity has a significant effect on dose and risk conversion coefficients, which cannot be predicted using body diameter alone. SSDE-normalized organ dose is not independent of diameter. SSDE itself generally overestimates organ dose for obese patients.