Introduction

In vertebrates, cerebellar activity underlies proper posture, defined as the relative orientation of body parts in space ^1–8. The cerebellum integrates sensory information from vestibular (balance), visual, and proprioceptive systems ³. These sensations are transformed into precise and coordinated adjustments in muscle tone and contraction allowing animals to control posture ⁹. Disruptions to mature cerebellar function lead to instability, unsteady gait and a compromised sense of balance ¹⁰. Development of the cerebellum coincides with postural maturation, and early development of the cerebellum has been extensively studied ^11–13. Notably, changes to morphology and activity of the output neurons of the cerebellar cortex, Purkinje cells, are thought to underlie the gradual refinement of motor control ¹⁴. To date, the contributions of developing Purkinje cells to postural control remains poorly understood.

Kinematic quantification by pose estimation in rodents ^15,16has opened a window into cerebellar contributions to postural behaviors in health and disease ^8,17,18. However, terrestrial gait and locomotion are complex and especially difficult to study during development; tracking and analysis is often limited to measures such as the time of head elevation or the duration of walking stance ¹⁹. In contrast, the biophysical challenges of maintaining posture underwater are straightforward to define ^20–22. For example, larval zebrafish balance in the pitch axis (nose-up/nose-down) by timing locomotion to countermand gravity-induced destabilization ^23,24 and by coordinated use of paired appendages (fins) and axial musculature (trunk) ²⁵. The small size and rapid development of the larval zebrafish allow high-throughput measurements of postural control (i.e. pitch-axis kinematics and fin/trunk coordination) from freely swimming subjects ²⁶.

The larval zebrafish is a powerful model to investigate cerebellar development and function ²⁷. Anatomically, the zebrafish cerebellum shares the same circuit structure as the mammalian cerebellum ²⁸. The zebrafish cerebellum is compartmentalized into regions with distinct response properties and output targets ^28–32. Multimodal representations were found in both cerebellar granule cells ^33,34and Purkinje cells ^31,35. Functional assays established a role for the larval zebrafish cerebellum in motor control, sensorimotor integration and predictive neural processing, particularly in response to visual input ^36–44. Finally, brain-wide imaging studies have established balance-relevant sensitivity in the cerebellum, identifying neurons that encode body angle and velocity ⁴⁵ and neurons responsive to direct inner-ear stimulation ⁴⁶. Overwhelmingly, this work has been done in reduced or restrained preparations, limiting insight into the cerebellar contribution to natural behaviors.

Powerful new opto- and chemogenetic ⁴⁷ approaches allow control of particular cerebellar cell types, reviewed in ⁴⁸. Recent work used such activation/inhibition to investigate cerebellar contributions to sensorimotor ^49–52 and non-sensorimotor behaviors ^53–56in health and disease ^57–59. Both approaches come with technical hurdles: optogenetics requires targeting light to the cerebellum, a particular challenge when untethered animals can move freely in depth, while chemogenetics uses bioactive co-factors ⁶⁰. A chemogenetic approach to cerebellar control with a non-bioactive ligand would be a welcome advance, particularly to study posture without visual interference (i.e. in the dark). One validated path forward is to express the rat non-selective cation channel TRPV1 and its ligand capsaicin in zebrafish ⁶¹. The endogenous zebrafish TRPV1 channel is capsaicin-insensitive ⁶², so targeted expression of rat TRPV1 allows cell-type specific control: low-doses of capsaicin can activate sensory and hypothalamic neurons while high-doses are excitotoxic ⁶¹. Capsaicin can be dissolved in water and is readily absorbed by freely-swimming larval zebrafish, sidestepping invasive procedures and the need for visible light. Finally, the conductance of a TRP channel is ∼1000x that of a channelrhodopsin ⁶³ suggesting that even low levels of TRPV1 expression will be biologically effective.

Here we used the TRPV1/capsaicin system to investigate the contribution of cerebellar Purkinje cells to postural behaviors as larval zebrafish swam freely in depth. Both activation and ablation of Purkinje cells could induce changes in pitch axis posture. Ablation in older larvae resulted in bigger disruptions to posture, allowing inference of the functional consequences of cerebellar development. Furthermore, ablation of Purkinje cells in older larvae disrupted the coordination of trunk and paired appendages (fins), impairing vertical navigation. Finally, we could reliably decode pitch-tilt direction from patterns of Purkinje cell activity. Taken together our results establish a clear role for the cerebellum in larval zebrafish postural control, even during the earliest stages of development. More broadly, our work establishes a new method to manipulate cerebellar output while performing quantitative high-throughput measures of unconstrained posture and locomotion. Our data are therefore a step towards defining an ancestral role for the highly-conserved cerebellum in postural control.

Results

A new tool for chemogenetic activation or ablation of Purkinje cells

We used a new tool to control Purkinje cells: the transgenic line Tg(aldoca:TRPV1-tagRFP). Fish in this line express rat TRPV1, a capsaicin-sensitive non-selective cation channel, exclusively in all cerebellar Purkinje cells (Figures 1A and 1B) ⁶⁴.

Figure 1:

Endogenous zebrafish TRPV1 channels are insensitive to capsaicin ⁶². Previous descriptions of rat TRPV1 in zebrafish sensory and hypothalamic neurons establish dose-dependent chemogenetic manipulation ⁶¹. We expect low-doses of capsaicin to depolarize Purkinje cells (Figure 1C, left), while high-doses should be excitotoxic (Figure 1C, right).

First, we assayed capsaicin concentrations and incubation times to identify a dose that would achieve long-term depolarization without cell death. We co-expressed a nuclear-targeted calcium indicator, GCaMP6f (Figure 1D) in all neurons using the Tg(elavl3:h2B-GCaMP6f) line for longitudinal imaging of neuronal activity. Previous work used 1 µM of capsaicin for long-term activation ⁶¹. We therefore imaged the cerebellum of Tg(aldoca:TRPV1-tagRFP);Tg(elavl3:h2B-GCaMP6f) fish prior to and 3, 6, and 9 hours after 1 µM capsaicin treatment (Figure 1D). We screened fish for comparable brightness and selected fish that had clearly visible expression but were not overly bright (Methods).

Prolonged exposure to a low dose of capsaicin increased cerebellar activity (Figures S1A and 1E). At each timepoint, TRPV1-expressing cells showed increased intensity relative to a pre-capsaicin baseline, while TRPV1-negative cells did not (Figure 1E, 3/6 hours post 1 µM capsaicin: 28%/20%/TRPV1+ cells F/F₀ > 2; 40 cells from 3 fish vs. 0%/0%/TRPV1-cells F/F₀ > 2; 44 cells from 4 fish; activated cells after 6 h of capsaicin treatment: 0/44 TRPV1-vs. 8/40 TRPV1+; Fisher’s exact test: p = 0.0018)

