Community-based benchmarking improves spike rate inference from two-photon calcium imaging data

A machine learning toolbox for the analysis of sharp-wave ripples reveals common waveform features across species

The study of sharp-wave ripples has advanced our understanding of memory function, and their alteration in neurological conditions such as epilepsy is considered a biomarker of dysfunction. Sharp-wave ripples exhibit diverse waveforms and properties that cannot be fully characterized by spectral methods alone. Here, we describe a toolbox of machine-learning models for automatic detection and analysis of these events. The machine-learning architectures, which resulted from a crowdsourced hackathon, are able to capture a wealth of ripple features recorded in the dorsal hippocampus of mice across awake and sleep conditions. When applied to data from the macaque hippocampus, these models are able to generalize detection and reveal shared properties across species. We hereby provide a user-friendly open-source toolbox for model use and extension, which can help to accelerate and standardize analysis of sharp-wave ripples, lowering the threshold for its adoption in biomedical applications.

On the contrast response function of adapted neural populations

The magnitude of neural responses in sensory cortex depends on the intensity of a stimulus and its probability of being observed within the environment. How these two variables combine to influence the overall response of cortical populations remains unknown. Here we show that, in primary visual cortex, the vector magnitude of the population response is described by a separable power law that factors the intensity of a stimulus and its probability. Moreover, the discriminability between two contrast levels in a cortical population is proportional to the logarithm of the contrast ratio.NEW & NOTEWORTHY The magnitude of neural responses in sensory cortex depends on the intensity of a stimulus and its probability of being observed within the environment. The authors show that, in primary visual cortex, the vector magnitude of the population response is described by a separable power law that factors the intensity of a stimulus and its probability.

A power law describes the magnitude of adaptation in neural populations of primary visual cortex

How do neural populations adapt to the time-varying statistics of sensory input? We used two-photon imaging to measure the activity of neurons in mouse primary visual cortex adapted to different sensory environments, each defined by a distinct probability distribution over a stimulus set. We find that two properties of adaptation capture how the population response to a given stimulus, viewed as a vector, changes across environments. First, the ratio between the response magnitudes is a power law of the ratio between the stimulus probabilities. Second, the response direction to a stimulus is largely invariant. These rules could be used to predict how cortical populations adapt to novel, sensory environments. Finally, we show how the power law enables the cortex to preferentially signal unexpected stimuli and to adjust the metabolic cost of its sensory representation to the entropy of the environment.

On the contrast response function of adapted neural populations

The magnitude of neural responses in sensory cortex depends on the intensity of a stimulus and its probability of being observed within the environment. How these two variables combine to influence the overall response of cortical populations remains unknown. Here we show that, in primary visual cortex, the vector magnitude of the population response is described by a separable power-law that factors the intensity of a stimulus and its probability.

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.

All-optical recreation of naturalistic neural activity with a multifunctional transgenic reporter mouse

Determining which features of the neural code drive behavior requires the ability to simultaneously read out and write in neural activity patterns with high precision across many neurons. All-optical systems that combine two-photon calcium imaging and targeted photostimulation enable the activation of specific, functionally defined groups of neurons. However, these techniques are unable to test how patterns of activity across a population contribute to computation because of an inability to both read and write cell-specific firing rates. To overcome this challenge, we make two advances: first, we introduce a genetic line of mice for Cre-dependent co-expression of a calcium indicator and a potent soma-targeted microbial opsin. Second, using this line, we develop a method for read-out and write-in of precise population vectors of neural activity by calibrating the photostimulation to each cell. These advances offer a powerful and convenient platform for investigating the neural codes of computation and behavior.

A machine learning toolbox for the analysis of sharp-wave ripples reveal common features across species

The study of sharp-wave ripples (SWRs) has advanced our understanding of memory function, and their alteration in neurological conditions such as epilepsy and Alzheimer's disease is considered a biomarker of dysfunction. SWRs exhibit diverse waveforms and properties that cannot be fully characterized by spectral methods alone. Here, we describe a toolbox of machine learning (ML) models for automatic detection and analysis of SWRs. The ML architectures, which resulted from a crowdsourced hackathon, are able to capture a wealth of SWR features recorded in the dorsal hippocampus of mice. When applied to data from the macaque hippocampus, these models were able to generalize detection and revealed shared SWR properties across species. We hereby provide a user-friendly open-source toolbox for model use and extension, which can help to accelerate and standardize SWR research, lowering the threshold for its adoption in biomedical applications.

Identifying properties of pattern completion neurons in a computational model of the visual cortex

Neural ensembles are found throughout the brain and are believed to underlie diverse cognitive functions including memory and perception. Methods to activate ensembles precisely, reliably, and quickly are needed to further study the ensembles' role in cognitive processes. Previous work has found that ensembles in layer 2/3 of the visual cortex (V1) exhibited pattern completion properties: ensembles containing tens of neurons were activated by stimulation of just two neurons. However, methods that identify pattern completion neurons are underdeveloped. In this study, we optimized the selection of pattern completion neurons in simulated ensembles. We developed a computational model that replicated the connectivity patterns and electrophysiological properties of layer 2/3 of mouse V1. We identified ensembles of excitatory model neurons using K-means clustering. We then stimulated pairs of neurons in identified ensembles while tracking the activity of the entire ensemble. Our analysis of ensemble activity quantified a neuron pair's power to activate an ensemble using a novel metric called pattern completion capability (PCC) based on the mean pre-stimulation voltage across the ensemble. We found that PCC was directly correlated with multiple graph theory parameters, such as degree and closeness centrality. To improve selection of pattern completion neurons in vivo, we computed a novel latency metric that was correlated with PCC and could potentially be estimated from modern physiological recordings. Lastly, we found that stimulation of five neurons could reliably activate ensembles. These findings can help researchers identify pattern completion neurons to stimulate in vivo during behavioral studies to control ensemble activation.

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.

