Abstract
Fish commonly use their lateral line system to detect moving bodies such as prey and predators. A remarkable case is the Mexican blind cavefish Astyanax fasciatus, which evolved the ability to detect nonmoving obstacles. The swimming body of A. fasciatus generates fluid disturbances, the alteration of which by an obstacle can be sensed by the fish's lateral line system. It is generally accepted that these alterations can provide information on the distance to the obstacle. We observed that A. fasciatus swimming in an unfamiliar environment open and close their mouths at high frequency (0.7–4.5 Hz) in order to generate suction flows. We hypothesized that repeated mouth suction generates a hydrodynamic velocity field, which is altered by an obstacle, inducing pressure gradients in the neuromasts of the lateral line and corresponding strong lateral line stimuli. We observed that the frequency and rate of mouthopening events varied with the fish's distance to obstacles, a hallmark of pulsebased navigation mechanisms such as echolocation. We formulated a mathematical model of this hitherto unrecognized mechanism of obstacle detection and parameterized it experimentally. This model suggests that suction flows induce lateral line stimuli that are weakly dependent on the fish's speed, and may be an order of magnitude stronger than the correspondent stimuli induced by the fish's gliding body. We illustrate that A. fasciatus can navigate nonvisually using a combination of two deeply ancestral and highly conserved mechanisms of rayfinned fishes: the mechanism of sensing water motion by the lateral line system and the mechanism of generating water motion by mouth suction.
INTRODUCTION
The Mexican blind cavefish, Astyanax fasciatus (Cuvier 1819) lost its eyes over the course of evolution following multiple independent colonization events of underground caves. It is nonetheless able to successfully avoid obstacles and navigate by utilizing hydrodynamic cues created by its own motion (Dijkgraaf, 1933; Dijkgraaf, 1947; Dijkgraaf, 1963). A moving fish creates hydrodynamic disturbances, namely fluid velocity and pressure fields, which vary with the distance to an obstacle. It has been hypothesized that to form the hydrodynamic image of an obstacle, the fish performs mapping of the hydrodynamic fields into the distance to the obstacle (Campenhausen et al., 1981; Dijkgraaf, 1933; Dijkgraaf, 1947; Dijkgraaf, 1963).
Like other rayfinned fishes, A. fasciatus can sense hydrodynamic disturbances using its lateral line, a specialized system of mechanoreceptors located along the fish's head and body (Bleckmann, 2007; Montgomery et al., 1995). The lateral line consists of mechanoreceptors protruding from the skin (superficial neuromasts) and mechanoreceptors located in canals beneath the skin (canal neuromasts). Superficial neuromasts measure the velocity of the flow very close to the skin, while canal neuromasts measure the water pressure difference between adjacent pores in the canal (Denton and Gray, 1983; Denton and Gray, 1988; Denton and Gray, 1989; Denton and Gray, 1982; Kroese and Schellart, 1992). Although surface and canal neuromasts have overlapping functions, blind cave fish can still locate objects even with disabled surface neuromasts (Montgomery et al., 2001), whereas the ability of these fish to navigate with disabled canal neuromasts is impaired or even lost (AbdelLatif et al., 1990). The pressure difference δp between two adjacent pores of the lateral line is generally considered as a hydrodynamic signal, which also varies with the distance to a wall h (Windsor et al., 2010). For a fish moving steadily in an unbounded fluid, δp is constant. In the presence of an obstacle, δp varies in time. As long as the fish is able to sense and analyze these variations, it is able to detect an obstacle.
Given that the pressure created by the motion of a body in fluid is proportional to the body's square speed, faster Mexican blind cavefish should detect an obstacle sooner and start the avoidance maneuver earlier, at a greater distance from the obstacle. However, turning distances of A. fasciatus from a wall were uncorrelated with its swimming speed (Teyke, 1985; Windsor et al., 2008). The apparent disparity was explained as being the result of a fish's analysis of the relative, rather than the absolute, magnitude of the pressure difference (Windsor et al., 2008). In the scientific literature, the gliding body of a fish is considered as the only source of variation in δp with distance to the obstacle. These pressure variations are denoted here as δp_{B}. To the best of our knowledge, alternative hydrodynamic mechanisms of obstacle detection, in which the stimuli that trigger an obstacleavoiding maneuver are weakly dependent on a fish's speed, are currently unknown.
Virtually all fish use mouth suction, an evolutionarily conserved and deeply ancestral method for prey capture and transport among teleosts (Lauder, 1980; Lauder, 1982; Lauder, 1985). To generate suction, fish rapidly open the mouth and expand the buccal cavity, generating fast flows and steep pressure gradients that extend in front of the mouth (Day et al., 2005; Day et al., 2007; Higham et al., 2006; Holzman et al., 2008; Weihs, 1980). The mouth is then rapidly closed, and the time from the onset of mouth opening until closing (hereafter ‘mouthopening event’) takes 10–100 ms (Gibb and FerryGraham, 2005). The same biomechanical mechanism is used to generate respiratory flows, which are characterized by much slower flows and weak accelerations (Brainerd and FerryGraham, 2006). In the context of suction feeding, it has been previously shown that the fluid velocity fields change with the distance to a wall (Nauwelaerts et al., 2007; Van Wassenbergh and Aerts, 2009). In other words, mouth suction could provide an independent source of variation in δp with distance to the obstacle. Such pressure differences are denoted here as δp_{S}. We hypothesize that mouth
List of symbols and abbreviations
 ∇
 differential operator of gradient
 a, b, c
 semiaxes of an equivalent ellipsoid
 BL
 body lengths
 G
 instantaneous mouth gape
 G_{max}
 peak mouth gape
 h
 instantaneous distance from the plane of the mouth aperture to a plane wall
 h_{0}
 initial distance from the plane of the mouth aperture to a plane wall
 h_{E}
 distance from the plane of the mouth to the wall when the suctioninduced stimulus attains its maximum
 h_{min}
 minimal distance at which fish starts the wallavoiding maneuver
 I
 symbolic notation of mirror image of fish
 L
 fish length
 p
 fluid pressure
 PIV
 particle image velocimetry
 flow rate through mouth aperture
 r
 equivalent radius of the mouth aperture
 R
 symbolic notation of fish
 SNR
 signaltonoise ratio
 t
 time
 T
 time to mouth peak gape
 t_{0}
 initial time of water suction
 U
 fish speed
 V
 fluid velocity vector
 v_{max}
 maximum velocity of peak intake flow through the mouth aperture
 xyz
 coordinates of the fishbodyfixed orthogonal coordinate system O_{1}xyz
 δp
 pressure difference between two adjacent canal pores
 δp_{B}
 pressure difference between two adjacent canal pores induced by body motion
 δp_{S}
 pressure difference between two adjacent canal pores induced by mouth suction
 δs
 distance between two adjacent pores ν water kinematic viscosity
 ξηζ
 coordinates of the mirrorimagefixed orthogonal coordinate system Oξηζ
 ρ
 water density
 φ_{B}
 the fluid velocity potential induced by the motion of the fish's body
 φ_{S}
 the fluid velocity potential induced by mouth suction
 the euclidean distance from the central axis of the disc of the mouth to the fluid point
