Abstract:The brain produces diverse functions, from perceiving sounds to producing arm reaches, through the collective activity of populations of many neurons. Determining if and how the features of these exogenous variables (e.g., sound frequency, reach angle) are reflected in population neural activity is important for understanding how the brain operates. Often, high-dimensional neural population activity is confined to low-dimensional latent spaces. However, many current methods fail to extract latent spaces that are clearly structured by exogenous variables. This has contributed to a debate about whether or not brains should be thought of as dynamical systems or representational systems. Here, we developed a new latent process Bayesian regression framework, the orthogonal stochastic linear mixing model (OSLMM) which introduces an orthogonality constraint amongst time-varying mixture coefficients, and provide Markov chain Monte Carlo inference procedures. We demonstrate superior performance of OSLMM on latent trajectory recovery in synthetic experiments and show superior computational efficiency and prediction performance on several real-world benchmark data sets. We primarily focus on demonstrating the utility of OSLMM in two neural data sets: μ ECoG recordings from rat auditory cortex during presentation of pure tones and multi-single unit recordings form monkey motor cortex during complex arm reaching. We show that OSLMM achieves superior or comparable predictive accuracy of neural data and decoding of external variables (e.g., reach velocity). Most importantly, in both experimental contexts, we demonstrate that OSLMM latent trajectories directly reflect features of the sounds and reaches, demonstrating that neural dynamics are structured by neural representations. Together, these results demonstrate that OSLMM will be useful for the analysis of diverse, large-scale biological time-series datasets. Extracting insight from data into the dynamic processes that produce observed phenomena is a ubiquitous challenge in biology. For example, brain functions are generated by the collective activity of many neurons. Extracting insight from biological time-series data can be challenging because the number of observations can be large, activities can be noisy, and collective dynamics can be complex. Here, we developed a new latent processes Bayesian regression model, the orthogonal stochastic linear mixing model (OSLMM) to address these challenges. Compared to competing methods, we demonstrate that OSLMM has superior recovery performance of latent trajectories in synthetic experiments and is computationally efficient. In two diverse neural data sets, we find that the latent spaces extracted from OSLMM are directly structured by task parameters. Our results demonstrate that latent dynamics can be structured by representations and suggest that OSLMM will be useful for data-driven discovery in large-scale biological time-series data.

A Flexible Bayesian Clustering of Dynamic Subpopulations in Neural Spiking Activity

Bayesian Bi-clustering of Neural Spiking Activity with Latent Structures

Flexible Bayesian Nonparametric Product Mixtures for Multi-scale Functional Clustering

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Dynamic Modeling of Spike Count Data With Conway-Maxwell Poisson Variability

Efficient Inference for Dynamic Flexible Interactions of Neural Populations

Clustering action potential spikes: Insights on the use of overfitted finite mixture models and Dirichlet process mixture models

Bayesian nonparametric analysis for the detection of spikes in noisy calcium imaging data

Bayesian inference of structured latent spaces from neural population activity with the orthogonal stochastic linear mixing model

Bayesian model-based clustering for populations of network data

Conditional Finite Mixtures of Poisson Distributions for Context-Dependent Neural Correlations

Spatiotemporal Clustering with Neyman-Scott Processes via Connections to Bayesian Nonparametric Mixture Models

Bayesian estimation of clustered dependence structures in functional neuroconnectivity

A functional clustering algorithm for the analysis of dynamic network data

An unbiased method to partition diverse neuronal responses into functional ensembles reveals interpretable population dynamics during innate social behavior

A novel method for sparse dynamic functional connectivity analysis from resting-state fMRI

Sparse Bayesian Hierarchical Modeling of High-dimensional Clustering Problems

Clustering of heterogeneous populations of networks

A Coarse-Grained Framework for Spiking Neuronal Networks: Between Homogeneity and Synchrony

A sticky Poisson Hidden Markov Model for spike data

Covariate-guided Bayesian mixture model for multivariate time series