Fast online deconvolution of calcium imaging data

Abstract

Fluorescent calcium indicators are a popular means for observing the spiking activity of large neuronal populations, but extracting the activity of each neuron from raw fluorescence calcium imaging data is a nontrivial problem. We present a fast online active set method to solve this sparse non-negative deconvolution problem. Importantly, the algorithm 3progresses through each time series sequentially from beginning to end, thus enabling real-time online estimation of neural activity during the imaging session. Our algorithm is a generalization of the pool adjacent violators algorithm (PAVA) for isotonic regression and inherits its linear-time computational complexity. We gain remarkable increases in processing speed: more than one order of magnitude compared to currently employed state of the art convex solvers relying on interior point methods. Unlike these approaches, our method can exploit warm starts; therefore optimizing model hyperparameters only requires a handful of passes through the data. A minor modification can further improve the quality of activity inference by imposing a constraint on the minimum spike size. The algorithm enables real-time simultaneous deconvolution of O(105) traces of whole-brain larval zebrafish imaging data on a laptop.

Full text

Its licence does not allow us to show the full text. Read it at PubMed Central:

References 50

1. 1.

   Grienberger C, Konnerth C. Imaging calcium in neurons. Neuron. 2012;73(5):862–885.

   Not in our library10.1016/j.neuron.2012.02.011

2. 2.

   Grewe BF, Langer D, Kasper H, Kampa BM, Helmchen F. High-speed in vivo calcium imaging reveals neuronal network activity with near-millisecond precision. Nat Methods. 2010;7(5):399–405.

   Not in our library10.1038/nmeth.1453

3. 3.

   Yaksi E, Friedrich RW. Reconstruction of firing rate changes across neuronal populations by temporally deconvolved Ca2+ imaging. Nat Methods. 2006;3(5):377–383.

   Not in our library10.1038/nmeth874

4. 4.

   Holekamp TF, Turaga D, Holy TE. Fast three-dimensional fluorescence imaging of activity in neural populations by objective-coupled planar illumination microscopy. Neuron. 2008;57(5):661–672.

   Not in our library10.1016/j.neuron.2008.01.011

5. 5.

   Vogelstein JT, Packer AM, Machado TA, Sippy T, Babadi B, Yuste R, et al Fast nonnegative deconvolution for spike train inference from population calcium imaging. J Neurophysiol. 2010;104(6):3691–3704.

   Not in our library10.1152/jn.01073.2009

6. 6.

   Vogelstein JT, Watson BO, Packer AM, Yuste R, Jedynak B, Paninski L. Spike inference from calcium imaging using sequential Monte Carlo methods. Biophys J. 2009;97(2):636–655.

   Not in our library10.1016/j.bpj.2008.08.005

7. 7.

   Pnevmatikakis, EA, Merel, J, Pakman, A, Paninski, L. Bayesian spike inference from calcium imaging data. Asilomar Conference on Signals, Systems and Computers. 2013;.

   Not in our library

8. 8.
   Accurate spike estimation from noisy calcium signals for ultrafast three-dimensional imaging of large neuronal populations in vivo · 2016

   No notice
9. 9.

   Sasaki T, Takahashi N, Matsuki N, Ikegaya Y. Fast and accurate detection of action potentials from somatic calcium fluctuations. J Neurophysiol. 2008;100(3):1668–1676.

   Not in our library10.1152/jn.00084.2008

10. 10.

    Theis L, Berens P, Froudarakis E, Reimer J, Rosón MR, Baden T, et al Benchmarking spike rate inference in population calcium imaging. Neuron. 2016;90(3):471–482.

    Not in our library10.1016/j.neuron.2016.04.014

11. 11.

    Mishchencko Y, Vogelstein JT, Paninski L. A Bayesian approach for inferring neuronal connectivity from calcium fluorescent imaging data. Ann Appl Stat. 2011; p. 1229–1261.

    Not in our library10.1214/09-aoas303