A power law of cortical adaptation

How do neural populations adapt to the time-varying statistics of sensory input? To investigate, we measured the activity of neurons in primary visual cortex adapted to different environments, each associated with a distinct probability distribution over a stimulus set. Within each environment, a stimulus sequence was generated by independently sampling form its distribution. We find that two properties of adaptation capture how the population responses to a given stimulus, viewed as vectors, are linked across environments. First, the ratio between the response magnitudes is a power law of the ratio between the stimulus probabilities. Second, the response directions are largely invariant. These rules can be used to predict how cortical populations adapt to novel, sensory environments. Finally, we show how the power law enables the cortex to preferentially signal unexpected stimuli and to adjust the metabolic cost of its sensory representation to the entropy of the environment.

Learnable latent embeddings for joint behavioural and neural analysis

Mapping behavioural actions to neural activity is a fundamental goal of neuroscience. As our ability to record large neural and behavioural data increases, there is growing interest in modelling neural dynamics during adaptive behaviours to probe neural representations1-3. In particular, although neural latent embeddings can reveal underlying correlates of behaviour, we lack nonlinear techniques that can explicitly and flexibly leverage joint behaviour and neural data to uncover neural dynamics3-5. Here, we fill this gap with a new encoding method, CEBRA, that jointly uses behavioural and neural data in a (supervised) hypothesis- or (self-supervised) discovery-driven manner to produce both consistent and high-performance latent spaces. We show that consistency can be used as a metric for uncovering meaningful differences, and the inferred latents can be used for decoding. We validate its accuracy and demonstrate our tool's utility for both calcium and electrophysiology datasets, across sensory and motor tasks and in simple or complex behaviours across species. It allows leverage of single- and multi-session datasets for hypothesis testing or can be used label free. Lastly, we show that CEBRA can be used for the mapping of space, uncovering complex kinematic features, for the production of consistent latent spaces across two-photon and Neuropixels data, and can provide rapid, high-accuracy decoding of natural videos from visual cortex.

Associations between in vitro, in vivo and in silico cell classes in mouse primary visual cortex

The brain consists of many cell classes yet in vivo electrophysiology recordings are typically unable to identify and monitor their activity in the behaving animal. Here, we employed a systematic approach to link cellular, multi-modal in vitro properties from experiments with in vivo recorded units via computational modeling and optotagging experiments. We found two one-channel and six multi-channel clusters in mouse visual cortex with distinct in vivo properties in terms of activity, cortical depth, and behavior. We used biophysical models to map the two one- and the six multi-channel clusters to specific in vitro classes with unique morphology, excitability and conductance properties that explain their distinct extracellular signatures and functional characteristics. These concepts were tested in ground-truth optotagging experiments with two inhibitory classes unveiling distinct in vivo properties. This multi-modal approach presents a powerful way to separate in vivo clusters and infer their cellular properties from first principles.

Monocular deprivation during the critical period alters neuronal tuning and the composition of visual circuitry

Abnormal visual experience during a developmental critical period degrades cortical responsiveness. Yet how experience-dependent plasticity alters the response properties of individual neurons and composition of visual circuitry is unclear. Here, we measured with calcium imaging in alert mice how monocular deprivation (MD) during the developmental critical period affects tuning for binocularity, orientation, and spatial frequency for neurons in primary visual cortex. MD of the contralateral eye did not uniformly shift ocular dominance (OD) of neurons towards the fellow ipsilateral eye but reduced the number of monocular contralateral neurons and increased the number of monocular ipsilateral neurons. MD also impaired matching of preferred orientation for binocular neurons and reduced the percentage of neurons responsive at most spatial frequencies for the deprived contralateral eye. Tracking the tuning properties for several hundred neurons before and after MD revealed that the shift in OD is complex and dynamic, with many previously monocular neurons becoming binocular and binocular neurons becoming monocular. Binocular neurons that became monocular were more likely to lose responsiveness to the deprived contralateral eye if they were better matched for orientation prior to deprivation. In addition, the composition of visual circuitry changed as population of neurons more responsive to the deprived eye were exchanged for neurons with tuning properties more similar to the network of responsive neurons altered by MD. Thus, plasticity during the critical period adapts to recent experience by both altering the tuning of responsive neurons and recruiting neurons with matching tuning properties.

Associations between in vitro , in vivo and in silico cell classes in mouse primary visual cortex

The brain consists of many cell classes yet in vivo electrophysiology recordings are typically unable to identify and monitor their activity in the behaving animal. Here, we employed a systematic approach to link cellular, multi-modal in vitro properties from experiments with in vivo recorded units via computational modeling and optotagging experiments. We found two one-channel and six multi-channel clusters in mouse visual cortex with distinct in vivo properties in terms of activity, cortical depth, and behavior. We used biophysical models to map the two one- and the six multi-channel clusters to specific in vitro classes with unique morphology, excitability and conductance properties that explain their distinct extracellular signatures and functional characteristics. These concepts were tested in ground-truth optotagging experiments with two inhibitory classes unveiling distinct in vivo properties. This multi-modal approach presents a powerful way to separate in vivo clusters and infer their cellular properties from first principles.

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.