Different cells showed increased activity at the 3,6, and 9 hour timepoints, and the same cells were differentially active at different timepoints. We interpret this as evidence that 1 µM of capsaicin could sporadically activate subsets of Purkinje cells. Notably, in one fish that had particularly strong tagRFP expression we observed a small number of neurons at the 9 h timepoint with bright, speckled fluorescence suggestive of cell death (Figure S1C). We did not observe any signs of cell death at the 6 h timepoints (6 TRPV1+ fish at 6 h post 1 µM capsaicin). We therefore set an upper limit of 6 h of exposure to 1 µM capsaicin for activation experiments.

Induced activation was reversible, even after prolonged exposure to 1 µM of capsaicin. We tested whether the elevated patterns of neuronal activity that we observed in the presence of capsaicin would return to baseline by imaging cerebellar Purkinje cells in Tg(elavl3:h2B-GCaMP6f) before exposure, after 6 h of 1 µM capsaicin, and 40 min after washout. Relative to baseline, fluorescent intensities increased after 6 h, as in Figure 1E. Importantly, fluorescence returned to baseline levels after 40 min of washout (Figure S1B, 6h post 1 µM capsaicin: 40.9%/TRPV1+ cells F/F₀ > 2; washout: 0%/ TRPV1+ cells F/F₀ > 2; 22 cells from 3 fish; activated cells before 1 µM capsaicin treatment vs. after washout: Fisher’s exact test: p = 1). We conclude that capsaicin-induced activation is reversible after washout.

Exposure to high doses of capsaicin caused rapid axonal degeneration and cell death. We developed a protocol for Purkinje cell lesion: Tg(aldoca:TRPV1-tagRFP) larvae (without GCaMP6f) were imaged at 7 dpf, at 8 dpf after 1 h of 10 µM capsaicin treatment and again at 9 dpf (Figure 1G). Timelapse imaging of the Purkinje cell axons showed rapid degeneration already 15 min after capsaicin treatment started (Figure 1F). Cell numbers rapidly declined after 1 h of 10 µM capsaicin treatment and did not show any signs of recovery at 9 dpf (Figures 1G and 1H) (median [inter-quartile range]; 7 dpf: control 213 [195 – 271] cells vs. pre lesion 282 [255 – 366] cells; 9 dpf: control 218 [201 – 250] cells vs. post lesion 68 [52 – 70] cells; 3 control and 3 lesioned fish).

Consistent with prior work in other cell populations ⁶¹, we found that chemogenetic use of the capsaicin/TRPV1 system can be used to reversibly activate or rapidly ablate cerebellar Purkinje cells in larval zebrafish.

Purkinje cells regulate postural control in the pitch axis

We used our Scalable Apparatus to Measure Posture and Locomotion (SAMPL) to measure posture and locomotion in freely swimming zebrafish ²⁶. SAMPL is a high-throughput videographic approach that measures kinematic parameters of posture and locomotion from fish swimming in a predominantly vertical arena that encourages navigation in depth (Figures 2A and 2B) allowing us to analyze postural changes in the pitch (nose-up/nose-down) axis. Larval zebrafish locomote in discrete bouts of rapid translation (Figure 2B, grey lines). To navigate up/down, fish sequence these bouts while maintaining a nose-up/nose-down pitch ⁶⁵. Notably, climb/dive bouts are defined relative to the trajectory of the bout. Climb/dive bouts can therefore be initiated from either nose-up (positive) or nose-down (negative) postures.

Figure 2:

Nose-up “climb” bouts (Figure 2E) engage both axial musculature of the body and the fins to produce a net upward trajectory while nose-down “dive” bouts (Figure 2J) rely on axial musculature alone and have a net downward trajectory ²⁵. Notably, postural angles after either climb or dive bouts tend to increase, a consequence of restorative rotations that counteract destabilizing torques ²⁴. SAMPL’s automated and high-throughput design yields data with large numbers of observations. To ensure a focus on only the most meaningful differences, we adopted two stringent criteria for significance: p-values <0.05, and an effect size of ≥ 15%. All p-values and effect sizes are reported in Tables 1 to 5.

Table 1.

We used the timing and capsaicin concentrations we had previously validated (Figure 1) to design two behavioral paradigms: one to activate and one to ablate cerebellar Purkinje cells. Experiments were done from 7–9 dpf, and began with a single day without perturbations. Activation was then achieved by exposing Tg(aldoca:TRPV1-tagRFP); Tg(elavl3:h2B-GCaMP6f) fish to two 6h periods of 1 µM capsaicin while they swam freely in the dark (Figure 2C). Alternatively, Purkinje cells were ablated by exposing Tg(aldoca:TRPV1-tagRFP) fish to 10 µM of capsaicin for 1h (Figure 2D). All fish were screened before experiments for comparable levels of tagRFP fluorescence and control and experimental groups were randomly selected. A single experimental repeat consisted of 1-3 apparatus run in parallel with fish from a single clutch of embryos (i.e. siblings). To maintain consistency with genotypes used for validation, the activation and ablation experiments had different backgrounds (i.e. the presence/absence of the elavl3:h2B-GCaMP6f) allele). Because of background variation ²⁶, all comparisons were restricted to control vs. experimental groups within an experimental paradigm over the same time period. We focused our analysis on the pitch axis as the current version of the SAMPL apparatus ²⁶ is not optimized for quantification of roll axis behavior. Across our datasets (Tables 1 to 3) we did not observe meaningful differences between the control and experimental groups in the pre-manipulation period. To avoid adding noise to our estimates of effect size, we therefore report comparisons between control and experimental groups after perturbation. We did not observe global consequences for swimming: swim speed, swim frequency and bout duration were unaffected during Purkinje cell activation (Figures S2A and S2B) or after Purkinje cell lesion. Similarly, raw bout numbers were not different between the control and activation (median [inter-quartile range] 1256 [571 – 1454] bouts vs. 656 [444 – 1485] bouts; p-value 0.46) or lesion groups (2004 [1507 – 2471] bouts vs. 1913 [1556 – 2416] bouts; p-value 0.84, Tables 1 and 2).

Table 2.

Climbing postures were perturbed after both activation and ablation of Purkinje cells. During activation, fish adopted more nose-up postures before and throughout climb bouts. We observed a shift towards more positive values across the distribution of postures before fish initiated a climb bout (Figure 2F). The average climb posture of fish during depolarization was 29% higher than in control fish (Figure 2G, median [95% confidence interval]: 14.7° [14.0 – 15.4°] vs. 19.0° [18.5 – 19.7°], p-value < 0.001, effect size: 29%). Similarly, after Purkinje cell lesion, the average climb bout postural angle increased 36% relative to controls (Figures 2H and 2I, 10.0° [9.5 – 10.7°] vs. 13.6° [13.1 – 14.3°], p-value < 0.001, effect size: 36%).

