Compression Detects Changes in Spiking Neural Data from Cortical Lesions

The complexity of neural data changes as the brain processes information during events. Universal lossless compression algorithms, which are broadly applicable and grounded in information theory, identify and exploit redundancies in data in order to compress it to essentially-optimal sizes regardless of underlying statistics. These algorithms may be used to conveniently and efficiently estimate a given signal's Shannon entropy rate, a biologically relevant measure of the complexity of a signal. It is therefore natural to explore their effectiveness in the analysis of spiking neural data.

This work focuses on using compression to analyze recordings (96-channel Utah arrays) taken from motor cortex of animals performing reaching tasks for three days before and three days after administering electrolytic lesions (Subject U: 4 lesions, H: 3). In particular, we use the inverse compression ratio (ICR), which compares the sizes of compressed and uncompressed data to estimate the amount of statistically unique information. We calculate ICR with temporally-independent lossless compression (gzip) and temporally-dependent lossy compression (H.264, MPEG-2). Compression-based ICR was compared to single-neuron measures used to understand spiking data, such as average firing rates and Fano factor. Compression is also compared to common dimensionality reduction techniques, principal component analysis (PCA) and factor analysis (FA).

Statistical tests on aggregate data comparing each metric before and after lesioning reveal that ICR is able to significantly (Mann-Whitney U test, p < 0.01) detect lesions with higher accuracy than single-neuron metrics, but not dimensionality reduction (ICR methods: 85.7%, single-neuron methods: 78.6%, dimensionality reduction: 100%). Additionally, statistical results on the same data show that ICR metrics remain more stable than single-neuron methods after lesion. The bitrate parameter of lossy compression algorithms is swept to better understand the effect of information rates and "optimal" compression on lesion detection performance. Our conclusions are confirmed by the same analyses performed on several different simulated neural datasets.

These results suggest that compression algorithms may be a useful tool to detect and better understand perturbations to the underlying structure of neural data. Information-theoretic analyses may complement techniques like dimensionality reduction and firing rate tuning as a convenient and useful tool to characterize neural data.