Mice and primates use distinct strategies for visual segmentation

The rodent visual system has attracted great interest in recent years due to its experimental tractability, but the fundamental mechanisms used by the mouse to represent the visual world remain unclear. In the primate, researchers have argued from both behavioral and neural evidence that a key step in visual representation is 'figure-ground segmentation', the delineation of figures as distinct from backgrounds. To determine if mice also show behavioral and neural signatures of figure-ground segmentation, we trained mice on a figure-ground segmentation task where figures were defined by gratings and naturalistic textures moving counterphase to the background. Unlike primates, mice were severely limited in their ability to segment figure from ground using the opponent motion cue, with segmentation behavior strongly dependent on the specific carrier pattern. Remarkably, when mice were forced to localize naturalistic patterns defined by opponent motion, they adopted a strategy of brute force memorization of texture patterns. In contrast, primates, including humans, macaques, and mouse lemurs, could readily segment figures independent of carrier pattern using the opponent motion cue. Consistent with mouse behavior, neural responses to the same stimuli recorded in mouse visual areas V1, RL, and LM also did not support texture-invariant segmentation of figures using opponent motion. Modeling revealed that the texture dependence of both the mouse's behavior and neural responses could be explained by a feedforward neural network lacking explicit segmentation capabilities. These findings reveal a fundamental limitation in the ability of mice to segment visual objects compared to primates.

A deep learning framework for inference of single-trial neural population dynamics from calcium imaging with subframe temporal resolution

In many areas of the brain, neural populations act as a coordinated network whose state is tied to behavior on a millisecond timescale. Two-photon (2p) calcium imaging is a powerful tool to probe such network-scale phenomena. However, estimating the network state and dynamics from 2p measurements has proven challenging because of noise, inherent nonlinearities and limitations on temporal resolution. Here we describe Recurrent Autoencoder for Discovering Imaged Calcium Latents (RADICaL), a deep learning method to overcome these limitations at the population level. RADICaL extends methods that exploit dynamics in spiking activity for application to deconvolved calcium signals, whose statistics and temporal dynamics are quite distinct from electrophysiologically recorded spikes. It incorporates a new network training strategy that capitalizes on the timing of 2p sampling to recover network dynamics with high temporal precision. In synthetic tests, RADICaL infers the network state more accurately than previous methods, particularly for high-frequency components. In 2p recordings from sensorimotor areas in mice performing a forelimb reach task, RADICaL infers network state with close correspondence to single-trial variations in behavior and maintains high-quality inference even when neuronal populations are substantially reduced.

Multiregion neuronal activity: the forest and the trees

The past decade has witnessed remarkable advances in the simultaneous measurement of neuronal activity across many brain regions, enabling fundamentally new explorations of the brain-spanning cellular dynamics that underlie sensation, cognition and action. These recently developed multiregion recording techniques have provided many experimental opportunities, but thoughtful consideration of methodological trade-offs is necessary, especially regarding field of view, temporal acquisition rate and ability to guarantee cellular resolution. When applied in concert with modern optogenetic and computational tools, multiregion recording has already made possible fundamental biological discoveries - in part via the unprecedented ability to perform unbiased neural activity screens for principles of brain function, spanning dozens of brain areas and from local to global scales.

Conceptual framework for neuronal ensemble identification and manipulation related to behavior using calcium imaging

Significance: The identification and manipulation of spatially identified neuronal ensembles with optical methods have been recently used to prove the causal link between neuronal ensemble activity and learned behaviors. However, the standardization of a conceptual framework to identify and manipulate neuronal ensembles from calcium imaging recordings is still lacking. Aim: We propose a conceptual framework for the identification and manipulation of neuronal ensembles using simultaneous calcium imaging and two-photon optogenetics in behaving mice. Approach: We review the computational approaches that have been used to identify and manipulate neuronal ensembles with single cell resolution during behavior in different brain regions using all-optical methods. Results: We proposed three steps as a conceptual framework that could be applied to calcium imaging recordings to identify and manipulate neuronal ensembles in behaving mice: (1) transformation of calcium transients into binary arrays; (2) identification of neuronal ensembles as similar population vectors; and (3) targeting of neuronal ensemble members that significantly impact behavioral performance. Conclusions: The use of simultaneous two-photon calcium imaging and two-photon optogenetics allowed for the experimental demonstration of the causal relation of population activity and learned behaviors. The standardization of analytical tools to identify and manipulate neuronal ensembles could accelerate interventional experiments aiming to reprogram the brain in normal and pathological conditions.

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.

Sustained deep-tissue voltage recording using a fast indicator evolved for two-photon microscopy

Genetically encoded voltage indicators are emerging tools for monitoring voltage dynamics with cell-type specificity. However, current indicators enable a narrow range of applications due to poor performance under two-photon microscopy, a method of choice for deep-tissue recording. To improve indicators, we developed a multiparameter high-throughput platform to optimize voltage indicators for two-photon microscopy. Using this system, we identified JEDI-2P, an indicator that is faster, brighter, and more sensitive and photostable than its predecessors. We demonstrate that JEDI-2P can report light-evoked responses in axonal termini of Drosophila interneurons and the dendrites and somata of amacrine cells of isolated mouse retina. JEDI-2P can also optically record the voltage dynamics of individual cortical neurons in awake behaving mice for more than 30 min using both resonant-scanning and ULoVE random-access microscopy. Finally, ULoVE recording of JEDI-2P can robustly detect spikes at depths exceeding 400 μm and report voltage correlations in pairs of neurons.

