Multi-scale approaches for high-speed imaging and analysis of large neural populations

Optical constraints on two-photon voltage imaging

Significance

Genetically encoded voltage indicators (GEVIs) are a valuable tool for studying neural circuits in vivo, but the relative merits and limitations of one-photon (1P) versus two-photon (2P) voltage imaging are not well characterized.

Aim

We consider the optical and biophysical constraints particular to 1P and 2P voltage imaging and compare the imaging properties of commonly used GEVIs under 1P and 2P excitation.

Approach

We measure the brightness and voltage sensitivity of voltage indicators from commonly used classes under 1P and 2P illumination. We also measure the decrease in fluorescence as a function of depth in the mouse brain. We develop a simple model of the number of measurable cells as a function of reporter properties, imaging parameters, and desired signal-to-noise ratio (SNR). We then discuss how the performance of voltage imaging would be affected by sensor improvements and by recently introduced advanced imaging modalities.

Results

Compared with 1P excitation, 2P excitation requires ∼ 10 4 -fold more illumination power per cell to produce similar photon count rates. For voltage imaging with JEDI-2P in the mouse cortex with a target SNR of 10 (spike height to baseline shot noise), a measurement bandwidth of 1 kHz, a thermally limited laser power of 200 mW, and an imaging depth of > 300    μ m , 2P voltage imaging using an 80-MHz source can record from no more than ∼ 12 neurons simultaneously.

Conclusions

Due to the stringent photon-count requirements of voltage imaging and the modest voltage sensitivity of existing reporters, 2P voltage imaging in vivo faces a stringent tradeoff between shot noise and tissue photodamage. 2P imaging of hundreds of neurons with high SNR at a depth of > 300    μ m will require either major improvements in 2P GEVIs or qualitatively new approaches to imaging.

Optical constraints on two-photon voltage imaging

Significance

Genetically encoded voltage indicators (GEVIs) are a valuable tool for studying neural circuits in vivo, but the relative merits and limitations of one-photon (1P) vs. two-photon (2P) voltage imaging are not well characterized.

Aim

We consider the optical and biophysical constraints particular to 1P and 2P voltage imaging and compare the imaging properties of commonly used GEVIs under 1P and 2P excitation.

Approach

We measure brightness and voltage sensitivity of voltage indicators from commonly used classes under 1P and 2P illumination. We also measure the decrease in fluorescence as a function of depth in mouse brain. We develop a simple model of the number of measurable cells as a function of reporter properties, imaging parameters, and desired signal-to-noise ratio (SNR). We then discuss how the performance of voltage imaging would be affected by sensor improvements and by recently introduced advanced imaging modalities.

Results

Compared to 1P excitation, 2P excitation requires ~104-fold more illumination power per cell to produce similar photon count rates. For voltage imaging with JEDI-2P in mouse cortex with a target SNR of 10 (spike height:baseline shot noise), a measurement bandwidth of 1 kHz, a thermally limited laser power of 200 mW, and an imaging depth of > 300 μm, 2P voltage imaging using an 80 MHz source can record from no more 12 cells simultaneously.

Conclusions

Due to the stringent photon-count requirements of voltage imaging and the modest voltage sensitivity of existing reporters, 2P voltage imaging in vivo faces a stringent tradeoff between shot noise and tissue photodamage. 2P imaging of hundreds of neurons with high SNR at depth > 300 μm will require either major improvements in 2P GEVIs or qualitatively new approaches to imaging.

Barcoding of episodic memories in the hippocampus of a food-caching bird

The hippocampus is critical for episodic memory. Although hippocampal activity represents place and other behaviorally relevant variables, it is unclear how it encodes numerous memories of specific events in life. To study episodic coding, we leveraged the specialized behavior of chickadees-food-caching birds that form memories at well-defined moments in time whenever they cache food for subsequent retrieval. Our recordings during caching revealed very sparse, transient barcode-like patterns of firing across hippocampal neurons. Each "barcode" uniquely represented a caching event and transiently reactivated during the retrieval of that specific cache. Barcodes co-occurred with the conventional activity of place cells but were uncorrelated even for nearby cache locations that had similar place codes. We propose that animals recall episodic memories by reactivating hippocampal barcodes. Similarly to computer hash codes, these patterns assign unique identifiers to different events and could be a mechanism for rapid formation and storage of many non-interfering memories.

maskNMF: A denoise-sparsen-detect approach for extracting neural signals from dense imaging data

