Computational processing of neural recordings from calcium imaging data

Functional ultrasound imaging of the human spinal cord

Utilizing the first in-human functional ultrasound imaging (fUSI) of the spinal cord, we demonstrate the integration of spinal functional responses to electrical stimulation. We record and characterize the hemodynamic responses of the spinal cord to a neuromodulatory intervention commonly used for treating pain and increasingly used for the restoration of sensorimotor and autonomic function. We found that the hemodynamic response to stimulation reflects a spatiotemporal modulation of the spinal cord circuitry not previously recognized. Our analytical capability offers a mechanism to assess blood flow changes with a new level of spatial and temporal precision in vivo and demonstrates that fUSI can decode the functional state of spinal networks in a single trial, which is of fundamental importance for developing real-time closed-loop neuromodulation systems. This work is a critical step toward developing a vital technique to study spinal cord function and effects of clinical neuromodulation.

Learning enhances representations of taste-guided decisions in the mouse gustatory insular cortex

Learning to discriminate overlapping gustatory stimuli that predict distinct outcomes-a feat known as discrimination learning-can mean the difference between ingesting a poison or a nutritive meal. Despite the obvious importance of this process, very little is known about the neural basis of taste discrimination learning. In other sensory modalities, this form of learning can be mediated by either the sharpening of sensory representations or the enhanced ability of "decision-making" circuits to interpret sensory information. Given the dual role of the gustatory insular cortex (GC) in encoding both sensory and decision-related variables, this region represents an ideal site for investigating how neural activity changes as animals learn a novel taste discrimination. Here, we present results from experiments relying on two-photon calcium imaging of GC neural activity in mice performing a taste-guided mixture discrimination task. The task allows for the recording of neural activity before and after learning induced by training mice to discriminate increasingly similar pairs of taste mixtures. Single-neuron and population analyses show a time-varying pattern of activity, with early sensory responses emerging after taste delivery and binary, choice-encoding responses emerging later in the delay before a decision is made. Our results demonstrate that, while both sensory and decision-related information is encoded by GC in the context of a taste mixture discrimination task, learning and improved performance are associated with a specific enhancement of decision-related responses.

NeuroSeg-III: efficient neuron segmentation in two-photon Ca2+ imaging data using self-supervised learning

Two-photon Ca2+ imaging technology increasingly plays an essential role in neuroscience research. However, the requirement for extensive professional annotation poses a significant challenge to improving the performance of neuron segmentation models. Here, we present NeuroSeg-III, an innovative self-supervised learning approach specifically designed to achieve fast and precise segmentation of neurons in imaging data. This approach consists of two modules: a self-supervised pre-training network and a segmentation network. After pre-training the encoder of the segmentation network via a self-supervised learning method without any annotated data, we only need to fine-tune the segmentation network with a small amount of annotated data. The segmentation network is designed with YOLOv8s, FasterNet, efficient multi-scale attention mechanism (EMA), and bi-directional feature pyramid network (BiFPN), which enhanced the model's segmentation accuracy while reducing the computational cost and parameters. The generalization of our approach was validated across different Ca2+ indicators and scales of imaging data. Significantly, the proposed neuron segmentation approach exhibits exceptional speed and accuracy, surpassing the current state-of-the-art benchmarks when evaluated using a publicly available dataset. The results underscore the effectiveness of NeuroSeg-III, with employing an efficient training strategy tailored for two-photon Ca2+ imaging data and delivering remarkable precision in neuron segmentation.

Theta-frequency medial septal nucleus deep brain stimulation increases neurovascular activity in MK-801-treated mice

Introduction

Deep brain stimulation (DBS) has shown remarkable success treating neurological and psychiatric disorders including Parkinson's disease, essential tremor, dystonia, epilepsy, and obsessive-compulsive disorder. DBS is now being explored to improve cognitive and functional outcomes in other psychiatric conditions, such as those characterized by reduced N-methyl-D-aspartate (NMDA) function (i.e., schizophrenia). While DBS for movement disorders generally involves high-frequency (>100 Hz) stimulation, there is evidence that low-frequency stimulation may have beneficial and persisting effects when applied to cognitive brain networks.

Methods

In this study, we utilize a novel technology, functional ultrasound imaging (fUSI), to characterize the cerebrovascular impact of medial septal nucleus (MSN) DBS under conditions of NMDA antagonism (pharmacologically using Dizocilpine [MK-801]) in anesthetized male mice.

Results

Imaging from a sagittal plane across a variety of brain regions within and outside of the septohippocampal circuit, we find that MSN theta-frequency (7.7 Hz) DBS increases hippocampal cerebral blood volume (CBV) during and after stimulation. This effect was not present using standard high-frequency stimulation parameters [i.e., gamma (100 Hz)].

Discussion

These results indicate the MSN DBS increases circuit-specific hippocampal neurovascular activity in a frequency-dependent manner and does so in a way that continues beyond the period of electrical stimulation.

Memory encoding and retrieval by retrosplenial parvalbumin interneurons are impaired in Alzheimer's disease model mice

Memory deficits in Alzheimer's disease (AD) show a strong link with GABAergic interneuron dysfunctions.1,2,3,4,5,6,7 The ensemble dynamics of GABAergic interneurons represent memory encoding and retrieval,8,9,10,11,12 but how GABAergic interneuron dysfunction affects inhibitory ensemble dynamics in AD is unknown. As the retrosplenial cortex (RSC) is critical for episodic memory13,14,15,16 and is affected by β-amyloid accumulation in early AD,17,18,19,20,21 we address this question by performing Ca2+ imaging in RSC parvalbumin (PV)-expressing interneurons during a contextual fear memory task in healthy control mice and the 5XFAD mouse model of AD. We found that populations of PV interneurons responsive to aversive electric foot shocks during contextual fear conditioning (shock-responsive) significantly decreased in the 5XFAD mice, indicating dysfunctions in the recruitment of memory-encoding PV interneurons. In the control mice, ensemble activities of shock-responsive PV interneurons were selectively upregulated during the freezing epoch of the contextual fear memory retrieval, manifested by synaptic potentiation of PV interneuron-mediated inhibition. However, such changes in ensemble dynamics during memory retrieval and synaptic plasticity were both absent in the 5XFAD mice. Optogenetic silencing of PV interneurons during contextual fear conditioning in the control mice mimicked the memory deficits in the 5XFAD mice, while optogenetic activation of PV interneurons in the 5XFAD mice restored memory retrieval. These results demonstrate the critical roles of contextual fear memory-encoding PV interneurons for memory retrieval. Furthermore, synaptic dysfunction of PV interneurons may disrupt the recruitment of PV interneurons and their ensemble dynamics underlying contextual fear memory retrieval, subsequently leading to memory deficits in AD.

Modeling the diverse effects of divisive normalization on noise correlations

