The brain’s capabilities rely on both the molecular properties of individual cells and their interactions across brain-wide networks. However, relating gene expression to activity in individual neurons across the entire brain remains elusive. Here we developed an experimental-computational platform, WARP, for whole-brain imaging of neuronal activity during behavior, expansion-assisted spatial transcriptomics, and cellular-level registration of these two modalities. Through joint analysis of whole-brain neuronal activity during multiple behaviors, cellular gene expression, and anatomy, we identified functions of molecularly defined populations — including luminance coding in a cckb-pou4f2 midbrain population and task-structured activity in pvalb7-eomesa hippocampal-like neurons — and defined over 2,000 other function-gene-anatomy subpopulations. Analysis of this unprecedented multimodal dataset also revealed that most gene-matched neurons showed stronger activity correlations, highlighting a brain-wide role for gene expression in functional organization. WARP establishes a foundational platform and open-access dataset for cross-experiment discovery, high-throughput function-to-gene mapping, unification of cell biology and systems neuroscience, and scalable circuit modeling at the whole-brain scale.

Animals rapidly adapt to changing circumstances by shifting how they perceive, integrate, and act. Such flexibility is often attributed to transitions between internal states that exert widespread influence across the brain. Yet the mechanisms that drive state transitions and how they reconfigure brainwide computation remain unclear. Larval zebrafish, when actions are rendered futile by decoupling visual flow feedback from swimming in virtual reality, enter a temporary passive, energy-preserving state. In this state, astrocyte calcium levels are elevated, and swim reinitiation requires greater accumulated visual motion. Using whole-brain, cellular-resolution activity imaging, we observed widespread circuit alterations underlying this disengaged state: neuronal visual responses weakened, visual motion integration over time became dramatically leakier, motor inhibition increased, and motor preparation slowed, together suppressing conversion of sensory evidence into action. Astrocyte calcium rose during futile swimming, tracked the emergence and resolution of these brainwide changes, and was both necessary and sufficient to drive them. Thus, astrocytes orchestrate internal states that profoundly reshape neural computations, most powerfully at intermediate integrative processing stages, to meet changing demands.

Animals reprioritize behavioral goals in response to internal physiological states. Using larval zebrafish, we investigated whether engagement with a visuomotor task, the optomotor response (OMR), is coupled to cardiac dynamics. We discovered that threats lead to tachycardia that is synchronized with behavioral suppression. The change in heart rate is represented in the activity of specific neuronal populations. Severing the input to the sympathetic ganglia or ablating the vagus nerve revealed that the threat-related changes to behavioral state do not require interoceptive pathways. Direct tachycardic optopacing of the heart similarly suppressed the OMR response, but by reducing cardiac filling during diastole, thereby impacting oxygen delivery to the CNS. Optopacing also changed the activity of specific brain regions but in neurons distinct from those associated with threat-induced tachycardia. These cardiac function-associated central changes may have relevance to autonomic imbalances in anxiety, stress, and orthostatic disorders.

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 in utero 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, in vivo, 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.

Unveiling the genetic profiles of spatially distinguished cells is an important aspect in many areas of brain research, as the genetic identity contains information about a cell’s physiological properties and internal state. On top of this, knowledge of the genetic details of each cell can reveal structural organization within tissue. As image-based spatial transcriptomics moves toward applications in tissues with dense cellular packing, accurate assignment of detected mRNA transcripts ("spots") to correct segmented cells becomes increasingly difficult, rendering simple methods insufficient with many incorrect assignments to neighboring cells. Here we introduce SpotDMix, a statistical model for assigning spots to cells by modeling spots as coming from a mixture model of distributions matching segmented cell shapes, with assignment probabilities and shape parameters optimized using the Expectation Maximization algorithm. Performance is assessed and compared against several simple methods in various scenarios on both surrogate data and larval zebrafish data. In all tested scenarios SpotDMix outperforms the simple methods on all evaluated metrics, including individual transcript assignment accuracy, total assigned number of spots per cell error and cell type classification. Further, SpotDMix produces a higher degree of exclusivity between genes which are known to not or rarely co-express.

Sleep is regulated by a homeostatic process and associated with an increased arousal threshold, but the genetic and neuronal mechanisms that implement these essential features of sleep remain poorly understood. To address these fundamental questions, we performed a zebrafish genetic screen informed by human genome-wide association studies. We found that mutation of serine/threonine kinase 32a (stk32a) results in increased sleep and impaired sleep homeostasis in both zebrafish and mice, and that stk32a acts downstream of neurotensin signaling and the serotonergic raphe in zebrafish. stk32a mutation reduces phosphorylation of neurofilament proteins, which are co-expressed with stk32a in neurons that regulate motor activity and in lateral line hair cells that detect environmental stimuli, and ablating these cells phenocopies stk32a mutation. Neurotensin signaling inhibits specific sensory and motor populations, and blocks stimulus-evoked responses of neurons that relay sensory information from hair cells to the brain. Our work thus shows that stk32a is an evolutionarily conserved sleep regulator that links neuropeptidergic and neuromodulatory systems to homeostatic sleep drive and changes in arousal threshold, which are implemented through suppression of specific sensory and motor systems.

