Faculty Bibliography

Pyramidal cells in the rodent hippocampus often exhibit clear spatial tuning in navigation. Although it has been long suggested that pyramidal cell activity may underlie a topological code rather than a topographic code, it remains unclear whether an abstract spatial topology can be encoded in the ensemble spiking activity of hippocampal place cells. Using a statistical approach developed previously, we investigate this question and related issues in greater detail. We recorded ensembles of hippocampal neurons as rodents freely foraged in one- and two-dimensional spatial environments and used a "decode-to-uncover" strategy to examine the temporally structured patterns embedded in the ensemble spiking activity in the absence of observed spatial correlates during periods of rodent navigation or awake immobility. Specifically, the spatial environment was represented by a finite discrete state space. Trajectories across spatial locations ("states") were associated with consistent hippocampal ensemble spiking patterns, which were characterized by a state transition matrix. From this state transition matrix, we inferred a topology graph that defined the connectivity in the state space. In both one- and two-dimensional environments, the extracted behavior patterns from the rodent hippocampal population codes were compared against randomly shuffled spike data. In contrast to a topographic code, our results support the efficiency of topological coding in the presence of sparse sample size and fuzzy space mapping. This computational approach allows us to quantify the variability of ensemble spiking activity, examine hippocampal population codes during off-line states, and quantify the topological complexity of the environment.

A fundamental task in neuroscience is to understand how neural ensembles represent information. Population decoding is a useful tool to extract information from neuronal populations based on the ensemble spiking activity. We propose a novel Bayesian decoding paradigm to decode unsorted spikes in the rat hippocampus. Our approach uses a direct mapping between spike waveform features and covariates of interest and avoids accumulation of spike sorting errors. Our decoding paradigm is nonparametric, encoding model-free for representing stimuli, and extracts information from all available spikes and their waveform features. We apply the proposed Bayesian decoding algorithm to a position reconstruction task for freely behaving rats based on tetrode recordings of rat hippocampal neuronal activity. Our detailed decoding analyses demonstrate that our approach is efficient and better utilizes the available information in the nonsortable hash than the standard sorting-based decoding algorithm. Our approach can be adapted to an online encoding/decoding framework for applications that require real-time decoding, such as brain-machine interfaces.

Behavioral state is known to influence interactions between thalamus and cortex, which are important for sensation, action, and cognition. The thalamic reticular nucleus (TRN) is hypothesized to regulate thalamo-cortical interactions, but the underlying functional architecture of this process and its state dependence are unknown. By combining the first TRN ensemble recording with psychophysics and connectivity-based optogenetic tagging, we found reticular circuits to be composed of distinct subnetworks. While activity of limbic-projecting TRN neurons positively correlates with arousal, sensory-projecting neurons participate in spindles and show elevated synchrony by slow waves during sleep. Sensory-projecting neurons are suppressed by attentional states, demonstrating that their gating of thalamo-cortical interactions is matched to behavioral state. Bidirectional manipulation of attentional performance was achieved through subnetwork-specific optogenetic stimulation. Together, our findings provide evidence for differential inhibition of thalamic nuclei across brain states, where the TRN separately controls external sensory and internal limbic processing facilitating normal cognitive function. PAPERFLICK:

Sparse Bayesian inference methods are applied to decode three-dimensional (3D) reach to grasp movement based on recordings of primary motor cortical (M1) ensembles from rhesus macaque. For three linear or nonlinear models tested, variational Bayes (VB) inference in combination with automatic relevance determination (ARD) is used for variable selection to avoid overfitting. The sparse Bayesian linear regression model achieved the overall best performance across objects and target locations. We assessed the sensitivity of M1 units in decoding and evaluated the proximal and distal representations of joint angles in population decoding. Our results suggest that the M1 ensembles recorded from the precentral gyrus area carry more proximal than distal information.

Neural spike train analysis is an important task in computational neuroscience which aims to understand neural mechanisms and gain insights into neural circuits. With the advancement of multielectrode recording and imaging technologies, it has become increasingly demanding to develop statistical tools for analyzing large neuronal ensemble spike activity. Here we present a tutorial overview of Bayesian methods and their representative applications in neural spike train analysis, at both single neuron and population levels. On the theoretical side, we focus on various approximate Bayesian inference techniques as applied to latent state and parameter estimation. On the application side, the topics include spike sorting, tuning curve estimation, neural encoding and decoding, deconvolution of spike trains from calcium imaging signals, and inference of neuronal functional connectivity and synchrony. Some research challenges and opportunities for neural spike train analysis are discussed.

Rodent hippocampal population codes represent important spatial information of 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. Specifically, we develop an infinite hidden Markov model (iHMM) and variational Bayes (VB) inference method to analyze rat hippocampal ensemble spike activity. We demonstrate the effectiveness of our approach using an open field navigation example and discuss the significance/implications of our results.

Respiratory sinus arrhythmia (RSA) is largely mediated by the autonomic nervous system through its modulating influence on the heart beats. We propose a robust algorithm for quantifying instantaneous RSA as applied to heart beat intervals and respiratory recordings under dynamic breathing patterns. The blood volume pressure-derived heart beat series (pulse intervals, PIs) are modeled as an inverse Gaussian point process, with the instantaneous mean PI modeled as a bivariate regression incorporating both past PIs and respiration values observed at the beats. A point process maximum likelihood algorithm is used to estimate the model parameters, and instantaneous RSA is estimated via a frequency domain transfer function evaluated at instantaneous respiratory frequency where high coherence between respiration and PIs is observed. The model is statistically validated using Kolmogorov-Smirnov goodness-of-fit analysis, as well as independence tests. The algorithm is applied to subjects engaged in meditative practice, with distinctive dynamics in the respiration patterns elicited as a result. The presented analysis confirms the ability of the algorithm to track important changes in cardiorespiratory interactions elicited during meditation, otherwise not evidenced in control resting states, reporting statistically significant increase in RSA gain as measured by our paradigm.