We observed an unexpected decrease in the climb bout postural angles for control fish in the post-lesion period (from 18.0° [17.6 – 18.4°]to 10.0° [9.5 – 10.7°], Table 2). We do not have an explanation for this particular change and we have confirmed that it is not due to a single outlier experiment (detected outliers: 0/15 pre-lesion control experiments; 0/14 post-lesion control experiments). Notably, if we assess the effect of adding 10 µM capsaicin by comparing the magnitude of the relative difference between pre- and post-lesion periods, normalized to the pre-lesion period, we still see a significant difference (control vs. lesion: -46% vs.-26%). We conclude that, even when accounting for observed changes between control fish at 7 vs. 8 dpf, Purkinje cell ablation modifies climb postures.

Dive bout postures were similarly perturbed after activation, but not ablation of Purkinje cells. Fish adopted more nose-down postural angles before and throughout dive bouts with a leftward shift of the distribution of postures before dive bouts (Figure 2K). Average dive bout posture was 24% more negative than in control fish (Figure 2L, median [95% confidence interval]: -16.6° [-16.9 – -16.1°] vs. -20.5° [-20.9 – -20.1°], p-value < 0.001, effect size = 24%). Purkinje cell lesions at 7 dpf did not shift the average posture for dive bouts (Figures 2M and 2N -11.7° [-11.9 – -11.5°] vs. -11.2° [-11.4 – -11.0°], p-value = 0.002, effect size = -4%).

We interpret these data as evidence that Purkinje cell activity is crucial to ensure that posture during climbs and dives is maintained within a normal range.

Loss of Purkinje cells in older fish results in more global deficits to posture

Over the first two weeks of life, larval zebrafish morphology and postural control strategies develop considerably ²³. These changes are matched by similarly pronounced cerebellar growth ⁶⁶ (Figures 3A and 3B). We observed that the number of Purkinje cells labelled in Tg(aldoca:TRPV1-tagRFP) roughly doubled between 7 and 14 dpf (Figure 3C, median [inter-quartile range] 7 dpf: 282 [255 – 336]; 14 dpf: 662 [591 – 733]). The increase in cell numbers is also evidence that the aldoca promoter continued to drive expression at later stages, allowing us to perform comparative experiments.

Figure 3:

Similar to lesions at 7 dpf, we did not observe any differences in swim speed, frequency or bout duration (Table 3). At 14 dpf, the effects of Purkinje cell lesions on postural angles were more widespread than at 7 dpf and also affected dive postures. We repeated our previous ablation experiments (Figure 2D) between 14–16 dpf, and analyzed climb (Figures S3A and 3D) and dive bouts (Figures S3B and 3E). Loss of Purkinje cells created more widespread behavioral deficits. Specifically, climb bout posture was increased by 20% after Purkinje cell lesion (median [95% confidence interval]: 14.3° [13.8 – 14.8°] vs. 17.1° [16.2 – 17.8°]; p-value < 0.001; effect size: 20%). At 14 dpf we also observed an effect on dive bout postures. After lesion dive bouts postures were 26% more negative (-9.8° [-10.1 – -9.5°] vs. -12.3° [-12.6 – -11.9°]; p-value < 0.001; effect size: 26%).

Table 3.

We conclude that, consistent with morphological growth, Purkinje cells of the cerebellum play a broader role in postural control at 14 dpf than at younger ages.

Purkinje cells regulate speed-dependent fin engagement

To climb, larval zebrafish coordinate fin movements that generate lift with axial rotations that direct thrust (Figure 4A). The greater the axial rotation, the stronger the lift-producing fin movements; this relationship increases as larvae develop ²⁵. Our previous work suggested that Purkinje cells were necessary for such fin-body coordination ²⁵. Here, we observed that fin engagement is speed-dependent, with faster bouts producing greater lift for a given axial rotation (Figure 4B, left, Spearman’s correlation coefficient: 0.2193; p = <0.001; of lift/rotation ratio [mm/deg] versus speed [mm/s] Figure S4A).

Figure 4:

After Purkinje cell ablation, 14 dpf fish produced less lift than expected when they swam fast. We divided swim bouts into three different bins according to their peak speed (slow: 5–7.5mm/s; medium: 7.5–15mm/s; fast >15mm/s) for both control and fish treated with 10 µM capsaicin. We parameterized the relationship between upward rotation and lift by fitting a line to swim bouts for each speed. After capsaicin exposure, the slopes of the medium and fast speed bins were significantly lower (Figures S4 and 4C), reflecting a loss of speed-dependent modulation (median [95% confidence interval]: slope slow: 0.029 [0.028 – 0.031 mm/°] vs. 0.033 [0.032 – 0.034 mm/°], p-value = 0.341, effect size: 6%; slope medium: 0.041 [0.040 – 0.044 mm/°] vs. 0.018 [0.018 – 0.019 mm/°], p-value <0.001, effect size: -34%; slope fast: 0.068 [0.067 – 0.071 mm/°] vs. 0.026 [0.026 – 0.027 mm/°], p-value <0.001, effect size: -62%;). The correlation between speed and fin-lift / rotation ratio is reduced after capsaicin exposure (Figure S4; Spearman correlation coefficient: control 0.2193; 10µM capsaicin: 0.0397; Z-test after z-transformation: p < 0.001) When analyzing fin lift and upward rotation across the three speed bins separately, we only observed significant differences for fast swim bouts, specifically showing a reduction in fin lift and an increase in upward rotationTable 4.

Table 4.

Next, to determine if lift was fin-dependent, we amputated the fins ^25,26 and repeated our experiments. A detailed explanation of how fin amputation affects swim kinematics can be found here ^25,26. We observed a near total loss of lift at all speeds; regardless of the speed bin, the slope of the relationship between upward rotation and lift was indistinguishable from zero (slope slow: 0.036 [0.035 – 0.037 mm/°] vs. -0.005 [-0.005 – -0.004 mm/°], p-value <0.00.1; effect size: -60%; slope medium: 0.055 [0.054 – 0.056] mm/° vs. -0.005 [-0.005 – -0.005 mm/°], p-value <0.001; effect size: -88%; slope fast: 0.068 [0.067 – 0.069 mm/°] vs. 0.013 [0.013 – 0.013 mm/°], p-value <0.001; effect size: -81%). Finally, we examined fin-body coordination in our 7 dpf activation and ablation datasets. In contrast to older larvae, we observed no meaningful changes after activation of Purkinje cells at 7 dpf. For Purkinje cell lesions at 7 dpf we found only the fin body coordination at fast bouts to be affected Tables 1 and 2.