Our objective was to test whether mouth suction can be used as an additional hydrodynamic mechanism of obstacle detection in Mexican blind cavefish. Specifically, we asked whether mouth suction in A. fasciatus is performed in the absence of food, and whether the rate of mouthopening events is modulated as a function of the environment. Lastly, we quantified the magnitude of the pressure gradient δp_{S} resulting from mouth suction and compared it with that resulting from the translating body δp_{B}.
RESULTS
Experimental results
Swimming freely in a familiar aquarium with no food and under wellventilated conditions (10.6–10.1 mg O_{2}l^{−1}), A. fasciatus had a mean (±s.d.) speed U of 0.75±0.23 body lengths (BL) s^{−1} (equivalent to 52.5±16.2 mm s^{−1}). In the familiar environment, fish opened and closed their mouths at a mean frequency of 0.2±0.23 Hz (N=30, maximum 1.08 Hz). After the obstacles in the tank were shifted to arbitrary locations, mouthopening frequency increased significantly (Fig. 1; ttest, d.f.=60, P<0.001) by ~3.5fold to a mean of 0.7±0.65 Hz (N=31, maximum 1.82 Hz) while swimming speed decreased significantly to 0.55±0.26 BL s^{−1} (equivalent to 38.5±18.6 mm s^{−1}; ttest, d.f.=53, P<0.006). Within treatments, variance in mouthopening frequency was not significantly different between individuals (Bartlett's test for homogeneity of variance, d.f.=2, P>0.8).
In a separate set of experiments we focused on mouthopening events that occurred when the fish swam towards the corner wall of the aquarium (Fig. 2). Fourier analysis revealed that the frequency of these events at a distance h of 110–0 mm from the corner wall was 4.5 Hz (Fig. 2, Fig. 3A). The distribution of mouthopening events with respect to the distance from the corner, perpendicular to the direction of motion, was significantly different from uniform (χ^{2}=61.31, d.f.=10, P<0.001). Mouthopening events were infrequent at a distance greater than 70 mm from the wall; their frequency doubled when the fish were at a distance of 70>h>10 mm, and was more than 6fold higher at the closest distance (10 mm) to the corner (Fig. 3B).
Highspeed videos taken during particle image velocimetry (PIV) experiments (1000 frames s^{−1}) indicate that the maximal mouth gape size G_{max} averaged 2.2±0.7 mm and time to peak mouth opening T averaged 78±51 ms (N=14; Fig. 4A, Table 1). The mouth returned to a closed state after an additional 85±61 ms. Highspeed PIV (supplementary material Movie 2) indicated that peak flow speed at the center of the mouth v_{max} averaged 66±22 mm s^{−1} (N=14; range 22–110 mm s^{−1}; Fig. 4B, Table 1).
Results of the mathematical modeling
Fig. 5A illustrates that the pressure difference δp_{B} induced by the gliding body increases monotonically with decreasing distance to the wall h, as expected (e.g. MilneThomson, 1968) [similar results for a fishlike body of revolution are shown in Windsor et al. (Windsor et al., 2010)]. Fig. 5A also illustrates that in each suction event the pressure difference induced by mouth opening oscillates with the distance to the wall h and reaches positive and negative extrema δp_{S}(h_{E}) at certain distances from the wall, h=h_{E}. The absolute values of the extrema also grow with decreasing distance to the wall. The larger of the extrema may be two orders of magnitude higher than the pressure difference induced by the gliding body in the range of swimming speeds we observed (10–100 mm s^{−1}). In other words, the walldetection stimuli induced by mouth suction may be much higher than those induced by the gliding body (Fig. 5A,B).
It follows from the Euler equations of water motion used to simulate the fluid motion (see Eqn 2 in Materials and methods) that the pressure difference induced by the gliding body δp_{B} is proportional to U^{2} (e.g. MilneThomson, 1968) [similar results for a fishlike body of revolution are shown in Windsor et al. (Windsor et al., 2010)]. The results of mathematical modeling also show that, for the flow parameters adopted here, the pressure difference induced by fast mouth opening δp_{S} is weakly dependent on a fish's speed. Therefore, the ratio of maximum δp_{S}(h_{E}) to δp_{B}(h_{E}) is proportional to U^{−2}. From Fig. 5B it is clear that for relatively low swimming speed (U≈20 mm s^{−1}) the detection signal induced by mouth suction is much stronger than the corresponding bodyinduced signal.
Additionally, in 14 mouthopening events for which suction flow speed, swimming speed and distance were available from PIV measurements (Table 1), the detection signals due to mouth suction were ~50fold stronger than those due to the gliding body (Fig. 5C; range 6–233, ttest, t=3.18, d.f.=13, P<0.0071).
DISCUSSION
In this paper we propose a new mechanism for nonvisual navigation in A. fasciatus. We suggest that Mexican blind cavefish repeatedly generate mouth suction flows (Figs 1, 2 and 3) to produce pressure signals that can be used to detect nonmoving obstacles (Fig. 5). Using experimental observations and mathematical modeling we illustrate that for A. fasciatus approaching a wall, the pressure difference δp_{S} between two adjacent canal pores due to mouth suction may be an order of magnitude stronger than the pressure difference δp_{B} due to the gliding body (Fig. 5). These results depend on swimming speed and the magnitude of suctions flows. For instance, for a fish swimming slowly (~20 mm s^{−1}) at a given distance from the wall, the bodyinduced walldetecting signal measured in a neuromast may be an order of magnitude weaker than the suctioninduced walldetecting signal. However, for faster swimming speeds (~100 mm s^{−1}), the two signals may be of the same order of magnitude.
The correlation between swimming speed and the bodyinduced pressure difference implies an innate tradeoff between the strength of the detection signal and the fish's reaction time to avoid an obstacle; speeding up increases the detection signal but reduces the corresponding reaction time. However, the mouthinduced walldetecting signal depends weakly on the speed of a fish, and therefore it can be useful even when the fish is not moving. This novel finding may provide an alternative explanation to the lack of correlation between swimming speed and turning distance in A. fasciatus (Teyke, 1985; Windsor et al., 2008).
In the Earthfixed frame of references the gliding body of a fish and its mouth suction generate water flow in opposite directions. In the presence of an obstacle and in the frame of reference of the fish, the flow velocity induced by the motion of the body varies slowly in time, whereas that induced by the mouth suction varies rapidly. However, the pressure gradient measured by the canal neuromasts is proportional to the flow acceleration. Thus, even if the two velocities are of the same order of magnitude, their accelerations are not. In this respect, a fish that uses its own gliding body for obstacle detection generates a monotonic and slowly varying disturbance. Multiple mouth suction events generate oscillating signals (Fig. 2 and Fig. 5A). Interestingly, canal neuromasts of A. fasciatus are relatively insensitive to slow flow variations but sensitive to oscillating flows (Montgomery et al., 2009; van Netten, 2006). The mouth gape and the related velocity of the suction flows oscillate in time. These oscillating flows can interact with an obstacle to induce an oscillating pressure difference in the canal neuromast. The dominant frequency of these oscillations is of the order of 5–20 Hz (Fig. 5, Table 1), which is close to the lower end of the sensitivity range of the canal neuromasts [10–100 Hz (Montgomery et al., 2009; van Netten, 2006)]. Thus, the stronger signal produced by mouth suction is also more likely to generate a strong neural response due to the unsteady flows it produces. It should be stressed that the advantage of the bodyinduced walldetecting signal is that it acts continuously, whereas the mouthinduced signal is intermittent and sporadic.
The suction flow creates fluid disturbances in a certain steradian angle, which can be derived from Eqn 6 given in the Appendix. The central axis of the angle is perpendicular to the mouth plane. For a fish moving with its head perpendicular to a wall, the central axis of the steradian angle is also perpendicular to the wall. For small ratios of the mouth diameter to the distance to the obstacle, this angle is small. Therefore, our model of mouth suction flow takes into account only the fluid velocity along the longitudinal axis of the steradian angle. Cases when the fish swim toward the wall at an angle very different from 90 deg require reformulation of the problem, which is beyond the scope of the present paper.
To detect underwater objects, cetaceans generate highfrequency pressure disturbances that propagate in water as acoustic waves with length much smaller than the distance to a detected object (Friedlander, 1958). Using mouth suction, A. fasciatus generates lowfrequency pressure oscillation with a wavelength much larger than the distance to the obstacle. In such a case the speed of sound can be considered infinite and the water incompressible (Friedlander, 1958). Thus, the suggested mechanism of obstacle detection using mouth suction by transmitting oscillation signals is clearly not acoustic. However, it does share certain features with echolocation. First, the highfrequency suction flow generated by mouth suction in A. fasciatus can be seen as an active emission of obstacledetection signals, one of the hallmarks of pulsebased navigation mechanisms such as echolocation. Second, a salient feature of echolocation is that animals increase the signal frequency in unfamiliar environments and when approaching a target (Au, 1993; Busnel and Fish, 1980; Simmons et al., 1979; Yovel et al., 2010). Such modulation of mouthopening frequency was observed in our experiments with Mexican blind cavefish (Fig. 1, Fig. 3B).
The frequency of mouth opening in familiar and unfamiliar environments (0.2–0.7 Hz) is well within the reported range for mouth opening for respiration (Brainerd and FerryGraham, 2006; Hughes, 1970; Hughes, 1960). However, our data suggest that the increase in mouthopening frequency in unfamiliar environments is not due to increased metabolic rates.
First, mouthopening frequencies near the corner (4.5 Hz; Fig. 3A) were much higher than those observed during slow swimming in other species (Brainerd and FerryGraham, 2006; Hughes, 1970; Hughes, 1960). Second, in our experiments, mouthopening events were 6fold more likely to occur very close to the corner than at a distance of >7 cm from it, and this modulation is difficult to reconcile with variable metabolic demands over such short distances. Lastly, although we observed a decrease in swimming speed in the unfamiliar environment, A. fasciatus was previously observed to swim in an unfamiliar environment at a speed that was 20–30% higher than that observed in a familiar environment (Burt de Perera, 2004; Teyke, 1985). Based on studies in other fish species, such an increase in swimming speed should result in a maximal increase of ~20% in mouth ventilation rate, and sometimes even a decrease in that rate (Altimiras and Larsen, 2000; Clark and Seymour, 2006; Webb, 1971). It is therefore unlikely that the moderate increase in swimming speed previously observed in A. fasciatus would result in a 3.5fold increase in the rate of mouthopening events, as observed in our experiments.
From a biomechanical perspective, mouthopening events that generate respiratory and suctionfeeding flows are similar (Brainerd and FerryGraham, 2006; Westneat, 2006). Both behaviors consist of buccal expansion in order to generate unidirectional water flow through the mouth and out of the operculum. Hydrodynamically, both behaviors generate unsteady flows that can be modeled as passive flow into an orifice. Thus, both are mouth suction behaviors. Evolutionarily, both behaviors are ancestral and conserved across fishes, and have a shared biomechanical and neurological basis (Brainerd and FerryGraham, 2006; Westneat, 2006). Consequently, we cannot determine whether using mouth suction to detect obstacles evolved from respiratory or predatory flows. The three behaviors (feeding, respiration and obstacle detection), however, impose specific functional demands on mouth suction flows. In respiratory flows, the efficiency of gas exchange is higher under low flow speeds and weak accelerations (Brainerd and FerryGraham, 2006). In contrast, the functional demands for prey capture and obstacle detection by mouth suction are similar, as both depend on fast flows and steep accelerations (Holzman et al., 2012; Holzman et al., 2007; Wainwright et al., 2007). We do not rule out, however, that the mouth suction flows used for obstacle detection are modified respiratory flows, an intermediate behavior between respiratory and preycapturing flows.
The suggested mechanism of obstacle detection in A. fasciatus presents a combination of two ubiquitous mechanisms of rayfinned fishes: the ability to perceive the fluid motion using the lateral line system and the ability to generate unidirectional flows into their mouths. The lateral line system is an ancestral feature of teleosts (Philip et al., 2012), present also in amphibians and Elasmobranchii (Dijkgraaf, 1963). Similarly, the ability to generate mouth suction for prey capture, transport and respiration is ancestral to fishes, and is shared by rayfinned fish, lobefinned fish, amphibians and elasmobranchs (Lauder, 1982; Lauder, 1985; Westneat, 2006). This constitutes a remarkable example of the evolution of a novel mechanism from a combination of ancient mechanisms originally adopted for other functions. It also implies that this mechanism could be widespread, potentially used by other blind fishes, deepsea and nocturnal species.
MATERIALS AND METHODS
Filming of fish
Fish were obtained from the pet trade and housed in 400×300×100 mm aquaria for 4–6 weeks before being transferred to the experimental aquaria. Holding aquaria contained a sponge filter, with no additional obstacles. Fish were fed daily with commercial Tetra flakes. The experiments described below complied with IACUC approved guidelines for the use and care of animals in research at Tel Aviv University, Israel.
To characterize mouth opening in A. fasciatus, fish (total length L=40–70 mm) were introduced into an experimental aquarium (400×300×100 mm) and filmed using a pair of highspeed (125 frames s^{−1}) synchronized digital video cameras (1280×1024 pixels CMOS, Optronics GmBh, Germany) equipped with 60 mm/f2.8 and 24 mm/f1.8 lenses (Nikkor, Japan). The first camera filmed the aquarium from above while the second had a lateral view of one of the aquarium walls and corners (supplementary material Movie 1). Frame rate was selected such that mouth opening was captured by >4 consecutive frames. The cameras were set to cover an approximately 11×11 cm area of the aquarium, near one of the corners. For the first camera, a grid was placed under the aquarium bottom to calibrate the distances. For the second, a ruler was placed on the aquarium wall near the image boundary. Fish were introduced into the aquarium and immediately filmed when swimming along one of the walls, towards the aquarium corner. From the recorded videos we selected sequences in which the fish was continuously visible through its progress towards the wall (5–20 s long, depending of the fish's swimming speed). For each sequence analyzed, we recorded the time of each mouthopening event and the distance of the fish from the aquarium wall in front of it. Overall, 123 mouthopening events were analyzed.
To test whether the frequency of mouth opening in A. fasciatus is modulated in unfamiliar environments, we allowed three A. fasciatus individuals (total length L=40–70 mm) to acclimate for 48 h in a 400×300×100 mm aquarium. The aquarium contained five cylindrical obstacles of diameter 20–40 mm. Obstacles were made of machined polycarbonate, with a heavy steel base. Oxygen was monitored in the aquarium using YSI ProODO oxygen optode (YSI, Yellow Springs, OH, USA). Following the 48 h acclimation period, fish were filmed from the side at 60 frames s^{−1} for 15 min using a commercial Sony HDRCX550 camcorder (Sony, Japan). After filming, the obstacles were shifted to arbitrary locations by moving a large magnet under the aquarium, and the fish were filmed again. From the two filming periods (before and after moving the obstacles) we selected short (2–10 s) intervals in which a fish's mouth was continuously observed in the field of view. The number of mouthopening events in the filming period was quantified, and mouthopening rate was calculated for each sequence. Overall, we analyzed 30 sequences taken before obstacle shifting and 31 after. The duration of filming was selected to ensure a large enough sample size for statistical analysis.
We repeated this experiment, using a camera positioned above the aquarium, to quantify the swimming speed of the fish in familiar and unfamiliar environments. Movies were recorded using a GoPro Hero3 camera (San Mateo, CA, USA) recording at 60 frames s^{−1}. We randomly sampled 27 short clips (~1 s) during 5 min prior to obstacle shifting, and 27 short clips (~1 s) during 5 min immediately after obstacle shifting. We determined the swimming speed in each clip by digitizing the fish's head and dividing the cumulative swimming distance by clip duration.
PIV
To characterize the flows produced in front of the mouth due to rapid mouth opening, we used a flow visualization technique termed PIV. The details of this method are described elsewhere (Holzman and Wainwright, 2009; Raffel et al., 1998) and are discussed here only in brief. In general, this technique is used to obtain instantaneous velocity measurements and derived properties in fluids. The fluid is seeded with small, neutrally buoyant particles and is illuminated such that the particles are visible. The motion of the particles is recorded using a highspeed camera, and is analyzed to calculate the speed and direction of the flow in the field of view.
Fish were allowed to acclimate to their aquaria and to the laser sheet described below for a week prior to the PIV experiments. At the onset of each trial, the obstacles were moved to arbitrary locations in the tank. We focused on an area of 5×5 cm near one of the corners of the tank. Fish that swam voluntarily into this area were filmed and their suction flows were analyzed. A Coherent Magnum II 665 nm, 1500 mW solidstate continuous wave laser (fan angle of 10 deg; Coherent Inc., Santa Clara, CA, USA) was used to produce a laser sheet in the experimental aquarium. The laser sheet, ~5 cm wide and 0.1 mm thick, was parallel to the long wall of the aquarium and ~15 mm distant from it. We focused on an area of approximately 50×50 mm, bordering one of the aquarium's corners. The plane of the sheet coincided with the centerline of fish that swam parallel to the wall. To visualize the flow, the water was seeded with 12 μm silvercoated, hollow glass beads with specific gravity of 1.05 (Potter Industries Inc., Carlstadt, NJ, USA). Fish were filmed in lateral view using a highspeed digital video camera (1000 frames s^{−1}, Photron SA3, Japan) equipped with a 105 mm/f2.8 lens (Nikkor, Japan). Additionally, a camcorder recording at 120 frames s^{−1} (Sony, Japan) captured anterior views of the swimming fish, which were used to verify the orientation and location of the fish within the laser sheet. Sequential images taken during mouth opening, treated as image pairs, were analyzed using a crosscorrelation algorithm in MatPIV (http://folk.uio.no/jks/matpiv/index2.html), an open software toolbox for PIV analysis in MATLAB (MathWorks Inc., MA, USA). Image pairs were analyzed using a windows shifting technique, starting with 64×64 pixel interrogation areas and ending with 16×16 pixel areas (with 50% overlap) after six passes. The crosscorrelation algorithm returned a twodimensional grid of vertical and horizontal velocities and signaltonoise ratio (SNR) for each image pair analyzed.
We extracted data on the magnitude of flow speed at the center of the mouth at each time point. Velocity values with an SNR value lower than 2 were omitted (<10% of the cases). We define peak flow speed as the maximum flow speed observed during a mouthopening event. We analyzed only sequences in which the laser sheet intersected with the midsagittal plane of the fish, as verified with the anterior view camera. In addition to the calculation of flow speed, we determined for each frame the longitudinal and transverse coordinates of the anteriormost points on the fish's upper and lower jaws, using the MATLAB free package DLTdv5 (Hedrick, 2008). We used these landmarks to calculate gape distance and gape angle, the angle between the upper and lower jaw. We also calculated the time to mouth opening, defined as the time it takes the fish to open its mouth from 20% to 95% of the maximal gape observed during the mouthopening event (Holzman et al., 2008). Overall, we analyzed 14 mouthopening events for the three fish.
Statistical analysis
We used Bartlett's test for homogeneity of variance implemented in the software R statistics (R Development Core Team, 2009) to determine whether variation in mouthopening frequency differed significantly between individuals. A fast Fourier transform (FFT) algorithm implemented in MATLAB was used on sequences of repeated mouth opening to estimate the dominant frequency of these events. A chisquare test was used to test whether the frequency of mouthopening events was higher at different distances from the corner by comparing the observed frequencies to a uniform distribution (same number of events in each distance category). Unless otherwise stated, statistical tests were carried out using the software R statistics (R Development Core Team, 2009). Throughout the text, means are shown ±s.d.
A mathematical model of obstacle detection using mouth suction
The primary aim of our modeling was to clarify mathematically which of the physical processes, the steadily gliding body of a fish or its unsteady mouth suction, generates stronger walldetection signals that can be sensed by the lateral line. For this purpose we use the premises of inviscid incompressible fluid, the closedform expressions for the hydrodynamic field generated by a moving ellipsoid (MilneThomson, 1968) and the closedform expressions for the hydrodynamic field generated by a disk of sinks (Lamb, 1932). Using the method of hydrodynamic images, we calculated the pressure that is induced by the wall in two adjacent canal pores for a fish moving perpendicular to the wall. We then compared the magnitude of the pressure differences in two adjacent canal pores generated by the two mechanisms.
The parameters used for our mathematical modeling are the following: fish length L, the speed of fish U, peak mouth opening G_{max}, time to peak gape T, the peak velocity of intake flows at the center of the mouth aperture v_{max} and the minimal distance h_{min} at which a Mexican blind fish starts the wallavoiding maneuver. The typical characteristic scales were adopted here for numerical calculations: L=50 mm, U=20–1000 mm s^{−1}, h_{min}=5 mm (e.g. Windsor et al., 2008). For numerical estimates of T, G_{max} and v_{max} we used the values measured in our filming and PIV measurements (see Results, ‘Experimental results’): G_{max}=2.2 mm, v_{max}=66 mm s^{−1} and T=0.078 s (see Results).
For the range of swimming speeds adopted here, the Reynolds numbers Re=UL/v=2000–10,000 is sufficiently low to assume that a fish's body boundary layer is laminar and thin (Schlichting, 1979; Vogel, 1994). In such a case, according to the basic assumptions of the boundary layer theory, the flow outside the boundary layer can be considered as inviscid and irrotational (potential) (Schlichting, 1979). The distance in front of the mouth aperture, where the viscosity influences the velocity field due to mouth suction √(vT)≈0.3 mm (Tuck, 1970), is much smaller than h_{min}=5 mm. Therefore, we also considered the flow generated by the mouth aperture as inviscid and potential.
Given that pores of canal neuromasts are spaced at relatively small intervals δs with respect to the fish's length (Schemmel, 1967), the pressure difference in a neuromast can be approximated as: (1) where δs is the vector connecting two adjacent pores and ∇p is the pressure gradient. In the Earthbound frame of reference the equation of motion of inviscid incompressible fluid can be described by the Euler and continuity equations (MilneThomson, 1968): (2) (3) where ρ is the fluid density and t is time.
If a body moves in an unbounded fluid with constant velocity U then, in the frame of reference attached to the body, the fluid motion is steady, and the local fluid acceleration ∂V/∂t on the lefthand side of Eqn 2 vanishes. If a body approaches a wall, the geometry of the fluid domain varies in time, the fluid velocity becomes time dependent, and the fluid acceleration in Eqn 2 becomes nonzero. It is important to note that the fluid acceleration may also become nonzero as a result of unsteady mouth suction.
Consider a fish moving with constant velocity in the direction normal to an infinite plane wall. Assume that the velocity field created by a fish in unbounded fluid is known. To satisfy the conditions of impermeability on the wall exactly and on the fish body approximately, we used the method of mirror images (MilneThomson, 1968), where the hydrodynamic combination ‘fish–wall’ is replaced by two identical fish ‘R+I’ swimming in unbounded fluid along the same line with the same velocity U but in opposite directions, where I is a mirror image of R with respect to the wall (Fig. 6). We further assume that fish R is fixed in space and that the water moves with constant velocity U directed from the fish's head to its tail; whereas the fish image I moves with velocity 2U with respect to R and with velocity U with respect to water. The moving fish I generates water disturbances by its gliding body or by its mouth suction. These disturbances create pressure gradients ∇p in the neuromasts of R, which are identical to ∇p created in the same neuromasts by the wall in the combination ‘R+wall’. Once ∇p is defined, the pressure difference in two adjacent pores is also defined by Eqn 2.
We further used two orthogonal coordinate systems, O_{1}xyz fixed in nonmoving fish R and Oξηζ fixed in moving fish I (O_{1} xOξ, O_{1}yOη and O_{1}zOζ). The axes O_{1} x and Oξ are collinear with U. The relationship between the two coordinate systems is ξ=x−2Ut+h_{0}, where h_{0} is the distance between O and O_{1} at the initial moment of time t=0. Assume that the fluid velocity of fluid disturbances generated by the fish image I in an unbounded domain is known in the fixed coordinate system Oξηζ. The total vector of the fluid in the fixed inspace coordinate system O_{1}xyz can be represented as: (4) where V_{B} is the fluid velocity of fluid disturbances pertaining to the gliding body and V_{S} to the mouth suction. Assume that the fish is a slender and streamlined body which generates weak fluid disturbances V_{B,S}<<U. Then, substituting Eqn 4 into Eqn 2, and neglecting higher order terms, we can decompose the pressure gradients into the pressure gradients pertaining to a fish's gliding body, and to that pertaining to the mouth suction.
To calculate the fluid disturbances generated by a fish's gliding body, we replaced it with a threedimensional ellipsoid ξ^{2}/a^{2}+η^{2}/b^{2}+ζ^{2}/c^{2}=1 moving in unbounded fluid with velocity U directed along the axis Oξ (supplementary material Fig. S1).
The algorithm for calculating the corresponding fluid velocities that result from the motion of an ellipsoid (MilneThomson, 1968) is provided in the Appendix.
We consider the flow created by mouth suction as that created by a disk of sinks of uniform density over the plane area contained by a circle η^{2}+ζ^{2}=r^{2}(t) (Lamb, 1932), where the radius of the disk r varies in time as: (5) (6) (7) τ=(t–t_{0})/T, t_{0} is the initial time of suction and α≈5 is an experimental fitting parameter. The total flow rate through the circle can be estimated using PIV measurements of the suction velocity v_{max}f(t) in the center of the mouth aperture: (8) The fluid velocity on the central axis of the disk Oξ can be written as: (9) where ξ^{+} is the distance from the plane of mouth to the point of a fluid normalized with the radius r(t) (see Appendix). Eqn 9 captures the main spatial and temporal characteristics of the velocity field created by mouth suction (Day et al., 2005; Day et al., 2007; Holzman et al., 2008; Muller et al., 1982) (supplementary material Fig. S2). Once the fluid components of the diskinduced flow are known, the diskinduced velocity vector and its derivative with respect to time are also defined.
ACKNOWLEDGEMENTS
We thank Anatoli Zapolski for his help with filming experiments. We are indebted to Amatzia Genin, Matt McHenry, Victor Mishkevich, Peter Wainwright and Yossi Yovel for discussions and comments on the manuscript.
APPENDIX
To calculate the fluid velocities, we used the theory of potential flow (MilneThomson, 1968) where the fluid velocity is calculated as a gradient of the fluid velocity potential. Once the fluid velocity potential is known, the fluid velocity field is also known. The fluid velocity potential of an ellipsoid moving along its longitudinal axis is given as (MilneThomson, 1968): (A1) where: (A2) λ=λ(ξ, y, ζ) is the positive real root of an algebraic equation: (A3) k=−U/(2−A_{0}) and A_{0} is the value of A for λ=0. The fluid velocity components induced by the motion of the ellipsoid can be bound from Eqns A1, A2 as: (A4) The derivatives ∂λ/∂ξ, ∂λ/∂η and ∂λ/∂ζ can be obtained by differentiating Eqn A3 with respect to ξ, η, ζ. For instance, differentiation with respect to ξ gives the following linear equation with respect to ∂λ/∂ξ: (A5) the solution of which is straightforward. Thus, calculating the fluid velocity components induced by image I on the surface of R is reduced to a numeric solution of an algebraic Eqn A3 and numeric calculation of Eqn A2.
We consider the flow created by mouth suction as that created by a disk of sinks of uniform density over the plane area contained by a circle η^{2}+ζ^{2}=r^{2}(t), where the radius r varies in time. The fluid velocity potential φ_{S} induced by the disk of sinks can be written as (Lamb, 1932): (A6) where ξ_{1} is the distance from the disk to the field point, and and J_{0,1} are the Bessel functions of the first kind of a real argument. Differentiation of Eqn A6 with respect to ξ_{1}, η and ζ yields the fluid velocity components induced by the disk of sources. For instance, given that J_{0}(0)=1, the fluid velocity on the central axis of the disk can be written as: (A7) where ξ^{+}=ξ_{1}/r (supplementary material Fig. S2).
FOOTNOTES

Author contributions
R.H. and G.Z. conceived and designed the research; R.H., G.Z. and S.P.F. conducted the experiments; R.H., S.P.F. and G.Z. analyzed the results; G.Z. carried out the mathematical modelling; R.H. and G.Z. wrote the paper.

Competing interests
The authors declare no competing financial interests.

Funding
The study was funded by Israeli Science Foundation (ISF) 158/11 to R.H.

Supplementary material
Supplementary material available online at http://jeb.biologists.org/lookup/suppl/doi:10.1242/jeb.098384//DC1
 © 2014. Published by The Company of Biologists Ltd