Calcium imaging has been widely adopted for its ability to record from large neuronal populations. To summarize the time course of neural activity, dimensionality reduction methods, which have been applied extensively to population spiking activity, may be particularly useful. However, it is unclear if the dimensionality reduction methods applied to spiking activity are appropriate for calcium imaging. We thus carried out a systematic study of design choices based on standard dimensionality reduction methods. We also developed a novel method to perform deconvolution and dimensionality reduction simultaneously (termed CILDS). CILDS most accurately recovered the single-trial, low-dimensional time courses from calcium imaging that would have been recovered from spiking activity. CILDS also outperformed the other methods on calcium imaging recordings from larval zebrafish and mice. More broadly, this study represents a foundation for summarizing calcium imaging recordings of large neuronal populations using dimensionality reduction in diverse experimental settings.

Biophysically accurate multicompartmental models of individual neurons have significantly advanced our understanding of the input-output function of single cells. These models depend on a large number of parameters that are difficult to estimate. In practice, they are often hand-tuned to match measured physiological behaviors, thus raising questions of identifiability and interpretability. We propose a statistical approach to the automatic estimation of various biologically relevant parameters, including 1) the distribution of channel densities, 2) the spatiotemporal pattern of synaptic input, and 3) axial resistances across extended dendrites. Recent experimental advances, notably in voltage-sensitive imaging, motivate us to assume access to: i) the spatiotemporal voltage signal in the dendrite and ii) an approximate description of the channel kinetics of interest. We show here that, given i and ii, parameters 1-3 can be inferred simultaneously by nonnegative linear regression; that this optimization problem possesses a unique solution and is guaranteed to converge despite the large number of parameters and their complex nonlinear interaction; and that standard optimization algorithms efficiently reach this optimum with modest computational and data requirements. We demonstrate that the method leads to accurate estimations on a wide variety of challenging model data sets that include up to about 10(4) parameters (roughly two orders of magnitude more than previously feasible) and describe how the method gives insights into the functional interaction of groups of channels.

The early formation of sensorimotor circuits is essential for survival. While the development and function of exteroceptive circuits and their associated motor pathways are well characterized, far less is known about the circuits that convey viscerosensory inputs to the brain and transmit visceromotor commands from the central nervous system to internal organs. Technical limitations, such as the development of viscerosensory and visceromotor circuits and the invasiveness of procedures required to access them, have hindered studies of their functional development in mammals. Using larval zebrafish-which are genetically accessible and optically transparent-we tracked, , how cardiosensory and cardiomotor neural circuits assemble and begin to function. We uncovered a staged program. First, a minimal efferent circuit suffices for heart-rate control: direct brain-to-heart vagal motor innervation is required, intracardiac neurons are not, and heart rate is governed exclusively by the motor vagus nerve. Within the hindbrain, we functionally localize a vagal premotor population that drives this early efferent control. Second, sympathetic innervation arrives and enhances the dynamics and amplitude of cardiac responses, as neurons in the most anterior sympathetic ganglia acquire the ability to drive cardiac acceleration. These neurons exhibit proportional, integral, and derivative-like relationships to heart rate, consistent with controller motifs that shape gain and dynamics. Third, vagal sensory neurons innervate the heart. Distinct subsets increase activity when heart rate falls or rises, and across spontaneous fluctuations, responses to aversive stimuli, and optogenetically evoked cardiac perturbations, their dynamics are captured by a single canonical temporal kernel with neuron-specific phase offsets, supporting a population code for heart rate. This temporally segregated maturation isolates three experimentally tractable regimes-unidirectional brain-to-heart communication, dual efferent control, and closed-loop control after sensory feedback engages-providing a framework for mechanistic dissection of organism-wide heart-brain circuits.

Calcium imaging with protein-based indicators is widely used to follow neural activity in intact nervous systems, but current protein sensors report neural activity at timescales much slower than electrical signalling and are limited by trade-offs between sensitivity and kinetics. Here we used large-scale screening and structure-guided mutagenesis to develop and optimize several fast and sensitive GCaMP-type indicators. The resulting 'jGCaMP8' sensors, based on the calcium-binding protein calmodulin and a fragment of endothelial nitric oxide synthase, have ultra-fast kinetics (half-rise times of 2 ms) and the highest sensitivity for neural activity reported for a protein-based calcium sensor. jGCaMP8 sensors will allow tracking of large populations of neurons on timescales relevant to neural computation.