Our data show that loss of Purkinje cells disrupts the speed-dependent increase in fin-mediated lift in older, and to a lesser degree in younger fish. We interpret this finding as evidence that Purkinje cells are indispensable for normal coordination of the fins and body.

Purkinje cells encode pitch direction at both individual and population levels

Our experiments establish that manipulations of Purkinje cells interfere with balance in the pitch axis. We therefore hypothesized that Purkinje cell activity would be modulated by nose-up/nose-down body tilts. We used Tilt In Place Microscopy (TIPM) ²⁶ to measure the response of individual Purkinje cells (Figure 5A) to rapid pitch tilts. Briefly, fish are mounted on a mirror galvanometer and rapidly rotated to eccentric angles (Figure 5B, ±30°nose-up/nose-down).

Figure 5:

We used Tg(aldoca:GAL4);Tg(UAS:GCaMP6s) to label Purkinje cells in the lateral parts of the cerebellum thought to receive vestibular input ^30,39,67. To facilitate identification of the same cells from volumes imaged at both ±30°, we used doubly mono-allelic fish and screened for sparse expression of Purkinje cells. In total, we imaged 43 Purkinje cells from 8 fish. Of those, 31 cells could reliably identified at ±30° and were included in the analysis.

We calculated a directionality index (DI) for all cells and categorized cells as either tuned (|DI| > 0.35) or as untuned (|DI| < 0.35). Some cells showed, on average, higher responses to one direction (Figure 5F), but due to highly variable responses (median [inter-quartile range] Up response range: 23.3 [18.3 – 29.5]; Down response range: 43 [23.7 – 49.2]), they did not exhibit consistent tuning as determined by the directionality index. Individual Purkinje cells showed either directionally-tuned (Figure 5C, n=18) or untuned (Figure 5D, n=13) patterns of responses. Tuned cells were distributed throughout the lateral cerebellum (Figure 5E), and showed a slight preference for nose-down stimuli (12 vs. 6, Figure 5F). We did not observe any systematic differences in the response properties across each experiment from untuned cells (Figure 5G).

While untuned cells did not show overt directional preferences, pooling their responses allowed decoding of stimulus direction. We applied principal component analysis (PCA) to visualize the data and assess whether the trial types even for untuned cells exhibited distinct separation based on their trial identity (nose-up or nose-down). Principal component analysis of the integral of the full responses on each trial from untuned neurons showed near-complete segregation of trial types (Figure 5H). To assay whether there was indeed directional information we trained a decoder (support vector machine) and tested its accuracy on pseudo-populations of untuned cells of different sizes ranging from 3 - 13 cells (Figure 5I). We focused only on untuned cells, as including even a single tuned cell for the population coding will lead to excellent results. Training and test trials were different to avoid over-fitting. Pseudo-populations with more than 3 cells achieved accurate decoding well above chance levels (determined by shuffling trial identity)(median [inter-quartile range] accuracy: 3/5/7/10/13 cells: 0.78 [0.68 – 0.91] / 0.88 [0.70 – 0.88] / 1 [0.84 – 1] / 1 [0.97 – 1] / 1 [1 – 1]). To ensure that the choice of cutoff for the directionality index (0.35) does not bias the decoding results, we also calculated the decoding accuracy using the 3, 5, and 7 least directionally tuned cells based on their directionality index. We again found that more than 3 cells achieve accurate decoding above chance level (3/5/7 least tuned cells: 0.67 [0.64 – 0.78] / 0.86 [0.73 – 0.88] / 0.88 [0.86 – 0.92]).

Older larvae showed additional changes to dive postures after Purkinje cell lesions. We therefore tested if: (1) Purkinje cells in older larvae exhibited differences in the numbers or direction of tuned cells or (2) if population-level decoding accuracy changed. We performed longitudinal TIPM, sampling from zebrafish larvae at 7 and 14 dpf. To improve throughput, we recorded the responses upon return from ±19° stimuli (Figure S5A) (7 dpf: 138/11 cells/fish; 14 dpf: 90/7 cells/fish; of those 23/3 cells/fish were imaged at both timepoints); previous work established that responses upon return to baseline are highly correlated with the response at the eccentric position ⁶⁸. We observed increased fluorescence relative to baseline values in 7 dpf and 14 dpf Purkinje cells upon return from ±19° steps (Figure S5B and Figure S5C). To analyze directional tuning we compared the maximum fluorescence in the first second after return to baseline. The relative number of tuned cells per fish was comparable between 7 and 14 dpf larvae (Figure S5D; median [inter-quartile range] at 7 dpf: 7 [6 – 34%]; 14 dpf: 8 [2 – 19%]; p-value = 0.7763). While most cells were not directionally selective, the preferred direction of tuned cells was different at 7 and 14 dpf: at 7 dpf more Purkinje cells were nose-down tuned (2/31 up/down) but at 14 dpf more cells were nose-up tuned (11/3 up/down Figure S5E; Fisher’s exact test, p-value < 0.001).

We next assayed accuracy of directional encoding of simultaneously recorded untuned cells. We performed principal component analysis for untuned cells at 7 (Figure S5F) and 14 dpf (Figure S5G) and tested decoding accuracy on the untuned cells of individual fish at 7 and 14 dpf. We did not observe differences in decoding accuracy between 7 and 14 dpf larvae (Figure S5H; median [inter-quartile range] 7 dpf: 0.68 [0.63 – 0.83]; 14 dpf: 0.73 [0.65 – 0.79]; p-value = 0.9468). We conclude that cerebellar Purkinje cells can encode pitch direction both at the single neuron and population levels with similar encoding accuracy in young and older larvae.

Discussion

We used a validated chemogenetic tool to investigate the role of cerebellar Purkinje cells in postural behavior as larval zebrafish swam freely in depth. Activation of Purkinje cells could induce changes in pitch axis (nose-up/nose-down) posture. Purkinje cell ablation changed posture, with broader effects in older larvae. Ablation disrupted fin-body coordination responsible for proper climbing. Finally, we could reliably decode pitch-tilt direction from patterns of Purkinje cell activity. We did not observe developmental changes in population coding of direction but found a shift in the tuning direction of Purkinje cells. Taken together our results establish a role for the cerebellum in postural control even during the earliest stages of larval zebrafish development. Our work establishes a new method that combines bidirectional manipulation of cerebellar output and quantitative high-throughput measures of unconstrained posture and locomotion.

Contributions of Purkinje cells to posture