12. 12.

    Picardo MA, Merel J, Katlowitz KA, Vallentin D, Okobi DE, Benezra SE, et al Population-Level Representation of a Temporal Sequence Underlying Song Production in the Zebra Finch. Neuron. 2016;90(4):866–876.

    Not in our library10.1016/j.neuron.2016.02.016

13. 13.

    Pnevmatikakis EA, Soudry D, Gao Y, Machado TA, Merel J, Pfau D, et al Simultaneous Denoising, Deconvolution, and Demixing of Calcium Imaging Data. Neuron. 2016;89(2):285–299.

    Not in our library10.1016/j.neuron.2015.11.037

14. 14.

    Grosenick L, Marshel JH, Deisseroth K. Closed-loop and activity-guided optogenetic control. Neuron. 2015;86(1):106–139.

    Not in our library10.1016/j.neuron.2015.03.034

15. 15.

    Rickgauer JP, Deisseroth K, Tank DW. Simultaneous cellular-resolution optical perturbation and imaging of place cell firing fields. Nat Neurosci. 2014;17(12):1816–1824.

    Not in our library10.1038/nn.3866

16. 16.

    Packer AM, Russell LE, Dalgleish HW, Häusser M. Simultaneous all-optical manipulation and recording of neural circuit activity with cellular resolution in vivo. Nat Methods. 2015;12(2):140–146.

    Not in our library10.1038/nmeth.3217

17. 17.

    Clancy KB, Koralek AC, Costa RM, Feldman DE, Carmena JM. Volitional modulation of optically recorded calcium signals during neuroprosthetic learning. Nat Neurosci. 2014;17(6):807–809.

    Not in our library10.1038/nn.3712

18. 18.

    Lewi J, Butera R, Paninski L. Sequential optimal design of neurophysiology experiments. Neural Comput. 2009;21(3):619–687.

    Not in our library10.1162/neco.2008.08-07-594

19. 19.

    Park M, Pillow JW. Bayesian active learning with localized priors for fast receptive field characterization. In: Adv Neural Inf Process Syst; 2012 p. 2348–2356.

    Not in our library

20. 20.

    Shababo B, Paige B, Pakman A, Paninski L. Bayesian inference and online experimental design for mapping neural microcircuits. In: Adv Neural Inf Process Syst; 2013 p. 1304–1312.

    Not in our library

21. 21.

    Ahrens MB, Orger MB, Robson DN, Li JM, Keller PJ. Whole-brain functional imaging at cellular resolution using light-sheet microscopy. Nat Methods. 2013;10(5):413–420.

    Not in our library10.1038/nmeth.2434

22. 22.

    Vladimirov N, Mu Y, Kawashima T, Bennett DV, Yang CT, Looger LL, et al Light-sheet functional imaging in fictively behaving zebrafish. Nat Methods. 2014;.

    Not in our library10.1038/nmeth.3040

23. 23.

    Potra FA, Wright SJ. Interior-point methods. J Comput Appl Math. 2000;124(1):281–302.

    Not in our library10.1016/s0377-0427(00)00433-7

24. 24.

    Ayer M, Brunk HD, Ewing GM, Reid WT, Silverman E. An empirical distribution function for sampling with incomplete information. Ann Math Stat. 1955;26(4):641–647.

    Not in our library10.1214/aoms/1177728423

25. 25.

    Barlow RE, Bartholomew DJ, Bremner JM, Brunk HD. Statistical inference under order restrictions: The theory and application of isotonic regression. Wiley New York; 1972.

    Not in our library

26. 26.

    van Eeden C. Testing and Estimating Ordered Parameters of Probability Distributions. PhD thesis, University of Amsterdam; 1958.

    Not in our library

27. 27.

    Miles RE. The complete amalgamation into blocks, by weighted means, of a finite set of real numbers. Biometrika. 1959;46(3/4):317–327.

    Not in our library10.1093/biomet/46.3-4.317

