Calcium Imaging Analysis Guide

Calcium Imaging Analysis Guide

Purpose

This skill encodes expert methodological knowledge for analyzing calcium imaging data from fluorescent genetically encoded calcium indicators (GECIs). It covers the domain-specific decisions that a general-purpose programmer or data scientist would get wrong without specialized training in optical neurophysiology: choosing deconvolution parameters based on indicator kinetics, correcting neuropil contamination, handling modality-specific preprocessing, and interpreting fluorescence signals as neural activity.

When to Use This Skill

• Designing an analysis pipeline for two-photon, one-photon/miniscope, or fiber photometry data
• Choosing motion correction strategy for awake behaving animal recordings
• Selecting ROI detection method for your preparation density and imaging modality
• Setting deconvolution parameters matched to your calcium indicator
• Evaluating whether extracted signals reflect genuine single-neuron activity
• Troubleshooting common artifacts: photobleaching, neuropil contamination, motion

Research Planning Protocol

Before executing the domain-specific steps below, you MUST:

1. State the research question — What neural activity question is this calcium imaging analysis addressing?
2. Justify the method choice — Why calcium imaging (not electrophysiology, fMRI, etc.)? What alternatives were considered?
3. Declare expected outcomes — What activity patterns would support vs. refute the hypothesis?
4. Note assumptions and limitations — What does this method assume? Where could it mislead (e.g., indicator kinetics, neuropil)?
5. Present the plan to the user and WAIT for confirmation before proceeding.

For detailed methodology guidance, see the research-literacy skill.

⚠️ Verification Notice

This skill was generated by AI from academic literature. All parameters, thresholds, and citations require independent verification before use in research. If you find errors, please open an issue.

Decision Tree: Which Pipeline for Your Modality

What is your imaging modality?
 |
 +-- Two-photon (2P) microscopy
 | |
 | +-- Sparse labeling (Cre-dependent, cell-type specific)?
 | | --> Suite2P or CaImAn with standard CNMF
 | | Neuropil coefficient ~0.7 (Chen et al., 2013)
 | |
 | +-- Dense labeling (pan-neuronal)?
 | --> Suite2P with increased max_overlap (>=0.75)
 | or CaImAn CNMF with careful merge parameters
 | Consider soma-targeted indicators (Chen et al., 2020)
 |
 +-- One-photon (1P) / miniscope
 | |
 | --> CNMF-E (Zhou et al., 2018) or MIN1PIPE (Lu et al., 2018)
 | Standard CNMF will FAIL: 1P has large structured background
 | that requires explicit background modeling
 | CaImAn supports 1P via CNMF-E mode
 |
 +-- Fiber photometry (population-level)
 |
 --> No single-cell extraction needed
 Use isosbestic channel (405-415 nm) for motion/bleaching correction
 IRLS regression preferred over OLS (Lerner et al., 2015)
 Compute dF/F or z-scored signal

Core Pipeline Steps

Step 1: Motion Correction

Motion correction must precede all other analysis. Uncorrected motion creates false transients and blurs cellular signals.

Domain judgment:

• Always use non-rigid registration for awake behaving animals. Brain tissue deforms non-uniformly due to heartbeat, respiration, and locomotion (Dombeck et al., 2007).
• For 1P/miniscope data, apply spatial high-pass filtering before estimating motion vectors, because the large fluctuating background confounds rigid correlation-based registration (Giovannucci et al., 2019).
• Inspect registration quality by computing the correlation of each frame to the reference. Frames with correlation below 0.5 indicate severe motion and should be flagged or excluded (expert consensus).
• Edge pixels that shift out of the FOV during registration should be cropped; ROIs near edges may have unreliable signals.

Step 2: ROI Detection

Domain judgment:

• Set expected cell diameter to match your preparation. Cortical pyramidal somata are ~10-15 um; Purkinje cells are ~25 um; granule cells are ~5 um. Converting to pixels depends on your magnification and pixel size.
• For dense labeling, increase max_overlap to 0.75-1.0 (Suite2P) or adjust merge thresholds (CaImAn). Default overlap rejection discards valid overlapping neurons.
• Set connected=False when detecting dendrites or axonal boutons, which are not spatially contiguous in 2D projections.
• Always manually curate detected ROIs. Automated classifiers achieve near-human performance (Giovannucci et al., 2019) but are not perfect, especially for non-standard preparations.

Step 3: Neuropil Contamination Correction

Why this matters: In 2P imaging, each ROI's fluorescence contains signal from the surrounding neuropil (dense mesh of dendrites and axons). Without correction, you will observe artificial correlations between neurons, inflated response amplitudes, and obscured cell-specific tuning (Chen et al., 2013).

Correction formula:

F_corrected = F_raw - r * F_neuropil

Domain judgment:

• The coefficient r = 0.7 is an empirically derived average for 2P cortical imaging. It was estimated by measuring fluorescence in blood vessels (which should have zero true signal) relative to surrounding neuropil (Chen et al., 2013).
• Suite2P estimates r iteratively together with spike deconvolution, minimizing the deconvolution residual. This data-driven approach is preferred over a fixed coefficient.
• For 1P data, neuropil correction is handled differently: CNMF-E models the background as a low-rank spatiotemporal component rather than a per-ROI annulus (Zhou et al., 2018).
• Over-subtraction (r too high) produces negative fluorescence values. If you see many negative dF/F values, reduce r.
• For soma-targeted indicators (e.g., soma-GCaMP6f; Chen et al., 2020), neuropil contamination is reduced but not eliminated.

Step 4: dF/F Computation

Baseline estimation methods:

Formula:

dF/F = (F(t) - F0) / F0

Domain judgment:

• The sliding window for rolling percentile should be 30-60 seconds to capture slow baseline fluctuations without tracking fast transients (expert consensus).
• Using the mean as baseline inflates dF/F for highly active neurons. Use a low percentile (8th-10th) or the mode instead (Dombeck et al., 2007).
• Photobleaching causes an exponential decay in baseline fluorescence. If uncorrected, early time points will have artificially low dF/F and late time points artificially high dF/F. Fit and subtract an exponential before computing dF/F (Giovannucci et al., 2019).

Step 5: Deconvolution / Spike Inference

Deconvolution estimates the underlying spike train from the slow calcium fluorescence signal.

Critical: The deconvolution kernel decay constant (tau) must match your calcium indicator. See references/indicator-parameters.md for the full table.

Domain judgment:

• If tau is wrong, deconvolution is systematically biased. Too long produces sluggish traces; too short produces overshoot artifacts.
• The kernel is fixed to tau in Suite2P; it is not fit to the data. Always verify your indicator and set tau explicitly.
• Do NOT threshold deconvolved events to create binary spike trains. This discards magnitude information: a 1-spike event and a 10-spike burst become identical (Pachitariu et al., 2017).
• Deconvolved amplitudes are in arbitrary units. There is an unknown scaling factor between fluorescence and spike count. Do not interpret amplitudes as absolute firing rates.
• CASCADE (Rupprecht et al., 2021) can produce calibrated spike rate estimates by resampling ground truth data to match your recording's noise level, but still has systematic errors at low SNR.