# The mechanism for directional hearing in fish

Johannes Veith, Thomas Chaigne, Ana Svanidze, Lena Elisa Dressler, Maximilian Hoffmann, Ben Gerhardt, Benjamin Judkewitz

► To cite this version:

Johannes Veith, Thomas Chaigne, Ana Svanidze, Lena Elisa Dressler, Maximilian Hoffmann, et al. The mechanism for directional hearing in fish. Nature, 2024, 631 (8019), pp.118-124. 10.1038/s41586-024-07507-9. hal-04741474

HAL Id: hal-04741474

https://hal.science/hal-04741474v1

Submitted on 17 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# The mechanism for directional hearing in fish

https://doi.org/10.1038/s41586-024-07507-9

Received: 16 January 2024

Accepted: 2 May 2024

Published online: 19 June 2024

Open access

<figure><img src="image_2.png" /></figure>

Check for updates

Johannes Veith<fnref n="1" />,²,⁶, Thomas Chaigne<fnref n="1" />,³,⁶, Ana Svanidze<fnref n="1" />, Lena Elisa Dressler<fnref n="1" />,⁴, Maximilian Hoffmann<fnref n="1" />,⁵, Ben Gerhardt<fnref n="1" />,² & Benjamin Judkewitz<fnref n="1" />,⁹

Locating sound sources such as prey or predators is critical for survival in many vertebrates. Terrestrial vertebrates locate sources by measuring the time delay and intensity difference of sound pressure at each ear<fnref n="1" />⁻⁵. Underwater, however, the physics of sound makes interaural cues very small, suggesting that directional hearing in fish should be nearly impossible⁶. Yet, directional hearing has been confirmed behaviourally, although the mechanisms have remained unknown for decades. Several hypotheses have been proposed to explain this remarkable ability, including the possibility that fish evolved an extreme sensitivity to minute interaural differences or that fish might compare sound pressure with particle motion signals⁷,⁸. However, experimental challenges have long hindered a definitive explanation. Here we empirically test these models in the transparent teleost Danionella cerebrum, one of the smallest vertebrates⁹,¹⁰. By selectively controlling pressure and particle motion, we dissect the sensory algorithm underlying directional acoustic startles. We find that both cues are indispensable for this behaviour and that their relative phase controls its direction. Using micro-computed tomography and optical vibrometry, we further show that D. cerebrum has the sensory structures to implement this mechanism. D. cerebrum shares these structures with more than 15% of living vertebrate species, suggesting a widespread mechanism for inferring sound direction.

Whether it is the song of a whale, the rustling of a mouse or the quiet sneaking of a cat, when sound emanates from a source, it travels through the surrounding medium as an oscillation of motion and pressure. The ability to detect its direction turns hearing into a spatial sense that is pivotal for survival.

Terrestrial vertebrates sense direction by sampling pressure at the ears (that is, two positions far apart; Fig. 1a). In the late 1940s, Lloyd Jeffress speculated that a combination of delay lines and coincidence detector neurons could measure the interaural time difference (ITD)<fnref n="1" />. An implementation of the Jeffress model was later found by Carr and Konishi in their seminal work on the barn owl brainstem², making it one of the most well-known canonical circuits in vertebrates and a textbook example for the convergence of theory, behaviour, biophysics and physiology³,⁴. In addition to the ITD, the difference in the pressure amplitude across the ears, called the interaural level difference (ILD), also carries information about the sound direction. This quantity is amplified by the head, which reflects airborne sound and casts an acoustic shadow onto the far ear⁵.

The vertebrate sense of hearing evolved from a common fish ancestor, and yet these well-established models for directional hearing falter when applied to underwater environments. As sound travels approximately five times faster in water than in air, ITDs are reduced to very small levels (Fig. 1b). As biological tissues have acoustic impedances similar to water, ILDs are substantially smaller, too. This presents a

conundrum: according to prevailing models, fish should not be able to localize sound. Yet, behavioural evidence shows that fishes such as the Atlantic cod⁷,⁸,¹¹⁻¹⁴, the plainfin midshipman¹⁵⁻¹⁹, herring²⁰⁻²³ and goldfish²⁴⁻³⁰ can determine the direction of sound sources.

What cues are available to fish that might enable directional hearing? Fish have two distinct peripheral auditory pathways³¹,³². First, otolithic end organs of the inner ear, in addition to their vestibular function, also act as particle motion sensors for nanometre to micrometre displacements. Owing to the morphology of hair cells, this direct hearing pathway is inherently spatial, but it allows animals to tell only the axis of sound propagation, not the direction—a limitation termed the 180° ambiguity problem of directional hearing⁶. The second, so-called indirect, hearing pathway relies on the swim bladder, which is filled with compressible gas that oscillates in a pressure field. In Otophysi, a large superorder containing about 66% of freshwater fish species³³, a series of bones (Weberian ossicles) transmits this motion to the inner ear³⁴⁻³⁹.

In 1975, Arie Schuijf proposed a model⁷,⁸ suggesting that fish could theoretically deduce the direction of sound if they were able to separately sense and compare its motion and pressure components. Alternative hypotheses include the possibility that fish evolved an extreme sensitivity to minuscule ITDs and ILDs or use their mechanosensory lateral line organ for directional inference through an unknown mechanism. Despite almost a century of careful work on this topic

<footnote n="1">¹Einstein Center for Neurosciences, Charité – Universitätsmedizin Berlin, Berlin, Germany. ²Institut für Biologie, Humboldt-Universität zu Berlin, Berlin, Germany. ³Aix Marseille Univ, CNRS, Centrale Med, Institut Fresnel, Marseille, France. ⁴Present address: Museum für Naturkunde Berlin, Berlin, Germany. ⁵Present address: Rockefeller University, New York, NY, USA. ⁶These authors contributed equally: Johannes Veith, Thomas Chaigne. ⁷e-mail: benjamin.judkewitz@charite.de</footnote>
a

<figure><img src="image_3.png" /></figure>

$$ \Delta t = \sin(\alpha) \cdot 500 \mu s $$

b

<figure><img src="image_4.png" /></figure>

$$ \Delta t = \sin(\alpha) \cdot 1 \mu s $$

c

<figure><img src="image_5.png" /></figure>

Camera: 120 fps

d

<figure><img src="image_6.png" /></figure>

e

<figure><img src="image_7.png" /></figure>

f

<figure><img src="image_8.png" /></figure>

g

<figure><img src="image_9.png" /><figcaption>**Fig. 1** | Sounds elicit a directional startle reflex in *Danionella cerebrum*.</figcaption></figure>

a, Schematic of a pressure wave, arriving at the auditory organs with a detectable ITD in humans. b, ITDs are heavily diminished underwater (value approximated for *D. cerebrum*). c, Behavioural setup (Methods). d, Playback paradigm. Before the experiment, sound pressure and particle acceleration are calibrated at multiple points inside the inner tank (top left, orange crosses; see also Extended Data Fig. 1c–e). Playback is triggered if three conditions are met: the fish swims into the trigger zone (top right, dotted green rectangle), the fish is oriented ≤45° to the y axis, and ≥5 s have passed since the last playback (Methods). e, Startles are detected by a speed threshold after sound playback (Methods; see Extended Data Fig. 4 for details on startle dynamics). Top: centred trajectories after playback for startles (n = 1,415) and non-startles (n = 2,383) across all fish (n = 65). Bottom: average fish speed for startles and non-startles, aligned to sound trigger at t = 0. f, Centred startle trajectories in two sound configurations show a directional escape away from the left (81% of n = 125 startle trials across 58 fish; two-sided binomial test: *P* = 2 × 10⁻¹²) or right (79% of n = 115 startle trials across 56 fish; two-sided binomial test: *P* = 2 × 10⁻¹⁰) speaker. g, Pooled centred trajectories from f with flipped trajectories for the right speaker stimulus summarize the directional escape away from the single speaker (80% of n = 240 startle trials across 63 fish; two-sided binomial test: *P* = 1 × 10⁻²¹). f,g, The heat maps are normalized and smoothed two-dimensional histograms over endpoint positions of the trajectories (grey for single-stimulus data; blue for pooled data). Scale bars, 10 cm (d) and 5 mm (f,g).

(Supplementary Table 1), starting with Karl von Frisch's studies on minnows in 1935⁴⁰, the mechanism of directional hearing in fish is debated to this date.

Empirical tests of the directional hearing mechanism have been fraught with several difficulties. First, reliably controlling sound stimuli underwater is challenging owing to echoes and sound reverberations, as well as near-field effects in small tanks⁴¹⁻⁴⁵. Second, although fish are sensitive to particle motion, many studies control only for pressure, so several authors have urged experimenters to control both quantities⁴¹,⁴⁶,⁴⁷. Third, the sense of hearing, mediated by the inner ear, must be distinguished from the lateral line sense, which is sensitive to low-frequency water flow (<200 Hz)⁴⁸⁻⁵¹. Fourth, behavioural paradigms should exclude the possibility of klinotaxis (that is, a sequential gradient sampling strategy) instead of true directional hearing. Finally, the study of sound transduction to the fish's inner ear³⁶,³⁸,³⁹ is complicated by tissue opacity.

Here we address these challenges and systematically test hypotheses of directional hearing mechanisms using the transparent fish *Danionella cerebrum*, one of the smallest known vertebrates⁹,¹⁰,⁵²⁻⁵⁹. Its small size makes *D. cerebrum* well suited for high-throughput experiments under controlled laboratory conditions. Moreover, *D. cerebrum* communicate acoustically, underlining the importance of hearing for their behaviour. With an inner ear separation of less than 1 mm, *D. cerebrum* put interaural comparison mechanisms to their ultimate test (Fig. 1b).

We find that *D. cerebrum* perform directional startles away from a sound source and that this ability is independent of the lateral line. We then present an extensive set of controlled particle motion and pressure stimuli that lead to differential predictions for directional responses, depending on seven alternative hypotheses for the mechanism of directional hearing. Finally, we carry out laser-scanning vibrometry of auditory structures across the transparent body to determine the physical cues available to *D. cerebrum*. Together, the findings of our experiments lead us to reject all but one of the proposed hypotheses for directional hearing, and they provide strong support for Schuifj's model that fish compare the phase between pressure and motion to tell the direction of a sound source.

## Directional startle responses

To test whether *D. cerebrum* can hear sound direction, we tracked their motion in an aquarium surrounded by underwater speakers (Fig. 1c, Methods and Extended Data Fig. 1a,b). We played back transient sounds (about 12 ms duration, about 0.7 ms rise time, 780 Hz centre frequency; Methods and Extended Data Figs. 2 and 3) and quantified the direction of their startle reflex. The experiment was carried out with one fish at a time, and playback depended on the fish's orientation and position, to test left-right directional hearing and cancel position-dependent echoes (Fig. 1d and Methods).

Shortly after sound onset (within 17 ms or 2 video frames), *D. cerebrum* performed a characteristic startle reflex involving fast sideward displacement (see Fig. 1d,e, Methods and Extended Data Figs. 4 and 5 for startle dynamics, probabilities and habituation). We used the relative displacement along the left-right speaker axis 50 ms after startle onset as a readout of directional response and found that *D. cerebrum* startle away from the speaker, irrespective of the speaker location (left or right; Fig. 1f). Consequently, left and right responses were pooled in a metric for directional escape away from the speaker (Fig. 1g). These single-speaker sound playbacks resulted in directional escapes for 80% of all startles (191 of n = 240 startles across 63 fish with at least one startle, two-sided binomial test: *P* = 1 × 10⁻²¹), an effect present in both sexes (Extended Data Fig. 6a). A replication of the experiment in lateral line-ablated fish showed equivalent results of directional hearing (Extended Data Figs. 7b and 9b). Thus, female and male *D. cerebrum* exhibit directional hearing independent of lateral line function.
<figure><img src="image_10.png" /><figcaption>**Fig. 2** | The relative phase between pressure and particle motion predicts startle direction. **a**, Schematic of *D. cerebrum* hearing apparatus anatomy. **b**, Schuijf's model: a plane wave from the left differs from a plane wave from the right in terms of the phase relationship between pressure and directed particle velocity. **c**, Illustration of stimulus configurations. Plus, minus and arrows inside speaker symbols refer to speaker signal polarity. L, left speaker; R, right speaker; P, only pressure; M, only motion. The schematics of resulting traces illustrate the pressure and particle velocity relationship (see Extended Data Figs. 2 and 3 for actual traces). Speaker schematics show a simplified configuration. Active echo cancellation typically involved three speakers (Supplementary Table 2). **d–f**, Centred startle trajectories after sound playback. All statistical tests are two-sided binomial tests. **d**, Directional escapes away from a single speaker for a positive polarity signal (same data as Fig. 1g) and negative polarity signal (80% escapes, *P* = 7 × 10⁻²³, *n* = 258 startle trials to this sound in 61 fish). **e**, Absence of significant directional bias in the positive polarity condition with only pressure (44% of *n* = 90 to right, not significant; *P* > 0.05) and only particle motion (46% of *n* = 192 to right, not significant; *P* > 0.05). **f**, Single-speaker playbacks pooled over both polarities evoke a directional escape (80% escapes, *P* = 4 × 10⁻⁴³, *n* = 498 startle trials in 64 fish). Selective inversion of pressure polarity by an additional pair of speakers along the orthogonal axis inverts the relative polarity between pressure and particle velocity, which tricks the fish into performing startles approaching the active speaker (67% approach, *P* = 1 × 10⁻⁹, *n* = 331 startle trials in 61 fish). Scale bars, 500 µm (**a**) and 5 mm (**d–f**).</figcaption></figure>

