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.

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.

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.

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/ .

Hidden Markov modeling for maximum probability neuron reconstruction.

Recent advances in brain clearing and imaging have made it possible to image entire mammalian brains at sub-micron resolution. These images offer the potential to assemble brain-wide atlases of neuron morphology, but manual neuron reconstruction remains a bottleneck. Several automatic reconstruction algorithms exist, but most focus on single neuron images. In this paper, we present a probabilistic reconstruction method, ViterBrain, which combines a hidden Markov state process that encodes neuron geometry with a random field appearance model of neuron fluorescence. ViterBrain utilizes dynamic programming to compute the global maximizer of what we call the most probable neuron path. We applied our algorithm to imperfect image segmentations, and showed that it can follow axons in the presence of noise or nearby neurons. We also provide an interactive framework where users can trace neurons by fixing start and endpoints. ViterBrain is available in our open-source Python package brainlit.

Is Neuroscience FAIR? A Call for Collaborative Standardisation of Neuroscience Data.

In this perspective article, we consider the critical issue of data and other research object standardisation and, specifically, how international collaboration, and organizations such as the International Neuroinformatics Coordinating Facility (INCF) can encourage that emerging neuroscience data be Findable, Accessible, Interoperable, and Reusable (FAIR). As neuroscientists engaged in the sharing and integration of multi-modal and multiscale data, we see the current insufficiency of standards as a major impediment in the Interoperability and Reusability of research results. We call for increased international collaborative standardisation of neuroscience data to foster integration and efficient reuse of research objects.

Fitting Splines to Axonal Arbors Quantifies Relationship Between Branch Order and Geometry.

Neuromorphology is crucial to identifying neuronal subtypes and understanding learning. It is also implicated in neurological disease. However, standard morphological analysis focuses on macroscopic features such as branching frequency and connectivity between regions, and often neglects the internal geometry of neurons. In this work, we treat neuron trace points as a sampling of differentiable curves and fit them with a set of branching B-splines. We designed our representation with the Frenet-Serret formulas from differential geometry in mind. The Frenet-Serret formulas completely characterize smooth curves, and involve two parameters, curvature and torsion. Our representation makes it possible to compute these parameters from neuron traces in closed form. These parameters are defined continuously along the curve, in contrast to other parameters like tortuosity which depend on start and end points. We applied our method to a dataset of cortical projection neurons traced in two mouse brains, and found that the parameters are distributed differently between primary, collateral, and terminal axon branches, thus quantifying geometric differences between different components of an axonal arbor. The results agreed in both brains, further validating our representation. The code used in this work can be readily applied to neuron traces in SWC format and is available in our open-source Python package brainlit: http://brainlit.neurodata.io/.

A 7 Tesla Amygdalar-Hippocampal Shape Analysis of Lithium Response in Bipolar Disorder.

Research to discover clinically useful predictors of lithium response in patients with bipolar disorder has largely found them to be elusive. We demonstrate here that detailed neuroimaging may have the potential to fill this important gap in mood disorder therapeutics. Lithium treatment and bipolar disorder have both been shown to affect anatomy of the hippocampi and amygdalae but there is no consensus on the nature of their effects. We aimed to investigate structural surface anatomy changes in amygdala and hippocampus correlated with treatment response in bipolar disorder. Patients with bipolar disorder (N = 14) underwent lithium treatment, were classified by response status at acute and long-term time points, and scanned with 7 Tesla structural MRI. Large Deformation Diffeomorphic Metric Mapping was applied to detect local differences in hippocampal and amygdalar anatomy between lithium responders and non-responders. Anatomy was also compared to 21 healthy comparison participants. A patch of the ventral surface of the left hippocampus was found to be significantly atrophied in non-responders as compared to responders at the acute time point and was associated at a trend-level with long-term response status. We did not detect an association between response status and surface anatomy of the right hippocampus or amygdala. To the best of our knowledge, this is the first shape analysis of hippocampus and amygdala in bipolar disorder using 7 Tesla MRI. These results can inform future work investigating possible neuroimaging predictors of lithium response in bipolar disorder.

Multiscale Coupling of an Agent-Based Model of Tissue Fibrosis and a Logic-Based Model of Intracellular Signaling.

Wound healing and fibrosis following myocardial infarction (MI) is a dynamic process involving many cell types, extracellular matrix (ECM), and inflammatory cues. As both incidence and survival rates for MI increase, management of post-MI recovery and associated complications are an increasingly important focus. Complexity of the wound healing process and the need for improved therapeutics necessitate a better understanding of the biochemical cues that drive fibrosis. To study the progression of cardiac fibrosis across spatial and temporal scales, we developed a novel hybrid multiscale model that couples a logic-based differential equation (LDE) model of the fibroblast intracellular signaling network with an agent-based model (ABM) of multi-cellular tissue remodeling. The ABM computes information about cytokine and growth factor levels in the environment including TGFß, TNFα, IL-1ß, and IL-6, which are passed as inputs to the LDE model. The LDE model then computes the network signaling state of individual cardiac fibroblasts within the ABM. Based on the current network state, fibroblasts make decisions regarding cytokine secretion and deposition and degradation of collagen. Simulated fibroblasts respond dynamically to rapidly changing extracellular environments and contribute to spatial heterogeneity in model predicted fibrosis, which is governed by many parameters including cell density, cell migration speeds, and cytokine levels. Verification tests confirmed that predictions of the coupled model and network model alone were consistent in response to constant cytokine inputs and furthermore, a subset of coupled model predictions were validated with in vitro experiments with human cardiac fibroblasts. This multiscale framework for cardiac fibrosis will allow for systematic screening of the effects of molecular perturbations in fibroblast signaling on tissue-scale extracellular matrix composition and organization.