Quantitative estimation of calcium dynamics from ratiometric measurements: a direct, nonratioing method

A simple method for getting standard error on the ratiometric calcium estimator

The ratiometric fluorescent calcium indicator Fura-2 plays a fundamental role in the investigation of cellular calcium dynamics. Despite of its widespread use in the last 30 years, only one publication (Joucla et al., 2010)) proposed a way of obtaining confidence intervals on fitted calcium dynamic model parameters from single ‘calcium transients’. Shortcomings of this approach are its requirement for a ‘3 wavelengths’ protocol (excitation at 340 and 380 nm as usual plus at 360 nm, the isosbestic point) as well as the need for an autofluorence / background fluorescence model at each wavelength. Here, we propose a simpler method that eliminates both shortcommings:1.

a precise estimation of the standard errors of the raw data is obtained first,

2.

the standard error of the ratiometric calcium estimator (a function of the raw data values) is derived using both the propagation of uncertainty and a Monte-Carlo method.

Once meaningful standard errors for calcium estimates are available, standard errors on fitted model parameters follow directly from the use of nonlinear least-squares optimization algorithms.

Estimating background-subtracted fluorescence transients in calcium imaging experiments: a quantitative approach

Calcium imaging has become a routine technique in neuroscience for subcellular to network level investigations. The fast progresses in the development of new indicators and imaging techniques call for dedicated reliable analysis methods. In particular, efficient and quantitative background fluorescence subtraction routines would be beneficial to most of the calcium imaging research field. A background-subtracted fluorescence transients estimation method that does not require any independent background measurement is therefore developed. This method is based on a fluorescence model fitted to single-trial data using a classical nonlinear regression approach. The model includes an appropriate probabilistic description of the acquisition system's noise leading to accurate confidence intervals on all quantities of interest (background fluorescence, normalized background-subtracted fluorescence time course) when background fluorescence is homogeneous. An automatic procedure detecting background inhomogeneities inside the region of interest is also developed and is shown to be efficient on simulated data. The implementation and performances of the proposed method on experimental recordings from the mouse hypothalamus are presented in details. This method, which applies to both single-cell and bulk-stained tissues recordings, should help improving the statistical comparison of fluorescence calcium signals between experiments and studies.

Photon shot noise limits on optical detection of neuronal spikes and estimation of spike timing

Optical approaches for tracking neural dynamics are of widespread interest, but a theoretical framework quantifying the physical limits of these techniques has been lacking. We formulate such a framework by using signal detection and estimation theory to obtain physical bounds on the detection of neural spikes and the estimation of their occurrence times as set by photon counting statistics (shot noise). These bounds are succinctly expressed via a discriminability index that depends on the kinetics of the optical indicator and the relative fluxes of signal and background photons. This approach facilitates quantitative evaluations of different indicators, detector technologies, and data analyses. Our treatment also provides optimal filtering techniques for optical detection of spikes. We compare various types of Ca(2+) indicators and show that background photons are a chief impediment to voltage sensing. Thus, voltage indicators that change color in response to membrane depolarization may offer a key advantage over those that change intensity. We also examine fluorescence resonance energy transfer indicators and identify the regimes in which the widely used ratiometric analysis of signals is substantially suboptimal. Overall, by showing how different optical factors interact to affect signal quality, our treatment offers a valuable guide to experimental design and provides measures of confidence to assess optically extracted traces of neural activity.

Making neurophysiological data analysis reproducible: why and how?

Reproducible data analysis is an approach aiming at complementing classical printed scientific articles with everything required to independently reproduce the results they present. "Everything" covers here: the data, the computer codes and a precise description of how the code was applied to the data. A brief history of this approach is presented first, starting with what economists have been calling replication since the early eighties to end with what is now called reproducible research in computational data analysis oriented fields like statistics and signal processing. Since efficient tools are instrumental for a routine implementation of these approaches, a description of some of the available ones is presented next. A toy example demonstrates then the use of two open source software programs for reproducible data analysis: the "Sweave family" and the org-mode of emacs. The former is bound to R while the latter can be used with R, Matlab, Python and many more "generalist" data processing software. Both solutions can be used with Unix-like, Windows and Mac families of operating systems. It is argued that neuroscientists could communicate much more efficiently their results by adopting the reproducible research paradigm from their lab books all the way to their articles, thesis and books.

Fast nonnegative deconvolution for spike train inference from population calcium imaging

Fluorescent calcium indicators are becoming increasingly popular as a means for observing the spiking activity of large neuronal populations. Unfortunately, extracting the spike train of each neuron from a raw fluorescence movie is a nontrivial problem. This work presents a fast nonnegative deconvolution filter to infer the approximately most likely spike train of each neuron, given the fluorescence observations. This algorithm outperforms optimal linear deconvolution (Wiener filtering) on both simulated and biological data. The performance gains come from restricting the inferred spike trains to be positive (using an interior-point method), unlike the Wiener filter. The algorithm runs in linear time, and is fast enough that even when simultaneously imaging >100 neurons, inference can be performed on the set of all observed traces faster than real time. Performing optimal spatial filtering on the images further refines the inferred spike train estimates. Importantly, all the parameters required to perform the inference can be estimated using only the fluorescence data, obviating the need to perform joint electrophysiological and imaging calibration experiments.