Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, August 13 until 12:00 AM ET on Friday, August 14 due to maintenance. We apologize for the inconvenience.


Title: Robust point-process Granger causality analysis in presence of exogenous temporal modulations and trial-by-trial variability in spike trains
Assessing directional influences between neurons is instrumental to understand how brain circuits process information. To this end, Granger causality, a technique originally developed for time-continuous signals, has been extended to discrete spike trains. A fundamental assumption of this technique is that the temporal evolution of neuronal responses must be due only to endogenous interactions between recorded units, including self-interactions. This assumption is however rarely met in neurophysiological studies, where the response of each neuron is modulated by other exogenous causes such as, for example, other unobserved units or slow adaptation processes. Here, we propose a novel point-process Granger causality technique that is robust with respect to the two most common exogenous modulations observed in real neuronal responses: within-trial temporal variations in spiking rate and between-trial variability in their magnitudes. This novel method works by explicitly including both types of modulations into the generalized linear model of the neuronal conditional intensity function (CIF). We then assess the causal influence of neuron i onto neuron j by measuring the relative reduction of neuron j ’s point process likelihood obtained considering or removing neuron i . CIF’s hyper-parameters are set on a per-neuron basis by minimizing Akaike’s information criterion. In synthetic data sets, generated by means of random processes or networks of integrate-and-fire units, the proposed method recovered with high accuracy, sensitivity and robustness the underlying ground-truth connectivity pattern. Application of presently available point-process Granger causality techniques produced instead a significant number of false positive connections. In real spiking responses recorded from neurons in the monkey pre-motor cortex (area F5), our method revealed many causal relationships between neurons as well as the temporal structure of their interactions. Given its robustness our method can be effectively applied to real neuronal data. Furthermore, its explicit estimate of the effects of unobserved causes on the recorded neuronal firing patterns can help decomposing their temporal variations into endogenous and exogenous components.  more » « less
Award ID(s):
1942585 1755780 2226123
PAR ID:
10220612
Author(s) / Creator(s):
; ;
Editor(s):
Morrison, Abigail
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
17
Issue:
1
ISSN:
1553-7358
Page Range / eLocation ID:
e1007675
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Extracting directional connectivity in a neuronal ensemble from spiking observations is a key challenge in understanding the circuit mechanisms of brain function. Existing methods proceed in two stages, by first estimating the latent processes that govern spiking, followed by characterizing connectivity using said estimates. As such, the extracted networks in the second stage are highly sensitive to the accuracy of the estimates in the first stage. In this work, we introduce a framework to directly extract Granger causal links from spiking observations, without requiring intermediate time-domain estimation, by explicitly modeling the endogenous and exogenous latent processes that underlie spiking activity. Our proposed method integrates several techniques such as point processes, state-space modeling and Pólya-Gamma augmentation. We demonstrate the utility of our proposed approach using simulated data and application to real data from the rat brain during sleep. 
    more » « less
  2. null (Ed.)
    Neuronal activity correlations are key to understanding how populations of neurons collectively encode information. While two-photon calcium imaging has created a unique opportunity to record the activity of large populations of neurons, existing methods for inferring correlations from these data face several challenges. First, the observations of spiking activity produced by two-photon imaging are temporally blurred and noisy. Secondly, even if the spiking data were perfectly recovered via deconvolution, inferring network-level features from binary spiking data is a challenging task due to the non-linear relation of neuronal spiking to endogenous and exogenous inputs. In this work, we propose a methodology to explicitly model and directly estimate signal and noise correlations from two-photon fluorescence observations, without requiring intermediate spike deconvolution. We provide theoretical guarantees on the performance of the proposed estimator and demonstrate its utility through applications to simulated and experimentally recorded data from the mouse auditory cortex. 
    more » « less
  3. The organization of neurons into functionally related assemblies is a fundamental feature of cortical networks, yet our understanding of how these assemblies maintain distinct identities while sharing members remains limited. Here we analyze how spike timing-dependent plasticity (STDP) shapes the formation and stability of overlapping neuronal assemblies in recurrently coupled networks of spiking neuron models. Using numerical simulations and an associated mean-field theory, we demonstrate that the temporal structure of the STDP rule, specifically its degree of causality, critically determines whether assemblies that share neurons maintain segregation or merge together after training is completed. We find that causal STDP rules, where potentiation or depression occurs strictly when presynaptic spikes precede or follow postsynaptic spikes, allow assemblies to remain distinct even with substantial overlap in membership. This stability arises because causal STDP effectively cancels the symmetric correlations introduced by common inputs from shared neurons. In contrast, acausal STDP rules lead to assembly fusion when overlap exceeds a critical threshold, due to unchecked growth of common input correlations. Our results provide theoretical insight into how spike timing-dependent learning rules can support distributed representations where individual neurons participate in multiple assemblies while maintaining functional specificity. 
    more » « less
  4. Abstract: Neurons in higher-order visual areas integrate information through a canonical computation called normalization. The strength of normalization is highly heterogeneous across neurons, and this heterogeneity correlates with attention-mediated modulations in neural responses. However, the circuit mechanism underlying the heterogeneous normalization strength is unclear. In this work, we study normalization in a spiking neuron network model of visual cortex. Our model reveals that the heterogeneity of normalization strength is highly correlated with the inhibitory current each neuron receives. The correlation between inhibition and other synaptic inputs explains the experimentally observed dependence of spike count correlations on normalization strength. Further, we find that neurons with stronger normalization encode information more efficiently, and that networks with more heterogeneity in normalization have more information and higher capacity in encoding visual stimuli. Together, our model provides a mechanistic explanation of heterogeneous normalization strengths in the visual cortex, and sheds light on the computational benefits of neuronal heterogeneity. 
    more » « less
  5. Neurons in higher-order visual areas integrate information through a canonical computation called normalization. The strength of normalization is highly heterogeneous across neurons, and this heterogeneity correlates with attention-mediated modulations in neural responses. However, the circuit mechanism underlying the heterogeneous normalization strength is unclear. In this work, we study normalization in a spiking neuron network model of visual cortex. Our model reveals that the heterogeneity of normalization strength is highly correlated with the inhibitory current each neuron receives. The correlation between inhibition and other synaptic inputs explains the experimentally observed dependence of spike count correlations on normalization strength. Further, we find that neurons with stronger normalization encode information more efficiently, and that networks with more heterogeneity in normalization encode visual stimuli with higher information and capacity. Together, our model provides a mechanistic explanation of heterogeneous normalization strengths in the visual cortex and sheds light on the computational benefits of neuronal heterogeneity. 
    more » « less