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.
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).
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.
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).
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.
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).
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.
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).
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.
Multiple comparison of mean degree of nodes from Figure 3
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.
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).
A 3-dimensional rotation of the dorsal view shown in Figure 2A.
A 3-dimensional rotation of the dorsal view shown in Figure 2B.
A 3-dimensional rotation of the dorsal view shown in Figure 2C.
A 3-dimensional rotation of the dorsal view shown in Figure 7C.
A 3-dimensional rotation of the dorsal view shown in Figure 7D.
Tail beat frequency during flow stimulation
Correlation of motor activity to flow stimulation
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