Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Featured ArticleResearch Articles, Systems/Circuits

Brain-Wide Mapping of Water Flow Perception in Zebrafish

Gilles Vanwalleghem, Kevin Schuster, Michael A. Taylor, Itia A. Favre-Bulle and Ethan K. Scott
Journal of Neuroscience 20 May 2020, 40 (21) 4130-4144; DOI: https://doi.org/10.1523/JNEUROSCI.0049-20.2020
Gilles Vanwalleghem
1The Queensland Brain Institute, The University of Queensland, St. Lucia, Queensland 4072, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Gilles Vanwalleghem
Kevin Schuster
2School of Biomedical Sciences, The University of Queensland, St. Lucia, Queensland 4072, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael A. Taylor
3Australian Institute of Bioengineering and Nanotechnology, The University of Queensland, St. Lucia, Queensland 4072, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Itia A. Favre-Bulle
1The Queensland Brain Institute, The University of Queensland, St. Lucia, Queensland 4072, Australia
4School of Mathematics and Physics, The University of Queensland, St. Lucia, Queensland 4072, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ethan K. Scott
1The Queensland Brain Institute, The University of Queensland, St. Lucia, Queensland 4072, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ethan K. Scott
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Information about water flow, detected by lateral line organs, is critical to the behavior and survival of fish and amphibians. While certain aspects of water flow processing have been revealed through electrophysiology, we lack a comprehensive description of the neurons that respond to water flow and the network that they form. Here, we use brain-wide calcium imaging in combination with microfluidic stimulation to map out, at cellular resolution, neuronal responses involved in perceiving and processing water flow information in larval zebrafish. We find a diverse array of neurons responding to head-to-tail (h-t) flow, tail-to-head (t-h) flow, or both. Early in this pathway, in the lateral line ganglia, neurons respond almost exclusively to the simple presence of h-t or t-h flow, but later processing includes neurons responding specifically to flow onset, representing the accumulated displacement of flow during a stimulus, or encoding the speed of the flow. The neurons reporting on these more nuanced details are located across numerous brain regions, including some not previously implicated in water flow processing. A graph theory-based analysis of the brain-wide water flow network shows that a majority of this processing is dedicated to h-t flow detection, and this is reinforced by our finding that details like flow velocity and the total accumulated flow are only encoded for the h-t direction. The results represent the first brain-wide description of processing for this important modality, and provide a departure point for more detailed studies of the flow of information through this network.

SIGNIFICANCE STATEMENT In aquatic animals, the lateral line is important for detecting water flow stimuli, but the brain networks that interpret this information remain mysterious. Here, we have imaged the activity of individual neurons across the entire brains of larval zebrafish, revealing all response types and their brain locations as water flow processing occurs. We find neurons that respond to the simple presence of water flow, and others attuned to the direction, speed, and duration of flow, or the accumulated displacement of water that has passed during the stimulus. With this information, we modeled the underlying network, describing a system that is nuanced in its processing of water flow simulating head-to-tail motion but rudimentary in processing flow in the tail-to-head direction.

  • calcium imaging
  • lateral line
  • light sheet microscopy
  • microfluidic
  • sensory
  • zebrafish

Introduction

The lateral line system provides water flow sensation in fish and amphibians, informing behaviors as varied as noncontact assessment, predator evasion, hunting, and rheotaxis (Partridge and Pitcher, 1980; Montgomery et al., 1997; McHenry et al., 2009; Olszewski et al., 2012; Suli et al., 2012; Butler and Maruska, 2015). Its sensory neuromasts are spread across the body of the animal and are composed of bundles of hair cells whose polarized organization makes them direction selective either in the anterior–posterior or dorsal–ventral axis (Dijkgraaf, 1963; Ghysen and Dambly-Chaudière, 2004). Signals from neuromasts are carried to the lateral line ganglia (LLGs), and from there into the neuropil of the medial octavolateral nucleus (MON).

Each afferent neuron in the LLGs selectively innervates hair cells of a given polarity (Nagiel et al., 2008; Faucherre et al., 2009; Ji et al., 2018). As a result, although they may receive input from hair cells spanning multiple neuromasts, afferent neurons are direction selective (Nagiel et al., 2008; Faucherre et al., 2009; Ji et al., 2018). The lateral line system can also be topographically divided into the following two regions: the anterior lateral line, comprising cranial neuromasts whose afferent neurons are found in the anterior LLG (Raible and Kruse, 2000), and the posterior lateral line, comprising trunk and tail neuromasts innervated by the posterior lateral line ganglion (pLLG) afferent neurons (Metcalfe et al., 1985; Raible and Kruse, 2000; Schuster et al., 2010).

Tract-tracing studies and single-cell labeling have shown that, in addition to the MON and Mauthner neurons, the afferent neurons from the lateral line ganglia project to the eminentia granularis (McCormick, 1989; Alexandre and Ghysen, 1999; Bleckmann, 2008; Liao and Haehnel, 2012), which further projects to the torus semicircularis (TS), homologous to the inferior colliculus, and the deep layers of the optic tectum (OT), homologous to the superior colliculus (McCormick and Hernandez, 1996; Bleckmann, 2008). Lateral line information is finally conveyed to the telencephalon from the TS via the lateral preglomerular nucleus (Wullimann, 1997; Bleckmann, 2008). In goldfish, neurons in the MON respond at the onset of the stimuli or show sustained responses for the duration of the stimuli (Mogdans et al., 1997; Ali et al., 2010; Künzel et al., 2011), and similar responses have been described in the TS (Plachta et al., 1999; Ali et al., 2010).

The algorithm underlying rheotaxis has recently been described theoretically: larvae appear to integrate the difference of flow velocity across their sides to direct turning behavior into the flow of water (Oteiza et al., 2017). This simple decision-making rule aligns them, after a few swim bouts, within the center of the flow velocity gradient. Existing anatomic and physiological descriptions of lateral line processing, described above, are insufficient to explain how this basic vectorial information is encoded in the brain, and how the adaptive decision to turn toward the direction of the flow is generated. Indeed, our current understanding of the brain-wide flow-processing network is insufficient to explain how animals execute a host of behaviors that rely on the lateral line.

In this study, we have used a microfluidic device to deliver water flow stimuli to immobilized zebrafish larvae while imaging a genetically encoded calcium indicator throughout the brain. This has produced whole-brain, cellular-resolution datasets of neural activity in larvae exposed to different speed and orientations of water flow. Our data reveal functional categories of neurons with direction-selective responses to water flow, flow onset, and the accumulation of water flow through time, as well as direction-nonselective neurons responding to water flow or flow onset in either direction. We observe lateral line processing in numerous regions across the brain, including but not limited to those that have been described previously. The results inform a new model of brain-wide water flow processing in larval zebrafish.

Materials and Methods

Zebrafish larvae and calcium imaging.

All procedures were performed with approval from the University of Queensland Animal Welfare Unit in accordance with approval SBMS/378/16/ARC.