## Physical cues and hypotheses

We then inquired about possible algorithms that may underlie this directional hearing behaviour. As directionality has to be inferred from physical acoustic cues, we gathered those cues potentially available to *D. cerebrum* to list compatible hypotheses for a directional hearing algorithm. Adults have an inner ear separation (about 0.6 mm; Fig. 2a) that is orders of magnitude smaller than the wavelength of the sounds recorded in their natural habitat (≥150 mm for sounds up to 10 kHz). In addition, the mismatch between the characteristic acoustic impedance

of water (1.5 MPa s m⁻¹) and biological tissue (about 1.6–1.7 MPa s m⁻¹)⁶⁰ is small, unlike the approximately 4,000-fold mismatch between biological tissue and air (0.0004 MPa s m⁻¹). Hence, *D. cerebrum* can hardly break left–right symmetry by casting an acoustic shadow. However, close to a monopole sound source, ILDs could occur owing to the steep decay of the sound field with distance (Extended Data Fig. 8a). In *D. cerebrum*, the level drop between both inner ears at 3 cm distance from a monopole sound source can be as large as about 2% for pressure, irrespective of frequency, and about 4% for particle velocity at frequencies below 4.2 kHz. Hence, sound direction could be inferred
through two pressure sensors (P-ILD; hypothesis 1) or two particle motion sensors (M-ILD; hypothesis 2). These two strategies would work best for nearby sounds coming from a direction along the axis of the sensor pair.

The other binaural cue may be time differences in pressure (P-ITD, hypothesis 3) or particle motion (M-ITD, hypothesis 4). The maximal ITD between *D. cerebrum*'s inner ears is 0.4 µs, orders of magnitude smaller than in terrestrial vertebrates⁴. An ITD mechanism would thus point to an extreme sensitivity to minute time differences.

Each inner ear hair cell deflects along a preferred axis. Continuous sinusoidal sounds from opposing directions would stimulate the same hair cells, leading to the 180° ambiguity problem. If, however, all ecologically relevant sounds started with compression, the animal could interpret initial particle motion as motion directed away from the source. This sense would require just a single particle motion sensor for each axis. We call this possibility M-polarity hearing for either positive or negative polarity (hypotheses 5 and 6).

Finally, Schuijf's model for directional hearing resolves the 180° ambiguity problem by using pressure as a reference quantity (hypothesis 7). The idea is most easily illustrated for a plane wave (Fig. 2b), but its applicability is not limited to it (Extended Data Fig. 8b): at the phase of high compression, particles move away from the source, no matter if it started with push or pull. Thus, the pressure signal can act as a reference signal to resolve the 180° ambiguity of particle motion.

To distinguish between these hypotheses, we took advantage of *D. cerebrum*'s small size and precise control over stimuli under laboratory conditions. This allowed us to create stimuli that would normally not occur in nature (for example, sounds with pure pressure and no particle motion component) and to behaviourally dissect the algorithm used by *D. cerebrum* to tell sound direction.

## The directional hearing algorithm

The observed directional startle responses (Fig. 1f,g) are a reaction to naturalistic two-component sound consisting of pressure and particle motion signals (Fig. 2c(i)). To investigate how the isolated components of this sound field affect *D. cerebrum*, we created a pure pressure stimulus and a pure particle motion stimulus. The pure pressure stimulus was realized by driving opposing speakers in phase (Fig. 2c(iv)), and the pure motion stimulus was realized by driving them out of phase (Fig. 2c(v)), creating standing waves with nodes at the animal's location. Our approach additionally took into account echo cancellation and the spatially mapped frequency response of our recording tanks to present controlled and reproducible stimuli (Methods and Extended Data Figs. 2 and 3). We found that either component alone (pure pressure or pure motion) can elicit startles (Extended Data Fig. 7). However, neither component by itself elicited directionally biased responses, suggesting that both are necessary to tell sound direction (Fig. 2e and Extended Data Fig. 7).

The unbiased responses to the pure particle motion stimulus were the first evidence against a directional hearing algorithm based on initial particle motion polarity (M-polarity). To further test this hypothesis, we presented amplitude-inverted waveforms, for which the M-polarity algorithm would predict a reversed startle response. However, *D. cerebrum* performed startles away from the speaker for both polarities (Fig. 2d; pooled single-speaker data in Fig. 2f: 80% escape of *n* = 498 startles across 64 fish, two-sided binomial test: *P* = 4 × 10⁻⁴³). Hence, we ruled out particle motion polarity (M-polarity) as the sole cue for directional hearing.

To test whether it may be the relationship between pressure and motion that determines the startle direction, we selectively inverted the pressure signal while leaving particle motion unchanged (similar to refs. 8,14). As we have seen, activating the left speaker (for example) evokes startles towards the right, away from the speaker, for both positive or negative waveform polarity (configurations shown

in Fig. 2c(i,ii)). By introducing an additional pressure source through in-phase activation of two speakers orthogonal to the left-right axis (Fig. 2c(vi)), we were able to selectively invert pressure, creating stimuli with pressure and particle motion cues akin to a sound originating from the right (Fig. 2c(iii)), despite the right speaker being inactive (Methods). In this 'trick condition', Schuijf's model predicts reversed startle behaviour (that is, 'escape' towards a speaker).

Indeed, *D. cerebrum* could be tricked: following pressure inversion, *D. cerebrum* performed startles towards the active speaker rather than away (Fig. 2f, 67% approach, two-sided binomial test: *P* = 1 × 10⁻⁹, *n* = 331 startle trials in 61 fish). This result held true for both sexes, within individual fish, for both pressure polarities and in lateral line-ablated fish (Extended Data Figs. 6b-d, 7 and 9d).

To check whether any binaural mechanism explains startles towards the speaker in the trick condition, we estimated the amplitudes and signs of P-ITD, M-ITD, P-ILD and M-ILD on the basis of theory and pressure measurements on both sides of the fish (see the section of the Methods entitled Estimation of binaural cues). We found that the signs of M-ITD and M-ILD remain unchanged when creating the trick condition and cannot explain the reversal of *D. cerebrum*'s escape direction. This suggests that neither interaural time nor level comparisons are sufficient to infer sound direction from particle motion.

In the single-speaker experiments, the distance-dependent decay of absolute pressure is 4.4 Pa over 600 µm, which is 2% of the pressure amplitude and potentially large enough to be detected. When we selectively invert pressure to realize the trick condition, we also effectively invert the sign of the level gradient along the horizontal x axis, as well as the sign of the phase delay (see Extended Data Fig. 8d for a geometrical explanation). Therefore, P-ILD and P-ITD are the other mechanisms that are in agreement with the behavioural response of *D. cerebrum*.

In summary, *D. cerebrum* can be tricked into startling towards the speaker rather than away (see also Supplementary Video 1). Among all seven hypotheses considered, this behaviour is consistent with only three (Fig. 4): Schuijf's model, which relies on the phase comparison between pressure and particle motion; and the P-ILD and P-ITD mechanisms, which both rely on sensing pressure level at two positions in space.

To determine which of these remaining three hypotheses is compatible with *D. cerebrum*'s sensory anatomy, we next asked whether *D. cerebrum* possess sensory organs that can detect pressure and particle motion, as required by Schuijf's model, or sensory organs that can detect a difference in pressure amplitude along the azimuth, as required by the P-ILD and P-ITD mechanisms.

## The hearing apparatus

On the basis of micro-computed tomography (micro-CT) and optical vibrometry, we characterized *D. cerebrum*'s auditory organs to narrow down the candidate mechanisms for directional hearing. We started by visualizing the anatomy of *D. cerebrum*'s hearing apparatus with micro-CT imaging of a phosphomolybdic acid-stained sample (Methods). We segmented the main components of the hearing apparatus such as the swim bladder, the labyrinths, the otoliths (lapillus, asteriscus and sagitta) and otolithic end organs (utricle, lagena and saccule), the ossicles of the Weberian apparatus (including tripus and scaphium), and the lymphatic chambers that engulf the lagena and the saccule (Fig. 3a and Supplementary Video 2).

To study the motion of these auditory structures in response to the sound stimuli, we built an imaging vibrometer based on laser-scanning confocal reflectance microscopy: by time-gated acquisition of each pixel with respect to a continuous sinusoidal sound playback, it was possible to infer the relative phase of the motion of auditory structures in two dimensions (Fig. 3b, Methods and Extended Data Fig. 10).

First, we investigated whether the hearing apparatus of *D. cerebrum* is capable of detecting pressure. To this end, we drove two opposite
<figure><img src="image_11.png" /><figcaption>**Fig. 3** | The *D. cerebrum* hearing apparatus can detect pressure and particle motion separately. a, Segmentation of the hearing apparatus based on micro-CT. The fields of view for vibrometry phase maps in c-e are indicated. b, Illustration of the vibrometry method. A laser-scanning confocal reflectance microscope is used to image the motion of auditory structures in a sound field. Each pixel is sampled at four sound phases to record sound-induced motion in the x-y plane. See Methods and Extended Data Fig. 10 for details. c-e, Particle velocity phase (φ) and amplitude (A) maps for motion along the left-right speaker axis; phase colour wheel shown in b, normalized to the amplitude Amax indicated above the maps. c, A pressure stimulus causes anti-phase motion of the tips of the tripus and scaphium along the medial-lateral axis (indicated with red and cyan arrows). d, Along the dorso-ventral axis, a pressure stimulus creates phase-lagged motion of the scaphium and sagitta with the surrounding tissue. e, A particle motion stimulus results in relative motion between the</figcaption></figure>

speakers in phase, thereby creating a 'pressure-only' stimulus. In this
condition, we observed anti-phase motion of the left and right tripus
and scaphium along the medial-lateral axis, periodically compressing
lymphatic space (Fig. 3c,d). This motion stems from swim bladder
compression and expansion in a pressure field. In addition, we observed
rotational motion of the sagitta, part of the saccular end organ that lies
embedded in an adjacent second lymphatic space (Fig. 3d). By contrast,
we did not observe relative motion between the surrounding tissue and
the lapillus or the asteriscus (Extended Data Fig. 11). We concluded that
an indirect pressure sensing pathway exists in *D. cerebrum* and that

lapillus and surrounding tissue along the medial-lateral axis. Pooling over pixels within the region of interest, the displacement amplitude of the surrounding tissue is 1.24 µm ± 0.14 µm, and that of the lapillus is 0.95 µm ± 0.05 µm. At a relative phase of 0.14 π ± 0.05 π, the relative displacement can be as large as 0.56 µm (that is, 45% of the surrounding tissue motion). f, Indirect pathway for pressure sensing: in a pressure field, the swim bladder oscillates, which moves the tripus and the scaphium inwards and outwards. Lymphatic spaces probably couple the motion of the scaphium to the sagitta<sup>73,74</sup>. g, Direct pathway for particle motion sensing: the tissue of the fish couples to the particle motion of water, but the denser lapillus lags in phase and moves with lower amplitude, leading to a relative displacement that could be sensed by utricular hair cells underneath the lapillus. Observations shown in c-e were repeated once, four times and once in other fish, respectively, with similar results. Arrows in f,g indicate direction of motion. Scale bars, 250 µm (a) and 100 µm (c-e).

the saccule may be its main end organ, in agreement with findings in
other fishes³⁸,³⁹,⁶¹,⁶².

Second, we studied particle motion sensing by driving two opposite
speakers in anti-phase. In this condition, with large particle motion but
low-pressure signal, the lapillus, the otolith of the utricular end organ,
moved with a phase lag (about 0.14 π ± 0.05 π) and at a lower amplitude
(about 76% ± 10%) than the surrounding tissue (Fig. 3e). This creates a
relative motion that is expected to stimulate the underlying hair cell
epithelium (direct pathway). We did not detect such relative motion for
the tripus, the sagitta and the asteriscus. Note that we used stimuli with
<figure><img src="image_12.png" /><figcaption>**Fig. 4 | Evidence for Schuijf's model of phase comparison.** We consider seven models for directional hearing that depend on different sensory structures and predict different behaviours (an eighth lateral line (LL)-based mechanism can be ruled out as we observe directional behaviour despite lateral line ablation). Interaural time difference (M-ITD) and level difference (M-ILD) mechanisms based on particle motion can be rejected on the basis of the behavioural data in the trick configuration. A strategy that is based on escaping positive (M-polarity (+)) or negative (M-polarity (-)) initial motion can be rejected as inverted polarity waveforms fail to invert the startle direction. Finally, a mechanism based on sensing pressure level or time differences (P-ILD, P-ITD) – consistent with behavioural data – can be rejected as *D. cerebrum* possess only a single pressure sensor. This leaves Schuijf's model as the one that correctly predicts an inversion of startle direction in the trick configuration and that is based on sensory cues that *D. cerebrum* is able to sense.</figcaption></figure>