Optical recording of intricate molecular dynamics is becoming an indispensable technique for biological studies, accelerated by the development of new or improved biosensors and microscopy technology. This creates major computational challenges to extract and quantify biologically meaningful spatiotemporal patterns embedded within complex and rich data sources, many of which cannot be captured with existing methods. Here, we introduce activity quantification and analysis (AQuA2), a fast, accurate, and versatile data analysis platform built upon advanced machine-learning techniques. It decomposes complex live-imaging-based datasets into elementary signaling events, allowing accurate and unbiased quantification of molecular activities and identification of consensus functional units. We demonstrate applications across a wide range of biosensors, cell types, organs, animal models, microscopy techniques, and imaging approaches. As exemplar findings, we show how AQuA2 identified drug-dependent interactions between neurons and astroglia, as well as distinct sensorimotor signal propagation patterns in the mouse spinal cord.

Preprint: https://doi.org/10.1101/2024.05.02.592259

Understanding how neural circuits give rise to behavior requires comprehensive knowledge of neuronal morphology, connectivity, and function. Atlas platforms play a critical role in enabling the visualization, exploration, and dissemination of such information. Here, we present FishExplorer, an interactive and expandable community platform designed to integrate and analyze multimodal brain data from larval zebrafish. FishExplorer supports datasets acquired through light microscopy (LM), electron microscopy (EM), and X-ray imaging, all co-registered within a unified spatial coordinate system which enables seamless comparison of neuronal morphologies and synaptic connections. To further assist circuit analysis, FishExplorer includes a suite of tools for querying and visualizing connectivity at the whole-brain scale. By integrating data from recent large-scale EM reconstructions (presented in companion studies), FishExplorer enables researchers to validate circuit models, explore wiring principles, and generate new hypotheses. As a continuously evolving resource, FishExplorer is designed to facilitate collaborative discovery and serve the growing needs of the teleost neuroscience community.

Differentiable simulations of optical systems can be combined with deep learning-based reconstruction networks to enable high performance computational imaging via end-to-end (E2E) optimization of both the optical encoder and the deep decoder. This has enabled imaging applications such as 3D localization microscopy, depth estimation, and lensless photography via the optimization of local optical encoders. More challenging computational imaging applications, such as 3D snapshot microscopy which compresses 3D volumes into single 2D images, require a highly non-local optical encoder. We show that existing deep network decoders have a locality bias which prevents the optimization of such highly non-local optical encoders. We address this with a decoder based on a shallow neural network architecture using global kernel Fourier convolutional neural networks (FourierNets). We show that FourierNets surpass existing deep network based decoders at reconstructing photographs captured by the highly non-local DiffuserCam optical encoder. Further, we show that FourierNets enable E2E optimization of highly non-local optical encoders for 3D snapshot microscopy. By combining FourierNets with a large-scale multi-GPU differentiable optical simulation, we are able to optimize non-local optical encoders 170× to 7372× larger than prior state of the art, and demonstrate the potential for ROI-type specific optical encoding with a programmable microscope.

When a behavior repeatedly fails to achieve its goal, animals often give up and become passive, which can be strategic for preserving energy or regrouping between attempts. It is unknown how the brain identifies behavioral failures and mediates this behavioral-state switch. In larval zebrafish swimming in virtual reality, visual feedback can be withheld so that swim attempts fail to trigger expected visual flow. After tens of seconds of such motor futility, animals became passive for similar durations. Whole-brain calcium imaging revealed noradrenergic neurons that responded specifically to failed swim attempts and radial astrocytes whose calcium levels accumulated with increasing numbers of failed attempts. Using cell ablation and optogenetic or chemogenetic activation, we found that noradrenergic neurons progressively activated brainstem radial astrocytes, which then suppressed swimming. Thus, radial astrocytes perform a computation critical for behavior: they accumulate evidence that current actions are ineffective and consequently drive changes in behavioral states.

The ability to measure synaptic connectivity and properties is essential for understanding neuronal circuits. However, existing methods that allow such measurements at cellular resolution are laborious and technically demanding. Here, we describe a system that allows such measurements in a high-throughput way by combining two-photon optogenetics and volumetric Ca2+ imaging with whole-cell recording. We reveal a circuit motif for generating fast undulatory locomotion in zebrafish.

Nonvisual photosensation enables animals to sense light without sight. However, the cellular and molecular mechanisms of nonvisual photobehaviors are poorly understood, especially in vertebrate animals. Here, we describe the photomotor response (PMR), a robust and reproducible series of motor behaviors in zebrafish that is elicited by visual wavelengths of light but does not require the eyes, pineal gland, or other canonical deep-brain photoreceptive organs. Unlike the relatively slow effects of canonical nonvisual pathways, motor circuits are strongly and quickly (seconds) recruited during the PMR behavior. We find that the hindbrain is both necessary and sufficient to drive these behaviors. Using in vivo calcium imaging, we identify a discrete set of neurons within the hindbrain whose responses to light mirror the PMR behavior. Pharmacological inhibition of the visual cycle blocks PMR behaviors, suggesting that opsin-based photoreceptors control this behavior. These data represent the first known light-sensing circuit in the vertebrate hindbrain.