Abstract
In motor systems, higher-order neurons provide commands to lower-level central pattern generators (CPGs) that autonomously produce rhythmic motor patterns. Such hierarchical organization is often thought to be inherent in the anatomical position of the neurons. Here, however, we report that a neuron that is member of a CPG in one species acts as a higher-order neuron in another species. In the nudibranch mollusc, Melibe leonina, swim interneuron 1 (Si1) is in the CPG underlying swimming, firing rhythmic bursts of action potentials as part of the swim motor pattern. We found that its homolog in another nudibranch, Dendronotus iris, serves as a neuromodulatory command neuron for the CPG of a homologous swimming behavior. In Dendronotus, Si1 fired irregularly throughout the swim motor pattern. The burst and spike frequencies of Dendronotus swim CPG neurons correlated with Si1 firing frequency. Si1 activity was both necessary and sufficient for the initiation and maintenance of the swim motor pattern. Each Si1 was electrically coupled to all of the CPG neurons and made monosynaptic excitatory synapses with both Si3s. Si1 also bilaterally potentiated the excitatory synapse from Si3 to Si2. “Virtual neuromodulation” of both Si3-to-Si2 synapses using dynamic clamp combined with depolarization of both Si3s mimicked the effects of Si1 stimulation on the swim motor pattern. Thus, in Dendronotus, Si1 is a command neuron that turns on, maintains, and accelerates the motor pattern through synaptic and neuromodulatory actions, thereby differing from its homolog in Melibe in its functional position in the motor hierarchy.
SIGNIFICANCE STATEMENT Cross-species comparisons of motor system organization can provide fundamental insights into their function and origin. Central pattern generators (CPGs) are lower in the functional hierarchy than the neurons that initiate and modulate their activity. This functional hierarchy is often reflected in neuroanatomical organization. This paper definitively shows that an identified cerebral ganglion neuron that is a member of a CPG underlying swimming in one nudibranch species serves as a command neuron for the same behavior in another species. We describe and test the synaptic and neuromodulatory mechanisms by which the command neuron initiates and accelerates rhythmic motor patterns. Thus, the functional position of neurons in a motor hierarchy can shift from one level to another over evolutionary time.
Introduction
A hierarchical organization implies directionality, with higher-level components issuing commands that are carried out by lower levels (Kupfermann and Weiss, 1978). This organization is often considered to be inherent in the neuroanatomy of motor systems. For example, brainstem neurons provide commands to locomotor central pattern generators (CPGs) in the spinal cord (Dubuc et al., 2008; Roberts et al., 2008; Kiehn, 2016). Similarly, in invertebrates, commands for locomotion and feeding come from neurons outside of the CPG and often in other ganglia (Gillette et al., 1982; Gamkrelidze et al., 1995; Panchin et al., 1995; Frost and Katz, 1996; Stein, 2009; Puhl et al., 2012). Here, we show that the hierarchical position of a neuron differs in two species even though the neuron is involved in the production of the same behavior in those species.
The nudibranch molluscs, Melibe leonina and Dendronotus iris exhibit characteristic swimming behaviors that consist of a series of laterally-directed body flexions in which the whole body alternately bends toward the left and right sides (Sakurai et al., 2011). A phylogenetic analysis indicates that the most recent common ancestor of these species likely swam in this manner, making the swimming behaviors homologous (Goodheart et al., 2015; Sakurai and Katz, 2017).
In addition to the behavioral homology, a bilaterally represented neuron was shown to be homologous. First identified in Melibe, the neuron was named Swim interneuron 1 (Si1) because of its involvement in swim CPG (Thompson and Watson, 2005). It was previously established that Si1 in Dendronotus is homologous to Si1 in Melibe based on its conserved neuronal anatomy (soma location in the cerebral ganglion, proximity to a cluster of serotonergic neurons, unique axon shape; Fig. 1A1,B1), and neurochemistry (FMRFamide immunoreactivity), which uniquely identify this neuron in several nudibranch species (Sakurai et al., 2011; Newcomb et al., 2012; Gunaratne et al., 2017). Homologous neurons are a classic feature of gastropod nervous systems (Croll, 1987; Bulloch and Ridgway, 1995; Katz and Quinlan, 2019).
In Melibe, Si1 is electrically-coupled to the ipsilateral Si2 and the contralateral Si4 in the pedal ganglia (Fig. 1A2; Sakurai et al., 2014). Each neuron inhibits its contralateral counterpart to form the main half-center oscillator. The burst period of the main half-center oscillator is regulated by inhibitory synapses from Si3 (Fig. 1A2; Sakurai et al., 2014). Si1 also synapses onto Si3 on both sides, which plays roles in setting the phase relationship and promoting the transition of Si3 bursting from one side to the other (Sakurai et al., 2014). As it is part of the CPG, tonic depolarization or hyperpolarization of one Si1 interferes with rhythmic activity (Figs. 1A3, 4; Sakurai et al., 2011).
In contrast to Melibe, in Dendronotus, Si1 has no reciprocal inhibition and is not a member of the CPG (Fig. 1B2); each Si1 is electrically coupled to both left and right Si2 and fires irregularly throughout the swim motor pattern (Sakurai et al., 2011; Sakurai and Katz, 2016, 2017). The swim CPG is composed of Si2 and Si3, which form the half-center oscillator through reciprocal inhibition of their contralateral counterparts (Fig. 1B2). A homolog of Si4 has not been identified. It was shown that tonic depolarization of one Si1 speeds up bursting, whereas hyperpolarization of Si1 slows it down (Figs. 1B3,4; Sakurai et al., 2011). However, the mechanism of such modulatory action and its roles in the motor pattern generation were not previously elucidated.
Here, we find that unlike its homolog in Melibe, Si1 in Dendronotus is a command neuron for swimming; its spiking activity is necessary for the production of the motor pattern and sufficient to elicit it. Depolarizing Si1 initiates, maintains, and accelerates the swim motor pattern through its synaptic and neuromodulatory actions on the members of the swim CPG. Thus, homologous neurons occupy different positions in the motor hierarchies of Dendronotus and Melibe.
Materials and Methods
Animals.
Adult specimens of Melibe leonina and Dendronotus iris were collected on the west coast of North America by Monterey Abalone Company, Marinus Scientific, and Living Elements. Animals were shipped to Atlanta overnight and kept in artificial seawater (Instant Ocean) at 10–12°C with a 12 h light/dark cycle. Elemental composition of the Instant Ocean is as follows (in mm): 444 Na, 9.5 K, 9.2 Ca, 49.4 Mg, 519 Cl, 26 SO4, 2.3 HCO3, 0.4 H3BO3, 0.3 Br (Creswell, 1993).
Brain preparation.
The animal was first anesthetized by injecting 0.33 m MgCl2 solution into the body cavity. Then, the left body wall near the esophagus was cut open and the brain, consisting of the left and right pedal and fused cerebral-pleural ganglia, was removed, together with a portion of the esophagus. The tissue was placed in a Petri dish lined with Sylgard 184 (Dow Corning) and filled with artificial seawater kept at 4°C. The composition of the artificial sea water was as follows (in mm): 420 NaCl, 10 KCl, 10 CaCl2, 50 MgCl2, 11 d-glucose, and 10 HEPES, pH 7.6. Connective tissue and brain sheath were manually dissected from the dorsal side of the cerebral ganglia and pedal ganglia with forceps and scissors. After desheathing, all nerve roots except for pedal nerve 3 (PdN3) were crushed with forceps to avoid spontaneous movement of any muscle tissue attached to the brain. The desheathed brain was then continuously superfused with artificial seawater at a rate of 0.5–1.0 ml/min at 10°C and allowed to rest for >1 h before starting electrophysiological experiments.
In some experiments, spontaneous neuronal activity was suppressed by applying high divalent cation (Hi-Di) saline, which raises the threshold for neuronal spiking. The composition of the Hi-Di saline was as follows (in mm): 285 NaCl, 10 KCl, 25 CaCl2, 125 MgCl2, 11 d-glucose, and 10 HEPES, pH 7.6. Hi-Di saline was used to examine direct monosynaptic and electrical connections between neurons. For zero calcium experiments, CaCl2 was replaced by MgCl2 to block chemical synaptic transmission.
Electrophysiology.
For intracellular recording of neuronal membrane potential, sharp microelectrodes were made from borosilicate glass capillaries (outer diameter, 1.0 mm; inner diameter, 0.78 mm) with a micropipette puller (Model P-97, Sutter Instrument). The electrodes were filled with 2 m potassium acetate and 0.1 m potassium chloride. The resistance of the microelectrodes was within the range from 20 to 60 MΩ. The membrane potential recordings of the swim interneurons were made using AxoClamp 2B amplifier (Molecular Devices). The output bandwidth of the amplifier was set to 3 kHz. All the output signals from these amplifiers were digitized at a sampling frequency of >2 kHz using a Micro1401 A/D converter (Cambridge Electronic Design) and acquired on a desktop computer running Windows 10 with Spike2 software v8 (Cambridge Electric Design).
Neurons were stimulated by injecting current steps (0.5–3 nA) through a bridge-balanced electrode. To measure electrical connections between two neurons, one of them was impaled with two microelectrodes, one for membrane potential recording and the other for current injection and negative current pulses of varying amplitude (−1 to −4 nA) was injected through the second electrode. The current/voltage relationship was always linear in the neuron being injected with the current pulses of varying amplitude when the voltage change was measured at the peak. To evoke the swim motor pattern, either left or right PdN3 was stimulated with a train of voltage pulses (5 V, 1 ms) at 5 Hz for 1–3 s applied through a suction pipette made from polyethylene tube. A train of stimulus pulses were generated by MASTER-8 pulse stimulator (AMPI) and fed to the electrode through a stimulus isolator (ISO-Flex, AMPI).
Dynamic clamp.
To experimentally boost the strength of the Si3-to-Si2 synapse, we created artificial synaptic currents simultaneously by running Dynamic Clamp software StdpC (Kemenes et al., 2011) on a second Windows computer. Each Si2 was impaled with two microelectrodes, one for membrane potential recording and the other for current injection. Each Si3 was impaled with one microelectrode and the amplifier was set to the discontinuous current-clamp mode. The artificial synaptic current, ISyn, was calculated as described previously (Destexhe et al., 1994; Sharp et al., 1996; Kemenes et al., 2011): where S(t) is the instantaneous synaptic activation, gSyn is the maximum synaptic conductance, and VSyn is the reversal potential of the synapse. In this study VSyn was set to 0 mV for the Si3-to-Si2 synapses based on previous measurements under voltage-clamp (Sakurai and Katz, 2017).
The instantaneous activation, S(t) is given by the differential equation: where S∞ is the steady-state synaptic activation and τSyn is the time constant for synaptic decay. VPre is the presynaptic membrane potential of Si3 and VThresh is the threshold potential for the release of neurotransmitter, which was set to be 50% of the amplitude of the smallest Si3 action potentials. VSlope is the synaptic slope parameter of the activation curve. In this study, τSyn and VSlope was set to 40 ms and 25 mV, respectively, as described previously (Sakurai and Katz, 2016, 2017). gSyn was varied between 40 and 160 nS.
Experimental design and statistical analysis.
A train of Si2 spikes was counted as a burst if the spike interval was <0.5 s and the duration of the spike train was >0.5 s. The burst frequency was measured by the interval between the median spike in each burst of Si1 or Si2 and averaged over 30 s. When bursting in left-right alternation was not detected as described previously (Sakurai and Katz, 2016), then the burst frequency was considered as 0 Hz and the average spike frequency in the duration of >30 s was used as the intraburst spike frequency. Normalized EPSP/C amplitudes in the Figure 7C and D graph were calculated by dividing the amplitude of EPSP or EPSC by the averaged amplitude before the Si1 stimulation.
Statistical comparisons were performed using SigmaPlot v12.5 (Jandel Scientific) for Student's t test (one-tailed), paired t test [see Fig. 7Cii (N = 4), Dii (N = 9)], Pearson correlation [Fig. 2C (N = 34 bursts), D (N = 35 bursts)], and one-way repeated-measures ANOVA [see Figs. 5B,C (N = 11 preparations), 8D1,D2 (N = 8 preparations), 9D (N = 11 preparations), E (N = 7 preparations)] with Holm–Sidak pairwise multiple-comparison procedures. Shapiro–Wilk test was used to test normality of data structure. When the normality test failed, Wilcoxon signed rank test was used instead of a paired t test [Figs. 2B (N = 40), 6C (N = 25)]. Results are expressed as the mean in Figure 7Ci,Di, the mean ± SD in Figures 2B, 5B,C, 8D1,D2, and 9D,E, or plotted in a box plot in Figure 6C. The experimental manipulations of varying strength were given with random order. Missing data were due to the loss of intracellular recording from the neurons.
Results
Si1 firing frequency is correlated with the intensity of the motor pattern
Si1 was previously shown to be extrinsic to swim CPG in Dendronotus (Sakurai et al., 2011). Although it is active throughout the motor program, Si1 does not fire bursts of action potentials like the CPG neurons do (Fig. 2Ai). In 90% of preparations, Si1 was silent before the onset of the swim motor pattern (N = 36 of 40). After triggering the motor pattern with electrical stimulation of a body wall nerve, Si1 started firing action potentials or significantly increased its ongoing spike rate, which peaked quickly before gradually decreasing (Fig. 2Aii). There was a significant difference between the average spike frequency of Si1 before and during the swim motor pattern (N = 40 preparations; Fig. 2B).
In parallel with the increase in Si1 spiking, the burst frequency and intraburst spike frequency of Si2 and Si3 also rapidly increased at the beginning of the swim motor bout and then slowly declined over the course of the bout (Fig. 2Aiii,Aiv). There was a significant correlation between the Si2 burst frequency and Si1 spike frequency (Fig. 2C). Significant positive correlations were seen in 88% of preparations (29 of 33) with the average Pearson correlation coefficient of 0.57 ± 0.18. Similarly, the intraburst spike frequencies of Si2 and Si3 also correlated with Si1 spiking (Fig. 2D). Significant positive correlations were seen in all 33 preparations with the average Pearson correlation coefficient of 0.70 ± 0.15.
Si1 spiking is necessary and sufficient for initiation and maintenance of the swim motor pattern
The swim motor pattern could be terminated by halting Si1 activity. Hyperpolarizing current injection into both Si1s caused a reduction in spiking in Si2 and Si3, leading to a cessation of bursting (Fig. 3A). When the Si1s were released from hyperpolarization, the motor pattern resumed. Similar results were seen in all preparations examined (N = 14). In addition, in 40% of preparations, suppression of spiking in just one Si1 terminated the motor pattern (N = 6 of 15). These results indicate that spiking of at least one Si1 was necessary for the maintenance of the swim motor pattern.
In quiescent preparations, injecting superthreshold depolarizing current into one Si1 (N = 38) or both Si1s (N = 21) induced bursting in the swim CPG neurons (Fig. 3B; N = 53 total preparations). Bursting activity slowly dissipated after termination of the current pulse. Thus, although not a member of the CPG itself, Si1 has a command neuron-like function: it is sufficient to initiate bursting and necessary to maintain it.
Si1 provides excitatory synaptic drive to Si3
Previously established connectivity of Si1 is summarized in Figure 1B. In this study we determine that Si1 has additional synaptic connections (compare Fig. 10A). It was previously shown that each Si1 is electrically coupled to both left and right Si2s (Sakurai et al., 2011). We found that Si1 also has chemical and electrical synapses with both Si3s. In quiescent preparations, when Si1 spiking was induced by injecting depolarizing current, we observed EPSPs in Si3 that corresponded one-for-one with Si1 spikes in all preparations examined (Fig. 3C; N = 43).
To determine whether the synapse from Si1 to Si3 was monosynaptic, we applied Hi-Di saline, which raises the threshold for action potential firing and thus minimizes spontaneous activity and polysynaptic pathways (Fig. 4). Si1 evoked fast and slow EPSPs in Si3 (Fig. 4A1). The fast EPSPs in both left and right Si3s followed the Si1 spikes one-for-one at fixed latencies, indicating that they are monosynaptic (Fig. 4A2). The average amplitude of the Si1-evoked EPSP in the Si3 ipsilateral to Si1 was 2.61 ± 1.57 mV (N = 22) and 2.73 ± 1.94 mV (N = 23) in the contralateral Si3. There was no significant difference in the Si1-evoked EPSPs between left and right Si3s (p = 0.87 by Mann–Whitney rank sum test). The fast EPSPs rode on top of a slow depolarizing potential that lasted for >10 s (Fig. 4A1). The amplitude of the slow component evoked by a Si3 spike train (5 Hz, 5 s) was 5.61 ± 2.72 mV (N = 8). The decay time to 50% was 5.47 ± 2.44 s (N = 8). Both the fast and slow components of Si1-evoked EPSPs in Si3 were greatly diminished in Ca2+-free saline, suggesting that they are mediated by chemical synapses (Fig. 4B; N = 11).
The Ca2+-free saline also revealed an electrotonic component of the excitatory connection between Si1 and both Si3s (Fig. 4B1; N = 11). The electrical synapse was symmetrical; hyperpolarizing current injected into one Si3 produced electrotonic potentials in the contralateral Si3 and both Si1s (Fig. 4C). There was no significant difference in the strength of electrical coupling between the ipsilateral and contralateral neurons; the coupling coefficient between Si1 and the contralateral Si3 was 0.023 ± 0.011 (N = 9) and 0.021 ± 0.008 (N = 7) between Si1 and the ipsilateral Si3 (one-tailed, p = 0.32 by Student's t test). In Ca2+-free saline, Si1 spikes produced small depolarizing potentials (Fig. 4B2), which were mediated by electrical connections. The average amplitude of such depolarizing potentials in the Si3 ipsilateral to Si1 was 0.27 ± 0.20 mV (N = 9) and 0.34 ± 0.18 mV (N = 8) in the contralateral Si3. There was no significant difference between left and right Si3s in the size of Si1-evoked electrotonic potentials (one-tailed, p = 0.23 by Student's t test). Thus, the excitatory action from Si1 to Si3 that is mediated by electrical coupling is approximately one-tenth the size of the chemical excitatory synaptic action. These results suggest that each Si1 equally drives both the left and right Si3 via excitatory synapses and electrical coupling, but the excitatory action is mediated mainly through chemical synapses (Fig. 4D).
Depolarization of Si3 does not account for the Si1-evoked acceleration of the motor pattern
Si1 stimulation accelerates an ongoing swim motor pattern (Fig. 1B2; Sakurai et al., 2011). To determine whether this effect on the swim motor program is caused by the synaptic actions of Si1 on Si3, we mimicked the Si1-induced synaptic depolarization by depolarizing both Si3s with positive current injection. Bilateral depolarization of Si3 did not consistently mimic the effect of Si1 stimulation, increasing burst frequency in 46% of preparations (Fig. 5A1; N = 5 of 11) and decreasing it in 36% of preparations (Fig. 5A2; N = 4). In two preparations, the effect changed from trial to trial. Overall, there was no clear effect of bilateral Si3 depolarization on the swim cycle frequency (Fig. 5B). However, current injection caused a significant increase in the intraburst spike frequency of Si3 (Fig. 5C; p < 0.001 by one-way repeated-measures ANOVA, F(6,36) = 19.5, N = 11 preparations). These results suggest that Si1-evoked depolarization of Si3 increases the firing rate within every burst but does not change the frequency of the bursts.
Si1 potentiates the Si3-to-Si2 synapse
Previously it was shown that the synapse from Si3 to the contralateral Si2 is necessary for the CPG to produce the swim motor pattern (Sakurai and Katz, 2016). Here, we found that in all preparations examined, this synaptic potential increased in amplitude during the swim motor pattern (Fig. 6A–C; N = 25 preparations). Before the occurrence of the swim motor pattern, sparse spiking of Si3 evoked small EPSPs in Si2 (Fig. 6A1), which can be seen when overlaid (Fig. 6A2). During the swim motor pattern, the Si3-evoked EPSPs in Si2 increased in amplitude (Fig. 6B1), corresponding one-to-one with the Si3 spikes (Fig. 6B2). Overall, the amplitude of the initial EPSP during the swim motor pattern was 8.9 ± 4.3 mV (N = 25), which was significantly larger than those evoked before a swim motor pattern 3.9 ± 1.9 mV (Fig. 6C; N = 25, W = 325, p < 0.001 by Wilcoxon signed rank test, N = 25). This increase might have been due at least in part to activity-dependent facilitation of the Si3 synapse. However, because Si1 activity was tightly correlated with the swim motor pattern (compare Fig. 2), it was also feasible that Si1 had a neuromodulatory action on the synapse. Therefore, we examined the direct effect of Si1 onto the Si3 synaptic strength while suppressing spontaneous activity with Hi-Di saline (Fig. 7).
To measure Si3 synaptic strength, Si2 was voltage-clamped at −50 mV in Hi-Di saline while Si3 was stimulated to fire one action potential every 5 s (Fig. 7A). The Si3-evoked EPSCs maintained a constant amplitude until Si1 was stimulated to fire a train of action potentials (5 Hz, 5 s), causing a large increase in the amplitude of the EPSCs (Fig. 7A–C). The magnitude of the increase varied across preparations; the average EPSC amplitude 5 s before Si1 stimulation was 0.73 ± 0.45 nA, which increased to 1.22 ± 0.56 nA 5 s after the stimulation (p = 0.004 by a two-tailed paired t test, N = 4). Similar results were obtained by measuring Si3-evoked EPSPs in Si2 (Fig. 7D). The average EPSP amplitude 5 s before Si1 stimulation was 3.90 ± 1.25 mV, which increased to 6.06 ± 2.15 mV (p = 0.001 by a two-tailed paired t test, N = 9). In both cases, the potentiation lasted ∼1 min (Fig. 7C,D). The Si1 stimulation also caused a small increase in the input resistance of Si2 (4.4 ± 1.7%, N = 5 preparations). Although the increased input resistance of the postsynaptic Si2 may have contributed to the enhancement of EPSPs, the transient large increase in synaptic current suggests that Si1 directly enhances transmitter release from Si3 or modifies the response of receptors in Si2.
Si1 accelerates the swim motor pattern by potentiating the Si3-to-Si2 synapses
To determine the effects of increasing the strength of the Si3-to-Si2 synapses on the swim motor pattern, we performed a “virtual neuromodulation” experiment by artificially increasing the synaptic conductance using the dynamic clamp (Fig. 8). For the dynamic clamp, the amplitude and time course of the synaptic current was calculated in real time (see Materials and Methods) and injected into the postsynaptic Si2 through a second electrode in response to spikes in Si3.
When the Si3-to-Si2 synapse was artificially boosted, Si2 became more likely to fire an action potential in response to Si3 spiking (Fig. 8B; Sakurai and Katz, 2016). Artificial enhancement of both Si3-to-Si2 synapses accelerated the swim motor pattern (Fig. 8C). The burst frequency (Fig. 8D1) and the intraburst spike frequency (Fig. 8D2) of Si2 increased as the artificial synaptic conductance was increased (N = 8 preparations). This demonstrates that increasing the strength of the Si3-to-Si2 synapse alone is sufficient to increase the burst frequency and the intraburst spike frequency.
Combined depolarization of Si3 and potentiation of Si3 synapse is sufficient to evoke the swim motor pattern
In a previous study, we demonstrated that depolarization of both Si2s caused them to fire bursts of long duration that irregularly alternated between the left and right sides. Activation of spiking in the Si3 pair alone did not induce bursting unless the Si2s were also recruited by the excitatory Si3-to-Si2 synapse (Sakurai and Katz, 2016). Here, in a similar experiment, we show that spiking activity in the Si3s alone did not initiate or maintain the swim motor pattern (Fig. 9A1) unless it was strong enough to recruit spiking in the Si2s (Fig. 9A2). Thus, although the Si3 pair is reciprocally inhibitory, it is not sufficient to produce rhythmic bursting; the CPG requires the participation of Si2s firing with Si3s to produce rhythmic bursting activity. Because Si1 has neuromodulatory actions that potentiate Si3-to-Si2 synapses in addition to depolarizing Si3, we examined the effect of synaptic potentiation of generation of the swim motor pattern.
When the Si3s were both depolarized by bilateral current injection, boosting the Si3-to-Si2 synapses by dynamic clamping induced the swim motor pattern by allowing the Si3s to more easily recruit the Si2s (Fig. 9B1). With dynamic clamping, as little as 0.5 nA injected into both Si3s evoked bursting in 100% of preparations (Fig. 9C, black squares). Injection of 2 nA current alone induced rhythmic bursting in >50% of preparations (Fig. 9A2,C, white circles). Once initiated, burst frequency did not increase as a function of current injected into Si3 either with or without dynamic clamping (Fig. 9B2,D). Without dynamic clamping, the Si2 intraburst spike frequency showed no significant change with larger injected currents (N = 11 preparations); however, with dynamic clamping it increased with larger injected current (N = 7 preparations; Fig. 9E). Thus, Si1 acts via a combination of actions to initiate and maintain the swim motor pattern: depolarization of Si3 and enhancement of Si3 synaptic strength. The depolarization of Si3 is crucial for the initiation of spiking activity whereas the synaptic potentiation allows Si3 to reliably recruit robust Si2 spiking.
Discussion
Despite being homologous in Dendronotus and Melibe, Si1 has different functions in the generation of homologous swimming behaviors in these two nudibranch species (Sakurai et al., 2011). In Dendronotus, Si1 is not part of the swim CPG; it does not fire bursts of action potentials in phase with the motor pattern and it lacks reciprocal inhibition across the midline (Sakurai et al., 2011). Here, we showed that Si1 acts as an extrinsic command-like neuron whose non-rhythmic activity is necessary and sufficient for the swim CPG to produce a rhythmic motor pattern.
In Melibe, it was previously shown that Si1 acts as a member of the swim CPG, exhibiting bursting activity during the swim motor pattern, which is maintained through reciprocal inhibition with the contralateral Si1 and Si2 and strong electrical coupling with the ipsilateral Si2 (Thompson and Watson, 2005; Sakurai et al., 2011, 2014). Each Si1 synapses bilaterally onto the Si3 pair, but their actions are opposite between left and right Si3; Si1 makes an excitatory synapse onto the ipsilateral Si3 and a biphasic synapse onto the contralateral Si3 (Fig. 1A1; Sakurai et al., 2014). These asymmetric synapses in Melibe play a role in promoting the transition of Si3 bursting from one side to the other and thus are important for setting the phase relationship between two half-center kernels (Sakurai et al., 2014). Thus, homologous neurons in these two species differ in their functional positions in the motor system hierarchy; Si1 is part of the swim CPG in Melibe, but serves as a higher-level command neuron in Dendronotus.
Dual actions by Si1 activate the motor pattern
In Dendronotus, each Si1 has two complementary bilateral actions (Fig. 10A). The first action of Si1 is to evoke fast and slow excitatory potentials in both Si3s, which summate to cause them to fire action potentials. The second action of Si1 is to heterosynaptically potentiate the Si3-to-Si2 synapses. This potentiation is necessary to allow Si3 to excite Si2 above threshold for firing action potentials. Mimicking these two actions by depolarization of Si3s and dynamic clamping of the Si3-to-Si2 synapses elicits a swim-like motor pattern. Thus, the dual effects of Si1, depolarizing Si3 and enhancing its synaptic strength, are sufficient for the CPG to produce rhythmic bursting activity.
Because Si1 has electrical coupling with all other neurons involved in the circuit, it is likely that the electrical coupling also contributes to generation of the swim motor pattern to some extent. However, based on the facts that (1) the chemical synaptic potentials are ∼10 times larger than the electrotonic potentials and that (2) the swim motor pattern can be blocked by reducing the strength of the Si3-to-Si2 chemical synapse either by pharmacological manipulation (Sakurai and Katz, 2017) or by dynamic clamping (Sakurai and Katz, 2016), we conclude that Si1 initiates and maintains the swim motor pattern mainly through its chemical synaptic action onto Si3 and the heterosynaptic potentiation of Si3 synapses.
Activation of Si1 also accelerates the on-going swim motor pattern. In this study, we demonstrated that artificially boosting the Si3-to-Si2 synapses by dynamic clamping mimicked the acceleratory effect of Si1, but simple bilateral depolarization of Si3s did not. In the Dendronotus swim CPG, it was previously shown that neither Si2 nor Si3 possesses any endogenous bursting properties, but action potential firing in both neurons is essential to generate rhythmic bursting (Sakurai and Katz, 2016). It is still unclear how strengthening the excitatory synapse from Si3 to Si2 accelerates the oscillatory activity. We hypothesize that Si3 not only causes Si2 to fire, but also that it may boost dynamic membrane properties of Si2 and/or Si3. In a half-center configuration, the dynamic properties of neuronal components play crucial roles in prompting phase transitions from one side to the other (Friesen and Stent, 1978; Arbas and Calabrese, 1987; Friesen, 1994).
Si1 as a command neuron
Since the early description of “command neurons” by Wiersma and Ikeda (1964), there have been many reports describing the functions of command neurons (Kupfermann and Weiss, 1978; Puhl et al., 2012; Bouvier et al., 2015; Hampel et al., 2015; Kim et al., 2017). Extrinsic command neurons are a common feature of motor systems such as leech swimming (Kristan et al., 2005), insect leg movements (Berg et al., 2015), and vertebrate locomotion (Dubuc et al., 2008; McLean et al., 2008; Hägglund et al., 2010; Kimura et al., 2013; Bouvier et al., 2015; Grillner and El Manira, 2015).
Si1 fulfills the criteria for a command neuron (cf. Kupfermann and Weiss, 1978): (1) Si1 always fired during the swim motor pattern. (2) Si1 is necessary for production of the swim motor pattern; suppression of both Si1s blocked the swim motor pattern. (3) Activation of both Si1s is sufficient to initiate the swim motor pattern.
The synaptic and neuromodulatory mechanism underlying the Si1 command
The Dendronotus organization closely resembles that originally proposed for a half-center oscillator with excitatory drive arising from a command neuron (Fig. 10B,C; Friesen, 1994). In the vertebrate brainstem and spinal cord, respiratory and locomotor CPGs convert the tonic brainstem input into phasic motor patterns (Dubuc et al., 2008; Hägglund et al., 2010; Kimura et al., 2013; Grillner and El Manira, 2015; but see Soffe et al., 2009; Li et al., 2010). However, little is known about how extrinsic tonic input turns on an emergent oscillatory circuit with no intrinsically bursting neurons.
Here, we determined that the excitatory synaptic inputs onto the CPG neurons and the potentiation of their synapses mediate the command-like action of Si1. Moreover, mimicking the dual bilateral actions of Si1 by dynamic clamping induced rhythmic bursting of Si2 and Si3 that resembled the swim motor pattern. These experiments directly showed that the synaptic and neuromodulatory actions of Si1 are sufficient to induce the swim motor pattern.
Dendronotus command mechanism differs from that of Tritonia
In another nudibranch, Tritonia diomedea, the mechanism for commanding a swim motor pattern separates the excitatory synaptic drive from the neuromodulation. The Dorsal Ramp Interneuron (DRI) is the command neuron for dorsal-ventral flexion swimming (Frost and Katz, 1996). It evokes large EPSPs in a set of swim CPG neurons (Dorsal Swim Interneurons, DSIs) similar to the large EPSPs evoked by Si1 in Si3. However, unlike Si1, DRI does not enhance the strength of synapses made by the DSIs. Instead, the DSIs themselves enhance the strength of synapses made by other CPG neurons by releasing serotonin (Katz et al., 1994; Katz and Frost, 1995a,b; Sakurai and Katz, 2003). Thus, the functionality of excitation and neuromodulation is divided between two neurons in Tritonia, one of which is intrinsic to the swim CPG, whereas in Dendronotus, a single neuron, Si1, combines both actions.
The Dendronotus command is separate from the swim CPG
The line between command neuron and CPG neuron can be blurred in some cases. Extrinsic modulatory neurons can receive rhythmic drive from the CPG causing them to oscillate when the CPG is activated (Coleman and Nusbaum, 1994; Coleman et al., 1995). In Tritonia, DRI also receives feedback from the swim CPG and has its spiking shaped by the motor pattern (Frost and Katz, 1996). In this study we found no evidence that Si1 receives feedback from the CPG neurons except for electrical coupling, which is ∼10× weaker than the ipsilateral Si1-to-Si2 coupling in Melibe (Sakurai et al., 2011). Thus, in Dendronotus, there is a distinct separation of command and CPG function.
Anatomical organization of the hierarchical structure
It may be significant that the soma of Si1 resides in the cerebral ganglion. In gastropod molluscs, command neurons tend to be in the cerebral ganglia which receive inputs from visual, tactile, chemical, and vestibular senses. For example, several identified neurons in the cerebral ganglia of the pteropod Clione project axons to the pedal ganglia to control the swim motor pattern (Panchin et al., 1995; Satterlie and Norekian, 1995; Deliagina et al., 1998). Similarly, a group of neurons in the cerebral ganglion of Aplysia were shown to be necessary and sufficient to elicit for locomotor activity from the pedal ganglia (Fredman and Jahan-Parwar, 1983; Jing et al., 2008). Command neurons for the feeding motor network in the buccal ganglia are also found in the cerebral ganglia in Aplysia (Jing et al., 2004), Pleurobranchaea (Gillette et al., 1978, 1982), Lymnaea (Norekian and Satterlie, 1993), and Clione (Alania et al., 2008). The somata of the CPG neurons Si2 and Si3 are located in the pedal ganglia, which is where efferent neurons in gastropods generally reside.
The presence of Si1 in the cerebral ganglion suggests that Dendronotus might represent the ancestral state with Si1 having a command role. Si1 in Melibe might have been co-opted into a rhythmic CPG role (Sakurai et al., 2014). It would also be interesting to know how the motor hierarchy is organized in other species, such as Flabellina iodinea, which evolved the swimming behavior independently of Dendronotus and Melibe and has a rhythmically active Si1 (Gunaratne et al., 2017). More species need to be sampled to better ascertain the direction of evolutionary change.
Footnotes
This work was supported by NSF Grant IOS-1455527 to P.S.K. We thank Dr. Andrey Shilnikov for helpful discussions, and Trevor Fay and Andrew Kim at the Monterey Abalone Company for their efforts in collecting and shipping animals to meet our needs.
The authors declare no competing financial interests.
- Correspondence should be addressed to Akira Sakurai at akira{at}gsu.edu