An Information-Theoretic Framework to Measure the Dynamic Interaction Between Neural Spike Trains

Objective: While understanding the interaction patterns among simultaneous recordings of spike trains from multiple neuronal units is a key topic in neuroscience, existing methods either do not consider the inherent point-process nature of spike trains or are based on parametric assumptions. This work presents an information-theoretic framework for the model-free, continuous-time estimation of both undirected (symmetric) and directed (Granger-causal) interactions between spike trains. Methods: The framework computes the mutual information rate (MIR) and the transfer entropy rate (TER) for two point processes X and Y, showing that the MIR between X and Y can be decomposed as the sum of the TER along the directions X. Y and Y. X. We present theoretical expressions and introduce strategies to estimate efficiently the two measures through nearest neighbor statistics. Results: Using simulations of independent and coupled point processes, we show the accuracy of MIR and TER to assess interactions even for weakly coupled and short realizations, and demonstrate the superiority of continuous-time estimation over the standard discrete-time approach. We also apply the MIR and TER to real-world data, specifically, recordings from in-vitro preparations of spontaneously-growing cultures of cortical neurons. Using this dataset, we demonstrate the ability of MIR and TER to describe how the functional networks between recording units emerge over the course of the maturation of the neuronal cultures. Conclusion and Significance: the proposed framework provides principled measures to assess undirected and directed spike train interactions with more efficiency and flexibility than previous discrete-time or parametric approaches, opening new perspectives for the analysis of point-process data in neuroscience and many other fields.

Automated summary

This study introduces an information-theoretic framework for model-free, continuous-time estimation of interactions between spike trains, demonstrating its accuracy for weakly coupled and short realizations compared to the standard discrete-time approach. By applying this framework to real-world data, it shows the ability to describe how functional networks between neuronal units emerge over time, providing more efficient and flexible measures for assessing spike train interactions in neuroscience.