Introduction

Developing sensorimotor reflex circuits must precisely connect functional subtypes of neurons to ensure appropriate behavior. For example, withdrawal from noxious stimuli requires maturation of a sensorimotor circuit that uses subtypes of spinal interneurons to transform noxious stimulation into activation of both ipsilateral flexor and contralateral extensor motor neurons 1. Work over the past 40 years has highlighted motor partner populations as possible orchestrators of connectivity in pre-motor reflex circuits 26, but controversy remains about the nature of their role. In the spinal cord, molecular perturbations of motor neuron identity have provided evidence both for 711 and against 1216 an instructive role in establishing connectivity. Part of this controversy stems from the wide variety of inputs to spinal motor neurons 17, the molecular and functional heterogeneity of pre-motor interneurons 14,15, and their complex roles in gait and posture 18. Further, transcription factors play multivariate and redundant roles in spinal motor neuron development 19,20, such that the effects of molecular perturbations of identity can be masked.

The sensorimotor circuit for vertical gaze stabilization offers a simple framework to evaluate whether and how motor neurons shape pre-motor circuit fate and connectivity. The vertebrate vestibulo-ocular reflex circuit consists of three neuron types – peripheral sensory, central projection, and extraocular motor neurons – that stabilize gaze after head/body tilts (Figure 1A) 21. Subtype fate, anatomical connectivity, and function are inextricably linked: directionally-tuned sensory neurons innervate nose-up/nose-down subtypes of projection neurons, which in turn innervate specific motor neurons that selectively control either eyes-down or eyes-up muscles 3,2229. As both the recipients and origin of directional information, projection neuron fate specification is paramount to proper circuit assembly. Recent work has established the vertical vestibulo-ocular reflex circuit in zebrafish as a model to uncover determinants of fate and connectivity 2426,29 given the ease of optical imaging, abundant tools for genetic perturbations, rapid development, and robust evolutionary conservation.

phox2a loss-of-function mutants fail to develop nIII/nIV motor neurons and vertical eye rotation behavior.

Associated with Figure S2.

(A) Schematic of vestibulo-ocular reflex circuitry and the genetic loss-of-function approach used to perturb motor-derived signals.

(B) Fluorescent in situ hybridization showing phox2a transcript expression in statoacoustic ganglion sensory afferents (left), central projection neurons in the tangential nucleus (middle), or nIII/nIV extraocular motor neurons (right) at 5 days post-fertilization (dpf). Top: probe only, nuclei outlined with dashed lines. Bottom: probe (green) merged with somata, labeled by Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) (sensory, central) or Tg(isl1:GFP) (motor).

(C) Schematic of CRISPR/Cas9 mutagenesis approach. Top: Red star shows location of guides against phox2a DNA. Bottom: RNA sequence in wildtype and phox2ad22 alleles. Red dashed lines show deleted sequence; “STOP” box shows predicted premature stop codon due to deletion. Right shows predicted protein sequence.

(D) Transmitted light image of a 5 dpf wildtype (top) and phox2a null mutant (bottom). White arrows point to a normally inflated (top) or absent (bottom) swim bladder.

(E) Images of nIII/nIV motor neurons in one hemisphere, labeled by Tg(isl1:GFP), in wildtype siblings (left) and phox2a null mutants (right) at 5 dpf. Scale bar, 20 µm.

(F) Quantification of the number of Tg(isl1:GFP)+ neurons in nIII/nIV from N=6 wildtype siblings and N=10 phox2a null mutants.

The current model for vestibulo-ocular reflex circuit development was motivated by pioneering work in chick 3,30 and formalized by Hans Straka: “[circuit assembly] is accomplished by a specification process that retrogradely transmits post-synaptic target identities to pre-synaptic neurons.” 31. In its strongest form, this “retrograde” model posits a causal role for extraocular motor neurons in specifying the fate (sensory selectivity) of central projection neurons. This key prediction – that loss of motor neurons would disrupt sensory selectivity in pre-synaptic projection neurons – remains untested. In zebrafish, extraocular motor neurons are temporally poised for such a role. Motor neurons are organized into spatial pools, and though synaptogenesis at ocular muscle targets begins late in development 32, motor neuron fate (muscle target and pool location) is determined early 33. Projection neurons are born at roughly the same time as motor neurons and extend axons shortly afterwards, poising them to receive deterministic signals that could retrogradely specify their sensory selectivity 29.

Here, we adopted a loss-of-function approach to determine whether motor partner populations specify identity or instruct connectivity across an entire vestibular reflex circuit in zebrafish. We generated a new mutant allele for the phox2a gene to eliminate the extraocular motor neurons used for vertical gaze stabilization. Combining functional, anatomical, and sequencing approaches, we then demonstrated that motor neurons are dispensable for three aspects of pre-motor reflex circuit assembly: (i) directionally-appropriate connectivity between sensory and projection neurons, (ii) assembly of projection neurons with motor partners, and (iii) the transcriptional profiles of projection neurons. The current model of vestibulo-ocular reflex circuit development must there-fore be revised: up/down projection neuron subtype fate cannot be retrogradely established by a motor partner-derived signal. Instead, the signals that specify fate must lie elsewhere. More broadly, our work argues against a deterministic role of motor neurons in specifying the fate and sensory connectivity of pre-motor circuit components.

Results

Constitutive loss of phox2a prevents extraocular motor neuron specification and impairs vertical gaze stabilization behavior

Extraocular motor neurons for vertical/torsional gaze stabilization are located in cranial nuclei III (nIII) and IV (nIV). To eliminate nIII/nIV motor neurons and by extension, any secreted signals, we used a genetic loss-of-function approach (Figure 1A). A single highly-conserved transcription factor, phoxa, specifies nIII/nIV fate 3437. In the vestibulo-ocular reflex circuit, phox2a is exclusively expressed in nIII/nIV motor neurons but not its upstream partners (Figure 1B). Therefore, phox2a is an ideal genetic target to eliminate motor-derived signals without compromising evaluations of upstream functional development.

Prior mutagenesis established a phox2a loss-of-function allele in zebrafish 34, but the line has since been lost. Here, we generated three new phox2a loss-of-function alleles using CRISPR/Cas9 mutagenesis (Figure 1C) (one allele shown here; additional alleles described in Methods). Consistent with prior reports and human mutations 38, both eyes in phox2a null mutants were exotropic (rotated towards the ears) reflecting a loss of motor neurons in nIII/nIV. phox2a mutants failed to hatch from their chorions without manual intervention and did not inflate their swim bladders by 5 days post-fertilization (dpf) (Figure 1D), phenotypes not previously reported 34. Consequently, null mutants do not survive past 7 dpf. We did not observe these morphological phenotypes in wildtype and heterozygous siblings (Figure 1D). As vestibulo-ocular reflex circuit architecture and behavior is established by 5 dpf 24,25,29, premature lethality did not preclude further measurements of circuit development.

To validate phox2a loss-of-function, we leveraged a downstream transcription factor: isl1 39. The Tg(isl1:GFP) line 40 labels all nIII/nIV motor neurons except inferior oblique neurons 33, which comprise one of four pools for upwards eye rotations. We first quantified changes in the number of labeled nIII/nIV neurons (Figure 1E-Figure 1F). In phox2a mutants, we observed an expected and near-total loss of isl1 expression (WT: 298±19 neurons across both hemispheres; null: 19±11 neurons; Wilcoxon rank sum test, p=2.5×10−4) at 5 dpf, well-after nIII/nIV differentiation is complete 33. As a secondary measure of nIII/nIV motor neuron fate, we performed fluorescent in situ hybridization against vachta, a marker of cholinergic motor neurons, when differentiation is complete at 2 dpf. Unlike sibling controls, nIII/nIV neurons in phox2a mutants failed to express vachta, validating that loss of phox2a perturbs at least two key transcriptional markers of nIII/nIV fate (Figure S1). Unexpectedly, we also observed slightly fewer neurons in phox2a heterozygotes (heterozygote: 229±20 neurons; Wilcoxon rank sum test against WT, p=6.7×10−4). In heterozygotes, loss of isl1 fluorescence was restricted to the medial domain of dorsal nIII, which contains some of the earliest-born neurons in nIII/nIV (Figure S2A-Figure S2C) 33. Globally, this manifested as a rostral and ventral shift in the positions of all neurons mapped (Figure S2D) (two-sample, two-tailed KS test, WT vs. heterozygotes: mediolateral axis, p=0.13; rostrocaudal: p=4.0×10−29; dorsoventral: p=2.5×10−9). This region contains two motor pools that control the inferior (IR) and medial rectus (MR) muscles 33. We conclude that phox2a acts in a dose- and birthdate-dependent manner to specify nIII motor pool fate.

Together, these observations validate our phox2a loss of function alleles as a selective means to disrupt nIII/nIV motor neuron fate specification and vertical eye rotation behavior.

Peripheral-to-central circuit assembly does not require motor partners

Vertical gaze stabilization requires that (1) peripheral VIIIth nerve sensory afferents relay tilt sensation (nose-up/nose-down) directly to projection neurons in the tangential nucleus, and (2) projection neurons innervate appropriate nIII/nIV pools (eyes-up/eyes-down). For proper circuit function, appropriate connectivity must first develop across up/down circuit subtypes. The “retrograde” model predicts that motor partners specify circuit assembly. Therefore, in the absence of motor neurons, projection neurons should fail to respond selectively to directional tilt sensations – either due to loss of their fate, the fate of upstream sensory afferents, or sensory-to-central connectivity.