A number of calcium imaging methods have been developed to monitor the activity of large populations of neurons. One particularly promising approach, Bessel imaging, captures neural activity from a volume by projecting within the imaged volume onto a single imaging plane, therefore effectively mixing signals and increasing the number of neurons imaged per pixel. These signals must then be computationally demixed to recover the desired neural activity. Unfortunately, currently-available demixing methods can perform poorly in the regime of high imaging density (i.e., many neurons per pixel). In this work we introduce a new pipeline (maskNMF) for demixing dense calcium imaging data. The main idea is to first denoise and temporally sparsen the observed video; this enhances signal strength and reduces spatial overlap significantly. Next we detect neurons in the sparsened video using a neural network trained on a library of neural shapes. These shapes are derived from segmented electron microscopy images input into a Bessel imaging model; therefore no manual selection of "good" neural shapes from the functional data is required here. After cells are detected, we use a constrained non-negative matrix factorization approach to demix the activity, using the detected cells' shapes to initialize the factorization. We test the resulting pipeline on both simulated and real datasets and find that it is able to achieve accurate demixing on denser data than was previously feasible, therefore enabling faithful imaging of larger neural populations. The method also provides good results on more "standard" two-photon imaging data. Finally, because much of the pipeline operates on a significantly compressed version of the raw data and is highly parallelizable, the algorithm is fast, processing large datasets faster than real time.

Barcoding of episodic memories in the hippocampus of a food-caching bird

Episodic memory, or memory of experienced events, is a critical function of the hippocampus1-3. It is therefore important to understand how hippocampal activity represents specific events in an animal's life. We addressed this question in chickadees - specialist food-caching birds that hide food at scattered locations and use memory to find their caches later in time4,5. We performed high-density neural recordings in the hippocampus of chickadees as they cached and retrieved seeds in a laboratory arena. We found that each caching event was represented by a burst of firing in a unique set of hippocampal neurons. These 'barcode-like' patterns of activity were sparse (<10% of neurons active), uncorrelated even for immediately adjacent caches, and different even for separate caches at the same location. The barcode representing a specific caching event was transiently reactivated whenever a bird later interacted with the same cache - for example, to retrieve food. Barcodes co-occurred with conventional place cell activity6,7, as well as location-independent responses to cached seeds. We propose that barcodes are signatures of episodic memories evoked during memory recall. These patterns assign a unique identifier to each event and may be a mechanism for rapid formation and storage of many non-interfering memories.

An entorhinal-like region in food-caching birds

The mammalian entorhinal cortex routes inputs from diverse sources into the hippocampus. This information is mixed and expressed in the activity of many specialized entorhinal cell types, which are considered indispensable for hippocampal function. However, functionally similar hippocampi exist even in non-mammals that lack an obvious entorhinal cortex or, generally, any layered cortex. To address this dilemma, we mapped extrinsic hippocampal connections in chickadees, whose hippocampi are used for remembering numerous food caches. We found a well-delineated structure in these birds that is topologically similar to the entorhinal cortex and interfaces between the hippocampus and other pallial regions. Recordings of this structure revealed entorhinal-like activity, including border and multi-field grid-like cells. These cells were localized to the subregion predicted by anatomical mapping to match the dorsomedial entorhinal cortex. Our findings uncover an anatomical and physiological equivalence of vastly different brains, suggesting a fundamental nature of entorhinal-like computations for hippocampal function.

Glioma epileptiform activity and progression are driven by IGSF3-mediated potassium dysregulation

Seizures are a frequent pathophysiological feature of malignant glioma. Recent studies implicate peritumoral synaptic dysregulation as a driver of brain hyperactivity and tumor progression; however, the molecular mechanisms that govern these phenomena remain elusive. Using scRNA-seq and intraoperative patient ECoG recordings, we show that tumors from seizure patients are enriched for gene signatures regulating synapse formation. Employing a human-to-mouse in vivo functionalization pipeline to screen these genes, we identify IGSF3 as a mediator of glioma progression and dysregulated neural circuitry that manifests as spreading depolarization (SD). Mechanistically, we discover that IGSF3 interacts with Kir4.1 to suppress potassium buffering and found that seizure patients exhibit reduced expression of potassium handlers in proliferating tumor cells. In vivo imaging reveals that dysregulated synaptic activity emanates from the tumor-neuron interface, which we confirm in patients. Our studies reveal that tumor progression and seizures are enabled by ion dyshomeostasis and identify SD as a driver of disease.

An entorhinal-like region in food-caching birds

