Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 651.
Published online 2017 September 21. doi:  10.1038/s41467-017-00310-3
PMCID: PMC5608914

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.


To survive and thrive, motile organisms use sensory cues to navigate towards environments where they are more likely to avoid predators, obtain food or find mates. Efficient goal-directed locomotion requires closed-loop coordination between motor action and sensory perception. Each movement induces a new sensory signal, which in turn modulates the forthcoming motor output. This mechanism is at play in a number of goal-directed behaviours, in organisms ranging from bacteria1 to nematodes2, 3 and insects4, 5, but also among humans6. Numerous models have been proposed to account for this complex-coordinated motion, but to date no data are available to understand how these behavioural strategies might be implemented at the circuit level in the vertebrate brain.

Here, we take advantage of the accessibility of zebrafish larvae to whole-brain imaging79 to investigate the neural sensorimotor computation underlying phototaxis. This behaviour, which drives the animal towards illuminated regions, is already present at 5 days post-fertilisation and is thus likely to be hard-wired in the larval brain1012. To obtain information regarding the direction of a light source, larvae use two complementary strategies: stereo-visual comparison12 and spatio-temporal sampling13. The first mode relies on the difference in instantaneous perceived intensity at two different angles to infer the direction of the source. In the second mode, the illumination spatial gradient is extracted from two successive samplings obtained before and after a gaze-reorienting movement. This latter approach requires that the visual and motor information are integrated in a timely fashion.

We thus postulated the existence of a central neural circuit that would drive spontaneous gaze shift while integrating unilateral and bilateral changes in illumination so as to bias reorienting bouts towards a light source. We first establish, through behavioural assays, that the gaze and turning-bout orientations are robustly coordinated, and that the statistics of gaze orientation are biased towards illuminated regions. We then use volumetric functional imaging and optogenetic activation to identify a bilaterally distributed neuronal ensemble in the rostral hindbrain that appears to be a strong candidate for the control of spontaneous gaze dynamics. We further investigate how the self-oscillatory dynamics of this population are modulated by bilateral and unilateral changes in illumination. Finally, we develop a rate model of this circuit that reproduces most of our observations and provides the first comprehensive description of how phototaxis can be implemented in the vertebrate brain.


Visually induced modulation of gaze shift dynamics

Larvae can redirect their gaze by triggering coherent angular excursions of both eyes—a process called a saccade—or through whole-body reorientation. We examined the endogenous temporal structure of both processes by monitoring tail deflections and eye orientations in the absence of visual cues in larvae partially restrained in agarose. Spontaneous ocular saccades in zebrafish larvae are highly stereotyped (Fig. 1a–d; Supplementary Fig. 1), with a quasi-periodic alternation between leftward and rightward gaze periods. This property is reflected by the quasi-bimodal distribution of gaze orientations and a distribution of delays between successive reorientations peaked at ~12 s (12.2 ± 3.5 s, mean ± std, N = 29). The orientations of the gaze and the turning bouts were found to be strongly correlated (Fig. 1e–i; Supplementary Fig. 1 and Supplementary Movie 1). We observed that over 85% of the tail-beats that could be unequivocally classified as an attempted turning bout were oriented in the direction of the gaze. Moreover, a large fraction (60.2%) of these bouts occurred within the first second following the onset of an ipsiversive saccade. This coordination mechanism was true for both large-amplitude reorienting saccades that drove the right–left alternation of gaze direction, and for low-amplitude secondary saccades that maintained the gaze to the right or to the left by compensating for slow ocular drift (Supplementarty Fig. 1d–g).

Fig. 1
Behavioural assays of ocular-saccade–turning-bout coordination and light-induced gaze bias. a Definition of the eye and tail kinematic parameters. The larva’s tail and eyes are free. The gaze angle, θ gaze, is defined as the mean ...

These observations suggest that a unique command circuit may drive ocular saccades and set the direction of successive tail-beats, such that the gaze sequence may be viewed as a proxy for fish reorienting dynamics. We thus posited that phototaxis should manifest as a bias in the animal gaze-angle distribution towards illuminated regions. We tested this hypothesis by replicating, in agarose-restrained animals, two visual stimulation protocols known to evoke phototaxis in freely swimming larvae while monitoring gaze dynamics12, 13. Visual stimulation was obtained by projecting two LEDs onto a screen below the animal, each delivering uniform illumination to one eye’s field of view (Supplementary Methods).

We first submitted the larva to a 1 min period of whole-field (bilateral) illumination followed by a step decrement of the intensity to one eye, whereas illumination to the other eye was held constant (Fig. 1j; Supplementary Fig. 2). The unilateral illumination period was maintained for 30 s, and the sequence was then repeated. In 67% of the fish (18 of 27), we observed that the gaze orientation was significantly biased towards the more illuminated side during unilateral illumination periods (Fig. 1k). This process was reflected in the trial-averaged post-extinction gaze sequence, which exhibited a transient drift toward the illuminated region (Supplementary Fig. 2f). This bias was associated with a 27 ± 3.4% (mean ± std) increase in fixation time towards the illuminated area. In total, a significant phototactic behaviour was observed in 77% of the fish, yet 10% displayed negative phototaxis, in accord with the experiments performed in freely swimming configurations12. We confirmed that, owing to gaze–tail coordination, this light-induced gaze bias resulted in a turning bias of the animal toward the light by performing similar experiments while simultaneously monitoring both the eyes and tail movements (Supplementary Fig. 2k–m).