To evaluate upstream circuit formation, we measured tilt-evoked responses in projection neurons using Tilt-In-Place Microscopy (TIPM) 41 (Figure 2A-Figure 2B). Peripherally, tilts activate utricular VIIIth nerve sensory inputs to projection neurons 29,41. We used a galvanometer to deliver tonic nose-up and nose-down pitch tilts to phox2a null larvae and sibling controls. We then measured the activity of a calcium indicator, GCaMP6s 42, in projection neurons. We performed experiments at 5 dpf, when nearly all projection neurons are selective for one tilt direction 29, circuit architecture is stable 25, and gaze stabilization behavior is directionally-appropriate 24.

Motor neurons are dispensable for proper connectivity between utricular sensory afferents and projection neurons.

Associated with Table 1.

(A) Schematic of pitch vestibulo-ocular reflex circuitry. Dashed lines outline projection neurons as calcium imaging target. Nose-down/eyes-up channel represented with blue; orange, nose-up/eyes-down.

(B) Schematic of tonic pitch-tilt stimulus delivered with Tilt-In-Place Microscopy (TIPM). Shaded regions show calcium imaging windows when fish were oriented horizontally immediately following tilts. Inset shows timecourse of the rapid step to restore horizontal position after tilts. Imaging experiments used larvae from the Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:GCaMP6s) line.

(C) Proportion of subtypes observed in sibling controls and phox2a null mutants. Blue: nose-down. Orange: nose-up. Grey: Neurons without directional tuning (criteria in Methods).

(D/E) Soma position of nose-down (blue) and nose-up (orange) neurons in sibling controls (left) and phox2a null mutants (right). Soma size scaled by strength of directional selectivity (min=0; max=1; see Methods).

(F/I) Heatmaps showing example tilt responses from nose-down (F) or nose-up (I) neurons in sibling controls (top) and phox2a null mutants (bottom). n=10 neurons with strongest ΔF/F responses to tilts shown. Each row shows an individual neuron. Shaded bars show calcium imaging window immediately following restoration from eccentric position. Black arrow points to first second of tilt response used for analyses.

(G/J) Distributions of maximum ΔF/F responses to tilts for nose-down (G) or nose-up (J) neurons in sibling controls (black) and phox2a null mutants (red). Solid and shaded lines show mean and standard deviation, respectively, of bootstrapped data (Methods)

(H/K) Distributions of directional tuning score to tilts for nose-down (H) or nose-up (K) neurons in sibling controls (black) and phox2a null mutants (red). Tuning score ranges from 0 (equal responses to both tilt directions, no tuning) to 1 (responses to one tilt direction only); criteria detailed in Methods. Solid and shaded lines show mean and standard deviation, respectively, of bootstrapped data.

Projection neuron responses and topography were strikingly unchanged in phox2a mutants compared to controls. There was no change in the number of projection neurons observed across N=5 sibling controls and phox2a mutants (mean control number: 42±4 neurons per hemisphere; mean in phox2a mutants: 39±3 neurons; two-tailed Wilcoxon rank sum test, p=0.22). We next recorded the activity of n=297 neurons from N=16 phox2a mutants and n=440 neurons from N=21 sibling controls (Methods and Table 1 split by genotype). We observed comparable ratios of projection neuron subtypes (sib: 46% nose-down, 47% nose-up, 7% untuned; phox2a: 49% nose-down, 44% nose-up, 7% untuned) (Figure 2C). Next, we evaluated their topography (Figure 2D-Figure 2E). Projection neurons are topographically organized along the dorso-ventral axis by their directional selectivity 29. Global spatial separation between subtypes remained significant in phox2a mutants (one-way multivariate analysis of variance, p=0.004). We also compared the topography of nose-up and nose-down neurons separately across phox2a genotypes (Methods and Ta-ble 2 split by genotype). Nose-down neurons were comparably distributed between null and control larvae (one-way multivariate analysis of variance, p=0.15). We observed a minor lateral shift to nose-up neurons in null mutants (median mediolateral position, sib: 15.2µm from medial edge; phox2a, 13.2µm; two-tailed, two-sample KS test, p=3.0×10−4) but no changes in other spatial axes (dorsoventral: p=0.16; rostrocaudal: p=0.56). The small medial deviation (2µm across a 40µm space) is within the limits of our registration error. We conclude that projection neuron topography is established independently of motor partners.

Statistical comparisons of tilt responses across genotypes.

WT (sampled) refers to an n=125 neuron subset, sampled with replacement from a reference dataset of wildtype projection neurons. Data shown is mean/standard deviation unless otherwise noted. p val generated from a one-way ANOVA with multiple comparisons. Associated with Figure 2 and Figure 3.

Projection neuron sensitivity and selectivity also developed comparably between phox2a mutants and siblings (Figure 2F-Figure 2K). Projection neurons responded to tilt sensations with comparable magnitudes (Figure 2G-Figure 2J) (nose-down mean ΔF/F, sib: 1.86±1.69; phox2a: 2.07±1.48; two-tailed Wilcoxon rank sum test, p=0.98; nose-up mean ΔF/F, sib: 1.24±1.23; phox2a: 1.02±0.89; p=0.18). Previously, we defined a metric to describe a neuron’s selectivity for one tilt direction (0 = equal responsesto up/down; 1 = maximally selective) 29. Directional selectivity remained unchanged in phox2a mutants (Figure 2H-Figure 2K) (nose-down mean index, sib: 0.73±0.29; phox2a: 0.68±0.29; two-tailed Wilcoxon rank sum test, p=0.85; nose-up mean index, sib: 0.85±0.26; phox2a: 0.81±0.29; p=0.12). Collectively, this demonstrates that the functional responses of projection neurons and, by inference, connectivity with utricular afferents are not shaped by motor partners.

Ventral projection neurons receive additional input from the semicircular canals 29, which encode phasic (fast) tilt sensation. To activate sensory afferents from the semicircular canals, we used TIPM to deliver two impulses of angular rotation (Figure 3A-Figure 3B) 29,41. We observed no changes in phox2a mutants (Table 1 and Table 2 for genotype splits). Projection neurons responded to impulses in comparable ratios (Figure 3E) (sib: 58% responsive; phox2a: 71% responsive). Responsive projection neu-rons remained localized to the ventral nucleus (dorsoventral axis: two-tailed, two-sample KS test, p=0.99). We observed minor positional changes in the rostrocaudal and mediolateral axes (Table 2), which are not topographically organized by impulse responsivity 29; this deviation is again within our margin of registration error. Lastly, the functional properties of projection neurons were unchanged (Figure 3F-Figure 3G). We observed no change in calcium response magnitudes (Figure 3F) (mean ΔF/F, sib: 0.33±0.29; phox2a: 0.36±0.40; two-tailed Wilcoxon rank sum test, p=0.85) or lack of directional selectivity (Figure 3G) (mean index, sib: 0.07±0.41; phox2a: 0.06±0.40; p=0.39). Therefore, fate and connectivity between phasic sensory afferents and projection neurons must not require motor partners.

Statistical comparisons of projection neuron topography across genotypes.

WT (sampled) refers to an n=125 neuron subset, sampled with replacement from a reference dataset of wildtype projection neurons. Data shown is the median/standard deviation distance from the ventro-lateral and rostral edges of the tangential nucleus (total size: 40 µm across each spatial axis). p val from one-way analysis of variance (individual spatial axes) or multivariate analysis of variance (global organization), respectively.

Motor neurons are dispensable for proper connectivity between semicircular canal sensory afferents and projection neurons.

Associated with Table 1.

(A) Schematic of impulse tilt experiments. Yellow dashed lines outline projection neurons as calcium imaging target. Impulse-responsive neurons (ventrally-localized) shown with purple; unresponsive neurons, grey.

(B) Schematic of impulse stimuli delivered with TIPM. Shaded regions show calcium imaging windows at horizontal immediately following impulses. Inset shows timecourse of impulse stimulus. Imaging experiments used larvae from the Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:GCaMP6s) line.

(C) Proportion of impulse-responsive (purple) and unresponsive (grey) neurons observed in sibling controls and phox2a null mutants.

(D) Soma position of impulse-responsive neurons in sibling controls (left) and phox2a null mutants (right). Soma size scaled by strength of calcium response (ΔF/F), normalized by max observed ΔF/F.

(E) Heatmaps showing example impulse responses from neurons in sibling controls (left) and phox2a null mutants (right). n=10 example neurons shown. Each row shows an individual neuron. Shaded bars show calcium imaging window immediately following impulse delivery. Black arrow points to first second of tilt response used for analyses. Note smaller scale (0-0.75) of impulse responses relative to Figures 2F and 2I.

(F) Distributions of maximum ΔF/F responses to impulses in sibling controls (black) and phox2a null mutants (red). Solid and shaded lines show mean and standard deviation, respectively, from bootstrapped data.

(G) Distributions of directional tuning score to impulses in sibling controls (black) and phox2a null mutants (red). Tuning score ranges from 0 (equal responses to both tilt directions, no tuning) to 1 (responses to one tilt direction only); criteria detailed in Methods. Solid and shaded lines show mean and standard deviation, respectively, from bootstrapped data.