Spontaneous Activity Predicts Survival of Developing Cortical Neurons

Spontaneous activity plays a crucial role in brain development by coordinating the integration of immature neurons into emerging cortical networks. High levels and complex patterns of spontaneous activity are generally associated with low rates of apoptosis in the cortex. However, whether spontaneous activity patterns directly encode for survival of individual cortical neurons during development remains an open question. Here, we longitudinally investigated spontaneous activity and apoptosis in developing cortical cultures, combining extracellular electrophysiology with calcium imaging. These experiments demonstrated that the early occurrence of calcium transients was strongly linked to neuronal survival. Silent neurons exhibited a higher probability of cell death, whereas high frequency spiking and burst behavior were almost exclusively detected in surviving neurons. In local neuronal clusters, activity of neighboring neurons exerted a pro-survival effect, whereas on the functional level, networks with a high modular topology were associated with lower cell death rates. Using machine learning algorithms, cell fate of individual neurons was predictable through the integration of spontaneous activity features. Our results indicate that high frequency spiking activity constrains apoptosis in single neurons through sustained calcium rises and thereby consolidates networks in which a high modular topology is reached during early development.

Toward an integrative neurovascular framework for studying brain networks

Brain functional connectivity based on the measure of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals has become one of the most widely used measurements in human neuroimaging. However, the nature of the functional networks revealed by BOLD fMRI can be ambiguous, as highlighted by a recent series of experiments that have suggested that typical resting-state networks can be replicated from purely vascular or physiologically driven BOLD signals. After going through a brief review of the key concepts of brain network analysis, we explore how the vascular and neuronal systems interact to give rise to the brain functional networks measured with BOLD fMRI. This leads us to emphasize a view of the vascular network not only as a confounding element in fMRI but also as a functionally relevant system that is entangled with the neuronal network. To study the vascular and neuronal underpinnings of BOLD functional connectivity, we consider a combination of methodological avenues based on multiscale and multimodal optical imaging in mice, used in combination with computational models that allow the integration of vascular information to explain functional connectivity.

Dynamic Distortion of Orientation Representation after Learning in the Mouse Primary Visual Cortex

Learning is an essential cognitive mechanism allowing behavioral adaptation through adjustments in neuronal processing. It is associated with changes in the activity of sensory cortical neurons evoked by task-relevant stimuli. However, the exact nature of those modifications and the computational advantages they may confer are still debated. Here, we investigated how learning an orientation discrimination task alters the neuronal representations of the cues orientations in the primary visual cortex (V1) of male and female mice. When comparing the activity evoked by the task stimuli in naive mice and the mice performing the task, we found that the representations of the orientation of the rewarded and nonrewarded cues were more accurate and stable in trained mice. This better cue representation in trained mice was associated with a distortion of the orientation representation space such that stimuli flanking the task-relevant orientations were represented as the task stimuli themselves, suggesting that those stimuli were generalized as the task cues. This distortion was context dependent as it was absent in trained mice passively viewing the task cues and enhanced in the behavioral sessions where mice performed best. Those modifications of the V1 population orientation representation in performing mice were supported by a suppression of the activity of neurons tuned for orientations neighboring the orientations of the task cues. Thus, visual processing in V1 is dynamically adapted to enhance the reliability of the representation of the learned cues and favor generalization in the task-relevant computational space.SIGNIFICANCE STATEMENT Performance improvement in a task often requires facilitating the extraction of the information necessary to its execution. Here, we demonstrate the existence of a suppression mechanism that improves the representation of the orientations of the task stimuli in the V1 of mice performing an orientation discrimination task. We also show that this mechanism distorts the V1 orientation representation space, leading stimuli flanking the task stimuli orientations to be generalized as the task stimuli themselves.

Designing AAV Vectors for Monitoring the Subtle Calcium Fluctuations of Inferior Olive Network in vivo

Adeno-associated viral (AAV) vectors, used as vehicles for gene transfer into the brain, are a versatile and powerful tool of modern neuroscience that allow identifying specific neuronal populations, monitoring and modulating their activity. For consistent and reproducible results, the AAV vectors must be engineered so that they reliably and accurately target cell populations. Furthermore, transgene expression must be adjusted to sufficient and safe levels compatible with the physiology of studied cells. We undertook the effort to identify and validate an AAV vector that could be utilized for researching the inferior olivary (IO) nucleus, a structure gating critical timing-related signals to the cerebellum. By means of systematic construct generation and quantitative expression profiling, we succeeded in creating a viral tool for specific and strong transfection of the IO neurons without adverse effects on their physiology. The potential of these tools is demonstrated by expressing the calcium sensor GCaMP6s in adult mouse IO neurons. We could monitor subtle calcium fluctuations underlying two signatures of intrinsic IO activity: the subthreshold oscillations (STOs) and the variable-duration action potential waveforms both in-vitro and in-vivo. Further, we show that the expression levels of GCaMP6s allowing such recordings are compatible with the delicate calcium-based dynamics of IO neurons, inviting future work into the network dynamics of the olivo-cerebellar system in behaving animals.

Striatal neuronal ensembles reveal differential actions of amantadine and clozapine to ameliorate mice L-DOPA-induced dyskinesia