Divisive normalization, a prominent descriptive model of neural activity, is employed by theories of neural coding across many different brain areas. Yet, the relationship between normalization and the statistics of neural responses beyond single neurons remains largely unexplored. Here we focus on noise correlations, a widely studied pairwise statistic, because its stimulus and state dependence plays a central role in neural coding. Existing models of covariability typically ignore normalization despite empirical evidence suggesting it affects correlation structure in neural populations. We therefore propose a pairwise stochastic divisive normalization model that accounts for the effects of normalization and other factors on covariability. We first show that normalization modulates noise correlations in qualitatively different ways depending on whether normalization is shared between neurons, and we discuss how to infer when normalization signals are shared. We then apply our model to calcium imaging data from mouse primary visual cortex (V1), and find that it accurately fits the data, often outperforming a popular alternative model of correlations. Our analysis indicates that normalization signals are often shared between V1 neurons in this dataset. Our model will enable quantifying the relation between normalization and covariability in a broad range of neural systems, which could provide new constraints on circuit mechanisms of normalization and their role in information transmission and representation.

Advances in In Vitro Models of Neuromuscular Junction: Focusing on Organ-on-a-Chip, Organoids, and Biohybrid Robotics

The neuromuscular junction (NMJ) is a peripheral synaptic connection between presynaptic motor neurons and postsynaptic skeletal muscle fibers that enables muscle contraction and voluntary motor movement. Many traumatic, neurodegenerative, and neuroimmunological diseases are classically believed to mainly affect either the neuronal or the muscle side of the NMJ, and treatment options are lacking. Recent advances in novel techniques have helped develop in vitro physiological and pathophysiological models of the NMJ as well as enable precise control and evaluation of its functions. This paper reviews the recent developments in in vitro NMJ models with 2D or 3D cultures, from organ-on-a-chip and organoids to biohybrid robotics. Related derivative techniques are introduced for functional analysis of the NMJ, such as the patch-clamp technique, microelectrode arrays, calcium imaging, and stimulus methods, particularly optogenetic-mediated light stimulation, microelectrode-mediated electrical stimulation, and biochemical stimulation. Finally, the applications of the in vitro NMJ models as disease models or for drug screening related to suitable neuromuscular diseases are summarized and their future development trends and challenges are discussed.

Dissociable roles of thalamic nuclei in the refinement of reaches to spatial targets

Reaches are complex movements that are critical for survival, and encompass the control of different aspects such as direction, speed, and endpoint precision. Complex movements have been postulated to be learned and controlled through distributed motor networks, of which the thalamus is a highly connected node. Still, the role of different thalamic circuits in learning and controlling specific aspects of reaches has not been investigated. We report dissociable roles of two distinct thalamic nuclei - the parafascicular (Pf) and ventroanterior/ventrolateral (VAL) nuclei - in the refinement of spatial target reaches in mice. Using 2-photon calcium imaging in a head-fixed joystick task where mice learned to reach to a target in space, we found that glutamatergic neurons in both areas were most active during reaches early in learning. Reach-related activity in both areas decreased late in learning, as movement direction was refined and reaches increased in accuracy. Furthermore, the population dynamics of Pf, but not VAL, covaried in different subspaces in early and late learning, but eventually stabilized in late learning. The neural activity in Pf, but not VAL, encoded the direction of reaches in early but not late learning. Accordingly, bilateral lesions of Pf before, but not after learning, strongly and specifically impaired the refinement of reach direction. VAL lesions did not impact direction refinement, but instead resulted in increased speed and target overshoot. Our findings provide new evidence that the thalamus is a critical motor node in the learning and control of reaching movements, with specific subnuclei controlling distinct aspects of the reach early in learning.

Unraveling the mechanisms of deep-brain stimulation of the internal capsule in a mouse model

Deep-brain stimulation (DBS) is an effective treatment for patients suffering from otherwise therapy-resistant psychiatric disorders, including obsessive-compulsive disorder. Modulation of cortico-striatal circuits has been suggested as a mechanism of action. To gain mechanistic insight, we monitored neuronal activity in cortico-striatal regions in a mouse model for compulsive behavior, while systematically varying clinically-relevant parameters of internal-capsule DBS. DBS showed dose-dependent effects on both brain and behavior: An increasing, yet balanced, number of excited and inhibited neurons was recruited, scattered throughout cortico-striatal regions, while excessive grooming decreased. Such neuronal recruitment did not alter basic brain function such as resting-state activity, and only occurred in awake animals, indicating a dependency on network activity. In addition to these widespread effects, we observed specific involvement of the medial orbitofrontal cortex in therapeutic outcomes, which was corroborated by optogenetic stimulation. Together, our findings provide mechanistic insight into how DBS exerts its therapeutic effects on compulsive behaviors.

Interneuron-Targeted Disruption of SYNGAP1 Alters Sensory Representations in the Neocortex and Impairs Sensory Learning

SYNGAP1 haploinsufficiency in humans leads to severe neurodevelopmental disorders characterized by intellectual disability, autism, epilepsy, and sensory processing deficits. However, the circuit mechanisms underlying these disorders are not well understood. In mice, a decrease of SynGAP levels results in cognitive deficits by interfering with the development of excitatory glutamatergic connections. Recent evidence suggests that SynGAP also plays a crucial role in the development and function of GABAergic inhibitory interneurons. Nevertheless, it remains uncertain whether and to what extent the expression of SYNGAP1 in inhibitory interneurons contributes to cortical circuit function and related behaviors. The activity of cortical neurons has not been measured simultaneously with behavior. To address these gaps, we recorded from layer 2/3 neurons in the primary whisker somatosensory cortex (wS1) of mice while they learned to perform a whisker tactile detection task. Our results demonstrate that mice with interneuron-specific SYNGAP1 haploinsufficiency exhibit learning deficits characterized by heightened behavioral responses in the absence of relevant sensory input and premature responses to unrelated sensory stimuli not associated with reward acquisition. These behavioral deficits are accompanied by specific circuit abnormalities within wS1. Interneuron-specific SYNGAP1 haploinsufficiency increases detrimental neuronal correlations directly related to task performance and enhances responses to irrelevant sensory stimuli unrelated to the reward acquisition. In summary, our findings indicate that a reduction of SynGAP in inhibitory interneurons impairs sensory representation in the primary sensory cortex by disrupting neuronal correlations, which likely contributes to the observed cognitive deficits in mice with pan-neuronal SYNGAP1 haploinsufficiency.SIGNIFICANCE STATEMENT SYNGAP1 haploinsufficiency leads to severe neurodevelopmental disorders. The exact nature of neural circuit dysfunction caused by SYNGAP1 haploinsufficiency remains poorly understood. SynGAP plays a critical role in the function of GABAergic inhibitory interneurons as well as glutamatergic pyramidal neurons in the neocortex. Whether and how decreasing SYNGAP1 level in inhibitory interneurons disrupts a behaviorally relevant circuit remains unclear. We measure neural activity and behavior in mice learning a perceptual task. Mice with interneuron-targeted disruption of SYNGAP1 display increased detrimental neuronal correlations and elevated responses to irrelevant sensory inputs, which are related to impaired task performance. These results show that cortical interneuron dysfunction contributes to sensory deficits in SYNGAP1 haploinsufficiency with important implications for identifying therapeutic targets.

Social odor discrimination and its enhancement by associative learning in the hippocampal CA2 region