While activation and ablation manipulations both produced biologically meaningful behavior changes, the two experiments were run with different genetic backgrounds and on different generations of the SAMPL apparatus. Conse-quentially, our ability to define precisely what role Purkinje cells play in balance behaviors in larval zebrafish is limited. Activation experiments are particularly laborious as they require thorough pre-screening to ensure adequate brightness levels to achieve sufficient depolarization without excitotoxicity. To maintain consistency in depolarization effects, experiments involving different genetic backgrounds or developmental stages would require additional calcium imaging to confirm that 1 µM capsaicin elicits comparable responses. Hence, we restricted our experiments to 7 days post-fertilization (dpf) and did not extend the study to other developmental time points. Given that the primary purpose of this series of experiments was to establish TRPV1-mediated manipulation of Purkinje cells as a means to investigate postural control, it is beyond the scope of the work to repeat the experiments. Nonetheless, we consider the findings individually below in the context of prior work.

Purkinje cell ablations modifies postural stability. Importantly, the differences we observed were more widespread in older larvae, underscoring the developmental importance of Purkinje cells for balance. Purkinje cell output is inhibitory ^69,70, Purkinje cells in the lateral cerebellum project to vestibular nuclei ^30,39, and Purkinje cells are tonically active ^71,72. We propose that the net effect of Purkinje cell loss would be disinhibition of target nuclei responsible for encoding posture and parameterizing corrective pitch-axis behaviors. While the precise nature of the transformation between larval zebrafish pitch and posture control kinematics is not yet known, loss of cerebellar-targeted nuclei can disrupt postural behaviors ^25,73.

The effects of ablations became more widespread in older larvae. During early development, larval zebrafish grow in volume by roughly an order of magnitude and shift their postural control strategies to better climb/dive as they navigate in depth ^23,25. Unlike climb bouts, changes to postural stability during dives only emerge at 14 dpf. As activation of Purkinje cells produced meaningful changes during dives at 7 dpf, we infer that the delayed emergence of ablation effects does not reflect incomplete integration of Purkinje cells into dive-control circuits. Instead, we propose that this delay reflects the gradual, functional emergence of Purkinje cell control over dives as development progresses. This aligns with previous studies showing that Purkinje cell circuits mature at different rates during development ¹⁴.

Notably, the basal posture during dive bouts decreases in older control fish (Figure 3E) — ablation shifts the posture comparable to its younger state. Future work with our system enables testing of the hypothesis that Purkinje cell output plays a role in setting the postures older fish adopt during dives.

Purkinje cell activation also modifies postural stability. Intriguingly, activation broadened the distribution of observed postures in the same way as ablation. Our imaging assay established that 1 µM of capsaicin would stochastically activate subsets of Purkinje cells. This stochasticity could reflect normal fluctuations in basal levels of activity, or it could arise from cells going in and out of depolarization block ⁷⁴. Synchronized/precisely-timed Purkinje cell output is thought to shape movements ^75–80, though perhaps not for all behaviors ⁸¹. Our imaging suggests that the set of Purkinje cells activated at any one moment in time is limited and random. We therefore propose that the net effect of 1 µM of capsaicin is ultimately disruptive to Purkinje cell synchrony, and thus likely disruptive. Future work could test this hypothesis by intracellular recording from cerebello-recipient populations like the vestibular nuclei ^82,83

Previously, we reported that larval zebrafish coordinate their fins and trunk to climb effectively ²⁵. The relationship between trunk-mediated changes to trajectory (upward rotation) and fin-mediated lift depends on locomotor speed. Here we observed that after Purkinje cell loss, speed-dependent increases in lift with greater trunk rotation are disrupted (Figures S4A to 4C). In the fastest speed bin, we observed an increase in upward rotation and a decrease in average fin lift. In contrast, the medium-speed bin showed no significant changes in average fin lift or upward rotation, yet already displayed coordination deficits. Based on these observations, we argue that Purkinje cell lesions primarily affect coordination, rather than simply reducing one specific parameter such as lift or rotation. As we did not observe any change to locomotor speed after ablation (Tables 1 to 3), we infer that Purkinje cell loss disrupts speed-dependent coordination for climbing. These results extend our original report where a lower-throughput method (photoablation) suggested that Purkinje cell loss impacted the fin-trunk relationship ²⁵. In larval zebrafish, the neuronal substrates for axial speed control ^84–89 and fin engagement ⁹⁰ are known. The potential for whole-brain imaging in larval zebrafish ⁴¹, particularly with high-speed voltage indicators ⁹¹ and cutting-edge modeling approaches ⁹², stands to reveal how Purkinje cell activity comes to coordinate body and fin movements. Importantly, since our behavioral data suggest that Purkinje cell activity impacts fin-trunk coordination more strongly in older larvae, longitudinal approaches will be key to understanding the developmental changes to cerebellar signaling that underlie effective coordination of trunk and limbs.

Encoding strategies for body tilt stimuli

Purkinje cell activity reflects both sensory and motor inputs. One limitation of TIPM is that larvae are immobilized in agarose during tilts. Consequentially, our measurements of Purkinje cell activity are artificially constrained. Nonetheless, a subset of Purkinje cells were unambiguously direction-selective, and a simple decoder could use the activity of non-selective cells to differentiate tilt direction. We infer that vestibular information directly related to pitch axis posture is represented by the Purkinje cell population targeted in our ablation/activation experiments, consistent with broader imaging of cerebellar responses to body tilt ^45,46. Similar to the behavior results, we observed an asymmetry in the tuning direction of Purkinje cells at 7 dpf, with more cells being tuned to the nose-down direction. This asymmetry shifted between 7 and 14 dpf, suggesting developmental changes in how navigation in the pitch axis is processed in the cerebellum. These changes underscore the importance of longitudinal measurements of Purkinje cell activity across early development to understand emergent control of posture.

The ability to decode tilt direction from the collective activity of “untuned” Purkinje cells suggests a role for population coding. Such mechanisms have been proposed for head/body motion ⁹³ and eye movements ^81,94in the primate cerebellum. Population coding requires that multiple Purkinje cells converge onto downstream targets, which is well-established in cerebellar target nuclei ^76,95. In larval zebrafish, Purkinje cells involved in locomotion converge on eurydendroid cells; electrophysiological recordings confirm a many-to-one convergence scheme that could similarly support population coding ³¹. Vestibular-sensitive cells are located in the lateral cerebellum ^45,46, which projects to hindbrain regions that contain vestibular nuclei ⁶⁶. Comparing activity of vestibular nucleus neurons involved in tilt-driven behaviors ^82,96,97before/after TRPV1-mediated ablation would speak to the collective contributions of Purkinje cells.

TRPV1/capsaicin as a tool to study cerebellar contributions to behavior

