Spike inference from calcium imaging using sequential Monte Carlo methods. Vogelstein, J., T., Watson, B., O., Packer, A., M., Yuste, R., Jedynak, B., M., & Paninski, L. Biophysical Journal, 97(2):636-655, 7, 2009.
abstract   bibtex   
As recent advances in calcium sensing technologies facilitate simultaneously imaging action potentials in neuronal populations, complementary analytical tools must also be developed to maximize the utility of this experimental paradigm. Although the observations here are fluorescence movies, the signals of interest--spike trains and/or time varying intracellular calcium concentrations--are hidden. Inferring these hidden signals is often problematic due to noise, nonlinearities, slow imaging rate, and unknown biophysical parameters. We overcome these difficulties by developing sequential Monte Carlo methods (particle filters) based on biophysical models of spiking, calcium dynamics, and fluorescence. We show that even in simple cases, the particle filters outperform the optimal linear (i.e., Wiener) filter, both by obtaining better estimates and by providing error bars. We then relax a number of our model assumptions to incorporate nonlinear saturation of the fluorescence signal, as well external stimulus and spike history dependence (e.g., refractoriness) of the spike trains. Using both simulations and in vitro fluorescence observations, we demonstrate temporal superresolution by inferring when within a frame each spike occurs. Furthermore, the model parameters may be estimated using expectation maximization with only a very limited amount of data (e.g., approximately 5-10 s or 5-40 spikes), without the requirement of any simultaneous electrophysiology or imaging experiments.
@article{
 title = {Spike inference from calcium imaging using sequential Monte Carlo methods.},
 type = {article},
 year = {2009},
 identifiers = {[object Object]},
 keywords = {animals,biological,calcium,calcium metabolism,cytology/metabolism,fluorescence,inbred c57bl,intracellular space,intracellular space metabolism,metabolism,mice,models,monte carlo method,neurons,neurons cytology,neurons metabolism,probability,time factors},
 pages = {636-655},
 volume = {97},
 websites = {http://dx.doi.org/10.1016/j.bpj.2008.08.005},
 month = {7},
 institution = {Department of Neuroscience, The Johns Hopkins School of Medicine, Baltimore, Maryland, USA. joshuav@jhu.edu},
 id = {ea71e181-8857-33ea-b063-d3d46eb47138},
 created = {2015-06-19T07:45:09.000Z},
 file_attached = {false},
 profile_id = {182bbbf9-24a3-3af3-9ed6-563e8f89259b},
 group_id = {8d229673-0aec-3014-b0f6-eda47f83e147},
 last_modified = {2015-06-19T07:45:23.000Z},
 read = {false},
 starred = {false},
 authored = {true},
 confirmed = {true},
 hidden = {false},
 citation_key = {smc-oopsi},
 source_type = {article},
 abstract = {As recent advances in calcium sensing technologies facilitate simultaneously imaging action potentials in neuronal populations, complementary analytical tools must also be developed to maximize the utility of this experimental paradigm. Although the observations here are fluorescence movies, the signals of interest--spike trains and/or time varying intracellular calcium concentrations--are hidden. Inferring these hidden signals is often problematic due to noise, nonlinearities, slow imaging rate, and unknown biophysical parameters. We overcome these difficulties by developing sequential Monte Carlo methods (particle filters) based on biophysical models of spiking, calcium dynamics, and fluorescence. We show that even in simple cases, the particle filters outperform the optimal linear (i.e., Wiener) filter, both by obtaining better estimates and by providing error bars. We then relax a number of our model assumptions to incorporate nonlinear saturation of the fluorescence signal, as well external stimulus and spike history dependence (e.g., refractoriness) of the spike trains. Using both simulations and in vitro fluorescence observations, we demonstrate temporal superresolution by inferring when within a frame each spike occurs. Furthermore, the model parameters may be estimated using expectation maximization with only a very limited amount of data (e.g., approximately 5-10 s or 5-40 spikes), without the requirement of any simultaneous electrophysiology or imaging experiments.},
 bibtype = {article},
 author = {Vogelstein, Joshua T. and Watson, Brendon O and Packer, Adam M and Yuste, Rafael and Jedynak, Bruno M and Paninski, Liam},
 journal = {Biophysical Journal},
 number = {2}
}

Downloads: 0