In a second approach, we developed a virtual-reality phototaxis assay in which the animal had access to spatial cues through spatio-temporal probing only. Both eyes were exposed to the same whole-field illumination at an intensity that was locked in real-time to the gaze orientation (Fig. 1l; Supplementary Fig. 2i–j). The illumination was thus at its maximum when the gaze was oriented to one side, and at its minimum for the other side. A significant light-induced bias of the gaze orientation was observed in 85% of the fish (p < 0.05, N = 11 of 13). All responsive larvae showed an increased fixation time towards the brighter direction, consistent with a positive phototactic behaviour (Fig. 1m).

In both assays, the imposed spatial light gradient did not trigger a systematic and immediate gaze shift. Instead, it induced a statistical bias of the spontaneous saccadic dynamics whose net result was to increase the relative duration of gaze fixations towards the brightest region.

Functional mapping of gaze-tuned neuronal populations

These observations suggested that the saccadic dynamics play a central role in phototaxis by enabling a sequential sampling of illumination gradients through spontaneous gaze shifts, while driving the reorienting dynamics via gaze–tail coordination. We thus sought to identify the saccade command circuit. Ocular saccades in zebrafish have been studied mostly as a model of neural integration1416. The accessibility of this vertebrate model to optogenetic techniques allows one to examine with unmatched precision how a transient neural command that rapidly brings the eyes into a desired position can be transformed into a sustained signal to maintain a fixed eye position17. However, in contrast with mammals18, the upstream command circuit that drives the rhythmic alternation between leftward and rightward saccades in zebrafish remains elusive19.

Our approach to identify the saccade command neuronal population consisted of recording large fractions of the brain using light-sheet functional imaging in animals expressing calcium reporter in nearly every neuron (Tg(elavl3:GCaMP6f)) while monitoring eye dynamics (Fig. 2a; Supplementary Movie 2). We restricted our analysis to the hindbrain and the caudal midbrain, where saccade-induced distortions of the brain tissues remained small enough to enable consistent signal extraction. We implemented a multilinear regression approach to classify individual voxels based on their tuning with respect to the angular position and velocity signals of the eyes20 (Supplementary Methods). We then applied a non-rigid registration method to align these functional maps onto a single reference brain obtained by averaging four different samples imaged using the one-photon light-sheet microscope at high spatial resolution. A zebrafish brain atlas21 was finally used to anatomically localise the various neuronal clusters that displayed robust tuning with the gaze position and velocity (Fig. 2b, c; Supplementary Movie 3).

