Magnetic actuation of otoliths allows behavioral and brain-wide neuronal exploration of vestibulo-motor processing in larval zebrafish.

The vestibular system in the inner ear plays a central role in sensorimotor control by informing the brain about the orientation and acceleration of the head. However, most experiments in neurophysiology are performed using head-fixed configurations, depriving animals of vestibular inputs. To overcome this limitation, we decorated the utricular otolith of the vestibular system in larval zebrafish with paramagnetic nanoparticles. This procedure effectively endowed the animal with magneto-sensitive capacities: applied magnetic field gradients induced forces on the otoliths, resulting in robust behavioral responses comparable to those evoked by rotating the animal by up to 25°. We recorded the whole-brain neuronal response to this fictive motion stimulation using light-sheet functional imaging. Experiments performed in unilaterally injected fish revealed the activation of a commissural inhibition between the brain hemispheres. This magnetic-based stimulation technique for larval zebrafish opens new perspectives to functionally dissect the neural circuits underlying vestibular processing and to develop multisensory virtual environments, including vestibular feedback.

Emergence of time persistence in a data-driven neural network model.

Establishing accurate as well as interpretable models of network activity is an open challenge in systems neuroscience. Here, we infer an energy-based model of the anterior rhombencephalic turning region (ARTR), a circuit that controls zebrafish swimming statistics, using functional recordings of the spontaneous activity of hundreds of neurons. Although our model is trained to reproduce the low-order statistics of the network activity at short time scales, its simulated dynamics quantitatively captures the slowly alternating activity of the ARTR. It further reproduces the modulation of this persistent dynamics by the water temperature and visual stimulation. Mathematical analysis of the model unveils a low-dimensional landscape-based representation of the ARTR activity, where the slow network dynamics reflects Arrhenius-like barriers crossings between metastable states. Our work thus shows how data-driven models built from large neural populations recordings can be reduced to low-dimensional functional models in order to reveal the fundamental mechanisms controlling the collective neuronal dynamics.

Neural assemblies uncovered by generative modeling explain whole-brain activity statistics and reflect structural connectivity.

Patterns of endogenous activity in the brain reflect a stochastic exploration of the neuronal state space that is constrained by the underlying assembly organization of neurons. Yet, it remains to be shown that this interplay between neurons and their assembly dynamics indeed suffices to generate whole-brain data statistics. Here, we recorded the activity from ∼40,000 neurons simultaneously in zebrafish larvae, and show that a data-driven generative model of neuron-assembly interactions can accurately reproduce the mean activity and pairwise correlation statistics of their spontaneous activity. This model, the compositional Restricted Boltzmann Machine (cRBM), unveils ∼200 neural assemblies, which compose neurophysiological circuits and whose various combinations form successive brain states. We then performed in silico perturbation experiments to determine the interregional functional connectivity, which is conserved across individual animals and correlates well with structural connectivity. Our results showcase how cRBMs can capture the coarse-grained organization of the zebrafish brain. Notably, this generative model can readily be deployed to parse neural data obtained by other large-scale recording techniques.

A scalable assay for chemical preference of small freshwater fish.

Sensing the chemical world is of primary importance for aquatic organisms, and small freshwater fish are increasingly used in toxicology, ethology, and neuroscience by virtue of their ease of manipulation, tissue imaging amenability, and genetic tractability. However, precise behavioral analyses are generally challenging to perform due to the lack of knowledge of what chemical the fish are exposed to at any given moment. Here we developed a behavioral assay and a specific infrared dye to probe the preference of young zebrafish for virtually any compound. We found that the innate aversion of zebrafish to citric acid is not mediated by modulation of the swim but rather by immediate avoidance reactions when the product is sensed and that the preference of juvenile zebrafish for ATP changes from repulsion to attraction during successive exposures. We propose an information-based behavioral model for which an exploration index emerges as a relevant behavioral descriptor, complementary to the standard preference index. Our setup features a high versatility in protocols and is automatic and scalable, which paves the way for high-throughput preference compound screening at different ages.

Evolutionary divergence of locomotion in two related vertebrate species.

