Matrix Inversion and Subset Selection (MISS): A pipeline for mapping of diverse cell types across the murine brain.

The advent of increasingly sophisticated imaging platforms has allowed for the visualization of the murine nervous system at single-cell resolution. However, current experimental approaches have not yet produced whole-brain maps of a comprehensive set of neuronal and nonneuronal types that approaches the cellular diversity of the mammalian cortex. Here, we aim to fill in this gap in knowledge with an open-source computational pipeline, Matrix Inversion and Subset Selection (MISS), that can infer quantitatively validated distributions of diverse collections of neural cell types at 200-µm resolution using a combination of single-cell RNA sequencing (RNAseq) and in situ hybridization datasets. We rigorously demonstrate the accuracy of MISS against literature expectations. Importantly, we show that gene subset selection, a procedure by which we filter out low-information genes prior to performing deconvolution, is a critical preprocessing step that distinguishes MISS from its predecessors and facilitates the production of cell-type maps with significantly higher accuracy. We also show that MISS is generalizable by generating high-quality cell-type maps from a second independently curated single-cell RNAseq dataset. Together, our results illustrate the viability of computational approaches for determining the spatial distributions of a wide variety of cell types from genetic data alone.

Emergence of directional bias in tau deposition from axonal transport dynamics.

Defects in axonal transport may partly underpin the differences between the observed pathophysiology of Alzheimer's disease (AD) and that of other non-amyloidogenic tauopathies. Particularly, pathological tau variants may have molecular properties that dysregulate motor proteins responsible for the anterograde-directed transport of tau in a disease-specific fashion. Here we develop the first computational model of tau-modified axonal transport that produces directional biases in the spread of tau pathology. We simulated the spatiotemporal profiles of soluble and insoluble tau species in a multicompartment, two-neuron system using biologically plausible parameters and time scales. Changes in the balance of tau transport feedback parameters can elicit anterograde and retrograde biases in the distributions of soluble and insoluble tau between compartments in the system. Aggregation and fragmentation parameters can also perturb this balance, suggesting a complex interplay between these distinct molecular processes. Critically, we show that the model faithfully recreates the characteristic network spread biases in both AD-like and non-AD-like mouse tauopathy models. Tau transport feedback may therefore help link microscopic differences in tau conformational states and the resulting variety in clinical presentations.

Neural connectivity predicts spreading of alpha-synuclein pathology in fibril-injected mouse models: Involvement of retrograde and anterograde axonal propagation.

In Parkinson's disease, some of the first alpha-synuclein aggregates appear in the olfactory system and the dorsal motor nucleus of the vagus nerve before spreading to connected brain regions. We previously demonstrated that injection of alpha-synuclein fibrils unilaterally into the olfactory bulb of wild type mice leads to widespread synucleinopathy in brain regions directly and indirectly connected to the injection site, consistently, over the course of periods longer than 6 months. Our previously reported observations support the idea that alpha-synuclein inclusions propagates between brain region through neuronal networks. In the present study, we further defined the pattern of propagation of alpha-synuclein inclusions and developed a mathematical model based on known mouse brain connectivity. Using this model, we first predicted the pattern of alpha-synuclein inclusions propagation following an injection of fibrils into the olfactory bulb. We then analyzed the fitting of these predictions to our published histological data. Our results demonstrate that the pattern of propagation we observed in vivo is consistent with axonal transport of alpha-synuclein aggregate seeds, followed by transsynaptic transmission. By contrast, simple diffusion of alpha-synuclein fits very poorly our in vivo data. We also found that the spread of alpha-synuclein inclusions appeared to primarily follow neural connections retrogradely until 9 months after injection into the olfactory bulb. Thereafter, the pattern of spreading was consistent with anterograde propagation mathematical models. Finally, we applied our mathematical model to a different, previously published, dataset involving alpha-synuclein fibril injections into the striatum, instead of the olfactory bulb. We found that the mathematical model accurately predicts the reported progressive increase in alpha-synuclein neuropathology also in that paradigm. In conclusion, our findings support that the progressive spread of alpha-synuclein inclusions after injection of protein fibrils follows neural networks in the mouse connectome.

Directionality bias underpins divergent spatiotemporal progression of Alzheimer-related tauopathy in mouse models.

Mounting evidence implicates trans-synaptic connectome-based spread as a shared mechanism behind different tauopathic conditions, yet also suggests there is divergent spatiotemporal progression between them. A potential parsimonious explanation for this apparent contradiction could be that different conditions incur differential rates and directional biases in tau transmission along fiber tracts. In this meta-analysis we closely examined this hypothesis and quantitatively tested it using spatiotemporal tau pathology patterns from 11 distinct models across 4 experimental studies. For this purpose, we extended a network-based spread model by incorporating net directionality along the connectome. Our data unambiguously supports the directional transmission hypothesis. First, retrograde bias is an unambiguously better predictor of tau progression than anterograde bias. Second, while spread exhibits retrograde character, our best-fitting biophysical models incorporate the mixed effects of both retrograde- and anterograde-directed spread, with notable tau-strain-specific differences. We also found a nontrivial association between directionality bias and tau strain aggressiveness, with more virulent strains exhibiting less retrograde character. Taken together, our study implicates directional transmission bias in tau transmission along fiber tracts as a general feature of tauopathy spread and a strong candidate explanation for the diversity of spatiotemporal tau progression between conditions. This simple and parsimonious mechanism may potentially fill a critical gap in our knowledge of the spatiotemporal ramification of divergent tauopathies.