Lastly, we considered whether loss of one subtype of nIII/nIV neurons might alter connectivity. For example, loss of eyes-down motor pools could impair wiring between their corresponding nose-up sensory and projection neuron partners. Here, we leveraged phox2a heterozygotes, which lack a subtype of nIII neurons (IR/MR) that contribute to downwards eye rotations (Figure S2). We observed no differences in tonic tilt responses between phox2a wildtype, heterozygote, and null larvae, though we did note a minor decrease in response strength to impulses (statistics in Table 1). We note that phox2a heterozygotes do not lack all motor pools for downwards eye rotations. Nevertheless, we conclude that individual motor pools do not meaningfully contribute to connectivity between sensory and projection neurons.

Taken together, these experiments demonstrate intact directional selectivity for two peripheral sensory inputs – utricular and semi-circular canal VIIIth nerve afferents – and appropriate connectivity with projection neurons. We conclude that functional sensory-to-central circuit formation is established independently of motor partners.

Projection neurons remain competent to assemble with appropriate motor targets

Motor partners could secrete signals that initiate pre-motor axon outgrowth, target arriving axons to specific motor pools, or trigger synaptogenesis 3. Motor pool topography in nIII/nIV reflects ocular muscle targets: dorsal pools innervate downward-rotating muscles (superior oblique and inferior rectus), while ventral pools target the converse (eyes-up, superior rectus and inferior oblique) 33,43. In turn, projection neuron somatic and axonal organization mirrors motor pool topography 28,29, which could facilitate directionally-selective circuit assembly. We reasoned that projection neurons may fail to initiate axon outgrowth, target spatially-appropriate motor pools, and/or form synapses in phox2a mutants. To test this hypothesis, we measured changes in projection neuron anatomy at 5 dpf, when axonal arbors are established and stable 25.

To test whether projection neurons establish gross, long-range (hindbrain to midbrain) axonal outgrowths, we performed optical retrograde labeling 44 using a photolabile protein, Kaede. We targeted the medial longitudinal fasciculus at the midbrain-hindbrain boundary, which contains projection neuron axons 25,29 (Figure 4A). In both phox2a mutants and sibling controls, we observed retrograde photolabeling of projection neuron soma (Figure 4B), supporting that initial axon outgrowth does not require motor partner-derived signals.

Projection neurons are anatomically and molecularly poised to assemble with motor neuron partners in phox2a mutants.

(A) Schematic of retrograde photofill experiments. Projection neuron axons expressing the photolabile protein Kaede are targeted at the midbrain-hindbrain boundary with ultraviolet light. Converted protein (magenta) retrogradely diffuses to the soma, while the unconverted nucleus remains green.

(B) Projection neuron somata in sibling controls (left) and phox2a null mutants (right) after retrograde photolabeling. Experiments performed at 5 dpf. Neurons visualized in Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:E1b-Kaede).

(C) Top two panels: Projection neuron axons at the hindbrain (inset) and midbrain-hindbrain boundary in sibling controls (top) and phox2a null mutants (bottom). Axons visualized using Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:E1b-Kaede). White dashed outline (circle) shows arborization fields in nIII/nIV. Dashed box over axons shows location of two bottom panels. MHB and yellow dashed line, midbrain-hindbrain boundary. nucMLF: nucleus of the longitudinal fasciculus. Bottom two panels: Zoom of axons (dashed rectangle above). Spatial segregation between early-born (magenta+green) and late-born (green only) axons. White dashed line reflects separation between dorsal (nose-up, early-born) and ventral (nose-down, late-born) axon bundles. Image at 5 dpf in sagittal view.

(D) Projection neuron axon bundle in a phox2a null mutant at 3 dpf. White arrows point to single collateral to two nIII/nIV neurons that were not eliminated following phox2a knockout.

(E) Fluorescent in situ hybridization against RNA for three pre-synaptic markers: synaptophysin a (sypa; left), synaptic vesicle glycoprotein 2 (sv2, middle), and synapsin I (syn1, right). Top row, sibling controls. Bottom row, phox2a null mutants. For each panel set, left images show in situ probe expression (green) and right images show merge with projection neurons labeled in Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:E1b-Kaede). Dashed lines outline the projection nucleus. Cell and transcript expression outside the projection nucleus is removed for visual clarity. Images taken at 5 dpf in sagittal mount. All scale bars, 20 µm.

Next, we evaluated whether projection neuron axons remain capable of wiring with spatially-appropriate motor partners. Projection neuron axons segregate along the dorsal (nose-up) and ventral (nose-down) axes according to their birth order (early/late born, respectively) 28,29 and the pool topography of their motor targets 33. To test whether projection neurons retain this topography, we optically labeled the axons of early-born (before 30 hpf) projection neurons 29. In phox2a mutants, axons remained dorsoventrally segregated at midbrain targets (Figure 4C, inset). Typically, projection neurons robustly collateralize to nIII/nIV targets at the midbrain-hindbrain boundary. We did not observe collaterals to nIII/nIV in phox2a mutants (Figure 4C). However, projection neurons still robustly arborized to more rostral, spinal-projecting targets in the nucleus of the medial longitudinal fasciculus, suggesting they retain the machinery necessary to collateralize. Consistent with this hypothesis, we observed that projection neurons formed occasional, small collaterals in phox2a mutants with few (1-5%) nIII/nIV neurons remaining (Figure 4D). We conclude that projection neurons remain competent to assemble with spatially-appropriate targets.

If motor neurons are required to initiate synaptogenesis, then projection neurons should fail to develop pre-synaptic machinery. To test this hypothesis, we performed fluorescent in situ hybridization against common pre-synaptic transcripts: synaptophysin a (sypa), synaptic vesicle glycoprotein (sv2), and synapsin I (syn1). In both phox2a mutants and controls, we observed robust transcript expression in projection neuron somata at 5 dpf (Figure 4E), well-after synaptogenesis onset in wildtype larvae 29. Motor partner-derived signals are thus not required for projection neurons to develop the necessary components for synaptogenesis.

Though motor neurons may play later roles in selecting and/or refining pre-motor input specificity, our data supports that projection neurons remain anatomically and molecularly poised to assemble with appropriate targets. We predict that absent collaterals and synapses reflect a lack of adhesive contact necessary to stabilize 45,46, but not instruct the formation of nascent structures.

The transcriptional profiles of projection neurons are intact in the absence of motor partners

We found that functional and anatomical connectivity, from peripheral sensors to motor targets, develop independently of motor partners. Fate in the vestibulo-ocular reflex circuit follows from connectivity 21, but neuronal fate can also be defined with respect to unique transcriptional signatures. Previously, we developed a sequencing pipeline to discover transcription factors that specify functional subtypes of spinal motor neurons and evaluate the consequences of perturbations on transcriptional fate 47. We adapted this approach to determine if loss of motor-derived signals changed the transcriptional profiles of projection neurons.

We compared the transcriptional profiles of projection neurons in phox2a mutants and sibling controls (Figure 5A) using bulk RNA sequencing. Sibling controls included both wildtype and heterozygous phox2a larvae (Figure S2). We performed sequencing experiments at 72 hours post-fertilization (hpf), after projection neuron differentiation is complete and synaptogenesis to motor targets has peaked 29. We sequenced projection neurons labeled by Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) 24,25,48,49 (Methods, Figure 5B, Figure S3). Neurons labeled in this line include, but are not exclusive to the projection neurons in the tangential nucleus used for vertical gaze stabilization. Therefore, we evaluated our bulk RNA sequencing dataset in the context of a single-cell reference atlas derived from the same transgenic line (Methods, Figure S4) to minimize noise from other labeled populations. We used evx2 50 as a reference, as it was expressed in all projection neurons (Figure S4D) and highly detected (36%) in singly-profiled projection neurons.

Motor neurons are dispensable for normal transcriptional profiles of projection neurons.

Associated with Figure S3-Figure S6, Table 3.

(A) Schematic of sequencing approach. Central projection neurons (Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:E1b-Kaede)) are harvested from 3 dpf larvae. Flow cytometry is used to exclude neurons not labelled by Tg(−6.7Tru.Hcrtr2:GAL4-VP16). Bulk RNA sequencing is performed to compare the profiles of projection neurons in siblings and phox2a null mutants.

(B) Example of projection neurons before (left) and after (right) harvesting. Neurons visualized with Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:E1b-Kaede). Dashed lines outline projection neurons in the tangential nucleus; dotted lines, medial vestibular nucleus. Yellow region shows margin of harvesting error: non-projection neurons that may be included in bulk sequencing dataset.

(C) Number of differentially expressed genes in projection neurons at 3 dpf after applying progressive filters based on gene expression in a reference single-cell dataset. Data shown on logarithmic scale. Solid, dashed, and dotted lines represent differentially-expressed gene with p adjusted<0.5, p adjusted<0.01, or p adjusted<0.001 significance, respectively.

(D) Volcano plot showing differentially expressed genes in projection neurons between control and phox2a null larvae at 3 dpf. Dashed lines represent significance cutoffs: horizontal line, p adjusted>0.05; vertical line, Log2 Fold Change > 2.0. Each circle is a gene. Genes to the left and right of 0 on the horizontal axis show downregulated and upregulated genes, respectively. Colors indicate percent of reference cells that express a given gene. Grey-colored genes are below both significance thresholds.

