Inferring Neuronal Dynamics from Calcium Imaging Data Using Biophysical Models and Bayesian Inference

33Citations
Citations of this article
144Readers
Mendeley users who have this article in their library.

Abstract

Calcium imaging has been used as a promising technique to monitor the dynamic activity of neuronal populations. However, the calcium trace is temporally smeared which restricts the extraction of quantities of interest such as spike trains of individual neurons. To address this issue, spike reconstruction algorithms have been introduced. One limitation of such reconstructions is that the underlying models are not informed about the biophysics of spike and burst generations. Such existing prior knowledge might be useful for constraining the possible solutions of spikes. Here we describe, in a novel Bayesian approach, how principled knowledge about neuronal dynamics can be employed to infer biophysical variables and parameters from fluorescence traces. By using both synthetic and in vitro recorded fluorescence traces, we demonstrate that the new approach is able to reconstruct different repetitive spiking and/or bursting patterns with accurate single spike resolution. Furthermore, we show that the high inference precision of the new approach is preserved even if the fluorescence trace is rather noisy or if the fluorescence transients show slow rise kinetics lasting several hundred milliseconds, and inhomogeneous rise and decay times. In addition, we discuss the use of the new approach for inferring parameter changes, e.g. due to a pharmacological intervention, as well as for inferring complex characteristics of immature neuronal circuits.

Figures

  • Fig 1. The generativemodel and non-hybrid quadratic integrate-and-fire (family of QGIF) models. (A) Graph representing the generative model and its inversion, which are comprised of evolution (i.e., a neuron model) and observation equations. The illustrated hierarchy in the graph displays how neuronal dynamics relate to fluorescence traces. (B andC) Sample voltage traces of (B) QGIF and (C) bursting-QGIF models. The traces show the rhythmic activity of these models in response to sustained depolarizations. Note that the family of QGIF models does not require any reset condition for spike/burst generation. The persistent Na+ and M-type K+ currents (red and green lines, respectively) of bursting-QGIF model are in units of [μA / cm2]. Parameters for the simulations: (B) Iapp = 0.2 μA / cm 2, (C) Irep ¼ 80 mA mS=mF, Iapp = 1 μA / cm2. See Table 1 for the rest of parameter values.
  • Fig 2. Assessment of biophysical aspects of the bursting-QGIF model. (A-C) Simulated voltage traces of the bursting-QGIF model in response to (A) sustained (not shown) and (B) brief square (magenta line) positive current pulses, and (C) a sinusoidal input current (magenta line). (A) The increment in persistent Na+ current enhances the burstiness, similarly to [57]: As gNaP increases, the number of spikes within each burst is increased, and the interspike intervals become shorter. Note that as expected biophysically the model exhibits a tonic repetitive spiking pattern for weak gNaP densities. (B) M-type K + current governs the recovery mechanism for membrane potential, similarly to [57]: Blocking the M-type K+ channel by setting gM = 0 leads to a prolonged burst in response to a short-duration depolarized current pulse (magenta line). This indicates that the activation of this channel is important for terminating the bursts. (C) The bursting-QGIF model is selective in the input slope, similarly to [53]: the periodic burst response of the model to a sinusoidal input current shows that the burst are mainly initiated on the positive slope of the input (magenta line) thus signalling the input slope. The conductances and currents acrossA-C are in units of [mS / cm2] and [μA / cm2], respectively; input currents are in arbitrary units. Parameters for the simulation: (A-C) Irep ¼ 80 mA mS=mF, (A) Iapp = 1, (B) Iapp = 0.6, (C) Iapp = 2.5sin(0.03t) μA / cm2. See Table 1 for the rest of parameter values.
  • Table 1. Parameter values.
  • Table 1. (Continued)
  • Fig 3. Field of view and low frequency temporal drifts. (A) Sample average (over frames) Oregon Green BAPTA 1 (OGB-1) fluorescence image of a neuronal population from the CA3 area of a hippocampal slice. (B) An in vitroOGB-1 fluorescence trace before (blue line) and after (red line) removal of its slowly varying components, by using the fourth degree polynomial detrending method.
  • Table 2. Prior densities.
  • Table 3. Generative models.
  • Fig 4. Results of the proposed approach for synthetic single-spike-evoked transients with fast rise times. (A): QGIF model, (B): FHNmodel with high SNR, (C): FHNmodel with low SNR. (A-C) Inferring neuronal dynamics from transients with fast rise times (e.g. 3ms). The fluorescence traces (first row) and [Ca2+] kinetics (second row) were simulated by using (A) QGIF and (B andC) FHN generative models. Inverting each model for its corresponding trace by

References Powered by Scopus

A quantitative description of membrane current and its application to conduction and excitation in nerve

16308Citations
N/AReaders
Get full text

Impulses and Physiological States in Theoretical Models of Nerve Membrane

4907Citations
N/AReaders
Get full text

Ultrasensitive fluorescent proteins for imaging neuronal activity

4567Citations
N/AReaders
Get full text

Cited by Powered by Scopus

Calcium imaging and dynamic causal modelling reveal brain-wide changes in effective connectivity and synaptic dynamics during epileptic seizures

46Citations
N/AReaders
Get full text

Electrophysiology read-out tools for brain-on-chip biotechnology

38Citations
N/AReaders
Get full text

Insights Into Spinal Dorsal Horn Circuit Function and Dysfunction Using Optical Approaches

30Citations
N/AReaders
Get full text

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Cite

CITATION STYLE

APA

Rahmati, V., Kirmse, K., Marković, D., Holthoff, K., & Kiebel, S. J. (2016). Inferring Neuronal Dynamics from Calcium Imaging Data Using Biophysical Models and Bayesian Inference. PLoS Computational Biology, 12(2). https://doi.org/10.1371/journal.pcbi.1004736

Readers over time

‘16‘17‘18‘19‘20‘21‘22‘23‘24010203040

Readers' Seniority

Tooltip

PhD / Post grad / Masters / Doc 71

64%

Researcher 32

29%

Professor / Associate Prof. 7

6%

Lecturer / Post doc 1

1%

Readers' Discipline

Tooltip

Neuroscience 45

48%

Agricultural and Biological Sciences 19

20%

Computer Science 15

16%

Engineering 14

15%

Article Metrics

Tooltip
Mentions
News Mentions: 1
Social Media
Shares, Likes & Comments: 127

Save time finding and organizing research with Mendeley

Sign up for free
0