A preprocessing toolbox for 2-photon subcellular calcium imaging.

Recording the spiking activity from subcellular compartments of neurons such as axons and dendrites during behavior with 2-photon calcium imaging is increasingly common yet remains challenging due to low signal-to-noise, inaccurate region-of-interest (ROI) identification, movement artifacts, and difficulty in grouping ROIs from the same neuron. To address these issues, we present a computationally efficient pre-processing pipeline for subcellular signal detection, movement artifact identification, and ROI grouping. For subcellular signal detection, we capture the frequency profile of calcium transient dynamics by applying Fast Fourier Transform (FFT) on smoothed time-series calcium traces collected from axon ROIs. We then apply band-pass filtering methods (e.g. 0.05 to 0.12 Hz) to select ROIs that contain frequencies that match the power band of transients. To remove motion artifacts from z-plane movement, we apply Principal Component Analysis on all calcium traces and use a Bottom-Up Segmentation change-point detection model on the first principal component. After removing movement artifacts, we further identify calcium transients from noise by analyzing their prominence and duration. Finally, ROIs with high activity correlation are grouped using hierarchical or k-means clustering. Using axon ROIs in the CA1 region, we confirm that both clustering methods effectively determine the optimal number of clusters in pairwise correlation matrices, yielding similar groupings to "ground truth" data. Our approach provides a guideline for standardizing the extraction of physiological signals from subcellular compartments during behavior with 2-photon calcium imaging.