Single-cell RNAseq analysis of spinal locomotor circuitry in larval zebrafish.

Identification of the neuronal types that form the specialized circuits controlling distinct behaviors has benefited greatly from the simplicity offered by zebrafish. Electrophysiological studies have shown that in addition to connectivity, understanding of circuitry requires identification of functional specializations among individual circuit components, such as those that regulate levels of transmitter release and neuronal excitability. In this study, we use single-cell RNA sequencing (scRNAseq) to identify the molecular bases for functional distinctions between motoneuron types that are causal to their differential roles in swimming. The primary motoneuron, in particular, expresses high levels of a unique combination of voltage-dependent ion channel types and synaptic proteins termed functional 'cassettes.' The ion channel types are specialized for promoting high-frequency firing of action potentials and augmented transmitter release at the neuromuscular junction, both contributing to greater power generation. Our transcriptional profiling of spinal neurons further assigns expression of this cassette to specific interneuron types also involved in the central circuitry controlling high-speed swimming and escape behaviors. Our analysis highlights the utility of scRNAseq in functional characterization of neuronal circuitry, in addition to providing a gene expression resource for studying cell type diversity.

Single cell RNA-seq analysis of spinal locomotor circuitry in larval zebrafish.

Identification of the neuronal types that form the specialized circuits controlling distinct behaviors has benefited greatly from the simplicity offered by zebrafish. Electrophysiological studies have shown that additional to connectivity, understanding of circuitry requires identification of functional specializations among individual circuit components, such as those that regulate levels of transmitter release and neuronal excitability. In this study we use single cell RNA sequencing (scRNAseq) to identify the molecular bases for functional distinctions between motoneuron types that are causal to their differential roles in swimming. The primary motoneuron (PMn) in particular, expresses high levels of a unique combination of voltage-dependent ion channel types and synaptic proteins termed functional 'cassettes'. The ion channel types are specialized for promoting high frequency firing of action potentials and augmented transmitter release at the neuromuscular junction, both contributing to greater power generation. Our transcriptional profiling of spinal neurons further assigns expression of this cassette to specific interneuron types also involved in the central circuitry controlling high speed swimming and escape behaviors. Our analysis highlights the utility of scRNAseq in functional characterization of neuronal circuitry, in addition to providing a gene expression resource for studying cell type diversity.

Primary and secondary motoneurons use different calcium channel types to control escape and swimming behaviors in zebrafish.

The escape response and rhythmic swimming in zebrafish are distinct behaviors mediated by two functionally distinct motoneuron (Mn) types. The primary (1°Mn) type depresses and has a large quantal content (Qc) and a high release probability (Pr). Conversely, the secondary (2°Mn) type facilitates and has low and variable Qc and Pr. This functional duality matches well the distinct associated behaviors, with the 1°Mn providing the strong, singular C bend initiating escape and the 2°Mn conferring weaker, rhythmic contractions. Contributing to these functional distinctions is our identification of P/Q-type calcium channels mediating transmitter release in 1°Mns and N-type channels in 2°Mns. Remarkably, despite these functional and behavioral distinctions, all ∼15 individual synapses on each muscle cell are shared by a 1°Mn bouton and at least one 2°Mn bouton. This blueprint of synaptic sharing provides an efficient way of controlling two different behaviors at the level of a single postsynaptic cell.

Zebrafish neuromuscular junction: The power of N.

In the early 1950s, Katz and his colleagues capitalized on the newly developed intracellular microelectrode recording technique to investigate synaptic transmission. For study they chose frog neuromuscular junction (NMJ), which was ideally suited due to the accessibility and large size of the muscle cells. Paradoxically, the large size precluded the use of next generation patch clamp technology. Consequently, electrophysiological study of synaptic function shifted to small central synapses made amenable by patch clamp. Recently, however, the unique features offered by zebrafish have rekindled interest in the NMJ as a model for electrophysiological study of synaptic transmission. The small muscle size and synaptic simplicity provide the singular opportunity to perform in vivo spinal motoneuron-target muscle patch clamp recordings. Additional incentive is provided by zebrafish lines harboring mutations in key synaptic proteins, many of which are embryonic lethal in mammals, but all of which are able to survive well past synapse maturation in zebrafish. This mini-review will highlight features that set zebrafish NMJs apart from traditional NMJs. We also draw into focus findings that offer the promise of identifying features that define release sites, which serve to set the upper limit of transmitter release. Since its conception several candidates representing release sites have been proposed, most of which are based on distinctions among vesicle pools in their state of readiness for release. However, models based on distinctions among vesicles have become enormously complicated and none adequately account for setting an upper limit for exocytosis in response to an action potential (AP). Specifically, findings from zebrafish NMJ point to an alternative model, positing that elements other than vesicles per se set the upper limits of release.

