Models and simulation of 3D neuronal dendritic trees using Bayesian networks

NeuroSuites: An online platform for running neuroscience, statistical, and machine learning tools

Nowadays, an enormous amount of high dimensional data is available in the field of neuroscience. Handling these data is complex and requires the use of efficient tools to transform them into useful knowledge. In this work we present NeuroSuites, an easy-access web platform with its own architecture. We compare our platform with other software currently available, highlighting its main strengths. Thanks to its defined architecture, it is able to handle large-scale problems common in some neuroscience fields. NeuroSuites has different neuroscience-oriented applications and tools to integrate statistical data analysis and machine learning algorithms commonly used in this field. As future work, we want to further expand the list of available software tools as well as improve the platform interface according to user demands.

Computational synthesis of cortical dendritic morphologies

Neuronal morphologies provide the foundation for the electrical behavior of neurons, the connectomes they form, and the dynamical properties of the brain. Comprehensive neuron models are essential for defining cell types, discerning their functional roles, and investigating brain-disease-related dendritic alterations. However, a lack of understanding of the principles underlying neuron morphologies has hindered attempts to computationally synthesize morphologies for decades. We introduce a synthesis algorithm based on a topological descriptor of neurons, which enables the rapid digital reconstruction of entire brain regions from few reference cells. This topology-guided synthesis generates dendrites that are statistically similar to biological reconstructions in terms of morpho-electrical and connectivity properties and offers a significant opportunity to investigate the links between neuronal morphology and brain function across different spatiotemporal scales. Synthesized cortical networks based on structurally altered dendrites associated with diverse brain pathologies revealed principles linking branching properties to the structure of large-scale networks.

Efficient metadata mining of web-accessible neural morphologies

Advancements in neuroscience research have led to steadily accelerating data production and sharing. The online community repository of neural reconstructions NeuroMorpho.Org grew from fewer than 1000 digitally traced neurons in 2006 to more than 140,000 cells today, including glia that now constitute 10.1% of the content. Every reconstruction consists of a detailed 3D representation of branch geometry and connectivity in a standardized format, from which a collection of morphometric features is extracted and stored. Moreover, each entry in the database is accompanied by rich metadata annotation describing the animal subject, anatomy, and experimental details. The rapid expansion of this resource in the past decade was accompanied by a parallel rise in the complexity of the available information, creating both opportunities and challenges for knowledge mining. Here, we introduce a new summary reporting functionality, allowing NeuroMorpho.Org users to efficiently download digests of metadata and morphometrics from multiple groups of similar cells for further analysis. We demonstrate the capabilities of the tool for both glia and neurons and present an illustrative statistical analysis of the resulting data.

Comparing the Electrophysiology and Morphology of Human and Mouse Layer 2/3 Pyramidal Neurons With Bayesian Networks

Pyramidal neurons are the most common neurons in the cerebral cortex. Understanding how they differ between species is a key challenge in neuroscience. We compared human temporal cortex and mouse visual cortex pyramidal neurons from the Allen Cell Types Database in terms of their electrophysiology and dendritic morphology. We found that, among other differences, human pyramidal neurons had a higher action potential threshold voltage, a lower input resistance, and larger dendritic arbors. We learned Gaussian Bayesian networks from the data in order to identify correlations and conditional independencies between the variables and compare them between the species. We found strong correlations between electrophysiological and morphological variables in both species. In human cells, electrophysiological variables were correlated even with morphological variables that are not directly related to dendritic arbor size or diameter, such as mean bifurcation angle and mean branch tortuosity. Cortical depth was correlated with both electrophysiological and morphological variables in both species, and its effect on electrophysiology could not be explained in terms of the morphological variables. For some variables, the effect of cortical depth was opposite in the two species. Overall, the correlations among the variables differed strikingly between human and mouse neurons. Besides identifying correlations and conditional independencies, the learned Bayesian networks might be useful for probabilistic reasoning regarding the morphology and electrophysiology of pyramidal neurons.

Comparing basal dendrite branches in human and mouse hippocampal CA1 pyramidal neurons with Bayesian networks

Pyramidal neurons are the most common cell type in the cerebral cortex. Understanding how they differ between species is a key challenge in neuroscience. A recent study provided a unique set of human and mouse pyramidal neurons of the CA1 region of the hippocampus, and used it to compare the morphology of apical and basal dendritic branches of the two species. The study found inter-species differences in the magnitude of the morphometrics and similarities regarding their variation with respect to morphological determinants such as branch type and branch order. We use the same data set to perform additional comparisons of basal dendrites. In order to isolate the heterogeneity due to intrinsic differences between species from the heterogeneity due to differences in morphological determinants, we fit multivariate models over the morphometrics and the determinants. In particular, we use conditional linear Gaussian Bayesian networks, which provide a concise graphical representation of the independencies and correlations among the variables. We also extend the previous study by considering additional morphometrics and by formally testing whether a morphometric increases or decreases with the distance from the soma. This study introduces a multivariate methodology for inter-species comparison of morphology.