Although the hippocampus is crucial for social memory, how social sensory information is combined with contextual information to form episodic social memories remains unknown. Here, we investigated the mechanisms for social sensory information processing using two-photon calcium imaging from hippocampal CA2 pyramidal neurons (PNs)-which are crucial for social memory-in awake head-fixed mice exposed to social and non-social odors. We found that CA2 PNs represent social odors of individual conspecifics and that these representations are refined during associative social odor-reward learning to enhance the discrimination of rewarded compared with unrewarded odors. Moreover, the structure of the CA2 PN population activity enables CA2 to generalize along categories of rewarded versus unrewarded and social versus non-social odor stimuli. Finally, we found that CA2 is important for learning social but not non-social odor-reward associations. These properties of CA2 odor representations provide a likely substrate for the encoding of episodic social memory.

Distinct topographic organization and network activity patterns of corticocollicular neurons within layer 5 auditory cortex

The auditory cortex (AC) modulates the activity of upstream pathways in the auditory brainstem via descending (corticofugal) projections. This feedback system plays an important role in the plasticity of the auditory system by shaping response properties of neurons in many subcortical nuclei. The majority of layer (L) 5 corticofugal neurons project to the inferior colliculus (IC). This corticocollicular (CC) pathway is involved in processing of complex sounds, auditory-related learning, and defense behavior. Partly due to their location in deep cortical layers, CC neuron population activity patterns within neuronal AC ensembles remain poorly understood. We employed two-photon imaging to record the activity of hundreds of L5 neurons in anesthetized as well as awake animals. CC neurons are broader tuned than other L5 pyramidal neurons and display weaker topographic order in core AC subfields. Network activity analyses revealed stronger clusters of CC neurons compared to non-CC neurons, which respond more reliable and integrate information over larger distances. However, results obtained from secondary auditory cortex (A2) differed considerably. Here CC neurons displayed similar or higher topography, depending on the subset of neurons analyzed. Furthermore, specifically in A2, CC activity clusters formed in response to complex sounds were spatially more restricted compared to other L5 neurons. Our findings indicate distinct network mechanism of CC neurons in analyzing sound properties with pronounced subfield differences, demonstrating that the topography of sound-evoked responses within AC is neuron-type dependent.

Effective and efficient neural networks for spike inference from in vivo calcium imaging

Calcium imaging provides advantages in monitoring large populations of neuronal activities simultaneously. However, it lacks the signal quality provided by neural spike recording in traditional electrophysiology. To address this issue, we developed a supervised data-driven approach to extract spike information from calcium signals. We propose the ENS2 (effective and efficient neural networks for spike inference from calcium signals) system for spike-rate and spike-event predictions using ΔF/F0 calcium inputs based on a U-Net deep neural network. When testing on a large, ground-truth public database, it consistently outperformed state-of-the-art algorithms in both spike-rate and spike-event predictions with reduced computational load. We further demonstrated that ENS2 can be applied to analyses of orientation selectivity in primary visual cortex neurons. We conclude that it would be a versatile inference system that may benefit diverse neuroscience studies.

Quantifying uncertainty in spikes estimated from calcium imaging data

In recent years, a number of methods have been proposed to estimate the times at which a neuron spikes on the basis of calcium imaging data. However, quantifying the uncertainty associated with these estimated spikes remains an open problem. We consider a simple and well-studied model for calcium imaging data, which states that calcium decays exponentially in the absence of a spike, and instantaneously increases when a spike occurs. We wish to test the null hypothesis that the neuron did not spike-i.e., that there was no increase in calcium-at a particular timepoint at which a spike was estimated. In this setting, classical hypothesis tests lead to inflated Type I error, because the spike was estimated on the same data used for testing. To overcome this problem, we propose a selective inference approach. We describe an efficient algorithm to compute finite-sample $p$-values that control selective Type I error, and confidence intervals with correct selective coverage, for spikes estimated using a recent proposal from the literature. We apply our proposal in simulation and on calcium imaging data from the $\texttt{spikefinder}$ challenge.

Coregistration of heading to visual cues in retrosplenial cortex

Spatial cognition depends on an accurate representation of orientation within an environment. Head direction cells in distributed brain regions receive a range of sensory inputs, but visual input is particularly important for aligning their responses to environmental landmarks. To investigate how population-level heading responses are aligned to visual input, we recorded from retrosplenial cortex (RSC) of head-fixed mice in a moving environment using two-photon calcium imaging. We show that RSC neurons are tuned to the animal's relative orientation in the environment, even in the absence of head movement. Next, we found that RSC receives functionally distinct projections from visual and thalamic areas and contains several functional classes of neurons. While some functional classes mirror RSC inputs, a newly discovered class coregisters visual and thalamic signals. Finally, decoding analyses reveal unique contributions to heading from each class. Our results suggest an RSC circuit for anchoring heading representations to environmental visual landmarks.

FPGA-Based In-Vivo Calcium Image Decoding for Closed-Loop Feedback Applications

Miniaturized calcium imaging is an emerging neural recording technique that has been widely used for monitoring neural activity on a large scale at a specific brain region of rats or mice. Most existing calcium-image analysis pipelines operate offline. This results in long processing latency, making it difficult to realize closed-loop feedback stimulation for brain research. In recent work, we have proposed an FPGA-based real-time calcium image processing pipeline for closed-loop feedback applications. It can perform real-time calcium image motion correction, enhancement, fast trace extraction, and real-time decoding from extracted traces. Here, we extend this work by proposing a variety of neural network based methods for real-time decoding and evaluate the tradeoff among these decoding methods and accelerator designs. We introduce the implementation of the neural network based decoders on the FPGA, and show their speedup against the implementation on the ARM processor. Our FPGA implementation enables the real-time calcium image decoding with sub-ms processing latency for closed-loop feedback applications.

Restructuring of olfactory representations in the fly brain around odor relationships in natural sources

A core challenge of olfactory neuroscience is to understand how neural representations of odor are generated and progressively transformed across different layers of the olfactory circuit into formats that support perception and behavior. The encoding of odor by odorant receptors in the input layer of the olfactory system reflects, at least in part, the chemical relationships between odor compounds. Neural representations of odor in higher order associative olfactory areas, generated by random feedforward networks, are expected to largely preserve these input odor relationships1-3. We evaluated these ideas by examining how odors are represented at different stages of processing in the olfactory circuit of the vinegar fly D. melanogaster. We found that representations of odor in the mushroom body (MB), a third-order associative olfactory area in the fly brain, are indeed structured and invariant across flies. However, the structure of MB representational space diverged significantly from what is expected in a randomly connected network. In addition, odor relationships encoded in the MB were better correlated with a metric of the similarity of their distribution across natural sources compared to their similarity with respect to chemical features, and the converse was true for odor relationships encoded in primary olfactory receptor neurons (ORNs). Comparison of odor coding at primary, secondary, and tertiary layers of the circuit revealed that odors were significantly regrouped with respect to their representational similarity across successive stages of olfactory processing, with the largest changes occurring in the MB. The non-linear reorganization of odor relationships in the MB indicates that unappreciated structure exists in the fly olfactory circuit, and this structure may facilitate the generalization of odors with respect to their co-occurence in natural sources.

A Simple Ca2+-Imaging Approach of Network-Activity Analyses for Human Neurons

Rapid advances in light microscopy and development of all-optical electrophysiological imaging tools have greatly leveraged the speed and the depth of neurobiology studies. Calcium imaging is a common method that is useful for measuring calcium signals in cells and has been used as a functional proxy for neuronal activity. Here I describe a simple, stimulation-free approach that measures neuronal network activity and single-neuron dynamics in human neurons. This protocol provides the experimental workflow that includes step-wise illustrations of sample preparations, data processing, and analyses that can be used for quick phenotypical assessment and serves as a quick functional readout for mutagenesis or screen effort for neurodegenerative studies.