Astrocytic modulation of excitatory synaptic signaling in a mouse model of Rett syndrome.

Studies linking mutations in Methyl CpG Binding Protein 2 (MeCP2) to physiological defects in the neurological disease, Rett syndrome, have focused largely upon neuronal dysfunction despite MeCP2 ubiquitous expression. Here we explore roles for astrocytes in neuronal network function using cortical slice recordings. We find that astrocyte stimulation in wild-type mice increases excitatory synaptic activity that is absent in male mice lacking MeCP2 globally. To determine the cellular basis of the defect, we exploit a female mouse model for Rett syndrome that expresses wild-type MeCP2-GFP in a mosaic distribution throughout the brain, allowing us to test all combinations of wild-type and mutant cells. We find that the defect is dependent upon MeCP2 expression status in the astrocytes and not in the neurons. Our findings highlight a new role for astrocytes in regulation of excitatory synaptic signaling and in the neurological defects associated with Rett syndrome.

A Gradient in Synaptic Strength and Plasticity among Motoneurons Provides a Peripheral Mechanism for Locomotor Control.

The recruitment of motoneurons during force generation follows a general pattern that has been confirmed across diverse species [1-3]. Motoneurons are recruited systematically according to synaptic inputs and intrinsic cellular properties and corresponding to movements of different intensities. However, much less is known about the output properties of individual motoneurons and how they affect the translation of motoneuron recruitment to the strength of muscle contractions. In larval zebrafish, spinal motoneurons are recruited in a topographic gradient according to their input resistance (Rin) at different swimming strengths and speeds. Whereas dorsal, lower-Rin primary motoneurons (PMns) are only activated during behaviors that involve strong and fast body bends, more ventral, higher-Rin secondary motoneurons (SMns) are recruited during weaker and slower movements [4-6]. Here we perform in vivo paired recordings between identified spinal motoneurons and skeletal muscle cells in larval zebrafish. We characterize individual motoneuron outputs to single muscle cells and show that the strength and reliability of motoneuron outputs are inversely correlated with motoneuron Rin. During repetitive high-frequency motoneuron drive, PMn synapses undergo depression, whereas SMn synapses potentiate. We monitor muscle cell contractions elicited by single motoneurons and show that the pattern of motoneuron output strength and plasticity observed in electrophysiological recordings is reflected in muscle shortening. Our findings indicate a link between the recruitment pattern and output properties of spinal motoneurons that can together generate appropriate intensities for muscle contractions. We demonstrate that motoneuron output properties provide an additional peripheral mechanism for graded locomotor control at the neuromuscular junction.

Fatigue in Rapsyn-Deficient Zebrafish Reflects Defective Transmitter Release.

Rapsyn-deficient myasthenic syndrome is characterized by a weakness in voluntary muscle contraction, a direct consequence of greatly reduced synaptic responses that result from poorly clustered acetylcholine receptors. As with other myasthenic syndromes, the general muscle weakness is also accompanied by use-dependent fatigue. Here, we used paired motor neuron target muscle patch-clamp recordings from a rapsyn-deficient mutant line of zebrafish to explore for the first time the mechanisms causal to fatigue. We find that synaptic responses in mutant fish can follow faithfully low-frequency stimuli despite the reduced amplitude. This is in part helped by a compensatory increase in the number of presynaptic release sites in the mutant fish. In response to high-frequency stimulation, both wild-type and mutant neuromuscular junctions depress to steady-state response levels, but the latter shows exaggerated depression. Analysis of the steady-state transmission revealed that vesicle reloading and release at individual release sites is significantly slower in mutant fish during high-frequency activities. Therefore, reductions in postsynaptic receptor density and compromised presynaptic release collectively serve to reduce synaptic strength to levels that fall below the threshold for muscle action potential generation, thus accounting for use-dependent fatigue. Our findings raise the possibility that defects in motor neuron function may also be at play in other myasthenic syndromes that have been mapped to mutations in muscle-specific proteins. SIGNIFICANCE STATEMENT: Use-dependent fatigue accompanies many neuromuscular myasthenic syndromes, including muscle rapsyn deficiency. Here, using a rapsyn-deficient line of zebrafish, we performed paired motor neuron target muscle patch-clamp recordings to investigate the mechanisms causal to this phenomenon. Our findings indicate that the reduced postsynaptic receptor density resulting from defective rapsyn contributes to weakness, but is not solely responsible for use-dependent fatigue. Instead, we find unexpected involvement of altered transmitter release from the motor neuron. Specifically, slowed reloading of vesicle release sites leads to augmented synaptic depression during repeated action potentials. Even at moderate stimulus frequencies, the depression levels for evoked synaptic responses fall below the threshold for the generation of muscle action potentials. The associated contraction failures are manifest as use-dependent fatigue.

