Precise hyperacuity estimation of spike timing from calcium imaging

Two-photon imaging is a major recording technique in neuroscience, but it suffers from several limitations, including a low sampling rate, the nonlinearity of calcium responses, the slow dynamics of calcium dyes and a low signal-to-noise ratio, all of which impose a severe limitation on the application of two-photon imaging in elucidating neuronal dynamics with high temporal resolution. Here, we developed a hyperacuity algorithm (HA_time) based on an approach combining a generative model and machine learning to improve spike detection and the precision of spike time inference. First, Bayesian inference estimates the calcium spike model by assuming the constancy of the spike shape and size. A support vector machine employs this information and detects spikes with higher temporal precision than the sampling rate. Compared with conventional thresholding, HA_time improved the precision of spike time estimation up to 20-fold for simulated calcium data. Furthermore, the benchmark analysis of experimental data from different brain regions and simulation of a broader range of experimental conditions showed that our algorithm was among the best in a class of hyperacuity algorithms. We encourage experimenters to use the proposed algorithm to precisely estimate hyperacuity spike times from two-photon imaging.