NeuroSeg-II: A deep learning approach for generalized neuron segmentation in two-photon Ca2+ imaging

The development of two-photon microscopy and Ca²⁺ indicators has enabled the recording of multiscale neuronal activities in vivo and thus advanced the understanding of brain functions. However, it is challenging to perform automatic, accurate, and generalized neuron segmentation when processing a large amount of imaging data. Here, we propose a novel deep-learning-based neural network, termed as NeuroSeg-II, to conduct automatic neuron segmentation for in vivo two-photon Ca²⁺ imaging data. This network architecture is based on Mask region-based convolutional neural network (R-CNN) but has enhancements of an attention mechanism and modified feature hierarchy modules. We added an attention mechanism module to focus the computation on neuron regions in imaging data. We also enhanced the feature hierarchy to extract feature information at diverse levels. To incorporate both spatial and temporal information in our data processing, we fused the images from average projection and correlation map extracting the temporal information of active neurons, and the integrated information was expressed as two-dimensional (2D) images. To achieve a generalized neuron segmentation, we conducted a hybrid learning strategy by training our model with imaging data from different labs, including multiscale data with different Ca²⁺ indicators. The results showed that our approach achieved promising segmentation performance across different imaging scales and Ca²⁺ indicators, even including the challenging data of large field-of-view mesoscopic images. By comparing state-of-the-art neuron segmentation methods for two-photon Ca²⁺ imaging data, we showed that our approach achieved the highest accuracy with a publicly available dataset. Thus, NeuroSeg-II enables good segmentation accuracy and a convenient training and testing process.

Stability and flexibility of odor representations in the mouse olfactory bulb

Dynamic changes in sensory representations have been basic tenants of studies in neural coding and plasticity. In olfaction, relatively little is known about the dynamic range of changes in odor representations under different brain states and over time. Here, we used time-lapse in vivo two-photon calcium imaging to describe changes in odor representation by mitral cells, the output neurons of the mouse olfactory bulb. Using anesthetics as a gross manipulation to switch between different brain states (wakefulness and under anesthesia), we found that odor representations by mitral cells undergo significant re-shaping across states but not over time within state. Odor representations were well balanced across the population in the awake state yet highly diverse under anesthesia. To evaluate differences in odor representation across states, we used linear classifiers to decode odor identity in one state based on training data from the other state. Decoding across states resulted in nearly chance-level accuracy. In contrast, repeating the same procedure for data recorded within the same state but in different time points, showed that time had a rather minor impact on odor representations. Relative to the differences across states, odor representations remained stable over months. Thus, single mitral cells can change dynamically across states but maintain robust representations across months. These findings have implications for sensory coding and plasticity in the mammalian brain.

Mapping of structural arrangement of cells and collective calcium transients: an integrated framework combining live cell imaging using confocal microscopy and UMAP-assisted HDBSCAN-based approach

Live cell calcium (Ca2+) imaging is one of the important tools to record cellular activity during in vitro and in vivo preclinical studies. Specially, high-resolution microscopy can provide valuable dynamic information at the single cell level. One of the major challenges in the implementation of such imaging schemes is to extract quantitative information in the presence of significant heterogeneity in Ca2+ responses attained due to variation in structural arrangement and drug distribution. To fill this gap, we propose time-lapse imaging using spinning disk confocal microscopy and machine learning-enabled framework for automated grouping of Ca2+ spiking patterns. Time series analysis is performed to correlate the drug induced cellular responses to self-assembly pattern present in multicellular systems. The framework is designed to reduce the large-scale dynamic responses using uniform manifold approximation and projection (UMAP). In particular, we propose the suitability of hierarchical DBSCAN (HDBSCAN) in view of reduced number of hyperparameters. We find UMAP-assisted HDBSCAN outperforms existing approaches in terms of clustering accuracy in segregation of Ca2+ spiking patterns. One of the novelties includes the application of non-linear dimension reduction in segregation of the Ca2+ transients with statistical similarity. The proposed pipeline for automation was also proved to be a reproducible and fast method with minimal user input. The algorithm was used to quantify the effect of cellular arrangement and stimulus level on collective Ca2+ responses induced by GPCR targeting drug. The analysis revealed a significant increase in subpopulation containing sustained oscillation corresponding to higher packing density. In contrast to traditional measurement of rise time and decay ratio from Ca2+ transients, the proposed pipeline was used to classify the complex patterns with longer duration and cluster-wise model fitting. The two-step process has a potential implication in deciphering biophysical mechanisms underlying the Ca2+ oscillations in context of structural arrangement between cells.

Reward expectation extinction restructures and degrades CA1 spatial maps through loss of a dopaminergic reward proximity signal

Hippocampal place cells support reward-related spatial memories by forming a cognitive map that over-represents reward locations. The strength of these memories is modulated by the extent of reward expectation during encoding. However, the circuit mechanisms underlying this modulation are unclear. Here we find that when reward expectation is extinguished in mice, they remain engaged with their environment, yet place cell over-representation of rewards vanishes, place field remapping throughout the environment increases, and place field trial-to-trial reliability decreases. Interestingly, Ventral Tegmental Area (VTA) dopaminergic axons in CA1 exhibit a ramping reward-proximity signal that depends on reward expectation and inhibiting VTA dopaminergic neurons largely replicates the effects of extinguishing reward expectation. We conclude that changing reward expectation restructures CA1 cognitive maps and determines map reliability by modulating the dopaminergic VTA-CA1 reward-proximity signal. Thus, internal states of high reward expectation enhance encoding of spatial memories by reinforcing hippocampal cognitive maps associated with reward.

Review on data analysis methods for mesoscale neural imaging in vivo

Significance: Mesoscale neural imaging in vivo has gained extreme popularity in neuroscience for its capacity of recording large-scale neurons in action. Optical imaging with single-cell resolution and millimeter-level field of view in vivo has been providing an accumulated database of neuron-behavior correspondence. Meanwhile, optical detection of neuron signals is easily contaminated by noises, background, crosstalk, and motion artifacts, while neural-level signal processing and network-level coordinate are extremely complicated, leading to laborious and challenging signal processing demands. The existing data analysis procedure remains unstandardized, which could be daunting to neophytes or neuroscientists without computational background. Aim: We hope to provide a general data analysis pipeline of mesoscale neural imaging shared between imaging modalities and systems. Approach: We divide the pipeline into two main stages. The first stage focuses on extracting high-fidelity neural responses at single-cell level from raw images, including motion registration, image denoising, neuron segmentation, and signal extraction. The second stage focuses on data mining, including neural functional mapping, clustering, and brain-wide network deduction. Results: Here, we introduce the general pipeline of processing the mesoscale neural images. We explain the principles of these procedures and compare different approaches and their application scopes with detailed discussions about the shortcomings and remaining challenges. Conclusions: There are great challenges and opportunities brought by the large-scale mesoscale data, such as the balance between fidelity and efficiency, increasing computational load, and neural network interpretability. We believe that global circuits on single-neuron level will be more extensively explored in the future.

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.

Task specificity in mouse parietal cortex