particle motion along the mediolateral axis relevant for the left-right startle behaviour. To support directional hearing in three dimensions, *D. cerebrum* may have further direct motion sensing pathways along additional axes, in line with particle motion tuning in saccular<sup>63,64</sup> or lagenar<sup>63,65</sup> afferents in other species. In summary, we did not find two pressure sensors that could detect a pressure level difference (P-ILD) or time difference (P-ITD) along the binaural axis, but rather a single pressure sensor (Fig. 3f) and a set of particle motion sensors (Fig. 3g).

We consider whether *D. cerebrum* could have another pressure sensor that the vibrometry measurements did not detect. Like in other cyprinids, their swim bladder is divided into an anterior and a posterior part. To rule out *D. cerebrum* using a pressure difference between these divisions, we repeated our behavioural analysis for only those startles in which the anterior-posterior axis of the fish was near-orthogonal to the axis of sound presentation, giving equivalent results (Extended Data Fig. 12a). Theoretically, *D. cerebrum* might possess other, unknown pressure sensors to implement P-ILD or P-ITD. However, all known sound pressure sensors in fish are based on compressible gas-filled organs. As gas-filled structures have a high micro-CT contrast, hypothetical pressure sensing organs would have to be either microscopic to evade detection, or based on an unknown principle of sound pressure transduction without compressible gas. Neither of these options is supported by our current knowledge of fish biology and physics of sound.

We therefore reject P-ILD and P-ITD as plausible mechanisms. Instead, the *D. cerebrum* anatomy is well suited for implementing Schuijf's model for directional hearing.

## Discussion

In this study, we report direct experimental evidence for directional hearing in fish and identify its underlying biophysical mechanism. Directional startles were acoustically elicited in male and female *D. cerebrum*, irrespective of lateral line ablation. A selective inversion of