Locomotion exists in diverse forms in nature; however, little is known about how closely related species with similar neuronal circuitry can evolve different navigational strategies to explore their environments. Here, we investigate this question by comparing divergent swimming pattern in larval Danionella cerebrum (DC) and zebrafish (ZF). We show that DC displays long continuous swimming events when compared with the short burst-and-glide swimming in ZF. We reveal that mesencephalic locomotion maintenance neurons in the midbrain are sufficient to cause this increased swimming. Moreover, we propose that the availability of dissolved oxygen and timing of swim bladder inflation drive the observed differences in the swim pattern. Our findings uncover the neural substrate underlying the evolutionary divergence of locomotion and its adaptation to their environmental constraints.

An analysis pipeline to compare explorative locomotion across fish species.

Recently, we introduced a powerful approach that leverages differences in swimming behaviors of two closely related fish species to identify previously unreported locomotion-related neuronal correlates. Here, we present this analysis approach applicable for any species of fish to compare their short and long timescale swimming kinematics. We describe steps for data collection and cleaning, followed by the calculation of short timescale kinematics using half tail beats and the analysis of long timescale kinematics using mean square displacement and heading decorrelation. For complete details on the use and execution of this protocol, please refer to Rajan et al. (2022).1.

Thermal modulation of Zebrafish exploratory statistics reveals constraints on individual behavioral variability.

BACKGROUND: Variability is a hallmark of animal behavior. It contributes to survival by endowing individuals and populations with the capacity to adapt to ever-changing environmental conditions. Intra-individual variability is thought to reflect both endogenous and exogenous modulations of the neural dynamics of the central nervous system. However, how variability is internally regulated and modulated by external cues remains elusive. Here, we address this question by analyzing the statistics of spontaneous exploration of freely swimming zebrafish larvae and by probing how these locomotor patterns are impacted when changing the water temperatures within an ethologically relevant range. RESULTS: We show that, for this simple animal model, five short-term kinematic parameters - interbout interval, turn amplitude, travelled distance, turn probability, and orientational flipping rate - together control the long-term exploratory dynamics. We establish that the bath temperature consistently impacts the means of these parameters, but leave their pairwise covariance unchanged. These results indicate that the temperature merely controls the sampling statistics within a well-defined kinematic space delineated by this robust statistical structure. At a given temperature, individual animals explore the behavioral space over a timescale of tens of minutes, suggestive of a slow internal state modulation that could be externally biased through the bath temperature. By combining these various observations into a minimal stochastic model of navigation, we show that this thermal modulation of locomotor kinematics results in a thermophobic behavior, complementing direct gradient-sensing mechanisms. CONCLUSIONS: This study establishes the existence of a well-defined locomotor space accessible to zebrafish larvae during spontaneous exploration, and quantifies self-generated modulation of locomotor patterns. Intra-individual variability reflects a slow diffusive-like probing of this space by the animal. The bath temperature in turn restricts the sampling statistics to sub-regions, endowing the animal with basic thermophobicity. This study suggests that in zebrafish, as well as in other ectothermic animals, ambient temperature could be used to efficiently manipulate internal states in a simple and ethological way.

Blind deconvolution for spike inference from fluorescence recordings.

The parallel developments of genetically-encoded calcium indicators and fast fluorescence imaging techniques allows one to simultaneously record neural activity of extended neuronal populations in vivo. To fully harness the potential of functional imaging, one needs to infer the sequence of action potentials from fluorescence traces. Here we build on recently proposed computational approaches to develop a blind sparse deconvolution (BSD) algorithm based on a generative model for inferring spike trains from fluorescence traces. BSD features, (1) automatic (fully unsupervised) estimation of the hyperparameters, such as spike amplitude, noise level and rise and decay time constants, (2) a novel analytical estimate of the sparsity prior, which yields enhanced robustness and computational speed with respect to existing methods, (3) automatic thresholding for binarizing spikes that maximizes the precision-recall performance, (4) super-resolution capabilities increasing the temporal resolution beyond the fluorescence signal acquisition rate. BSD also uniquely provides theoretically-grounded estimates of the expected performance of the spike reconstruction in terms of precision-recall and temporal accuracy for each recording. The performance of the algorithm is established using synthetic data and through the SpikeFinder challenge, a community-based initiative for spike-rate inference benchmarking based on a collection of joint electrophysiological and fluorescence recordings. Our method outperforms classical sparse deconvolution algorithms in terms of robustness, speed and/or accuracy and performs competitively in the SpikeFinder challenge. This algorithm is modular, easy-to-use and made freely available. Its novel features can thus be incorporated in a straightforward way into existing calcium imaging packages.