(E) Same data as Figure 5D. Colored genes show eight candidates evaluated with fluorescent in situ hybridization: red, upregulated; blue, downregulated; yellow, highly-expressed controls (evx2.

(F) Fluorescent in situ hybridization against candidate genes that met projection neuron filter criteria. Top row shows sibling controls; bottom row, phox2a null mutants. For each gene, left panels show RNA probe (green) and right panels show merge with projection neurons labeled by Tg(−6.7Tru.Hcrtr2:GAL4-VP16) (grey). Dashed lines outline the projection nucleus. Cell and transcript expression outside the projection nucleus is masked for visual clarity. Arrows denote whether genes are upregulated (red), downregulated (blue), or not significantly changed (yellow). Percentage refers to fraction of cells in a single-cell RNA sequencing reference atlas (Methods) with detected transcript. Candidates: itga9 (log2 fold change=23.0, p adj.=3.9×10−6), twf1b (log2 fold change=5.9, p adj.=0.024, p4hb (log2 fold change=5.1, p adj.=0.04), mapk6 (log2 fold change=5.1, p adj.=0.06), rxfp2a (log2 fold change=-8.5, p adj.=1.1×10−5), satb1a ((log2 fold change=-3.0, p adj.=0.001), evx2 (log2 fold change=0.46, p adj.=0.99), myt1la ((log2 fold change=2.6, p adj.=0.44). All scale bars, 20 µm.

There were strikingly few differentially-expressed genes in projection neurons between phox2a siblings and null mutants (Figure 5C, Table 3). All candidate differentially-expressed genes were lowly-expressed (detected in <10% of reference projection neurons, Figure 5D). To determine if any candidates were differentially expressed in projection neurons, we used a fluorescent in situ hybridization method 51 in which fluorescence intensity correlates with detected transcript expression reliably across individual larvae (Figure S5). We evaluated 8 candidate genes (Figure 5E-Figure 5F); four with the highest detection levels in projection neurons (satb1a, rxfp2a, mapk6, and p4hb), two with high significance and fold change thresholds (itga9, twf1b), and two highly-detected controls (evx2, myt1la). Qualitatively, we observed no differences in expression patterns between phox2a mutants and siblings.

Differentially expressed genes in projection neurons.

Star indicates a gene was evaluated using fluorescent in situ hybridization. # symbol indicates a gene was also differentially expressed in adjacent phox2a-expressing medial vestibular neurons (see Figure S7). “% of projection neurons with expression” refers to detection in a filtered subset of projection neurons from a single-cell reference atlas of neurons labeled in Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) (Methods, Figure S4). Genes sorted by p adjusted value. Data associated with Figure 5.

We considered that our inability to detect differentially-expressed genes could arise from our exclusion of candidates based on their expression in our reference single-cell atlas. Therefore, we repeated our analyses in unfiltered bulk sequencing data. The top 50 highest-expressed genes in phox2a siblings were highly detected in singly-profiled neurons labeled by Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede), including projection neurons (Table 4) (mean detection in projection neurons: 55%±21%, min=0%, max=87%). This suggests that our dissections adequately captured our target population. However, we again identified few differentially-expressed genes in our unfiltered data (Figure S6A-Figure S6C), with substantial decreases as significance stringency increased. In situ hybridization validated that top candidates remained lowly expressed in projection neurons in both phox2a siblings and mutants (Figure S6D), Importantly, nearly all candidates had low detection across all neurons in our reference single-cell atlas and had predicted expression in populations such as glia and the caudal hindbrain (Methods, Table 5). Notably, some candidates were highly expressed in a subset of medial vestibular nucleus neurons, which lies on the medial edge of the tangential nucleus and expressed phox2a (Methods, Figure S7, and Table 6). Together, we conclude that any differential gene expression in our data either reflects noise or contamination from other labeled populations, but not projection neurons in the tangential nucleus.

Top 50 expressed genes in an unfiltered bulk RNA sequencing dataset of phox2a siblings.

“% of unfiltered 10x neurons” refers to gene detection in a single-cell atlas of neurons labeled in Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) (n=1,468 neurons). “% of projection neurons” refers to gene detection in a subset of the single-cell atlas containing projection neurons in the tangential nucleus (n=159 neurons). Data associated with Figure 5.

Top 50 differentially expressed genes in an unfiltered bulk RNA sequencing dataset of phox2a siblings and null mutants.

One star indicates a gene was retained in a filtered subset of projection neurons; %, evaluated using fluorescent in situ hybridization. “% of unfiltered 10x neurons” refers to gene detection in an unfiltered single-cell reference atlas of neurons labeled in Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) (n=1,468 neurons). Putative origin inferred from gene expression in the annotated 10x dataset (Methods, Figure S4)). Genes sorted by p adjusted value. Data associated with Figure 5.

Differentially expressed genes in phox2a-expressing medial vestibular neurons.

Star indicates a gene was evaluated in projection neurons using fluorescent in situ hybridization. # indicates a gene was significantly differentially expressed in projection neurons. “% of medial vestibular neurons” refers to detection in a subset of phox2a-expressing medial vestibular neurons in a single-cell reference atlas (Methods, Figure S4). “% of projection neurons with expression” refers to detection in a filtered subset of projection neurons. Gene sorted by p adjusted value. Data associated with Figure 5

We acknowledge the possibility that our in situ method is insufficiently quantitative to detect subtle differences in expression. Similarly, despite using both bulk and single-cell RNA sequencing approaches, we may lack the resolution to uncover differential gene expression within projection neurons, particularly given the inclusion of heterozygous phox2a siblings in controls. Nevertheless, consistent with functional and anatomical characterization, our sequencing data argues that projection neurons acquire the correct transcriptional profiles in the absence of motor partner-derived signals. Our findings are reminiscent of recent reports that the molecular signatures of spinal interneurons develop independently of motor partners 15.

Discussion

Here, we show that motor neurons are dispensable for fate specification in a canonical sensorimotor circuit. We first demonstrated that peripheral sensory and central projection neurons develop appropriate, directionally-selective connectivity and topography independently of their motor partners. Next, we established that projection neurons remain anatomically and molecularly competent to assemble with motor partners. Lastly, we show that loss of motor neurons does not meaningfully alter the transcriptional signatures of their pre-motor projection neuron partners. By providing causal evidence against an instructional role of motor partners for sensory connectivity, our work forces a revision of the current model for vestibulo-ocular reflex circuit formation. As proper connectivity across multiple synapses is foundational for proper function, our work speaks to mechanisms responsible for sensorimotor circuit assembly.

Transcriptional influences on motor neuron fate specification

While the primary focus of our work was circuit assembly, we found that, unexpectedly, phox2a acts in a dose-dependent manner to specify extraocular motor pool fate. Key evidence comes from phox2a heterozygotes, in which the earliest-born dorsal neurons in nIII are lost but later-born neurons in nIII/nIV are intact. This observation extends prior characterizations of phox2a mutations in zebrafish 34, chick 36, and human 37,38. Prior work hypothesized that phox2a dosage may regulate midbrain motor neuron differentiation into visceral and somatic types 36. In other systems, transcription 5254, growth 55 and axon guidance factors 56 can act in such a graded manner to regulate coarse cell type specification and wiring specificity. We extend these ideas to show that phox2a dose-dependency acts both over closely-related subtypes (pools within a single cranial nucleus) and along a temporal axis, where partial dosage preferentially targets the earliest-born neurons 33. Specifically, if phox2a is expressed in neural progenitor cells that give rise to nIII/nIV, then the earliest-born motor neurons would have the shortest exposure to phox2a.

Molecular insight into ocular motor neuron pool specification is sparse but would be welcome given the strong links between genetic development and ocular motor disease 38,5759. For example, subpopulation markers could resolve the topography of pools within dorsal nIII; whether IR/MR pools are spatially segregated or intermingled 33,60; whether the medial/lateral axis reflects functional differences among motor neuron subtypes; and whether/how local interactions between motor neuron pools contributes to fate specification 61. In spinal circuits, the rich molecular understanding of motor pool specification 47,6265 has enabled targeted perturbations of pool identity, allowing for major discoveries of their roles in circuit assembly 9,13,15,66. Our findings thus represent a step forward towards understanding how developmental deficits may contribute to ocular motor disorders 67.

Motor neurons: active or passive architects of pre-motor connectivity?

Our discoveries advance outstanding controversies over whether motor neurons actively or passively shape pre-motor connectivity. We find that extraocular motor neuron axons do not serve as “pioneers” 6872, with pre-motor axon targeting following passively from motor-derived pathfinding signals 2. Such a model predicts that projection neuron targeting would be entirely ablated after constitutive loss of extraocular motor neurons and their secreted signals 7377. Instead, we observed that projection neurons still establish long-range (hindbrain to midbrain) axonal projections, with appropriate spatial segregation that matches the topography of their motor partners 28,29,33. Our findings complement reports in spinal circuits that pre-motor targeting is grossly appropriate after manipulating the spatial source of, but not ablating, potential pathfinding signals 13, and that the transcriptional fate of premotor projection neurons similarly develops independently 15. We point to the late development of ocular musculature 78,79 compared to spinal musculature 2 as a potential source of the dispensability of muscle-derived signals.

Our work is also inconsistent with the strongest form of the “retrograde” hypothesis for vestibulo-ocular reflex circuit assembly. Originally, the retrograde model posited that motor neurons release a diffusable or cell-surface available signal that instructs pre-motor collaterals to sprout and then innervate specific pools, enabling behavioral specificity 3,31. Here, the proper spatial and temporal segregation of projection neuron axons suggests they remain poised to wire with spatially-appropriate (dorsal/ventral pools) targets. Additional evidence comes from incomplete phox2a knockouts (1-5% of nIII/nIV remaining), where projection neurons still form collaterals, though not robustly or reliably. We predict that projection neuron axons do not require a target-derived cue to grow, search, and synapse onto motor targets, and simply lack the adhesive contact necessary to stabilize nascent structures (reviewed in 45,46).

