Clonally related, Notch-differentiated spinal neurons integrate into distinct circuits
Abstract
Shared lineage has diverse effects on patterns of neuronal connectivity. In mammalian cortex, excitatory sister neurons assemble into shared microcircuits. In Drosophila, in contrast, sister neurons with different levels of Notch expression (NotchON/NotchOFF) develop distinct identities and diverge into separate circuits. Notch-differentiated sister neurons have been observed in vertebrate spinal cord and cerebellum, but whether they integrate into shared or distinct circuits remains unknown. Here, we evaluate how sister V2a (NotchOFF)/V2b (NotchON) neurons in the zebrafish integrate into spinal circuits. Using an in vivo labeling approach, we identified pairs of sister V2a/b neurons born from individual Vsx1+ progenitors and observed that they have somata in close proximity to each other and similar axonal trajectories. However, paired whole-cell electrophysiology and optogenetics revealed that sister V2a/b neurons receive input from distinct presynaptic sources, do not communicate with each other, and connect to largely distinct targets. These results resemble the divergent connectivity in Drosophila and represent the first evidence of Notch-differentiated circuit integration in a vertebrate system.
Editor's evaluation
This is an important article that describes the connectivity of sibling neurons in the zebrafish spinal cord, where one sibling receives Notch signaling (Notch-ON) and the other does not (Notch-OFF). They find that V2a and V2b siblings have different morphology, inputs, outputs, and fail to connect to each other; this provides new insight into the role of lineage in specifying neuronal connectivity. The experiments are convincing and the conclusions are supported by the data presented.
https://doi.org/10.7554/eLife.83680.sa0eLife digest
The brain is populated by neurons which are generated during embryonic development from cells called progenitors. Neurons that come from the same progenitor cell are considered to be ‘sisters’. In certain brain regions of mice, sister neurons are often wired into shared networks, meaning they are more likely to receive input from the same neurons and connect with each other than non-sister cells.
In contrast, in invertebrate animals, like the fruit fly, sister neurons often have different identities and are less likely to connect with each other. This may be because sister neurons in fruit flies often have varied levels of a protein called Notch, which plays an important role in establishing the identity of cells. Vertebrate and invertebrate animals are different in many respects, and it remained unclear whether Notch levels dictate which sister neurons connect together in vertebrates as they do in fruit flies.
To investigate, Bello-Rojas and Bagnall studied two neurons in the spinal cord of zebrafish embryos which come from the same type of progenitor cell: the V2a neuron which has low levels of Notch, and the V2b neuron which has high levels of Notch. Fish, like humans, are vertebrates; however, their embryos are mostly transparent, making it easier to track how their neurons make connections during development using a microscope. This enabled Bello-Rojas and Bagnall to monitor whether V2a and V2b sister neurons joined the same network, like in other vertebrates, or different networks, akin to sister neurons in fruit flies which also have differing levels of Notch.
Bello-Rojas and Bagnall found that sister V2a and V2b neurons stayed close to one another and seemed to connect through similar paths. However, closer investigation revealed that the sister neurons did not receive input from the same source. They also did not connect to each other or the same output neuron, suggesting that V2a and V2b sister neurons are part of different networks.
This is the first time Notch levels have been shown to regulate which network a neuron will join in a vertebrate species. Since the V2a and V2b neurons are involved in controlling body movement, future work should determine whether adding progenitor cells that produce these neurons into the spinal cord could help the neuron network recover after injury or disease.
Introduction
How does shared lineage affect neuronal circuitry? Neurons arising from common progenitors are more likely to exhibit stereotypic patterns of connectivity, in two models from vertebrate and invertebrate systems. In mouse cortex, clonally related excitatory sister neurons preferentially form connections within a shared microcircuit (Xu et al., 2014; Yu et al., 2009). In contrast, clonally related sister neurons in Drosophila form distinct NotchON and NotchOFF hemilineages which innervate distinct targets and often express different neurotransmitters (Artavanis-Tsakonas et al., 1999; Endo et al., 2007; Harris et al., 2015; Lacin et al., 2019; Lacin and Truman, 2016; Mark et al., 2021; Pinto-Teixeira and Desplan, 2014; Skeath and Doe, 1998).
Notch-differentiated clonally related sister neurons have been observed in the vertebrate spinal cord and cerebellum (Kimura et al., 2008; Peng et al., 2007; Zhang et al., 2021), but it remains unknown whether these clonally related neurons integrate into shared circuits. In ventral spinal cord, motor neurons and interneurons develop from five progenitor domains (p0, p1, p2, pMN, and p3) (Goulding, 2009; Goulding and Lamar, 2000; Jessell, 2000). Progenitors in the p2 domain transiently express the transcription factor Vsx1 (Kimura et al., 2008; Passini et al., 1998). Each p2 progenitor makes a final paired division into an excitatory V2a (NotchOFF) and an inhibitory V2b (NotchON) neuron via asymmetric expression of Delta ligands and subsequent Notch-mediated lateral inhibition (Del Barrio et al., 2007; Francius et al., 2016; Kimura et al., 2008; Okigawa et al., 2014; Peng et al., 2007).
Although both V2a and V2b neurons project axons ipsilaterally and caudally, these neuron classes differ in other aspects. V2a interneurons express vsx2 (referred to as chx10 in this paper for clarity) and provide glutamatergic drive onto motor populations (Kimura et al., 2006), whereas V2b interneurons express gata3 and provide glycinergic and GABAergic inhibition onto motor populations (Andrzejczuk et al., 2018; Callahan et al., 2019). V2b neurons also support flexor/extensor alternation through reciprocal inhibition in limb circuits (Britz et al., 2015; Zhang et al., 2014). Given their shared origin but divergent cellular identities, it remains unknown whether these V2a/b sister neurons integrate into shared or distinct functional spinal circuits.
We investigated whether V2a/b sister neurons in zebrafish spinal cord preferentially integrate into shared circuits, as with clonally related cortical neurons, or distinct circuits, as with Notch-differentiated hemilineages in Drosophila. Using a sparse labeling approach, we directly observed and identified individual pairs of sister V2a/b neurons arising from a single progenitor. Our morphological and electrophysiological analyses reveal that although sister V2a/b neurons share anatomical characteristics, these sister neurons diverge into separate circuits, with largely distinct presynaptic and postsynaptic partners. To the best of our knowledge, this is the first assessment of circuit integration of Notch-differentiated clonally related neurons in vertebrate models.
Results
Micro-injection of vsx1 plasmid allows for clonal pair tracking in vivo
In both zebrafish and mice, vsx1+ progenitors give rise to two distinct daughter populations, V2a and V2b neurons (Kimura et al., 2008; Peng et al., 2007). Using transgenic zebrafish, individual vsx1+ progenitors have been shown to undergo a final paired division into one V2a (NotchOFF) and one V2b neuron (NotchON) (Figure 1A; Kimura et al., 2008). We aimed to develop a protocol to label and identify individual clonal pairs resulting from this division in vivo. To label individual pairs, we micro-injected titrated amounts of a bacterial artificial chromosome (BAC) construct, vsx1:GFP, into fertilized zebrafish embryos at the single-cell stage (Figure 1B). At the 21-somite stage, larval zebrafish were screened for vsx1 GFP+ progenitors and then imaged every 5 min to capture the progenitor division (Figure 1C). Progenitors became elongated before dividing into two distinct cells.
When the fish become free swimming at 4 days post fertilization (dpf), vsx1 GFP+ pairs were assessed for co-expression of known V2a/b transcription factors (chx10/gata3) to verify their neuronal identities, using transgenic fish Tg(chx10:lox-dsRed-lox:GFP) (Kimura et al., 2006) or Tg(gata3:lox-dsRed-lox:GFP) (Callahan et al., 2019; Figure 1D). For simplicity, these fish lines will be referred to as chx10:Red and gata3:Red. We assayed these in separate experiments due to overlap in fluorescence from reporter lines. Figure 1D presents example images of vsx1 GFP+ pairs in which one of the two neurons in the pair expresses the appropriate marker: a clonal pair (green) where one neuron co-expresses the V2a marker Chx10 (left), and a different clonal pair in which one neuron co-expresses the V2b marker Gata3 (middle). Based on previous work, we expected that every vsx1 GFP+ pair would consist of one V2a and one V2b neuron (Kimura et al., 2008). However, among clonal pairs imaged in the chx10:Red background, only 61/92 (66.3%) of vsx1 GFP+ pairs included one identified V2a neuron (Figure 1E). In contrast, in the gata3:Red background, 35/38 (92.1%) of vsx1 GFP+ pairs included one identified V2b neuron (Figure 1E). Rarely, both vsx1 GFP+ neurons in a pair expressed both Chx10 or Gata3 markers (<10%). However, in 25% of vsx1 GFP+ pairs in chx10:Red fish, neither neuron expressed the V2a marker.
A possible explanation for the lower rate of V2a marker expression could be under-labeling in the fluorescent reporter line. Alternatively, there is at least one additional population of neurons to emerge from the V2 domain. In zebrafish, the V2s population is glycinergic and expresses Sox1a (Gerber et al., 2019). V2s neurons resemble V2c neurons in mice in that they both express Sox1a and arise after V2a/b development; however, V2c neurons are GABAergic while V2s neurons are purely glycinergic (Gerber et al., 2019; Panayi et al., 2010). Using the Tg(sox1a:dmrt3a-gata2a:EFP(ka705)) reporter line, here referred to as sox1:GFP, (Gerber et al., 2019), we assayed the presence of sox1+/vsx1+ neurons by injecting a vsx1:mCherry BAC in embryos at the single-cell stage. In 16/65 (24.6%) of vsx1 mCherry+ pairs, one of the two sister neurons co-labeled with sox1a (Figure 1D, right), and in 49/65 (75.4%) of vsx1 mCherry+ pairs, neither neuron co-labeled with sox1a (Figure 1E). These results suggest that not all vsx1+ progenitors differentiate into V2a/b pairs. Instead, approximately 75% of vsx1 progenitors divide into V2a/b pairs while the remainder divide into V2b/s pairs. We did not see any vsx1+ triplets or singlets in co-label experiments (0/195), suggesting that vsx1 progenitors only undergo a single, terminal division. Based on these results, we conclude that our stochastic labeling approach successfully labeled clonally related V2 neurons, but required a chx10 co-label to properly identify vsx1 pairs as V2a/b neurons in vivo.
Sister V2a/b neurons remain in close proximity to each other
Immediately after progenitor division around 1 dpf, sister V2a/b neurons are located in close proximity to each other (Kimura et al., 2008), but they have not been followed out to 3–5 dpf when the spinal circuit transitions from spontaneous coiling during embryonic stages to the beat-and-glide locomotion at the larval stage. To assess somatic relationships between sister V2 neurons at larval stages, we measured inter-soma distances among sister and non-sister pairs. In an example fish (Figure 2A), the inter-soma distance between the GFP-only sister neuron (presumed V2b) to the GFP/Chx10:Red co-labeled V2a neuron (dark red arrow) was shorter than the distance between non-sister neurons (white arrows). Across animals, sister V2 neurons were usually closer neighbors than non-sister V2 neurons, whether measured in the V2a or V2b reporter lines (Figure 2B). Sister vsx1+ neurons remained in close proximity to each other throughout embryonic and larval development. Beginning at 24 hpf, we embedded embryos in low melting point agarose, imaged, and then re-imaged at 48 hpf. At 24 hpf, sister neuron centers were ~7 µm apart, or effectively adjacent. By 48 hpf, this inter-soma distance increased slightly to ~9 µm (Figure 2C). In a separate set of experiments, we tracked vsx1+ sister neurons from 48 to 96 hpf by embedding fish for imaging at 48 hpf, freeing from agarose after imaging, and re-embedding at 96 hpf. The distance between somata increased slightly, but still remained relatively short (Figure 2C). Because V2a/b somata are ~10 µm in size (Callahan et al., 2019; Kimura et al., 2006; Menelaou et al., 2014), our data suggest that sister V2 neurons usually remain adjacent. Lastly, restricting our analysis to sister V2a/b neurons using chx10:Red fish, we found that V2b neurons were typically positioned more dorsally than their sister V2a counterparts (Figure 2D), consistent with previous work showing inhibitory populations are located more dorsally than excitatory neurons in spinal cord (Kimura et al., 2006; McLean et al., 2007). Altogether, our data demonstrate that sister V2a/b neurons develop and remain close to each other during larval stages. As a result, in subsequent experiments we inferred that sparsely labeled vsx1 GFP+ neurons located close to each other at 3–4 dpf represented sister pairs.
Though V2a axons are consistently longer, sister V2a/b axons travel along similar trajectories
As V2a/b neurons both project descending, ipsilateral axons, we next assessed whether the axons of clonally related V2a/b neurons exhibited any consistent morphological characteristics. Vsx1 GFP+ pairs were labeled in chx10:Red fish using a vsx1:GFP plasmid and later imaged on a confocal microscope. V2a/b axons were reconstructed (Figure 3A), and the descending axon length of each clonal V2a/b neuron was measured. Sister V2a neurons exhibited axons that were on average 61% longer than their V2b counterparts (Figure 3B), consistent with work showing that Notch expression attenuates axon growth (Mark et al., 2021; Mizoguchi et al., 2020). There was no relationship between the length of the axons and their location along the rostral-caudal axis of the fish (Figure 3B). To measure axon proximity, the shortest distance between the V2b and V2a axon was calculated along each point of the V2b neuron, beginning at the axon hillock (Figure 3C, inset). The fraction of those inter-axon distances within 5 µm was calculated for each pair. Indeed, clonally related V2a/b neurons send axons along a similar trajectory, with a median of 37.1% of the V2b axon length within 5 µm of the V2a axon (Figure 3C). Because the axons follow similar paths, these results suggest a possibility for sister V2a/b neurons to contact shared synaptic targets.
Sister V2a/b neurons receive input from distinct synaptic circuits
Work in hippocampus has shown that sister neurons are more likely to receive synaptic input from shared presynaptic partners than non-sister neurons (Xu et al., 2014). In contrast, NotchON and NotchOFF sister neurons in Drosophila integrate into separate hemilineages that segregate spatially, although whether they receive shared input is not known (Harris et al., 2015). To evaluate whether sister V2a/b neurons receive input from shared or distinct presynaptic partners in vivo, we performed paired whole-cell electrophysiology in voltage clamp from clonally related pairs of V2a/b neurons identified as above (Figure 4A and B; Bagnall and McLean, 2014). Both sister neurons were held at –80 mV, the chloride reversal potential, to isolate excitatory postsynaptic currents (EPSCs), while a bright-field stimulus was used to elicit fictive swim (Figure 4C). The timing of EPSCs arriving in each neuron of the pair was asynchronous, as exemplified by an overlay of several hundred EPSCs from either a V2a/b and the associated EPSC-triggered average in its sister neuron (Figure 4D and E). A summary of the amplitudes of detected EPSCs and associated EPSC-triggered averages for this example neuron is shown in Figure 4F. Across recordings from 13 clonally related pairs in vivo, we consistently saw little to no synchronous synaptic input (Figure 4G).
Lastly, we wanted to compare whether this asynchrony in EPSC input was present in non-sister V2a/b neurons from the same segment. Using the same analysis, it appeared that non-sister V2a/b neurons from the same spinal segment receive input from distinct synaptic sources as well (Figure 4H). The asynchronous timing of these inputs suggests that they cannot be arriving from a shared presynaptic source, but rather, different presynaptic sources which fire at different times (Bagnall and McLean, 2014). Altogether, these data show that not only sister V2a/b neurons, but any V2a-V2b pair, clonal or non-clonal, from the same segment receives input from distinct presynaptic sources during light-evoked locomotion at slow to medium locomotor speeds. We cannot rule out the possibility that sister neurons receive shared inputs from circuits for fast locomotion or specialized behaviors.
Sister V2a/b neurons do not form synaptic connections with each other
Clonal pair analysis in cortex has shown that sister neurons preferentially form synapses onto each other (Yu et al., 2009; Zhang et al., 2017). To identify whether sister V2a/b neurons form synaptic connections with each other, paired in vivo whole-cell recordings were performed in chx10:Red fish as described above (Figure 5A). Spiking was elicited by depolarizing current steps in either the V2a or the V2b neuron while the other neuron was held in voltage clamp to measure synaptic responses (Vhold of –80 mV in V2b neurons to measure EPSCs, Vhold of 0 mV in V2a neurons to measure IPSCs). In both cases, there were no detectable evoked currents, showing that sister V2a/b neurons do not connect with each other (Figure 5C). Similarly, non-clonally related V2a/b neurons exhibited no interconnectivity (Figure 5D). Therefore, V2a/b neurons in the same segment do not form direct synapses with each other. Any connectivity among V2a and V2b neurons likely occurs between neurons in different segments (Sengupta and Bagnall, 2022).
Sister V2a/b neurons provide asymmetric input onto downstream neurons in spinal cord
Research in cortex has shown that clonally related inhibitory interneurons form synaptic connections with shared downstream targets (Zhang et al., 2017), although this claim is disputed (Mayer et al., 2016). Given the proximity of sister V2a/b axons (Figure 3), it was plausible that they share common downstream targets. To address this question, we micro-injected a vsx1:Gal4 BAC and a UAS:CoChR2-tdTomato plasmid (Antinucci et al., 2020; Schild and Glauser, 2015) in embryos at the single-cell stage to drive stochastic expression of this channelrhodopsin variant in vsx1 sister neurons (Figure 6A) for selective optical stimulation of sister neurons. To distinguish presumed V2a/b pairs from V2b/s pairs, we screened CoChR2-tdTomato+ vsx1 sister pairs for morphological characteristics (‘Methods’). We validated that the optical stimuli effectively elicited spiking in CoChR2-tdTomato+ vsx1 sister neurons by performing cell-attached recordings while providing a 10 ms light pulse (Figure 6B). All CoChR2-tdTomato+ vsx1 neurons fired action potentials in response to optical stimulation (Figure 6C, n = 13 neurons from 12 fish). Spiking was elicited in both V2a and V2b CoChR2-tdTomato+ neurons (Figure 6D). Similar experiments were performed on nearby CochR2-tdTomato(-) neurons to ensure that optical stimuli evoked spiking only in neurons expressing CoChR2-tdTomato. All CoChR2-tdTomato(-) neurons remained inactive during the light stimulus (Figure 6C, n = 22 neurons from 18 fish). In CoChR2-tdTomato+ neurons, most light-evoked spikes were observed throughout the duration of the stimulus, with some spiking following the stimulus window (Figure 6C and E). This prolonged activity is most likely due to the long inactivation kinetics of the CoChR2 variant (Antinucci et al., 2020). Altogether, our optogenetic approach is a feasible method for assessing downstream connectivity of sister V2a/b neurons.
Having validated our optogenetic approach, we proceeded to perform whole-cell patch-clamp recordings on known V2a/b neuron downstream spinal targets (i.e., motor neurons, V1, V2a, V2b neurons), which were located 1–4 segments caudal to the V2a/b sister pair in voltage-clamp mode using a cesium-based internal solution (Figure 7A; Bagnall and McLean, 2014; Callahan et al., 2019; Kimura et al., 2006; Menelaou and McLean, 2019). Because sister vsx1+ neurons are close to each other, our optogenetic stimulus would activate both neurons simultaneously. However, by clamping the target neuron at different reversal potentials, we could isolate either evoked EPSCs or inhibitory postsynaptic currents (IPSCs) (Figure 7B–D). In most recorded neurons, optical stimuli evoked neither EPSCs nor IPSCs, consistent with sparse connectivity in the spinal cord (n = 85/99; Figure 7E and F). In six target neurons, we recorded evoked EPSCs (Vhold –80 mV) but not IPSCs, demonstrating that the target neuron received synaptic input from the CoChR2-labeled V2a neuron but not the V2b (Figure 7B and E). In another six target neurons, we detected evoked IPSCs (Vhold 0 mV) but no EPSCs, demonstrating that the target neuron received synaptic input from the V2b but not the V2a neuron (Figure 7C and E). In a subset of experiments, NBQX/APV or strychnine were used to block responses and confirm glutamatergic or glycinergic connections, respectively (Figure 7B, n = 2; Figure 7C, n = 4). In two instances, a target neuron received both evoked EPSCs and IPSCs, with the magnitude of IPSCs approximately fivefold larger than the magnitude of the EPSCs, suggesting an asymmetric connection from sister V2a/b neurons (Figure 7E).
In 10 neurons, we detected a slow depolarizing current when target neurons were held at –80 mV (Figure 7D and F, gray), but not at 0 mV. This evoked current had a lower amplitude and longer rise time than fast evoked EPSCs (Figure 7G). This slow excitatory current may be caused by a weak disynaptic electrical connection (Menelaou and McLean, 2019), but we were not able to eliminate it with gap junction blockers (carbenoxylone and 18-β-glycyrrhetinic acid). We summarize the identities of target neurons receiving synaptic input from sister V2a/b neurons in Table 1. Target neurons were evenly divided between motor neurons (early and late born), and excitatory and inhibitory interneurons. Overall, these results demonstrate that clonally related V2a/b neurons do not preferentially form synaptic connections with shared targets.
Discussion
In this study, we showed that clonally related V2a/b neurons exhibit similar morphological characteristics, but form synapses with and receive information from largely distinct neuronal partners. Through our use of plasmid injections and time-lapse imaging, we definitively identified individual pairs of clonally related V2a/b neurons born from a single vsx1+ progenitor cell in vivo (Figure 7H). Additionally, some vsx1+ progenitors appear to divide into V2b/s pairs. Within V2a/b pairs, we saw that sister neuron somata remain in close proximity to each other and send their axons along similar trajectories. However, our electrophysiological data showed that these sister neurons integrate into distinct circuits. Clonally related V2a/b neurons do not communicate with each other, do not receive input from similar sources, and infrequently connect to the same downstream target. This connectivity pattern resembles circuitry seen in Drosophila Notch-differentiated hemilineages (Figure 7I). Our results represent the first evidence of Notch-differentiated circuit integration in a vertebrate system, and may reflect a means of cell-type and circuit diversification in earlier evolved neural structures.
Notch determines cellular identity of vsx1+ sister neurons
Notch is an important regulator in V2a/b differentiation, and during vsx1+ progenitor division, differences in Notch expression result in the onset of V2a (NotchOFF) or V2b (NotchON) programming (Batista et al., 2008; Debrulle et al., 2020; Kimura et al., 2008; Mizoguchi et al., 2020; Okigawa et al., 2014). However, it remains unknown whether Notch plays a role in sister V2a/b development beyond initiating cellular identity or whether it functions as an intermediary step before other molecular factors determine cellular morphology postmitotically (Kozak et al., 2020; Mizoguchi et al., 2020). Our morphological analysis showed differences in V2a/b axon lengths and dorsoventral position (Figures 2 and 3). Further experiments are needed to evaluate whether these differences are a result of Notch signaling or intrinsic to postmitotic cellular identity. Constitutive manipulation of Notch levels results in a skewing of V2a/b numbers (Mizoguchi et al., 2020). Therefore, a temporally controlled manipulation is needed to address the role of Notch in postmitotic morphological and functional development of sister V2a/b neurons.
Similarly, the recently discovered V2s population relies on Notch signaling for its development with Notch KO mutants showing a decrease in sox1a+ neurons (Gerber et al., 2019). We speculate that some vsx1+ progenitors give rise to V2b/s sister pairs in addition to the previously described V2a/b pairs. Our experiments in reporter lines (Figure 1) suggested that ~75% of vsx1 GFP+ progenitors divided into V2a/b sister neurons, whereas ~25% resulted in V2b/s neuron pairs. V2s neurons arise later than the initial wave of V2a/b pairs (Gerber et al., 2019). Because Notch has been shown to exhibit different effects on cellular identity during early and late development, we suggest that delayed Notch activity causes some later born vsx1+ sister neurons to adopt a V2b/s identity, which are both NotchON (Jacobs et al., 2022). Similarly, only early cerebellar progenitors appear to undergo Notch differentiation into distinct cell types, the Purkinje and granule cells (Zhang et al., 2021). Notch overexpression experiments could have biased differentiation in favor of V2b/s pairs earlier in development, accounting for the increase in V2b and decrease in V2a numbers (Mizoguchi et al., 2020). However, these experiments have not looked at changes to V2s numbers, so selective evaluation of later born V2 progenitors is needed to identify whether V2b/s clonal pairs exist and if so whether they are temporally delayed relative to V2a/b pairs.
Notch-mediated development influences circuit formation
Lineage pathfinding and innervation differences in Drosophila are well documented, and Notch-differentiated sister neurons in these organisms develop different axon trajectories, presumably connecting to different downstream targets (Harris et al., 2015; Truman et al., 2010). Our data show that vsx1+ sister neurons in spinal cord have similar descending trajectories albeit different axon lengths (Figure 3). Analysis of Notch-differentiated lineages in vertebrate cerebellum have shown that Notch mediates cerebellar progenitor differentiation into excitatory and inhibitory cerebellar cell types (Zhang et al., 2021), but it is not yet known whether the resulting neurons integrate into shared or distinct circuits. Our results are consistent with a framework in which the progeny of Notch-differentiated divisions preferentially integrate into distinct networks in both invertebrates and vertebrates. In contrast with cortical lineages, the divergent cellular identities of sister V2a/b neurons appear to determine that they participate in distinct circuits. We speculate that earlier evolved neural structures rely on Notch-differentiated divisions as a means to diversify neuronal populations during development. The presence of Notch-differentiated sister neurons in both cerebellum and spinal cord could represent an efficient mechanism to generate diverse cell types early in development, in contrast to cortical reliance on dedicated streams of excitatory and inhibitory neural progenitors (Goulding, 2009; Leto et al., 2016; Ma et al., 2018). This would allow for the development of several neuronal cell types, each governed by their own intrinsic molecular cues.
Shared vsx1+ progenitor birthdates do not lead to shared integration
Developmental timing allows for proper integration of neurons into functional speed-dependent locomotor circuits. In zebrafish, motor neurons and interneurons born during similar developmental windows are active and recruited at similar speeds (McLean et al., 2007; McLean and Fetcho, 2009). These birthdate-dependent microcircuits emerge in larvae and persist into adulthood (Ampatzis et al., 2014). By definition, vsx1+ sister neurons share a birthdate, suggesting that both neurons are likely recruited at similar speeds and therefore might integrate into shared microcircuits. However, our work shows that vsx1 sister neurons neither synapse onto each other, receive similar inputs, nor frequently target the same neurons. One possible explanation is that sister V2a/b divergence in cellular identity may cause integration into different hemilineage temporal cohorts, similar to Drosophila, which then determine their neuronal connectivity (Mark et al., 2021). Additionally, V2b neurons, whose recruitment patterns have not yet been described, may participate in different behaviors than V2a neurons. This separation of pathways driving excitatory and inhibitory neurons would allow for independent activation (accelerator) or inactivation (brake) of movement (Callahan et al., 2019; Eklöf-Ljunggren et al., 2012). It is worth noting that we measured synaptic inputs during fictive locomotion induced by bright-field stimuli, and that the possibility remains sister vsx1 neurons do receive similar inputs under different behavioral paradigms, such as turns or escapes.
Lastly, the sister vsx1 neurons infrequently connected to the same downstream targets (Figure 7). Because we saw two examples of targets receiving input from both the V2a and V2b neuron of a clonal pair, it is unclear whether sister neurons are explicitly discouraged from sharing downstream targets or whether it is simply random. In either case, the observed connectivity divergence might function to coordinate antagonistic components during locomotion. Spinal V1 interneurons target different populations of neurons along the rostral-caudal length of the spinal cord (Sengupta et al., 2021). Even if non-clonally related V2a and V2b neurons generally form synaptic contacts onto the same populations, such as motor neurons, they may exhibit different connectivity patterns in the longitudinal axis, preventing clonally related pairs from sharing downstream targets. Mapping the rostrocaudal connectivity of V2a and V2b populations would address this hypothesis.
Materials and methods
Experimental model and subject details
Request a detailed protocolAll fish used for experiments were 0 to 6 dpf, before the onset of sexual maturation. All experiments and procedures were approved by the Animal Studies Committee at Washington University and adhere to NIH guidelines. Adult zebrafish (Danio rerio) were maintained at 28.5°C with a 14:10 light:dark cycle in the Washington University Zebrafish Facility up to 1 year following standard care procedures. Larval zebrafish used for experiments were kept in Petri dishes in system water or housed with system water flow.
To target V2a and V2b neurons, the Tg(chx10:loxP-dsRed-loxP:GFP) (Kimura et al., 2006) (ZDB-ALT-061204-4) and Tg(gata3:loxP-dsRed-loxP:GFP) (Callahan et al., 2019) (ZDB-ALT-190724-4) lines were used. We visualized V2s neurons in Tg(sox1a:dmrt3a-gata2a:EFP(ka705)) (ZDB-ALT-191113-2) (Gerber et al., 2019), a gift from Dr. Uwe Strähle.
Stochastic single-cell labeling by microinjections
Request a detailed protocolTg(chx10:loxP-dsRed-loxP:GFP) and Tg(gata3:loxP-dsRed-loxP:GFP) were injected with a vsx1:GFP bacterial artificial chromosome (BAC) at a final concentration of 5 ng/µL (a gift from Dr. Shin-ichi Higashijima). Tg(sox1a:dmrt3a-gata2a:EFP(ka705)) were injected with a vsx1:mCherry BAC at 15 ng/µL (generated from BAC CH211-67N1 by VectorBuilder, Inc). To label clonal pairs with an optogenetic activator, wild-type embryos were injected with a vsx1:Gal4 BAC (VectorBuilder) and UAS:CoChR2-tdTomato plasmid (Addgene Cat# 124233) at 20 ng/µL and 25 ng/µL, respectively. The embryos were transferred to system water to develop. Embryos were screened between 1 and 4 dpf for sparse expression of Red/GFP fluorophores and selected for confocal imaging and electrophysiology. In optogenetic experiments, it was not feasible to use Tg(chx10:Red) animals to identify V2a neurons because of fluorophore overlap. Instead, we screened CoChR2-tdTomato+ vsx1 sister neurons for distinguishable V2a morphology, specifically the presence of an ascending collateral, which is characteristic of V2a but not V2b or V2s neurons (Callahan et al., 2019; Gerber et al., 2019; Menelaou et al., 2014).
Confocal imaging
Request a detailed protocolAt 18–24 hpf, larvae were anesthetized in 0.02% MS-222 and embedded in low-melting point agarose (0.7%) in a 10 mm FluoroDish (WPI). Spinal segments with sparsely labeled progenitors were imaged with a time-lapse approach, consisting of one Z-stack every 5 min, under a spinning disk confocal microscope (Crest X-Light V2; laser line 470 nm; upright Scientifica microscope; 40× objective; imaged with Photometrics BSI Prime camera). After progenitor division, larvae were kept in the FluoroDish inside of an incubator and reimaged at a higher-resolution at 2 dpf with a laser confocal (Olympus FV1200, 488 nm laser, XLUMPlanFl-20× W/0.95 NA water immersion objective).
Larvae imaged beginning at 2 dpf were anesthetized in 0.02% MS-222 and embedded in low melting point agarose (1.5%) in a 10 mm FluoroDish (WPI). Images were acquired on an Olympus FV1200 confocal microscope as above. A transmitted light image was obtained along with laser scanning fluorescent images to identify spinal segments. Sequential scanning used for multi-wavelength images. Fish were unembedded from the agarose and placed separately in labeled Petri dishes and later reimaged at 4 dpf as described above. In some cases, fish were only imaged at 4 dpf using the embedding methods described above. Transcription factor co-expression was evaluated visually.
Image analysis
Request a detailed protocolConfocal images were analyzed using Imaris (9.8, Bitplane) and ImageJ (1.53q, FIJI) (Schindelin et al., 2012). For axon tracing, stitched projection images were made with the Pairwise stitching (Preibisch et al., 2009) ImageJ plugin. The overlap of the fused image was smoothed with linear blending and was registered based on the fill channel or the average of all channels. Three-dimensional (3D) images were reconstructed and analyzed using Imaris. Axon length measurements of each reconstructed neuron were obtained using the Filament function to trace over the 3D rendering. Axon length includes only the descending branches of the neuron, starting at the axon hillock. 3D axon coordinates of descending projections were exported from Imaris, and separation of axon distances was calculated as the shortest distance between sister V2b to sister V2a axons. Muscle segment number was counted under differential interference contrast (DIC). Distances between soma centers were measured in three dimensions using the Points function in Imaris. Normalized dorsoventral soma position was calculated by measuring the distance of the soma from the notochord and dividing by the total height of the spinal cord, with 0 as the ventral-most point.
Electrophysiological recordings
Request a detailed protocolCell-attached recordings were targeted to stochastically labeled WT fish with vsx1:Gal4 BAC and UAS:CoChR2-tdTomato plasmids to calibrate firing of vsx1:Gal4;UAS:CoChR2-tdTomato vsx1+ pairs. Whole-cell patch-clamp recordings were performed in Tg(chx10:loxP-dsRed-loxP:GFP) injected with vsx1:GFP and Tg(chx10:GFP;gata3:loxP-dsRed-loxP:GFP) larvae at 4–6 dpf for paired clonal V2a/b and non-clonal V2a/b recordings, respectively. Additionally, whole-cell patch-clamp recordings were performed in stochastically labeled WT fish with vsx1:Gal4 BAC and UAS:CoChR2-tdTomato in downstream targets. Larvae were immobilized with 0.1% α-bungarotoxin and fixed to a Sylgard lined Petri dish with custom-sharpened tungsten pins. Each larva was then transferred to a microscope (Scientifica SliceScope Pro) equipped with infrared differential interface contrast optics, epifluorescence, and immersion objectives (Olympus: 40×, 0.8 NA). One muscle segment overlaying the spinal cord was removed (segments 7–17) using a blunt-end glass electrode and suction (Wen and Brehm, 2010). The bath solution consisted of (in mM) 134 NaCl, 2.9 KCl, 1.2 MgCl2, 10 HEPES, 10 glucose, and 2.1 CaCl2. Osmolarity was adjusted to ~295 mOsm and pH to 7.5.
Patch pipettes (5–15 MΩ) were filled with internal solution for voltage and current clamp and cell-attached composed of (in mM) 125 K gluconate, 2 MgCl2, 4 KCl, 10 HEPES, 10 EGTA, and 4 Na2ATP. Whole-cell optogenetic and some paired recordings were performed using internal solution composed of (in mM) 122 cesium methanesulfonate, 1 tetraethylammonium-Cl, 3 MgCl2, 1 QX-314 Cl, 10 HEPES, 10 EGTA, and 4 Na2ATP. Additionally, Alexa Fluor 647 hydrazide 0.05–0.1 mM or sulforhodamine (0.02%) was included to visualize morphology of recorded cells post hoc. Osmolarity was adjusted to ~285 mOsm and KOH or CsOH was used to bring the pH to 7.5. Patch recordings were made in whole-cell configuration using a Multiclamp 700B, filtered at 10 kHz (current clamp) or 2 kHz (voltage clamp). All recordings were digitized at 100 kHz with a Digidata 1440 (Molecular Devices) and acquired with pClamp 10 (Molecular Devices). The following drugs were bath applied where noted: strychnine (10 µM), NBQX (10 µM), APV (100 µM), 18-beta-glycyrrhetinic acid (150 µM), and carbenoxolone disodium salt (500 µM).
During paired electrophysiology recordings, fictive swimming sometimes occurred spontaneously and in other instances was elicited by white light illumination of the animal. In optogenetic experiments examining channelrhodopsin firing and V2a/b targeting, light stimulation was provided with high-intensity epifluorescent illumination (CoolLED pE-300), 10% intensity with a 40× (0.8 NA) water immersion objective for 10 ms. The objective was positioned over a single spinal segment prior to stimulus delivery.
Electrophysiology data were imported into Igor Pro 6.37 (Wavemetrics) using NeuroMatic (Rothman and Silver, 2018). The detection algorithm was based on the event detection instantiated in the SpAcAn environment for Igor Pro (Rousseau et al., 2012) and as previously described (Bagnall and McLean, 2014). All events detected were additionally screened manually to exclude spurious noise artifacts. EPSCs were analyzed using custom-written codes in Igor and MATLAB. For synchronized input evaluation, recordings were excluded from analysis if we could not induce robust fictive locomotion.
Statistics
Statistical tests were performed using MATLAB (R2018a, MathWorks). Due to the non-normal distribution of physiological results, we used nonparametric statistics and tests for representations and comparisons. Details of statistical tests, p-values used, and sample sizes are described in the corresponding figure legends.
Data availability
All analyses generated during this study are included in the manuscript source data file. Computational analyses were performed using code available at https://github.com/bagnall-lab/Event-detection, (copy archived at swh:1:rev:b2fbdb9c9acb263116c2f041caf89d28d1ed1c5f) (Bagnall, 2022). Raw data will be available on Dryad at the following DOI: https://doi.org/10.5061/dryad.4xgxd25dh.
-
Dryad Digital RepositoryClonally related, Notch-differentiated spinal neurons integrate into distinct circuits EPhys raw data.https://doi.org/10.5061/dryad.4xgxd25dh
References
-
Vsx1 and chx10 paralogs sequentially secure V2 interneuron identity during spinal cord developmentCellular and Molecular Life Sciences 77:4117–4131.https://doi.org/10.1007/s00018-019-03408-7
-
Vsx1 transiently defines an early intermediate V2 interneuron precursor compartment in the mouse developing spinal cordFrontiers in Molecular Neuroscience 9:145.https://doi.org/10.3389/fnmol.2016.00145
-
Neuronal patterning: making stripes in the spinal cordCurrent Biology 10:R565–R568.https://doi.org/10.1016/s0960-9822(00)00615-1
-
Circuits controlling vertebrate locomotion: moving in a new directionNature Reviews. Neuroscience 10:507–518.https://doi.org/10.1038/nrn2608
-
Neuronal specification in the spinal cord: inductive signals and transcriptional codesNature Reviews. Genetics 1:20–29.https://doi.org/10.1038/35049541
-
Consensus paper: cerebellar developmentCerebellum 15:789–828.https://doi.org/10.1007/s12311-015-0724-2
-
Neural lineage tracing in the mammalian brainCurrent Opinion in Neurobiology 50:7–16.https://doi.org/10.1016/j.conb.2017.10.013
-
Differences in the morphology of spinal v2a neurons reflect their recruitment order during swimming in larval zebrafishThe Journal of Comparative Neurology 522:1232–1248.https://doi.org/10.1002/cne.23465
-
Sox1 is required for the specification of a novel p2-derived interneuron subtype in the mouse ventral spinal cordThe Journal of Neuroscience 30:12274–12280.https://doi.org/10.1523/JNEUROSCI.2402-10.2010
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Paired patch clamp recordings from motor-neuron and target skeletal muscle in zebrafishJournal of Visualized Experiments 2:2351.https://doi.org/10.3791/2351
Article and author information
Author details
Funding
National Institute on Deafness and Other Communication Disorders (R01 DC016413)
- Martha W Bagnall
Pew Charitable Trusts
- Martha W Bagnall
McKnight Endowment Fund for Neuroscience
- Martha W Bagnall
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. Rebecca Callahan for initial contribution in experimental design and topic development, Dr. Mohini Sengupta for thoughtful critiques of the paper, Dr. Shin-ichi Higashijima for the vsx1:GFP BAC construct, and Dr. Uwe Strähle for the sox1:GFP fish line. We are grateful to Drs. Andreas Burkhalter and Haluk Lacin for insightful comments on manuscript. We also acknowledge the Washington University Zebrafish Facility for fish care and Washington University Center for Cellular Imaging (WUCCI) for supporting the confocal imaging experiments. This work was supported by funding through the National Institute of Health (NIH) R01 DC016413 (MWB). MWB is a Pew Biomedical Scholar and a McKnight Foundation Scholar.
Ethics
All experiments and procedures were approved by the Animal Studies Committee at Washington University within IACUC protocol 20-0367 and adhere to NIH guidelines.
Copyright
© 2022, Bello-Rojas and Bagnall
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
-
- 862
- views
-
- 130
- downloads
-
- 10
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Time estimation is an essential prerequisite underlying various cognitive functions. Previous studies identified ‘sequential firing’ and ‘activity ramps’ as the primary neuron activity patterns in the medial frontal cortex (mPFC) that could convey information regarding time. However, the relationship between these patterns and the timing behavior has not been fully understood. In this study, we utilized in vivo calcium imaging of mPFC in rats performing a timing task. We observed cells that showed selective activation at trial start, end, or during the timing interval. By aligning long-term time-lapse datasets, we discovered that sequential patterns of time coding were stable over weeks, while cells coding for trial start or end showed constant dynamism. Furthermore, with a novel behavior design that allowed the animal to determine individual trial interval, we were able to demonstrate that real-time adjustment in the sequence procession speed closely tracked the trial-to-trial interval variations. And errors in the rats’ timing behavior can be primarily attributed to the premature ending of the time sequence. Together, our data suggest that sequential activity maybe a stable neural substrate that represents time under physiological conditions. Furthermore, our results imply the existence of a unique cell type in the mPFC that participates in the time-related sequences. Future characterization of this cell type could provide important insights in the neural mechanism of timing and related cognitive functions.
-
- Neuroscience
Granule cells of the cerebellum make up to 175,000 excitatory synapses on a single Purkinje cell, encoding the wide variety of information from the mossy fibre inputs into the cerebellar cortex. The granule cell axon is made of an ascending portion and a long parallel fibre extending at right angles, an architecture suggesting that synapses formed by the two segments of the axon could encode different information. There are controversial indications that ascending axon (AA) and parallel fibre (PF) synapse properties and modalities of plasticity are different. We tested the hypothesis that AA and PF synapses encode different information, and that the association of these distinct inputs to Purkinje cells might be relevant to the circuit and trigger plasticity, similar to the coincident activation of PF and climbing fibre inputs. Here, by recording synaptic currents in Purkinje cells from either proximal or distal granule cells (mostly AA and PF synapses, respectively), we describe a new form of associative plasticity between these two distinct granule cell inputs. We show for the first time that synchronous AA and PF repetitive train stimulation, with inhibition intact, triggers long-term potentiation (LTP) at AA synapses specifically. Furthermore, the timing of the presentation of the two inputs controls the outcome of plasticity and induction requires NMDAR and mGluR1 activation. The long length of the PFs allows us to preferentially activate the two inputs independently, and despite a lack of morphological reconstruction of the connections, these observations reinforce the suggestion that AA and PF synapses have different coding capabilities and plasticity that is associative, enabling effective association of information transmitted via granule cells.