Parietal cortex is implicated in a variety of behavioral processes, but it is unknown whether and how its individual neurons participate in multiple tasks. We trained head-fixed mice to perform two visual decision tasks involving a steering wheel or a virtual T-maze and recorded from the same parietal neurons during these two tasks. Neurons that were active during the T-maze task were typically inactive during the steering-wheel task and vice versa. Recording from the same neurons in the same apparatus without task stimuli yielded the same specificity as in the task, suggesting that task specificity depends on physical context. To confirm this, we trained some mice in a third task combining the steering wheel context with the visual environment of the T-maze. This hybrid task engaged the same neurons as those engaged in the steering-wheel task. Thus, participation by neurons in mouse parietal cortex is task specific, and this specificity is determined by physical context.

Two-photon calcium imaging of neuronal activity

In vivo two-photon calcium imaging (2PCI) is a technique used for recording neuronal activity in the intact brain. It is based on the principle that, when neurons fire action potentials, intracellular calcium levels rise, which can be detected using fluorescent molecules that bind to calcium. This Primer is designed for scientists who are considering embarking on experiments with 2PCI. We provide the reader with a background on the basic concepts behind calcium imaging and on the reasons why 2PCI is an increasingly powerful and versatile technique in neuroscience. The Primer explains the different steps involved in experiments with 2PCI, provides examples of what ideal preparations should look like and explains how data are analysed. We also discuss some of the current limitations of the technique, and the types of solutions to circumvent them. Finally, we conclude by anticipating what the future of 2PCI might look like, emphasizing some of the analysis pipelines that are being developed and international efforts for data sharing.

Visual-area-specific tonic modulation of GABA release by endocannabinoids sets the activity and coordination of neocortical principal neurons

Perisomatic inhibition of pyramidal neurons (PNs) coordinates cortical network activity during sensory processing, and this role is mainly attributed to parvalbumin-expressing basket cells (BCs). However, cannabinoid receptor type 1 (CB1)-expressing interneurons are also BCs, but the connectivity and function of these elusive but prominent neocortical inhibitory neurons are unclear. We find that their connectivity pattern is visual area specific. Persistently active CB1 signaling suppresses GABA release from CB1 BCs in the medial secondary visual cortex (V2M), but not in the primary visual cortex (V1). Accordingly, in vivo, tonic CB1 signaling is responsible for higher but less coordinated PN activity in the V2M than in the V1. These differential firing dynamics in the V1 and V2M can be captured by a computational network model that incorporates visual-area-specific properties. Our results indicate a differential CB1-mediated mechanism controlling PN activity, suggesting an alternative connectivity scheme of a specific GABAergic circuit in different cortical areas.

•

CB1⁺ basket cells exhibit visual-area-specific morphology and connectivity patterns

•

Tonic CB1 signaling underlies high pyramidal neurons (PN) activity in V2M but not V1

•

Tonic CB1 signaling differentially modulates PN-correlated activity in V1 and V2M

•

Numerical simulations capture specific CB1-dependent firing dynamics of V1 and V2M

Inferring monosynaptic connections from paired dendritic spine Ca2+imaging and large-scale recording of extracellular spiking

Techniques to identify monosynaptic connections between neurons have been vital for neuroscience research, facilitating important advancements concerning network topology, synaptic plasticity, and synaptic integration, among others. Here, we introduce a novel approach to identify and monitor monosynaptic connections using high-resolution dendritic spine Ca²⁺imaging combined with simultaneous large-scale recording of extracellular electrical activity by means of high-density microelectrode arrays (HD-MEAs). We introduce an easily adoptable analysis pipeline that associates the imaged spine with its presynaptic unit and test it onin vitrorecordings. The method is further validated and optimized by simulating synaptically-evoked spine Ca²⁺transients based on measured spike trains in order to obtain simulated ground-truth connections. The proposed approach offers unique advantages asi) it can be used to identify monosynaptic connections with an accurate localization of the synapse within the dendritic tree,ii) it provides precise information of presynaptic spiking, andiii) postsynaptic spine Ca²⁺signals and, finally,iv) the non-invasive nature of the proposed method allows for long-term measurements. The analysis toolkit together with the rich data sets that were acquired are made publicly available for further exploration by the research community.

ICoRD: Iterative correlation-based ROI detection method for the extraction of neural signals in calcium imaging

OBJECTIVE

In vivo calcium imaging is a standard neuroimaging technique that allows selective observation of target neuronal activities. In calcium imaging, neuron activation signals provide key information for the investigation of neural circuits. For efficient extraction of the calcium signals of neurons, selective detection of the region of interest (ROI) pixels corresponding to the active subcellular region of the target neuron is essential. However, current ROI detection methods for calcium imaging data exhibit a relatively low signal extraction performance from neurons with a low signal-to-noise power ratio (SNR). This is problematic because a low SNR is unavoidable in many biological experiments.

APPROACH

Therefore, we propose an iterative correlation-based ROI detection (ICoRD) method that robustly extracts the calcium signal of the target neuron from a calcium imaging series with severe noise.

MAIN RESULTS

ICoRD extracts calcium signals closer to the ground-truth calcium signal than the conventional method from simulated calcium imaging data in all low SNR ranges. Additionally, this study confirmed that ICoRD robustly extracts activation signals against noise, even within in vivo environments.

SIGNIFICANCE

ICoRD showed reliable detection from neurons with a low SNR and sparse activation, which were not detected by conventional methods. ICoRD will facilitate our understanding of neural circuit activity by providing significantly improved ROI detection in noisy images.

Computational Methods for Neuron Segmentation in Two-Photon Calcium Imaging Data: A Survey

Calcium imaging has rapidly become a methodology of choice for real-time in vivo neuron analysis. Its application to large sets of data requires automated tools to annotate and segment cells, allowing scalable image segmentation under reproducible criteria. In this paper, we review and summarize the most recent methods for computational segmentation of calcium imaging. The contributions of the paper are three-fold: we provide an overview of the main algorithms taxonomized in three categories (signal processing, matrix factorization and machine learning-based approaches), we highlight the main advantages and disadvantages of each category and we provide a summary of the performance of the methods that have been tested on public benchmarks (with links to the public code when available).

Axial motion estimation and correction for simultaneous multi-plane two-photon calcium imaging

Two-photon imaging in behaving animals is typically accompanied by brain motion. For functional imaging experiments, for example with genetically encoded calcium indicators, such brain motion induces changes in fluorescence intensity. These motion-related intensity changes or motion artifacts can typically not be separated from neural activity-induced signals. While lateral motion, within the focal plane, can be corrected by computationally aligning images, axial motion, out of the focal plane, cannot easily be corrected. Here, we developed an algorithm for axial motion correction for non-ratiometric calcium indicators taking advantage of simultaneous multi-plane imaging. Using temporally multiplexed beams, recording simultaneously from at least two focal planes at different z positions, and recording a z-stack for each beam as a calibration step, the algorithm separates motion-related and neural activity-induced changes in fluorescence intensity. The algorithm is based on a maximum likelihood optimisation approach; it assumes (as a first order approximation) that no distortions of the sample occurs during axial motion and that neural activity increases uniformly along the optical axis in each region of interest. The developed motion correction approach allows axial motion estimation and correction at high frame rates for isolated structures in the imaging volume in vivo, such as sparse expression patterns in the fruit fly brain.