NeuroPath2Path: Classification and elastic morphing between neuronal arbors using path-wise similarity

Neuron shape and connectivity affect function. Modern imaging methods have proven successful at extracting morphological information. One potential path to achieve analysis of this morphology is through graph theory. Encoding by graphs enables the use of high throughput informatic methods to extract and infer brain function. However, the application of graph-theoretic methods to neuronal morphology comes with certain challenges in term of complex subgraph matching and the difficulty in computing intermediate shapes in between two imaged temporal samples. Here we report a novel, efficacious graph-theoretic method that rises to the challenges. The morphology of a neuron, which consists of its overall size, global shape, local branch patterns, and cell-specific biophysical properties, can vary significantly with the cell's identity, location, as well as developmental and physiological state. Various algorithms have been developed to customize shape based statistical and graph related features for quantitative analysis of neuromorphology, followed by the classification of neuron cell types using the features. Unlike the classical feature extraction based methods from imaged or 3D reconstructed neurons, we propose a model based on the rooted path decomposition from the soma to the dendrites of a neuron and extract morphological features from each constituent path. We hypothesize that measuring the distance between two neurons can be realized by minimizing the cost of continuously morphing the set of all rooted paths of one neuron to another. To validate this claim, we first establish the correspondence of paths between two neurons using a modified Munkres algorithm. Next, an elastic deformation framework that employs the square root velocity function is established to perform the continuous morphing, which, as an added benefit, provides an effective visualization tool. We experimentally show the efficacy of NeuroPath2Path, NeuroP2P, over the state of the art.

An open repository for single-cell reconstructions of the brain forest

NeuroMorpho.Org was launched in 2006 to provide unhindered access to any and all digital tracings of neuronal morphology that researchers were willing to share freely upon request. Today this database is the largest public inventory of cellular reconstructions in neuroscience with a content of over 80,000 neurons and glia from a representative diversity of animal species, anatomical regions, and experimental methods. Datasets continuously contributed by hundreds of laboratories worldwide are centrally curated, converted into a common non-proprietary format, morphometrically quantified, and annotated with comprehensive metadata. Users download digital reconstructions for a variety of scientific applications including visualization, classification, analysis, and simulations. With more than 1,000 peer-reviewed publications describing data stored in or utilizing data retrieved from NeuroMorpho.Org, this ever-growing repository can already be considered a mature resource for neuroscience.

Dendritic branching angles of pyramidal cells across layers of the juvenile rat somatosensory cortex

The characterization of the structural design of cortical microcircuits is essential for understanding how they contribute to function in both health and disease. Since pyramidal neurons represent the most abundant neuronal type and their dendritic spines constitute the major postsynaptic elements of cortical excitatory synapses, our understanding of the synaptic organization of the neocortex largely depends on the available knowledge regarding the structure of pyramidal cells. Previous studies have identified several apparently common rules in dendritic geometry. We study the dendritic branching angles of pyramidal cells across layers to further shed light on the principles that determine the geometric shapes of these cells. We find that the dendritic branching angles of pyramidal cells from layers II-VI of the juvenile rat somatosensory cortex suggest common design principles, despite the particular morphological and functional features that are characteristic of pyramidal cells in each cortical layer. J. Comp. Neurol. 524:2567-2576, 2016. © 2016 Wiley Periodicals, Inc.

Tree Topology Estimation

Tree-like structures are fundamental in nature, and it is often useful to reconstruct the topology of a tree - what connects to what - from a two-dimensional image of it. However, the projected branches often cross in the image: the tree projects to a planar graph, and the inverse problem of reconstructing the topology of the tree from that of the graph is ill-posed. We regularize this problem with a generative, parametric tree-growth model. Under this model, reconstruction is possible in linear time if one knows the direction of each edge in the graph - which edge endpoint is closer to the root of the tree - but becomes NP-hard if the directions are not known. For the latter case, we present a heuristic search algorithm to estimate the most likely topology of a rooted, three-dimensional tree from a single two-dimensional image. Experimental results on retinal vessel, plant root, and synthetic tree data sets show that our methodology is both accurate and efficient.

Bayesian networks in neuroscience: a survey