Amantadine and clozapine have proved to reduce abnormal involuntary movements (AIMs) in preclinical and clinical studies of L-DOPA-Induced Dyskinesias (LID). Even though both drugs decrease AIMs, they may have different action mechanisms by using different receptors and signaling profiles. Here we asked whether there are differences in how they modulate neuronal activity of multiple striatal neurons within the striatal microcircuit at histological level during the dose-peak of L-DOPA in ex-vivo brain slices obtained from dyskinetic mice. To answer this question, we used calcium imaging to record the activity of dozens of neurons of the dorsolateral striatum before and after drugs administration in vitro. We also developed an analysis framework to extract encoding insights from calcium imaging data by quantifying neuronal activity, identifying neuronal ensembles by linking neurons that coactivate using hierarchical cluster analysis and extracting network parameters using Graph Theory. The results show that while both drugs reduce LIDs scores behaviorally in a similar way, they have several different and specific actions on modulating the dyskinetic striatal microcircuit. The extracted features were highly accurate in separating amantadine and clozapine effects by means of principal components analysis (PCA) and support vector machine (SVM) algorithms. These results predict possible synergistic actions of amantadine and clozapine on the dyskinetic striatal microcircuit establishing a framework for a bioassay to test novel antidyskinetic drugs or treatments in vitro.

Representing the dynamics of high-dimensional data with non-redundant wavelets

A crucial question in data science is to extract meaningful information embedded in high-dimensional data into a low-dimensional set of features that can represent the original data at different levels. Wavelet analysis is a pervasive method for decomposing time-series signals into a few levels with detailed temporal resolution. However, obtained wavelets are intertwined and over-represented across levels for each sample and across different samples within one population. Here, using neuroscience data of simulated spikes, experimental spikes, calcium imaging signals, and human electrocorticography signals, we leveraged conditional mutual information between wavelets for feature selection. The meaningfulness of selected features was verified to decode stimulus or condition with high accuracy yet using only a small set of features. These results provide a new way of wavelet analysis for extracting essential features of the dynamics of spatiotemporal neural data, which then enables to support novel model design of machine learning with representative features.

•

WCMI can extract meaningful information from high-dimensional data

•

Extracted features from neural signals are non-redundant

•

Simple decoders can read out these features with superb accuracy

One of the essential questions in data science is to extract meaningful information from high-dimensional data. A useful approach is to represent data using a few features that maintain the crucial information. The leading property of spatiotemporal data is foremost ever-changing dynamics in time. Wavelet analysis, as a classical method for disentangling time series, can capture temporal dynamics with detail. Here, we leveraged conditional mutual information between wavelets to select a small subset of non-redundant features. We demonstrated the efficiency and effectiveness of features using various types of neuroscience data with different sampling frequencies at the level of the single cell, cell population, and coarse-scale brain activity. Our results shed new insights into representing the dynamics of spatiotemporal data using a few fundamental features extracted by wavelet analysis, which may have wide implications to other types of data with rich temporal dynamics.

Opposite forms of adaptation in mouse visual cortex are controlled by distinct inhibitory microcircuits

Sensory processing in the cortex adapts to the history of stimulation but the mechanisms are not understood. Imaging the primary visual cortex of mice we find here that an increase in stimulus contrast is not followed by a simple decrease in gain of pyramidal cells; as many cells increase gain to improve detection of a subsequent decrease in contrast. Depressing and sensitizing forms of adaptation also occur in different types of interneurons (PV, SST and VIP) and the net effect within individual pyramidal cells reflects the balance of PV inputs, driving depression, and a subset of SST interneurons driving sensitization. Changes in internal state associated with locomotion increase gain across the population of pyramidal cells while maintaining the balance between these opposite forms of plasticity, consistent with activation of both VIP->SST and SST->PV disinhibitory pathways. These results reveal how different inhibitory microcircuits adjust the gain of pyramidal cells signalling changes in stimulus strength.

The authors describe the role of inhibitory microcircuits in the visual cortex of mice in adaptation to contrast. They show how external stimuli and internal state interact to adjust processing in the visual cortex.

Vision-dependent specification of cell types and function in the developing cortex

The role of postnatal experience in sculpting cortical circuitry, while long appreciated, is poorly understood at the level of cell types. We explore this in the mouse primary visual cortex (V1) using single-nucleus RNA sequencing, visual deprivation, genetics, and functional imaging. We find that vision selectively drives the specification of glutamatergic cell types in upper layers (L) (L2/3/4), while deeper-layer glutamatergic, GABAergic, and non-neuronal cell types are established prior to eye opening. L2/3 cell types form an experience-dependent spatial continuum defined by the graded expression of ∼200 genes, including regulators of cell adhesion and synapse formation. One of these genes, Igsf9b, a vision-dependent gene encoding an inhibitory synaptic cell adhesion molecule, is required for the normal development of binocular responses in L2/3. In summary, vision preferentially regulates the development of upper-layer glutamatergic cell types through the regulation of cell-type-specific gene expression programs.

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.

MTNeuro: A Benchmark for Evaluating Representations of Brain Structure Across Multiple Levels of Abstraction

There are multiple scales of abstraction from which we can describe the same image, depending on whether we are focusing on fine-grained details or a more global attribute of the image. In brain mapping, learning to automatically parse images to build representations of both small-scale features (e.g., the presence of cells or blood vessels) and global properties of an image (e.g., which brain region the image comes from) is a crucial and open challenge. However, most existing datasets and benchmarks for neuroanatomy consider only a single downstream task at a time. To bridge this gap, we introduce a new dataset, annotations, and multiple downstream tasks that provide diverse ways to readout information about brain structure and architecture from the same image. Our multi-task neuroimaging benchmark (MTNeuro) is built on volumetric, micrometer-resolution X-ray microtomography images spanning a large thalamocortical section of mouse brain, encompassing multiple cortical and subcortical regions. We generated a number of different prediction challenges and evaluated several supervised and self-supervised models for brain-region prediction and pixel-level semantic segmentation of microstructures. Our experiments not only highlight the rich heterogeneity of this dataset, but also provide insights into how self-supervised approaches can be used to learn representations that capture multiple attributes of a single image and perform well on a variety of downstream tasks. Datasets, code, and pre-trained baseline models are provided at: https://mtneuro.github.io/.