Large-scale two-photon calcium imaging in freely moving mice

We developed a miniaturized two-photon microscope (MINI2P) for fast, high-resolution, multiplane calcium imaging of over 1,000 neurons at a time in freely moving mice. With a microscope weight below 3 g and a highly flexible connection cable, MINI2P allowed stable imaging with no impediment of behavior in a variety of assays compared to untethered, unimplanted animals. The improved cell yield was achieved through a optical system design featuring an enlarged field of view (FOV) and a microtunable lens with increased z-scanning range and speed that allows fast and stable imaging of multiple interleaved planes, as well as 3D functional imaging. Successive imaging across multiple, adjacent FOVs enabled recordings from more than 10,000 neurons in the same animal. Large-scale proof-of-principle data were obtained from cell populations in visual cortex, medial entorhinal cortex, and hippocampus, revealing spatial tuning of cells in all areas.

•

We made a light-weight 2-photon miniscope for calcium imaging in freely moving mice

•

Stable high-quality imaging was observed during a wide spectrum of behaviors

•

Activity can be monitored in volumes of over 1,000 visual or entorhinal-cortex cells

•

A custom-designed z-scanning module allows fast imaging across multiple planes

Calcium Imaging and Electrophysiology of hippocampal Activity under Anesthesia and natural Sleep in Mice

The acute effects of anesthesia and their underlying mechanisms are still not fully understood. Thus, comprehensive analysis and efficient generalization require their description in various brain regions. Here we describe a large-scale, annotated collection of 2-photon calcium imaging data and multi-electrode, extracellular electrophysiological recordings in CA1 of the murine hippocampus under three distinct anesthetics (Isoflurane, Ketamine/Xylazine and Medetomidine/Midazolam/Fentanyl), during natural sleep, and wakefulness. We cover several aspects of data quality standardization and provide a set of tools for autonomous validation, along with analysis workflows for reuse and data exploration. The datasets described here capture various aspects of neural activity in hundreds of pyramidal cells at single cell resolution. In addition to relevance for basic biological research, the dataset may find utility in computational neuroscience as a benchmark for models of anesthesia and sleep.

To deconvolve, or not to deconvolve: Inferences of neuronal activities using calcium imaging data

BACKGROUND

With the increasing popularity of calcium imaging in neuroscience research, choosing the right methods to analyze calcium imaging data is critical to address various scientific questions. Unlike spike trains measured using electrodes, fluorescence intensity traces provide an indirect and noisy measurement of the underlying neuronal activities. The observed calcium traces are either analyzed directly or deconvolved to spike trains to infer neuronal activities. When both approaches are applicable, it is unclear whether deconvolving calcium traces is a necessary step.

METHODS

In this article, we compare the performance of using calcium traces or their deconvolved spike trains for three common analyses: clustering, principal component analysis (PCA), and population decoding.

RESULTS

We found that (1) the two approaches lead to diverging results; (2) estimated spike trains, when smoothed or binned appropriately, usually lead to satisfactory performances, such as more accurate estimation of cluster membership; (3) although estimate spike train produce results more similar to true spike data than trace data, we found that the PCA results from trace data might better reflect the underlying neuronal ensembles (clusters); and (4) for both approaches, decobability can be improved by using denoising or smoothing methods.

COMPARISON WITH EXISTING METHODS

Our simulations and applications to real data suggest that estimated spike data outperform trace data in cluster analysis and give comparable results for population decoding. In addition, the decobability of estimated spike data can be slightly better than that of calcium trace data with appropriate filtering / smoothing methods.

CONCLUSION

We conclude that spike detection might be a useful pre-processing step for certain problems such as clustering; however, the continuous nature of calcium imaging data provides a natural smoothness that might be helpful for problems such as dimensional reduction.

FARCI: Fast and Robust Connectome Inference

The inference of neuronal connectome from large-scale neuronal activity recordings, such as two-photon Calcium imaging, represents an active area of research in computational neuroscience. In this work, we developed FARCI (Fast and Robust Connectome Inference), a MATLAB package for neuronal connectome inference from high-dimensional two-photon Calcium fluorescence data. We employed partial correlations as a measure of the functional association strength between pairs of neurons to reconstruct a neuronal connectome. We demonstrated using in silico datasets from the Neural Connectomics Challenge (NCC) and those generated using the state-of-the-art simulator of Neural Anatomy and Optimal Microscopy (NAOMi) that FARCI provides an accurate connectome and its performance is robust to network sizes, missing neurons, and noise levels. Moreover, FARCI is computationally efficient and highly scalable to large networks. In comparison with the best performing connectome inference algorithm in the NCC, Generalized Transfer Entropy (GTE), and Fluorescence Single Neuron and Network Analysis Package (FluoroSNNAP), FARCI produces more accurate networks over different network sizes, while providing significantly better computational speed and scaling.

Fluorescent Biosensors Based on Silicon Nanowires

Nanostructures are arising as novel biosensing platforms promising to surpass current performance in terms of sensitivity, selectivity, and affordability of standard approaches. However, for several nanosensors, the material and synthesis used make the industrial transfer of such technologies complex. Silicon nanowires (NWs) are compatible with Si-based flat architecture fabrication and arise as a hopeful solution to couple their interesting physical properties and surface-to-volume ratio to an easy commercial transfer. Among all the transduction methods, fluorescent probes and sensors emerge as some of the most used approaches thanks to their easy data interpretation, measure affordability, and real-time in situ analysis. In fluorescent sensors, Si NWs are employed as substrate and coupled with several fluorophores, NWs can be used as quenchers in stem-loop configuration, and have recently been used for direct fluorescent sensing. In this review, an overview on fluorescent sensors based on Si NWs is presented, analyzing the literature of the field and highlighting the advantages and drawbacks for each strategy.

Reliable Sensory Processing in Mouse Visual Cortex through Cooperative Interactions between Somatostatin and Parvalbumin Interneurons

Intrinsic neuronal variability significantly limits information encoding in the primary visual cortex (V1). However, under certain conditions, neurons can respond reliably with highly precise responses to the same visual stimuli from trial to trial. This suggests that there exists intrinsic neural circuit mechanisms that dynamically modulate the intertrial variability of visual cortical neurons. Here, we sought to elucidate the role of different inhibitory interneurons (INs) in reliable coding in mouse V1. To study the interactions between somatostatin-expressing interneurons (SST-INs) and parvalbumin-expressing interneurons (PV-INs), we used a dual-color calcium imaging technique that allowed us to simultaneously monitor these two neural ensembles while awake mice, of both sexes, passively viewed natural movies. SST neurons were more active during epochs of reliable pyramidal neuron firing, whereas PV neurons were more active during epochs of unreliable firing. SST neuron activity lagged that of PV neurons, consistent with a feedback inhibitory SST→PV circuit. To dissect the role of this circuit in pyramidal neuron activity, we used temporally limited optogenetic activation and inactivation of SST and PV interneurons during periods of reliable and unreliable pyramidal cell firing. Transient firing of SST neurons increased pyramidal neuron reliability by actively suppressing PV neurons, a proposal that was supported by a rate-based model of V1 neurons. These results identify a cooperative functional role for the SST→PV circuit in modulating the reliability of pyramidal neuron activity.SIGNIFICANCE STATEMENT Cortical neurons often respond to identical sensory stimuli with large variability. However, under certain conditions, the same neurons can also respond highly reliably. The circuit mechanisms that contribute to this modulation remain unknown. Here, we used novel dual-wavelength calcium imaging and temporally selective optical perturbation to identify an inhibitory neural circuit in visual cortex that can modulate the reliability of pyramidal neurons to naturalistic visual stimuli. Our results, supported by computational models, suggest that somatostatin interneurons increase pyramidal neuron reliability by suppressing parvalbumin interneurons via the inhibitory SST→PV circuit. These findings reveal a novel role of the SST→PV circuit in modulating the fidelity of neural coding critical for visual perception.

