A Nonparametric Bayesian Approach to Uncovering Rat Hippocampal Population Codes During Spatial Navigation

Rodent hippocampal population codes represent important spatial information about the environment during navigation. Several computational methods have been developed to uncover the neural representation of spatial topology embedded in rodent hippocampal ensemble spike activity. Here we extend our previous work and propose a nonparametric Bayesian approach to infer rat hippocampal population codes during spatial navigation. To tackle the model selection problem, we leverage a nonparametric Bayesian model. Specifically, to analyze rat hippocampal ensemble spiking activity, we apply a hierarchical Dirichlet process-hidden Markov model (HDP-HMM) using two Bayesian inference methods, one based on Markov chain Monte Carlo (MCMC) and the other based on variational Bayes (VB). We demonstrate the effectiveness of our Bayesian approaches on recordings from a freely-behaving rat navigating in an open field environment. We find that MCMC-based inference with Hamiltonian Monte Carlo (HMC) hyperparameter sampling is flexible and efficient, and outperforms VB and MCMC approaches with hyperparameters set by empirical Bayes

Advances in point process modeling: feature selection, goodness-of-fit and novel applications

The research contained in this thesis extends multivariate marked point process modeling methods for neuroscience, generalizes goodness-of-fit techniques for the class of marked point processes, and introduces the use of a general history-dependent point process model to the domain of sleep apnea. Our first project involves further development of a modeling tool for spiking data from neural populations using the theory of marked point processes. This marked point process model uses features of spike waveforms as marks in order to estimate a state variable of interest. We examine the informational content of geometric features as well as principal components of the waveforms at hippocampal place cell activity by comparing decoding accuracies of a rat's position along a track. We determined that there was additional information available beyond that contained in traditional geometric features used for decoding in practice. The expanded use of this marked point process model in neuroscience necessitates corresponding goodness-of-fit protocols for the marked case. In our second project, we develop a generalized time-rescaling method for marked point processes that produces uniformly distributed spikes under a proper model. Once rescaled, the ground process then behaves as a Poisson process and can be analyzed using traditional point process goodness-of-fit methods. We demonstrate the method's ability to detect quality and manner of fit through both simulation and real neural data analysis. In the final project, we introduce history-dependent point process modeling as a superior method for characterizing severe sleep apnea over the current clinical standard known as the apnea-hypopnea index (AHI). We analyze model fits using combinations of both clinical covariates and event observations themselves through functions of history. Ultimately, apnea onset times were consistently estimated with significantly higher accuracy when history was incorporated alongside sleep stage. We present this method to the clinical audience as a means to gain detailed information on patterns of apnea and to provide more customized diagnoses and treatment prescriptions. These separate yet complementary projects extend existing point process modeling methods and further demonstrate their value in the neurosciences, sleep sciences, and beyond

IST Austria Thesis

The hippocampus is a key brain region for spatial memory and navigation and is needed at all stages of memory, including encoding, consolidation, and recall. Hippocampal place cells selectively discharge at specific locations of the environment to form a cognitive map of the space. During the rest period and sleep following spatial navigation and/or learning, the waking activity of the place cells is reactivated within high synchrony events. This reactivation is thought to be important for memory consolidation and stabilization of the spatial representations. The aim of my thesis was to directly test whether the reactivation content encoded in firing patterns of place cells is important for consolidation of spatial memories. In particular, I aimed to test whether, in cases when multiple spatial memory traces are acquired during learning, the specific disruption of the reactivation of a subset of these memories leads to the selective disruption of the corresponding memory traces or through memory interference the other learned memories are disrupted as well. In this thesis, using a modified cheeseboard paradigm and a closed-loop recording setup with feedback optogenetic stimulation, I examined how the disruption of the reactivation of specific spiking patterns affects consolidation of the corresponding memory traces. To obtain multiple distinctive memories, animals had to perform a spatial task in two distinct cheeseboard environments and the reactivation of spiking patterns associated with one of the environments (target) was disrupted after learning during four hours rest period using a real-time decoding method. This real-time decoding method was capable of selectively affecting the firing rates and cofiring correlations of the target environment-encoding cells. The selective disruption led to behavioural impairment in the memory tests after the rest periods in the target environment but not in the other undisrupted control environment. In addition, the map of the target environment was less stable in the impaired memory tests compared to the learning session before than the map of the control environment. However, when the animal relearned the task, the same map recurred in the target environment that was present during learning before the disruption. Altogether my work demonstrated that the reactivation content is important: assembly-related disruption of reactivation can lead to a selective memory impairment and deficiency in map stability. These findings indeed suggest that reactivated assembly patterns reflect processes associated with the consolidation of memory traces

Low-rank graphical models and Bayesian inference in the statistical analysis of noisy neural data

We develop new methods of Bayesian inference, largely in the context of analysis of neuroscience data. The work is broken into several parts. In the first part, we introduce a novel class of joint probability distributions in which exact inference is tractable. Previously it has been difficult to find general constructions for models in which efficient exact inference is possible, outside of certain classical cases. We identify a class of such models that are tractable owing to a certain "low-rank" structure in the potentials that couple neighboring variables. In the second part we develop methods to quantify and measure information loss in analysis of neuronal spike train data due to two types of noise, making use of the ideas developed in the first part. Information about neuronal identity or temporal resolution may be lost during spike detection and sorting, or precision of spike times may be corrupted by various effects. We quantify the information lost due to these effects for the relatively simple but sufficiently broad class of Markovian model neurons. We find that decoders that model the probability distribution of spike-neuron assignments significantly outperform decoders that use only the most likely spike assignments. We also apply the ideas of the low-rank models from the first section to defining a class of prior distributions over the space of stimuli (or other covariate) which, by conjugacy, preserve the tractability of inference. In the third part, we treat Bayesian methods for the estimation of sparse signals, with application to the locating of synapses in a dendritic tree. We develop a compartmentalized model of the dendritic tree. Building on previous work that applied and generalized ideas of least angle regression to obtain a fast Bayesian solution to the resulting estimation problem, we describe two other approaches to the same problem, one employing a horseshoe prior and the other using various spike-and-slab priors. In the last part, we revisit the low-rank models of the first section and apply them to the problem of inferring orientation selectivity maps from noisy observations of orientation preference. The relevant low-rank model exploits the self-conjugacy of the von Mises distribution on the circle. Because the orientation map model is loopy, we cannot do exact inference on the low-rank model by the forward backward algorithm, but block-wise Gibbs sampling by the forward backward algorithm speeds mixing. We explore another von Mises coupling potential Gibbs sampler that proves to effectively smooth noisily observed orientation maps

Transductive neural decoding for unsorted neuronal spikes of rat hippocampus

Neural decoding is an important approach for extracting information from population codes. We previously proposed a novel transductive neural decoding paradigm and applied it to reconstruct the rat's position during navigation based on unsorted rat hippocampal ensemble spiking activity. Here, we investigate several important technical issues of this new paradigm using one data set of one animal. Several extensions of our decoding method are discussed.National Institutes of Health (U.S.) (NIH Grant DP1-OD003646)National Institutes of Health (U.S.) (NIH Grant MH061976