Establishing neuroanatomical correspondences across mouse and marmoset brain structures.

Interest in the common marmoset is growing due to evolutionarily proximity to humans compared to laboratory mice, necessitating a comparison of mouse and marmoset brain architectures, including connectivity and cell type distributions. Creating an actionable comparative platform is challenging since these brains have distinct spatial organizations and expert neuroanatomists disagree. We propose a general theoretical framework to relate named atlas compartments across taxa and use it to establish a detailed correspondence between marmoset and mice brains. Contrary to conventional wisdom that brain structures may be easier to relate at higher levels of the atlas hierarchy, we find that finer parcellations at the leaf levels offer greater reconcilability despite naming discrepancies. Utilizing existing atlases and associated literature, we created a list of leaf- level structures for both species and establish five types of correspondence between them. One-to-one relations were found between 43% of the structures in mouse and 47% in marmoset, whereas 25% of mouse and 10% of marmoset structures were not relatable. The remaining structures show a set of more complex mappings which we quantify. Implementing this correspondence with volumetric atlases of the two species, we make available a computational tool for querying and visualizing relationships between the corresponding brains. Our findings provide a foundation for computational comparative analyses of mesoscale connectivity and cell type distributions in the laboratory mouse and the common marmoset.

Establishing neuroanatomical correspondences across mouse and marmoset brain structures.

Interest in the common marmoset is growing due to evolutionarily proximity to humans compared to laboratory mice, necessitating a comparison of mouse and marmoset brain architectures, including connectivity and cell type distributions. Creating an actionable comparative platform is challenging since these brains have distinct spatial organizations and expert neuroanatomists disagree. We propose a general theoretical framework to relate named atlas compartments across taxa and use it to establish a detailed correspondence between marmoset and mice brains. Contrary to conventional wisdom that brain structures may be easier to relate at higher levels of the atlas hierarchy, we find that finer parcellations at the leaf levels offer greater reconcilability despite naming discrepancies. Utilizing existing atlases and associated literature, we created a list of leaf-level structures for both species and establish five types of correspondence between them. One-to-one relations were found between 43% of the structures in mouse and 47% in marmoset, whereas 25% of mouse and 10% of marmoset structures were not relatable. The remaining structures show a set of more complex mappings which we quantify. Implementing this correspondence with volumetric atlases of the two species, we make available a computational tool for querying and visualizing relationships between the corresponding brains. Our findings provide a foundation for computational comparative analyses of mesoscale connectivity and cell type distributions in the laboratory mouse and the common marmoset.

Spatial cell-type enrichment predicts mouse brain connectivity.

A fundamental neuroscience topic is the link between the brain's molecular, cellular, and cytoarchitectonic properties and structural connectivity. Recent studies relate inter-regional connectivity to gene expression, but the relationship to regional cell-type distributions remains understudied. Here, we utilize whole-brain mapping of neuronal and non-neuronal subtypes via the matrix inversion and subset selection algorithm to model inter-regional connectivity as a function of regional cell-type composition with machine learning. We deployed random forest algorithms for predicting connectivity from cell-type densities, demonstrating surprisingly strong prediction accuracy of cell types in general, and particular non-neuronal cells such as oligodendrocytes. We found evidence of a strong distance dependency in the cell connectivity relationship, with layer-specific excitatory neurons contributing the most for long-range connectivity, while vascular and astroglia were salient for short-range connections. Our results demonstrate a link between cell types and connectivity, providing a roadmap for examining this relationship in other species, including humans.

The effects of microglia on tauopathy progression can be quantified using Nexopathy in silico (Nexis) models.

The prion-like transsynaptic propagation of misfolded tau along the brain's connectome has previously been modeled using connectome-based network diffusion models. In addition to the connectome, interactions between the general neurological "milieu" in the neurodegenerative brain and proteinopathic species can also contribute to pathology propagation. Such a molecular nexopathy framework posits that the distinct characteristics of neurodegenerative disorders stem from interactions between the network and surrounding molecular players. However, the effects of these modulators remain unquantified. Here, we present Nexopathy in silico ("Nexis"), a quantitative model of tau progression augmenting earlier models by including parameters of pathology propagation defined by the molecular modulators of connectome-based spread. Our Nexis:microglia model provides the first quantitative characterization of this effect on the whole brain by expanding previous models of neuropathology progression by incorporating microglial influence. We show that Trem2, but not microglial homeostasis genes, significantly improved the model's predictive power. Trem2 appears to reduce tau accumulation rate while increasing its interregional spread from the hippocampal seed area, causing higher tau burden in the striatum, pallidum, and contralateral hippocampus. Nexis provides an improved understanding and quantification of microglial contribution to tau propagation and can be flexibly modified to include other modulators of progressive neurodegeneration.