Nonequivalent release sites govern synaptic depression.

Synaptic depression is prominent among synapses, but the underlying mechanisms remain uncertain. Here, we use paired patch clamp recording to study neuromuscular transmission between the caudal primary motor neuron and target skeletal muscle in zebrafish. This synapse has an unusually low number of release sites, all with high probabilities of release in response to low-frequency stimulation. During high-frequency stimulation, the synapse undergoes short-term depression and reaches steady-state levels of transmission that sustain the swimming behavior. To determine the release parameters underlying this steady state, we applied variance analysis. Our analysis revealed two functionally distinct subclasses of release sites differing by over 60-fold in rates of vesicle reloading. A slow reloading class requires seconds to recover and contributes to depression onset but not the steady-state transmission. By contrast, a fast reloading class recovers within tens of milliseconds and is solely responsible for steady-state transmission. Thus, in contrast to most current models that assign levels of steady-state depression to vesicle availability, our findings instead assign this function to nonuniform release site kinetics. The duality of active-site properties accounts for the highly nonlinear dependence of steady-state depression levels on frequency.

Zebrafish CaV2.1 calcium channels are tailored for fast synchronous neuromuscular transmission.

The CaV2.2 (N-type) and CaV2.1 (P/Q-type) voltage-dependent calcium channels are prevalent throughout the nervous system where they mediate synaptic transmission, but the basis for the selective presence at individual synapses still remains an open question. The CaV2.1 channels have been proposed to respond more effectively to brief action potentials (APs), an idea supported by computational modeling. However, the side-by-side comparison of CaV2.1 and CaV2.2 kinetics in intact neurons failed to reveal differences. As an alternative means for direct functional comparison we expressed zebrafish CaV2.1 and CaV2.2 α-subunits, along with their accessory subunits, in HEK293 cells. HEK cells lack calcium currents, thereby circumventing the need for pharmacological inhibition of mixed calcium channel isoforms present in neurons. HEK cells also have a simplified morphology compared to neurons, which improves voltage control. Our measurements revealed faster kinetics and shallower voltage-dependence of activation and deactivation for CaV2.1. Additionally, recordings of calcium current in response to a command waveform based on the motorneuron AP show, directly, more effective activation of CaV2.1. Analysis of calcium currents associated with the AP waveform indicate an approximately fourfold greater open probability (PO) for CaV2.1. The efficient activation of CaV2.1 channels during APs may contribute to the highly reliable transmission at zebrafish neuromuscular junctions.

Synchronous and asynchronous modes of synaptic transmission utilize different calcium sources.

Asynchronous transmission plays a prominent role at certain synapses but lacks the mechanistic insights of its synchronous counterpart. The current view posits that triggering of asynchronous release during repetitive stimulation involves expansion of the same calcium domains underlying synchronous transmission. In this study, live imaging and paired patch clamp recording at the zebrafish neuromuscular synapse reveal contributions by spatially distinct calcium sources. Synchronous release is tied to calcium entry into synaptic boutons via P/Q type calcium channels, whereas asynchronous release is boosted by a propagating intracellular calcium source initiated at off-synaptic locations in the axon and axonal branch points. This secondary calcium source fully accounts for the persistence following termination of the stimulus and sensitivity to slow calcium buffers reported for asynchronous release. The neuromuscular junction and CNS neurons share these features, raising the possibility that secondary calcium sources are common among synapses with prominent asynchronous release. DOI: http://dx.doi.org/10.7554/eLife.01206.001.