Fig. 2
Regression-based identification of gaze-tuned neuronal populations. a Schematic of the experimental setup and regression analysis. Volumetric recordings on GCaMP6f-expressing larvae were performed using one-photon light-sheet imaging (20 sections per ...

The most prominent velocity-tuned neuronal assemblies were located in two bilaterally symmetric clusters in rhombomeres 4–6 (rh 4–6, Fig. 2d–f, region 1). These populations correspond to the previously identified saccade generator burst neurons22 (SGBN). Each lateral cluster triggers ipsiversive saccades through direct activation of two oculomotor nuclei (abducens and oculomotor nucleus III). Position-tuned neurons were found in the caudal hindbrain in a region consistent with the previously described velocity-position neural integrator20 (VPNI, Fig. 2d–f, region 2). These analysis allowed us to further discover the existence of position-tuned clusters in the rostral hindbrain (rh 2–3) in the form of four stripes bilaterally distributed on each side of the midline (Fig. 2d–f, region 3). These clusters systematically displayed increased and prolonged activity after every ipsiversive saccade, whether reorienting or secondary (Supplementary Fig. 3). We also found two symmetric clusters, lying at the rostral border of the hindbrain (rh 1) ventral to the cerebellum, that were strongly correlated with contraversive gaze orientation (Fig. 2d–f, region 4). Notice that the same gaze-tuned neuronal clusters in the rostral hindbrain were independently identified by Alex Ramirez and Emre Aksay in a yet unpublished study.

Identifying a putative saccade command circuit

The morphology and anatomical location of this newly identified gaze-tuned ensemble in the rostral hindbrain (the four stripes located in rh 2–3) appeared very similar to the self-oscillating circuit identified by Ahrens et al. in eyes-fixed experiments based on activity correlation analysis23 (the so-called hindbrain oscillator, or HBO). Consistent with these findings, we observed that this particular neuronal population displayed sustained antiphasic dynamics, with each side activated in alternation, in recordings where the saccades were abolished through eye fixation in agar (Fig. 3a–g). To properly delineate the neuronal population engaged in these self-generated oscillations, which we refer to as the HBO, we first approximated the oscillatory signal using a small subset of neurons and then computed the tuning of each voxel to this reference trace (Supplementary Methods; Supplementary Fig. 4). This functional criterion yields a neuronal map that essentially encompasses the gaze-tuned circuit in rh 1-3 as formerly identified through regression with the gaze signal (Fig. 3a, b; Supplementary Fig. 5). In more caudal regions of the hindbrain (SGBN and VPNI areas, rh4-7), a small number of neurons also participated in these self-oscillatory dynamics. We analysed the temporal structure of the HBO’s endogenous oscillations in eye-fixed conditions by running the Hilbert transform on the differential signal (left minus right circuit activity, Fig. 3e–g). We found the period distribution to be peaked at ~20 s (19.9 ± 4.8 s mean ± std, N = 8 fish), i.e. close to the spontaneous saccadic period.

Fig. 3
Self-oscillatory dynamics and optogenetically evoked saccades. a Dorso-ventral projection of the gaze-tuned (top-half, N = 7 fish) and self-oscillatory hindbrain population (bottom-half, N = 8 eye-fixed fish). In the latter, ...

The quasi-periodic nature of the saccadic dynamics suggests the existence of a circuit, akin to a central pattern generator, capable of maintaining oscillations with a period on the order of 20 s in the absence of rhythmic sensory or proprioceptive inputs. The HBO constitutes the best candidate for this saccade command circuit. This putative role was found to be consistent with optogenetic assays using larvae expressing ChR2 pan-neurally (Fig. 3h–l; Supplementary Fig. 6). We delivered series of 2.75 s-long light pulses over 47 × 57 μm2 brain areas, alternately on either side of the midline with a 18 s period (one pulse every 9 s), while monitoring saccades (Fig. 3j). Each pair of bilaterally symmetric regions was tested five times before the targeted regions were relocated to eventually probe a large fraction of the hindbrain encompassing both the HBO and the SGBN. As previously reported22, activation of the SGBN (rh 4–6) consistently entrained the saccadic dynamics (8 fish of 19, see Supplementary Methods for details on the quantification of the response). For these responsive fish, activation of the hindbrain in the rh 2–3 region evoked ipsiversive saccades with comparable efficiency (Fig. 3l). These responses could not result from direct activation of the oculomotor nuclei, which lie immediately rostral to rh 1 but drive contraversive saccades. In contrast, optogenetic activation of the two most rostral regions of the oscillator, located in rh 1, failed to evoke conjugated saccades, which might indicate that these two assemblies lie downstream of the HBO’s gaze-driving neurons. Finally, a small area in the caudal midbrain was found to trigger ipsiversive saccades in two larvae.

Careful examination of the functional imaging data in eyes-free preparations offers hints regarding a putative connectivity map between the HBO and the SGBN. We identified two small clusters in rh 7 systematically activated 2–3 s prior to the saccade onset (Supplementary Fig. 3a–c), which are thus analogous to the so-called long-lead burst neurons in the mammalian saccadic circuit24. Such a functional trait can be understood within the assumption that they receive inhibitory inputs from the contralateral HBO circuit, whose spike-rate consistently decays before the saccade onset. In contrast, the ipsilateral HBO circuit remained silent until the saccade onset, which precludes the possibility that it drives, via excitatory inputs, these long-lead burst neurons. Our observations further indicate that the HBO may receive ascending ipsilateral signals from saccadic pre-motor centres (efference copy). When two successive saccades occurred in the same direction, the second one was systematically followed by a rebound of the HBO active module (Supplementary Fig. 3d, e). This process cannot result from the reciprocal inhibitory coupling with the contralateral HBO module25, as the latter remains silent during such phases. Although more work will be needed to definitely establish that the HBO drives spontaneous saccades, these various observations suggest a neural architecture of the saccadic command circuit as sketched in Supplementary Fig. 3f.

The hindbrain oscillator responds to visual stimuli

Dunn et al.25 recently reported that the HBO left and right activity biased the orientation of fictive turning bouts to the right and the left, respectively. This finding is in line with our hypothesis that a unique command circuit controls both turning-bout orientation and spontaneous ocular saccades. The resulting gaze–tail coordination process, analogous to the well-documented head–eye coordination in vertebrates26, 27, explains why the orientation of successive turning bouts exhibits a ~10 s autocorrelation decay in freely swimming larvae13, as this value corresponds to the typical period between successive reorienting saccades.

To evaluate the HBO’s putative role in driving the phototactic behaviour, we sought to investigate whether its endogenous dynamics could be visually biased. We thus successively probed its response to asymmetric (unilateral illumination of one eye) and symmetric (similar illumination on both eyes) visual stimuli (Fig. 4). To prevent visual perturbations by the light-sheet, these experiments were carried out using two-photon excitation, i.e. with an infrared source that is invisible to the fish28. In a first series of experiments, we alternately illuminated each eye for 15 s while performing brain-wide functional imaging. The associated response maps indicated that the most prominent visually sensitive clusters in the hindbrain corresponded to the HBO (Fig. 4b–f; Supplementary Fig. 5b). The right HBO, whose activity is associated with a rightward reorienting bias, was activated when the right eye was stimulated, and vice-versa. These unilateral visual responses are consistent with a stereo-visual positive phototaxis mechanism. In 1 of 11 fish, however, we found a reverse response, i.e. the visual stimulation elicited activity in the contralateral HBO. This observation parallels, and may explain, the 10% of fish found to display negative phototaxis in the behavioural assays.

Fig. 4
Response of the HBO to asymmetric and symmetric visual stimuli. a Schematic of the two-photon light-sheet imaging setup with stereo-visual stimulation. b Comparison of the gaze-tuned and visual response projection maps. Top: the red and blue colours encode ...

Phototaxis has been shown to implicate OFF retinal ganglion cells, i.e. neurons activated upon light decrement12. We thus asked whether the visually evoked transient of the left or right HBO was driven by the stimulus onset (light-on) on the ipsilateral eye, the stimulus offset (light-off) on the contralateral eye, or both. To answer this question, we monitored the HBO dynamics while exposing one eye to alternating on–off illumination, the other being maintained at constant illumination. We found that light-on and light-off stimuli induced a sharp increase in activity in the ipsilateral and contralateral sub-populations, respectively. This indicates that the ON and OFF pathways both contribute to driving the HBO (Supplementary Fig. 7, N = 8 fish).

In a second experiment, we investigated the HBO response to bilateral illumination by delivering a series of 100 ms-long flashes simultaneously to both eyes with a 10 s inter-flash period (Fig. 4g, h). Both left and right HBO displayed a significant response as evidenced by the trial-averaged peri-stimulus signal (Fig. 4i, n = 1311 flashes, N = 12 fish, Supplementary Fig. 7d, e), but the across-trials dispersion appeared surprisingly large (coefficient of variation 170%). The visual responses of the left and right HBO were strongly anti-correlated (r = −0.5) such that a given flash seemed to stochastically evoke a response in either one of them (Fig. 4j). To evaluate whether this apparent stochasticity may reflect a state-dependent sensitivity of the HBO, we examined how the stimulus-evoked transient depended on the particular phase of the oscillation at which the stimulus was delivered. We used the Hilbert transform to extract the phase [var phi](t) of the HBO oscillation. We then computed the derivative of the left and right HBO signals at all times corresponding to a given phase, separating those at which a stimulus had been delivered from the times when the oscillation was running free. We then subtracted the means of both groups to extract the sole contribution of the stimulus to the HBO transient dynamics. This analysis yielded a phase-dependent response curve for each left and right HBO (Fig. 4k). It revealed that each subpopulation is responsive to the visual stimulus within a particular phase range and is essentially unresponsive for the rest of the cycle. More specifically, each flash selectively reinforces the cluster that is already active or about to become active after a period of rest. The transition between left–right HBO activity being associated with gaze shift events, such a phase-dependent response curve can be interpreted as a motor-related gating mechanism. In this view, each subpopulation selectively integrates light increments generated by ipsiversive saccades as expected during spatio-temporal sampling.

HBO entrainment by periodic uniform visual stimuli

In biological and physical oscillators, phase-dependent responses can be revealed through phase-locking and frequency entrainment by periodic forcing2931. Consistently, we observed that the HBO oscillation exhibited a mean period of 20 s in the symmetric stimulation experiments, i.e. precisely twice the stimulation period T stim = 10 s (Fig. 4l, T osc = 20.5 ± 1.4 s, N = 12 fish). Furthermore, the phase offset between the HBO and the stimulus signal, defined as δφ(t) = φ(t) - π*tTstim, displayed a marked bimodality, with two peaks located at 0 and π (Fig. 4m, Rayleigh test, p < 0.01, Supplementary Fig. 8a). This indicates that visual stimuli induce a phase-shift in the HBO oscillations such that the flashes tend to coincide with successive right-to-left and left-to-right transitions (Fig. 4m inset).

We suspected that this period-doubling mechanism ensued from the particular choice of the stimulation period, which is close to half the HBO endogenous period T endo. We sought to explore how this synchronization process varied with the stimulation T stim. To guarantee that the mean illumination was independent of the stimulation period, we exposed the fish to a series of alternating ON–OFF illuminations at six frequencies (20 periods each, 10 s < T stim < 60 s) separated by 100 s-long periods of constant illumination (Fig. 5a). For each sequence, we extracted the mean oscillation period T osc and the phase-offset distribution (Fig. 5b, c). In all tested visually responsive fish (N = 5 of 7), this analysis revealed the existence of an entrainment plateau at intermediate stimulation periods (20 s < T stim < 40 s) for which the HBO frequency was controlled by the stimulation period (T osc/T stim [congruent with] 1) and the phase-offset distribution was strongly non-uniform (Fig. 5d–f, Rayleigh test, p < 0.01). At higher and lower frequencies, the entrainment ratio T osc/T stim was close to 2 and 0.5, respectively. The same experiments were replicated for asymmetric stimuli and yielded qualitatively similar results, albeit with a more pronounced phase-locking and a more extended entrainment plateau (Fig. 5g–i, N = 11, Supplementary Fig. 8).

Fig. 5
Phase-locking and frequency entrainment of the HBO by periodic visual stimuli. a Example HBO signals during symmetric visual stimulus presentation consisting of alternating light-ON and light-OFF periods of similar duration. Difference of the left and ...

A stochastic neural model of phototaxis

A rate model approach32 was used to model the oscillatory dynamics of the HBO and its visually evoked responses based on the connectivity architecture suggested by our experimental observations (Fig. 6a). We described the HBO as a half-centred oscillator consisting of two symmetric modules with recurrent excitation, reciprocal inhibition, and adaptation currents (Supplementary Methods). Such a simple model, derived from a spiking, conductance-based model, has been proposed for cortical circuits to account for the phenomenology of perceptual bi-stability33, 34. This description is consistent with the neurotransmitter identity of the different clusters that participate in the self-oscillatory dynamics25. The medial clusters of the HBO were identified as glutamatergic and the lateral clusters as primarily GABAergic, reflecting the hindbrain columnar organisation of the neurotransmitter classes35.

Fig. 6
A rate model of the HBO’s visually entrained dynamics leading to phototaxis. a, b Schematic of HBO network connectivity. Visual stimuli are relayed to the HBO via the ON and OFF pathways, as detailed in (b). The HBO activity controls gaze shifts ...

The model parameters were adjusted to best reproduce the self-sustained quasi-periodic oscillations observed in the absence of visual input (Fig. 6c; Supplementary Fig. 9). Visual inputs were assumed to be relayed to the HBO in the simplest way consistent with the proposed connectivity (Fig. 6b). A step increase in the illumination of the right eye was associated with a burst of current in the right circuit, followed by a plateau, whereas a step decrease in the illumination induced a transient current in the left circuit to account for the OFF pathway contribution. We then numerically implemented the stimulation protocols used in the functional imaging experiments. The left and right circuit responded to ipsilateral asymmetric visual stimuli, but we further observed a phase-dependent intensity response curve to symmetric flashes in qualitative agreement with its experimental counterpart (Figs. 4k and 6d, e). The model allowed us to reinterpret the latter result in terms of a phase-response curve (PRC). Within this scheme, the effect of a stimulus is to evoke a phase offset, i.e. a transient slowing down or speeding up of the oscillatory dynamics. For each phase, we thus computed the visually induced phase delay Δ[var phi] by comparing the phase evolution after stimulation and in the free-running (without stimulation) regime (Fig. 6f). The resulting π-periodic PRC provides a direct interpretation of the phase-locking mechanism around the phases π/2 and 3π/2, where the PRC displays a negative slope30. Consistently, a phase-locking mechanism was observed with an entrainment plateau around the endogenous oscillation frequency (Fig. 6g). However, for this process to quantitatively compare with the experimental data, the current noise had to be set at a relatively low value (Supplementary Fig. 10).

We sought to evaluate whether this minimal neural model could account for the orientational phototactic behaviour under the assumption that an imbalance in the HBO left vs. right activity distribution may be viewed as a proxy for the statistical orientational bias in the animal locomotion25. We considered configurations in which a distant light source was oriented at an angle θ s with respect to the animal’s rostro-caudal axis. The total illumination impinging on each eye was taken as a slowly varying function of θ s, reflecting the large visual field of view of the zebrafish, with a maximum at θ s = 80° for the right eye (Fig. 6h; Supplementary Methods). We successively considered the stereo-visual comparison and spatio-temporal sampling phototactic strategies. In the first case, we examined the sole effect of a continuous unbalanced bilateral illumination. For each angle θ s, the illumination on each eye was converted into a constant input current on the ipsilateral subcircuit. We then numerically computed the HBO dynamics and compared the mean left vs. right activity. We found the HBO oscillations to be statistically biased towards the subnetwork ipsilateral to the light source. This was assessed by comparing both the time fraction during which one circuit was more active than the other (Fig. 6h) and the mean activities of the two circuits (Supplementary Fig. 11). Within the assumption that such an imbalance reflects a reorienting bias25, this result appears consistent with a positive phototactic mechanism.

In a second approach, we focused on the spatio-temporal sampling strategy (Fig. 6i). We imposed identical illumination on both eyes at all times so as to provide no explicit spatial cue to the circuit. The stimulus intensity was assumed to be a slowly decreasing function of the gaze orientation relative to the light source. We posited that each left-to-right or right-to-left transition of the HBO corresponded to a rightward or leftward gaze shift of amplitude 30°, respectively. This ±15° gaze oscillation in turn elicited a series of abrupt changes in the perceived stimulus. This procedure yielded a left vs. right activity bias that depended on the light source azimuth (Fig. 6i; Supplementary Fig. 11) and was again consistent with a positive phototactic behaviour.


A central goal in neuroscience is to decipher the neural computation underlying the execution of natural tasks. In this respect, whole-brain functional imaging constitutes a game-changing technique because it allows one to systematically map entire neuronal populations that respond to a given sensory stimulus7, 28, 36 or correlate with a particular motor pattern20, 25. Here, this strategy was taken one step further: we used three distinct functional properties as joint criteria to identify a sensorimotor hub likely to play a key role for phototaxis. On the basis of behavioural assays, we reasoned that there should exist a circuit that both sets the pace of spontaneous saccades and integrates visual stimuli. We thus successively delineated three neuronal populations that (i) correlated with the gaze signal during spontaneous saccades, (ii) displayed self-sustained oscillating activity in eye-fixed configurations, and (iii) responded to unilateral visual stimulation. We then used a morphological registration algorithm to demonstrate the existence of a well-defined intersecting neuronal ensemble in the rostral hindbrain, identified as the HBO, that combines all three functionalities (Supplementary Movie 4).

Using a rate model, we showed that the HBO’s phase-dependent visual sensitivity suffices to explain the tendency of larvae to reorient towards illuminated areas using both stereo-visual comparison and spatio-temporal sampling. The proposed sensorimotor processing leading to phototaxis can be summarised as follows: (i) the HBO self-sustained oscillations drive a quasi-periodic sequence of gaze shifts that allow the animal to actively explore the light angular gradient; (ii) asymmetric visual stimuli enhance the activity of the HBO clusters ipsilateral to the eye that receives the larger illumination, which in turn biases the fish reorientation towards the light source; and (iii) a gaze shift-induced light increment reinforces the active HBO subpopulation (ipsilateral to the new gaze direction), thus delaying further left–right transition. Over multiple cycles, this process yields a statistical bias of the fish turning probability towards the brightest region.

In foveate species, ocular saccades enable detailed visual exploration by sequentially bringing different regions of the visual scene onto a small area of the retina (the fovea) where the resolving power is maximal37. The present study suggests that saccades in afoveate species such as zebrafish may subserve a non-image forming visual function by allowing the animal to sequentially sample the light level at two distinct angles without the need to execute a whole-body reorientation. Within this scheme, the large stereotyped gaze shifts that characterize spontaneous saccades in zebrafish might be understood as a way to maximise the detected contrast in the context of weak illumination gradients.

The HBO presents striking similarities with central pattern generators (CPG). These specialized networks drive cyclic motor patterns underlying rhythmic behaviours such as locomotion, chewing, or breathing38. Although they can operate in the absence of external sensory or proprioceptive inputs, the latter have an important role in regulating the circuit dynamics to produce behaviourally relevant motor output39. Their action generally varies with the particular phase within the cycle at which they occur40. The HBO differs from standard CPGs in the sense that it drives a pseudo-periodic motor pattern—the saccadic period displays significant variations from cycle to cycle—but also controls the sequence of discrete turning-bout orientations. By modulating the rhythmic HBO oscillations, the visual inputs thus modulate the temporal pattern of reorienting movements, effectively biasing the animal’s trajectory towards brightest areas over multiple cycles. It has been recently demonstrated, in Drosophila larvae, that sensorimotor neurons can be engaged in distinct goal-directed behaviours41. It will thus be interesting to see whether other sensory information—chemical cues, water flow, temperature gradients, etc.—may also be relayed to the HBO, and how these multi-sensory cues are weighted in order for the larva to navigate towards the most favourable environments.

Beyond phototaxis, the sensorimotor processing operated by the HBO may also be examined in the broader context of active sensing. In vertebrates, including mammals, sensory perception—olfaction, touch, vision, etc.—generally rely on conjugated motor routines42, 43 such as sniffing, whisking, eye saccades, etc. The rhythms of these motor patterns have been shown to entrain low-frequency neuronal oscillations in primary sensory cortices, which in turn modulate their sensitivity according to specific motor phases44. This sensorimotor coordination has a key role in sensory processing. Given the generic nature of this neural mechanism, it will be exciting to further investigate whether the HBO oscillation may similarly drive, through afferent projections, rhythmic changes in sensitivity in sensory centres such that the olfactory bulb or the optic tectum.

Data availability

All data and codes used for the analysis are available from the authors on request.

Electronic supplementary material


We are very grateful to Claire Wyart for allowing us to use her optogenetic setup. We would like to thank the IBPS fish facility staff, and in particular Stéphane Tronche and Alex Bois for their invaluable help. We are also grateful to Carounagarane Dore for his important contribution to the design of the light-sheet microscope. This work was supported by C’Nano Ile de France, École des Neurosciences de Paris Ile-de-France (ENP), Fondation pour la Recherche Médicale (FRM: FDT20160435670 and FDT20150532639) and Fondation Pierre-Gilles de Gennes (FPGG).

Author contributions

Author contributions

S.W., A.D., V.B., R.C., R.M. and G.D. designed the project. S.W., T.B., V.B., R.C., S.K. and G.D. carried out the behavioural and neuroimaging experiments and analysis. D.G.C.H. and I.H.B. generated the elavl3:GCaMP6f transgenic zebrafish line. A.D. and R.M. developed the model and performed data analysis. S.W. and U.L.B. performed the optogenetic assays. All the authors contributed to writing the manuscript.


Competing interests

The authors declare no competing financial interests.


Sébastien Wolf and Alexis M. Dubreuil contributed equally to this work.

Electronic supplementary material

Supplementary Information accompanies this paper at doi:10.1038/s41467-017-00310-3.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


1. Berg HC, Brown DA. Chemotaxis in Escherichia coli analysed by three-dimensional tracking. Nature. 1972;239:500–504. doi: 10.1038/239500a0. [PubMed] [Cross Ref]
2. Chalasani SH, et al. Dissecting a circuit for olfactory behaviour in Caenorhabditis elegans. Nature. 2007;450:63–70. doi: 10.1038/nature06292. [PubMed] [Cross Ref]
3. Ward A, Liu J, Feng Z, Xu XZS. Light-sensitive neurons and channels mediate phototaxis in C. elegans. Nat. Neurosci. 2008;11:916–922. doi: 10.1038/nn.2155. [PMC free article] [PubMed] [Cross Ref]
4. Gomez-Marin A, Stephens GJ, Louis M. Active sampling and decision making in Drosophila chemotaxis. Nat. Commun. 2011;2:441. doi: 10.1038/ncomms1455. [PMC free article] [PubMed] [Cross Ref]
5. Gepner R, et al. Computations underlying Drosophila photo-taxis, odor-taxis, and multi-sensory integration. Elife. 2015;4:e06229. doi: 10.7554/eLife.06229. [PMC free article] [PubMed] [Cross Ref]
6. Porter J, et al. Mechanisms of scent-tracking in humans. Nat. Neurosci. 2007;10:27–29. doi: 10.1038/nn1819. [PubMed] [Cross Ref]
7. Ahrens MB, et al. Brain-wide neuronal dynamics during motor adaptation in zebrafish. Nature. 2012;485:471–477. [PMC free article] [PubMed]
8. Panier T, et al. Fast functional imaging of multiple brain regions in intact zebrafish larvae using selective plane illumination microscopy. Front. Neural Circuits. 2013;7:65. doi: 10.3389/fncir.2013.00065. [PMC free article] [PubMed] [Cross Ref]
9. Feierstein CE, Portugues R, Orger MB. Seeing the whole picture: a comprehensive imaging approach to functional mapping of circuits in behaving zebrafish. Neuroscience. 2015;296:26–38. doi: 10.1016/j.neuroscience.2014.11.046. [PubMed] [Cross Ref]
10. Brockerhoff SE, et al. A behavioral screen for isolating zebrafish mutants with visual system defects. Proc. Natl Acad. Sci. USA. 1995;92:10545–10549. doi: 10.1073/pnas.92.23.10545. [PubMed] [Cross Ref]
11. Orger MB, Baier H. Channeling of red and green cone inputs to the zebrafish optomotor response. Vis. Neurosci. 2005;22:275–281. doi: 10.1017/S0952523805223039. [PubMed] [Cross Ref]
12. Burgess HA, Schoch H, Granato M. Distinct retinal pathways drive spatial orientation behaviors in zebrafish navigation. Curr. Biol. 2010;20:381–386. doi: 10.1016/j.cub.2010.01.022. [PMC free article] [PubMed] [Cross Ref]
13. Chen X, Engert F. Navigational strategies underlying phototaxis in larval zebrafish. Front. Syst. Neurosci. 2014;8:39. doi: 10.3389/fnsys.2014.00039. [PMC free article] [PubMed] [Cross Ref]
14. Miri A, et al. Spatial gradients and multidimensional dynamics in a neural integrator circuit. Nat. Neurosci. 2011;14:1150–1159. doi: 10.1038/nn.2888. [PMC free article] [PubMed] [Cross Ref]
15. Daie K, Goldman MS, Aksay ERF. Spatial patterns of persistent neural activity vary with the behavioral context of short-term memory. Neuron. 2015;85:847–860. doi: 10.1016/j.neuron.2015.01.006. [PMC free article] [PubMed] [Cross Ref]
16. Lee MM, Arrenberg AB, Aksay ERF. A structural and genotypic scaffold underlying temporal integration. J. Neurosci. 2015;35:7903–7920. doi: 10.1523/JNEUROSCI.3045-14.2015. [PMC free article] [PubMed] [Cross Ref]
17. Joshua M, Lisberger SG. A tale of two species: neural integration in zebrafish and monkeys. Neuroscience. 2015;296:80–91. doi: 10.1016/j.neuroscience.2014.04.048. [PMC free article] [PubMed] [Cross Ref]
18. Sparks DL. The brainstem control of saccadic eye movements. Nat. Rev. Neurosci. 2002;3:952–964. doi: 10.1038/nrn986. [PubMed] [Cross Ref]
19. Gonçalves PJ, Arrenberg AB, Hablitzel B, Baier H, Machens CK. Optogenetic perturbations reveal the dynamics of an oculomotor integrator. Front. Neural Circuits. 2014;8:10. doi: 10.3389/fncir.2014.00010. [PMC free article] [PubMed] [Cross Ref]
20. Miri A, Daie K, Burdine RD, Aksay E, Tank DW. Regression-based identification of behavior-encoding neurons during large-scale optical imaging of neural activity at cellular resolution. J. Neurophysiol. 2011;105:964–980. doi: 10.1152/jn.00702.2010. [PMC free article] [PubMed] [Cross Ref]
21. Randlett O, et al. Whole-brain activity mapping onto a zebrafish brain atlas. Nat. Methods. 2015;12:1039–1046. doi: 10.1038/nmeth.3581. [PMC free article] [PubMed] [Cross Ref]
22. Schoonheim PJ, Arrenberg AB, Del Bene F, Baier H. Optogenetic localization and genetic perturbation of saccade-generating neurons in zebrafish. J. Neurosci. 2010;30:7111–7120. doi: 10.1523/JNEUROSCI.5193-09.2010. [PMC free article] [PubMed] [Cross Ref]
23. Ahrens MB, et al. Whole-brain functional imaging at cellular resolution using light-sheet microscopy. Nat. Methods. 2013;10:413–420. doi: 10.1038/nmeth.2434. [PubMed] [Cross Ref]
24. Scudder CA, Kaneko CS, Fuchs AF. The brainstem burst generator for saccadic eye movements: a modern synthesis. Exp. Brain Res. 2002;142:439–462. doi: 10.1007/s00221-001-0912-9. [PubMed] [Cross Ref]
25. Dunn TW, et al. Brain-wide mapping of neural activity controlling zebrafish exploratory locomotion. Elife. 2016;5:e12741. doi: 10.7554/eLife.12741. [PMC free article] [PubMed] [Cross Ref]
26. Berthoz A, Grantyn A. Neuronal mechanisms underlying eye-head coordination. Prog. Brain Res. 1986;64:325–343. doi: 10.1016/S0079-6123(08)63427-5. [PubMed] [Cross Ref]
27. Einhäuser W, et al. Eye-head coordination during free exploration in human and cat. Ann. N. Y. Acad. Sci. 2009;1164:353–366. doi: 10.1111/j.1749-6632.2008.03709.x. [PubMed] [Cross Ref]
28. Wolf S, et al. Whole-brain functional imaging with two-photon light-sheet microscopy. Nat. Methods. 2015;12:379–380. doi: 10.1038/nmeth.3371. [PubMed] [Cross Ref]
29. Brown TG. On the nature of the fundamental activity of the nervous centres; together with an analysis of the conditioning of rhythmic activity in progression, and a theory of the evolution of function in the nervous system. J. Physiol. 1914;48:18–46. doi: 10.1113/jphysiol.1914.sp001646. [PubMed] [Cross Ref]
30. Winfree AT. Biological rhythms and the behavior of populations of coupled oscillators. J. Theor. Biol. 1967;16:15–42. doi: 10.1016/0022-5193(67)90051-3. [PubMed] [Cross Ref]
31. Roenneberg T, Dragovic Z, Merrow M. Demasking biological oscillators: properties and principles of entrainment exemplified by the Neurospora circadian clock. Proc. Natl Acad. Sci. USA. 2005;102:7742–7747. doi: 10.1073/pnas.0501884102. [PubMed] [Cross Ref]
32. Dayan, P. & Abbott, L. F. Theoretical Neuroscience (MIT Press, 2001).
33. Laing CR, Chow CC. A spiking neuron model for binocular rivalry. J. Comput. Neurosci. 2002;12:39–53. doi: 10.1023/A:1014942129705. [PubMed] [Cross Ref]
34. Vattikuti S, et al. Canonical cortical circuit model explains rivalry, intermittent rivalry, and rivalry memory. PLoS Comput. Biol. 2016;12:e1004903. doi: 10.1371/journal.pcbi.1004903. [PMC free article] [PubMed] [Cross Ref]
35. Kinkhabwala A, et al. A structural and functional ground plan for neurons in the hindbrain of zebrafish. Proc. Natl Acad. Sci. 2011;108:1164–1169. doi: 10.1073/pnas.1012185108. [PubMed] [Cross Ref]
36. Portugues R, Feierstein CE, Engert F, Orger MB. Whole-brain activity maps reveal stereotyped, distributed networks for visuomotor behavior. Neuron. 2014;81:6. doi: 10.1016/j.neuron.2014.01.019. [PMC free article] [PubMed] [Cross Ref]
37. Martinez-Conde S, Macknik SL. Fixational eye movements across vertebrates: comparative dynamics, physiology, and perception. J. Vis. 2008;8:28. doi: 10.1167/8.14.28. [PubMed] [Cross Ref]
38. Marder E, Calabrese RL. Principles of rhythmic motor pattern generation. Physiol. Rev. 1996;76:687–717. [PubMed]
39. Blitz DM, Nusbaum MP. Neural circuit flexibility in a small sensorimotor system. Curr. Opin. Neurobiol. 2011;21:544–552. doi: 10.1016/j.conb.2011.05.019. [PMC free article] [PubMed] [Cross Ref]
40. Pearson KG. Generating the walking gait: role of sensory feedback. Prog. Brain Res. 2004;143:123–129. doi: 10.1016/S0079-6123(03)43012-4. [PubMed] [Cross Ref]
41. Tastekin I, et al. Role of the subesophageal zone in sensorimotor control of orientation in Drosophila larva. Curr. Biol. 2015;25:1448–1460. doi: 10.1016/j.cub.2015.04.016. [PubMed] [Cross Ref]
42. Schroeder CE, Wilson DA, Radman T, Scharfman H, Lakatos P. Dynamics of active sensing and perceptual selection. Curr. Opin. Neurobiol. 2010;20:172–176. doi: 10.1016/j.conb.2010.02.010. [PMC free article] [PubMed] [Cross Ref]
43. Ranade S, Hangya B, Kepecs A. Multiple modes of phase locking between sniffing and whisking during active exploration. J. Neurosci. 2013;33:8250–8256. doi: 10.1523/JNEUROSCI.3874-12.2013. [PMC free article] [PubMed] [Cross Ref]
44. Ito J, Maldonado P, Singer W, Grün S. Saccade-related modulations of neuronal excitability support synchrony of visually elicited spikes. Cereb. Cortex. 2011;21:2482–2497. doi: 10.1093/cercor/bhr020. [PMC free article] [PubMed] [Cross Ref]

Articles from Nature Communications are provided here courtesy of Nature Publishing Group