Nevertheless, extraocular motor neurons might still play an active or passive role in selecting and/or refining input specificity from their projection neuron partners. In spinal circuits, motor pool position passively imposes geometric constrains on pre-motor axon targeting 13,14, and manipulating the dendritic structure of motor neuron axons transforms input specificity 9,11. Genetic perturbations of nIII/nIV motor neuron position selectively compromise ocular responses to directional visual stimuli 61, though the circuit-level origin of such impairments is unclear. For the vestibulo-ocular reflex circuit, transforming all motor pools to the same fate or genetically “scrambling” pool position could resolve whether motor input specificity is truly hard-wired in projection neurons, or whether projection neurons instead target gross spatial domains irrespective of partner identity 13. Motor neuron-derived signals are of course capable of shaping their input by strengthening/weakening their inputs. Importantly, our results suggest that such signals will not define the fate of projection neurons, and by extension, circuit architecture.

We note that our study does not eliminate one additional source of post-synaptic partner signals to projection neurons. As in primates 80, projection neurons also contact neurons in the interstitial nucleus of Cajal, also known as the nucleus of the medial longitudinal fasciculus (INC/nMLF) 24. INC/nMLF neurons project early in development 81 to spinal circuits used for postural stabilization during swimming 8286. Notably, ablation of projection neurons disrupts postural stability 50. As we did not observe postural deficits in phox2a mutants, we infer that projection neuron connectivity to INC/nMLF targets is present and functional. Correspondingly, the development of projection neuron collaterals and synapses to INC/nMLF neurons appeared qualitatively normal in phox2a mutants, supporting our interpretation that projection neurons retain the capacity to properly assemble with post-synaptic targets even though similar structures to extraocular motor neurons are absent. In the future, if a similarly specific marker like phox2a is identified that labels the INC/nMLF, it will be possible to test whether these neurons play a role in vestibulo-ocular reflex circuit development.

Alternative mechanisms for fate specification and sensory input specificity in projection neurons

What is the origin of signals that govern projection neuron fate and sensory input specificity, if not motor-derived? In comparable systems, fate signals can be intrinsically-expressed or originate from extrinsic sources. For example, intrinsic genetic mechanisms assemble laminar connectivity in visual circuits 87,88 and facilitate sensorimotor matching in spinal circuits 9,8991. In directionally-selective retinal circuits, subtype fate is established in a similar manner 92,93. In “intrinsic” models, synaptic specificity arises from molecular matching between subtypes 9496. Alternatively, in somatosensory and auditory circuits, transcriptional fate depends on extrinsic signals such as growth factors 97 and sensation 98, respectively. In spinal circuits, positional fate, which constrains connectivity 13, is established by extrinsic codes such as morphogen gradients in early development 99 and Hox factors 63,64. In “extrinsic” models, early inputs are often erroneous and refined by activity 100,101 or molecular factors 102. Collectively, these findings offer two alternative models for how vestibulo-ocular reflex circuit assembly emerges.

The tight links between birth order, somatic position, and stimulus selectivity 28,29,33,85,103,104 across vestibulo-ocular reflex circuit populations support an “intrinsic” determination model. Further, neurogenesis and initial axon targeting develops contempora-neously for sensory afferents 55,103,105,106, projection neurons 3,29,107109, and extraocular motor neurons 3,33,39,110112, suggesting that neurons are poised to assemble with targets as early as their time of differentiation. Importantly, an “intrinsic specification” model makes a testable prediction about how and when sensory selectivity should emerge across the circuit: projection neurons and extraocular motor neurons should be directionally selective as soon as pre-synaptic input is established. Such evidence would justify future molecular inquiries into the underlying genetic factors, expanding early characterizations of the mechanisms that shape hindbrain topography 113,114, recent molecular profiling of the zebrafish hindbrain 115,116, and reports of molecular match-ing between extraocular motor neurons and muscle 117,118. Operationally, the present study lays a foundation for molecular explorations of projection neuron subtype determinants by establishing bulk- and single-cell transcriptomic profiling and in situ validation pipelines.

Conversely, evidence that stimulus selectivity emerges gradually would suggest that sensory afferents and/or projection neurons initially wire indiscriminately and that circuit connectivity is refined in time by extrinsic forces. Prior work in the vestibulo-ocular reflex circuit has proposed developmental roles for sensory-derived trophic factors 119 and activity-dependent refinement 120,121, though sensory afferents develop typically in the absence of utricular input 122 and ocular motor behavior does not depend on stimulus-driven activity 123. Here, an “extrinsic” determination model would predict that connectivity is established by an anterogradely-transmitted signal – that is, from sensory afferents to ocular muscles. If so, then future investigations might constitutively ablate sensory afferents to eliminate activity-driven, diffusible, or cell-surface instructional signals, similar to the present study. The directional bias in opsin-evoked activity in projection neurons 25, together with their transcriptional profiles established here, offer a clear readout of the role of sensory-derived factors. However, genetic targets exclusive to vestibular sensory afferents for gaze stabilization have not been identified, and tissue-specific genetic ablations remain limited in zebrafish. Looking ahead, resolving when and how stimulus selectivity emerges across the vestibulo-ocular reflex circuit will be key to understanding whether connectivity with pre- and/or post-synaptic partners instructs subtype fate, or whether subtype fate instructs connectivity.

Conclusion

Here, we discovered that motor partners do not determine pre-motor fate and sensory connectivity for the projection neurons that stabilize gaze. Our results overturn the current model that stimulus selectivity and connectivity are retrogradely specified, a major step towards understanding the origin, and eventually nature, of mechanisms that assemble an archetypal sensorimotor reflex circuit. Instead, our data support and extend recent models in spinal systems that motor partners do not actively construct sensory-to-interneuron reflex circuit architecture, but may later refine their inputs. By defining the contribution of motor neurons to specification and sensory connectivity of gaze-stabilizing central projection neurons, our work speaks to principles of sensorimotor circuit assembly.

Materials and methods

Resource availability

Lead Contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, David Schoppik (schoppik@gmail.com).

Materials Availability

Mutant fish lines generated in this study will be deposited to the Zebrafish International Resource Center (ZIRC).

Data and code availability

  • All data and code are deposited at the Open Science Framework and are publicly available at DOI: 10.17605/OSF.IO/93V6E

Experimental model and subject details

Fish care

All protocols and procedures involving zebrafish were approved by the New York University Langone School of Medicine Institutional Animal Care & Use Committee (IACUC). All larvae were raised at 28.5°C at a density of 20-50 larvae in 25-40 ml of buffered E3 (1mM HEPES added). Larvae used for photofill experiments were raised in constant darkness; all other fish were raised on a standard 14/10h light/dark cycle. Larvae for experiments were between 3-5 days post-fertilization (dpf).

Transgenic lines

Experiments were conducted on the mifta-/- background to remove pigment. All experiments used larvae from the F3 generation or older of a newly-created line of phox2a mutants (described below) on the following backgrounds: Tg(isl1:GFP) 40 to validate phox2a loss-of-function; Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16) 25,49 to drive UAS reporter expression; Tg(UAS-E1b:Kaede) 48 for anatomical imaging experiments; and Tg(UAS:GCaMP6s) 83 for calcium imaging experiments. All larvae were selected for brightness of fluorescence relative to siblings. Mendelian ratios were observed, supporting that selected larvae were homozygous for fluorescent reporter alleles.

Generation of phox2a mutants

phox2a mutant lines were generated using CRISPR/Cas9 mutagenesis. Two guide RNAs (gRNAs) were designed using the Benchling CRISPR Guide RNA Design Tool (see: key resources, Table 7). gRNAs were located towards the 5’ region of exon 1 to minimize the size of any translated protein. gRNAs were incubated with Cas9 protein before co-injection into Tg(isl1:GFP) embryos at the single cell stage. Injected embryos were screened for anatomical phenotypes (reduction in isl1-positive nIII/nIV motor neurons). Phenotypic embryos (F0) and their embryos were raised and genotyped via sequencing to identify and validate germline mutations. Three founders were identified and used for experiments: (1) phox2ad22 has a 22 bp deletion from base pairs 249 to 270, (2) phox2ad19 has a 19 bp deletion from base pairs 262 to 280, and (3) phox2ai2 has a 2 bp insertion (AG) from base pairs 261 to 262. Each mutation created a nonsense mutation, causing a predicted premature stop codon at the beginning of the homeobox. All alleles were validated using complementation assays, and larvae from all three alleles were used in experiments. For brevity, only one allele (phox2ad22) is shown in Figure 1.

Key Resources Table (associated with Methods).

Maintenance of phox2a adults

phox2a null larvae do not survive past 7 dpf. Sibling embryos (phox2a+/+ or phox2a+/-) were raised and genotyped to identify heterozygotes for line propagation. Primers for genotyping are listed in the Key Resources table (Table 7). Genomic DNA was amplified using a polymerase (DreamTaq PCR Master Mix 2X, Thermo Fisher Scientific K1071), 60°annealing temperature, 30 second elongation time, and 35 cycles of PCR. PCR generates a 169 bp product (wildtype), 147 bp product (phox2ad22), 150 bp product (phox2ad19), or 171 bp product (phox2ai2). phox2ad22 and phox2ad19 DNA was evaluated using gel electrophoresis; phox2ai2 was assessed via sequencing with the reverse primer (Genewiz, Azenta Life Sciences, South Plainfield, New Jersey).