Zebrafish calls for reinterpretation for the roles of P/Q calcium channels in neuromuscular transmission.

A long-held tenet of neuromuscular transmission is that calcium-dependent neurotransmitter release is mediated by N-type calcium channels in frog but P/Q-type channels in mammals. The N-type assignment in frog is based principally on pharmacological sensitivity to ω-conotoxin GVIA. Our studies show that zebrafish neuromuscular transmission is also sensitive to ω-conotoxin GVIA. However, positional cloning of a mutant line with compromised neuromuscular function identified a mutation in a P/Q- rather than N-type channel. Cloning and heterologous expression of this P/Q-type channel confirmed a block by ω-conotoxin GVIA raising the likelihood that all vertebrates, including frog, use the P/Q-type calcium channel for neuromuscular transmission. In addition, our P/Q defective mutant line offered a means of testing the ability of roscovitine, known to potentiate frog neuromuscular transmission, to mediate behavioral and functional rescue. Acute treatment led to rapid improvement of both, pointing to potential therapeutic benefit for myasthenic disorders involving calcium channel dysfunction.

Zebrafish model for congenital myasthenic syndrome reveals mechanisms causal to developmental recovery.

Mutations in muscle ACh receptors cause slow-channel syndrome (SCS) and Escobar syndrome, two forms of congenital myasthenia. SCS is a dominant disorder with mutations reported for all receptor subunits except γ. Escobar syndrome is distinct, with mutations located exclusively in γ, and characterized by developmental improvement of muscle function. The zebrafish mutant line, twister, models SCS in terms of a dominant mutation in the α subunit (α(twi)) but shows the behavioral improvement associated with Escobar syndrome. Here, we present a unique electrophysiological study into developmental improvement for a myasthenic syndrome. The embryonic α(twi)ßδγ receptor isoform produces slowly decaying synaptic currents typical of SCS that transit to a much faster decay upon the appearance of adult ε, despite the α(twi) mutation. Thus, the continued expression of α(twi) into adulthood is tolerated because of the ε expression and associated recovery, raising the likelihood of unappreciated myasthenic cases that benefit from the γ-ε switch.

Acetylcholine receptor gating in a zebrafish model for slow-channel syndrome.

Slow-channel syndrome (SCS) is an autosomal-dominant disease resulting from mutations in muscle acetylcholine (ACh) receptor subunits. The associated fatigue and muscle degeneration are proposed to result from prolonged synaptic responses that overload intracellular calcium. Single-channel studies on reconstituted receptors bearing human mutations indicate that the prolonged responses result from an increase in receptor open duration and, in some cases, increased sensitivity to ACh. We show that both of these aberrant receptor properties are recapitulated in heterozygotic zebrafish bearing an L258P mutation in the α subunit, thus affording the unique opportunity to compare the single-channel properties of mutant receptors to the synaptic currents in vivo. Whole-cell recordings revealed synaptic currents that decayed along a multiexponential time course, reflecting receptors containing mixtures of wild-type and mutant α subunits. Treatment with quinidine, an open-channel blocker used to treat the human disorder, restored fast synaptic current kinetics and the ability to swim. Quinidine block also revealed that mutant receptors generate a large steady-state current in the absence of ACh. The spontaneous openings reflected a destabilization of the closed state, leading to an apparent increase in the sensitivity of these receptors to ACh. The effective block by quinidine on synaptic currents as well as nonliganded openings points to dual sources for the calcium-dependent myopathy in certain forms of SCS.

An acetylcholine receptor lacking both γ and ε subunits mediates transmission in zebrafish slow muscle synapses.