the pressure component tricked *D. cerebrum* into performing a startle towards the active speaker rather than away, consistent with the hypothesis that fish compare pressure and particle motion signals to infer direction (Schuijf's model). On the basis of anatomical and behavioural data, we rejected all known alternative models, including binaural models for directional hearing. Using optical vibrometry, we confirmed the existence of a direct pathway for particle motion sensing and an indirect pathway for pressure sensing in *D. cerebrum*. Hence, the *D. cerebrum* hearing apparatus supports a dual sense of hearing that allows for a comparison between pressure and particle motion. Together, these findings suggest that Schuijf's model is actually implemented in nature (Fig. 4).

Our work builds on a large body of pioneering work on fish hearing, summarized in Supplementary Table 1. To provide an overview of past evidence for directional hearing in fish, we categorized publications into five study types. We indicated to what degree authors controlled acoustic variables and whether they could judge their importance to directional hearing. These previous studies faced a trade-off: they were carried out either in open water, in which reverberations are negligible, but experiments are challenging so that only a few fish were tested<sup>66</sup>, or in the laboratory, in which more fish could be tested but pressure and particle motion were not fully controlled, and the lateral line function could not be ruled out as a near-field sense that aids directional hearing. Here we addressed these limitations by measuring the impulse response of our speakers, actively cancelling reverberations and precisely controlling natural and unnatural stimuli.

Hearing in humans refers to the perception of pressure oscillations. We have shown that *D. cerebrum* instead have a dual sense of hearing comprising pressure and particle motion sensing, which is used for directional hearing. It could be that the acoustic world is much richer for fish than for humans, with acoustic events carrying stereotypic dual pressure-motion signatures<sup>67,68</sup> that may even reveal their distance<sup>69</sup>. As *D. cerebrum* is a vocal species, directional hearing may also have a social function.

Schuijf's model has recently been extended through a proposal that fish compute the time-averaged product of pressure and motion (the acoustic intensity vector)<sup>6</sup>. This theory can account for phonotaxis of plainfin midshipmen towards monopole and dipole sources<sup>17,18</sup>. Future neurophysiological work may test whether *D. cerebrum* implements Schuijf's model this way.

*D. cerebrum* shares the Weberian apparatus with other otophysans, a superorder that includes 66% of living freshwater fish species and 15% of all living vertebrate species<sup>33,70</sup>. Even in fishes lacking the Weberian apparatus, otolithic end organs may still inherit pressure-induced swim bladder motion<sup>38,62,71,72</sup>. Hence, Schuijf's model may have widespread applicability.

## Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-024-07507-9.

1. Jeffress, L. A. A place theory of sound localization. *J. Comp. Physiol. Psychol.* **41**, 35–39 (1948).

2. Carr, C. & Konishi, M. A circuit for detection of interaural time differences in the brain stem of the barn owl. *J. Neurosci.* **10**, 3227–3246 (1990).

3. Konishi, M. Coding of auditory space. *Annu. Rev. Neurosci.* **26**, 31–55 (2003).

4. Knudsen, E. I. & Konishi, M. Mechanisms of sound localization in the barn owl (*Tyto alba*). *J. Comp. Physiol. A* **133**, 13–21 (1979).

5. Kandel, E. R., Koester, J. D., Mack, S. H. & Siegelbaum, S. A. in *Principles of Neural Science* (eds Kandel, E. R. et al.) 651–681 (McGraw Hill, 2021).

6. Sisneros, J. A. & Rogers, P. H. in *Fish Hearing and Bioacoustics* Vol. 877 (ed. Sisneros, J. A.) 121–155 (Springer, 2016).

7. Schuijf, A. Directional hearing of cod (*Gadus morhua*) under approximate free field conditions. *J. Comp. Physiol. A* **98**, 307–332 (1975).
8. Schuijf, A. & Buwalda, R. J. A. On the mechanism of directional hearing in cod (*Gadus morhua* L.). *J. Comp. Physiol.* **98**, 333-343 (1975).

9. Schulze, L. et al. Transparent *Danionella translucida* as a genetically tractable vertebrate brain model. *Nat. Methods* **15**, 977-983 (2018).

10. Britz, R., Conway, K. W. & Rüber, L. The emerging vertebrate model species for neurophysiological studies is *Danionella cerebrum*, new species (Teleostei: Cyprinidae). *Sci. Rep.* **11**, 18942 (2021).

11. Chapman, C. J. Field studies of hearing in teleost fish. *Helgoländer Wiss. Meeresunters.* **24**, 371-390 (1973).

12. Chapman, C. J. & Johnstone, A. D. F. Some auditory discrimination experiments on marine fish. *J. Exp. Biol.* **61**, 521-528 (1974).

13. Hawkins, A. D. & Sand, O. Directional hearing in the median vertical plane by the cod. *J. Comp. Physiol.* **A** **122**, 1-8 (1977).

14. Buwalda, R. J. A., Schuijf, A. & Hawkins, A. D. Discrimination by the cod of sounds from opposing directions. *J. Comp. Physiol.* **A** **150**, 175-184 (1983).

15. Ibara, R. M., Penny, L. T., Ebeling, A. W., Van Dykhuizen, G. & Cailliet, G. *in Predators and Prey in Fishes* Vol. 2 (eds Noakes, D. L. G. et al.) 205-212 (Springer, 1983).

16. McKibben, J. R. & Bass, A. H. Behavioral assessment of acoustic parameters relevant to signal recognition and preference in a vocal fish. *J. Acoust. Soc. Am.* **104**, 3520-3533 (1998).

17. Zeddies, D. G., Fay, R. R., Alderks, P. W., Shaub, K. S. & Sisneros, J. A. Sound source localization by the plainfin midshipman fish. *J. Acoust. Soc. Am.* **127**, 3104-3113 (2010).

18. Zeddies, D. G. et al. Local acoustic particle motion guides sound-source localization behavior in the plainfin midshipman fish, *Porichthys notatus*. *J. Exp. Biol.* **215**, 152-160 (2012).

19. Coffin, A. B. et al. Use of the swim bladder and lateral line in near-field sound source localization by fish. *J. Exp. Biol.* **217**, 2078-2088 (2014).

20. Olsen, K. *Directional Response in Herring to Sound and Noise Stimuli* (Int. Counc. Explor. Sea, 1969).

21. Blaxter, J. H. S., Gray, J. A. B. & Denton, E. J. Sound and startle responses in herring shoals. *J. Mar. Biol. Assoc. U. K.* **61**, 851-869 (1981).

22. Domenici, P. & Batty, R. S. Escape manoeuvres of schooling *Clupea harengus* in *J. Fish Biol.* **45**, 97-110 (1994).

23. Domenici, P. & Batty, R. S. Escape behaviour of solitary herring (*Clupea harengus*) and comparisons with schooling individuals. *Mar. Biol.* **128**, 29-38 (1997).

24. Moulton, J. M. & Dixon, R. H. in *Marine Bio-acoustics* Vol. 2 (ed. Tavolga, W. N.) 187-228 (Pergamon, 1967).

25. Eaton, R. C., Lavender, W. A. & Wieland, C. M. Identification of Mauthner-initiated response patterns in goldfish: evidence from simultaneous cinematography and electrophysiology. *J. Comp. Physiol.* **A** **144**, 521-531 (1981).

26. Eaton, R. C. & Emberley, D. S. How stimulus direction determines the trajectory of the Mauthner-initiated escape response in a teleost fish. *J. Exp. Biol.* **161**, 469-487 (1991).

27. Canfield, J. G. & Rose, G. J. Hierarchical sensory guidance of Mauthner-mediated escape responses in goldfish (*Carassius auratus*) and cichlids (*Haplochromis burtoni*). *Brain Behav. Evol.* **48**, 137-156 (1996).

28. Lewis, T. N. & Rogers, P. H. Directional acoustic startle response in the goldfish. *J. Acoust. Soc. Am.* **103**, 2863 (1998).

29. Preuss, T. & Faber, D. S. Central cellular mechanisms underlying temperature-dependent changes in the goldfish startle-escape behavior. *J. Neurosci.* **23**, 5617-5626 (2003).

30. Mirjany, M., Preuss, T. & Faber, D. S. Role of the lateral line mechanosensory system in directionality of goldfish auditory evoked escape response. *J. Exp. Biol.* **214**, 3358-3367 (2011).

31. Popper, A. N. & Coombs, S. Auditory mechanisms in teleost fishes: significant variations in both hearing capabilities and auditory structures are found among species of bony fishes. *Am. Sci.* **68**, 429-440 (1980).

32. Buwalda, R. J. A. & Van Der Steen, J. The sensitivity of the cod sacculus to directional and non-directional sound stimuli. *Comp. Biochem. Physiol.* **A** **64**, 467-471 (1979).

33. Nakatani, M., Miya, M., Mabuchi, K., Saitoh, K. & Nishida, M. Evolutionary history of Otophysi (Teleostei), a major clade of the modern freshwater fishes: Pangaean origin and Mesozoic radiation. *BMC Evol. Biol.* **11**, 177 (2011).

34. Fay, R. R. & Popper, A. N. Acoustic stimulation of the ear of the goldfish (*Carassius auratus*). *J. Exp. Biol.* **61**, 243-260 (1974).

35. Rogers, P. H., Popper, A. N., Hastings, M. C. & Saidel, W. M. Processing of acoustic signals in the auditory system of bony fish. *J. Acoust. Soc. Am.* **83**, 338-349 (1988).

36. Schulz-Mirbach, T. et al. In-situ visualization of sound-induced otolith motion using hard X-ray phase contrast imaging. *Sci. Rep.* **8**, 3121 (2018).

37. Vetter, B. J. & Sisneros, J. A. Swim bladder enhances lagena sensitivity to sound pressure and higher frequencies in female plainfin midshipman (*Porichthys notatus*). *J. Exp. Biol.* **223**, jeb225177 (2020).

38. Schulz-Mirbach, T. et al. Auditory chain reaction: effects of sound pressure and particle motion on auditory structures in fishes. *PLoS ONE* **15**, e0230578 (2020).

39. Maiditsch, I. P., Ladich, F., Heß, M., Schlepütz, C. M. & Schulz-Mirbach, T. Revealing sound-induced motion patterns in fish hearing structures in 4D: a standing wave tube-like setup designed for high-resolution time-resolved tomography. *J. Exp. Biol.* **225**, jeb243614 (2022).

40. von Frisch, K. & Dijkgraaf, S. Können Fische die Schallrichtung wahrnehmen? Z. Vgl. Physiol.* **22**, 641-655 (1935).

41. Gray, M., Rogers, P. H. & Zeddies, D. G. Acoustic particle motion measurement for bioacousticians: principles and pitfalls. *Proc. Mtgs. Acoust.* **27**, 010022 (2016).

42. van den Berg, A. V. & Schuijf, A. Acoustics of a standing wave tank for studying the hearing capacity of fish. *J. Acoust. Soc. Am.* **78**, 12-16 (1985).

43. Akamatsu, T., Okumura, T., Novarini, N. & Yan, H. Y. Empirical refinements applicable to the recording of fish sounds in small tanks. *J. Acoust. Soc. Am.* **112**, 3073-3082 (2002).

44. Duncan, A. J., Lucke, K., Erbe, C. & McCauley, R. D. Issues associated with sound exposure experiments in tanks. *Proc. Mtgs. Acoust.* **27**, 070008 (2016).

45. Mosharo, K. K. & Lobel, P. S. A comparison of underwater speakers for fish playback studies. *J. Acoust. Soc. Am.* **154**, 2365-2382 (2023).

46. Nedelec, S. L., Campbell, J., Radford, A. N., Simpson, S. D. & Merchant, N. D. Particle motion: the missing link in underwater acoustic ecology. *Methods Ecol. Evol.* **7**, 836-842 (2016).

47. Popper, A. N. & Hawkins, A. D. The importance of particle motion to fishes and invertebrates. *J. Acoust. Soc. Am.* **143**, 470-488 (2018).

48. Sand, O. in *Hearing and Sound Communication in Fishes* (eds Tavolga, W. N. et al.) 459-480 (Springer, 1981).

49. Mogdans, J. & Blalcmann, H. Peripheral lateral line responses to amplitude-modulated sinusoidal wave stimuli. *J. Comp. Physiol.* **A** **185**, 173-180 (1999).

50. Weeg, M. S. & Bass, A. H. Frequency response properties of lateral line superficial neuromasts in a vocal fish, with evidence for acoustic sensitivity. *J. Neurophysiol.* **88**, 1252-1262 (2002).

51. Higgs, D. M. & Radford, C. A. The contribution of the lateral line to 'hearing' in fish. *J. Exp. Biol.* https://doi.org/10.1242/jeb.078816 (2012).

52. Penalva, A. et al. Establishment of the miniature fish species *Danionella translucida* as a genetically and optically tractable neuroscience model. Preprint at bioRxiv https://doi.org/10.1101/444026 (2018).

53. Rajan, G. et al. Evolutionary divergence of locomotion in two related vertebrate species. *Cell Rep.* **38**, 110585 (2022).

54. Akbari, N. et al. Whole-brain optical access in a small adult vertebrate with two- and three-photon microscopy. *iScience* **25**, 105191 (2022).

55. Lam, P-Y. Longitudinal *in vivo* imaging of adult *Danionella cerebrum* using standard confocal microscopy. *Dis. Model. Mech.* **15**, dmm409753 (2022).

56. Lee, T. J. & Briggman, K. L. Visually guided and context-dependent spatial navigation in the translucent fish *Danionella cerebrum*. *Curr. Biol.* https://doi.org/10.1016/j.cub.2023.11.030 (2023).

57. Zada, D. et al. Development of neural circuits for social motion perception in schooling fish. Preprint at bioRxiv https://doi.org/10.1101/2023.10.25.563839 (2023).

58. Hoffmann, M., Henninger, J., Veith, J., Richter, L. & Judkewitz, B. Blazed oblique plane microscopy reveals scale-invariant inference of brain-wide population activity. *Nat. Commun.* **14**, 8019 (2023).

59. Cook, V. A. N. O. et al. Ultrafast sound production mechanism in one of the smallest vertebrates. *Proc. Natl Acad. Sci. USA* **121**, e2314017121 (2024).

60. Chan, V. & Perlas, A. in *Atlas of Ultrasound-Guided Procedures in Interventional Pain Management* (ed. Narouze, S. N.) 13-19 (Springer, 2011).

61. Schulz-Mirbach, T., Ladich, F., Plath, M. & Heß, M. Enigmatic ear stones: what we know about the functional role and evolution of fish otoliths: the role of fish otoliths in inner ear function. *Biol. Rev.* **94**, 457-482 (2019).

62. Colleye, O., Vetter, B. J., Mohr, R. A., Seeley, L. H. & Sisneros, J. A. Sexually dimorphic swim bladder extensions enhance the auditory sensitivity of female plainfin midshipman fish, *Porichthys notatus*. *J. Exp. Biol.* **222**, jeb204552 (2019).

63. Fay, R. R. The goldfish ear codes the axis of acoustic particle motion in three dimensions. *Science* **225**, 951-954 (1984).

64. Edds-Walton, P. L. in *Fish Hearing and Bioacoustics* Vol. 87 (ed. Sisneros, J. A.) 197-226 (Springer, 2016).

65. Lu, Z., Xu, Z. & Buchser, W. J. Acoustic response properties of lagena nerve fibers in the sleeper goby, *Dormitator latifrons*. *J. Comp. Physiol.* **A** **189**, 889-905 (2003).

66. Schuijf, A. in *Hearing and Sound Communication in Fishes* (eds Tavolga, W. N. et al.) 267-310 (Springer, 1981).

67. Fay, R. Soundscapes and the sense of hearing of fishes. *Integr. Zool.* **4**, 26-32 (2009).

68. Lara, R. A. & Vasconcelos, R. O. Characterization of the natural soundscape of zebrafish and comparison with the captive noise conditions. *Zebrafish* **16**, 152-164 (2019).

69. Schuijf, A. & Hawkins, A. D. Acoustic distance discrimination by the cod. *Nature* **302**, 143-144 (1983).

70. Bánki, O. et al. Catalogue of life checklist. *ChecklistBank*, https://doi.org/10.48580/df7lv (2023).

71. Schulz-Mirbach, T., Metscher, B. & Ladich, F. Relationship between swim bladder morphology and hearing abilities-a case study on Asian and African cichlids. *PLoS ONE* **7**, e42292 (2012).

72. Schulz-Mirbach, T., Heß, M., Metscher, B. D. & Ladich, F. A unique swim bladder-inner ear connection in a teleost fish revealed by a combined high-resolution microtomographic and three-dimensional histological study. *BMC Biol.* **11**, 75 (2013).

73. Ladich, F. & Schulz-Mirbach, T. Diversity in fish auditory systems: one of the riddles of sensory biology. *Front. Ecol. Evol.* **4**, https://doi.org/10.3389/fevo.2016.000028 (2016).

74. Finneran, J. J. & Hastings, M. C. A mathematical analysis of the peripheral auditory system mechanics in the goldfish (*Carassius auratus*). *J. Acoust. Soc. Am.* **108**, 1308-1321 (2000).

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

<figure><img src="image_13.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
### Animals

All animal experiments conformed to Berlin state, German federal and European Union animal welfare regulations and were approved by the LAGeSo, the Berlin authority for animal experiments. *D. cerebrum* were kept in commercial zebrafish aquaria (Tecniplast) with the following water parameters: pH 7.3, conductivity 350 µS cm⁻¹, temperature 27 °C. We used male and female adult fish between 4 and 11 months of age.

### Behavioural setup and protocol

The experimental setup comprised an inner 10 cm × 10 cm (length × width) tank with <200 µm thin optically opaque but acoustically transparent polypropylene sheet walls (cut out of plastic folders), surrounded by an outer tank with submerged speakers (4 × 3 Ekulit LSF-27M/SC 8 Ω in custom waterproof enclosures). Thus, the speakers were visually shielded from the fish inside the inner tank, and the fish were confined between the speakers (Fig. 1c provides a schematic; further details are provided in Extended Data Fig. 1a,b). The height of the water was 10 cm, and the transparent bottom of the inner tank was at 6.3 cm, leaving 3.7 cm to the water surface as a water column for the fish to swim in. The speakers were level with this water column, and all sounds were targeted for this water column. Infrared light-emitting diodes illuminated the fish from below. The inner tank was filmed with an overhead camera at 120 fps at 336 × 336 pixel resolution, and live tracking of the fish was carried out on a subset of frames at 15 fps. White light-emitting diodes lit the setup indirectly via reflections from the room walls. The room and water temperature was kept at 27 °C.

Each fish was tested once, and one fish was tested at a time. In the first minutes of the recording, a 10 cm × 10 cm acrylic plate with centimetre markings was placed in the inner tank to match the sound calibration grid to the video frame. Three minutes after placing the fish in the inner tank, playbacks were triggered for 45 min.

To probe into left-right directional hearing, playbacks from the front or back of the fish should be avoided. We prompted *D. cerebrum* to align with respect to the left-right speaker x axis to increase experimental throughput: previously, we observed that *D. cerebrum* swim closer to white than to black walls. By using two black plastic films as walls across the x axis and two white films across the y axis, we encouraged *D. cerebrum* to oscillate between the white walls, along the y axis (Extended Data Fig. 1c). Consequently, the ratio of distance covered along the y axis to the distance covered along the x axis was 1.6. Sound playback occurred only when fish were orthogonally oriented within a 45° angle measured from the orthogonal axis and within a 1.5 cm × 3 cm trigger zone at the centre, leaving at least 3.5 cm distance to the nearest wall (the typical startle displacement is mean ± s.d. = 1 cm ± 0.4 cm after the first 50 ms). The minimal delay between playbacks was set to ≥5 s with a minimum delay of 5 s plus a random delay, drawn from an exponential distribution with a mean of 5 s for each trial. This paradigm averaged to about 1.3 playbacks per minute in untreated fish and to about 0.6 playbacks per minute in lateral line-ablated fish.

Twelve target sounds were generated from a recorded pressure wave-form (see the section of the Methods entitled Sound stimulation wave-forms), targeted to the fish's current position to cancel reverberations (see the section of the Methods entitled Calibration and reverberation cancellation), and presented to the fish in random order following trigger events using custom-written code in Python 3.

The data in Figs. 1–4 and Extended Data Figs. 4–7 and 12a stem from 65 untreated fish (3,798 playbacks, 1,415 startles, about 37% startles). For each stimulus, we indicated the number of fish that responded with at least one startle. The same experiment, also comprising 12 sound configurations, was repeated with 74 lateral line-ablated fish (Extended Data Figs. 7 and 9; 2013 playbacks, 910 startles, 45% startles). A third sound playback experiment was carried out in the dark in 43 untreated fish, testing a subset of 4 sound configurations (Extended Data Fig. 12b).

### Behavioural analysis

**Tracking.** Pose tracking of *D. cerebrum*'s swimming behaviour was carried out with SLEAP⁷⁵. In total, 140 frames across nine random recordings of male and female fish were hand-labelled with a skeleton consisting of 7 equidistant nodes along the fish's body segments and 2 additional nodes, 1 for each eye. The 'single-animal' model was used for training. The model parameters and the trained model are available at the G-Node repository (see Data availability).

**Startle detection.** Plotting the fish's velocity against time around playback revealed a sharp increase in velocity after playback, clearly visible across all playbacks (Fig. 1e and Extended Data Fig. 4a). We defined a 25-ms time window around the time of peak velocity at which the speed distribution is bimodal and computed the average velocity in this time window for each trial to classify all playback trials with an average velocity above 17 cm s⁻¹ as startle trials (Extended Data Fig. 4b,c). The remaining ones were classified as non-startle trials. The decision criterion based on speed also resulted in a clear separation in terms of body bend (Extended Data Fig. 4d).

**Directional bias.** To classify startles into left or right, we measured the fish's x displacement during the first 50 ms after startle initiation (Fig. 1e). This duration was chosen because displacement heat maps at varying delays revealed that the initial, lateral displacement phase of the startle response peaks after 50 ms (Extended Data Fig. 4f). Wherever centred trajectories are shown, these initial 50 ms are depicted. The directional bias of the startle response is the fraction of startles to one indicated direction (left or right, away or towards speaker). This bias can be computed in two ways.

*Directional bias across trials and fish.* For each stimulus or set of stimuli, startle trials were pooled across all fish, and the fraction of startles in one direction was calculated. Using the two-sided binomial test, we calculated how likely a measured directional bias (approach or escape) would have been observed if the response was unbiased.

*Directional bias per fish.* In the analysis of bias across trials and fish, theoretically, all trials could stem from one performing animal (not of concern here; Extended Data Fig. 5a,b). To complement this measure, we also quantified the directional bias per fish. We had 12 sound configurations in each experiment and startles averaged to a total of about 22 startles per fish in an experiment; hence, a meaningful per-fish bias could be computed only on pooled sound configurations and for fish with many startles. To estimate the directional bias of individual fish, we filtered for fish that had ≥10 startles in both the single-speaker condition (pooled over 4 stimuli) and the trick condition (pooled over 4 stimuli; Extended Data Fig. 6c,d). Although the value reflects directional behaviour in the population and estimates fish-to-fish variability, it selects for fish that trigger many playbacks and startle often.

### Micro-CT

A 12-month-old male wild-type *D. cerebrum* was euthanized by ice shock and fixed with 4% paraformaldehyde in phosphate-buffered saline (PBS) at 4 °C overnight. The next day, the fish was washed for 15 min in PBS before being stained with 5% phosphomolybdic acid (Sigma Aldrich) solution in PBS at 4 °C overnight. After staining, the fish was washed in PBS for 15 min before embedding in 1% PBS-buffered agarose inside a cryo tube. The micro-CT scan was carried out at the ANATOMIX beamline at SOLEIL synchrotron by XPLORAYTION. The sample was placed into a 40-keV polychromatic (white) X-ray beam. A scan consisted of 3,200 projections collected at about ×10 optical magnification by a digital camera (Orca Flash 4.0 V2) with a sensor pixel size of 6.5 µm at 150 ms exposure time, yielding an effective pixel size of 0.6485 µm. The registered data were binned to 1.2970 µm voxel size. Key structures of the hearing apparatus were manually segmented. To this end, planes were hand-labelled using 3D Slicer⁷⁶ (v5.6, https://
slicer.org) and then interpolated using Biomedisa (v23)⁷⁷. FIJI ImageJ (v1.5)⁷⁸ was used to convert between different file types. The segments were turned into mesh grids and loaded into Blender for cleaning and rendering.

## Lateral line ablation and DASPEI staining

To rule out that the lateral line organ senses sound directionality in our experiments, we ablated the lateral line using neomycin⁷⁹. To ablate the neuromasts, fish were placed in a 200 µM neomycin solution for about 30 min. Afterwards, they were transferred to a beaker with tank water. Behavioural experiments started after ≥30 min. To confirm the reliability of the lateral line ablation protocol, we stained 30 neomycin-treated fish. After the behaviour experiment, they were transferred to a 100 µM DASPEI (2-[4-(dimethylamino)styryl]-1-ethylpyridinium iodide) solution and then to a beaker with tank water to wash out unbound DASPEI. Afterwards, the fish were euthanized with an ice shock and imaged with an epifluorescence microscope. Neuromasts were reliably stained in control fish but not in neomycin-treated fish, indicating reliable ablation (see Extended Data Fig. 9e,f for example images). As functional metrics we report an increase in number of wall contacts after startles (Extended Data Fig. 9g) and a decrease in foraging strikes in the dark (Extended Data Fig. 9h) in neomycin-treated fish.

## Vibrometry

**Confocal microscope.** The confocal reflectance microscope was based on a custom-built laser-scanning two-photon microscope (Extended Data Fig. 10a). The illumination source was a Ti:sapphire laser (MaiTai DeepSee; SpectraPhysics) operated at 810 nm (with or without mode-locking). Before entering a laser-scanning two-photon microscope, the beam passed through a 90:10 beam splitter (90% reflection, 10% transmission). The light back-scattered by the fish inner structures was descanned, reflected by the 90:10 beam splitter, and then focused by a lens (f= 50 mm) into a single-mode fibre (core diameter: 25 µm, numerical aperture: 0.1) acting as a confocal pinhole. The microscope was controlled by custom-written software (https://github.com/dationella/Ismaq).

**Acoustic stimulation.** Fish were anaesthetized in 120 mg l⁻¹ fish water-buffered MS-222. They were subsequently placed on a pre-formed agarose mould, which allowed the gill covers to move freely, and immobilized with 2% low-melting-point agarose (melting point 25 °C). A flow of aerated aquarium water (with anaesthetic) was delivered to their mouth through a glass capillary.

The fish was acoustically stimulated using two facing speakers sealed in custom-made waterproof enclosures. The diaphragms were exposed to water. The speakers were each placed about 1.3 cm away from the fish. They were driven using a DAQ card (National Instruments USB-6211), connected through audio amplifiers (Kemo M031N, 3.5 W). Pressures of up to about 176 dB (referenced to 1 µPa) were thus generated at the fish position in the pressure-only configuration and particle motion of up to about 8 mm s⁻¹ in the particle-motion-only configuration, consistent with the expected amplitude relationship between pressure and particle motion in the sound monopole near field.

**Motion phase maps.** The principle of the laser-scanning vibrometric measurement is illustrated in Fig. 3b and Extended Data Fig. 10b. The sample (Extended Data Fig. 10b(i)) was stimulated with an acoustic sinusoidal wave at frequency f<sub>stim</sub>, and imaged with a laser-scanning microscope with a line rate f<sub>scan</sub> (Extended Data Fig. 10b(ii)).

To reconstruct amplitudes and relative phases of sinusoidal object motion, we needed to measure each pixel under more than two different phases according to the Shannon-Nyquist sampling theorem. As noise can influence this measurement, we used four phase steps here, ensuring proper phase reconstruction while keeping acquisition sessions reasonably short.

To reconstruct the displacement of the moving structures inside the fish, each line of the image was repeatedly scanned nStep = 4 times, with a phase offset of π/2 between each line (Extended Data Fig. 10b(iii)). To this end, the stimulation frequency and the line rate must follow the relationship:

$$f_{\text{stim}} = (N + 1/\text{nStep})f_{\text{scan}}$$

with N being an integer. To maximize the line rate, we took N = floor(f<sub>stim</sub>/f<sub>scan</sub>).

This in turn set additional constraints on the various scanning parameters. We used f<sub>scan</sub> = 800 Hz and f<sub>stim</sub> = 1,000 Hz for the data presented in Fig. 3 and Extended Data Fig. 11.

To ensure repeatable measurements, the acoustic stimulation and the galvanometric scanning mirrors were synchronized so that each pixel was recorded at a known sound phase. This was achieved by triggering the sound generation on each single frame scan trigger.

Doing so, each pixel was stroboscopically probed at nStep = 4 different phases of the acoustic stimulation cycle. As sound propagates while scanning two consecutive pixels, the probed acoustic phase is shifted by 2π × pixelPeriod, which was taken into account in the motion reconstruction of the imaged structures (Extended Data Fig. 10c). These images were then reshaped to yield an (Nx,Ny,nStep) dataset (Extended Data Fig. 10b(iv)).

To analyse the motion of the inner structures of the fish, we used Matlab 2019b and a particle image velocimetry toolbox PIVlab⁸⁰, originally developed to characterize the motion of flowing particles for fluid mechanics. Essentially, the particle displacement is assessed by cross-correlating subregions with decreasing sizes of consecutive images (Extended Data Fig. 10b(v)). The contrast of the reflectance images was enhanced before the displacement analysis, and the results were curated in post-processing by removing outliers and interpolating detection gaps.

The motion detection yielded x- and y-displacement maps at each of the four phases in the acoustic stimulation period. The first Fourier component was computed for each pixel to extract the amplitude and phase of the local displacement (Extended Data Fig. 10b(vi)). The phase was finally corrected for the accumulating phase offset along the horizontal x direction due to the line scanning procedure (Extended Data Fig. 10c). Owing to the synchronization of the acoustic stimulation with the line scanning process, we could carry out this measurement in several planes and obtain a consistent volumetric complex map characterizing the motion response of the various inner structures to the acoustic stimulation. Maximum-amplitude projections across planes delivered the shown two-dimensional phase maps, one for motion along the speaker-speaker axis (x) and one for motion orthogonal to the speaker-speaker axis (y).

## Sound stimulation waveforms

We reasoned that *D. cerebrum* sense pressure and particle motion. Hence, our sound stimuli were defined in terms of three quantities: pressure, x acceleration and y acceleration, which were delivered to the fish’s current position by utilizing the frequency responses of speakers to cancel position-dependent reverberations (see the section of the Methods entitled Calibration and reverberation cancellation). y acceleration was always kept at zero, and only pressure and x acceleration were varied. In summary, 12 sounds were generated from a recorded pressure waveform and presented to the fish in a random sequence upon trigger events. The 12 sounds consisted of four single-speaker sounds (left or right × positive polarity or negative polarity), two sounds with only a pressure component (positive polarity or negative polarity), two sounds with only horizontal x-motion components (positive polarity or negative polarity) and four trick conditions, which exactly matched the four single-speaker target waveforms, but differed by the speakers that were active to realize these.
We observed that *D. cerebrum* startle when we drop a cylindrical piece of rubber into the water. We recorded the pressure waveform of this sound, high-pass filtered it at 100 Hz, and extracted a 12-ms snippet to serve as our pressure waveform template (note that conditioned sounds—that is, the actual speaker signals—were band-pass-filtered between 200 Hz and 1,200 Hz; see the following section). The target pressure amplitude was set to a peak sound pressure level of 167 dB (referenced to 1 µPa) by rescaling this pressure waveform accordingly. This amplitude was loud enough to elicit startles reliably and still supported by our small 2.7-cm-diameter speakers. The first peak's rise time (10% to 90% absolute amplitude) was 0.664 ms and the centre frequency of the pulse was about 780 Hz. The target horizontal particle acceleration waveform was computed from the pressure waveform using monopole theory for each Fourier component, as follows.

The pressure signal decays as $1/r$ with radial distance $r$ away from a sound monopole with amplitude $p_0$ at distance $r_0$

$$ \hat{p}(r, t) = \hat{p}_0 \frac{r_0 e^{ikr}}{r} e^{-i\omega t} $$

and with frequency $f$, $\omega = 2\pi f$, wavenumber $k = 2\pi/\lambda$, wavelength $\lambda$ and speed of sound $c = \lambda f$.

In a medium of density $\rho$, the radial particle velocity decays quadratically with distance in the near field ($kr \ll 1$, limit dependent on frequency):

$$ \hat{v}_r(r, t) = \left[ \frac{1}{\rho c} \left( 1 + \frac{i}{kr} \right) \hat{p}(r, t) \right] $$

By contrast, particle acceleration—the temporal derivative of particle velocity—decays quadratically with distance for nearby sounds ($r \ll 1$, limit independent of frequency):

$$ \hat{a}_r(r, t) = -i\omega\hat{v}_r(r, t) = \frac{1}{\rho}\left(\frac{1}{r} - ik\right)\hat{p}(r, t) $$

To compute the particle acceleration $a_r(r, t)$ at a distance $r$ to a sound monopole with pressure $p(r, t)$ for discrete signals of arbitrary waveform, we applied this equation separately for each Fourier component. Given a pressure waveform $\{\mathbf{p}_n\} := p_0 p_{1r} \cdots p_{N-1}$ with $N$ samples $p_{n'}$ spaced at $T=1/sr$ with sample rate $sr$, the particle acceleration $\{\mathbf{a}_n\} := a_0 a_{1r} \cdots a_{N-1}$ that would be observed at a distance $r = r_0$ from a sound monopole was calculated by carrying out the discrete Fourier transform $\{\mathbf{P}_l\} := P_0 P_{1r} \cdots P_{N-1}$

$$ P_l = \sum_{n=0}^{N-1} p_n e^{-i2\pi l \ln n} $$

and deriving particle acceleration for each Fourier component $\{\mathbf{A}_n\} := A_0 A_{1r} \cdots A_{N-1}$ independently. With corresponding frequencies $f_l \approx l/(NT)$, such that $k \approx 2\pi l/(NT)$, and the relationship between pressure and particle acceleration, $A_l$, is calculated as

$$ A_l = \frac{1}{\rho} \left( \frac{1}{r_0} - ik_l \right) P_l $$

which defines the radial particle acceleration through the inverse Fourier transform:

$$ a_n = \frac{1}{N} \sum_{l=0}^{N-1} A_l e^{i \frac{2\pi}{N} l n} $$

In the experiments, $r_0$ was set to 3 cm, thus simulating a monopole sound source at 3 cm, irrespective of *D. cerebrum*'s relative position to the speakers. This resulted in a peak particle acceleration of 7.59 m s⁻².

Other parameters used were $c = 1,500$ m s⁻¹, $\rho = 1,000$ kg m⁻³ and $sr = 51,200$ Hz. In terms of pressure, $x$ acceleration and $y$ acceleration ($p, a_x$ and $a_y$), there were eight different target configurations, with '+' indicating polarity of the template waveform and '-' indicating opposite polarity: four monopole configurations $(+,+,0), (-,-,0), (+,-,0)$ and $(-,+,0)$; two pressure configurations $(+,0,0)$ and $(-,0,0)$; and two motion configurations $(0,+0)$ and $(0,-0)$. Despite a total of eight target configurations, there were 12 sound configurations as the four monopole configurations can be realized in two ways, either with a single speaker or with three speakers (trick configuration; see the next section).

## Calibration and reverberation cancellation

Conducting experiments in small tanks presents challenges as both tank geometry and the receiver's position affect the sound amplitude and waveform sensed by the receiver (Extended Data Fig. 2c). By recording the speakers' impulse responses inside the inner tank in terms of pressure and particle acceleration (Extended Data Fig. 2b), speakers could be activated to precisely control pressure and particle acceleration components at the fish's location (Extended Data Fig. 2d).

### Pressure and acceleration measurements. Pressure.
Pressure was measured with a hydrophone (Aquarian Scientific AS-1, preamplifier: Aquarian Scientific PA-4, acquisition: NI-9231 sound and vibration module, National Instruments; Extended Data Fig. 1d). During repeated playback of the same sound, a single hydrophone was automatically moved across a 5 × 5 grid inside the inner tank, sampling with a spacing of 1.5 cm (Extended Data Fig. 1c,e). Hence, a 25-point pressure field was obtained for each sound configuration, spanning a 6 cm × 6 cm square at the centre of the inner tank between the speakers.

### Acceleration.
Particle acceleration was measured in two ways.

In the first method, particle acceleration was measured indirectly through the pressure gradient. Newton's second law of motion (pressure gradient force)

$$ \mathbf{a} = -\frac{1}{\rho} \nabla P $$

links the spatial pressure gradient to particle acceleration. In water, with density $\rho = 1,000$ kg m⁻³ and speed of sound $c = 1,500$ m s⁻¹, the following approximation holds for pressure signal frequencies $f \ll 100$ kHz, if the pressure gradient is sampled with step size $x_2 - x_1 = 1.5$ cm:

$$ a_x = - \frac{1}{\rho} \frac{p(x_2) - p(x_1)}{x_2 - x_1} $$

The approximation holds for all frequencies used in this experiment. For measuring gradients, moving a single hydrophone is preferred over a hydrophone array, as the gradient could be biased by small differences in hydrophone sensitivity and perturbations of the sound field by the presence of other hydrophones. We calculated $x$ and $y$ acceleration on the basis of the 25-point pressure field recorded with a single hydrophone. The pressure field included points outside the trigger zone to compute pressure gradients (that is, acceleration) across the trigger zone boundary.

In the second method, particle acceleration was additionally directly measured along all three axes with an acceleration sensor (Triaxial ICP - Model 356A45, PCB Piezotronics, acquired with NI-9231 sound and vibration module, National Instruments; Extended Data Fig. 1d). Like the hydrophone, the acceleration sensor was moved across all 5 × 5 grid positions during repeated playback of the same sound, giving measurements for $x, y$ and $z$ acceleration.

Whereas hydrophones are manufactured and calibrated for underwater use, the particle acceleration sensor is not made to measure particle acceleration underwater and is meant to be glued onto the
vibrating object. Owing to an acoustic impedance mismatch between metal and water, we expected the PCB sensor to underestimate particle acceleration.

We compared x and y acceleration waveforms for both measurement methods and found that the acceleration waveforms acquired through the direct method match the waveforms acquired through the indirect method after multiplication by a factor of about 2.4. The close match in rescaled waveforms confirms the validity of the gradient approximation in the indirect method.

Hence, in all experiments, x and y acceleration were measured through the indirect method, on the basis of spatial pressure gradients. The particle acceleration sensor still proved useful in measuring the vertical z acceleration in our setup.

**Impulse response-based sound targeting.** To create the same sounds at any position inside the inner tank, impulse responses for all 4 speakers were measured across 25 positions on a 5 × 5 grid with 1.5-cm spacing. In the following, the sound targeting method is described for one position.

Let $k_{i,p}$ be the pressure impulse response kernel, $k_{i,a_x}$ be the x acceleration impulse response kernel, and $k_{i,a_y}$ be the y acceleration impulse response kernel for the ith speaker. Using M speakers with signal $s_i$, pressure and acceleration can be predicted through convolution ($*$):

$$ 
\begin{aligned}
p &= \sum_{i=0}^{M-1} k_{i,p} * s_i \\
a_x &= \sum_{i=0}^{M-1} k_{i,a_x} * s_i \\
a_y &= \sum_{i=0}^{M-1} k_{i,a_y} * s_i
\end{aligned}
$$

In the Fourier domain, utilizing the convolution theorem, these become a system of equations for each Fourier component $l$.

$$ 
\begin{aligned}
P_l &= \sum_{i=0}^{M-1} K_{i,p,l} S_{i,l} \\
A_{x,l} &= \sum_{i=0}^{M-1} K_{i,a_x,l} S_{i,l} \\
A_{y,l} &= \sum_{i=0}^{M-1} K_{i,a_y,l} S_{i,l}
\end{aligned}
$$

On the basis of the Fourier components of the target waveforms (see the section of the Methods entitled Sound stimulation waveforms), $P_l$, $A_{x,l}$ and $A_{y,l}$, and the Fourier components of the impulse response kernel $K_{i,p,l}$, $K_{i,a_x,l}$ and $K_{i,a_y,l}$, the system of equations can be solved for the Fourier components of the speaker signals $S_{i,l}$ as long as $M \ge 3$ and the kernel components are non-zero and non-identical. The time-domain signal for the ith speaker is then given by the inverse Fourier transform using components $S_{i,l}$.

To increase robustness of the solutions (for example, to avoid speakers cancelling themselves unnecessarily and to limit speaker amplitude), speaker signal waveforms were forced to become similar to the target waveform. This was implemented by solving the system of equations with a least-square solver (scipyImprove.lsq linear) with bounds $-B_{i,l} < S_{i,l} < B_{i,l}$. The bound $B_{i,l}$ was computed as a rescaling of the absolute Fourier components of the target pressure waveform $P_l$

$$ B_l = \alpha_l y | P_l | $$

in which y is fixed and scales pressure to voltage and $\alpha_l$ is a rescaling parameter set independently for each speaker to give additional control over active speakers. We list our values for $\alpha_l$ used in different sound configurations in Supplementary Table 2.

After conditioning, all computed speaker signals were band-pass-filtered between 200 Hz and 1,200 Hz to avoid activating the lateral line.

To ensure that the trick configuration differed from the single-speaker configuration only by selective pressure inversion, a two-step sound conditioning was carried out. First, the speaker signals for the single-speaker configuration were calculated. Then, these signals were effectively fixed to closely resemble the single-speaker signal and only activations of the two speakers along the orthogonal axis were conditioned.

The above calculation was carried out for the 25 grid positions. The computed speaker signals accurately delivered the target waveforms to the target position (Extended Data Figs. 2d and 3 and Supplementary Table 2). To ensure consistency over experiments, the water level was kept at precisely 10 cm, and the pressure and acceleration fields inside the inner tank were measured several times (this includes before the first recording and after the last recording).

During the experiment, the fish's x-y position was detected at 15 Hz, and the loading for the speakers was linearly interpolated on the basis of targeted sounds at neighbouring grid positions.

In the section entitled Sound stimulation waveforms, we describe how we defined the pressure and particle motion target waveforms that were conditioned this way.

## Estimation of binaural cues (P-ILD, M-ILD, P-ITD, M-ITD)

**ILDs.** To estimate binaural cues in our behavioural experiment, we analysed the pressure and particle motion at sound calibration grid points 3 cm apart, ($x_0$) 1.5 cm to the left and ($x_1$) 1.5 cm to the right of the centre grid point. To estimate sign and peak amplitude of level differences, we calculated P-ILD (pressure ILD) as $max_t(\text{abs}(p(x_0, t))) - max_t(\text{abs}(p(x_1, t)))$ and M-ILD (particle motion ILD) as $max_t(\text{abs}(a_x(x_0, t))) - max_t(\text{abs}(a_x(x_1, t)))$. The level differences between these two points were divided by a factor 50 to estimate the level difference across the left-to-right inner ear axis of the fish (about 0.6 mm). Comparing the single-speaker configuration with the trick configuration, these data show that the sign of M-ILD remains the same (+0.11 m s⁻² versus +0.30 m s⁻²), but the sign of P-ILD is inverted (+4.4 Pa versus -4.6 Pa). For a geometrical illustration of the inversion of P-ILD, see Extended Data Fig. 8d.

**ITDs.** ITDs were estimated by calculating the phase propagation in different sound configurations under a monopole approximation (Extended Data Fig. 8c-f).

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

Fish trajectory, micro-CT and vibrometry data have been deposited in the G-Node repository at https://gin.g-node.org/danionella/veith_et_al_2024. Downsampled segmentation mesh-grid data have been deposited to MorphoSource at https://n2t.net/ark:/87602/m4/623999. Source data are provided with this paper.

## Code availability

Code used for data analysis is available on GitHub at https://github.com/danionella/veith_et_al_2024.

<footnote n="75">75. Pereira, T. D. et al. SLEAP: a deep learning system for multi-animal pose tracking. *Nat. Methods* **19**, 486–495 (2022).</footnote>

<footnote n="76">76. Fedorov, A. et al. 3D Slicer as an image computing platform for the quantitative imaging network. *Magn. Reson. Imaging* **30**, 1323–1341 (2012).</footnote>
77. Lösel, P. D. et al. Introducing Biomedisa as an open-source online platform for biomedical image segmentation. *Nat. Commun.* **11**, 5577 (2020).

78. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. *Nat. Methods* **9**, 676–682 (2012).

79. Harris, J. A. et al. Neomycin-induced hair cell death and rapid regeneration in the lateral line of zebrafish (*Danio rerio*). *J. Assoc. Res. Otolaryngol.* **4**, 219–234 (2003).

80. Thielicke, W. & Stamhuis, E. J. PIVlab – towards user-friendly, affordable and accurate digital particle image velocimetry in MATLAB. *J. Open Res. Softw.* **2**(1): e30, (2014).

**Acknowledgements** We thank M. Brecht, M. Larkum, J. Poulet, M. Lovett-Barron, T. Veith, B. Hesse, D. Markov, M. Kadobianskyi, V. Cook, L. Constien, J. Henninger, A. Gröneberg, U. Böhm, L. Kosina, F. O'Donovan, L. Schulze and J. Tiefenbacher for discussions and for critically reading our manuscript; and our fish facility team, A. Wrana, M. Renz, D. Krap, H. Schmidt, N. Kroworz, for fish care and experimental support. We acknowledge support by the European Research Council (ERC2021-CoG-101043615), the Einstein Foundation (EPP-2017-413), the German

Research Foundation (DFG, project EXC-2049-390688087) and the Alfred Krupp von Bohlen und Halbach Foundation. T.C. received a fellowship from the Alexander von Humboldt Foundation.

**Author contributions** Conceptualization: J.V., T.C. and B.J. Methodology: J.V., T.C., M.H. and B.G. Investigation: J.V., T.C., A.S., L.E.D. and M.H. Visualization: J.V. and T.C. Funding acquisition: B.J. Project administration: B.J. Supervision: T.C. and B.J. Writing: J.V., T.C. and B.J.

**Competing interests** The authors declare no competing interests.

**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41586-024-07507-9.

**Correspondence and requests for materials** should be addressed to Benjamin Judkewitz.

**Peer review information** Nature thanks Catherine Carr, Florian Engert and Michael Fine for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints.
<figure><img src="image_14.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>e</figcaption></figure>

Pressure measurement,
indirectly: 2D acceleration

<figure><img src="image_18.png" /></figure>

3D acceleration
measurement

<figure><img src="image_19.png" /></figure>

**Extended Data Fig. 1 | Setup design and sound calibration equipment.**

a, Photo of the setup. b, Schematic of the setup indicating size and material. The room and water temperature is kept at 27 °C; the room light is indirect light reflecting from white walls. Infrared (IR) light illuminates the setup from below. c, Schematic of sound calibration points inside the inner tank overlaid over a heatmap of fish positions inside the inner tank across all 65 recordings. *D. cerebrum* avoids the black walls and oscillates between the two white walls. Pressure and particle acceleration were measured at 25 points (x) on a grid spanning a 6 cm x 6 cm square, covering the trigger zone (green dotted

rectangle). Sounds were conditioned to deliver target waveforms to the fish's swimming position within the trigger zone based on each speaker's pressure and acceleration impulse responses for these 25 points. Points outside the trigger zone were measured to calculate pressure gradients across the edge of the trigger zone, see Methods. d, Left: Hydrophone, Aquarian Scientific, AS-1. Right: Triaxial acceleration sensor, ICP®-Model 356A45, PCB Piezotronics. e, The motorized arm moved the sensor to each grid position. Left: pressure measurement (was also used to compute acceleration indirectly). Right: direct acceleration measurement with the PCB sensor.
a

<figure><img src="image_20.png" /></figure>

b

<figure><img src="image_21.png" /></figure>

c

Same speaker activation for each fish position

<figure><img src="image_22.png" /></figure>

Different pressure and acceleration for each fish position

<figure><img src="image_23.png" /></figure>

d

Different speaker activation for each fish position

<figure><img src="image_24.png" /></figure>

Same pressure and acceleration for each fish position

<figure><img src="image_25.png" /></figure>
**Extended Data Fig. 2** | Sound targeting cancels reverb effects. **a**, Schematic of sound calibration points (x) inside the inner tank. Green dotted rectangle: trigger zone. Orange dotted square: 9 measurements and target locations shown in **c-d. b**, Impulse response kernel for the left speaker at the center calibration point for three sound cues: pressure, x acceleration, and y acceleration. Pressure and x/y acceleration impulse responses were recorded for each speaker and each position on the grid, totaling to 3 x 4 x 25 = 300 waveforms capturing the acoustic properties of the tank. **c**, Naive playback paradigm (not used): The same waveform is played back via the left speaker, and sound cues are recorded at nine positions. **Left:** constant speaker signals.

**Right:** Pressure and x/y acceleration measured at nine positions, spanning a 4.5 cm x 4.5 cm square at the center of the inner tank, vary considerably across this area, and the measured waveform does not resemble the playback waveform. **d**, Sound conditioning paradigm (used): Based on the impulse responses, different speaker signals are calculated for each target position. **Left:** speaker signals are different for each target position. Top and bottom speakers help to cancel y acceleration. **Right:** Pressure and x/y acceleration measured at the same nine positions. The waveforms resemble the defined target waveforms and are stable across positions.
a

<figure><img src="image_26.png" /></figure>

Row 1: Left speaker
Speaker activations
(Volts)

Row 2-4: x-, y-
acceleration,
pressure

prediction
target
measured

b

<figure><img src="image_27.png" /><figcaption>ii) Single left speaker -</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>iii) Single right speaker -</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>iv) Pressure configuration +</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>v) Particle motion configuration (1)</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>vi) Trick configuration, left -</figcaption></figure>