Bayesian networks are a type of probabilistic graphical models lie at the intersection between statistics and machine learning. They have been shown to be powerful tools to encode dependence relationships among the variables of a domain under uncertainty. Thanks to their generality, Bayesian networks can accommodate continuous and discrete variables, as well as temporal processes. In this paper we review Bayesian networks and how they can be learned automatically from data by means of structure learning algorithms. Also, we examine how a user can take advantage of these networks for reasoning by exact or approximate inference algorithms that propagate the given evidence through the graphical structure. Despite their applicability in many fields, they have been little used in neuroscience, where they have focused on specific problems, like functional connectivity analysis from neuroimaging data. Here we survey key research in neuroscience where Bayesian networks have been used with different aims: discover associations between variables, perform probabilistic reasoning over the model, and classify new observations with and without supervision. The networks are learned from data of any kind-morphological, electrophysiological, -omics and neuroimaging-, thereby broadening the scope-molecular, cellular, structural, functional, cognitive and medical- of the brain aspects to be studied.

Branching angles of pyramidal cell dendrites follow common geometrical design principles in different cortical areas

Unraveling pyramidal cell structure is crucial to understanding cortical circuit computations. Although it is well known that pyramidal cell branching structure differs in the various cortical areas, the principles that determine the geometric shapes of these cells are not fully understood. Here we analyzed and modeled with a von Mises distribution the branching angles in 3D reconstructed basal dendritic arbors of hundreds of intracellularly injected cortical pyramidal cells in seven different cortical regions of the frontal, parietal, and occipital cortex of the mouse. We found that, despite the differences in the structure of the pyramidal cells in these distinct functional and cytoarchitectonic cortical areas, there are common design principles that govern the geometry of dendritic branching angles of pyramidal cells in all cortical areas.

Neuronal morphology goes digital: a research hub for cellular and system neuroscience

The importance of neuronal morphology in brain function has been recognized for over a century. The broad applicability of "digital reconstructions" of neuron morphology across neuroscience subdisciplines has stimulated the rapid development of numerous synergistic tools for data acquisition, anatomical analysis, three-dimensional rendering, electrophysiological simulation, growth models, and data sharing. Here we discuss the processes of histological labeling, microscopic imaging, and semiautomated tracing. Moreover, we provide an annotated compilation of currently available resources in this rich research "ecosystem" as a central reference for experimental and computational neuroscience.

New insights into the classification and nomenclature of cortical GABAergic interneurons

A systematic classification and accepted nomenclature of neuron types is much needed but is currently lacking. This article describes a possible taxonomical solution for classifying GABAergic interneurons of the cerebral cortex based on a novel, web-based interactive system that allows experts to classify neurons with pre-determined criteria. Using Bayesian analysis and clustering algorithms on the resulting data, we investigated the suitability of several anatomical terms and neuron names for cortical GABAergic interneurons. Moreover, we show that supervised classification models could automatically categorize interneurons in agreement with experts' assignments. These results demonstrate a practical and objective approach to the naming, characterization and classification of neurons based on community consensus.

Improved estimates of visual field progression using bayesian linear regression to integrate structural information in patients with ocular hypertension

PURPOSE

To assess whether neuroretinal rim area (RA) measurements of the optic disc could be used to improve the estimate of the rate of change in visual field (VF) mean sensitivity in patients with ocular hypertension (OHT) using a Bayesian linear regression (BLR), compared to a standard ordinary least squares linear regression (OLSLR) of mean sensitivity (MS) measurements alone.

METHODS

MS and RA measurements were analyzed from a longitudinal series of 179 patients with OHT visiting Moorfields Eye Hospital between 1992 and 2000. For each patient, linear regression of RA was computed after an appropriate transformation to "scale" RA with MS measurements, and the slope coefficient from this regression was used as a prior for BLR of MS. The BLR then was compared with the OLSLR approach by evaluating how accurately each regression technique predicted future MS measurements.

RESULTS

On average, BLR was significantly more accurate than OLSLR for series up to 8 measurements long (root-mean-square prediction error [RMSPE] was 0.14 decibels [dB] smaller with BLR than OLSLR; P < 0.001, Wilcoxon signed-rank test), with OLSLR of VF data alone being more accurate for longer series (RMSPE was 0.06 dB smaller with OLSLR than BLR).

CONCLUSIONS

BLR provides a significantly more accurate estimate of the rate of change in MS than the standard OLSLR approach, especially in short time series, suggesting that structural measurements can be used successfully in statistical models to assist clinicians monitoring VF progression in patients with OHT. Further studies are necessary to validate the method in glaucoma patients.