Macrophage calcium reporter mice reveal immune cell communication in vitro and in vivo

Cell communication underlies emergent functions in diverse cell types and tissues. Recent evidence suggests that macrophages are organized in communicating networks, but new tools are needed to quantitatively characterize the resulting cellular conversations. Here, we infer cell communication from spatiotemporal correlations of intracellular calcium dynamics that are non-destructively imaged across cell populations expressing genetically encoded calcium indicators. We describe a hematopoietic calcium reporter mouse (Csf1rCreGCaMP5fl) and a computational analysis pipeline for inferring communication between reporter cells based on "excess synchrony." We observed signals suggestive of cell communication in macrophages treated with immune-stimulatory DNA in vitro and tumor-associated immune cells imaged in a dorsal window chamber model in vivo. Together, the methods described here expand the toolkit for discovery of cell communication events in macrophages and other immune cells.

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.

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.

Vision is required for the formation of binocular neurons prior to the classical critical period

Depth perception emerges from the development of binocular neurons in primary visual cortex. Vision is required for these neurons to acquire their mature responses to visual stimuli. The prevailing view is that vision does not influence binocular circuitry until the onset of the critical period, about a week after eye opening, and that plasticity of visual responses is triggered by increased inhibition. Here, we show that vision is required to form binocular neurons and to improve binocular tuning and matching from eye opening until critical period closure. Enhancing inhibition does not accelerate this process. Vision soon after eye opening improves the tuning properties of binocular neurons by strengthening and sharpening ipsilateral eye cortical responses. This progressively changes the population of neurons in the binocular pool, and this plasticity is sensitive to interocular differences prior to critical period onset. Thus, vision establishes binocular circuitry and guides binocular plasticity from eye opening.

A database and deep learning toolbox for noise-optimized, generalized spike inference from calcium imaging

Inference of action potentials ('spikes') from neuronal calcium signals is complicated by the scarcity of simultaneous measurements of action potentials and calcium signals ('ground truth'). In this study, we compiled a large, diverse ground truth database from publicly available and newly performed recordings in zebrafish and mice covering a broad range of calcium indicators, cell types and signal-to-noise ratios, comprising a total of more than 35 recording hours from 298 neurons. We developed an algorithm for spike inference (termed CASCADE) that is based on supervised deep networks, takes advantage of the ground truth database, infers absolute spike rates and outperforms existing model-based algorithms. To optimize performance for unseen imaging data, CASCADE retrains itself by resampling ground truth data to match the respective sampling rate and noise level; therefore, no parameters need to be adjusted by the user. In addition, we developed systematic performance assessments for unseen data, openly released a resource toolbox and provide a user-friendly cloud-based implementation.

Reinforcing neuron extraction and spike inference in calcium imaging using deep self-supervised denoising

Calcium imaging has transformed neuroscience research by providing a methodology for monitoring the activity of neural circuits with single-cell resolution. However, calcium imaging is inherently susceptible to detection noise, especially when imaging with high frame rate or under low excitation dosage. Here we developed DeepCAD, a self-supervised deep-learning method for spatiotemporal enhancement of calcium imaging data that does not require any high signal-to-noise ratio (SNR) observations. DeepCAD suppresses detection noise and improves the SNR more than tenfold, which reinforces the accuracy of neuron extraction and spike inference and facilitates the functional analysis of neural circuits.

Reconciling functional differences in populations of neurons recorded with two-photon imaging and electrophysiology

Extracellular electrophysiology and two-photon calcium imaging are widely used methods for measuring physiological activity with single-cell resolution across large populations of cortical neurons. While each of these two modalities has distinct advantages and disadvantages, neither provides complete, unbiased information about the underlying neural population. Here, we compare evoked responses in visual cortex recorded in awake mice under highly standardized conditions using either imaging of genetically expressed GCaMP6f or electrophysiology with silicon probes. Across all stimulus conditions tested, we observe a larger fraction of responsive neurons in electrophysiology and higher stimulus selectivity in calcium imaging, which was partially reconciled by applying a spikes-to-calcium forward model to the electrophysiology data. However, the forward model could only reconcile differences in responsiveness when restricted to neurons with low contamination and an event rate above a minimum threshold. This work established how the biases of these two modalities impact functional metrics that are fundamental for characterizing sensory-evoked responses.

Inferring spikes from calcium imaging in dopamine neurons

Calcium imaging has led to discoveries about neural correlates of behavior in subcortical neurons, including dopamine (DA) neurons. However, spike inference methods have not been tested in most populations of subcortical neurons. To address this gap, we simultaneously performed calcium imaging and electrophysiology in DA neurons in brain slices and applied a recently developed spike inference algorithm to the GCaMP fluorescence. This revealed that individual spikes can be inferred accurately in this population. Next, we inferred spikes in vivo from calcium imaging from these neurons during Pavlovian conditioning, as well as during navigation in virtual reality. In both cases, we quantitatively recapitulated previous in vivo electrophysiological observations. Our work provides a validated approach to infer spikes from calcium imaging in DA neurons and implies that aspects of both tonic and phasic spike patterns can be recovered.