Our use of TRPV1/capsaicin complements a modern suite of molecular tools to target cerebellar Purkinje cells ⁴⁸. In fish, different experiments have used opsins to excite / inhibit cerebellar Purkinje cells with exceptional temporal precision, establishing functional topography ³⁰ and an instructive role in learning ³¹. TRPV1/capsaicin is a well-validated approach ⁶¹ that permits parametric (i.e. dose-dependent) activation/ablation with a single transgenic line. It does not require light, facilitating dissociation of vestibular from visual contributions without requiring genetically-blind fish as in other studies using excitatory opsins ⁹⁸. Additionally, our approach of using TRPV1/capsaicin for cell ablation offers multiple advantages compared to other lesion methods. For example, Killer Red requires extended exposure to high light intensities and mounting of single fish ²⁵. Nitroreductase ablation, while effective, requires extended pro-drug exposure, leading to slower and less precise results ⁹⁹.

In contrast, TRPV1-mediated ablation is rapid and precise, with capsaicin triggering almost instantaneous cell death. This speed and simplicity make TRPV1 superior for experiments requiring quick, controlled ablation without the delays associated with other methods. Finally, chemogenetic approaches such as TRPV1/capsaicin permit prolonged experimentation in freely-moving animals, allowing us to collect large kinematic datasets necessary to rigorously study posture and locomotion.

Considerable progress has been made in recent years using new tools ^8,15–18and new perspectives ¹⁰⁰ to understand the cerebellar contributions to sensorimotor ^49–52and non-sensorimotor behaviors ^53–56in health and disease ^57–59. Underlying this considerable progress is an ever-improving ability to manipulate the cerebellum without compromising rigorous measures of behavior. Here — in support of similar goals — we validated a new chemogenetic approach (TRPV1/capsaicin-mediated activation and ablation) compatible with a high-throughput paradigm to measure behavior in freely swimming larval zebrafish (SAMPL). Our data uncover expected signatures of cerebellar contributions to posture and coordination, establishing the validity of our approach. Further, by comparing the impact of Purkinje cell ablation in time, we leverage the rapid maturation of the zebrafish to open a window into cerebellar control of posture and coordination across development. Our approach establishes a path forward for the larval zebrafish model to contribute to cerebellar mechanisms of postural control. The cerebellum emerged early in the evolution of vertebrates, when vertebrate life was underwater. Our work establishes a new tool to investigate ancient organizing principles of cerebellar function.

Materials and methods

Fish Care

All procedures involving zebrafish larvae (Danio rerio) were approved by the Institutional Animal Care and Use Committee of New York University. Fertilized eggs were collected and maintained at 28.5°C on a standard 14/10 hour light/dark cycle. Before 5 dpf, larvae were maintained at densities of 20-50 larvae per petri dish of 10 cm diameter, filled with 25-40 mL E3 with 0.5 ppm methylene blue. After 5 dpf, larvae were maintained at densities under 20 larvae per petri dish and fed cultured rotifers (Reed Mariculture) daily.

Fish Lines

To generate the Tg(aldoca:TRPV1-TagRFP) line, the 5-kbp aldolase Ca (aldoca) promoter ⁶⁴ and a gene cassette that includes TRPV1-Tag1RFP cDNA, rabbit beta-globin intron, and the SV40 polyadenylation signal (pAS) in pT2-4xUAS:TRPV1-RFPT ⁶¹ were subcloned into the Tol2 vector pT2KDest-RfaF ¹⁰¹ by the Gateway system (pT2K-aldoca-TRPV1-Tag1RFP-pAS). To establish stable transgenic lines, Tol2 plasmid and transposase mRNA (25 ng/µl each) were injected into one-cell-stage embryos.

The resulting Tg(aldoca:TRPV1-tagRFP) stable line allowed us to express the mammalian capsaicin-sensitive cation channel TRPV1 and the red fluorophore tagRFP in cerebellar Purkinje cells. Before exposure to capsaicin, fish were screened to ensure similar levels of tagRFP expression. When screening for transgene expression, we selected fish with clearly visible expression that was not excessively bright. The same criteria were applied when screening fish for GCaMP imaging and behavior experiments. Approximately a quarter of the fish that had aldoca:TRPV1-tagRFP expression had suitable expression levels for the activation experiment.

We measured neuronal activity using a genetically-encoded calcium indicator, Tg(UAS:GCaMP6s) ⁸⁶, driven by Tg(aldoca:GAL4) ¹⁰², or the Tg(elavl3:h2B-GCaMP6f) line ¹⁰³.

Confocal imaging of TRPV1-mediated activation / lesion