<figure><img src="image_32.png" /></figure>

**Extended Data Fig. 3** | Accuracy of six exemplary sound stimuli. a, Schematic of sound calibration points (x) inside the inner tank, indicating the sound target position and recording position for speaker signals and recordings shown in b, b, Pressure and acceleration measurements for sound configurations i) - vi) (see Fig. 2c for schematic), targeted for the inner tank's center calibration point and measured at the same position. Target waveforms (target) are defined as pressure, x accelerations, and y accelerations. Peak

pressure: 223.87 Pa, peak x acceleration: 7.59 m/s². The speaker signals (top panels in i-vi) are calculated based on impulse responses to deliver the target waveforms to the target position with the constraint that, e.g., the right speaker is inactive in the left speaker configurations (i & ii) and in the left speaker trick configurations (vi). For a given speaker signal, convolution with the respective impulse response kernels predicts (prediction) the measured waveforms.
<figure><img src="image_33.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_34.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_35.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_36.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_37.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_38.png" /></figure>

**Extended Data Fig. 4** | Startle detection and dynamics. **a**, Fish translational speed, shown for all n = 3798 playbacks across twelve stimulus configurations, and across 65 fish aligned to playback trigger at t = 0 ms (thin lines). The mean speed across playbacks is shown as a thick line. To trigger a playback, the fish had to swim into the trigger zone. This explains the small increase in mean speed prior to t = 0 ms (white arrow). **b**, We used a 25 ms time window (between vertical dashed lines) and classified trials as startle trials if the temporal average speed across this time window exceeded a threshold of 17 cm/s (horizontal line). The remaining ones were classified as non-startle trials. The threshold was set such that the average speed across non-startle trials increased slightly after playback (black arrow). This reflects a conservative choice for startle detection, classifying a few startles as non-startles rather than classifying non-startles as startles. **c**, Top: The average speeds computed