Multidimensional population activity in an electrically coupled inhibitory circuit in the cerebellar cortex

Inhibitory neurons orchestrate the activity of excitatory neurons and play key roles in circuit function. Although individual interneurons have been studied extensively, little is known about their properties at the population level. Using random-access 3D two-photon microscopy, we imaged local populations of cerebellar Golgi cells (GoCs), which deliver inhibition to granule cells. We show that population activity is organized into multiple modes during spontaneous behaviors. A slow, network-wide common modulation of GoC activity correlates with the level of whisking and locomotion, while faster (<1 s) differential population activity, arising from spatially mixed heterogeneous GoC responses, encodes more precise information. A biologically detailed GoC circuit model reproduced the common population mode and the dimensionality observed experimentally, but these properties disappeared when electrical coupling was removed. Our results establish that local GoC circuits exhibit multidimensional activity patterns that could be used for inhibition-mediated adaptive gain control and spatiotemporal patterning of downstream granule cells.

Calcium imaging analysis - how far have we come?

Techniques for calcium imaging were first demonstrated in the mid-1970s, whilst tools to analyse these markers of cellular activity are still being developed and improved today. For image analysis, custom tools were developed within labs and until relatively recently, software packages were not widely available between researchers. We will discuss some of the most popular methods for calcium imaging analysis that are now widely available and describe why these protocols are so effective. We will also describe some of the newest innovations in the field that are likely to benefit researchers, particularly as calcium imaging is often an inherently low signal-to-noise method. Although calcium imaging analysis has seen recent advances, particularly following the rise of machine learning, we will end by highlighting the outstanding requirements and questions that hinder further progress and pose the question of how far we have come in the past sixty years and what can be expected for future development in the field.

Calcium imaging analysis - how far have we come?

Techniques for calcium imaging were first achieved in the mid-1970s, whilst tools to analyse these markers of cellular activity are still being developed and improved. For image analysis, custom tools were developed within labs and until relatively recently, software packages were not widely available between researchers. We will discuss some of the most popular, alongside our preferred, methods for calcium imaging analysis that are now widely available and describe why these protocols are so effective. We will also describe some of the newest innovations in the field that are likely to benefit researchers, particularly as calcium imaging is often an inherently low signal-to-noise method. Although calcium imaging analysis has seen recent advances, particularly following the rise of machine learning, we will end by highlighting the outstanding requirements and questions that hinder further progress, and pose the question of how far we have come in the past sixty years and what can be expected for future development in the field.

MOD: A novel machine-learning optimal-filtering method for accurate and efficient detection of subthreshold synaptic events in vivo

BACKGROUND

To understand information coding in single neurons, it is necessary to analyze subthreshold synaptic events, action potentials (APs), and their interrelation in different behavioral states. However, detecting excitatory postsynaptic potentials (EPSPs) or currents (EPSCs) in behaving animals remains challenging, because of unfavorable signal-to-noise ratio, high frequency, fluctuating amplitude, and variable time course of synaptic events.

NEW METHOD

We developed a method for synaptic event detection, termed MOD (Machine-learning Optimal-filtering Detection-procedure), which combines concepts of supervised machine learning and optimal Wiener filtering. Experts were asked to manually score short epochs of data. The algorithm was trained to obtain the optimal filter coefficients of a Wiener filter and the optimal detection threshold. Scored and unscored data were then processed with the optimal filter, and events were detected as peaks above threshold.

RESULTS

We challenged MOD with EPSP traces in vivo in mice during spatial navigation and EPSC traces in vitro in slices under conditions of enhanced transmitter release. The area under the curve (AUC) of the receiver operating characteristics (ROC) curve was, on average, 0.894 for in vivo and 0.969 for in vitro data sets, indicating high detection accuracy and efficiency.

COMPARISON WITH EXISTING METHODS

When benchmarked using a (1 - AUC)^-1 metric, MOD outperformed previous methods (template-fit, deconvolution, and Bayesian methods) by an average factor of 3.13 for in vivo data sets, but showed comparable (template-fit, deconvolution) or higher (Bayesian) computational efficacy.

CONCLUSIONS

MOD may become an important new tool for large-scale, real-time analysis of synaptic activity.

Relationship between simultaneously recorded spiking activity and fluorescence signal in GCaMP6 transgenic mice

Fluorescent calcium indicators are often used to investigate neural dynamics, but the relationship between fluorescence and action potentials (APs) remains unclear. Most APs can be detected when the soma almost fills the microscope’s field of view, but calcium indicators are used to image populations of neurons, necessitating a large field of view, generating fewer photons per neuron, and compromising AP detection. Here, we characterized the AP-fluorescence transfer function in vivo for 48 layer 2/3 pyramidal neurons in primary visual cortex, with simultaneous calcium imaging and cell-attached recordings from transgenic mice expressing GCaMP6s or GCaMP6f. While most APs were detected under optimal conditions, under conditions typical of population imaging studies, only a minority of 1 AP and 2 AP events were detected (often <10% and ~20–30%, respectively), emphasizing the limits of AP detection under more realistic imaging conditions.