From behavior to circuit modeling of light-seeking navigation in zebrafish larvae.

Bridging brain-scale circuit dynamics and organism-scale behavior is a central challenge in neuroscience. It requires the concurrent development of minimal behavioral and neural circuit models that can quantitatively capture basic sensorimotor operations. Here, we focus on light-seeking navigation in zebrafish larvae. Using a virtual reality assay, we first characterize how motor and visual stimulation sequences govern the selection of discrete swim-bout events that subserve the fish navigation in the presence of a distant light source. These mechanisms are combined into a comprehensive Markov-chain model of navigation that quantitatively predicts the stationary distribution of the fish's body orientation under any given illumination profile. We then map this behavioral description onto a neuronal model of the ARTR, a small neural circuit involved in the orientation-selection of swim bouts. We demonstrate that this visually-biased decision-making circuit can capture the statistics of both spontaneous and contrast-driven navigation.

Whole-Brain Calcium Imaging during Physiological Vestibular Stimulation in Larval Zebrafish.

The vestibular apparatus provides animals with postural and movement-related information that is essential to adequately execute numerous sensorimotor tasks. In order to activate this sensory system in a physiological manner, one needs to macroscopically rotate or translate the animal's head, which in turn renders simultaneous neural recordings highly challenging. Here we report on a novel miniaturized, light-sheet microscope that can be dynamically co-rotated with a head-restrained zebrafish larva, enabling controlled vestibular stimulation. The mechanical rigidity of the microscope allows one to perform whole-brain functional imaging with state-of-the-art resolution and signal-to-noise ratio while imposing up to 25° in angular position and 6,000°/s2 in rotational acceleration. We illustrate the potential of this novel setup by producing the first whole-brain response maps to sinusoidal and stepwise vestibular stimulation. The responsive population spans multiple brain areas and displays bilateral symmetry, and its organization is highly stereotypic across individuals. Using Fourier and regression analysis, we identified three major functional clusters that exhibit well-defined phasic and tonic response patterns to vestibular stimulation. Our rotatable light-sheet microscope provides a unique tool for systematically studying vestibular processing in the vertebrate brain and extends the potential of virtual-reality systems to explore complex multisensory and motor integration during simulated 3D navigation.

Sensorimotor computation underlying phototaxis in zebrafish.

Animals continuously gather sensory cues to move towards favourable environments. Efficient goal-directed navigation requires sensory perception and motor commands to be intertwined in a feedback loop, yet the neural substrate underlying this sensorimotor task in the vertebrate brain remains elusive. Here, we combine virtual-reality behavioural assays, volumetric calcium imaging, optogenetic stimulation and circuit modelling to reveal the neural mechanisms through which a zebrafish performs phototaxis, i.e. actively orients towards a light source. Key to this process is a self-oscillating hindbrain population (HBO) that acts as a pacemaker for ocular saccades and controls the orientation of successive swim-bouts. It further integrates visual stimuli in a state-dependent manner, i.e. its response to visual inputs varies with the motor context, a mechanism that manifests itself in the phase-locked entrainment of the HBO by periodic stimuli. A rate model is developed that reproduces our observations and demonstrates how this sensorimotor processing eventually biases the animal trajectory towards bright regions.Active locomotion requires closed-loop sensorimotor co ordination between perception and action. Here the authors show using behavioural, imaging and modelling approaches that gaze orientation during phototaxis behaviour in larval zebrafish is related to oscillatory dynamics of a neuronal population in the hindbrain.

Rheotaxis of Larval Zebrafish: Behavioral Study of a Multi-Sensory Process.

Awake animals unceasingly perceive sensory inputs with great variability of nature and intensity, and understanding how the nervous system manages this continuous flow of diverse information to get a coherent representation of the environment is arguably a central question in systems neuroscience. Rheotaxis, the ability shared by most aquatic species to orient toward a current and swim to hold position, is an innate and robust multi-sensory behavior that is known to involve the lateral line and visual systems. To facilitate the neuroethological study of rheotaxic behavior in larval zebrafish we developed an assay for freely swimming larvae that allows for high experimental throughtput, large statistic and a fine description of the behavior. We show that there exist a clear transition from exploration to counterflow swim, and by changing the sensory modalities accessible to the fishes (visual only, lateral line only or both) and comparing the swim patterns at different ages we were able to detect and characterize two different mechanisms for position holding, one mediated by the lateral line and one mediated by the visual system. We also found that when both sensory modalities are accessible the visual system overshadows the lateral line, suggesting that at the larval stage the sensory inputs are not merged to finely tune the behavior but that redundant information pathways may be used as functional fallbacks.