Method details

Confocal imaging

Larvae were anesthetized in 0.2 mg/mL ethyl-3-aminobenzoic acid ethyl ester (MESAB, Sigma-Aldrich E10521, St. Louis, MO) prior to confocal imaging except where noted. Larvae were mounted dorsal side-up (axial view) or lateral side-up (sagittal view) in 2% low-melting point agarose (Thermo Fisher Scientific 16520) in E3. Images were collected on a Zeiss LSM800 confocal microscope with a 20x water-immersion objective (Zeiss W Plan-Apochromat 20x/1.0). Images of tangential nucleus soma and axons were acquired in a lateral mount with an 80×80 µm imaging window. Stacks spanned ∼30-40 µm, sampled every 1 µm. Images of nIII/nIV motor neurons were acquired in a dorsal mount with a 213×106 µm imaging window; stacks spanned approximately 90 µm, sampled every 1.5 µm. Images to validate nIII/nIV expression in a lateral mount were acquired using a 319×319 µm imaging window. Raw image stacks were analyzed using Fiji/ImageJ 124.

Identification of phox2a larvae

Prior to experiments, larvae were designated as phox2a mutants or sibling (wildtype/heterozygote) controls based on two criteria: gross loss of Tg(isl1:GFP) fluorescence in nIII/nIV at 2 dpf, visualized using a SugarCube LED Illuminator (Ushio America, Cypress CA) on a stereomicroscope (Leica Microsystems, Wetzlar, Germany) and absence of a swim bladder at 5 dpf. For anatomical and calcium imaging experiments, allele designations were validated using confocal imaging of nIII/nIV motor neurons: total or near-total loss of nIII/nIV neurons (null), selective loss of IR/MR neurons (heterozygote), or normal expression (wildtype). Designations were confirmed after experiments using genotyping. For RNA sequencing and fluorescent in situ experiments, sibling controls (wild-type/heterozygote) were combined.

Birthdating of nIII/nIV motor neurons

Early-born neurons in nIII/nIV were optically tagged using in vivo birthdating 29,33,125 on Tg(isl1:Kaede)ch103 larvae 126. Briefly, whole embryos were exposed to UV light for five minutes at experimenter-defined timepoints and subsequently raised in darkness to prevent background conversion. At 5 dpf, larvae were imaged on a confocal microscope. Neurons born before the time of photo-conversion expressed red, converted Kaede; neurons born after expressed only green, unconverted Kaede.

Fluorescent in situ hybridization and imaging

Experiments were performed using Hybridization Chain Reaction (HCR) for whole-mount zebrafish larvae 51,127. Probes were generated using the HCR 3.0 probe maker 128 using the sense sequence of the canonical gene cDNA from NCBI. All larvae were from the Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) background. Larvae were pre-identified as null mutants or siblings (wildtype or heterozygotes) and combined in equal ratios (8-10 larvae per condition, 16-20 larvae total) into a single 5 mL centrifuge tube for fixation and HCR. Larvae were fixed overnight with 4% PFA in PBS at 4° C and stored in 100% methanol at - 20° C. Subsequently, HCR was performed as described in 127, with adjustments to proteinase K incubation time based on age (2 dpf: 23 min incubation; 3 dpf: 30 min incubation; 5 dpf: 35 min incubation). HCR experiments used buffers and amplifiers from Molecular Instruments (Los Angeles, CA). DAPI staining was performed on some samples at 1:2000 from 5mg/ml stock solution, incubated overnight at 4° C. Samples were stored in 1x PBS at 4° C and imaged on a confocal microscope within four days. Prior to imaging, larvae were re-screened for Tg(isl1:GFP) fluorescence to identify null mutants and sibling controls. For each probe, imaging parameters were determined using a sibling control and kept constant for all subsequent larvae. Comparable settings (within 1% laser power) were used across probes.

Calcium imaging of tonic and impulse tilt stimuli responses

Experiments were performed as described in 29 using Tilt-In-Place Microscopy 41. All experiments used 5 dpf larvae from the Tg(isl1:GFP);(Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS:GCaMP6s) background. Briefly, larvae were mounted dorsal-up in 2% low-melt agarose in E3 onto a large beam diameter galvanometer system (ThorLabs GVS011). Tonic pitch-tilt stimuli were presented over a 65-second period in the following order: horizontal baseline (5 sec at 0°), nose-down tilt (15 sec at -19°), horizontal imaging (15 sec at 0°), nose-up tilt (15 sec at 19°), and horizontal imaging (15 sec at 0°). Impulse stimuli contained a 4 msec eccentric rotation, a 2 msec hold, and a 4 msec restoration step to horizontal and were presented twice over a 65-second imaging window: horizontal baseline (20 sec), impulse (10 msec), horizontal imaging (30 sec), impulse (10 msec), horizontal imaging (15 sec). Tonic and impulse stimuli were presented in alternating sets (impulse, then tonic) with a total of three stimulus set repeats.

Imaging was performed using a 20x water immersion objective (Olympus XLUMPLFLN20xW 20x/1.0), an infrared laser (Spectra-Physics MaiTai HP) at 920nm using 6.1-18.8 mW of power at the sample, and ThorLabs LS 3.0 software. Experiments were conducted in the dark. High-resolution anatomy scans of nIII/nIV motor neurons were performed for each experiment to validate allele designations. Scans used a 147×147 µm imaging window, a 90 µm stack sampled every 1.5 µm, and a 5.2 microsecond pixel dwell time. Anatomy scans of the tangential nucleus were acquired using a 148×91 µm imaging window as a 40-50 µm stack sampled every 1 µm. For stimulus imaging, the tangential nucleus was sampled every 3-6 µm based on cell density. 6-10 planes were sampled for each hemisphere. Ventral planes were imaged at higher magnification (112×68 µm imaging window) than dorsal planes (148×91 µm window) to avoid photomultiplier tube saturation from in-frame GFP fluorescence; magnification was corrected for in later analyses. Laser power was adjusted for each sampled plane due to the light scattering properties of zebrafish tissue. As greater power was required for ventral planes, imaging was always performed from ventral to dorsal to minimize photobleaching effects. Stimulus imaging was performed at 3 frames/second (2.2 µs pixel dwell time) with a total time of approximately two hours per fish.

The number of neurons sampled per fish for each genotype are as follows (n=mean/standard deviation per fish; total cells/total fish): (1) Wildtype reference dataset: n=22±14 neurons/fish, N=255/10 total; (2) phox2a+/+: 33±13 neurons/fish, N=76/5 total; (3) phox2a+/-: 43±9 neurons/fish, N=109/6 total; (4) phox2a-/-: 35±14neurons/fish, N=297/16 total.

Retrograde photolabeling of tangential nucleus neurons

Experiments were performed as described in 29 based on 44 on 5 dpf larvae from the Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) background. Briefly, experiments leveraged a photoconvertible protein, Kaede, which irreversibly converts from green to red with ultraviolet light. Larvae were raised in darkness to minimize background conversions. Larvae were mounted dorsal-up in 2% agarose under a confocal microscope. An imaging window was centered over the medial longitudinal fasciculus (MLF) and repeatedly scanned with a 405 nm laser for 30 seconds until fully converted (green to red). Off-target photo-conversion was assessed (e.g., conversion of projections lateral to the MLF). Larvae were unmounted, left to recover in E3 for 4 hours in darkness, and then re-mounted in a lateral mount. An imaging window was centered around the tangential nucleus (see: Confocal Imaging). Retrogradely-labeled soma were identified by their center-surround fluorescence appearance: red converted cytoplasm surrounding an unconverted green nucleus.

Neuron harvesting, dissociation, and flow cytometry

Experiments were performed on 72-74 hpf larvae from the Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) background. At 2 dpf, larvae were designated as null or sibling (wildtype/heterozygote) as described above. Three experimenters (D.G., K.R.H., and P.L) harvested neurons in parallel. Larvae were anesthetized in MESAB in Earle’s Balanced Salt Solution with calcium, magnesium, and phenol red (EBSS, Thermo Fisher Scientific 24010043) and mounted dorsal-up in 2% agarose. Fluorescence in tangential nucleus neurons was visualized using a SugarCube LED Illuminator (Ushio America, Cypress CA) using 10x eyepieces on a stereomicroscope (Leica Microsystems, Wetzlar, Germany). Neurons were harvested using a thin wall glass capillary tube (4 inch, OD 1.0 MM, World Precision Instruments) into EBSS in a non-stick Eppendorf tube and kept on ice until dissociation.

Neurons were dissociated in 20 units/mL of papain prepared in EBSS (Worthington Biochemical), 2000 units/mL of deoxyribonucleic prepared in EBSS (Worthington Biochemical), and 100 mg/mL of Type 1A Collagenase (Sigma Aldrich) prepared in Hanks Buffered Salt Solution without calcium/magnesium (HBSS, Thermo Fisher Scientific). Neurons were incubated for 45 minutes at 31.5°C with a gentle vortex every 10-15 min, then passed through a 20 µm filter and centrifuged for 10 mins at 300 x g. After removing supernatant, neurons were resuspended in L15 (Thermo Fisher Scientific) with 2% fetal bovine serum (Thermo Fisher Scientific). Cell health was evaluated using DAPI, applied at 0.5 µg/ml (Invitrogen) and incubated on ice for 30-45 mins prior to flow cytometry.