Fast and slow skeletal muscle types in larval zebrafish can be distinguished by a fivefold difference in the time course of their synaptic decay. Single-channel recordings indicate that this difference is conferred through kinetically distinct nicotinic acetylcholine receptor (AChR) isoforms. The underlying basis for this distinction was explored by cloning zebrafish muscle AChR subunit cDNAs and expressing them in Xenopus laevis oocytes. Measurements of single-channel conductance and mean open burst duration assigned α(2)ßδε to fast muscle synaptic current. Contrary to expectations, receptors composed of only αßδ subunits (presumed to be α(2)ßδ(2) receptors) recapitulated the kinetics and conductance of slow muscle single-channel currents. Additional evidence in support of γ/ε-less receptors as mediators of slow muscle synapses was reflected in the inward current rectification of heterologously expressed α(2)ßδ(2) receptors, a property normally associated with neuronal-type nicotinic receptors. Similar rectification was reflected in both single-channel and synaptic currents in slow muscle, distinguishing them from fast muscle. The final evidence for α(2)ßδ(2) receptors in slow muscle was provided by our ability to convert fast muscle synaptic currents to those of slow muscle by knocking down ε subunit expression in vivo. Thus, for the first time, muscle synaptic function can be ascribed to a receptor isoform that is composed of only three different subunits. The unique functional features offered by the α(2)ßδ(2) receptor likely play a central role in mediating the persistent contractions characteristic to this muscle type.

Paired patch clamp recordings from motor-neuron and target skeletal muscle in zebrafish.

Larval zebrafish represent the first vertebrate model system to allow simultaneous patch clamp recording from a spinal motor-neuron and target muscle. This is a direct consequence of the accessibility to both cell types and ability to visually distinguish the single segmental CaP motor-neuron on the basis of morphology and location. This video demonstrates the microscopic methods used to identify a CaP motor-neuron and target muscle cells as well as the methodologies for recording from each cell type. Identification of the CaP motor-neuron type is confirmed by either dye filling or by the biophysical features such as action potential waveform and cell input resistance. Motor-neuron recordings routinely last for one hour permitting long-term recordings from multiple different target muscle cells. Control over the motor-neuron firing pattern enables measurements of the frequency-dependence of synaptic transmission at the neuromuscular junction. Owing to a large quantal size and the low noise provided by whole cell voltage clamp, all of the unitary events can be resolved in muscle. This feature permits study of basic synaptic properties such as release properties, vesicle recycling, as well as synaptic depression and facilitation. The advantages offered by this in vivo preparation eclipse previous neuromuscular model systems studied wherein the motor-neurons are usually stimulated by extracellular electrodes and the muscles are too large for whole cell patch clamp. The zebrafish preparation is amenable to combining electrophysiological analysis with a wide range of approaches including transgenic lines, morpholino knockdown, pharmacological intervention and in vivo imaging. These approaches, coupled with the growing number of neuromuscular disease models provided by mutant lines of zebrafish, open the door for new understanding of human neuromuscular disorders.

Distinct roles for two synaptotagmin isoforms in synchronous and asynchronous transmitter release at zebrafish neuromuscular junction.

An obligatory role for the calcium sensor synaptotagmins in stimulus-coupled release of neurotransmitter is well established, but a role for synaptotagmin isoform involvement in asynchronous release remains conjecture. We show, at the zebrafish neuromuscular synapse, that two separate synaptotagmins underlie these processes. Specifically, knockdown of synaptotagmin 2 (syt2) reduces synchronous release, whereas knockdown of synaptotagmin 7 (syt7) reduces the asynchronous component of release. The zebrafish neuromuscular junction is unique in having a very small quantal content and a high release probability under conditions of either low-frequency stimulation or high-frequency augmentation. Through these features, we further determined that during the height of shared synchronous and asynchronous transmission these two modes compete for the same release sites.

Synaptic homeostasis in a zebrafish glial glycine transporter mutant.

Truncated escape responses characteristic of the zebrafish shocked mutant result from a defective glial glycine transporter (GlyT1). In homozygous GlyT1 mutants, irrigating brain ventricles with glycine-free solution rescues normal swimming. Conversely, elevating brain glycine levels restores motility defects. These experiments are consistent with previous studies that demonstrate regulation of global glycine levels in the CNS as a primary function of GlyT1. As GlyT1 mutants mature, their ability to mount an escape response naturally recovers. To understand the basis of this recovery, we assay synaptic transmission in primary spinal motor neurons by measuring stimulus-evoked postsynaptic potentials. At the peak of the motility defect, inhibitory synaptic potentials are both significantly larger and more prolonged indicating a prominent role for GlyT1 in shaping fast synaptic transmission. However, as GlyT1 mutants naturally regain their ability to swim, the amplitude of inhibitory potentials decreases to below wild-type levels. In parallel with diminishing synaptic potentials, the glycine concentration required to evoke the mutant motility defect increases 61-fold during behavioral recovery. Behavioral recovery is also mirrored by a reduction in the levels of both glycine receptor protein and transcript. These results suggest that increased CNS glycine tolerance and reduced glycine receptor expression in GlyT1 mutants reflect compensatory mechanisms for functional recovery from excess nervous system inhibition.