The habenula clock influences response to a stressor

The response of an animal to a sensory stimulus depends on the nature of the stimulus and on expectations, which are mediated by spontaneous activity. Here, we ask how circadian variation in the expectation of danger, and thus the response to a potential threat, is controlled. We focus on the habenula, a mediator of threat response that functions by regulating neuromodulator release, and use zebrafish as the experimental system. Single cell transcriptomics indicates that multiple clock genes are expressed throughout the habenula, while quantitative in situ hybridization confirms that the clock oscillates. Two-photon calcium imaging indicates a circadian change in spontaneous activity of habenula neurons. To assess the role of this clock, a truncated clocka gene was specifically expressed in the habenula. This partially inhibited the clock, as shown by changes in per3 expression as well as altered day-night variation in dopamine, serotonin and acetylcholine levels. Behaviourally, anxiety-like responses evoked by an alarm pheromone were reduced. Circadian effects of the pheromone were disrupted, such that responses in the day resembled those at night. Behaviours that are regulated by the pineal clock and not triggered by stressors were unaffected. We suggest that the habenula clock regulates the expectation of danger, thus providing one mechanism for circadian change in the response to a stressor.

State-Dependent Changes in Perception and Coding in the Mouse Somatosensory Cortex

An animal's behavioral state is reflected in the dynamics of cortical population activity and its capacity to process sensory information. To better understand the relationship between behavioral states and information processing, mice are trained to detect varying amplitudes of whisker-deflection under two-photon calcium imaging. Layer 2/3 neurons in the vibrissal primary somatosensory cortex are imaged across different behavioral states, defined based on detection performance (low to high-state) and pupil diameter. The neurometric curve in each behavioral state mirrors the corresponding psychometric performance, with calcium signals predictive of the animal's choice. High behavioral states are associated with lower network synchrony, extending over shorter cortical distances. The decrease in correlation across neurons in high state results in enhanced information transmission capacity at the population level. The observed state-dependent changes suggest that the coding regime within the first stage of cortical processing may underlie adaptive routing of relevant information through the sensorimotor system.

Neural anatomy and optical microscopy (NAOMi) simulation for evaluating calcium imaging methods

BACKGROUND

The past decade has seen a multitude of new in vivo functional imaging methodologies. However, the lack of ground-truth comparisons or evaluation metrics makes the large-scale, systematic validation vital to the continued development and use of optical microscopy impossible.

NEW-METHOD

We provide a new framework for evaluating two-photon microscopy methods via in silico Neural Anatomy and Optical Microscopy (NAOMi) simulation. Our computationally efficient model generates large anatomical volumes of mouse cortex, simulates neural activity, and incorporates optical propagation and scanning to create realistic calcium imaging datasets.

RESULTS

We verify NAOMi simulations against in vivo two-photon recordings from mouse cortex. We leverage this in silico ground truth to directly compare different segmentation algorithms and optical designs. We find modern segmentation algorithms extract strong neural time-courses comparable to estimation using oracle spatial information, but with an increase in the false positive rate. Comparison between optical setups demonstrate improved resilience to motion artifacts in sparsely labeled samples using Bessel beams, increased signal-to-noise ratio and cell-count using low numerical aperture Gaussian beams and nuclear GCaMP, and more uniform spatial sampling with temporal focusing versus multi-plane imaging.

COMPARISON WITH EXISTING METHODS

NAOMi is a first-of-its kind framework for assessing optical imaging modalities. Existing methods are either anatomical simulations or do not address functional imaging. Thus there is no competing method for simulating realistic functional optical microscopy data.

CONCLUSIONS

By leveraging the rich accumulated knowledge of neural anatomy and optical physics, we provide a powerful new tool to assess and develop important methods in neural imaging.

Direct extraction of signal and noise correlations from two-photon calcium imaging of ensemble neuronal activity

Neuronal activity correlations are key to understanding how populations of neurons collectively encode information. While two-photon calcium imaging has created a unique opportunity to record the activity of large populations of neurons, existing methods for inferring correlations from these data face several challenges. First, the observations of spiking activity produced by two-photon imaging are temporally blurred and noisy. Secondly, even if the spiking data were perfectly recovered via deconvolution, inferring network-level features from binary spiking data is a challenging task due to the non-linear relation of neuronal spiking to endogenous and exogenous inputs. In this work, we propose a methodology to explicitly model and directly estimate signal and noise correlations from two-photon fluorescence observations, without requiring intermediate spike deconvolution. We provide theoretical guarantees on the performance of the proposed estimator and demonstrate its utility through applications to simulated and experimentally recorded data from the mouse auditory cortex.

Cellular and Widefield Imaging of Sound Frequency Organization in Primary and Higher Order Fields of the Mouse Auditory Cortex

The mouse auditory cortex (ACtx) contains two core fields-primary auditory cortex (A1) and anterior auditory field (AAF)-arranged in a mirror reversal tonotopic gradient. The best frequency (BF) organization and naming scheme for additional higher order fields remain a matter of debate, as does the correspondence between smoothly varying global tonotopy and heterogeneity in local cellular tuning. Here, we performed chronic widefield and two-photon calcium imaging from the ACtx of awake Thy1-GCaMP6s reporter mice. Data-driven parcellation of widefield maps identified five fields, including a previously unidentified area at the ventral posterior extreme of the ACtx (VPAF) and a tonotopically organized suprarhinal auditory field (SRAF) that extended laterally as far as ectorhinal cortex. Widefield maps were stable over time, where single pixel BFs fluctuated by less than 0.5 octaves throughout a 1-month imaging period. After accounting for neuropil signal and frequency tuning strength, BF organization in neighboring layer 2/3 neurons was intermediate to the heterogeneous salt and pepper organization and the highly precise local organization that have each been described in prior studies. Multiscale imaging data suggest there is no ultrasonic field or secondary auditory cortex in the mouse. Instead, VPAF and a dorsal posterior (DP) field emerged as the strongest candidates for higher order auditory areas.

Prefrontal deep projection neurons enable cognitive flexibility via persistent feedback monitoring

Cognitive flexibility, the ability to alter strategy according to changing stimulus-response-reward relationships, is critical for updating learned behavior. Attentional set-shifting, a test of cognitive flexibility, depends on the activity of prefrontal cortex (PFC). It remains unclear, however, what role PFC neurons play to support set-shifting. Using optogenetics and two-photon calcium imaging, we demonstrate that medial PFC activity does not bias sensorimotor responses during set-shifting, but rather enables set-shifting by encoding trial feedback information, a role it has been known to play in other contexts. Unexpectedly, the functional properties of PFC cells did not vary with their efferent projection targets. Instead, representations of trial feedback formed a topological gradient, with cells more strongly selective for feedback information located further from the pial surface, where afferent input from the anterior cingulate cortex was denser. These findings identify a critical role for deep PFC projection neurons in enabling set-shifting through behavioral feedback monitoring.