Flow cytometry was performed using a Sony SH800z cell sorter (100 µm nozzle, 20 psi) to isolate single neurons (Figure S3). Three controls were run: (1) non-fluorescent wildtype neurons, (2) non-fluorescent neurons + DAPI, (3) fluorescent (green) neurons from Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) + DAPI. On average, 2% of neurons were DAPI-positive and excluded. Neurons were evaluated for positive (green) fluorescence. Fluorescence was not evaluated to separate Tg(UAS-E1b:Kaede) neurons from those labeled by Tg(isl1:GFP). Neurons were sorted into an Eppendorf tube containing 700 µl of lysis buffer (RNAqueous Micro Total RNA Isolation Kit, Thermo Fisher Scientific) for downstream bulk RNA sequencing.

Bulk RNA sequencing

RNA isolation was performed using an RNAqueous Micro Total RNA Isolation Kit (Thermo Fisher Scientific). RNA concentration and quality (RIN > 8.0) was evaluated using an RNA 6000 Pico Kit and a 2100 Bioanalyzer system (Agilent Technologies, Santa Clara, California). RNA sequencing was performed by the NYU Genome Technology Center. Libraries were prepared using the low-input Clontech SMART-Seq HT with Nxt HT kit (Takara Bio USA) and sequenced using an Illumina NovaSeq 6000 with an S1 100 Cycle Flow Cell (v1.5).

Quantification and statistical analysis

Cell counting and spatial mapping of nIII/nIV motor neurons

Analysis was performed in Fiji/ImageJ 124 using the Cell Counter plugin. Anatomical stacks of nIII/nIV were subdivided in the dorsoventral axis as described in 33 to facilitate localization. A point ROI was dropped over each neuron in the plane in which the soma was brightest (center). The number of neurons in each dorsoventral plane and their coordinates were recorded. Neuron coordinates were standardized relative to a (0,0) point, defined as one corner of a standard-sized rectangular box centered over the extent of nIII/nIV in a maximum intensity projection. Differences in spatial location across genotypes was evaluated separately for each spatial axis using a two-tailed, two-sample Kolmogorov-Smirnov test. Probability distributions for figures were generated using the mean and standard deviation from bootstrapped data (n=100 iterations) to ensure results were robust to data from single larva.

Analysis of calcium imaging experiments

Analysis methods are detailed in 29 and summarized briefly here. Regions of Interest (ROIs) were drawn around tangential nucleus neurons for each stimulus plane sampled and adjusted for minor movement (1-2 µm) between trials. Raw fluorescence traces were extracted using Matlab R2020b (MathWorks, Natick, Massachusetts) and normalized by ROI size to account for variation in magnification. A neuron’s response to tonic or impulse stimuli was defined as the change in fluorescence in the first second of restoration to horizontal following tilt delivery. Responses were normalized using a baseline period, defined as the mean fluorescence across the initial baseline window (5 sec) preceding the nose-down tilt (nose-down response) or the last 3 sec of the horizontal restoration following nose-down tilt (nose-up response). This was used to generate a ΔF/F value. A ΔF/F response was defined as significant if it was greater than two standard deviations above baseline. Directional selectivity was assigned by normalizing the difference in ΔF/F responses to each tilt by their sum. This generated a scale of values of ±1 (i.e., positive values represent nose-up selectivity; negative values, nose-down). Some neurons responded to both tilt directions with high similarity; we set a minimum threshold of abs(0.1) to distinguish neurons with a clear directional selectivity from untuned neurons.

Spatial mapping of tangential nucleus neurons

Analysis methods are detailed in 29 and summarized briefly here. All imaged neurons were manually registered to a reference framework using Adobe Illustrator (2021). Anatomy stacks from all experiments were aligned in the XY (rostrocaudal, mediolateral) axes using established anatomical landmarks (e.g., Mauthner cell body, medial longitudinal fasciculus, otic capsule). For Z-registration (dorsoventral axis), stacks were subdivided into eight sections using landmarks within and around the tangential nucleus (e.g., Mauthner cell body, neuropil). All registered images were verified by two independent observers (D.G. and S.H.). Neurons were localized to one dorsoventral section and a reference circle, representing a cell, was placed in Illustrator. Coordinates for each reference circle were recorded and standardized to an absolute (0,0) point (dorsomedial-most point of the tangential nucleus). Coordinates were imported into Matlab (R2020b) and used to generate a spatial map of imaged neurons.

Statistical analysis of differences in tilt responses across phox2a genotypes

Statistical comparisons of tonic and impulse tilt responses are summarized in Table 1. Analyses used a one-way analysis of variance with multiple comparisons. No significant differences (tonic tilt responses) or small differences (impulse responses) were observed across genotypes. Control data reported in Results and Figure 2-Figure 3 is an aggregate from wildtype, phox2a+/+, and phox2a+/- larvae.

Alignment, quality control, and differential expression analysis of bulk sequencing data

Initial alignment and analyses were performed by the Applied Bioinformatics Laboratories at the NYU School of Medicine (RRID:SCR_019178). Sequencing data was aligned to the GRCz11 zebrafish reference genome and two fluorescent markers (Kaede, GFP; NCBI). Eight datasets from four experimental repeats were aligned: four from phox2a mutants, and four from sibling controls. One experimental repeat had significantly higher variance in the first and second principal components, likely due to poor quality leading to extremely low transcript counts, and was excluded from downstream analyses. Number of cells/larvae sequenced and used in downstream analysis are as follows: Repeat 1, n=532/n=904 cells from N=28/N=28 phox2a null/control lar-vae; Repeat 2, n=802/n=683 cells from N=27/N=26 phox2a null/control larvae; Repeat 3, n=1000/n=1007 cells from N=41/N=40 phox2a null/control larvae; Repeat 4 (excluded): n=690/n=571 cells from N=33/N=33 phox2a null/control larvae Differential gene expression between conditions (phox2a mutants vs. sibling controls) was assessed using DESeq2129. Differentially-expressed candidate genes met two criteria: log2 fold change >|2| and p adjusted < 0.05.

Filtering of bulk sequencing data using a reference single-cell sequencing dataset

Analyses were performed in R. Detection of markers for motor neurons (isl1, isl2a, isl2b) 130,131 and neurons caudal (hoxd4a) 113,132 and lateral (barhl2) 114 to rhombomeres 4-6 supported that our dataset included other populations. We applied a filter to exclude erroneous gene expression from non-tangential nucleus populations.

Filtering was performed using an existing single-cell atlas of neurons labeled in Tg(−6.7Tru.Hcrtr2:GAL4-VP16);(Tg(UAS-E1b:Kaede), generated with 10x Genomics. The reference atlas was generated from four experimental samples using the harvest, dissociation, and flow cytometry method described above. The sequenced atlas contained 1,468 neurons (Figure S4A-Figure S4B)). Data was analyzed using Seurat v4.0 133. Cluster annotation was performed using a combination of fluorescent in situ hybridization as described above (Figure S4C-Figure S4E and other data not shown) and published molecular data of the zebrafish hindbrain 113. n=159 neurons (11%) were validated as excitatory projection neurons from the tangential nucleus.

Genes in the bulk dataset were only included in downstream analyses if they were expressed above threshold percent of reference projection neurons: 1%, 3%, 5%, 10%, 30%, or 50%. The most stringent filter (50%) was set using the transcription factor evx2, which is reported to be expressed in all tangential nucleus neurons 50 and was detected in 36% of reference projection neurons. Qualitatively, we found that gene detection with fluorescent in situ hybridization scaled with reference filter stringency (Figure S5). Analyses were performed separately for each threshold. The total number of genes included for downstream analyses for each threshold are as follows: 28,807 (no threshold), 11,278 (1% of reference neurons), 8,346 (3%), 6,701 (5%), 4,185 (10%), 1,105 (30%), 384 (50%). We used the following significance thresholds for differential gene expression in filtered datasets: adjusted p value < 0.05 and abs(log2FoldChange) > 2. The number of differentially expressed genes for each threshold was as follows: 100 (no threshold), 23 (1% of reference neurons), 6 (3%), 2 (5%), 0 (10%).

Projection neurons in the tangential nucleus lie in close proximity to the medial vestibular nucleus (3-5 µm) to the medial edge of the tangential nucleus and 10-20 µm dorsal to the rostral/dorsal edge of the tangential nucleus). Some MVN neurons express phox2a (Figure S7). Our reference single-cell atlas isolated a small cluster of phox2a-expressing neurons in r5-6 that likely originates from this population. To control for the possibility that some differentially expressed genes are localized to the MVN, and not projection neurons, we also evaluated differential gene expression in the phox2a-expressing subset of MVN neurons (n=45 neurons, 3% of reference dataset). Data is shown in Figure S7.

Generation of representative images for fluorescent in situ hybridization

Images were generated using Fiji/ImageJ 124. An anatomical template of the tangential nucleus was generated based on 29. Briefly, for sagittal view images, a 30-µm stack was centered over the tangential nucleus. For each plane, a region of interest (ROI) was drawn over all cells within the bounds of the tangential nucleus. Transcript expression outside the ROI was masked. Maximum intensity projections were generated. Minimal or no alterations to brightness/contrast were made for probe expression given the correlation between fluorescence intensity and detected transcript 51.

Additional statistics

Bias and variability in probability distributions were estimated by bootstrapping, or resampling the raw distributions with replacement 134. Data shown is the mean and standard deviation of 100 bootstrapped distributions. Topography data was evaluated using two-tailed, two-way Kolmogorov-Smirnov tests. Functional responses to tilts (i.e., calcium response strength, directionality index) were evaluated using two-tailed Wilcoxon rank sum tests. Differences in responses across genotypes were analyzed using one-way analysis of variance tests.