across a 25 ms time window before each playback are below the selection threshold of 17 cm/s. Bottom: The average speed distribution within the 25 ms time window indicated in **b**, after playback, is bimodal. Hence, startles can be readily identified by a speed threshold. **d**, Body bend: The average absolute angle taken between six edges along the fish's body axis after playback across startles increases sharply. It remains constant if averaged across non-startles. Startles are initiated by a strong body bend. **e**, Centered trajectories for all playbacks with startle reaction (n = 1415) across all 65 fish after aligning initial fish heading. The first 50 ms of startle responses are shown. **f**, Normalized heatmaps of the fish head position at consecutive times relative to the position at the start of the startle response. The top row shows snapshots of the trajectories shown in **e**. After a fast initial displacement to the left or right, *D. cerebrum* continues on a forward trajectory.
<figure><img src="image_39.png" /><figcaption>c</figcaption></figure>

Minimal interplayback interval: 5 s + X, X ~ Exp(λ = 5 s)

<figure><img src="image_40.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_41.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_42.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_43.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_44.png" /></figure>
**Extended Data Fig. 5 | Habituation of the startle response.** a–b, Startle probability of each fish across all sound configurations as a function of the number of playbacks triggered by the fish. a, For main experiment (in untreated fish). b, For experiment in neomycin-treated fish. c–f, Details on habituation, data from main experiment (in untreated fish). c, Startle probability (the fraction of playbacks with startles, computed per bin) increases with time since the previous playback (green curve, the 95% confidence interval is given by Wilson scores). d, Startle probability (the fraction of playbacks with startles, computed per bin) decreases throughout the recording (green curve, the 95% confidence interval is given by Wilson scores). e, Startle probability as 2D binned statistics to show the interaction between time in recording and time since the last playback with the number of previous playbacks (Nth playback); white: no data. f, Average speeds over the