Images were collected using a Zeiss LSM800 confocal microscope using a 20x 1.0NA water immersion objective. Larvae were mounted in 2% low melting point agar (catalog #16520, Thermo Fisher Scientific) in a dorsal up position. Anatomical images were acquired from fish anesthetized with 0.2 mg/ml ethyl-3-aminobenzoic acid ethyl ester (MESAB, catalog # E10521, Sigma-Aldrich). To activate TRPV1-expressing Purkinje cells, fish were treated with 1 µM capsaicin in 0.2% DMSO in E3. To lesion Purkinje cells, fish were exposed to 10 µM capsaicin in 0.2% DMSO in E3. Control fish were treated with 0.2% DMSO in E3. Agar was removed around the tip of the tail to facilitate drug delivery. Fish were mounted throughout functional imaging experiments and kept in temperature controlled incubators between timepoints. Time series images were acquired with a 2 photon microscope (Thorlabs Bergamo equipped with a Mai Tai HP laser tuned to 920nm) with a framerate of 7.9fps. Images were analyzed in Fiji ¹⁰⁴; ROIs were drawn on nuclei of randomly selected Purkinje cells, which were then re-identified at each time point. Fluorescence for each cell and time point was normalized to the pre-capsaicin value. Cells with a dF/F0 of > 2 were considered activated and analyzed at the pre capsaicin and +6h timpoint. A Fisher’s exact test was performed to determine significance of the activated cells.

To image the anatomy of Purkinje cells exposed to 10 µM of capsaicin across time, the cerebellum was imaged at 7 dpf from fish mounted as above. Fish were unmounted and kept in E3 medium until the next day (8 dpf). At 8 dpf, fish were placed in 0.2% DMSO in E3 (control) or 10 µM capsaicin in 0.2% DMSO in E3 for 40-60min, and imaged again after 1h of recovery in E3 post-treatment. Fish from both groups were imaged again at 9 dpf. Confocal images were analyzed in Fiji and Purkinje cell somata were counted in both hemispheres of the cerebellum. A conservative approach was taken for cell counting, with inclusion of any structures still resembling cells, regardless of potential non-functionality or signs of degradation. Consequently, the counts are likely an underestimate of the actual percentage of cell loss.

Zebrafish behavior recordings

All behavior was measured using the Scalable Apparatus for Measuring Posture and Locomotion (SAMPL) apparatus, consisting of a chamber where larvae could swim freely, an infrared illuminator, a camera, and software to process video in real time. A comprehensive description of the apparatus is contained in ²⁶. Here we briefly describe the specific details of our experiments. Larvae were transferred to chambers at densities of 3-8 fish per chamber for 7 dpf experiments or 1-4 fish per chamber for 14 dpf experiments containing 25-30 ml of E3 or 0.2% DMSO / 1 µM capsaicin for activation experiments. After 24 h, behavior recordings were paused for 30-60 minutes for feeding (feeding pause) and 1-2 ml of rotifer culture was added to each chamber. Larvae were removed from the apparatus after 48h.

To monitor behavior before/during Purkinje cell activation, 7 dpf larvae were placed in chambers with E3. At 8 & 9 dpf, control fish were placed in 0.2% DMSO in E3 and the condition fish were placed in 1 µM capsaicin in 0.2% DMSO in E3 for 6h. The recording started about 10-15min after adding the fish to the capsaicin solution. Fish were fed after the 6h activation period. Video was sampled at 40Hz in constant darkness. Control: 9626 bouts (63% climb bouts)/149 fish/8 experimental repeats; Activation: 9664 bouts (61% climb bouts)/155 fish/8 experimental repeats;

To monitor behavior before/after Purkinje cell lesions, 7 dpf/14 dpf larvae were placed in the chambers with E3. After feeding at 8 dpf/15 dpf, fish were placed in petri dishes with 0.2% DMSO in E3 (control) or 10 µM capsaicin in 0.2% DMSO in E3 for 40-60min. Fish were then returned to the chambers in E3 and behavior recording was started. Video was sampled at 160Hz in constant darkness. 7 dpf lesions: Control: 17895 bouts (61% climb bouts)/110 fish/14 experimental repeats; Lesion: 17819 bouts (57% climb bouts)/120 fish/14 experimental repeats; 14 dpf lesion: Control: 10666 bouts (58% climb bouts)/48 fish/7 experimental repeats; Lesion: 10708 bouts (54% climb bouts)/44 fish/7 experimental repeats;

Pectoral fin amputations were performed at 13 dpf. Two length-matched siblings were anesthetized in 0.2 mg/ml ethyl-3-aminobenzoic acid ethyl ester (MESAB, catalog # E10521, Sigma-Aldrich) simultaneously and mounted in 2% low-melting temperature agar. Visualized under a stereomicroscope (Leica M80, 20x/12 eyepieces, 1.0x objective), the two pectoral fins from one larva were removed by pulling the base of the fin at the scapulocoracoid laterally with #5 Dumont forceps. After amputation, both fish were freed from the agar and allowed to recover in E3 until the next day, at which point half of the amputated and control fish were randomly selected for Purkinje cell lesions. Lesions were performed as above and behavior recorded for 48h. Behavior was recorded at a sampling rate of 160Hz with a 14/10h light-dark cycle. Control: 1506/5090/5353 (slow/medium/fast) bouts/15 fish/8 experimental repeats; Purkinje cell le- sion: 1667/6166/4299 (slow/medium/fast) bouts/18 fish/8 experimental repeats; Fin amputation: 1935/6295/4911 (slow/medium/fast) bouts/17 fish/8 experimental repeats.

Behavior analysis

Comprehensive descriptions of behavioral kinematics and baseline data for different genetic backgrounds are detailed in ²⁶. Here we describe the specific parameters used for our experiments. Behavior data were analyzed using custom-written software in MATLAB (Mathworks, Natick MA), which extracted individual swim bouts from the raw data (x/z position and pitch angle as a function of time). Only bouts during the circadian day were analyzed. Experimental repeats consisted of data collected across multiple SAMPL boxes from a single clutch of fish; the number of fish available determined how many apparatus were used (1-3). For comparisons across conditions (e.g. activation/control), fish from one clutch were randomly split into control and condition groups. As bout number is the fundamental unit of kinematic analysis, and different numbers of fish available would yield different numbers of bouts, we bounded our experiments to allow comparison across repeats. Specifically, if an experimental repeat contained less than 650 bouts it was excluded.

Between 22-27% of lesion experimental repeats contained less than 650 bouts and were not included in the analysis. For the activation experiments 56% (10 of 18) of experimental repeats were excluded with the 650 bouts threshold due to shorter recording times a higher fraction of experiments contained less than the threshold number of bouts. In subsequent analyses, the number of analyzed bouts was matched from both groups for a given experimental repeat to ensure an identical representation of control and condition bouts. Individual bouts were aligned at the time of peak speed. Bouts were excluded if their peak speed was <5mm/s or the fish rotated more than 30°(120°/sec) during the acceleration. The fractions excluded were as follows: for 7 dpf ablation: ctrl 0.2% lesion 0.15%; 7 dpf activation: ctrl 1% activation 1.7%; 14 dpf ablation dark: ctrl 0.05% ablation 0.05%; 14 dpf ablation light: ctrl 0.02% ablation 0.02%. For each experiment between 0.02% and 1.7% of bouts were excluded based on those criteria. Data was recorded either at 40Hz (activation experiments) or 160Hz (all other experiments). Effect size was calculated as the difference between the control value and the condition value relative to the control value. For fin body slope effect size the control value of the fast bin (i.e. largest slope) was used for effect size calculations to avoid overestimation of changes due to small control values.

Kinematic analyses proceeded as in ²⁶; key parameters were defined as follows:

• Posture is the pitch angle of the fish (long axis of the body relative to the horizon) at -250ms relative to peak speed, just before swim bout initiation. Positive values are nose-up.

• Climb Bouts are bouts with a trajectory of > 0°at the peak speed of the swim bout.

• Dive Bouts are bouts with a trajectory of < 0°at the peak speed of the swim bout.

• Upward rotation refers to the rotation from -250ms to the peak angular velocity; only bouts with positive upward rotation were included in the analysis of fin-body coordination.

• Lift is the residual change in depth (z) across a bout after subtracting the change expected from the posture of the fish as detailed in ²⁵. Briefly, the expected change is calculated using the distance the fish moves in x from -100 to 100ms and the pitch angle at -100ms. Only bouts with positive lift were included in the analysis of fin-body coordination.

• Fin-lift/rotation coordination is defined as the slope of the best linear fit between upward rotation and lift across bouts. The goodness of fit, R² was used as a measure of how well the fins and trunk are coordinated to generate lift, after ²⁵.

Functional GCaMP imaging in Purkinje cells

All calcium imaging experiments were performed using Tilt In Place Microscopy (TIPM), described comprehensively in ⁶⁸. Briefly, 7 dpf fish were mounted in the center of the uncoated side of a mirror galvanometer (catalog #GVS0111, Thorlabs) in 2% low-melting-point agarose. E3 was placed over the agarose, and the galvanometer mirror was placed under the microscope.

A microscope (Thorlabs Bergamo) was used to measure fluorescence elicited by multiphoton excitation (920nm) from a pulsed infrared laser (Mai Tai HP). Fast volumetric scanning was achieved using a piezo actuator (catalog #PFM450E, Thorlabs) to move the objective. Each frame of the volume (224 × 96 pixels) was collected with a 0.6 µs pixel dwell time (19.1 frames/s) resulting in a sampling rate of 3.82 volumes/s. While this imaging rate might be too slow to distinguish single spikes, it is suitable to measure a difference in calcium transients upon pitch stimulation to nose-up or nose-down direction. To set the galvanometer to a specific angle, a corresponding voltage was applied. The total angular range of the galvanometer is 40°. For the 30°stimuli, the galvanometer was driven to either +15°or -15°and then rotated so that the mirror was horizontal which allowed allowing for a 30°deflection in one direction. For each cell, 21 trials were initially recorded for one stimulus direction. The galvanometer was then remounted to allow for a 30°stimulus in the opposite direction, and 21 trials were similarly recorded for this direction. The order of nose-up and nose-down blocks were alternated for different fish. After all 42 trials were recorded fish were anesthetized with 0.2 mg/ml MESAB; after 10min the baseline fluorescence at ±30°was recorded to establish a baseline that controlled for eccentricity. Analysis was done using Fiji and MATLAB. In total 43 Purkinje cells were imaged and 31 cells were kept from 8 fish. Only Purkinje cells that could be reliably identified at ±30°were analyzed. To map the anatomical locations of the recorded cells, we imaged overview stacks for each fish. These stacks were manually aligned in Illustrator, and the cells included in the analysis were identified and color-coded according to their tuning properties.

Regions of interest were drawn in Fiji and loaded into MATLAB to extract the intensity of fluorescence after motion correction was performed ¹⁰⁵. The integral of each stimulus was calculated and trials of the same direction were averaged as the tonic response to ±30°pitch. To extract cells with directional information the directionality index (DI) was calculated by dividing the difference of the up and down responses by the sum of it. Cells with a DI greater than ±0.35 were considered directionally tuned. Only Purkinje cells that were not directionally tuned were used for principal component analysis and subsequent support vector machine decoding analysis.

To classify trial identity in the dataset, we used a support vector machine (SVM) with a linear kernel. The SVM model was trained using k-fold cross-validation, which splits the data into k subsets (folds). At each iteration, the model was trained on k-1 folds and tested on the remaining fold, ensuring that the model performance was evaluated on unseen data in each fold. Permutations were performed on randomized trial identity as a null hypothesis (5-fold cross-validation; 100 shuffles for randomization). Accuracy was calculated as 1 minus the classification loss.

For calcium imaging in 7 and 14 dpf larvae, a horizontal imaging protocol was used. In total 11 fish were imaged at 7 dpf and 7 fish at 14 dpf. A total of 138/90 (7/14 dpf) cells were recorded. Cells were imaged while the fish was horizontal. For horizontal imaging, we used a ±19°stimulus, enabling us to alternate between up and down trials without the need to remount the galvanometer. Before each trial, a 15 sec period was recorded; the average activity during this time was used as the baseline. Fish were pitched nose-down (-19°) for 15s and rapidly returned to horizontal, where-upon calcium activity was measured. This stimulus was then repeated in the nose-up (+19°) direction. The maximum dFF of the first second upon return was analyzed. Cells were classified into directional or non-directional based on the directionality index as described above. PCA and decoder analyses were performed using activity from non-directional cells. Decoding accuracy was tested for each fish individually.

Statistics

All statistical testing was done in Matlab R2020a. Data across repeats was pooled for analysis. To assess the variability and determine whether pooling individual experimental repeats within each group was appropriate, we performed a two-way analysis of variance (ANOVA) on the interquartile ranges (IQRs) of the single experimental repeats for the 7 days post-fertilization (dpf) activation, the 7dpf lesion, and the 14dpf lesion experiments without excluding experimental repeats. The results of the ANOVAs and the IQRs for all experimental repeats are reported in Tables 6 to 11.

Table 5.

Table 6.

Table 7.

Table 8.

Table 9.

Table 10.

Table 11.

To estimate the spread of the data we resampled distributions 100 times with replacement from the data from each condition and computed the expected value for control and perturbed datasets. These permutations were then used to explicitly compute a p-value for fitted variables (slope and R² of fin body coordination).

For other variables two-sided Wilcoxon rank sum tests were performed. To correct for multiple testing the critical p-value was calculated based on α = 0.05 using Šidáks method. The critical p-value for each data set is reported in the respective table. Outliers were determined as deviating more than 3 times the scaled median absolute deviation (MAD) from the median. A scaling factor of 1.4826 was used to ensure that MAD-based outlier detection is consistent with other methods like Z-scores. Data is shown as median and 95% confidence interval of the median for measured parameters or as median with 25^th and 75^th percentile for bootstrapped variables. The 95% confidence intervals of the median were bootstrapped using 1000 samples. The medians with 95% confidence intervals for all parameters are reported in the tables. For linear fits a robust regression model (bisquare) was used and fitted variables (slope and R² of fin body coordination) were bootstrapped (using 100 samples). To test speed dependency of the fin-lift/rotation ratio Spearman’s rank correlation was computed and control and lesion values were compared. First, the correlation coefficients where transformed using Fisher’s z-transformation to enable direct comparison of the z-scores. The difference between the z-scores was divided by the standard error and a z-test was performed. Additionally, we only considered effect sizes of ≥15% to be biologically relevant.

Data & Code

All data, raw and analyzed, as well as code necessary to generate the figures is available at the following DOI: 10.17605/OSF.IO/9X57Z

Acknowledgements

Tg(aldoca:TRPV1-tagRFP) fish were generated using a plasmid that was a gift from David Prober’s laboratory. Research was supported by the National Institute on Deafness and Communication Disorders of the National Institutes of Health under award number R01DC017489. The authors would like to thank Martha Bagnall and Abigail Person along with the members of the Schoppik and Nagel lab for their valuable feedback and discussions.

Additional information

Author contributions

Conceptualization: FA and DS, Methodology: FA and DS, Investigation: FA and KN, Visualization: FA, Writing: FA, DS Editing: DS, Funding Acquisition: FA and DS, Supervision: DS, MH. Tg(aldoca:TRPV1-tagRFP) fish were generated by KM.