Acknowledgements

Research was supported by the National Institute on Deafness and Communication Disorders of the National Institutes of Health under award numbers R01DC017489, F31DC020910, and F31DC019554, the National Institute of Neurological Disorders and Stroke under award numbers F99NS129179, T32NS086750, and the National Cancer Institute P30CA016087. The authors would like to thank Hannah Gelnaw for assistance with fish care, and Jeremy Dasen, Claude Desplan, Katherine Nagel, Dan Sanes, along with the members of the Schoppik and Nagel labs for their valuable feedback and discussions. Finally, the authors gratefully acknowledge the late Hans Straka for his generous insights and encouragement throughout.

Author contributions

Conceptualization: DG and DS, Methodology: DG, BR, KRH, PL, and DS, Investigation: DG, BR, KRH, PL, ML, HP, HG, SH, and CQ, Visualization: DG, Writing: DG, Editing: DS, Funding Acquisition: DG and DS, Supervision: DS.

Author competing interests

The authors declare no competing interests.

Note

This reviewed preprint has been updated to correct the affiliations.

Supplementary figures

Extraocular motor neuron fate is perturbed in phox2a mutants

Associated with Figure 1.

(A) Fluorescent in situ hybridization images of nIII/nIV motor neurons, labeled by Tg(isl1:GFP), in wildtype larvae at 2 dpf. Panels from left to right show: (1) isl1+ motor neurons, (2) mRNA for vachta, a marker for cholinergic motor neurons, (3) DAPI, and (4) merge. White dashed lines show the extent of nIII/nIV.

(B) Loss of both isl1 and vachta identity in nIII/nIV motor neurons in phox2a mutants at 2 dpf. All scale bars, 20 µm.

phox2a specifies nIII motor neuron fate in a dose- and birthdate-dependent manner.

Associated with Figure 1.

(A) Images of nIII/nIV motor neurons, labeled in Tg(isl1:GFP), in wildtype siblings (left) and phox2a heterozygotes (middle) at 5 dpf. Wildtype image same as in Figure 1E. One hemisphere shown. White dashed lines outline the dorsal extent of nIII, which contains inferior rectus and medial rectus neurons 33. Scale bar, 20 µm.

(B) Location of the earliest-born neurons in nIII/nIV (left, magenta) against all nIII/nIV neurons labeled in Tg(isl1:Kaede) (right, grey). Larvae birthdated at 34 hpf (Methods). One hemisphere shown. White dashed lines outline the dorsal extent of nIII. Scale bar, 20 µm.

(C) Quantification of the number of Tg(isl1:GFP)+ neurons in nIII/nIV from N=6 wildtype siblings (grey) and N=8 phox2a heterozygotes (teal). Wildtype data same as Figure 1F.

(D) Distributions showing probability of nIII/nIV soma location across each spatial axis in wildtype (black) and heterozygous (teal) phox2a larvae. Solid and shaded lines show mean and standard deviation, respectively, from bootstrapped data. Data from same fish quantified in Figure S2C. ns, not significant; star, significant at the p<0.001 level.

Flow cytometry gating strategy to sort fluorescently-labeled neurons for bulk RNA sequencing.

Associated with Figure 5.

(A) Sequential gates used to sort fluorescent neurons labeled with Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede);Tg(isl1:GFP). Gate A excluded presumptive debris (small cells). Gate B isolated single cells and excluded large cells and doublets. Gate C excluded DAPI+ (dead or unhealthy) neurons. Gate D isolated fluorescent (GFP or Kaede+) neurons; neurons in this gate were sorted. Gates were set using negative controls (not shown; Methods). Gates shown for one of four experimental repeats.

Molecular identification of projection neurons using a reference single-cell RNA sequencing atlas.

Associated with Figure 5.

(A) UMAP visualization of a single-cell RNA sequencing atlas of n=1,468 neurons labeled in Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede), generated with 10x Genomics (Methods). Each circle is a single neuron. Neurons are clustered (colors) according to their transcriptional identity. Annotations are based on validated marker genes (data not shown). TAN, tangential nucleus; MVN, medial vestibular nucleus; r, rhombomere; DON, dorsal octavolateral nucleus. Note: Cluster 2 is broadly inclusive of lateral, excitatory vestibular projection neurons near the ventral base of the otic capsule. It includes, but is not exclusive to, the tangential nucleus.

(B) Heatmap showing genes unique to each annotated cluster. Each row is a gene; names unlisted for clarity. Columns show distinct clusters. Color bar on top reflects clusters in Figure S4A. Yellow and purple reflect stronger or weaker gene expression, respectively. Black outlines show two populations of interest: tangential nucleus neurons and medial vestibular nucleus neurons.

(C) Fluorescent in situ hybridization against three markers (pou4f1, penkb, chrna2b) that are negative for tangential nucleus projection neurons and positive for medial vestibular nucleus neurons. Top row shows RNA expression (green); bottom row, merge with neurons labeled in Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede). Dashed lines outline the tangential nucleus (TAN) and medial vestibular nucleus (MVN). Data from 72 hpf larvae. Images shown in an axial view.

(D) Fluorescent in situ hybridization against a positive marker (evx2) for both tangential nucleus and medial vestibular nucleus neurons. All scale bars, 20 µm.

Visualization of transcripts in siblings and phox2a null mutants with fluorescent in situ hybridization is (1) consistent across larvae and (2) scales with predicted detection in projection neurons.

Associated with Figure 5.

(A-A’) Fluorescent in situ hybridization against itga9 for three sibling (A) or phox2a null mutant (A’) larvae (72 hpf), imaged with identical conditions. Left column shows RNA (green); right column, merge with projection neurons visualized with Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) (grey). Dashed lines outline the projection nucleus. Cell and transcript expression outside the projection nucleus is removed for visual clarity. Percentage (1.9%) refers to fraction of cells in a single-cell RNA sequencing reference atlas (Methods) with detected transcript. All scale bars, 20 µm.

(B-B’) Fluorescent in situ hybridization against htt for three sibling (B) and phox2a mutant (B’) larvae (72 hpf).

(C-C’) Fluorescent in situ hybridization against evx2 for three sibling (C) and phox2a mutant (C’) larvae (72 hpf).

Differential gene expression in an unfiltered bulk sequencing dataset of siblings and phox2a mutants.

Associated with Figure 5.

(A) Volcano plot showing differentially expressed genes across an unfiltered bulk RNA sequencing dataset. Dashed lines represent significance cutoffs: horizontal line, p adjusted>0.05; vertical line, Log2 Fold Change > 2.0. Each circle is a gene. Genes to the left and right of 0 on the horizontal axis show downregulated and upregulated genes, respectively. Red color shows genes that are differentially expressed in a filtered subset of projection neurons (Figure 5). Grey-colored genes are below both significance thresholds.

(B) Same data as Figure S6A, now highlighting candidate genes evaluated by fluorescent in situ (Figure 5) with red. Two candidates (yellow) that did not meet projection neuron filter criteria (Methods) are shown in Figure S6D; remaining candidates (included in filtered data) shown in Figure 5F.

(C) Same data as Figure S6A-Figure S6B, showing the number of differentially expressed genes at progressive significance thresholds (p adjusted). Red and blue lines show the number of significantly upregulated and downregulated genes, respectively.

(D) Fluorescent in situ hybridization against two differentially expressed genes, slc22a7a (log2 fold change=10.2, p adj.=1.6×10−4) and bckdhbl (log2 fold change=-9.1, p adj.=0.001), that did not meet projection neuron filter criteria. Percentage refers to fraction of all singly-profiled hindbrain vestibular neurons with expression (Methods). Left columns show RNA (green); right columns, merge with projection neurons labeled with Tg(−6.7Tru.Hcrtr2:GAL4-VP16);Tg(UAS-E1b:Kaede) (grey). Dashed lines outline the projection nucleus. Cell and transcript expression outside the projection nucleus is removed for visual clarity. All scale bars, 20 µm.

phox2a expression in the medial vestibular nucleus may underscore differential gene expression phenotypes in bulk data.

Associated with Figure 5.

(A) Fluorescent in situ hybridization against phox2a in a 5 dpf larvae (axial view). Top panel shows phox2a RNA (green); bottom panel, merge with neurons visualized with Tg(isl1:GFP);Tg(−6.7Tru.Hcrtr2:GAL4-VP16(;Tg(UAS-E1b:Kaede) (grey). White dashed lines outline three nuclei of interest: projection neurons in the tangential nucleus (TAN), the medial vestibular nucleus (MVN), and the facial nucleus (nVII). All scale bars, 20 µm.

(B) Volanco plot showing differentially expressed genes in phox2a-expressing medial vestibular nucleus neurons between control and phox2a null larvae at 3 dpf. Dashed lines represent significance cutoffs: horizontal line, p>0.05; vertical line, Log2 Fold Change > 2.0. Each circle is a gene. Genes to the left and right of 0 on the horizontal axis show downregulated and upregulated genes, respectively. Colors indicate percent of reference phox2a-expressing medial vestibular neurons (Methods, Figure S4A-Figure S4B) that express a given gene. Grey-colored genes are below both significance thresholds. Data for all differentially expressed candidates in medial vestibular neurons shown in Table 6.

(C) Same data as Figure S7B. Color shows genes that are differentially expressed in both medial vestibular nucleus neurons and projection neurons.