The mammalian entorhinal cortex routes inputs from diverse sources into the hippocampus. This information is mixed and expressed in the activity of many specialized entorhinal cell types, which are considered indispensable for hippocampal function. However, functionally similar hippocampi exist even in non-mammals that lack an obvious entorhinal cortex, or generally any layered cortex. To address this dilemma, we mapped extrinsic hippocampal connections in chickadees, whose hippocampi are used for remembering numerous food caches. We found a well-delineated structure in these birds that is topologically similar to the entorhinal cortex and interfaces between the hippocampus and other pallial regions. Recordings of this structure revealed entorhinal-like activity, including border and multi-field grid-like cells. These cells were localized to the subregion predicted by anatomical mapping to match the dorsomedial entorhinal cortex. Our findings uncover an anatomical and physiological equivalence of vastly different brains, suggesting a fundamental nature of entorhinal-like computations for hippocampal function.

Review of data processing of functional optical microscopy for neuroscience

Functional optical imaging in neuroscience is rapidly growing with the development of optical systems and fluorescence indicators. To realize the potential of these massive spatiotemporal datasets for relating neuronal activity to behavior and stimuli and uncovering local circuits in the brain, accurate automated processing is increasingly essential. We cover recent computational developments in the full data processing pipeline of functional optical microscopy for neuroscience data and discuss ongoing and emerging challenges.

An improved platform for cultured neuronal network electrophysiology: multichannel optogenetics integrated with MEAs

Cultured neuronal networks (CNNs) are powerful tools for studying how neuronal representation and adaptation emerge in networks of controlled populations of neurons. To ensure the interaction of a CNN and an artificial setting, reliable operation in both open and closed loops should be provided. In this study, we integrated optogenetic stimulation with microelectrode array (MEA) recordings using a digital micromirror device and developed an improved research tool with a 64-channel interface for neuronal network control and data acquisition. We determined the ideal stimulation parameters including light intensity, frequency, and duty cycle for our configuration. This resulted in robust and reproducible neuronal responses. We also demonstrated both open and closed loop configurations in the new platform involving multiple bidirectional channels. Unlike previous approaches that combined optogenetic stimulation and MEA recordings, we did not use binary grid patterns, but assigned an adjustable-size, non-binary optical spot to each electrode. This approach allowed simultaneous use of multiple input-output channels and facilitated adaptation of the stimulation parameters. Hence, we advanced a 64-channel interface in that each channel can be controlled individually in both directions simultaneously without any interference or interrupts. The presented setup meets the requirements of research in neuronal plasticity, network encoding and representation, closed-loop control of firing rate and synchronization. Researchers who develop closed-loop control techniques and adaptive stimulation strategies for network activity will benefit much from this novel setup.

Assessing Local and Branch-specific Activity in Dendrites

Dendrites are elaborate neural processes which integrate inputs from various sources in space and time. While decades of work have suggested an independent role for dendrites in driving nonlinear computations for the cell, only recently have technological advances enabled us to capture the variety of activity in dendrites and their coupling dynamics with the soma. Under certain circumstances, activity generated in a given dendritic branch remains isolated, such that the soma or even sister dendrites are not privy to these localized signals. Such branch-specific activity could radically increase the capacity and flexibility of coding for the cell as a whole. Here, we discuss these forms of localized and branch-specific activity, their functional relevance in plasticity and behavior, and their supporting biophysical and circuit-level mechanisms. We conclude by showcasing electrical and optical approaches in hippocampal area CA3, using original experimental data to discuss experimental and analytical methodology and key considerations to take when investigating the functional relevance of independent dendritic activity.

High-fidelity estimates of spikes and subthreshold waveforms from 1-photon voltage imaging in vivo

The ability to probe the membrane potential of multiple genetically defined neurons simultaneously would have a profound impact on neuroscience research. Genetically encoded voltage indicators are a promising tool for this purpose, and recent developments have achieved a high signal-to-noise ratio in vivo with 1-photon fluorescence imaging. However, these recordings exhibit several sources of noise and signal extraction remains a challenge. We present an improved signal extraction pipeline, spike-guided penalized matrix decomposition-nonnegative matrix factorization (SGPMD-NMF), which resolves supra- and subthreshold voltages in vivo. The method incorporates biophysical and optical constraints. We validate the pipeline with simultaneous patch-clamp and optical recordings from mouse layer 1 in vivo and with simulated and composite datasets with realistic noise. We demonstrate applications to mouse hippocampus expressing paQuasAr3-s or SomArchon1, mouse cortex expressing SomArchon1 or Voltron, and zebrafish spines expressing zArchon1.