Zebrafish (Danio rerio) larvae, of either sex, carrying the transgene elavl3:H2B-GCaMP6s or elavl3:H2B-GCaMP6f (Chen et al., 2013) were maintained at 28.5°C on a 14 h ON/10 h OFF light cycle. Adult fish were maintained, fed, and mated as previously described in the study by Westerfield (2000). All experiments were conducted in nacre mutant larvae of the TL strain (Chen et al., 2013). Larvae at 6 days postfertilization (dpf) were immobilized in 2% low-melting point (LMP) agarose (Progen Biosciences) and imaged at 5 Hz on a custom-built selective plane illumination microscope (SPIM, Thompson et al., 2016; Taylor et al., 2018). The tail was freed below the swim bladder and the fish was paralyzed with 100 μm tubocurarine (tubocurarine hydrochloride pentahydrate, Sigma-Aldrich; CHEBI:9774). The flow speed was chosen to be below the startle threshold and was not strong enough to elicit passive movement of the tail. The fish were mounted in a custom-built microfluidics chamber (Fig. 1A) hooked up to a syringe pump (NE-1000X, New Era Pump Systems) that drove water flow, and the flow was triggered using an Arduino board to synchronize with the imaging. The microscope setup was isolated from vibrations on a micro-g lab table (catalog #63–534, TMC). In each larva, horizontal planes were imaged in either the dorsoventral or ventrodorsal direction, at 20 µm increments from the dorsal-most neurons in the brain to the deepest brain region that could be clearly imaged using SPIM. For most larvae, this resulted in a stack of images spanning ∼240 µm dorsoventrally and capturing the entire rostrocaudal and lateral extents of the brain. This means that most of the brain was robustly sampled, but that some of the deepest regions (composing the ventral-most 50 µm, approximately) may have been missed in some larvae and may therefore be underrepresented in our dataset. All image acquisition and stimulus presentation was controlled by the microManager software (Edelstein et al., 2010).

Motion artifacts caused by slow drift of the image or by spontaneous movements by the larva were corrected in Fiji (RRID:SCR_002285) using a rigid body transformation in StackReg (Thévenaz et al., 1998). We used the CaImAn toolbox for the fluorescence signal extraction, using 4000 components for the initialization (Pnevmatikakis et al., 2016).

Behavioral experiments.

The larvae were mounted as described above, but the imaging was with a 4 × 0.1 numerical aperture Olympus microscope objective (PLN 4×) placed below the microfluidic device, and a tube lens L4 projected onto a Basler aca1920 camera recording movements at 100 frames per second (fps; Favre-Bulle et al., 2017). Videos were processed using a median filter of 1 pixel radius and background subtraction in Fiji, the tail motion was measured in Fiji and further analyzed in MATLAB (RRID:SCR_001622). The motions were high-pass filtered with a 5 Hz filter to identify fast tail flicks, and we measured the frequency of tail beats during flow stimulation with the findpeaks function on the high-pass filtered data.

Microparticle imaging velocimetry.

Water flow profiles were measured using particle imaging velocimetry (PIV). The flow chamber and pump were set up in a configuration similar to that in imaging experiments, with water replaced by a suspension of 4.5-µm-diameter polystyrene latex microspheres (Thermo Fisher Scientific) in water at a concentration of 2.5 × 10−4 (by weight). The chamber was imaged at 100 fps, and the flow field was calculated from 800 sequential images recorded while the pump generated a flow of 10 µl/s (Fig. 1B). The images were analyzed using the PIVlab software package version 2.00 (Thielicke and Stamhuis, 2014).

Clustering and registration of calcium responses.

Linear regression, clustering, location, and quantitative analyses were performed in MATLAB, using custom scripts, which can be made available on request. For the first set of experiments, we built regressors, with an average GCaMP response, for each repeat of the flow stimulus onset and offset. The coefficient of determination (r2) of the linear regression models was used to select responsive neurons, and we chose a 0.1 threshold based on the r2 distribution of our models to allow for conservative filtering of the data.

Direction selectivity index (DSI) was calculated as in the study by Grama and Engert (2012). Embedded Image Resph-t is the response (z-scored) to head-to-tail (h-t) flow and Respt-h is the response to tail-to-head (t-h) flow.

The coordinates of individual neurons were normalized by creating a template from all imaged larvae using ANTs (Advanced Normalization Tools; Avants et al., 2010, 2011). The resulting template was then registered onto the Z-brain atlas (Marquart et al., 2017), as previously described (Favre-Bulle et al., 2018).

To generate the nodes for the network analysis, the region of interest (ROI) coordinates of each functional cluster were clustered using k-means with a squared Euclidean distance in an unsupervised manner. The number of spatial clusters was constrained by the necessity to have at least 10 ROIs from three different fish contributing to each node; otherwise, the number of spatial clusters was lowered until such a criterion was met. As a result of the light sheet coming from the right side of the fish, we observe a slight asymmetry in the distribution of the nodes, this is explained by the slightly better image quality (signal-to-noise ratio and out-of-focus light). We used the average z-scored response of each node to compute a pairwise correlation matrix, and this matrix gave us the correlation metric between each pair of nodes. From this metric, we computed an undirected binary graph where an edge would link two nodes if the correlation between them was above 0.6. Graph metrics were extracted using the brain connectivity toolbox (RRID:SCR_004841; Rubinov and Sporns, 2010). Specifically, we computed the degree (number of edges from each node) and participation coefficient (proportion of different functional clusters that each node is connected to) of each node belonging to a functional cluster and averaged them.

For the assessment of neural decoding of features, we used the Neural Decoding toolbox (https://github.com/KordingLab/Neural_Decoding/) with linear (Wiener filter), nonlinear (Wiener Cascade), and a dense neural network models (Glaser et al., 2017). We trained each model on the first third, tested on the second third, and validated on the last part of the integrators neuronal responses to the complex flow. We report the Pearson's correlation between the features predicted by each model and the actual features (speed, time, and distance) of the complex flow. To obtain replicates, we split the neuronal population into 10 random groups that were processed in the same way.

Experimental design and statistical analysis.

The design included 13 fish for the control group, and 3 fish for the neomycin-treated group in the first set of experiments, as well as 7 fish for the second set of experiments. Linear regressions were performed using the fitlm function of MATLAB with default parameters. The k-means clustering was performed on the filtered data (as described above, r2 > 0.1) with 50 clusters using the city block distance with five replicates.

We used a Kolmogorov–Smirnov test in MATLAB to compare the directional selectivity of each cluster against a normal distribution with the following parameters, µ of 0 (no selectivity) and σ of 0.2 corresponding to the mean σ of the experimental distributions.

We used an ANOVA (after testing for normality using Anderson–Darling tests) in GraphPad prism (RRID:SCR_002798) to compare the mean degrees and participation coefficients. The parameters were set to use the Geisser–Greenhouse correction for nonsphericity, assuming matched measures for each functional clusters and using the Tukey correction for multiple comparisons. We report the adjusted p value, as well as mean differences and mean SE of the differences.

We used Friedman tests in GraphPad Prism to compare the correlation and movement frequencies between baseline and stimuli. We used Dunn's correction for the comparison between areas under the curve and report the adjusted p values.

For the neural decoding, we used a Wilcoxon test to compare the correlation between time and distance in our pseudo-replicates, and we report the one-tailed p value.

Results

Combining microfluidics and light sheet microscopy to detect water flow responses

Brain-wide mapping of water flow networks requires the controlled application of flow stimuli in a setting where calcium imaging can take place. With this goal, we designed a 3D-printed microfluidic device compatible with our custom-built light sheet microscope (Fig. 1A; Favre-Bulle et al., 2018; Taylor et al., 2018). We embedded 6 dpf larvae in LMP agarose to immobilize them, cut the agarose away from their tails below the swim bladder, and inserted them into the device's larger channel. We then used a dual-syringe pump to deliver smooth consistent stimuli, pushing water through the small side channels and pulling it out of the center channel to produce h-t flow, or reversing this to produce t-h flow. To gauge the resulting water flow stimuli, we suspended polystyrene latex beads in water and imaged their movements while the pump was active (Fig. 1B). This preparation produced an orderly flow that peaked at 2 mm/s in the middle of the chamber where the tail of the fish was positioned.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

A microfluidics device for presenting larval zebrafish with water flow. A, 3D-printed device, with one large channel and two small channels through which water flows. The flow of water through the channels (for an h-t flow stimulus) is indicated by arrows, and the position and orientation of the larva is indicated with a black arrow. B, Micro-PIV analysis of the water flow in the chamber, fish is for illustration only. The traces (bottom) indicate the speed of the flow along the width of the channel. C, Raster plot of the 65,000 fluorescent traces extracted from the 13 fish (left), the direction and timing of the flow is indicated by the colored bars at the bottom (colors match B). D, The raster plots for 6510 flow responsive neurons (after r2 filtering), grouped by similarity to one another.

We first performed a broad unbiased search for water flow-responsive neurons across the brain. We started with a simple stimulus train of alternating h-t and t-h flow (10 s each, with 10 s of rest between stimuli) with a single consistent flow rate of 2 mm/s, well below the threshold for eliciting startles and within the bounds of larval swim speed (McHenry et al., 2009; Severi et al., 2014). Thirteen fish expressing a pan-neuronal nuclear-targeted GCaMP6s (Chen et al., 2013) were paralyzed with tubocurarine and exposed to the above stimulus in the microfluidics device, while whole-brain calcium imaging proceeded. Following motion correction and segmentation of our images into ROIs corresponding to neurons (see Materials and Methods), we obtained fluorescence traces for ∼63,000 ROIs (Fig. 1C). We registered these ROIs to the Z-brain atlas, a reference atlas for 6 dpf larval zebrafish in which brain regions were identified through both anatomic location and expression patterns, allowing responses to be registered spatially from animal to animal, and against the brain regions delineated in Z-brain (Randlett et al., 2015). As previously described (Favre-Bulle et al., 2018), we used linear regression to identify ∼6500 ROIs that were consistently and specifically responsive to flow stimuli (Fig. 1D).

A basic characterization of flow-responsive neurons across the larval zebrafish brain

Our next goal was to categorize the different types of responses that water flow elicits across the brain. From our 6500 water flow-responsive ROIs, we used k-means clustering with a city block distance to identify eight distinct functional clusters of water flow-responsive neuron (Figs. 2, 3). One broad category of flow-responsive ROIs corresponded to onset-specific neurons (Fig. 2A, Movie 1) and comprised the following three functional clusters: “bidirectional onset” neurons that had similar responses to flow in either direction; “h-t onset” neurons responsive selectively responsive to h-t flow; and “t-h onset” neurons specifically responsive to t-h flow. In each case, there was an initially strong response to the stimulus that quickly settled back to baseline, even as the flow stimulus persisted. Each also showed a weak response to the offset of the opposite direction stimulus, or to the cessation of stimuli of both directions in the case of bidirectional onset neurons. This could result from either a small amount or backflow in the chamber at the end of these stimuli or a weak sensitivity to relative changes in flow rate (such that the cessation of t-h flow is equivalent to the onset of h-t flow).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Whole-brain response profiles to the water flow stimuli. A, The first group of responses show a rapid response to the onset of the stimuli, indicated by colored rectangles (mean responses, SD shaded area). Some responses are weakly direction selective (top), the others are selective to h-t (middle) or t-h (bottom) water flow. Outlines indicate telencephalon (Telen), OT, TS (dotted), and MON. Scale bar, 100 µm. Rotation in Multimedia 1. B, The second group shows ongoing responses for the duration of the stimuli (mean responses, SD shaded area). These can again be subdivided by their direction selectivity: weakly h-t (top), strongly h-t (middle), or strongly t-h (bottom). Rotation in Multimedia 2. C, The third group shows a slow buildup of the GCaMP signal, peaking at the end of the stimuli. Only two such clusters fit into this category: one strongly selective for h-t flow and the other strongly selective for t-h flow. Rotation in Multimedia 3. D, Scatter plot of the average DSI of each cluster. DSI = 0.16 ± 0.13, 0.83 ± 0.07, −0.66 ± 0.23, 0.30 ± 0.19, 0.93 ± 0.03, −0.87 ± 0.07, 0.9 ± 0.07, −0.84 ± 0.1; mean ± SD, n = 13 fish). The p values of KS tests against normal distribution (µ = 0, σ = 0.2) = 2.4981e-02, 1.8228e-04, 5.9128e-12, 5.8589e-12, 7.6683e-13, 5.9769e-11, 5.9260e-12, 6.3809e-12. E, Scatter plot of the average time to maximum intensity of each cluster, both for the h-t (left) and t-h (right) stimuli (2.3 ± 1.1, 4.3 ± 0.9, 2.4 ± 0.7, 4.7 ± 0.5, 10.2 ± 1.2, 13.7 ± 3.4, 13.7 ± 4.6, 8.3 ± 5.6; mean ± SD, n = 13 fish). F, tSNE plot of all the members from the eight clusters, using correlation as a distance. 3-dimensional rotations of the dorsal views shown in A, B, and C can be found in Movies 1, 2, and 3, respectively.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Heat maps of the eight clusters of flow-responsive neurons. The z-scored fluorescent traces of all the ROIs composing each of the eight clusters from Figure 2. Colored rectangles at the bottom indicate the timing of the stimuli (magenta for h-t, green for t-h).

ROIs belonging to the three “onset” clusters were distributed broadly across the brain, and overlapped extensively with one another, especially in the hindbrain. There were, however, differences in the density of different clusters in different brain regions, with bidirectional onset neurons enriched in the TS (Figs. 2A, 4).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Distribution of the eight clusters across flow-responsive brain regions. Normalized distribution of each cluster from Figure 2 across the major flow-responsive brain regions.

A second broad category of ROIs responded to the presence of flow with activity that persisted for the duration of the stimulus (Fig. 2B, Movie 2). These included clusters of “bidirectional on,” “h-t on,” and “t-h on” neurons, paralleling the direction selectivity shown by the onset neurons described above. Again, the spatial distributions of these three clusters were highly overlapping, especially in the hindbrain and MON, but regional enrichments for specific clusters were also evident (Figs. 2B, 4). Notably, the directional-selective on clusters (of both orientations) were the nearly exclusive occupants of the pLLG, reinforcing the observations that the responses of these neurons reflect activity in the neuromasts (Pichler and Lagnado, 2019).

The final broad category of ROI, comprising the “h-t integrator” and “t-h integrator” clusters, showed persistent responses that grew in intensity throughout the duration of the stimulus, then gradually returned to baseline after the stimulus ended (Fig. 2C, Movie 3). Our analysis did not reveal bidirectional integrators. The h-t integrators were more numerous, but the distributions of h-t and t-h integrators were similar, with a majority of neurons in these clusters appearing in the tectum or the hindbrain.

Quantitative analyses of the response properties of these clusters bore out the initial observations presented above. Six of the eight clusters showed strong direction selectivity (according to the DSI; Eq. 1), with only the two bidirectional clusters failing to diverge significantly from 0 (Fig. 2D). In terms of the time that ROIs took to reach their peak responses to water flow stimuli, the three onset clusters showed fast precise responses, the on clusters had average peaks near the midpoint of the 10 s stimuli, and the integrator clusters peaked at or just after the termination of the stimulus (Fig. 2E).

Clustering is subject both to the inadvertent separation of equivalently responding ROIs (overclustering) and to the combining of distinctly responding ROIs (underclustering). While it is not possible to exclude the possibility of either of these types of error in our clustering, we performed a t-SNE (t-distributed stochastic neighbor embedding for long) dimension reduction, using the correlation distance and not the city block distance of our k-means, to visualize the members of the eight clusters (Fig. 2F) relative to one another. Reassuringly, the eight clusters show clear segregation, with the h-t- and t-h-selective clusters separated, and the bidirectional clusters in between the two. Within each of these direction-selective groups, the subclasses of functional responses are also separated. These observations provide support for the validity of our eight clusters as biologically relevant subtypes of flow-responsive neuron.

As a final control to confirm that we are observing bona fide lateral line responses, we treated larvae with neomycin to ablate their lateral line hair cells (Harris et al., 2003). Following lateral line ablation, we observed a drastic drop in flow responsiveness across the brain, as judged by the distribution of the coefficient of determination of all ROIs to the same simple stimulus train used above (Fig. 5A). In attempting to cluster the responses that did occur, we found very small numbers of neurons fitting into two noisy clusters approximately corresponding to the h-t on and t-h on clusters from our initial experiments, albeit with weaker response strength (Fig. 5). The small number of responsive ROIs, along with their poor response fidelity, indicates that the neomycin is blocking ∼90% of the flow responses. This suggests that the lateral line system rather than, for instance, proprioception, is responsible for the numerous, strong, and diverse responses that we describe above.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Neomycin-treated whole-brain responses to water flow. A, Normalized distribution of the r2 values from the linear regression of each fluorescence trace for every ROI to the water flow stimuli, WT in orange and neomycin-treated, with their lateral line ablated, fish in blue (n = 3). B, Clustered responses, after filtering on the r2 values as we did in Figure 2 from the neomycin-treated fish, average responses (left), and heatmaps (right). Colored rectangles at the bottom indicate the timing of the stimuli (magenta for h-t, green for t-h).

Correlation patterns across the brain-wide water flow network

To gain a better sense of the network organization of the observed responses, we applied the tools of graph theory to our results. We performed an unsupervised clustering of the spatial localization of all the responsive ROIs, for each functional cluster, ensuring that the nodes would be represented in all fish (Fig. 6D). This unsupervised approach confirmed the prevalence of mostly h-t and t-h on clusters in the pLLG (Fig. 6D, dashed outlines). We then generated an average correlation matrix between all the nodes across our fish, with correlation <0.6 ignored, to generate an undirected graph (Fig. 6A,E). The graph metrics of degree and participation across clusters were consistent at other thresholds of correlation, which shows that the 0.6 correlation threshold is not special and the following results are generalizable independent of the threshold (Fig. 6C). A null model was generated by applying an amplitude-adjusted Fourier transform to the average time series of each node, which preserves the amplitude and temporal structure of the data while shuffling the temporal indices (Fig. 6B; Theiler et al., 1992). This null model loses any structure and shows no correlation higher than our threshold (highest correlation, 0.33).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Graph analysis of flow-responsive network. A, Correlation matrix for the 126 functional nodes. Color bar indicates correlation coefficient. B, Amplitude adjusted Fourier transform null model correlation matrix. C, We computed the degree and participation of the graph from Figure 3 using lowering correlation thresholds from 0.9 (light gray) to 0.5 (black). D, Spatial localization of the representative 126 nodes from our eight functional clusters, colored as in Figure 2. The light sheet enters the brain from the right, which gives us as better signal on that side and explains the observed assymetry. E, Graph showing functional edges between nodes with correlation coefficients exceeding 0.6. F, The sizes of the nodes represent the degree (number of connections) of each node. Larger nodes have more connections. G, Participation coefficient between functional classes (fraction of different classes with which each node is connected), with larger nodes having stronger participation scores. H, Bar plots of mean degree and mean participation (±SD, n = 13 fish) across all nodes from each of the eight functional cluster. p Values and mean differences for the multiple comparisons can be found in Tables 1 and 2 for degrees and participation, respectively. I, Scatter plots of the difference of mean degree and mean participation from H, for each fish, comparing the h-t clusters to their equivalent clusters in other directions.

View this table:
  • View inline
  • View popup
Table 1.

Multiple comparison of mean degree of nodes from Figure 3

View this table:
  • View inline
  • View popup
Table 2.

Multiple comparison of mean participation of nodes from Figure 3

We observe that the network is biased toward the h-t on response cluster, which is significantly more connected as measured by the degree of each node (the number of nodes to which each node, regardless of function, is correlated; Fig. 6F,H,I, Table 1). We also compared the participation coefficient (the proportion of different functional clusters with which each node connects) and again the h-t-specific clusters had a significantly higher participation coefficient than t-h-specific ones (Fig. 6G,I, Table 2). Overall, this suggests a network that encodes various properties of flow stimuli in both directions, but with a greater network-wide weighting on the detection and processing of h-t flow.

Flow-responsive neurons encode specific properties of water flow stimuli

Our identification of “on,” “onset,” and “integrator” neurons, along with functional clusters within these categories responsive to h-t flow, t-h flow, or both reveals the cardinal response types to flow information in the larval zebrafish brain. Our simple stimulus train, however, limits our ability to gauge how these neurons encode key properties of flow stimuli, including their speed, duration, and overall magnitude. Specific questions include whether any of these neurons encode speed, whether onset neurons respond only to the onset of stimuli or to changes in speed while a stimulus is ongoing, and whether integrator neurons accumulate the duration or volume of flow. To explore these response properties further, we designed a more complex stimulus train (Fig. 7A, top) with stimuli of various strengths, directions, and durations, and we paired or separated these stimuli temporally to isolate responses to stimulus onset, change in flow speed, and stimulus offset. To increase the temporal resolution of our calcium imaging during these rapid transitions, we used GCaMP6f, a GCaMP variant with faster kinetics than the GCaMP6s used above (Chen et al., 2013). This will also allow us to validate that some of the response profiles were not artifacts of the slow rise and decay of GCaMP6s.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Brain-wide representations of the complex features of water flow. A, A more complex stimulus train was used to assess the representation of speed, acceleration, and volume through the whole brain. At the top is a visual representation of the stimulus train (forward simulated motion, or h-t flow, is positive in magenta; t-h flow is negative in green). The bottom shows the raster plot of the responses similar to the clusters of the simple stimulus from Figure 2. Colored rectangles indicate to which functional clusters from Figure 2 the ROIs belong to. B, Onset functional clusters do not respond to stepwise change in speed. Top, Average response for the bidirectional, t-h-specific, and h-t specific onset clusters. Bottom, Visualization of their localization (same colors as Fig. 2; blue is bidirectional, green is t-h, and red is h-t; n = 7 fish). Scale bar, 100 µm. C, The h-t integrators encode a combination of speed and duration. Top, Average response of integrators to the flow of increasing the simulated distance traveled; labels indicate the length of each stimulus. The third and fourth, as well as fifth and sixth are repeats of the stimulus (respectively, the fifth and last h-t stimuli, and the second and eighth from A). Below are the average areas under the curves for the different conditions (n = 7 fish, mean ± SD). The only significant differences were between the first or second stimulus and the last two or last three stimuli, respectively (adjusted p values for Dunn's test: first vs sixth = 0.0108; first vs seventh = 0.0444; second vs fifth = 0.00,444; second vs sixth = 0.0004; second vs seventh = 0.0022). Bottom, Visualization of their localization. Rotation in Multimedia 4. D, Speed encoding of the flow speed. Top, Average response of the two categories of speed-encoding neurons, either 1 mm/s specific (blue) or 2 mm/s specific (red). Bottom, Visualization of their localization. Rotation in Multimedia 5. 3-dimensional rotations of the dorsal views shown in C and D can be found in Movies 4 and 5, respectively.

We first used the eight previous functional clusters from our simple h-t/t-h stimulus train to predict responses to our more complex stimulus train (Fig. 7A, top), we obtained ROIs with similar profiles across our seven fish (Fig. 7A, raster plot). However, we failed to identify responses corresponding to the t-h integrators, which may indicate that these were an artifact of the GCaMP6s, of our clustering approach, or both.

Using the response profiles to this more complex stimulus train (Fig. 7A, top), we asked specific questions regarding the salient stimulus features responsible for activity in the onset neurons. Are they triggered specifically at the start of a flow stimulus, or are they sensitive to any changes in flow speed? To address this, our complex stimulus train included a 30 s stimulus that changed from 1 to 2 mm/s after 10 s, stayed at 2 mm/s for 10 s, and then dropped back to 1 mm/s for the final 10 s. Clustering of the onset responses for a restricted window of time surrounding these changes of speed did not reveal consistent responses to the change from 1 to 2 mm/s, indicating that the onset cells represent the start of flow, but not a change of speed (Fig. 7B).

Our next question involved the integrators and which of the following flow features they integrate: time, or a combination of time and speed, which would be equivalent to the volume of water flow experienced by the fish. As mentioned above, we failed to identify the t-h integrator responses, but the h-t integrator showed responses that correlated strongly to a combination of the time and speed of the flow, which may represent the amount of fluid displacement of the flow stimulus, or simulated motion from the fish frame of reference (Fig. 7C, Movie 4). Specifically, we measured the areas under the response curves during different combinations of flow speed and duration for the h-t flow. This showed that the area under the curve does not appear to increase simply for either time or speed, but that responses from integrators are impacted both by the time and speed of the stimulus (Fig. 7C). We speculate that the lack of significant difference between the 1 mm/s for 10 or 20 s stimulus could be an artifact of increased responses to the first stimulus presented in the stimulus train (the 10 s stimulus, as also observed for several clusters in Fig. 2A; Favre-Bulle et al., 2018), or could be explained by a floor effect. To validate which feature the integrators may best encode, we used the Neural Decoding toolbox (Glaser et al., 2017) and trained three models on the integrators neuronal responses (see Table 5). In all three models, the distance feature (combining speed and duration) best explained by the responses, ahead of either speed or duration individually. The difference between time and distance encoding was significant in all three models [Wiener filter (linear), p = 0.002; Wiener cascade (nonlinear), p = 0.002; and dense neural network, p = 0.001; one-tailed Wilcoxon test).

In larval zebrafish, different water flow speeds trigger different behaviors ranging from rheotaxis to startles (McHenry et al., 2009; Olszewski et al., 2012; Suli et al., 2012; Oteiza et al., 2017), suggesting that the speed of flow stimuli is represented in the brain. To search for speed-encoding ROIs, we observed responses that occurred during the pooled stimuli at 1 or 2 mm/s of flow in our complex stimulus train. In analyzing the resulting data, we identified putative speed encoding ROIs by selecting all responses that differed by >2 SDs between the two speeds, and then clustered these ROIs to separate noisy or inconsistent responses from those possibly encoding speed. This targeted analysis yielded two functional clusters (Fig. 7D, Movie 5), one with approximately threefold stronger responses to flow at 2 than 1 mm/s (680 ROIs across seven larvae), and the other with twofold stronger response to 1 than 2 mm/s (598 ROIs across seven larvae). Although these two clusters also responded to t-h flow, their responses to t-h flow were not speed sensitive. The ROIs belonging to these clusters were mostly found in the MON, the hindbrain, and the TS, with some scattered in other brain regions (Fig. 7D). The 1 mm/s responses were clustered in the middle of the MON, while the 2 mm/s responses were at the rostral and caudal edges of the MON.

Finally, to ensure that the observed responses were sensory and not motor responses, we imaged fish exposed to the complex stimulus. We embedded the fish in the same manner that they were embedded for the calcium-imaging experiment, but without tubocurarine treatment. We did not observe any significant changes in tail movements during stimulation compared with baseline (Table 3), and those movements were also not significantly correlated to the stimuli (Tables 4, 5). We speculate that the flow speed to which our larvae were exposed may be too slow to elicit significant movements. These results show that most of our observed water flow responses were sensory, rather than premotor or motor commands. An alternative explanation lies in the recent observation that the glia may silence futile behaviors in embedded larval zebrafish; the fish may eventually give up swimming during our settling period before imaging can begin (Mu et al., 2019).

Movie 1.

A 3-dimensional rotation of the dorsal view shown in Figure 2A.

Movie 2.

A 3-dimensional rotation of the dorsal view shown in Figure 2B.

Movie 3.

A 3-dimensional rotation of the dorsal view shown in Figure 2C.

Movie 4.

A 3-dimensional rotation of the dorsal view shown in Figure 7C.

Movie 5.

A 3-dimensional rotation of the dorsal view shown in Figure 7D.

View this table:
  • View inline
  • View popup
Table 3.

Tail beat frequency during flow stimulation

View this table:
  • View inline
  • View popup
Table 4.

Correlation of motor activity to flow stimulation

View this table:
  • View inline
  • View popup
Table 5.

Pearson's correlations of neural decoding of complex features by the forward integrators

Discussion

Overview of flow-responsive neurons across the brain

By performing brain-wide calcium imaging during controlled water flow stimulation, we have identified various categories of responsive neurons, the stimulus features that they are attuned to, and their distributions across the brain. This allows several observations about the structure and functional architecture of the water flow-sensing network of the brain. We have found that the direction of flow is already represented in the pLLG, but that none of the more complex stimulus properties (onset or total accumulated flow) are represented in this structure. This observation, which is consistent with past electrophysiological recordings of LLG neurons (Liao and Haehnel, 2012; Haehnel-Taguchi et al., 2014), suggests that these simple directional signals provide the initial input for all subsequent water flow processing in the brain.

The afferent neurons from the pLLG all project to the MON, where we found second-order neurons attuned to all of the flow stimulus properties (presence, direction, onset, and total volume) identified in our study. These diverse response types observed in the MON are intermingled, showing no spatial organization, consistent with the observations that these projections may be topographically, but not functionally, segregated (Faucherre et al., 2009; Liao, 2010; Liao and Haehnel, 2012; Pujol-Martí et al., 2012). Topographical organization of the flow-sensitive structures, with some units showing tight receptive fields and with individual neuromasts represented in restricted areas of the MON (Plachta et al., 1999), would not have been revealed with our setup, which stimulated all trunk neuromasts simultaneously.

We also observed nuanced responses to water flow stimuli in structures downstream of the MON, including the OT, TS, cerebellum, and telencephalon (Fig. 4). These downstream regions are well positioned to integrate water flow stimuli with other modalities such as visual, vestibular, or auditory (Thompson et al., 2016; Vanwalleghem et al., 2017; Favre-Bulle et al., 2018; Constantin et al., 2019). These telencephalic responses may also be involved in decision-making, perhaps through the algorithm described in the study by Oteiza et al. (2017) to align with the flow (rheotaxis). Another region rich in flow-responsive neurons, the anterior hindbrain, has been recently identified as a central location for evidence accumulation and decision-making (Bahl and Engert, 2020). Overall, this shows how the fish could use the downstream regions to integrate the water flow stimuli with other sensory stimuli to inform its behavioral choices. Future studies could combine stimuli from multiple modalities with different saliencies to assess the precise function that each region would play in their processing and integration.

Possible neural mechanisms for speed detection and integration

Our initial observation of integrator neurons was consistent with their representing the accumulation of flow across either time or distance, and our initial stimulus train did not allow us to distinguish between these possibilities. Using a more complex stimulus train, we found evidence for an increased GCaMP response for given neurons at particular speeds (Fig. 7C), which may indicate that the firing rate in these neurons encodes flow speed (Chen et al., 2013). As the fast-encoding ROIs respond (weakly) to the 1 mm/s stimuli, it seems unlikely that these are simply high-threshold ROIs, and this suggests that these signals may actually encode speed. Flow speed could also be represented through temporal coding or coincidence detectors (Mogdans and Bleckmann, 2012), but the temporal kinetics of GCaMP would not have permitted this observation in our preparation. In the future, fluorescent voltage sensors or traditional electrophysiology may allow such observations to be made (Akemann et al., 2010).

The integrators that we have described match nicely with the type of information that a recently developed rheotaxis algorithm would require, in which the accumulated flow across the two sides of the body is compared with drive-turning behavior (Oteiza et al., 2017). We showed that these integrators do not simply encode the duration of the stimulus, but a combination of time and flow speed. The size of the zebrafish larva and the design of our microfluidics chamber precluded the presentation of lateralized stimuli, but such an experiment in the future could test this rheotaxis model concretely (Oteiza et al., 2017).

Coding of detailed information about h-t water flow

By several measures, the brain-wide water flow network appears to be slanted toward the detection and complex processing of water flow like that experienced by a forward-swimming larva. This evidence includes a greater number of h-t-responsive neurons (Fig. 1), more pronounced integration of h-t flow stimuli (Figs. 2, 4), speed-encoding neurons for h-t flow (Fig. 7), and greater functional involvement of h-t-responsive neurons (as measured by degree and participation) in correlation-based network graphs (Fig. 6). Indeed, even the bidirectional response profiles that we have described were slanted toward h-t responsiveness (Fig. 2D).

Combined, these observations suggest that larvae are sensitive to the presence or the onset of water flow in the t-h direction, which may be sufficient, for instance, to detect the water flow produced by the attack of a suction predator. A more nuanced appreciation of the details of h-t flow may subserve not only evasive behaviors, but also the fine control of swimming and navigation. The h-t flow is also the perceived flow direction when fish swim to align with the flow during rheotaxis to remain in place (Olszewski et al., 2012; Oteiza et al., 2017). Viewed this way, the emphasis on h-t flow processing that we observe is unsurprising, since this type of flow would be more common, and likely more ethologically relevant, to larval zebrafish. Recent observations of specific silencing of the h-t flow response at the neuromast level by efferent copies during self-motion also support an ethological relevance to the h-t flow (Pichler and Lagnado, 2020). However, the paralyzing agent we used blocks the cholinergic receptors that mediate this silencing by the motor neurons, and we did not observe swim bouts as a result of our stimuli. How this silencing interacts with our observed responses and how the brain still manages to make the relevant computations would be an interesting topic for future study.

It is also possible that we are underestimating this bias toward perception of a simulated forward movement, given our imaging setup. Since we are using a head-embedded preparation, we are stimulating only the trunk lateral line. Because we are not stimulating the cranial neuromasts, we may be missing responses that would presumably, based on their position on the body of the animal, report on h-t flow. We may also be missing more specialized functions of the cranial neuromasts, such as the use of flow or vibrational information during predatory behavior (Pohlmann et al., 2004; Carrillo and McHenry, 2016).

The encoding and processing of the water flow stimuli may change with age, as the ratio of body size to stimulus features changes, and with the development of canal neuromasts (Van Trump and McHenry, 2008; Wada et al., 2014). As such, studying the ontogeny of the processing would be interesting, especially for the afferent neurons receiving information from the newly developed canal neuromasts, which have been linked to foraging in zebrafish (Carrillo et al., 2019). However, collecting such an ontological dataset would require imaging in ossified zebrafish, whose skulls and pigmentation preclude the use of one-photon light sheet microscopy. Danionella translucida, a promising new fish model, allows imaging of the adult brain and could thus be used for such a developmental study of the lateral line sensory network (Schulze et al., 2018).

To our knowledge, this is the first brain-wide study of lateral line information processing. By observing activity across the brain at cellular resolution, we have uncovered multiple classes of responses representing several features of the flow, and have located each of these response types to its position in the brain. While these results are, in a sense, a comprehensive accounting of lateral line processing in this system, such studies leave many open questions about the circuit-level mechanisms governing the network (Vanwalleghem et al., 2018). The use of GCaMP prevents our judging the order in which events occur across the network, making it difficult to infer the direction of information flow. We also do not know the structures or connectivity of the neurons whose signals we are detecting, and this prevents an anatomic validation of our correlation-based network that we have described using graph theory. In the future, registration of these response types to the morphologies of associated neurons, for instance using databases of neurons across the larval zebrafish brain (Kunst et al., 2019), would allow our activity-based models to be grounded in a plausible anatomic framework. Registration against electron microscopy data could allow us to infer connectivity within the network at the synaptic level (Hildebrand et al., 2017).

More sophisticated stimuli, targeting individual neuromasts or temporally controlled combinations of neuromasts, will be necessary to reveal topography or the computations that allow larvae to detect complex real-world flow patterns. Finally, future functional work, likely using optogenetics, will be necessary to test the resulting models of information flow concretely. Viewed from this perspective, the current study provides a first tantalizing glimpse of this brain-wide network and a departure point for targeted studies of its structural and mechanistic details.

Footnotes

  • The authors declare no competing financial interests.

  • Support for this research was provided by an National Health and Medical Research Council Project Grant (APP1066887) and three Australia Research Council Discovery Project Grants (DP140102036, DP110103612, and DP190103430) to E.K.S., and an EMBO Long-Term Fellowship to G.V.; and a Fellowship from the Human Frontier Science Program to M.A.T. Support was also provided by the Australian National Fabrication Facility, QLD node. We thank the Biological Resources Aquatics Team at the University of Queensland for animal care.

  • Correspondence should be addressed to Gilles Vanwalleghem at g.vanwalleghem{at}uq.edu.au or Ethan K. Scott at ethan.scott{at}uq.edu.au

References

  1. ↵
    1. Akemann W,
    2. Mutoh H,
    3. Perron A,
    4. Rossier J,
    5. Knöpfel T
    (2010) Imaging brain electric signals with genetically targeted voltage-sensitive fluorescent proteins. Nat Methods 7:643–649. doi:10.1038/nmeth.1479 pmid:20622860
    OpenUrlCrossRefPubMed
  2. ↵
    1. Alexandre D,
    2. Ghysen A
    (1999) Somatotopy of the lateral line projection in larval zebrafish. Proc Natl Acad Sci U S A 96:7558–7562. doi:10.1073/pnas.96.13.7558 pmid:10377454
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Ali R,
    2. Mogdans J,
    3. Bleckmann H
    (2010) Responses of medullary lateral line units of the goldfish, Carassius auratus, to amplitude-modulated sinusoidal wave stimuli. Int J Zool 2010:1–14. doi:10.1155/2010/762621
    OpenUrlCrossRef
  4. ↵
    1. Avants BB,
    2. Yushkevich P,
    3. Pluta J,
    4. Minkoff D,
    5. Korczykowski M,
    6. Detre J,
    7. Gee JC
    (2010) The optimal template effect in hippocampus studies of diseased populations. Neuroimage 49:2457–2466. doi:10.1016/j.neuroimage.2009.09.062 pmid:19818860
    OpenUrlCrossRefPubMed
  5. ↵
    1. Avants BB,
    2. Tustison NJ,
    3. Song G,
    4. Cook PA,
    5. Klein A,
    6. Gee JC
    (2011) A reproducible evaluation of ANTs similarity metric performance in brain image registration. Neuroimage 54:2033–2044. doi:10.1016/j.neuroimage.2010.09.025 pmid:20851191
    OpenUrlCrossRefPubMed
  6. ↵
    1. Bahl A,
    2. Engert F
    (2020) Neural circuits for evidence accumulation and decision making in larval zebrafish. Nat Neurosci 23:94–102. doi:10.1038/s41593-019-0534-9 pmid:31792464
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bleckmann H
    (2008) Peripheral and central processing of lateral line information. J Comp Physiol A Neuroethol Sens Neural Behav Physiol 194:145–158. doi:10.1007/s00359-007-0282-2
    OpenUrlCrossRefPubMed
  8. ↵
    1. Butler JM,
    2. Maruska KP
    (2015) The mechanosensory lateral line is used to assess opponents and mediate aggressive behaviors during territorial interactions in an African cichlid fish. J Exp Biol 218:3284–3294. doi:10.1242/jeb.125948 pmid:26491195
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Carrillo A,
    2. McHenry MJ
    (2016) Zebrafish learn to forage in the dark. J Exp Biol 219:582–589. doi:10.1242/jeb.128918 pmid:26889003
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Carrillo A,
    2. Van Le D,
    3. Byron ML,
    4. Jiang H,
    5. McHenry M
    (2019) Canal neuromasts enhance foraging in zebrafish (Danio rerio). Bioinspir Biomim 14:035003.
  11. ↵
    1. Chen TW,
    2. Wardill TJ,
    3. Sun Y,
    4. Pulver SR,
    5. Renninger SL,
    6. Baohan A,
    7. Schreiter ER,
    8. Kerr RA,
    9. Orger MB,
    10. Jayaraman V,
    11. Looger LL,
    12. Svoboda K,
    13. Kim DS
    (2013) Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499:295–300. doi:10.1038/nature12354 pmid:23868258
    OpenUrlCrossRefPubMed
  12. ↵
    1. Constantin L,
    2. Poulsen RE,
    3. Favre-Bulle IA,
    4. Taylor MA,
    5. Sun B,
    6. Goodhill GJ,
    7. Vanwalleghem GC,
    8. Scott
    (2019) Altered brain-wide auditory networks in mutant larval zebrafish. bioRxiv. doi:10.1101/722082. doi:10.1101/722082
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Dijkgraaf S
    (1963) The functioning and significance of the lateral-line organs. Biol Rev Camb Philos Soc 38:51–105. doi:10.1111/j.1469-185x.1963.tb00654.x pmid:14027866
    OpenUrlCrossRefPubMed
  14. ↵
    1. Edelstein A,
    2. Amodaj N,
    3. Hoover K,
    4. Vale R,
    5. Stuurman N
    (2010) Computer control of microscopes using microManager. Curr Protoc Mol Biol Chapter14:Unit14.20.
  15. ↵
    1. Faucherre A,
    2. Pujol-Martí J,
    3. Kawakami K,
    4. López-Schier H
    (2009) Afferent neurons of the zebrafish lateral line are strict selectors of hair-cell orientation. PLoS One 4:e4477. doi:10.1371/journal.pone.0004477 pmid:19223970
    OpenUrlCrossRefPubMed
  16. ↵
    1. Favre-Bulle IA,
    2. Stilgoe AB,
    3. Rubinsztein-Dunlop H,
    4. Scott EK
    (2017) Optical trapping of otoliths drives vestibular behaviours in larval zebrafish. Nat Commun 8:630. doi:10.1038/s41467-017-00713-2 pmid:28931814
    OpenUrlCrossRefPubMed
  17. ↵
    1. Favre-Bulle IA,
    2. Vanwalleghem G,
    3. Taylor MA,
    4. Rubinsztein-Dunlop H,
    5. Scott EK
    (2018) Cellular-resolution imaging of vestibular processing across the larval zebrafish brain. Curr Biol 28:3711.e3–3722.e3. doi:10.1016/j.cub.2018.09.060 pmid:30449665
    OpenUrlCrossRefPubMed
  18. ↵
    1. Ghysen A,
    2. Dambly-Chaudière C
    (2004) Development of the zebrafish lateral line. Curr Opin Neurobiol 14:67–73. doi:10.1016/j.conb.2004.01.012 pmid:15018940
    OpenUrlCrossRefPubMed
  19. ↵
    1. Glaser JI,
    2. Benjamin AS,
    3. Chowdhury RH,
    4. Perich MG,
    5. Miller LE,
    6. Kording KP
    (2017) Machine learning for neural decoding. arXiv. Advance online publication. Retrieved April 13, 2020.
  20. ↵
    1. Grama A,
    2. Engert F
    (2012) Direction selectivity in the larval zebrafish tectum is mediated by asymmetric inhibition. Front Neural Circuits 6:59. doi:10.3389/fncir.2012.00059 pmid:22969706
    OpenUrlCrossRefPubMed
  21. ↵
    1. Haehnel-Taguchi M,
    2. Akanyeti O,
    3. Liao JC
    (2014) Afferent and motoneuron activity in response to single neuromast stimulation in the posterior lateral line of larval zebrafish. J Neurophysiol 112:1329–1339. doi:10.1152/jn.00274.2014 pmid:24966296
    OpenUrlCrossRefPubMed
  22. ↵
    1. Harris JA,
    2. Cheng AG,
    3. Cunningham LL,
    4. MacDonald G,
    5. Raible DW,
    6. Rubel EW
    (2003) Neomycin-induced hair cell death and rapid regeneration in the lateral line of zebrafish (Danio rerio). J Assoc Res Otolaryngol 4:219–234. doi:10.1007/s10162-002-3022-x pmid:12943374
    OpenUrlCrossRefPubMed
  23. ↵
    1. Hildebrand DGC,
    2. Cicconet M,
    3. Torres RM,
    4. Choi W,
    5. Quan TM,
    6. Moon J,
    7. Wetzel AW,
    8. Scott Champion A,
    9. Graham BJ,
    10. Randlett O,
    11. Plummer GS,
    12. Portugues R,
    13. Bianco IH,
    14. Saalfeld S,
    15. Baden AD,
    16. Lillaney K,
    17. Burns R,
    18. Vogelstein JT,
    19. Schier AF,
    20. Lee W-CA, et al
    . (2017) Whole-brain serial-section electron microscopy in larval zebrafish. Nature 545:345–349. doi:10.1038/nature22356 pmid:28489821
    OpenUrlCrossRefPubMed
  24. ↵
    1. Ji YR,
    2. Warrier S,
    3. Jiang T,
    4. Wu DK,
    5. Kindt KS
    (2018) Directional selectivity of afferent neurons in zebrafish neuromasts is regulated by Emx2 in presynaptic hair cells. Elife 7:e35796. doi:10.7554/eLife.35796
    OpenUrlCrossRef
  25. ↵
    1. Kunst M,
    2. Laurell E,
    3. Mokayes N,
    4. Kramer A,
    5. Kubo F,
    6. Fernandes AM,
    7. Forster D,
    8. Dal Maschio M,
    9. Baier H
    (2019) A cellular-resolution atlas of the larval zebrafish brain. Neuron 103:21.e5–38.e5. doi:10.1016/j.neuron.2019.04.034 pmid:31147152
    OpenUrlCrossRefPubMed
  26. ↵
    1. Künzel S,
    2. Bleckmann H,
    3. Mogdans J
    (2011) Responses of brainstem lateral line units to different stimulus source locations and vibration directions. J Comp Physiol A Neuroethol Sens Neural Behav Physiol 197:773–787. doi:10.1007/s00359-011-0642-9 pmid:21479569
    OpenUrlCrossRefPubMed
  27. ↵
    1. Liao JC
    (2010) Organization and physiology of posterior lateral line afferent neurons in larval zebrafish. Biol Lett 6:402–405. doi:10.1098/rsbl.2009.0995 pmid:20181553
    OpenUrlCrossRefPubMed
  28. ↵
    1. Liao JC,
    2. Haehnel M
    (2012) Physiology of afferent neurons in larval zebrafish provides a functional framework for lateral line somatotopy. J Neurophysiol 107:2615–2623. doi:10.1152/jn.01108.2011 pmid:22338025
    OpenUrlCrossRefPubMed
  29. ↵
    1. Marquart GD,
    2. Tabor KM,
    3. Horstick EJ,
    4. Brown M,
    5. Geoca AK,
    6. Polys NF,
    7. Nogare DD,
    8. Burgess HA
    (2017) High-precision registration between zebrafish brain atlases using symmetric diffeomorphic normalization. Gigascience 6:1–15. doi:10.1093/gigascience/gix056 pmid:28873968
    OpenUrlCrossRefPubMed
  30. ↵
    1. McCormick CA
    (1989) Central lateral line mechanosensory pathways in bony fish. In: The mechanosensory lateral line: neurobiology and evolution. (Coombs S, Görner P, Münz H, eds), pp 341–364. New York: Springer.
  31. ↵
    1. McCormick CA,
    2. Hernandez DV
    (1996) Connections of octaval and lateral line nuclei of the medulla in the goldfish, including the cytoarchitecture of the secondary octaval population in goldfish and catfish. Brain Behav Evol 47:113–137. doi:10.1159/000113232 pmid:8680846
    OpenUrlCrossRefPubMed
  32. ↵
    1. McHenry MJ,
    2. Feitl KE,
    3. Strother JA,
    4. Van Trump WJ
    (2009) Larval zebrafish rapidly sense the water flow of a predator's strike. Biol Lett 5:477–479. doi:10.1098/rsbl.2009.0048 pmid:19324627
    OpenUrlCrossRefPubMed
  33. ↵
    1. Metcalfe WK,
    2. Kimmel CB,
    3. Schabtach E
    (1985) Anatomy of the posterior lateral line system in young larvae of the zebrafish. J Comp Neurol 233:377–389. doi:10.1002/cne.902330307 pmid:3980776
    OpenUrlCrossRefPubMed
  34. ↵
    1. Mogdans J,
    2. Bleckmann H
    (2012) Coping with flow: behavior, neurophysiology and modeling of the fish lateral line system. Biol Cybern 106:627–642. doi:10.1007/s00422-012-0525-3 pmid:23099522
    OpenUrlCrossRefPubMed
  35. ↵
    1. Mogdans J,
    2. Bleckmann H,
    3. Menger N
    (1997) Sensitivity of central units in the goldfish, Carassius auratus, to transient hydrodynamic stimuli. Brain Behav Evol 50:261–283. doi:10.1159/000113341 pmid:9360004
    OpenUrlCrossRefPubMed
  36. ↵
    1. Montgomery JC,
    2. Baker CF,
    3. Carton AG
    (1997) The lateral line can mediate rheotaxis in fish. Nature 389:960–963. doi:10.1038/40135
    OpenUrlCrossRef
  37. ↵
    1. Mu Y,
    2. Bennett DV,
    3. Rubinov M,
    4. Narayan S,
    5. Yang CT,
    6. Tanimoto M,
    7. Mensh BD,
    8. Looger LL,
    9. Ahrens MB
    (2019) Glia accumulate evidence that actions are futile and suppress unsuccessful behavior. Cell 178:27.e19–43.e19. doi:10.1016/j.cell.2019.05.050 pmid:31230713
    OpenUrlCrossRefPubMed
  38. ↵
    1. Nagiel A,
    2. Andor-Ardó D,
    3. Hudspeth AJ
    (2008) Specificity of afferent synapses onto plane-polarized hair cells in the posterior lateral line of the zebrafish. J Neurosci 28:8442–8453. doi:10.1523/JNEUROSCI.2425-08.2008 pmid:18716202
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Olszewski J,
    2. Haehnel M,
    3. Taguchi M,
    4. Liao JC
    (2012) Zebrafish larvae exhibit rheotaxis and can escape a continuous suction source using their lateral line. PLoS One 7:e36661. doi:10.1371/journal.pone.0036661 pmid:22570735
    OpenUrlCrossRefPubMed
  40. ↵
    1. Oteiza P,
    2. Odstrcil I,
    3. Lauder G,
    4. Portugues R,
    5. Engert F
    (2017) A novel mechanism for mechanosensory-based rheotaxis in larval zebrafish. Nature 547:445–448. doi:10.1038/nature23014 pmid:28700578
    OpenUrlCrossRefPubMed
  41. ↵
    1. Partridge BL,
    2. Pitcher TJ
    (1980) The sensory basis of fish schools: relative roles of lateral line and vision. J Comp Physiol 135:315–325. doi:10.1007/BF00657647
    OpenUrlCrossRef
  42. ↵
    1. Pichler P,
    2. Lagnado L
    (2019) The transfer characteristics of hair cells encoding mechanical stimuli in the lateral line of zebrafish. J Neurosci 39:112–124. doi:10.1523/JNEUROSCI.1472-18.2018 pmid:30413644
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Pichler P,
    2. Lagnado L
    (2020) Motor behavior selectively inhibits hair cells activated by forward motion in the lateral line of zebrafish. Curr Biol 30:150.e3–157.e3. doi:10.1016/j.cub.2019.11.020 pmid:31866371
    OpenUrlCrossRefPubMed
  44. ↵
    1. Plachta D,
    2. Mogdans J,
    3. Bleckmann H
    (1999) Responses of midbrain lateral line units of the goldfish, Carassius auratus, to constant-amplitude and amplitude-modulated water wave stimuli. J Comp Physiol A Neuroethol Sens Neural Behav Physiol 185:405–417. doi:10.1007/s003590050401
    OpenUrlCrossRef
  45. ↵
    1. Pnevmatikakis EA,
    2. Soudry D,
    3. Gao Y,
    4. Machado TA,
    5. Merel J,
    6. Pfau D,
    7. Reardon T,
    8. Mu Y,
    9. Lacefield C,
    10. Yang W,
    11. Ahrens M,
    12. Bruno R,
    13. Jessell TM,
    14. Peterka DS,
    15. Yuste R,
    16. Paninski L
    (2016) Simultaneous denoising, deconvolution, and demixing of calcium imaging data. Neuron 89:285–299. doi:10.1016/j.neuron.2015.11.037 pmid:26774160
    OpenUrlCrossRefPubMed
  46. ↵
    1. Pohlmann K,
    2. Atema J,
    3. Breithaupt T
    (2004) The importance of the lateral line in nocturnal predation of piscivorous catfish. J Exp Biol 207:2971–2978. doi:10.1242/jeb.01129 pmid:15277552
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Pujol-Martí J,
    2. Zecca A,
    3. Baudoin J-P,
    4. Faucherre A,
    5. Asakawa K,
    6. Kawakami K,
    7. López-Schier H
    (2012) Neuronal birth order identifies a dimorphic sensorineural map. J Neurosci 32:2976–2987. doi:10.1523/JNEUROSCI.5157-11.2012 pmid:22378871
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Raible DW,
    2. Kruse GJ
    (2000) Organization of the lateral line system in embryonic zebrafish. J Comp Neurol 421:189–198. doi:10.1002/(SICI)1096-9861(20000529)421:2<189::AID-CNE5>3.0.CO;2-K
    OpenUrlCrossRefPubMed
  49. ↵
    1. Randlett O,
    2. Wee CL,
    3. Naumann EA,
    4. Nnaemeka O,
    5. Schoppik D,
    6. Fitzgerald JE,
    7. Portugues R,
    8. Lacoste AM,
    9. Riegler C,
    10. Engert F,
    11. Schier AF
    (2015) Whole-brain activity mapping onto a zebrafish brain atlas. Nat Methods 12:1039–1046. doi:10.1038/nmeth.3581 pmid:26778924
    OpenUrlCrossRefPubMed
  50. ↵
    1. Rubinov M,
    2. Sporns O
    (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52:1059–1069. doi:10.1016/j.neuroimage.2009.10.003 pmid:19819337
    OpenUrlCrossRefPubMed
  51. ↵
    1. Schulze L,
    2. Henninger J,
    3. Kadobianskyi M,
    4. Chaigne T,
    5. Faustino AI,
    6. Hakiy N,
    7. Albadri S,
    8. Schuelke M,
    9. Maler L,
    10. Del Bene F,
    11. Judkewitz B
    (2018) Transparent Danionella translucida as a genetically tractable vertebrate brain model. Nat Methods 15:977–983. doi:10.1038/s41592-018-0144-6 pmid:30323353
    OpenUrlCrossRefPubMed
  52. ↵
    1. Schuster K,
    2. Dambly-Chaudière C,
    3. Ghysen A
    (2010) Glial cell line-derived neurotrophic factor defines the path of developing and regenerating axons in the lateral line system of zebrafish. Proc Natl Acad Sci U S A 107:19531–19536. doi:10.1073/pnas.1002171107 pmid:20974953
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Severi KE,
    2. Portugues R,
    3. Marques JC,
    4. O'Malley DM,
    5. Orger MB,
    6. Engert F
    (2014) Neural control and modulation of swimming speed in the larval zebrafish. Neuron 83:692–707. doi:10.1016/j.neuron.2014.06.032 pmid:25066084
    OpenUrlCrossRefPubMed
  54. ↵
    1. Suli A,
    2. Watson GM,
    3. Rubel EW,
    4. Raible DW
    (2012) Rheotaxis in larval zebrafish is mediated by lateral line mechanosensory hair cells. PLoS One 7:e29727. doi:10.1371/journal.pone.0029727 pmid:22359538
    OpenUrlCrossRefPubMed
  55. ↵
    1. Taylor MA,
    2. Vanwalleghem GC,
    3. Favre-Bulle IA,
    4. Scott EK
    (2018) Diffuse light-sheet microscopy for stripe-free calcium imaging of neural populations. J Biophotonics 11:e201800088. doi:10.1002/jbio.201800088 pmid:29920963
    OpenUrlCrossRefPubMed
  56. ↵
    1. Theiler J,
    2. Eubank S,
    3. Longtin A,
    4. Galdrikian B,
    5. Doyne Farmer J
    (1992) Testing for nonlinearity in time series: the method of surrogate data. Physica D 58:77–94. doi:10.1016/0167-2789(92)90102-S
    OpenUrlCrossRefPubMed
  57. ↵
    1. Thévenaz P,
    2. Ruttimann UE,
    3. Unser M
    (1998) A pyramid approach to subpixel registration based on intensity. IEEE Trans Image Process 7:27–41. doi:10.1109/83.650848 pmid:18267377
    OpenUrlCrossRefPubMed
  58. ↵
    1. Thielicke W,
    2. Stamhuis E
    (2014) PIVlab–towards user-friendly, affordable and accurate digital particle image velocimetry in MATLAB. J Open Res Softw 2:e30. doi:10.5334/jors.bl
    OpenUrlCrossRefPubMed
  59. ↵
    1. Thompson AW,
    2. Vanwalleghem GC,
    3. Heap LA,
    4. Scott EK
    (2016) Functional profiles of visual-, auditory-, and water flow-responsive neurons in the zebrafish tectum. Curr Biol 26:743–754. doi:10.1016/j.cub.2016.01.041 pmid:26923787
    OpenUrlCrossRefPubMed
  60. ↵
    1. Van Trump WJ,
    2. McHenry MJ
    (2008) The morphology and mechanical sensitivity of lateral line receptors in zebrafish larvae (Danio rerio). J Exp Biol 211:2105–2115. doi:10.1242/jeb.016204 pmid:18552300
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Vanwalleghem G,
    2. Heap LA,
    3. Scott EK
    (2017) A profile of auditory-responsive neurons in the larval zebrafish brain. J Comp Neurol 525:3031–3043. doi:10.1002/cne.24258 pmid:28599354
    OpenUrlCrossRefPubMed
  62. ↵
    1. Vanwalleghem GC,
    2. Ahrens MB,
    3. Scott EK
    (2018) Integrative whole-brain neuroscience in larval zebrafish. Curr Opin Neurobiol 50:136–145. doi:10.1016/j.conb.2018.02.004 pmid:29486425
    OpenUrlCrossRefPubMed
  63. ↵
    1. Wada H,
    2. Iwasaki M,
    3. Kawakami K
    (2014) Development of the lateral line canal system through a bone remodeling process in zebrafish. Dev Biol 392:1–14. doi:10.1016/j.ydbio.2014.05.004 pmid:24836859
    OpenUrlCrossRefPubMed
  64. ↵
    1. Westerfield M
    (2000) The zebrafish book. A guide for the laboratory use of zebrafish (Danio rerio), Ed 4. Eugene, OR: University of Oregon.
  65. ↵
    1. Wullimann MF
    (1997) The physiology of fishes, Ed 2. Abingdon, UK: Taylor and Francis.
Back to top

In this issue

The Journal of Neuroscience: 40 (21)
Journal of Neuroscience
Vol. 40, Issue 21
20 May 2020
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Brain-Wide Mapping of Water Flow Perception in Zebrafish
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Brain-Wide Mapping of Water Flow Perception in Zebrafish
Gilles Vanwalleghem, Kevin Schuster, Michael A. Taylor, Itia A. Favre-Bulle, Ethan K. Scott
Journal of Neuroscience 20 May 2020, 40 (21) 4130-4144; DOI: 10.1523/JNEUROSCI.0049-20.2020

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Brain-Wide Mapping of Water Flow Perception in Zebrafish
Gilles Vanwalleghem, Kevin Schuster, Michael A. Taylor, Itia A. Favre-Bulle, Ethan K. Scott
Journal of Neuroscience 20 May 2020, 40 (21) 4130-4144; DOI: 10.1523/JNEUROSCI.0049-20.2020
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • calcium imaging
  • lateral line
  • light sheet microscopy
  • microfluidic
  • sensory
  • zebrafish

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Axonal ER Ca2+ Release Selectively Enhances Activity-Independent Glutamate Release in a Huntington Disease Model
  • Sound improves neuronal encoding of visual stimuli in mouse primary visual cortex
  • Decoding of working memory contents in auditory cortex is not distractor-resistant
Show more Research Articles

Systems/Circuits

  • Axonal ER Ca2+ Release Selectively Enhances Activity-Independent Glutamate Release in a Huntington Disease Model
  • Sound improves neuronal encoding of visual stimuli in mouse primary visual cortex
  • Decoding of working memory contents in auditory cortex is not distractor-resistant
Show more Systems/Circuits
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.