Function of neuromuscular synapses in the zebrafish choline-acetyltransferase mutant bajan.

We have identified a zebrafish mutant line, bajan, in which compromised motility and fatigue result from a point mutation in the gene coding choline acetyltransferase (ChAT), the enzyme responsible for acetylcholine (ACh) synthesis. Although the mutation predicts loss of ChAT function, bajan inexplicably retains low levels of neuromuscular transmission. We exploited this residual activity and determined the consequences for synaptic function. The attenuated synaptic responses were a direct consequence of a decrease in both resting mean quantal size and quantal content. To replicate behavioral fatigue in swimming, motorneurons were stimulated at high frequencies. A prominent reduction in quantal content, reflecting vesicle depletion, was coincident with a small additional reduction in quantal size. In humans, defective ChAT leads to episodic apnea, a form of congenital myasthenic syndrome characterized by use-dependent fatigue. In contrast to bajan, however, afflicted individuals exhibit a normal resting quantal size and quantal content. The fatigue in humans results from a pronounced long-lasting drop in quantal size with little or no change in quantal content. These differences have important implications for interpreting fatigue as well as on understanding the impact of ACh availability on vesicle filling and recycling.

An electrically coupled network of skeletal muscle in zebrafish distributes synaptic current.

Fast and slow skeletal muscle types are readily distinguished in larval zebrafish on the basis of differences in location and orientation. Additionally, both muscle types are compact, rendering them amenable to in vivo patch clamp study of synaptic function. Slow muscle mediates rhythmic swimming, but it does so purely through synaptic drive, as these cells are unable to generate action potentials. Our patch clamp recordings from muscle pairs of zebrafish reveal a network of electrical coupling in slow muscle that allows sharing of synaptic current within and between segmental boundaries of the tail. The synaptic current exhibits slow kinetics (tau(decay) approximately 4 ms), which further facilitates passage through the low pass filter, a consequence of the electrically coupled network. In contrast to slow muscle, fast skeletal muscle generates action potentials to mediate the initial rapid component of the escape response. The combination of very weak electrical coupling and synaptic kinetics (tau(decay) <1 ms) too fast for the network low pass filter minimizes intercellular sharing of synaptic current in fast muscle. These differences between muscle types provide insights into the physiological role(s) of electrical coupling in skeletal muscle. First, intrasegmental coupling among slow muscle cells allows effective transfer of synaptic currents within tail segments, thereby minimizing differences in synaptic depolarization. Second, a fixed intersegmental delay in synaptic current transit, resulting from the low pass filter properties of the slow muscle network, helps coordinate the rostral-caudal wave of contraction.

Paired motor neuron-muscle recordings in zebrafish test the receptor blockade model for shaping synaptic current.

The transparent spinal cord and electrically compact fast muscle of zebrafish offer the first opportunity to perform simultaneous whole-cell patch-clamp recordings from both motor neuron and target skeletal muscle in situ. Our paired recordings reveal the fastest reported kinetics for both spontaneous and evoked synaptic currents at any synapse and a large quantal size that facilitates the resolution of spontaneous synaptic currents. We used this preparation to test the recent proposal that open channel block of the acetylcholine receptor by acetylcholine modulates the kinetics and timing of transmission between nerve and muscle in larval zebrafish (Legendre et al., 2000). Contrary to the predictions of this model, we find similar delay and onset kinetics of synaptic current at positive and negative muscle membrane potentials, even after inhibition of acetylcholinesterase. In contrast, blockade of motor neuron K channels by 4-aminopyridine prolonged the action potential and introduced a significant delay and slowing of evoked synaptic currents, demonstrating our ability to measured altered transmitter release with this system. We conclude that the kinetics of neuromuscular synaptic currents in zebrafish is not governed by receptor block.