28. 28.

    Mair P, Hornik K, de Leeuw J. Isotone optimization in R: pool-adjacent-violators algorithm (PAVA) and active set methods. J Stat Softw. 2009;32(5):1–24.

    Not in our library

29. 29.

    Best MJ, Chakravarti N. Active set algorithms for isotonic regression; a unifying framework. Math Prog. 1990;47(1–3):425–439.

    Not in our library10.1007/bf01580873

30. 30.

    Grotzinger SJ, Witzgall C. Projections onto order simplexes. Appl Math Optim. 1984;12(1):247–270.

    Not in our library10.1007/bf01449044

31. 31.

    Podgorski K, Haas K. Fast non-negative temporal deconvolution for laser scanning microscopy. J Biophotonics. 2013;6(2):153–162.

    Not in our library10.1002/jbio.201100133

32. 32.

    Brent RP. Algorithms for Minimization Without Derivatives. Courier Corporation; 1973.

    Not in our library

33. 33.

    Nocedal J, Wright S. Numerical optimization. Springer Science & Business Media; 2006.

    Not in our library

34. 34.

    Donoho DL. De-noising by soft-thresholding. IEEE Trans Inf Theory. 1995;41(3):613–627.

    Not in our library10.1109/18.382009

35. 35.

    Mallat SG, Zhang Z. Matching pursuits with time-frequency dictionaries. IEEE Trans Signal Processing. 1993;41(12):3397–3415.

    Not in our library10.1109/78.258082

36. 36.

    Chen TW, Wardill T, Sun Y, Pulver SR, Renninger SL, Baohan A, et al Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature. 2013;499(7458):295–300.

    Not in our library10.1038/nature12354

37. 37.

    Lawson CL, Hanson RJ. Solving least squares problems. vol. 15 SIAM; 1995.

    Not in our library

38. 38.

    Bro R, De Jong S. A fast non-negativity-constrained least squares algorithm. J Chemometrics. 1997;11(5):393–401.

    Not in our library10.1002/(sici)1099-128x(199709/10)11:5<393∷aid-cem483>3.3.co;2-c

39. 39.

    Diamond S, Boyd S. CVXPY: A Python-Embedded Modeling Language for Convex Optimization. J Mach Learn Res. 2016;17(83):1–5.

    Not in our library

40. 40.

    Domahidi A, Chu E, Boyd S. ECOS: An SOCP solver for embedded systems. In: European Control Conference (ECC); 2013. p. 3071–3076.

    Not in our library

Cited by papers we hold 13

• Glutamate indicators with improved activation kinetics and localization for imaging synaptic transmission · 2023
  No notice
• Facemap: a framework for modeling neural activity based on orofacial tracking · 2023
  No notice
• Evaluation of Recombinant Botulinum Neurotoxin Type A1 Efficacy in Peripheral Inflammatory Pain in Mice · 2022
  Not checked
• Minian, an open-source miniscope analysis pipeline · 2022
  Not checked
• Gut–brain circuits for fat preference · 2022
  No notice
• Large-scale two-photon calcium imaging in freely moving mice · 2022
  No notice
• Relationship between simultaneously recorded spiking activity and fluorescence signal in GCaMP6 transgenic mice · 2021
  No notice
• Calcium Imaging and the Curse of Negativity · 2021
  Not checked
• The importance of non-pharmaceutical interventions during the COVID-19 vaccine rollout. · 2021
  Not checked
• CaImAn an open source tool for scalable calcium imaging data analysis · 2019
  No notice
• Efficient and accurate extraction of in vivo calcium signals from microendoscopic video data · 2018
  No notice
• Robustness of Spike Deconvolution for Neuronal Calcium Imaging · 2018
  No notice
• Community-based benchmarking improves spike rate inference from two-photon calcium imaging data · 2018
  No notice

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.