Place Cells in Head-Fixed Mice Navigating a Floating Real-World Environment

The hippocampal place cell system in rodents has provided a major paradigm for the scientific investigation of memory function and dysfunction. Place cells have been observed in area CA1 of the hippocampus of both freely moving animals, and of head-fixed animals navigating in virtual reality environments. However, spatial coding in virtual reality preparations has been observed to be impaired. Here we show that the use of a real-world environment system for head-fixed mice, consisting of an air-floating track with proximal cues, provides some advantages over virtual reality systems for the study of spatial memory. We imaged the hippocampus of head-fixed mice injected with the genetically encoded calcium indicator GCaMP6s while they navigated circularly constrained or open environments on the floating platform. We observed consistent place tuning in a substantial fraction of cells despite the absence of distal visual cues. Place fields remapped when animals entered a different environment. When animals re-entered the same environment, place fields typically remapped over a time period of multiple days, faster than in freely moving preparations, but comparable with virtual reality. Spatial information rates were within the range observed in freely moving mice. Manifold analysis indicated that spatial information could be extracted from a low-dimensional subspace of the neural population dynamics. This is the first demonstration of place cells in head-fixed mice navigating on an air-lifted real-world platform, validating its use for the study of brain circuits involved in memory and affected by neurodegenerative disorders.

Online analysis of microendoscopic 1-photon calcium imaging data streams

In vivo calcium imaging through microendoscopic lenses enables imaging of neuronal populations deep within the brains of freely moving animals. Previously, a constrained matrix factorization approach (CNMF-E) has been suggested to extract single-neuronal activity from microendoscopic data. However, this approach relies on offline batch processing of the entire video data and is demanding both in terms of computing and memory requirements. These drawbacks prevent its applicability to the analysis of large datasets and closed-loop experimental settings. Here we address both issues by introducing two different online algorithms for extracting neuronal activity from streaming microendoscopic data. Our first algorithm, OnACID-E, presents an online adaptation of the CNMF-E algorithm, which dramatically reduces its memory and computation requirements. Our second algorithm proposes a convolution-based background model for microendoscopic data that enables even faster (real time) processing. Our approach is modular and can be combined with existing online motion artifact correction and activity deconvolution methods to provide a highly scalable pipeline for microendoscopic data analysis. We apply our algorithms on four previously published typical experimental datasets and show that they yield similar high-quality results as the popular offline approach, but outperform it with regard to computing time and memory requirements. They can be used instead of CNMF-E to process pre-recorded data with boosted speeds and dramatically reduced memory requirements. Further, they newly enable online analysis of live-streaming data even on a laptop.

Calcium imaging methods enable researchers to measure the activity of genetically-targeted large-scale neuronal subpopulations. Whereas previous methods required the specimen to be stable, e.g. anesthetized or head-fixed, new brain imaging techniques using microendoscopic lenses and miniaturized microscopes have enabled deep brain imaging in freely moving mice. However, the very large background fluctuations, the inevitable movements and distortions of imaging field, and the extensive spatial overlaps of fluorescent signals complicate the goal of efficiently extracting accurate estimates of neural activity from the observed video data. Further, current activity extraction methods are computationally expensive due to the complex background model and are typically applied to imaging data long after the experiment is complete. Moreover, in some scenarios it is necessary to perform experiments in real-time and closed-loop—analyzing data on-the-fly to guide the next experimental steps or to control feedback –, and this calls for new methods for accurate real-time processing. Here we address both issues by adapting a popular extraction method to operate online and extend it to utilize GPU hardware that enables real time processing. Our algorithms yield similar high-quality results as the original offline approach, but outperform it with regard to computing time and memory requirements. Our results enable faster and scalable analysis, and open the door to new closed-loop experiments in deep brain areas and on freely-moving preparations. Our algorithms can be used for newly enabled real-time analysis of streaming data, as well as swapped in directly to replace the computationally costly offline approach.

Localized semi-nonnegative matrix factorization (LocaNMF) of widefield calcium imaging data

Widefield calcium imaging enables recording of large-scale neural activity across the mouse dorsal cortex. In order to examine the relationship of these neural signals to the resulting behavior, it is critical to demix the recordings into meaningful spatial and temporal components that can be mapped onto well-defined brain regions. However, no current tools satisfactorily extract the activity of the different brain regions in individual mice in a data-driven manner, while taking into account mouse-specific and preparation-specific differences. Here, we introduce Localized semi-Nonnegative Matrix Factorization (LocaNMF), a method that efficiently decomposes widefield video data and allows us to directly compare activity across multiple mice by outputting mouse-specific localized functional regions that are significantly more interpretable than more traditional decomposition techniques. Moreover, it provides a natural subspace to directly compare correlation maps and neural dynamics across different behaviors, mice, and experimental conditions, and enables identification of task- and movement-related brain regions.