Calcium Imaging and the Curse of Negativity

The imaging of neuronal activity using calcium indicators has become a staple of modern neuroscience. However, without ground truths, there is a real risk of missing a significant portion of the real responses. Here, we show that a common assumption, the non-negativity of the neuronal responses as detected by calcium indicators, biases all levels of the frequently used analytical methods for these data. From the extraction of meaningful fluorescence changes to spike inference and the analysis of inferred spikes, each step risks missing real responses because of the assumption of non-negativity. We first show that negative deviations from baseline can exist in calcium imaging of neuronal activity. Then, we use simulated data to test three popular algorithms for image analysis, CaImAn, suite2p, and CellSort, finding that suite2p may be the best suited to large datasets. We also tested the spike inference algorithms included in CaImAn, suite2p, and Cellsort, as well as the dedicated inference algorithms MLspike and CASCADE, and found each to have limitations in dealing with inhibited neurons. Among these spike inference algorithms, FOOPSI, from CaImAn, performed the best on inhibited neurons, but even this algorithm inferred spurious spikes upon the return of the fluorescence signal to baseline. As such, new approaches will be needed before spikes can be sensitively and accurately inferred from calcium data in inhibited neurons. We further suggest avoiding data analysis approaches that, by assuming non-negativity, ignore inhibited responses. Instead, we suggest a first exploratory step, using k-means or PCA for example, to detect whether meaningful negative deviations are present. Taking these steps will ensure that inhibition, as well as excitation, is detected in calcium imaging datasets.

From Neurons to Cognition: Technologies for Precise Recording of Neural Activity Underlying Behavior

Understanding how brain activity encodes information and controls behavior is a long-standing question in neuroscience. This complex problem requires converging efforts from neuroscience and engineering, including technological solutions to perform high-precision and large-scale recordings of neuronal activity in vivo as well as unbiased methods to reliably measure and quantify behavior. Thanks to advances in genetics, molecular biology, engineering, and neuroscience, in recent decades, a variety of optical imaging and electrophysiological approaches for recording neuronal activity in awake animals have been developed and widely applied in the field. Moreover, sophisticated computer vision and machine learning algorithms have been developed to analyze animal behavior. In this review, we provide an overview of the current state of technology for neuronal recordings with a focus on optical and electrophysiological methods in rodents. In addition, we discuss areas that future technological development will need to cover in order to further our understanding of the neural activity underlying behavior.

A Systematic Evaluation of Interneuron Morphology Representations for Cell Type Discrimination

Quantitative analysis of neuronal morphologies usually begins with choosing a particular feature representation in order to make individual morphologies amenable to standard statistics tools and machine learning algorithms. Many different feature representations have been suggested in the literature, ranging from density maps to intersection profiles, but they have never been compared side by side. Here we performed a systematic comparison of various representations, measuring how well they were able to capture the difference between known morphological cell types. For our benchmarking effort, we used several curated data sets consisting of mouse retinal bipolar cells and cortical inhibitory neurons. We found that the best performing feature representations were two-dimensional density maps, two-dimensional persistence images and morphometric statistics, which continued to perform well even when neurons were only partially traced. Combining these feature representations together led to further performance increases suggesting that they captured non-redundant information. The same representations performed well in an unsupervised setting, implying that they can be suitable for dimensionality reduction or clustering.

Improved hyperacuity estimation of spike timing from calcium imaging

Two-photon imaging is a major recording technique used in neuroscience. However, it suffers from several limitations, including a low sampling rate, the nonlinearity of calcium responses, the slow dynamics of calcium dyes and a low SNR, all of which severely limit the potential of two-photon imaging to elucidate neuronal dynamics with high temporal resolution. We developed a hyperacuity algorithm (HA_time) based on an approach that combines a generative model and machine learning to improve spike detection and the precision of spike time inference. Bayesian inference was performed to estimate the calcium spike model, assuming constant spike shape and size. A support vector machine using this information and a jittering method maximizing the likelihood of estimated spike times enhanced spike time estimation precision approximately fourfold (range, 2–7; mean, 3.5–4.0; 2SEM, 0.1–0.25) compared to the sampling interval. Benchmark scores of HA_time for biological data from three different brain regions were among the best of the benchmark algorithms. Simulation of broader data conditions indicated that our algorithm performed better than others with high firing rate conditions. Furthermore, HA_time exhibited comparable performance for conditions with and without ground truths. Thus HA_time is a useful tool for spike reconstruction from two-photon imaging.

On optimal coupling of the 'electronic photoreceptors' into the degenerate retina

Objective

To restore sight in atrophic age-related macular degeneration, the lost photoreceptors can be replaced with electronic implants, which replicate their two major functions: (1) converting light into an electric signal, and (2) transferring visual information to the secondary neurons in the retinal neural network—the bipolar cells (BC). We study the selectivity of BC activation by subretinal implants and dynamics of their response to pulsatile waveforms in order to optimize the electrical stimulation scheme such that retinal signal processing with 'electronic photoreceptors' remains as close to natural as possible.

Approach

A multicompartmental model of a BC was implemented to simulate responses of the voltage-gated calcium channels and subsequent synaptic vesicle release under continuous and pulsatile stimuli. We compared the predicted response under various frequencies, pulse durations, and alternating gratings to the corresponding experimental measurements. In addition, electric field was computed for various electrode configurations in a 3-d finite element model to assess the stimulation selectivity via spatial confinement of the field.

Main results

The modeled BC-mediated retinal responses were, in general, in good agreement with previously published experimental results. Kinetics of the calcium pumps and of the neurotransmitter release in ribbon synapses, which underpin the BC's temporal filtering and rectifying functions, allow mimicking the natural BC response with high frequency pulsatile stimulation, thereby preserving features of the retinal signal processing, such as flicker fusion, adaptation to static stimuli and non-linear summation of subunits in receptive field. Selectivity of the BC stimulation while avoiding direct activation of the downstream neurons (amacrine and ganglion cells—RGCs) is improved with local return electrodes.

Significance

If the retinal neural network is preserved to a large extent in age-related macular degeneration, selective stimulation of BCs with proper spatial and temporal modulation of the extracellular electric field may retain many features of the natural retinal signal processing and hence allow highly functional restoration of sight.

Real-Time Neuron Detection and Neural Signal Extraction Platform for Miniature Calcium Imaging

Real-time neuron detection and neural activity extraction are critical components of real-time neural decoding. In this paper, we propose a novel real-time neuron detection and activity extraction system using a dataflow framework to provide real-time performance and adaptability to new algorithms and hardware platforms. The proposed system was evaluated on simulated calcium imaging data, calcium imaging data with manual annotation, and calcium imaging data of the anterior lateral motor cortex. We found that the proposed system accurately detected neurons and extracted neural activities in real time without any requirement for expensive, cumbersome, or special-purpose computing hardware. We expect that the system will enable cost-effective, real-time calcium imaging-based neural decoding, leading to precise neuromodulation.