Whisker Contact Detection of Rodents Based on Slow and Fast Mechanical Inputs.

Rodents use their whiskers to locate nearby objects with an extreme precision. To perform such tasks, they need to detect whisker/object contacts with a high temporal accuracy. This contact detection is conveyed by classes of mechanoreceptors whose neural activity is sensitive to either slow or fast time varying mechanical stresses acting at the base of the whiskers. We developed a biomimetic approach to separate and characterize slow quasi-static and fast vibrational stress signals acting on a whisker base in realistic exploratory phases, using experiments on both real and artificial whiskers. Both slow and fast mechanical inputs are successfully captured using a mechanical model of the whisker. We present and discuss consequences of the whisking process in purely mechanical terms and hypothesize that free whisking in air sets a mechanical threshold for contact detection. The time resolution and robustness of the contact detection strategies based on either slow or fast stress signals are determined. Contact detection based on the vibrational signal is faster and more robust to exploratory conditions than the slow quasi-static component, although both slow/fast components allow localizing the object.

A microfluidic device to study neuronal and motor responses to acute chemical stimuli in zebrafish.

Zebrafish larva is a unique model for whole-brain functional imaging and to study sensory-motor integration in the vertebrate brain. To take full advantage of this system, one needs to design sensory environments that can mimic the complex spatiotemporal stimulus patterns experienced by the animal in natural conditions. We report on a novel open-ended microfluidic device that delivers pulses of chemical stimuli to agarose-restrained larvae with near-millisecond switching rate and unprecedented spatial and concentration accuracy and reproducibility. In combination with two-photon calcium imaging and recordings of tail movements, we found that stimuli of opposite hedonic values induced different circuit activity patterns. Moreover, by precisely controlling the duration of the stimulus (50-500 ms), we found that the probability of generating a gustatory-induced behavior is encoded by the number of neurons activated. This device may open new ways to dissect the neural-circuit principles underlying chemosensory perception.

An amplitude modulation/demodulation scheme for whisker-based texture perception.

Whisking rodents can discriminate finely textured objects using their vibrissae. The biomechanical and neural processes underlying such sensory tasks remain elusive. Here we combine the use of model micropatterned substrates and high-resolution videography of rats' whiskers during tactile exploration to study how texture information is mechanically encoded in the whisker motion. A biomechanical modeling of the whisker is developed, which yields quantitative predictions of the spectral and temporal characteristics of the observed whisker kinetics, for any given topography. These texture-induced whisker vibrations are then replayed via a multiwhisker stimulator while recording neuronal responses in the barrel field of the primary somatosensory cortex (S1bf). These results provide a comprehensive description of the transduction process at play during fine texture sensing in rats. They suggest that the sensory system operates through a vibratory amplitude modulation/demodulation scheme. Fine textural properties are encoded in the time-varying envelope of the whisker-resonant vibrations. This quantity is then recovered by neural demodulation, as it effectively drives the spiking-rate signal of a large fraction of S1 cortical neurons. This encoding/decoding scheme is shown to be robust against variations in exploratory conditions, such as the scanning speed or pad-to-substrate distance, thus allowing for reliable tactile discrimination in realistic conditions.

Fast functional imaging of multiple brain regions in intact zebrafish larvae using selective plane illumination microscopy.

The optical transparency and the small dimensions of zebrafish at the larval stage make it a vertebrate model of choice for brain-wide in-vivo functional imaging. However, current point-scanning imaging techniques, such as two-photon or confocal microscopy, impose a strong limit on acquisition speed which in turn sets the number of neurons that can be simultaneously recorded. At 5 Hz, this number is of the order of one thousand, i.e., approximately 1-2% of the brain. Here we demonstrate that this limitation can be greatly overcome by using Selective-plane Illumination Microscopy (SPIM). Zebrafish larvae expressing the genetically encoded calcium indicator GCaMP3 were illuminated with a scanned laser sheet and imaged with a camera whose optical axis was oriented orthogonally to the illumination plane. This optical sectioning approach was shown to permit functional imaging of a very large fraction of the brain volume of 5-9-day-old larvae with single- or near single-cell resolution. The spontaneous activity of up to 5,000 neurons was recorded at 20 Hz for 20-60 min. By rapidly scanning the specimen in the axial direction, the activity of 25,000 individual neurons from 5 different z-planes (approximately 30% of the entire brain) could be simultaneously monitored at 4 Hz. Compared to point-scanning techniques, this imaging strategy thus yields a ≃20-fold increase in data throughput (number of recorded neurons times acquisition rate) without compromising the signal-to-noise ratio (SNR). The extended field of view offered by the SPIM method allowed us to directly identify large scale ensembles of neurons, spanning several brain regions, that displayed correlated activity and were thus likely to participate in common neural processes. The benefits and limitations of SPIM for functional imaging in zebrafish as well as future developments are briefly discussed.