All cells in an animal collectively ensure, moment-to-moment, the survival of the whole organism in the face of environmental stressors1,2. Physiology seeks to elucidate the intricate network of interactions that sustain life, which often span multiple organs, cell types, and timescales, but a major challenge lies in the inability to simultaneously record time-varying cellular activity throughout the entire body.We developed WHOLISTIC, a method to image second-timescale, time-varying intracellular dynamics across cell-types of the vertebrate body. By advancing and integrating volumetric fluorescence microscopy, machine learning, and pancellular transgenic expression of calcium sensors in transparent young Danio rerio (zebrafish) and adult Danionella, the method enables real-time recording of cellular dynamics across the organism. Calcium is a universal intracellular messenger, with a large array of cellular processes depending on changes in calcium concentration across varying time-scales, making it an ideal proxy of cellular activity3.Using this platform to screen the dynamics of all cells in the body, we discovered unexpected responses of specific cell types to stimuli, such as chondrocyte reactions to cold, meningeal responses to ketamine, and state-dependent activity, such as oscillatory ependymal-cell activity during periods of extended motor quiescence. At the organ scale, the method uncovered pulsating traveling waves along the kidney nephron. At the multi-organ scale, we uncovered muscle synergies and independencies, as well as muscle-organ interactions. Integration with optogenetics allowed us to all-optically determine the causal direction of brain-body interactions. At the whole-organism scale, the method captured the rapid brainstem-controlled redistribution of blood flow across the body.Finally, we advanced Whole-Body Expansion Microscopy4 to provide ground-truth molecular and ultrastructural anatomical context, explaining the spatiotemporal structure of activity captured by WHOLISTIC. Together, these innovations establish a new paradigm for systems biology, bridging cellular and organismal physiology, with broad implications for both fundamental research and drug discovery.

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.

We present a correlated light and electron microscopy (CLEM) dataset from a 7-day-old larval zebrafish, integrating confocal imaging of genetically labeled excitatory (vglut2a) and inhibitory (gad1b) neurons with nanometer-resolution serial section EM. The dataset spans the brain and anterior spinal cord, capturing >180,000 segmented soma, >40,000 molecularly annotated neurons, and 30 million synapses, most of which were classified as excitatory, inhibitory, or modulatory. To characterize the directional flow of activity across the brain, we leverage the synaptic and cell body annotations to compute region-wise input and output drive indices at single cell resolution. We illustrate the dataset’s utility by dissecting and validating circuits in three distinct systems: water flow direction encoding in the lateral line, recurrent excitation and contralateral inhibition in a hindbrain motion integrator, and functionally relevant targeted long-range projections from a tegmental excitatory nucleus, demonstrating that this resource enables rigorous hypothesis testing as well as exploratory-driven circuit analysis. The dataset is integrated into an open-access platform optimized to facilitate community reconstruction and discovery efforts throughout the larval zebrafish brain.

Preprint: https://www.biorxiv.org/content/early/2025/06/15/2025.06.10.658982

As animals adapt to new situations, neuromodulation is a potent way to alter behavior, yet mechanisms by which neuromodulatory nuclei compute during behavior are underexplored. The serotonergic raphe supports motor learning in larval zebrafish by visually detecting distance traveled during swims, encoding action effectiveness, and modulating motor vigor. We found that swimming opens a gate for visual input to cause spiking in serotonergic neurons, enabling encoding of action outcomes and filtering out learning-irrelevant visual signals. Using light-sheet microscopy, voltage sensors, and neurotransmitter/modulator sensors, we tracked millisecond-timescale neuronal input-output computations during behavior. Swim commands initially inhibited serotonergic neurons via GABA, closing the gate to spiking. Immediately after, the gate briefly opened: voltage increased consistent with post-inhibitory rebound, allowing swim-induced visual motion to evoke firing through glutamate, triggering serotonin secretion and modulating motor vigor. Ablating GABAergic neurons impaired raphe coding and motor learning. Thus, serotonergic neuromodulation arises from action-outcome coincidence detection within the raphe, suggesting the existence of similarly fast and precise circuit computations across neuromodulatory nuclei.

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