Abstract
Sensorimotor reflex circuits engage distinct neuronal subtypes, defined by precise connectivity, to transform sensation into compensatory behavior. Whether and how motor neuron populations specify the subtype fate and/or sensory connectivity of their pre-motor partners remains controversial. Here, we discovered that motor neurons are dispensable for proper connectivity in the vestibular reflex circuit that stabilizes gaze. We first measured activity following vestibular sensation in premotor projection neurons after constitutive loss of their extraocular motor neuron partners. We observed normal responses and topography indicative of unchanged functional connectivity between sensory neurons and projection neurons. Next, we show that projection neurons remain anatomically and molecularly poised to connect appropriately with their down-stream partners. Lastly, we show that the transcriptional signatures that typify projection neurons develop independently of motor partners. Our findings comprehensively overturn a long-standing model: that connectivity in the circuit for gaze stabilization is retrogradely determined by motor partner-derived signals. By defining the contribution of motor neurons to specification of an archetypal sensorimotor circuit, our work speaks to comparable processes in the spinal cord and advances our understanding of principles of neural development.
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 2–6, but controversy remains about the nature of their role. In the spinal cord, molecular perturbations of motor neuron identity have provided evidence both for 7–11 and against 12–16 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,22–29. 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 24–26,29 given the ease of optical imaging, abundant tools for genetic perturbations, rapid development, and robust evolutionary conservation.
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 34–37. 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.
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.
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.
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.
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.
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.
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.
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 52–54, 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,57–59. 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,62–65 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” 68–72, 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 73–77. 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 82–86. 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,89–91. 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 94–96. 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,107–109, and extraocular motor neurons 3,33,39,110–112, 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.
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.
Note
This reviewed preprint has been updated to correct the affiliations.
Supplementary figures
References
- 1.The Integrative Action of the Nervous SystemYale University Press
- 2.A critical period for the specification of motor pools in the chick lumbosacral spinal cordDevelopment 122
- 3.The development of vestibulo-ocular circuitry in the chicken embryoJournal of Physiology-Paris 97:17–25
- 4.Motor circuits in action: Specification, connectivity, and functionNeuron 74:975–989
- 5.Master or servant? emerging roles for motor neuron subtypes in the construction and evolution of locomotor circuitsCurrent Opinion in Neurobiology 42:25–32
- 6.Assembly of motor circuits in the spinal cord: Driven to function by genetic and experience-dependent mechanismsNeuron 56:270–283
- 7.Primacy of flexor locomotor pattern revealed by ancestral reversion of motor neuron identityCell 162
- 8.Spinal locomotor circuits develop using hierarchical rules based on motorneuron position and identityNeuron 87:1008–1021
- 9.Target-induced transcriptional control of dendritic patterning and connectivity in motor neurons by the ETS gene pea3Cell 127:1439–1452
- 10.Columnar-intrinsic cues shape premotor input specificity in locomotor circuitsCell Reports 21:867–877
- 11.Positional strategies for connection specificity and synaptic organization in spinal sensory-motor circuitsNeuron 102:1143–1156
- 12.The waiting period of sensory and motor axons in early chick hindlimb: Its role in axon pathfinding and neuronal maturationThe Journal of Neuroscience 20
- 13.Patterns of spinal sensory-motor connectivity prescribed by a dorsoventral positional templateCell 147:653–665
- 14.Spinal inhibitory interneuron diversity delineates variant motor microcircuitsCell 165:207–219
- 15.Origin and segmental diversity of spinal inhibitory interneuronsNeuron 97:341–355
- 16.Intrinsic control of neuronal diversity and synaptic specificity in a proprioceptive circuiteLife 9
- 17.Decoding the organization of spinal circuits that control locomotionNature Reviews Neuroscience 17
- 18.The dual nature of the efferent nervous system: A further study of the static and kinetic systems, their function and symptomatology contentsArchives of Neurology & Psychiatry 10
- 19.Lim homeodomain factors lhx3 and lhx4 assign subtype identities for motor neuronsCell 95
- 20.Islet1 and islet2 have equivalent abilities to promote motoneuron formation and to specify motoneuron subtype identityDevelopment 133
- 21.Pathways and synaptic articulation patterns connecting vestibular receptors and oculomotor nucleiIn The Oculomotor System Hoeber Medical Division, Harper and Row :205–223
- 22.Rhombomeric organization of vestibular pathways in larval frogsJournal of Comparative Neurology 437:42–55
- 23.The relationship between hodological and cytoarchitectonic organization in the vestibular complex of the 11-day chicken embryoThe Journal of Comparative Neurology 457:87–105
- 24.The tangential nucleus controls a gravito-inertial vestibulo-ocular reflexCurrent Biology 22:1285–1295
- 25.Gaze-stabilizing central vestibular neurons project asymmetrically to extraocular motoneuron poolsThe Journal of Neuroscience 37:11353–11365
- 26.Development of vestibular behaviors in zebrafishCurrent Opinion in Neurobiology 53:83–89
- 27.Segmental analysis of the vestibular nerve and the efferents of the vestibular complexThe Anatomical Record 302:472–484
- 28.Organization of the gravity-sensing system in zebrafishNature Communications 13
- 29.Neuronal birthdate reveals topography in a vestibular brainstem circuit for gaze stabilizationCurrent Biology 33:1265–1281
- 30.Neuroepithelial ‘compartments’ and the specification of vestibular projectionsIn Progress in Brain Research Elsevier :3–21
- 31.Ontogenetic rules and constraints of vestibulo-ocular reflex developmentCurrent Opinion in Neurobiology 20:689–695
- 32.α2-chimaerin regulates a key axon guidance transition during development of the oculomotor projectionJournal of Neuroscience 33:16540–16551
- 33.Extraocular motoneuron pools develop along a dorsoventral axis in zebrafish, Danio rerioJournal of Comparative Neurology 525:65–78
- 34.Development of noradrenergic neurons in the zebrafish hindbrain requires BMP, FGF8, and the homeodomain protein soulless/phox2aNeuron 24:555–566
- 35.Reciprocal gene replacements reveal unique functions for phox2 genes during neural differentiationThe EMBO Journal 24
- 36.PHOX2a regulation of oculomotor complex nucleogenesisDevelopment 137:1205–1213
- 37.Homozygous mutations in arix(phox2a) result in congenital fibrosis of the extraocular muscles type 2Nature Genetics 29
- 38.Neurological features of congenital fibrosis of the extraocular muscles type 2 with mutations in PHOX2aBrain 129:2363–2374
- 39.Differential expression of LIM homeobox genes among motor neuron subpopulations in the developing chick brain stemMolecular and Cellular Neuroscience 8:242–257
- 40.Visualization of cranial motor neurons in live transgenic zebrafish expressing green fluorescent protein under the control of the Islet-1 promoter/enhancerJournal of Neuroscience 20:206–218
- 41.Tilt in place microscopy: a simple, low-cost solution to image neural responses to body rotationsThe Journal of Neuroscience 43
- 42.Ultrasensitive fluorescent proteins for imaging neuronal activityNature 499:295–300
- 43.Extraocular motor nuclei: location, morphology and afferentsReviews of Oculomotor Research 2:81–117
- 44.Chronology-based architecture of descending circuits that underlie the development of locomotor repertoire after birtheLife 8
- 45.Cell adhesion molecules: signalling functions at the synapseNature Reviews Neuroscience 8
- 46.To stick or not to stick: The multiple roles of cell adhesion molecules in neural circuit assemblyFrontiers in Neuroscience 16
- 47.Determinants of motor neuron functional subtypes important for locomotor speedCell Reports 42
- 48.Targeting neural circuitry in zebrafish using GAL4 enhancer trappingNature Methods
- 49.A convergent and essential interneuron pathway for mauthner-cell-mediated escapesCurrent Biology 25:1526–1534
- 50.Biomechanics and neural circuits for vestibular-induced fine postural control in larval zebrafishNature Communications 14
- 51.Third-generation in situ hybridization chain reaction: multiplexed, quantitative, sensitive, versatile, robustDevelopment 145
- 52.Regulation of b lymphocyte and macrophage development by graded expression of pu.1Science 288
- 53.Induction of oligodendrocyte differentiation by olig2 and sox10: Evidence for reciprocal interactions and dosage-dependent mechanismsDevelopmental Biology 302
- 54.The level of the transcription factor pax6 is essential for controlling the balance between neural stem cell self-renewal and neurogenesisPLoS Genetics 5
- 55.A spatial and temporal gradient of fgf differentially regulates distinct stages of neural development in the zebrafish inner earPLoS Genetics 8
- 56.Graded expression of semaphorin-1a cell-autonomously directs dendritic targeting of olfactory projection neuronsCell 128
- 57.Loss of mafb function in humans and mice causes duane syndrome, aberrant extraocular muscle innervation, and inner-ear defectsThe American Journal of Human Genetics 98
- 58.Human cfeom1 mutations attenuate kif21a autoinhibition and cause oculomotor axon stallingNeuron 82
- 59.An inherited tubb2b mutation alters a kinesin-binding site and causes polymicrogyria, cfeom and axon dysinnervationHuman Molecular Genetics 21
- 60.Extra- and intracellular HRP analysis of the organization of extraocular motoneurons and internuclear neurons in the guinea pig and rabbitThe Journal of Comparative Neurology 262:429–445
- 61.Cadherins regulate nuclear topography and function of developing ocular motor circuitryeLife 9
- 62.Functionally related motor neuron pool and muscle sensory afferent subtypes defined by coordinate ets gene expressionCell 95
- 63.A hox regulatory network establishes motor neuron pool identity and target-muscle connectivityCell 123:477–491
- 64.Hox repertoires for motor neuron diversity and connectivity gated by a single accessory factor, FoxP1Cell 134:304–316
- 65.Parallel pbx -dependent pathways govern the coalescence and fate of motor columnsNeuron 91:1005–1020
- 66.Hox genes: Choreographers in neural development, architects of circuit organizationNeuron 80
- 67.Genetics of StrabismusSpringer International Publishing
- 68.Involvement of distinct pioneer neurons in the formation of layer-specific connections in the hippocampusThe Journal of Neuroscience 18
- 69.Subplate neurons pioneer the first axon pathway from the cerebral cortexScience 245
- 70.Pathfinding by peripheral pioneer neurons in grasshoppersScience 218
- 71.Pioneer neurones in an insect embryoNature 260
- 72.Pathfinding by zebrafish motoneurons in the absence of normal pioneer axonsDevelopment 114
- 73.Altered sensory projections in the chick hind limb following the early removal of motoneuronsDevelopmental Biology 118
- 74.Skin sensory innervation patterns in embryonic chick hindlimbs deprived of motoneuronsDevelopmental Biology 126
- 75.Sensory nerve routes in chick wing buds deprived of motor innervationDevelopment 95
- 76.Different subsets of axonal guidance cues are essential for sensory neurite outgrowth to cutaneous and muscle targets in the dorsal ramus of the embryonic chickJournal of Experimental Zoology 251
- 77.A transient population of neurons pioneers the olfactory pathway in the zebrafishThe Journal of Neuroscience 18
- 78.The development of eye movements in the zebrafish (Danio rerio)Developmental psychobiology 31:267–276
- 79.Differentiation of avian craniofacial muscles: I. patterns of early regulatory gene expression and myosin heavy chain synthesisDevelopmental Dynamics 216:96–112
- 80.Anatomical and physiological characteristics of vestibular neurons mediating the vertical vestibulo-ocular reflexes of the squirrel monkeyThe Journal of Comparative Neurology 264:571–594
- 81.Development of reticulospinal neurons of the zebrafish. ii. early axonal outgrowth and cell body positionJournal of Comparative Neurology 251
- 82.Neural control and modulation of swimming speed in the larval zebrafishNeuron 83:692–707
- 83.Descending control of swim posture by a midbrain nucleus in zebrafishNeuron 83:679–691
- 84.Selective responses to tonic descending commands by temporal summation in a spinal motor poolNeuron 83:708–721
- 85.Tiltable objective microscope visualizes discrimination of static and dynamic head movement originates at hair cellsbioRxiv
- 86.Brainstem circuits encoding start, speed, and duration of swimming in adult zebrafishNeuron 111:372–386
- 87.SidekicksCell 110
- 88.Dscam and sidekick proteins direct lamina-specific synaptic connections in vertebrate retinaNature 451
- 89.Ets gene er81 controls the formation of functional connections between group ia sensory afferents and motor neuronsCell 101
- 90.Specificity of sensory–motor connections encoded by sema3e–plxnd1 recognitionNature 459:842–846
- 91.Motor axon pathfindingCold Spring Harbor Perspectives in Biology 2
- 92.Type ii cadherins guide assembly of a direction-selective retinal circuitCell 158
- 93.The transcription factor tbx5 regulates direction-selective retinal ganglion cell development and image stabilizationCurrent Biology
- 94.Modular transcriptional programs separately define axon and dendrite connectivityeLife 8
- 95.Brain wiring determinants uncovered by integrating connectomes and transcriptomesCurrent Biology 33:3998–4005
- 96.Chapter 4 Transcriptional Networks in the Early Development of SensoryMotor CircuitsElsevier
- 97.The emergence of transcriptional identity in somatosensory neuronsNature 577
- 98.Sensory neuron diversity in the inner ear is shaped by activityCell 174:1229–1246
- 99.Snapshot: Spinal cord developmentCell 146:178–178
- 100.The molecular biology of axon guidanceScience 274
- 101.Developmental mechanisms that generate precise patterns of neuronal connectivityCell 72
- 102.Teneurin trans-axonal signaling prunes topographically missorted axonsCell Reports 42
- 103.The order and place of neuronal differentiation establish the topography of sensory projections and the entry points within the hindbrainJournal of Neuroscience 35:7475–7486
- 104.Distribution of neurosensory progenitor pools during inner ear morphogenesis unveiled by cell lineage reconstructioneLife 6
- 105.Dii reveals a prenatal arrival of efferents at the differentiating otocyst of miceHearing Research 65
- 106.Development of the Cranial Nerve Ganglia and Related Nuclei in the RatSpringer Berlin Heidelberg
- 107.Identification of early neurons in the brainstem and spinal cord: I. an autoradiographic study in the chickJournal of Comparative Neurology 192
- 108.Development of vestibular afferent projections into the hindbrain and their central targetsBrain Research Bulletin 60
- 109.Regional patterning of reticulospinal and vestibulospinal neurons in the hindbrain of mouse and rat embryosThe Journal of Comparative Neurology 411
- 110.Development of the brain stem in the rat. V. Thymidine-radiographic study of the time of origin of neurons in the midbrain tegmentumThe Journal of Comparative Neurology 198:677–716
- 111.Generation of the ocular motor nuclei and their cell types in the rabbitThe Journal of Comparative Neurology 200:69–82
- 112.Do oculomotor neuroblasts migrate across the midline in the fetal rat brain?Anatomy and Embryology 150
- 113.Constructing the hindbrain: Insights from the zebrafishDevelopmental Dynamics 224:1–17
- 114.A structural and functional ground plan for neurons in the hindbrain of zebrafishProceedings of the National Academy of Sciences 108:1164–1169
- 115.A single-cell transcriptome atlas for zebrafish developmentDevelopmental Biology 459:100–108
- 116.A single cell transcriptome atlas of the developing zebrafish hindbrainDevelopment
- 117.Axon guidance in the developing ocular motor system and duane retraction syndrome depends on semaphorin signaling via alpha2-chimaerinProceedings of the National Academy of Sciences 109
- 118.Patterning and axon guidance of cranial motor neuronsNature Reviews Neuroscience 8
- 119.The neuronal architecture and topography of the nucleus vestibularis tangentialis in the late chick embryoNeuroscience 2:189–IN2
- 120.Translating basic research to astronaut health in space: Nasa ames rodent specimen biobanking for the human research programGravitational and Space Research 4
- 121.Development of utricular otoliths, but not saccular otoliths, is necessary for vestibular function and survival in zebrafishJournal of Neurobiology 43:329–337
- 122.Delayed otolith development does not impair vestibular circuit formation in zebrafishJournal of the Association for Research in Otolaryngology 18:415–425
- 123.Development of functional hindbrain oculomotor circuitry independent of both vascularization and neuronal activity in larval zebrafishCurrent Neurobiology 7:62–73
- 124.Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682
- 125.In vivo birthdating by BAPTISM reveals that trigeminal sensory neuron diversity depends on early neurogenesisDevelopment 135:3259–3269
- 126.Vagus motor neuron topographic map determined by parallel mechanisms of hox5 expression and time of axon initiationCurrent Biology 27:3812–3825
- 127.A protocol for whole-mount immuno-coupled hybridization chain reaction (wichcr) in zebrafish embryos and larvaeSTAR Protocols 2
- 128.Segment number threshold determines juvenile onset of germline cluster expansion in platynereis dumeriliiJournal of Experimental Zoology Part B: Molecular and Developmental Evolution 338
- 129.Moderated estimation of fold change and dispersion for rna-seq data with deseq2Genome Biology 15
- 130.Requirement for lim homeobox gene isl1 in motor neuron generation reveals a motor neuron dependent step in interneuron differentiationCell 84
- 131.Molecular heterogeneity among primary motoneurons and within myotomes revealed by the differential mrna expression of novel islet-1 homologs in embryonic zebrafishDevelopmental Biology 171
- 132.Zebrafish hox genes: expression in the hindbrain region of wild-type and mutants of the segmentation gene, valentinoDevelopment 125
- 133.Integrated analysis of multimodal single-cell dataCell 184:3573–3587
- 134.Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracyStatistical Science 1:54–75
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Copyright
© 2024, Goldblatt et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 506
- downloads
- 26
- citations
- 2
Views, downloads and citations are aggregated across all versions of this paper published by eLife.