Fast and robust active neuron segmentation in two-photon calcium imaging using spatiotemporal deep learning

Calcium imaging records large-scale neuronal activity with cellular resolution in vivo. Automated, fast, and reliable active neuron segmentation is a critical step in the analysis workflow of utilizing neuronal signals in real-time behavioral studies for discovery of neuronal coding properties. Here, to exploit the full spatiotemporal information in two-photon calcium imaging movies, we propose a 3D convolutional neural network to identify and segment active neurons. By utilizing a variety of two-photon microscopy datasets, we show that our method outperforms state-of-the-art techniques and is on a par with manual segmentation. Furthermore, we demonstrate that the network trained on data recorded at a specific cortical layer can be used to accurately segment active neurons from another layer with different neuron density. Finally, our work documents significant tabulation flaws in one of the most cited and active online scientific challenges in neuron segmentation. As our computationally fast method is an invaluable tool for a large spectrum of real-time optogenetic experiments, we have made our open-source software and carefully annotated dataset freely available online.

Single-neuron perturbations reveal feature-specific competition in V1

The computations performed by local neural populations, such as a cortical layer, are typically inferred from anatomical connectivity and observations of neural activity. Here we describe a method-influence mapping-that uses single-neuron perturbations to directly measure how cortical neurons reshape sensory representations. In layer 2/3 of the primary visual cortex (V1), we use two-photon optogenetics to trigger action potentials in a targeted neuron and calcium imaging to measure the effect on spiking in neighbouring neurons in awake mice viewing visual stimuli. Excitatory neurons on average suppressed other neurons and had a centre-surround influence profile over anatomical space. A neuron's influence on its neighbour depended on their similarity in activity. Notably, neurons suppressed activity in similarly tuned neurons more than in dissimilarly tuned neurons. In addition, photostimulation reduced the population response, specifically to the targeted neuron's preferred stimulus, by around 2%. Therefore, V1 layer 2/3 performed feature competition, in which a like-suppresses-like motif reduces redundancy in population activity and may assist with inference of the features that underlie sensory input. We anticipate that influence mapping can be extended to investigate computations in other neural populations.

Analysis pipelines for calcium imaging data

Calcium imaging is a popular tool among neuroscientists because of its capability to monitor in vivo large neural populations across weeks with single neuron and single spike resolution. Before any downstream analysis, the data needs to be pre-processed to extract the location and activity of the neurons and processes in the observed field of view. The ever increasing size of calcium imaging datasets necessitates scalable analysis pipelines that are reproducible and fully automated. This review focuses on recent methods for addressing the pre-processing problems that arise in calcium imaging data analysis, and available software tools for high throughput analysis pipelines.

Neural data science: accelerating the experiment-analysis-theory cycle in large-scale neuroscience

Modern large-scale multineuronal recording methodologies, including multielectrode arrays, calcium imaging, and optogenetic techniques, produce single-neuron resolution data of a magnitude and precision that were the realm of science fiction twenty years ago. The major bottlenecks in systems and circuit neuroscience no longer lie in simply collecting data from large neural populations, but also in understanding this data: developing novel scientific questions, with corresponding analysis techniques and experimental designs to fully harness these new capabilities and meaningfully interrogate these questions. Advances in methods for signal processing, network analysis, dimensionality reduction, and optimal control-developed in lockstep with advances in experimental neurotechnology-promise major breakthroughs in multiple fundamental neuroscience problems. These trends are clear in a broad array of subfields of modern neuroscience; this review focuses on recent advances in methods for analyzing neural time-series data with single-neuronal precision.

CaImAn an open source tool for scalable calcium imaging data analysis

Advances in fluorescence microscopy enable monitoring larger brain areas in-vivo with finer time resolution. The resulting data rates require reproducible analysis pipelines that are reliable, fully automated, and scalable to datasets generated over the course of months. We present CaImAn, an open-source library for calcium imaging data analysis. CaImAn provides automatic and scalable methods to address problems common to pre-processing, including motion correction, neural activity identification, and registration across different sessions of data collection. It does this while requiring minimal user intervention, with good scalability on computers ranging from laptops to high-performance computing clusters. CaImAn is suitable for two-photon and one-photon imaging, and also enables real-time analysis on streaming data. To benchmark the performance of CaImAn we collected and combined a corpus of manual annotations from multiple labelers on nine mouse two-photon datasets. We demonstrate that CaImAn achieves near-human performance in detecting locations of active neurons.

