Abstract:Neuroscientific data analysis has classically involved methods for statistical signal and image processing, drawing on linear algebra and stochastic process theory. However, digitized neuroanatomical data sets containing labelled neurons, either individually or in groups labelled by tracer injections, do not fully fit into this classical framework. The tree-like shapes of neurons cannot mathematically be adequately described as points in a vector space. There is therefore a need for new approaches. Methods from computational topology and geometry are naturally suited to the analysis of neuronal shapes. Here we introduce methods from Discrete Morse Theory to extract tree-skeletons of individual neurons from volumetric brain image data, or to summarize collections of neurons labelled by localized anterograde tracer injections. Since individual neurons are topologically trees, it is sensible to summarize the collection of neurons labelled by a localized anterograde tracer injection using a consensus tree-shape. The algorithmic procedure includes an initial pre-processing step to extract a density field from the raw volumetric image data, followed by initial skeleton extraction from the density field using a discrete version of a 1-(un)stable manifold of the density field. Heuristically, if the density field is regarded as a mountainous landscape, then the 1-(un)stable manifold follows the "mountain ridges" connecting the maxima of the density field. We then simplify this skeleton-graph into a tree using a shortest-path approach and methods derived from persistent homology. The advantage of this approach is that it uses global information about the density field and is therefore robust to local fluctuations and non-uniformly distributed input signals. To be able to handle large data sets, we use a divide-and-conquer approach. The resulting software DiMorSC is available on Github.

Interpretation of autoencoder-learned collective variables using Morse–Smale complex and sublevelset persistent homology: An application on molecular trajectories

Analyzing multimodal probability measures with autoencoders

Structure-preserving visualization for single-cell RNA-Seq profiles using deep manifold transformation with batch-correction

Understanding recent deep-learning techniques for identifying collective variables of molecular dynamics

Autoencoders for dimensionality reduction in molecular dynamics: collective variable dimension, biasing and transition states

Nonlinear discovery of slow molecular modes using state-free reversible VAMPnets

Topological Skeletonization and Tree-Summarization of Neurons Using Discrete Morse Theory

Descriptors-free Collective Variables From Geometric Graph Neural Networks

Automated Global Analysis of Experimental Dynamics through Low-Dimensional Linear Embeddings

Effective Many-body Interactions in Reduced-Dimensionality Spaces Through Neural Network Models

Toward a structural identification of metastable molecular conformations

Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics

Improving Self-supervised Molecular Representation Learning using Persistent Homology

MLCV: Bridging Machine-Learning-Based Dimensionality Reduction and Free-Energy Calculation

Memory Kernel Minimization Based Neural Networks for Discovering Slow Collective Variables of Biomolecular Dynamics

Autoencoders for discovering manifold dimension and coordinates in data from complex dynamical systems

Analysis of the Geometrical Evolution in On-the-Fly Surface-Hopping Nonadiabatic Dynamics with Machine Learning Dimensionality Reduction Approaches: Classical Multidimensional Scaling and Isometric Feature Mapping.

Generalized Morse Theory of Distance Functions to Surfaces for Persistent Homology

Deep multi-modal learning for joint linear representation of nonlinear dynamical systems

Characterizing the Latent Space of Molecular Deep Generative Models with Persistent Homology Metrics

Molecular CT: Unifying Geometry and Representation Learning for Molecules at Different Scales