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

41. 41.

    Andersen ED, Andersen KD. The MOSEK interior point optimizer for linear programming: an implementation of the homogeneous algorithm In: High performance optimization. Springer; 2000 p. 197–232.

    Not in our library

42. 42.

    O’Donoghue B, Chu E, Parikh N, Boyd S. Conic Optimization via Operator Splitting and Homogeneous Self-Dual Embedding. J Optim Theory Appl. 2016; p. 1–27.

    Not in our library

43. 43.

    Gurobi Optimization Inc. Gurobi Optimizer Reference Manual; 2015. Available from: http://www.gurobi.com.

    Not in our library

44. 44.

    GENIE project, Janelia Research Campus, HHMI; Karel Svoboda (contact). Simultaneous imaging and loose-seal cell-attached electrical recordings from neurons expressing a variety of genetically encoded calcium indicators. CRCNS.org; 2015. Available from:.

    Not in our library10.6080/k02r3pmn

45. 45.

    Pologruto TA, Yasuda R, Svoboda K. Monitoring neural activity and [Ca2+] with genetically encoded Ca2+ indicators. J Neurosci. 2004;24(43):9572–9579.

    Not in our library10.1523/jneurosci.2854-04.2004

46. 46.

    Quan T, Lv X, Liu X, Zeng S. Reconstruction of burst activity from calcium imaging of neuronal population via Lq minimization and interval screening. Biomed Opt Express. 2016;7(6):2103–2117.

    Not in our library10.1364/boe.7.002103

47. 47.

    Dyer EL, Studer C, Robinson JT, Baraniuk RG. A robust and efficient method to recover neural events from noisy and corrupted data. In: Int IEEE/EMBS Conf Neural Eng (NER). 2013; p. 593–596.

    Not in our library

48. 48.

    Pachitariu M, Stringer C, Schröder S, Dipoppa M, Rossi LF, Carandini M, et al. Suite2p: beyond 10,000 neurons with standard two-photon microscopy. bioRxiv. 2016.

    Not in our library

49. 49.

    Giovannucci A, Friedrich J, Deverett B, Staneva V, Chklovskii D, Pnevmatikakis E. CaImAn: An open source toolbox for large scale calcium imaging data analysis on standalone machines Cosyne Abstracts 2017, Salt Lake City USA.

    Not in our library

50. 50.

    Friedrich J, Paninski L. Fast Active Set Methods for Online Spike Inference from Calcium Imaging. In: Adv Neural Inf Process Syst; 2016 p. 1984–1992.

    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.