course of the recording. Because fish that startle after playback have a higher speed before playback (see also Extended Data Fig. 4b), a decreased swimming speed could explain the decreasing startle probability over the course of the recording shown in d. Average speed is computed across 25 ms at a time 1 s before sound onset (before playback) and across 25 ms after sound onset (after playback), separately for startle and non-startle trials. The confidence interval indicates the standard error of the mean. Note the split y-axis. The speed of fish before playback stayed constant over the course of the recording. Hence habituation is not explained by a decrease in swimming activity over 50 min. a–f, Together, startle probability increases if the fish is in a fast swimming state and decreases, the less time has passed since the last playback and the more playbacks have been played previously (both signs of habituation).
a

<figure><img src="image_45.png" /></figure>

b

<figure><img src="image_46.png" /></figure>

c

<figure><img src="image_47.png" /></figure>

d

<figure><img src="image_48.png" /></figure>
**Extended Data Fig. 6 | Directional startles are present for both sexes and in individual fish.** Centered startle trajectories for single speaker configurations and trick configurations, same data as shown in Fig. 2f but separated by sex. Both male and female *D. cerebrum* startle away from single active speakers and startle towards the active speaker in the trick configuration. **a**, Single speaker configuration. Left: females (75% of 186 startles away from the speaker in 26 female fish that startled, two-sided binomial test: $p = 9 \times 10^{-12}$). Right: males (83% of 312 startles away from the speaker in 38 male fish that startled, two-sided binomial test: $p = 1 \times 10^{-33}$). **b**, Trick configuration. Left: females (72% of 107 startles towards the speaker in 23 female fish that startled, two-sided binomial test: $p = 6 \times 10^{-6}$). Right: males (64% of 224 startles towards the speaker in 38 fish that startled, two-sided binomial test: $p = 2 \times 10^{-5}$). **a-b**, Percentages indicate the fraction of startle displacements into the direction of the blue arrow. The heatmaps are normalized and smoothed 2D histograms

over endpoint positions of the shown trajectories. **c**, To estimate the directional bias in individual fish, we selected fish with at least one startle in both the single speaker condition and the trick condition (N = 46) to compute the number of fish with a mean directional bias away from the speaker (bias > 0.5): 41 of 46 fish in the single speaker configuration, two-sided binomial test: $p = 4 \times 10^{-8}$, and 10 of 46 in the trick configuration, two-sided binomial test: $p = 2 \times 10^{-4}$. Hence, individual fish escaped away from the single left-right speaker but approached it in the trick condition. **d**, The mean directional bias computed across fish with $\ge 10$ startles in both conditions (N = 5, a subset of **a**) amounted to a directional bias of mean $\pm$ std. = $84\% \pm 13\%$ away from the speaker in the single speaker configuration and to $72\% \pm 9\%$ towards the speaker in the trick configuration. **c-d**, Line thickness is proportional to the number of minimum startles in response to one of the two sets of stimuli. Blue line: average directional bias across fish directional biases.
a

<figure><img src="image_49.png" /></figure>

Without lateral line ablation: 65 fish tested

<figure><img src="image_50.png" /></figure>

b

<figure><img src="image_51.png" /></figure>

With lateral line ablation: 74 fish tested

**Extended Data Fig. 7** | Data summary on all twelve stimuli for both cohorts.

Centered startle trajectories and displacement heatmaps for each stimulus used in the experiments, startle trials are pooled across all tested fish. **a**, for fish without an ablated lateral line. **b**, for fish with an ablated lateral line. **a-b,**

Percentages indicate fraction of startles to the right. P-values report the probability that startles are unbiased in any direction (Two-sided binomial test, null hypothesis: p = 0.5, Bonferroni-corrected: n = 12 stimuli).
a

<figure><img src="image_52.png" /><figcaption>Sound monopole: 1 kHz</figcaption></figure>

<figure><img src="image_53.png" /></figure>

b

<figure><img src="image_54.png" /><figcaption>Sound monopole: 1 kHz</figcaption></figure>

$$ \dot{\tilde{p}}(r,t) = \hat{p}_0 \frac{r_0 e^{ikr}}{r} e^{-i\omega t} $$

$$ \dot{\tilde{v}}_r(r,t) = \left[ \frac{1}{\rho c} \left( 1 + \frac{i}{kr} \right) \right] \dot{\tilde{p}}(r,t) $$

c

<figure><img src="image_55.png" /><figcaption>i) Single speaker +</figcaption></figure>

<figure><img src="image_56.png" /></figure>

<figure><img src="image_57.png" /></figure>

d

<figure><img src="image_58.png" /><figcaption>vi) Trick configuration</figcaption></figure>

<figure><img src="image_59.png" /></figure>

<figure><img src="image_60.png" /></figure>

e

<figure><img src="image_61.png" /><figcaption>iv) Pressure configuration</figcaption></figure>

<figure><img src="image_62.png" /></figure>

<figure><img src="image_63.png" /></figure>

f

<figure><img src="image_64.png" /><figcaption>v) Particle motion configuration</figcaption></figure>

<figure><img src="image_65.png" /></figure>

<figure><img src="image_66.png" /></figure>
**Extended Data Fig. 8** | Sound monopoles and sound configurations. **a,** Pressure level (dashed line) falloff next to a sound monopole at several phase snapshots of a propagating wave. Left: Falloff of a 1 kHz wave over 1.5 m. Right: The falloff across the width of the fish at 3 cm distance to a monopole sound source stems from the level falloff with distance. **b,** Both amplitude ratio and relative phase between pressure and particle velocity change in a distance-dependent manner. Both sound directions (-x, x) stay separate in the relative phase between pressure and particle velocity. **c–f,** Results of a simple model used to illustrate level and phase of pressure and motion along the horizontal x-axis of the inner tank. The idealized speaker activations in the different sound configurations are modeled as sinusoidal monopole sound sources (see pressure equation in **b**) located 6 cm away from the origin with frequency f = 780 Hz and speed of sound c = 1500 m/s. Acceleration is calculated from the spatial pressure gradient along the x-axis. The top rows are phasor

representations of pressure or motion at five positions along the horizontal axis as indicated in the left cartoon. ILDs and ITDs are computed across a distance of 600 µm centered at the origin. See also Fig. 4 for a summary on ILD and ITD across sound configurations. **c,** In the single speaker configuration, P-ILD, P-ITD, M-ILD, and M-ITD could be interpreted as rightward cues by the fish. M-ITD is even smaller than P-ITD as motion phase propagates slower than pressure phase in the near field. **d,** In the trick configuration, P-ILD and P-ITD are inverted, while M-ILD and M-ITD remain unchanged as compared to the single speaker configuration. **c–d,** See Methods section on interaural cues for comparable P-ILD and M-ILD measurements in our setup. Note that we model monopoles in open water here, but the actual speakers are extended pressure sources in a tank. **e–f,** In both the pressure configuration and the particle motion configuration P-ILD, P-ITD, M-ILD, and M-ITD are zero or undefined.
<figure><img src="image_67.png" /></figure>

<figure><img src="image_68.png" /><figcaption>Extended Data Fig. 9 | See next page for caption.</figcaption></figure>
**Extended Data Fig. 9 | Directional behaviour is present after lateral line ablation.** *a-d*, Same as in Fig. 2c–f but for lateral line ablated (LLA) *D. cerebrum*. *a*, Illustration of six different stimulus configurations. The cartoons of resulting traces illustrate the pressure and particle velocity relationship, idealized for a plane wave scenario. *b-d*, Centered startle trajectories after sound playback. *b*, Directional escapes away from a single speaker for a positive polarity waveform (66% escapes, two-sided binomial test: *p* = 0.0003 across *n* = 142 startles in 67 fish with startle trials in this sound configuration) and negative polarity waveform (69% escapes, two-sided binomial test: *p* = 1 × 10⁻⁶ across *n* = 164 startles in 69 fish with startle trials in this sound configuration). *c*, Absence of directional bias in the positive polarity condition with only pressure (51%, two-sided binomial test: n.s.; *p* > 0.05) and only particle motion (48%, two-sided binomial test: n.s.; *p* > 0.05). See Extended Data Fig. 8 for opposite polarity results. *d*, Single speaker playbacks pooled over both polarities evoke a directional escape (67% escapes, two-sided binomial test: *p* = 1 × 10⁻⁹ across *n* = 306 startles in 72 fish with startle trials in this sound configuration). Selective inversion of pressure polarity by an additional opposing pair of speakers along the orthogonal axis inverts the relative polarity between pressure and particle velocity. This implements the trick condition, in which the fish is tricked into performing startles approaching the active speaker (68% approaches, two-sided binomial test: *p* = 2 × 10⁻⁹ across

*n* = 271 startles in 73 fish with startle trials in this sound configuration). *b-d*, Percentages indicate the fraction of startle displacements into the direction of the blue arrow. The heatmaps are normalized and smoothed 2D histograms over endpoint positions of the shown trajectories (gray for single stimulus data, blue for pooled data). *e-f, D. cerebrum* stained with DASPEI, imaged under an epifluorescence microscope. Images are constructed from several fields of view (gray dotted lines). Scale bar: 2 mm *e*, In untreated fish, neuromasts are visible as green dots (see exemplary white arrows). *f*, After neomycin treatment, no neuromasts are visible, confirming lateral line ablation. *g-h*, Functional indicators of lateral line ablation. *g*, Neomycin-treated fish have significantly more wall contacts within the secondary escape trajectory, i.e. within ~500 ms (independent one-sided t-test *p* = 1 × 10⁻⁵, after blind manual classification of the first ten startles in each recording into startles with and without wall contact. Some fish had less than ten startles, hence black dots are not always multiples of 10%). *h*, Left: In the light, untreated control shoals of mixed-sex adult *D. cerebrum* did not capture significantly more Artemia than neomycin-treated shoals (*n* = 4 shoals in each condition, independent one-sided t-test, *p* = 0.17). Right: In the dark, untreated shoals captured more Artemia (*p* = 0.024, independent one-sided t-test, *n* = 4 shoals in each condition), possibly using their intact lateral line sense. *g-h*, vertical gray stripe on violin plot indicates quartiles of data.
<figure><img src="image_69.png" /><figcaption>**b** Extracting motion amplitude and phase from line scan imaging</figcaption></figure>