Whisker encoding of mechanical events during active tactile exploration.

Rats use their whiskers to extract a wealth of information about their immediate environment, such as the shape, position or texture of an object. The information is conveyed to mechanoreceptors located within the whisker follicle in the form of a sequence of whisker deflections induced by the whisker/object contact interaction. How the whiskers filter and shape the mechanical information and effectively participate in the coding of tactile features remains an open question to date. In the present article, a biomechanical model was developed that provides predictions of the whisker dynamics during active tactile exploration, amenable to quantitative experimental comparison. This model is based on a decomposition of the whisker profile into a slow, quasi-static sequence and rapid resonant small-scale vibrations. It was applied to the typical situation of a rat actively whisking across a solid object. Having derived the quasi-static sequence of whisker deformation, the resonant properties of the whisker were analyzed, taking into account the boundary conditions imposed by the whisker/surface contact. We then focused on two elementary mechanical events that are expected to trigger significant neural responses, namely (1) the whisker/object first contact and (2) the whisker detachment from the object. Both events were found to trigger a deflection wave propagating upward to the mystacial pad at constant velocity of ≈3-5 m/s. This yielded a characteristic mechanical signature at the whisker base, in the form of a large peak of negative curvature occurring ≈4 ms after the event has been triggered. The dependence in amplitude and lag of this mechanical signal with the main contextual parameters (such as radial or angular distance) was investigated. The model was validated experimentally by comparing its predictions to high-speed video recordings of shock-induced whisker deflections performed on anesthetized rats. The consequences of these results on possible tactile encoding schemes are briefly discussed.

The role of exploratory conditions in bio-inspired tactile sensing of single topogical features.

We investigate the mechanism of tactile transduction during active exploration of finely textured surfaces using a tactile sensor mimicking the human fingertip. We focus in particular on the role of exploratory conditions in shaping the subcutaneous mechanical signals. The sensor has been designed by integrating a linear array of MEMS micro-force sensors in an elastomer layer. We measure the response of the sensors to the passage of elementary topographical features at constant velocity and normal load, such as a small hole on a flat substrate. Each sensor's response is found to strongly depend on its relative location with respect to the substrate/skin contact zone, a result which can be quantitatively understood within the scope of a linear model of tactile transduction. The modification of the response induced by varying other parameters, such as the thickness of the elastic layer and the confining load, are also correctly captured by this model. We further demonstrate that the knowledge of these characteristic responses allows one to dynamically evaluate the position of a small hole within the contact zone, based on the micro-force sensors signals, with a spatial resolution an order of magnitude better than the intrinsic resolution of individual sensors. Consequences of these observations on robotic tactile sensing are briefly discussed.

Effect of fingerprints orientation on skin vibrations during tactile exploration of textured surfaces.

In humans, the tactile perception of fine textures is mediated by skin vibrations when scanning the surface with the fingertip. These vibrations are encoded by specific mechanoreceptors, Pacinian corpuscules (PCs), located about 2 mm below the skin surface. In a recent article, we performed experiments using a biomimetic sensor which suggest that fingerprints (epidermal ridges) may play an important role in shaping the subcutaneous stress vibrations in a way which facilitates their processing by the PC channel. Here we further test this hypothesis by directly recording the modulations of the fingerpad/substrate friction force induced by scanning an actual fingertip across a textured surface. When the fingerprints are oriented perpendicular to the scanning direction, the spectrum of these modulations shows a pronounced maximum around the frequency v/lambda, where v is the scanning velocity and lambda the fingerprints period. This simple biomechanical result confirms the relevance of our previous finding for human touch.