A Guide to Emerging Technologies for Large-Scale and Whole-Brain Optical Imaging of Neuronal Activity

The mammalian brain is a densely interconnected network that consists of millions to billions of neurons. Decoding how information is represented and processed by this neural circuitry requires the ability to capture and manipulate the dynamics of large populations at high speed and high resolution over a large area of the brain. Although the use of optical approaches by the neuroscience community has rapidly increased over the past two decades, most microscopy approaches are unable to record the activity of all neurons comprising a functional network across the mammalian brain at relevant temporal and spatial resolutions. In this review, we survey the recent development in optical technologies for Ca2+ imaging in this regard and provide an overview of the strengths and limitations of each modality and its potential for scalability. We provide guidance from the perspective of a biological user driven by the typical biological applications and sample conditions. We also discuss the potential for future advances and synergies that could be obtained through hybrid approaches or other modalities.

Neural Coding With Bursts-Current State and Future Perspectives

Neuronal action potentials or spikes provide a long-range, noise-resistant means of communication between neurons. As point processes single spikes contain little information in themselves, i.e., outside the context of spikes from other neurons. Moreover, they may fail to cross a synapse. A burst, which consists of a short, high frequency train of spikes, will more reliably cross a synapse, increasing the likelihood of eliciting a postsynaptic spike, depending on the specific short-term plasticity at that synapse. Both the number and the temporal pattern of spikes in a burst provide a coding space that lies within the temporal integration realm of single neurons. Bursts have been observed in many species, including the non-mammalian, and in brain regions that range from subcortical to cortical. Despite their widespread presence and potential relevance, the uncertainties of how to classify bursts seems to have limited the research into the coding possibilities for bursts. The present series of research articles provides new insights into the relevance and interpretation of bursts across different neural circuits, and new methods for their analysis. Here, we provide a succinct introduction to the history of burst coding and an overview of recent work on this topic.

Efficient and accurate extraction of in vivo calcium signals from microendoscopic video data

In vivo calcium imaging through microendoscopic lenses enables imaging of previously inaccessible neuronal populations deep within the brains of freely moving animals. However, it is computationally challenging to extract single-neuronal activity from microendoscopic data, because of the very large background fluctuations and high spatial overlaps intrinsic to this recording modality. Here, we describe a new constrained matrix factorization approach to accurately separate the background and then demix and denoise the neuronal signals of interest. We compared the proposed method against previous independent components analysis and constrained nonnegative matrix factorization approaches. On both simulated and experimental data recorded from mice, our method substantially improved the quality of extracted cellular signals and detected more well-isolated neural signals, especially in noisy data regimes. These advances can in turn significantly enhance the statistical power of downstream analyses, and ultimately improve scientific conclusions derived from microendoscopic data.

FISSA: A neuropil decontamination toolbox for calcium imaging signals

In vivo calcium imaging has become a method of choice to image neuronal population activity throughout the nervous system. These experiments generate large sequences of images. Their analysis is computationally intensive and typically involves motion correction, image segmentation into regions of interest (ROIs), and extraction of fluorescence traces from each ROI. Out of focus fluorescence from surrounding neuropil and other cells can strongly contaminate the signal assigned to a given ROI. In this study, we introduce the FISSA toolbox (Fast Image Signal Separation Analysis) for neuropil decontamination. Given pre-defined ROIs, the FISSA toolbox automatically extracts the surrounding local neuropil and performs blind-source separation with non-negative matrix factorization. Using both simulated and in vivo data, we show that this toolbox performs similarly or better than existing published methods. FISSA requires only little RAM, and allows for fast processing of large datasets even on a standard laptop. The FISSA toolbox is available in Python, with an option for MATLAB format outputs, and can easily be integrated into existing workflows. It is available from Github and the standard Python repositories.

In vivo imaging of neural activity

Since the introduction of calcium imaging to monitor neuronal activity with single-cell resolution, optical imaging methods have revolutionized neuroscience by enabling systematic recordings of neuronal circuits in living animals. The plethora of methods for functional neural imaging can be daunting to the nonexpert to navigate. Here we review advanced microscopy techniques for in vivo functional imaging and offer guidelines for which technologies are best suited for particular applications.