<figure><img src="image_70.png" /><figcaption>**c** Correcting motion phase for sound propagation across field-of-view</figcaption></figure>

<figure><img src="image_71.png" /></figure>
**Extended Data Fig. 10 | Method for extracting 2D phase maps of tissue motion with laser scanning confocal reflectance microscopy. a,**

Experimental setup (see Methods section on vibrometry for details on the acoustic stimulation system). **b,** Illustration of phase map extraction: (i) A single bead is located within the field-of-view and (ii) oscillates in the horizontal x direction during acoustic stimulation. (iii) The bead is imaged with a laser scanning microscope, with each pixel being acquired at a different time. Line-scan and acoustic stimulation are synchronized to probe the bead at four different phases (blue: 0, red: $\pi/2$, yellow: $\pi$, green: $3\pi/2$). (iv) Data are reshaped to reconstruct the full movie of the bead motion. (v) The bead displacement is computed using a cross-correlation-based algorithm (see Methods). (vi) The

amplitude and phase of the first Fourier component of the bead displacement are extracted and plotted respectively with hue and color. **c,** (i) The sound phase and consequently the bead displacement phase is drifting when the pressure wave propagates along the horizontal x direction. Different motion phases are then detected for objects at different locations in the field-of-view, although being stimulated by the same sound wave. (ii) This additional phase $\Psi$ is subtracted to yield a phase map with free objects exhibiting the same phase for the same sound stimulation. The final phase relationship between various objects therefore only depends on the mechanical properties of the imaged structure.
a

Compare Fig. 3c

(dorsal, motion of tripus and
scaphium in pressure config.)

<figure><img src="image_72.png" /></figure>

Single loudspeaker
(pressure & p. motion)

<figure><img src="image_73.png" /></figure>

Two in-phase loudspeakers
(pressure)

<figure><img src="image_74.png" /></figure>

Two anti-phase loudspeakers
(particle motion)

<figure><img src="image_75.png" /></figure>

b

Compare Fig. 3d

(lateral, motion of sagitta and
scaphium in pressure config.)

<figure><img src="image_76.png" /></figure>

stimulation axis x

<figure><img src="image_77.png" /></figure>

(x 2)

<figure><img src="image_78.png" /></figure>

<figure><img src="image_79.png" /></figure>

y

<figure><img src="image_80.png" /></figure>

(x 2)

<figure><img src="image_81.png" /></figure>

<figure><img src="image_82.png" /></figure>

c

Compare Fig. 3e

(dorsal, relative motion of
lapillus in particle motion
config.)

<figure><img src="image_83.png" /></figure>

stimulation axis x

<figure><img src="image_84.png" /></figure>

(x 2)

<figure><img src="image_85.png" /></figure>

<figure><img src="image_86.png" /></figure>

y

<figure><img src="image_87.png" /></figure>

(x 2)

<figure><img src="image_88.png" /></figure>

<figure><img src="image_89.png" /></figure>
**Extended Data Fig. 11** | **2D optical vibrometry in three sound configurations.** Phase-amplitude particle velocity maps for 1 kHz component in response to 1 kHz speaker signals in the single speaker configuration (column 1), the pressure configuration (column 2), and the particle motion configuration (column 3). Depicted velocity amplitudes were scaled up by a factor of two in the single speaker condition to more readily compare responses to the pressure and motion configurations, where two speakers created a signal with approximately twice the pressure and motion amplitude. Across all three panels **a–c**, single-speaker maps resemble the addition of the

pressure and motion maps. Sound stimulation was along the x-axis in all cases. See Fig. 3a for the location of depicted views and Fig. 3c–e for names of anatomical regions. **a**, Tripus tip and scaphium contract along the fish's medial-lateral axis (x) in the pressure configuration but not in the particle motion configuration. **b**, The scaphium and the sagitta move at different phases along the fish's ventral-dorsal axis (x) relative to surrounding tissue in the pressure configuration but not in the particle motion configuration. **c**, The lapillus moves with a phase lag to surrounding tissue along the medial-lateral axis (x) in the particle motion configuration but not the pressure configuration.
<figure><img src="image_90.png" /></figure>

<figure><img src="image_91.png" /></figure>

**Extended Data Fig. 12** | Excluding rostro-caudal swim bladder mechanism and cryptic visual cues. **a**, To test the hypothesis that the pair of rostral and caudal swim bladder could implement a P-ILD or P-ITD sensor along the rostrocaudal axis, we filtered for playbacks where *D. cerebrum* was aligned orthogonal to the left-right axis during sound arrival at the fish location. Top: all data, same as Fig. 2f. Middle: filtered for ±10° alignment to y-axis, Single speaker (79% escapes, two-sided binomial test: $p = 6 \times 10^{-16}$, n = 183 startle trials in N = 59 fish), Trick configuration (72% escapes, two-sided binomial test: $p = 3 \times 10^{-7}$, n = 133 startle trials in N = 56 fish). Bottom: filtered for ±3° alignment to y-axis, Single speaker (74% escapes, two-sided binomial test: $p = 1 \times 10^{-4}$, n = 66 startle trials in N = 47 fish), Trick configuration (76% escapes, two-sided binomial test: $p = 8 \times 10^{-4}$, n = 45 startle trials in N = 45 fish). **b–c**, Replication of

the experiment in the dark to exclude the possibility of unknown visual cues. **b**, Top view of the modified setup for playback experiment in the dark. In the dark, the white inner tank walls could no longer prompt *D. cerebrum* to swim orthogonal to the left-right axis. To nevertheless trigger playbacks at the center of the tank with the fish being aligned orthogonally within a 45° cone, we added an additional constraining tank made from thin transparent plastic, thereby increasing the likelihood of triggering playback. **c**, Centered startle trajectory for experiment in the dark. Left: activation of the left speaker leads to rightward startles (71%, two-sided binomial test: $p = 9 \times 10^{-5}$, n = 92 startle trials in N = 43 fish). Right: activation of the right speaker leads to leftward startles (69%, two-sided binomial test: $p = 0.0004$, n = 93 startle trials in N = 43 fish).
# Reporting Summary

Nature Portfolio wishes to improve the reproducibility of the work that we publish. This form provides structure for consistency and transparency in reporting. For further information on Nature Portfolio policies, see our [Editorial Policies](#) and the [Editorial Policy Checklist](#).

## Statistics

For all statistical analyses, confirm that the following items are present in the figure legend, table legend, main text, or Methods section.

n/a | Confirmed

☐ ☑ The exact sample size (n) for each experimental group/condition, given as a discrete number and unit of measurement

☑ ☐ A statement on whether measurements were taken from distinct samples or whether the same sample was measured repeatedly

☐ ☑ The statistical test(s) used AND whether they are one- or two-sided
*Only common tests should be described solely by name; describe more complex techniques in the Methods section.*

☐ ☑ A description of all covariates tested

☑ ☐ A description of any assumptions or corrections, such as tests of normality and adjustment for multiple comparisons

☐ ☑ A full description of the statistical parameters including central tendency (e.g. means) or other basic estimates (e.g. regression coefficient) AND variation (e.g. standard deviation) or associated estimates of uncertainty (e.g. confidence intervals)

☐ ☑ For null hypothesis testing, the test statistic (e.g. F, t, r) with confidence intervals, effect sizes, degrees of freedom and *P* value noted
*Give *P* values as exact values whenever suitable.*

☑ ☐ For Bayesian analysis, information on the choice of priors and Markov chain Monte Carlo settings

☑ ☐ For hierarchical and complex designs, identification of the appropriate level for tests and full reporting of outcomes

☑ ☐ Estimates of effect sizes (e.g. Cohen's *d*, Pearson's *r*), indicating how they were calculated

*Our web collection on statistics for biologists contains articles on many of the points above.*

## Software and code

Policy information about [availability of computer code](#)

**Data collection** We used custom software written in Python 3 and MATLAB (v2019b) (https://github.com/danionella/veith_et_al_2024) as well as LSMAQ for laser scanning microscopy (https://github.com/danionella/lsmaq)

**Data analysis** We used custom software written in Python 3 and MATLAB (v2019b), as well as SLEAP (v1.3), BIOMEDISA (v23), FIJI (v1.5), 3D Slicer (v5.6)

For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Portfolio [guidelines for submitting code & software](#) for further information.

## Data

Policy information about [availability of data](#)

All manuscripts must include a [data availability statement](#). This statement should provide the following information, where applicable:

- Accession codes, unique identifiers, or web links for publicly available datasets

- A description of any restrictions on data availability

- For clinical datasets or third party data, please ensure that the statement adheres to our [policy](#)

Fish trajectory, micro-CT, and vibrometry data have been deposited in the G-Node repository at: gin.g-node.org/danionella/veith_et_al_2024
## Research involving human participants, their data, or biological material

Policy information about studies with human participants or human data. See also policy information about sex, gender (identity/presentation), and sexual orientation and race, ethnicity and racism.

<table><tr><td>Reporting on sex and gender</td><td>N/A</td></tr><tr><td>Reporting on race, ethnicity, or other socially relevant groupings</td><td>N/A</td></tr><tr><td>Population characteristics</td><td>N/A</td></tr><tr><td>Recruitment</td><td>N/A</td></tr><tr><td>Ethics oversight</td><td>N/A</td></tr></table>

Note that full information on the approval of the study protocol must also be provided in the manuscript.

## Field-specific reporting

Please select the one below that is the best fit for your research. If you are not sure, read the appropriate sections before making your selection.

<figure><img src="image_92.png" /></figure>

Life sciences

<figure><img src="image_93.png" /></figure>

Behavioural & social sciences

<figure><img src="image_94.png" /></figure>

Ecological, evolutionary & environmental sciences

For a reference copy of the document with all sections, see [nature.com/documents/nr-reporting-summary-flat.pdf](https://nature.com/documents/nr-reporting-summary-flat.pdf)

## Life sciences study design

All studies must disclose on these points even when the disclosure is negative.

<table><tr><td>Sample size</td><td>The sample size was chosen based on the standards in the field and estimated behavioral variability.</td></tr><tr><td>Data exclusions</td><td>No data were excluded.</td></tr><tr><td>Replication</td><td>We successfully replicated our findings in males (39), females (26), and lateral-line ablated fish (74). Observations shown in Figure 3 c-e were repeated once, four times, and once in other fish, respectively, with similar results.</td></tr><tr><td>Randomization</td><td>There was no group allocation which needed to be randomized.</td></tr><tr><td>Blinding</td><td>There was no group allocation which would have required blinding.</td></tr></table>

## Reporting for specific materials, systems and methods

We require information from authors about some types of materials, experimental systems and methods used in many studies. Here, indicate whether each material, system or method listed is relevant to your study. If you are not sure if a list item applies to your research, read the appropriate section before selecting a response.

### Materials & experimental systems

<table><tr><td>n/a</td><td>Involved in the study</td></tr><tr><td>☑</td><td>Antibodies</td></tr><tr><td>☑</td><td>Eukaryotic cell lines</td></tr><tr><td>☑</td><td>Palaeontology and archaeology</td></tr><tr><td>☐</td><td>Animals and other organisms</td></tr><tr><td>☑</td><td>Clinical data</td></tr><tr><td>☑</td><td>Dual use research of concern</td></tr><tr><td>☑</td><td>Plants</td></tr></table>

### Methods

<table><tr><td>n/a</td><td>Involved in the study</td></tr><tr><td>☑</td><td>ChIP-seq</td></tr><tr><td>☑</td><td>Flow cytometry</td></tr><tr><td>☑</td><td>MRI-based neuroimaging</td></tr></table>
## Animals and other research organisms

Policy information about studies involving animals; ARRIVE guidelines recommended for reporting animal research, and Sex and Gender in Research

**Laboratory animals** Wildtype adult (≥ 4 month old) Danionella cerebrum

**Wild animals** No wild animals were used in this study

**Reporting on sex** Findings apply to both sexes and sexes have been indicated explicitly

**Field-collected samples** No field-collected samples were used in this study

**Ethics oversight** All animal experiments conformed to Berlin state, German federal and European Union animal welfare regulations and were approved by the LAGeSo, the Berlin authority for animal experiments.

Note that full information on the approval of the study protocol must also be provided in the manuscript.

## Plants

**Seed stocks**

Report on the source of all seed stocks or other plant material used. If applicable, state the seed stock centre and catalogue number. If plant specimens were collected from the field, describe the collection location, date and sampling procedures.

**Novel plant genotypes**

Describe the methods by which all novel plant genotypes were produced. This includes those generated by transgenic approaches, gene editing, chemical/radiation-based mutagenesis and hybridization. For transgenic lines, describe the transformation method, the number of independent lines analyzed and the generation upon which experiments were performed. For gene-edited lines, describe the editor used, the endogenous sequence targeted for editing, the targeting guide RNA sequence (if applicable) and how the editor was applied.

**Authentication**

Describe any authentication procedures for each seed stock used or novel genotype generated. Describe any experiments used to assess the effect of a mutation and, where applicable, how potential secondary effects (e.g. second site T-DNA insertions, mosiacism, off-target gene editing) were examined.