Neurons, the cells that make up the nervous system, transmit information using electrical signals known as action potentials or spikes. Studying the spiking patterns of neurons in the brain is essential to understand perception, memory, thought, and behaviour. One way to do that is by recording electrical activity with microelectrodes. Another way to study neuronal activity is by using molecules that change how they interact with light when calcium binds to them, since changes in calcium concentration can be indicative of neuronal spiking. That change can be observed with specialized microscopes know as two-photon fluorescence microscopes. Using calcium indicators, it is possible to simultaneously record hundreds or even thousands of neurons. However, calcium fluorescence and spikes do not translate one-to-one.

In order to interpret fluorescence data, it is important to understand the relationship between the fluorescence signals and the spikes associated with individual neurons. The only way to directly measure this relationship is by using calcium imaging and electrical recording simultaneously to record activity from the same neuron. However, this is extremely challenging experimentally, so this type of data is rare.

To shed some light on this, Huang, Ledochowitsch et al. used mice that had been genetically modified to produce a calcium indicator in neurons of the visual cortex and simultaneously obtained both fluorescence measurements and electrical recordings from these neurons. These experiments revealed that, while the majority of time periods containing multi-spike neural activity could be identified using calcium imaging microscopy, on average, less than 10% of isolated single spikes were detectable. This is an important caveat that researchers need to take into consideration when interpreting calcium imaging results.

These findings are intended to serve as a guide for interpreting calcium imaging studies that look at neurons in the mammalian brain at the population level. In addition, the data provided will be useful as a reference for the development of activity sensors, and to benchmark and improve computational approaches for detecting and predicting spikes.

Signal-to-signal neural networks for improved spike estimation from calcium imaging data

Spiking information of individual neurons is essential for functional and behavioral analysis in neuroscience research. Calcium imaging techniques are generally employed to obtain activities of neuronal populations. However, these techniques result in slowly-varying fluorescence signals with low temporal resolution. Estimating the temporal positions of the neuronal action potentials from these signals is a challenging problem. In the literature, several generative model-based and data-driven algorithms have been studied with varied levels of success. This article proposes a neural network-based signal-to-signal conversion approach, where it takes as input raw-fluorescence signal and learns to estimate the spike information in an end-to-end fashion. Theoretically, the proposed approach formulates the spike estimation as a single channel source separation problem with unknown mixing conditions. The source corresponding to the action potentials at a lower resolution is estimated at the output. Experimental studies on the spikefinder challenge dataset show that the proposed signal-to-signal conversion approach significantly outperforms state-of-the-art-methods in terms of Pearson's correlation coefficient, Spearman's rank correlation coefficient and yields comparable performance for the area under the receiver operating characteristics measure. We also show that the resulting system: (a) has low complexity with respect to existing supervised approaches and is reproducible; (b) is layer-wise interpretable, and (c) has the capability to generalize across different calcium indicators.

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.

Identification and quantification of neuronal ensembles in optical imaging experiments

Recent technical advances in molecular biology and optical imaging have made it possible to record from up to thousands of densely packed neurons in superficial and deep brain regions in vivo, with cellular subtype specificity and high spatiotemporal fidelity. Such optical neurotechnologies are enabling increasingly fine-scaled studies of neuronal circuits and reliably co-active groups of neurons, so-called ensembles. Neuronal ensembles are thought to constitute the basic functional building blocks of brain systems, potentially exhibiting collective computational properties. While the technical framework of in vivo optical imaging and quantification of neuronal activity follows certain widely held standards, analytical methods for study of neuronal co-activity and ensembles lack consensus and are highly varied across the field. Here we provide a comprehensive step-by-step overview of theoretical, experimental, and analytical considerations for the identification and quantification of neuronal ensemble dynamics in high-resolution in vivo optical imaging studies.

A Data-Driven Approach to Predict and Classify Epileptic Seizures from Brain-Wide Calcium Imaging Video Data

The prediction of epileptic seizures has been an essential problem of epilepsy study. The calcium imaging video data images the whole brain-wide neurons activities with electrical discharge recorded by calcium fluorescence intensity (CFI). In this paper, using the zebrafish's brain-wide calcium image video data, we propose a data-driven approach to effectively detect the systemic change-point, and further predict the epileptic seizures. Our approach includes two phases: offline training and online testing. Specifically, during offline training, we extract features and confirm the existence of systemic change-point, then estimate the ratio of unchanged system duration to interictal period duration. For online testing, we implement a statistical model to estimate the change-point, and then predict the onset of epileptic seizure. The testing results show that our proposed approach could effectively predict the time range of future epileptic seizure. Furthermore, we explore the macroscopic patterns of epileptic and control cases, and extract features based on the pattern difference, then implement and compare the classification performance from four machine learning models. Based on the data structure, we also propose a new method to discretize related features, and combine with hierarchical clustering to better visualize and explain the pattern difference between epileptic and control cases.

Vision Changes the Cellular Composition of Binocular Circuitry during the Critical Period

High acuity stereopsis emerges during an early postnatal critical period when binocular neurons in the primary visual cortex sharpen their receptive field tuning properties. We find that this sharpening is achieved by dismantling the binocular circuit present at critical period onset and building it anew. Longitudinal imaging of receptive field tuning (e.g., orientation selectivity) of thousands of neurons reveals that most binocular neurons present in layer 2/3 at critical period onset are poorly tuned and are rendered monocular. In parallel, new binocular neurons are established by conversion of well-tuned monocular neurons as they gain matched input from the other eye. These improvements in binocular tuning in layer 2/3 are not inherited from layer 4 but are driven by the experience-dependent sharpening of ipsilateral eye responses. Thus, vision builds a new and more sharply tuned binocular circuit in layer 2/3 by cellular exchange and not by refining the original circuit.

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.