Abstract
Behaviorally adequate neuronal firing patterns are critically dependent on the specific types of ion channel expressed and on their subcellular localization. This study combines in situ electrophysiology with genetic and pharmacological intervention in larval Drosophila melanogaster of both sexes to address localization and function of L-type like calcium channels in motoneurons. We demonstrate that Dmca1D (Cav1 homolog) L-type like calcium channels localize to both the somatodendritic and the axonal compartment of larval crawling motoneurons. In situ patch-clamp recordings in genetic mosaics reveal that Dmca1D channels increase burst duration and maximum intraburst firing frequencies during crawling-like motor patterns in semi-intact animals. Genetic and acute pharmacological manipulations suggest that prolonged burst durations are caused by dendritically localized Dmca1D channels, which activate upon cholinergic synaptic input and amplify EPSPs, thus indicating a conserved function of dendritic L-type channels from Drosophila to vertebrates. By contrast, maximum intraburst firing rates require axonal calcium influx through Dmca1D channels, likely to enhance sodium channel de-inactivation via a fast afterhyperpolarization through BK channel activation. Therefore, in unmyelinated Drosophila motoneurons different functions of axonal and dendritic L-type like calcium channels likely operate synergistically to maximize firing output during locomotion.
SIGNIFICANCE STATEMENT Nervous system function depends on the specific excitabilities of different types of neurons. Excitability is largely shaped by different combinations of voltage-dependent ion channels. Despite a high degree of conservation, the huge diversity of ion channel types and their differential localization pose challenges in assigning distinct functions to specific channels across species. We find a conserved role, from fruit flies to mammals, for L-type calcium channels in augmenting motoneuron excitability. As in spinal cord, dendritic L-type channels amplify excitatory synaptic input. In contrast to spinal motoneurons, axonal L-type channels enhance firing rates in unmyelinated Drosophila motoraxons. Therefore, enhancing motoneuron excitability by L-type channels seems an old strategy, but localization and interactions with other channels are tuned to species-specific requirements.
Introduction
Motoneurons (MNs) are the relay output from the central pattern generating (CPG) networks to the muscles. Although CPGs generate the principle patterns of spiking output to the muscles (Kiehn and Kullander, 2004; Grillner et al., 2005; Marder et al., 2005), the final motor output depends not only on network connectivity but also on the intrinsic properties of MNs, because these compute synaptic input far beyond threshold summation (Kiehn, 1991; Kiehn et al., 2000; Heckman et al., 2003). Therefore, addressing the specific contributions of MN membrane currents to the production of rhythmic motor output is critical to understanding the control of movement in healthy and diseased animals.
L-type Ca2+ channels play important roles in regulating spinal MN excitability. Localized to dendrites, these channels may contribute to persistent inward current (PIC) and amplify synaptic input, especially in the presence of aminergic modulation (Heckman et al., 2003, 2005; Hultborn et al., 2013). The Drosophila homolog of vertebrate Cav1 L-type Ca2+ channels is Dmca1D (Zheng et al., 1995). Somatic patch-clamp recordings from larval Drosophila crawling MNs indicate a conserved somatodendritic localization of L-type channels (Worrell and Levine, 2008). However, the function of L-type channels in invertebrate neurons and the degree to which the regulation of MN excitability is conserved from insects to vertebrates remain incompletely understood.
This study shows directly by immunocytochemistry that Dmca1D localizes to patches of larval Drosophila MN dendritic membrane, but in contrast to spinal MNs, L-type channels also localize to the axons of Drosophila motoneurons. In situ patch-clamp recordings in genetic mosaics reveal that L-type like channels increase burst duration and maximum intraburst firing frequencies during locomotion. Combining imaging and current-clamp recordings with acute pharmacological manipulation indicates that dendritic channels function to amplify excitatory synaptic drive, as in vertebrates. By contrast, axonal Ca2+ influx through L-type channels promotes maximum MN firing rates by enhancing Na+ channel de-inactivation via a fast afterhyperpolarization that is caused by transient BK channel activation (Kadas et al., 2015).
Materials and Methods
Animals.
Drosophila melanogaster were reared in standard plastic vials with foam stoppers on a yeast-cornmeal-syrup-agar diet at 25°C with a 12 h light/dark regimen. Wandering third instar larvae of both sexes were used for all experiments. Canton-S (CS; Bloomington Stock Center, 64349; RRID:BDSC_64349) and w1118 (Bloomington Stock Center, 3605; RRID:BDSC_3605) larvae were used as wild-type or wild-type-like controls. The homozygous viable allele of Dmca1D, AR66 (Eberl et al., 1998; obtained from Dr. D. Eberl, University of Iowa, Iowa City, IA), contains a point mutation causing a hypomorphic phenotype (Ren et al., 1998) and was used to assess the consequences of reduced Dmca1D function in all cells.
We used mosaic analysis with a repressible cell marker (MARCM; Lee, Luo, 2001) to create animals with few Dmca1D-null mutant MNs in an otherwise heterozygous background. We used the null mutant X10 allele of Dmca1D (Bloomington Stock Center, 25141; RRID:BDSC_25141) to create a fly strain containing the X10 allele behind a flippase target recognition site (P {neoFRTry+}40A;Bloomington Stock Center, 5759; RRID:BDSC_5759) to obtain b1 Ca-α1DX10 P {neoFRTry+}40A/CyO P {ActGFP}JMR1. This strain was crossed to a fly strain containing a heat shock inducible flippase and the GAL4 inhibitor GAL80 under the control of the strong tubulin promoter behind the respective flippase recognition target site (P{hsFLP}1, y1w* P{UAS-mCD8::GFP.L}Ptp4ELL4; P{tubPGAL80}LL10 P{neoFRT}40A; P{tubP-GAL4}LL7; Bloomington Stock Center, 42725; RRID:BDSC_42725). The heat shock inducible flippase was activated by placing embryos 2 h after egg laying for 30 min from 25 to 34°C. Egg laying was permitted on fresh grape juice agar for 1 h.
To restrict Dmca1D-RNAi (VDRC, RRID:SCR_013805, stock 51491; and Bloomington Stock Center, 33413, HMS00294, RRID:BDSC_33413) knockdown to two identified MNs per ventral nerve cord hemisegment, expression was driven under the control of the even-skipped promoter (RN2-GAL4, Bloomington Stock Center, 7473; RRID:BDSC_7473; Fujioka et al., 2003). RN2 expresses in MN1-Ib, which innervates muscle 1 with big type I terminal boutons and in MNISN-Is, which innervates multiple dorsal muscles via the intersegmental nerve (ISN) and has small type I terminals (Hoang, Chiba, 2001). MN1-Ib (or short MN1b) and MNISN-Is (or short MN1s) are also named by their embryonic identities as aCC and RP2, respectively (Fujioka et al., 2003). In all RNAi experiments inclusion of extra Dicer-2 (UAS-Drc2; Bloomington Stock Center, 24650; RRID:BDSC_24650; Dietzl et al., 2007) was used to increase knockdown strength as previously reported (Ryglewski et al., 2012; Hutchinson et al., 2014).
Because RN2 is a relatively weak driver in the third instar, a flippase strategy was used to express the UAS-RNAi constructs in a mosaic fashion under the control of the particularly strong actin-GAL4 promoter. In this scheme, the flippase is driven by the weak RN2 promoter to remove a stop cassette to activate the strong actin-GAL4 promoter. This resulted in strong RNAi expression in only a small subset of aCC and RPs neurons (Hartwig et al., 2008; fly strain courtesy of Dr. S. Sanyal, BIOGEN-Idec). This scheme creates genetic mosaics, so that non-GFP labeled MNs could be used as internal controls to evaluate RNAi knockdown efficacy. Although the VDRC stock 51491 has previously been reported to yield 30–70% knockdown of Dmca1D current as recorded from larval MN somata (Worrell and Levine, 2008), we found that the TRiP RNAi construct (HMS00294, RRID:BDSC_33413) yielded a stronger and more reliable knockdown when tested by somatic patch-clamp recordings or by immunolabeling. Therefore, HMS00294 (full genotype: y1 sc* v1; P{TRiP.HMS00294}attP2) was used to address the function of Dmca1D in larval MNs during locomotion. TRiP.HMS00294 has been targeted with a VALIUM vector (Vermillion-attB-loxP-Intron-UAS-MCS) to an attP-landing site on the third chromosome of Drosophila (Ni et al., 2009). As genetic controls, we used animals with the same VALIUM vector inserted at the same attP landing site in the same genetic background (Bloomington Stock Center, 35786; RRID:BDSC_35786), but without siRNAi (Ni et al., 2008).
Larval preparation.
Third instar larvae were dissected in normal saline (composition in mm: 128 NaCl, 2 KCl, 1.8 CaCl2, 4 MgCl2, 5 HEPES, and ∼35 sucrose to adjust osmolality to 300 mOsm/kg−1; pH was adjusted to 7.25 with 1 m NaOH). Larvae were pinned dorsal side up in silicone elastomere (Sylgard 184, Dow Corning) lined Petri dishes with minute pins through the mouthhooks and the tail. Animals were dissected along the dorsal midline, and the dorsal cuticle and muscles were stretched laterally and pinned down with 2 min pins on each side. After removal of the gut and esophagus, the VNC was exposed, mounted onto an upright fixed stage Zeiss Axioskop 2 FS plus fluorescence microscope, and viewed with a 40× water-immersion objective. When constantly perfused with saline, ∼20% of the larval preparations spontaneously produced crawling-like motor patterns as characterized by peristaltic waves of rhythmic contractions, which propagate from posterior to anterior (Fox et al., 2006; Heckscher et al., 2012).
Electrophysiology.
Extracellular nerve root recordings were conducted as recently described (Kadas et al., 2015). To facilitate access to RP2 and aCC MNs with the recording electrode, the ganglionic sheath was digested by focally applying 1% protease (Protease type XIV from Streptomyces griseus; Sigma-Aldrich, catalog #P5147; CAS 9036-06-0) in normal saline with a large patch pipette (1 MΩ tip resistance). Recording patch pipettes with a tip resistance of 6–8 MΩ were pulled from borosilicate glass (1.5 mm outer diameter, 1.0 mm inner diameter, without filament, World Precision Instruments) with a vertical pipette puller (PC-10, Narishige). After protease treatment, the preparation was rinsed with 5 ml of normal saline for 2 min before patching onto the MNs. For current-clamp recordings electrodes contained internal solution of the following composition (in mm): 140 Kgluconate, 2 MgCl2, 2 Mg-ATP, 10 HEPES, 1.1 EGTA, glucose to adjust osmolality to 300 mOsm/kg−1, pH 7.25. External solution was normal saline (see above). For voltage-clamp recordings of Ca2+ currents TEA-Cl (Sigma-Aldrich, catalog #T2265, CAS 56-34-8) and 4-AP (Sigma-Aldrich, catalog #275875; CAS 504-24-5) were used to block K+ currents, and TTX (300 nm; Roth Chemicals, catalog #6973.1, CAS 4368-28-9) was used to block sodium current. The internal patch solution consisted of the following (in mm): 140 CsCl, 2 Mg-ATP, 0.5 CaCl2, 1.1 EGTA, 20 TEA-Br, 0.5 4-AP, 10 HEPES. The pH was adjusted to 7.24 with CsOH; osmolality was 323 mOsm/kg−1; osmolality of the respective external solution was adjusted to 320 mOsm/kg−1 with sucrose.
All recordings were conducted in a bath volume of 300 μl under constant saline flow (∼ 2 ml/min) at 24°C with an Axopatch 200B amplifier (Molecular Devices). Data were digitized with a Digidata 1322A (Molecular Devices) at a sampling rate of 10 kHz. Before approaching the cell, the offset potential was nulled manually. Upon obtaining a gigaohm seal, fast capacitive artifacts were compensated manually in patch mode at a holding potential of −70 mV. Whole-cell configuration was achieved with gentle negative pressure. Recordings were performed with serial resistances between 10 and 15 MΩ. For voltage-clamp recordings whole-cell capacitance was determined and compensated using the Axopatch 200B C-slow dial. Voltage errors caused by series resistance were predicted at 90%, and compensated for by 40–50% at a time constant of 2 μs. In voltage-clamp mode input resistance was calculated from the linear slope of the I–V relationship at subthreshold command potentials and subtracted off-line. In current-clamp input resistance was calculated identically but not corrected. Ca2+ currents were induced by voltage steps from −90 to +20 mV in 10 mV increments from a holding potential of −90 mV.
Intrinsic excitability of MNs was determined in current-clamp mode by somatic square-pulse current injections. At a resting membrane potential of −60 mV square-pulse currents of 400 ms duration were applied in 10 pA increments from −40 to +90 pA.
MN firing patterns during crawling-like movements were recorded in current-clamp in gap-free acquisition mode. Larval movements were monitored visually. Larval locomotion was not induced but all recordings shown were spontaneously occurring motor patterns. All electrophysiological data were acquired and analyzed with PClamp 10.4 software (Molecular Devices, RRID:SCR_011323).
MN postsynaptic responses to nicotinic acetylcholine receptor activation were recorded with somatic whole-cell recordings in current-clamp mode. The receptor agonist nicotine (2 × 10−5 m; Sigma-Aldrich, catalog #N3876; CAS, 54-11-5) was pressure applied focally to MN dendrites with a sharp microelectrode connected to a Picospritzer II (General Valve). Pulse duration was 3 ms and pressure was set to 20 psi.
Calcium imaging.
The genetically encoded calcium indicator GCaMP6s (20xUAS-GCaMP6s attP40; Bloomington Stock Center, 42746; RRID:BDSC_42746) was expressed in MNs under the control of RN2-GAL4. An Orca Flash 4.0 LT CMOS camera (C11440-42U; Hamamatsu Photonics) with HOKAWO 2.10 software was used for image acquisition. Exposure times were between 75 and 100 ms. Image series were streamed. Raw data were exported to Microsoft Excel 2010, and ΔF/F was calculated as previously described (Duch and Levine, 2002; Ryglewski et al., 2017).
Immunocytochemistry.
For immunostainings of GFP, the active zone marker bruchpilot (brp) and Dmca1D channels, larvae were dissected in saline, fixed in paraformaldehyde (4% in 0.1 m PBS) for 45 min, washed six times 20 min in PBS (0.1 m), and 4 times 30 min in PBS-TritonX (0.5%). Then preparations were incubated with primary antibodies (rabbit anti-GFP, 1:400, Invitrogen, A11122, RRID:AB_221569; mouse anti-Brp, nc82, DSHB, 1:200, RRID:AB_2314867; goat anti-Dmca1D, 1:400, Santa Cruz Biotechnology, sc-32083, RRID:AB_653056) overnight at 4°C in PBS-TritonX (0.1%). Specificity of the antibody for Drosophila Dmca1D channels was tested by Western blotting (Fig. 1C), immunocytochemistry on larval muscle (Fig. 1D), Dmca1DX10 mutant versus control embryos (Fig. 1Di), control compared with Dmca1D RNAi knockdown MN somata (Fig. 1E,F), and Dmca1DX10 mutant MNs in genetic mosaics (Fig. 1G,H). Following incubation with primary ABs preparations were washed six times 30 min in PBS and incubated in secondary ABs (donkey anti-rabbit AlexaFluor 488, Jackson ImmunoResearch, catalog #711-546-152, RRID:AB_2340619; donkey anti-mouse AlexaFluor 555, ThermoFisher, catalog #A-31570, RRID:AB_2536180; donkey anti-goat AlexaFluor 647, Jackson ImmunoResearch, catalog #705-605-147, RRID:AB_2340437; all 1:1000) overnight at 4°C. Then preparations were washed six times 30 min in PBS, dehydrated in an ascending ethanol series (50, 70, 90, 100%, 10 min each), cleared, and mounted in methylsalicylate.
Western blotting.
Flies were anesthetized on ice for 10 min. Then brains were removed quickly (<30 s per brain) and immediately homogenized in 85 μl sample buffer (2× concentrated: 125 mm 4× Tris-HCl/SDS, pH 6.8, 20% glycerol, 4% SDS, 20 mm dithiothreitol, 0.001% bromophenol blue (all Bio-Rad) in ddH2O). Samples were boiled for 3 min, spinned down for 30 s at 13,000 rpm, and stored at −28°C until used. Upon usage, samples were boiled for 3 min directly out of the freezer, spinned down for 1 min at 13,000 rpm and immediately loaded into gel pockets. The SDS-PAGE gel consisted of a 4% stacking gel, pH 6.8, and a 5% running gel, pH 8.8. Electrophoresis buffer consisted of: 25 mm Tris base, 190 mm glycine, and 0.1% SDS (all Bio-Rad) in ddH2O. Pockets were loaded with 85 μl samples containing 10 brains per lane. Band size was determined using Spectra Multicolor High Range Protein Ladder (Life Technologies, 26625), 40 μl per lane. The electrophoresis was run at room temperature (∼22°C). Large vertical electrophoresis chambers (40 ml, 15 × 100 μl pockets, with 1.5 mm thick spacers, Hoefer) guaranteed good separation of protein bands. Proteins were run through the stacking gel at a constant current of 20 mA (∼1.5 h). Upon reaching the running gel, current was increased to 30 mA until the bromophenol blue front reached the bottom of the gel (∼5 h). Proteins were transferred onto nitrocellulose membrane (Bio-Rad) overnight at 4°C at a constant voltage of 35 V using transfer buffer (25 mm Tris base, 192 mm glycine, both Bio-Rad), 15% methanol (Sigma-Aldrich) in ddH2O. The next day the nitrocellulose membrane was washed in ddH2O for 10 min. Nitrocellulose was washed 3 × 20 min in Tris-buffered saline with Tween 20 (TBST; 10 mm Tris pH 7.5, Biorad; 150 mm NaCl, 0.1% Tween 20, both Sigma-Aldrich), in ddH2O, then blocked for 2 h at room temperature in 10% milk powder (Sigma-Aldrich) in TBST, and washed again in TBST 3 × 20 min. Antibodies were prepared in 2.5% milk powder in TBST and not mixed to avoid cross-reaction. As loading control polyclonal rabbit α-heatshock protein 90 (hsp 90) antibody was used 1:1000 (Cell Signaling Technology, catalog #4874; RRID:AB_2121214). Polyclonal rabbit α-GFP IgG antibody (Invitrogen, A-11122; RRID:AB_221569) was used 1:400 to detect GFP-tagged Drosophila Dmca1A voltage gated Ca2+ channel (∼210–240 kDa), goat α-N-type Ca2+ channel α1B antibody (L-17; sc-32083, Santa Cruz Biotechnology; RRID:AB_653056) was used 1:400 to detect L-type like Drosophila Dmca1D voltage gated Ca2+ channel (∼270–280 kDa, dependent on splice variants). The nitrocellulose membrane was cut horizontally with a razor blade to separate the part of the membrane where the loading control was expected (at 84 kDa) from the rest of the gel to avoid antibody cross-reaction. Furthermore, lanes were separated vertically with a razor blade to incubate each lane separately with the desired antibody. Incubation was done in 14 ml Falcon tubes at 4°C overnight on a 3D-shaker. Then membranes were washed in TBST for 3 × 20 min. Secondary antibodies were goat α-rabbit IgG antibody, (H+L) HRP conjugate (Millipore, AP307P; RRID:AB_92641) and rabbit α-goat IgG antibody, HRP conjugate (Millipore, AP106P; RRID:AB_92411), both used at 1:5000 for 2 h at room temperature, shaking. Afterward the membrane was rinsed five times with TBST and washed three times for 20 min with TBST. Band detection was done after 5 min incubation using Immobilon Western Chemiluminescent HRP Substrate (Millipore, WBKLS0500). For detection with a chemiluminescence detection system (Fusion-SL 4.2) and Fusion Bio-1D software (both PEQLAB Biotechnologies) all pieces were reassembled in a transparent sheet protector so that the original membrane was restored.
Image acquisition and processing.
Immunolabeled whole-mount preparations were scanned with a Leica SP8 confocal laser scanning microscope with a 40 × oil-immersion lens (NA 1.2). Confocal setting for AlexaFluor 488, 555, and 647 dyes were as follows: argon laser excitation at 488 and detection between 500 and 530 nm. Solid-state laser excitation at 561 nm and detection between 570 and 600 nm. Helium-neon laser excitation at 633 and detection between 640 and 690 nm. AlexaFluor 488 and 647 were excited and scanned simultaneously, but AlexaFluor 555 was scanned separately to avoid signal crosstalk (SP8 sequential scan mode). Laser and detector settings were switched after each frame. Images were processed with AMIRA4.1 (RRID:SCR_014305), ImageJ (RRID:SCR_003070), and Corel Draw X7 (RRID:SCR_014235) software.
Experimental design and statistical analysis.
Analysis was conducted with Microsoft Excel 2010 and SPSS Statistics 22 (RRID:SCR_002865). Using the Shapiro–Wilk test, datasets were considered normally distributed at p > 0.1 and presented as mean and SD. Parametric testing was used for statistical comparison of experimental groups (Students t test for comparison of two groups, and ANOVA with LSD and/or Newman–Keuls post hoc testing for comparisons of >2 groups). If the criterion for normal distribution was not met, data were presented as boxplots (medians and quartiles in boxes and 10% and 90% values as error bars). Kruskal–Wallis ANOVA was used for statistical analysis of nonparametric data. Mann and Whitney U test was used for pairwise comparison of two groups, and for post hoc testing following detection of statistical significance with the Kruskal–Wallis test. Significance was accepted at *p < 0.05, **p < 0.01, ***p < 0.001. Sample size was not computed a priori because no information on expected effect sizes was available, and we did not use vertebrate animals. Instead, for key findings we report the level of significance, the variance, and the numbers of experiments, which together provide information on statistical power. In some cases we also provide information about effect sizes. Effect size for normally distributed data was calculated as Cohen's d (difference between the means of the experimental and the control group divided by the pooled SD). For nonparametric data effect size r was calculated by dividing the standard score by the square root of sample size (Cohen, 1988).
Results
Dmca1D channels localize along MN axons and to patches of MN dendrites
Dmca1D underlies L-type like Ca2+ current that can be recorded from the somata of larval RP2 (MNIs) and aCC (MN1b) MNs (Worrell and Levine, 2008). However, the localization of Dmca1D channels has not been tested directly in invertebrate neurons. Before analyzing Dmca1D localization in third instar larval MNs with extensively branched dendrites (∼3000 μm total length; Fig. 1A,B), we confirmed specificity of a commercially available polyclonal antibody (see Materials and Methods) for Dmca1D channels by Western blotting (Fig. 1C). Brain homogenate from wild-type flies yielded a single band at the predicted size of Dmca1D protein (274 kDa; Fig. 1C, α-Dmca1D, left lane). No bands were detected at the predicted sizes for the other two Drosophila VGCCs, namely DmαG (T-type like, 350 kDa) and Dmca1A (P/Q, N-type like, 210 kDa; Ryglewski et al., 2012). Moreover, the single band at the predicted size of Dmca1D (274 kDa) was also present in brain homogenate from DmαG-null mutants (Fig. 1C, α-Dmca1D, right lane). A knock-in scheme, pan-neuronal GFP-tagged UAS-Dmca1A transgene expression in a Dmca1A-null mutant background (Kawasaki et al., 2002), excluded possible detection of Dmca1A. Brain homogenate from this group was run separately with putative α-Dmca1D and α-GFP AB. Again, the putative α-Dmca1D AB yielded a single band at 274 kDa (Fig. 1C, α-Dmca1D, middle lane). By contrast, α-GFP yielded a signal at the predicted weight for Dmca1A (∼210–240 kDa; Fig. 1C, α-GFP, left lane, asterisk), but no signal in brain homogenate from Canton S controls. This excluded detection of Dmca1A channels, although a previous study suggested this AB to detect Dmca1A (Astorga et al., 2012).
Four additional lines of evidence for AB specificity were provided by immunocytochemistry. First, this AB yielded prominent immunolabel in larval muscle (Fig. 1D), which is known to express Dmca1D Ca2+ channels (Gielow et al., 1995; Singh and Wu, 1999). Second, the AB yielded immunopositive label in the CNS of heterozygous Dmca1DX10 mutant embryos, which were balanced over CyO P{Act-GFP}, and thus also GFP-positive (Fig. 1Di, left). By contrast no α-Dmca1D label and no GFP were present in homozygous Dmca1DX10 mutant embryos (Fig. 1Di, right). Third, mosaic RNAi knockdown of Dmca1D in subsets of MNs (see Materials and Methods) reduced the intensity of immunolabeling with this AB in MNs expressing Dmca1D-RNAi and GFP compared with internal control MNs without transgene expression (Fig. 1E). To rule out possible artificial effects of GFP on Dmca1D expression, we also created mosaic GFP expression in subsets of MNs without any RNAi. There, all MN somata revealed similar Dmca1D-immunopositive signal (Fig. 1F), independent of GFP expression. Finally, evidence for AB specificity was provided by the absence of label in Dmca1DX10-null mutant MNs. Because Dmca1D-null mutants are embryonic lethal, we used the MARCM technique (Lee and Luo, 2001) to create singled-out Dmca1DX10-null mutant but GFP-positive MNs in an otherwise heterozygous background (see Materials and Methods). No detectable AB label was found in GFP-positive null mutant MNs (Fig. 1G), but neighboring heterozygous neurons were labeled (Fig. 1G). By contrast, controls for mosaic GFP expression that did not carry the Dmca1DX10 mutation showed Dmca1D immunolabel in GFP-positive and in GFP-negative neurons (Fig. 1H). Together, these data provided multiple lines of evidence that this antibody was specific for Dmca1D channels.
Triple immunolabeling of larval RP2 and aCC MNs with α-GFP (RN2-GFP), neuropil regions with the active zone marker, α-Brp, and α-Dmca1D revealed localization of Dmca1D channels to multiple neuronal compartments (Fig. 2; 15 animals tested). Continuous Dmca1D label was found along all MN axons in all segmental nerve roots (Fig. 2C, white arrows). Dmca1D was also detected in axons of other neurons, which projected in CNS commissures or through segmental nerve roots (Fig. 2C, white arrowheads), indicating axonal Dmca1D localization in multiple different types of neurons. MN somata showed also immunopositive signal for Dmca1D channels (Fig. 2C, white asterisks), though it remained unclear whether somatic label reflected functional L-type channels in the membrane, or channel production. As suggested by previously published voltage-clamp recordings (Worrell and Levine, 2008), patches of Dmca1D label were also detected in MN dendrites (Fig. 2E–G, white arrowheads, asterisks for dendrite without Dmca1D). Therefore, Dmca1D channels localized along larval MN axons and clustered to some dendritic segments of the same MNs. The clustered appearance of Dmca1D to patches of MN dendrites is reminiscent of L-type channel localization in dendrites of spinal MNs (Heckman et al., 2003; Anelli et al., 2007), but the axonal localization has not been reported in vertebrate MNs.
Genetic mosaics allow probing Dmca1D function selectively in identified larval MNs
A general problem in addressing cell intrinsic ion channel function during network activity is that bath applied pharmacological blockers may affect all neurons in the network. Therefore, we used Drosophila genetics to target Dmca1D RNAi together with GFP as a reporter to subsets of larval crawling MNs (flippase strategy, see Materials and Methods). The comparison of Dmca1D immunolabel intensities in MN somata with and without RNAi (see above, 10 animals) indicated that Dmca1D-immunopositive label was reduced in MNs with GFP and RNAi expression (Fig. 1E, right), but not in internal control MNs in the same preparation (Fig. 1E, right). By contrast, in control animals with mosaic GFP but no RNAi expression, GFP-positive and GFP-negative MNs displayed similar Dmca1D immunolabel intensities (Fig. 1F; n = 10). This indicated that targeted Dmca1D RNAi knockdown to subsets of MNs was effective.
Voltage-clamp recordings confirmed that RNAi knockdown of Dmca1D in MNs reduced L-type-like current amplitude by ∼70% compared with internal controls (Fig. 3A–C). The remaining current was likely not the result of compensatory upregulation of other VGCCs, because blockers for Dmca1A (PLTXII) and DmαG (amiloride; Ryglewski et al., 2012) had no effect. In genetic mosaics, non-RNAi-expressing MNs displayed similar maximum current amplitudes at 0 mV as MNs in control animals for genetic background. By contrast, RNAi-expressing MNs displayed significantly reduced L-type current amplitudes, which were similar to those in heterozygous AR66 mutants (Fig. 3C). Therefore, our flippase RNAi scheme provided a useful tool to reduce Dmca1D expression selectively in a subset of identified MNs without affecting other parts of the motor network (but see below for potential indirect effects of genetic knockdown through development).
Dmca1D reduces MN firing responses to moderate amplitude somatic current injection
Representative MN1s firing responses to somatic square-pulse current injections are depicted in Figure 3D and E. For both control (Fig. 3D) and targeted Dmca1D RNAi knockdown (Fig. 3E) the voltage traces for two current injection amplitudes (40 and 50 pA) are superimposed. For better visualization additional four voltage responses for current injection amplitudes of 10, 20, 40, and 60 pA are shown for another control MN1s (Fig. 3G) and another Dmca1D RNAi knockdown (Fig. 3H). In both control (Fig. 3D,G) and following targeted Dmca1D knockdown (Fig. 3E,H), MN1s shows a characteristic biphasic depolarization upon square-pulse current injection that results in a marked delay to the first action potential and was previously described (Choi et al., 2004; Ping et al., 2011). Both the fast depolarization to ∼−40 mV followed by a slower depolarization to firing threshold and the resulting delay to the first action potential have been attributed to Shal (Kv4) channels, because RNAi knockdown of Shal eliminates the biphasic depolarization and the delay to firing (Schaefer et al., 2010). In agreement with earlier reports (Worrell and Levine, 2008) knockdown of Dmca1D significantly increased the frequency of MN firing responses to somatic square-pulse current injections (Fig. 3D–H). This was likely caused by a reduction in the amount of Ca2+ activated K+ outward currents (Worrell and Levine, 2008). However, we found increases in MN firing responses only for current injection amplitudes (20–70 pA), which resulted in moderate firing frequencies of 10–90 Hz (Fig. 3F–H). But note that Dmca1D RNAi knockdown did not cause significant increases of MN firing frequencies in responses to somatic current injection amplitudes that evoked locomotion-like firing frequencies (Fig. 3F; >90 Hz at current amplitudes ≥80 pA; see Fig. 7G–J). In square-pulse protocols, we did not increase the amplitude of injected current >90 pA because this often resulted in unstable recordings from the small somata of larval Drosophila MNs. However, ramp current injections that resulted in peak firing frequencies >100 Hz showed decreased firing frequencies upon blockade of Dmca1D (see Fig. 7D,G–K). Therefore, a reduction of Dmca1D channels increased MN responsiveness to somatic current injections in the lower frequency range (<90 Hz), did not affect MN firing rates in medium range locomotion-like frequencies (90–110 Hz), but decreased maximum intraburst firing frequencies (>120 Hz; see Fig. 7G,I). However, somatic current injections likely did not account for MN firing responses to normal network input. Therefore, we next recorded MNs in situ during crawling-like locomotor patterns in semi-intact preparations (Figs. 4, 5).
Dmca1D is not required for gross motor coordination during crawling
Larval forward crawling is characterized by rhythmical waves of segmental MN bursts that propagate from the posterior to the anterior abdominal segments (Fox et al., 2006; Heckscher et al., 2012; Kadas et al., 2015). Crawling-like motor patterns were displayed spontaneously by larvae that were dissected along the dorsal midline (Fig. 4A; Kadas et al., 2015). Characteristic delays between motor bursts in posterior and anterior segments could either be recorded extracellularly from segmental nerve roots (Fig. 4B), which contain ∼30 MN axons (Hoang and Chiba, 2001), or observed visually under the microscope. Cycle period (onset from 1 burst to the next in the same segment) was statistically similar in wild-type controls, AR66 Dmca1D hypomorphic mutants, targeted RNAi in MNs, and genetic controls for the mutant and RNAi strains (Fig. 4C). Similarly, rhythmicity and segmental delays were not significantly different in any of these genotypes. Therefore, Dmca1D channels were not required for coordinated crawling-like motor output, or for the overall central nervous regulation of locomotion speed.
Dmca1D enhances MN excitability during locomotion
To probe L-type Ca2+ channel function during locomotion specifically in MNs we used our flippase strategy. Dmca1D RNAi was targeted to random subsets of larval RP2 (MN1s) and aCC (MN1b) MNs, to then record either RNAi knockdown or control MNs (Figs. 4D,E, 5). In few instances we were able to conduct dual patch-clamp recordings (data not shown), either from two control MNs (Fig. 4D, GFP-positive aCC MNs), or from a control and a Dmca1D RNAi knockdown aCC MN in the same segment (Fig. 4E, for mosaic expression). The recordings indicated similar intrinsic excitabilities and similar spontaneous activity patterns of sister control aCC MNs (n = 4), but different excitabilities (n = 3) and different spontaneous activity patterns (n = 1) of control versus Dmca1D knockdown aCC MNs (data not shown). However, dual in situ recordings were technically challenging, resulting in insufficient data for thorough analysis, especially because only one mosaic animal displayed typical crawling-like locomotor patterns during dual recording. Therefore, for further analysis we compared individual recordings from MNs in seven mosaic RNAi knockdown animals and seven genetic controls (Fig. 5A–I).
Representative recordings of a control (Fig. 5A) and a Dmca1D RNAi knockdown MN1s (Fig. 5B) during crawling-like locomotor patterns indicated that RNAi knockdown of L-type like Ca2+ channels reduced MN firing responses to normal CPG activity. Quantification revealed that burst and synaptic drive potential duration were significantly reduced by Dmca1D knockdown (Fig. 5C–E). By contrast, the amplitude of the synaptic drive potential remained unaltered (Fig. 5C,F). MNs with reduced L-type like Ca2+ current fired significantly fewer action potentials per burst (Fig. 5G), though the mean intraburst firing rate remained unaltered (Fig. 5H). The maximum intraburst firing rates were significantly reduced following targeted Dmca1D knockdown (Fig. 5I). Therefore, Dmca1D channels likely affect two aspects of MN firing patterns during locomotion: they increase drive potential and burst duration, and they increase maximum intraburst firing rates. This likely enhances muscle contraction power for each segmental wave of MN excitation. Similarly, in vertebrates L-type channels enhance muscle force production at a given synaptic drive to MNs (Heckman et al., 2003).
Dendritic Dmca1D channels enhance MN responses to excitatory synaptic input
In contrast to dendritic L-type channels in vertebrate MNs, Dmca1D localizes to the somatodendritic and to the axonal domains of larval Drosophila MNs (Fig. 1). We next tried to obtain insight as to whether increased MN firing responses to CPG input were caused by axonal or dendritic Dmca1D channels, or both. In vertebrate MNs L-type Ca2+ channels are thought to contribute to PIC, which enhances synaptic drive to dendrites, and thus prolongs MN bursting (Heckman et al., 2003; Hultborn et al., 2013). Given that we observed Dmca1D localization to patches of larval Drosophila MN dendrites (Fig. 2), and that Dmca1D knockdown in postsynaptic MNs shortened burst and drive potential duration (Fig. 5D,E), we hypothesized dendritic amplification of EPSPs by L-type like Ca2+ channels also in Drosophila MNs. To test this we compared the responses of MNs with normal and reduced Dmca1D current to cholinergic synaptic input (Fig. 5J–L). To circumvent potential developmental effects that may be caused by genetic manipulation through development, we combined focal pressure application of the cholinergic agonist nicotine (2 × 10−5 m) with acute pharmacological blockade of Dmca1D with lanthanum (La3+; 1 μm). To prevent indirect effects via the network, we recorded EPSPs in response to nicotine puffs in the presence of the sodium channel antagonist TTX to block action potentials. Focal nicotine application to larval MN dendrites reliably induced EPSPs that could be recorded from the MN soma and typically contained a prolonged shoulder (Fig. 5J, arrowhead). Acute blockade of Dmca1D by bath application of La3+ eliminated this shoulder (Fig. 5J, red trace) and significantly shortened the EPSP (Fig. 5J,K). This indicated that dendritic L-type channels normally prolong MN postsynaptic responses to cholinergic input.
By contrast, EPSP amplitude was not reduced (Fig. 5L). Changes in EPSP amplitude as caused by Dmca1D knockdown were likely masked by a concomitant reduction in Ca2+ activated outward K+ current, since we have previously shown that BK channels reduce EPSP amplitude in these MNs (Kadas et al., 2015), and Ca2+ activated dendritic outward conductances are known to reduce the amplitude of prolonged EPSPs in many central neurons (Faber and Sah, 2003). We further confirmed dendritic Dmca1D channel activation upon focal nicotine application by Ca2+ imaging (Fig. 6). Nicotine puffs to larval MN dendrites caused relative increases in local dendritic calcium indicator fluorescence (GCaMP6s) by 60–80% (Fig. 6A,C). Following pharmacological blockade of Dmca1D channels by bath application of La3+ dendritic calcium signals were decreased to ∼20% (Fig. 6B,C). No run down in nicotine evoked local dendritic calcium responses was observed within 5 min in saline (Fig. 6D), but significant reductions were observed already 3 min after bath application of La3+ (Fig. 6D). These results were consistent with a role of dendritically localized Dmca1D channels in boosting synaptic drive to MN dendrites, and thus prolonging MN burst duration during locomotion.
Axonal Dmca1D channels enhance maximum MN firing rates
Dendritic Dmca1D function did not explain the reduction in maximum firing rates during locomotion upon Dmca1D RNAi in MNs (Fig. 5I), because synaptic drive potential amplitude remained unaltered (Fig. 5F). Therefore, decreased maximum MN intraburst firing rates in Dmca1D knockdown were not simply a consequence of reduced depolarization amplitude. We have recently shown that BK channels encoded by slowpoke (slo) promote maximum intraburst firing rates by increasing the amplitude of a fast afterhyperpolarization (fAHP), thus promoting de-inactivation of fast Na+ channels (Kadas et al., 2015). Similarly, in many fast-spiking mammalian neurons Kv3 channels augment AP repolarization and the recovery of Na+ channels from inactivation (Baranauskas et al., 2003). Based on these reports we hypothesized that axonally localized Dmca1D channels (Fig. 2C) may activate transient BK channels, and thus, indirectly enhance Na+ channel de-inactivation during high-frequency firing. Three sets of analysis were conducted to test this hypothesis:
First, we confirmed by Ca2+ imaging that action potentials caused axonal Ca2+ influx (Fig. 7A). Retrograde electrical stimulation of MN axons with a suction electrode placed en passant on the respective motor nerve evoked distinct axonal Ca2+ signals (Fig. 7A, middle trace). The amplitude of the axonal Ca2+ signals increased with the number of action potentials until saturation at 4–5 action potentials at 100 Hz (Fig. 7A, bottom trace). MN axonal Ca2+ signals were significantly reduced upon targeted expression of Dmca1D RNAi in MNs (Fig. 7Ai). Therefore, axonal Dmca1D channels as detected by immunocytochemistry were functional, not simply a reflection of axonal transport, and mediated activity dependent axonal Ca2+ influx.
Second, if Ca2+ influx through Dmca1D channels was indeed necessary for a BK-mediated fAHP, which in turn enhanced the release of fast Na+ channels from inactivation, RNAi knockdown of Dmca1D in MNs should have similar effects on firing rates and fAHP as loss of BK. Therefore, we compared maximum MN intraburst firing frequencies and fAHP amplitudes and durations during locomotor patterns in controls (Fig. 7B, top black traces), following selective Dmca1D RNAi in MNs (Fig. 7B, middle red traces), and in slo4 mutants (Fig. 7B, bottom blue traces). Compared with wild-type controls and to internal control MNs in genetic mosaics, maximum intraburst firing rates were significantly decreased in slo mutants (Kadas et al., 2015) and following Dmca1D RNAi knockdown in MNs (Fig. 7B–D). Similarly, compared with controls, the fAHP was smaller (Fig. 7B,C,E) and slower (Fig. 7B,C,F) in slo mutants and following Dmca1D RNAi knockdown in MNs. Because Dmca1D RNAi did not eliminate all L-type current (Fig. 3A–C) reductions in fAHP amplitude and maximum firing frequency were not as pronounced as in slo mutants, but statistically significant compared with internal control MNs in mosaic animals (p < 0.01, unpaired t test) or to wild-type MNs (Fig. 7D–F). This indicated that axonal L-type channels maximized intraburst firing rate by quickly releasing Na+ channels from inactivation via a fAHP mediated by BK channel activation.
However, results from Ca2+ channel knockdown through development must be viewed with caution, because Ca2+ influx may affect multiple aspects of neuronal differentiation (Spitzer, 2006), may cause compensatory changes in the expression of other ion channels (Marder and Goaillard, 2006), and L-type channels can link neural activity and gene expression (Flavell and Greenberg, 2008). Therefore, we next tried to gather additional support for our data by acute pharmacological blockade of Dmca1D. Larval MNs were recorded in current-clamp, and ramp current injections were adjusted to mimic MN crawling bursting patterns (Fig. 7G). In this setting, acute blockade of Dmca1D by bath application of La3+ (1 μm) reduced maximum firing rates (Fig. 7H,I) and the amplitude of the fAHP (Fig. 7H,J), and it increased the duration to maximum fAHP amplitude (Fig. 7H,K). Therefore, results from both sets of experiments, permanent RNAi knockdown and acute blockade, support the hypothesis that Dmca1D channels enhance and sharpen the fAHP amplitude to increase MN maximum firing rate during locomotion.
Discussion
A conserved role for L-type Ca2+ channels in augmenting MN excitability
By combining immunocytochemistry with targeted genetic and pharmacological manipulation, calcium imaging, and in situ patch-clamp recordings during locomotor-like behavior, this study shows that L-type Ca2+ current functions in Drosophila MNs to facilitate burst durations and maximum firing rates during rhythmical crawling behavior. Similarly, in vertebrate spinal cord MN burst duration and spike numbers are increased by L-type channels (Heckman et al., 2003), thus indicating a conserved role of VGCCs in augmenting MN excitability from fruit fly larvae (this study) to mammals (Heckman et al., 2003). However, differences exist between species in how Ca2+ channels operate in conjunction with other conductances in mediating this function (Bouhadfane et al., 2013). We find the Drosophila L-type homolog Dmca1D to be localized to the axonal and dendritic compartment of MNs and propose that dendritically located channels increase burst duration whereas axonal channels may promote maximum intraburst firing rates. The proposed dendritic function resembles that of spinal MN L-type channels, whereas the axonal function has not been described before. At present it remains unclear whether axonal and dendritic L-type channels are different splice variants of Dmca1D, possibly with different kinetics. Addressing these questions by focal voltage-clamp recordings from the soma versus the axon or by protein isoform isolation and analysis from isolated central versus peripheral parts of these MNs was beyond the scope of this study.
Dendritic Dmca1D channels boost excitatory synaptic drive to MNs
Spinal MN dendrites can generate PIC, which is either mediated directly by Ca2+ influx through dendritic L-type channels (Heckman et al., 2003), or indirectly by Ca2+-activated nonselective cation currents (Ican; Perrier and Hounsgaard, 1999; Bouhadfane et al., 2013). PIC boosts synaptic drive and thus enhances MN responses to a given synaptic drive (Heckman et al., 2004). Similarly, we found that L-type channels in Drosophila MNs enhance excitatory synaptic input to dendrites. This interpretation is supported by (1) shorter synaptic drive potentials and burst durations during locomotion in MNs with selective RNAi knockdown of Dmca1D, (2) by shorter durations of excitatory postsynaptic responses to focal nicotine puffs to MN dendrites following selective RNAi knockdown of Dmca1D in MNs, and (3) by reduced dendritic calcium elevations upon focal dendritic nAChR activation following pharmacological block of Dmca1D. Together these data indicate a highly conserved function of L-type channels in enhancing excitatory synaptic input to MN dendrites. In mammals, L-type Ca2+ current enhances MN firing rates and affects the pattern of MN recruitment (Collins et al., 2001, 2002).
Finally, in spinal MNs L-type current mediated PIC is contingent upon supraspinal aminergic control that facilitates locomotion (Heckman et al., 2003; Perrier and Delgado-Lezama, 2005; Hultborn et al., 2013). Amines also facilitate locomotion in invertebrates (Horvitz et al., 1982; Brembs et al., 2007; Pflüger and Duch, 2011), including Drosophila larvae (Fox et al., 2006), again pointing to conserved molecular strategies for the regulation of motor power output, though direct aminergic modulation of MN ion channels remains to be investigated in invertebrates.
Axonal Dmca1D channels enhance MN maximum firing rates
Evidence for L-type channels in nonsynaptic axonal membrane exists also for hippocampal neurons (Tippens et al., 2008), leech neurons (Beck et al., 2001), rat optic nerve, and spinal dorsal column white matter (Brown et al., 2001, 2003; Ouardouz et al., 2003), though the precise functions remain speculative. Because activity-dependent Ca2+ influx through VGCCs has numerous functions in axonal growth and synaptogenesis, it is difficult to distinguish between developmental and coding functions of axonal Ca2+ channels (Bucher and Goaillard, 2011). Our study provides evidence for axonal Ca2+ influx through L-type channels upon MN firing, thus indicating a coding function for MN axonal Dmca1D channels. Multiple sets of experiments are in agreement with the hypothesis that MN axonal L-type channels maximize intraburst firing rate by quickly releasing Na+ channels from inactivation via a fAHP mediated by BK channel activation. First, targeted Dmca1D RNAi knockdown in subsets of MNs reduced maximum intraburst firing rates during locomotion (Fig. 5I). Second, when mimicking MN crawling burst shape by somatic current injections, maximum intraburst firing rates were reduced by acute pharmacological blockade of Dmca1D (Fig. 7G,H). Therefore, Dmca1D enhances MN excitability independent of synaptic mechanisms. Third, we have previously reported a reduction in MN maximum intraburst firing rates upon acute blockade and/or genetic knockdown of slo encoded BK channels (Kadas et al., 2015). Slo increases MN firing rates by mediating a fAHP that augments the release of fast Na+ channels from inactivation (Kadas et al., 2015). We now report that RNAi knockdown of Dmca1D in MNs mimics the effects of slo on fAHP amplitude and duration (Fig. 7B–F), thus indicating that slo functions downstream of Ca2+ influx through Dmca1D. However, given that Ca2+ influx through L-type channels may affect multiple aspects of neural differentiation (Spitzer, 2006), or alter transcriptional programs (Flavell and Greenberg, 2008), our data from RNAi knockdown of Dmca1D must be interpreted with caution. Nonetheless, we judge it unlikely that Dmca1D affects the fAHP in MNs via indirect developmental effects, because acute pharmacological blockade of Dmca1D also decreased the amplitude and increased the duration of the fAHP (Fig. 7G–J). Therefore, reduced fAHP amplitude and firing rates in MNs were unlikely caused by indirect developmental effects of Dmca1D knockdown, but rather by reduced acute slo activation because of reduced activity-dependent Ca2+ influx.
In summary, our data indicate a conserved role for L-type VGCCs in enhancing MN excitability. Although synaptic input amplification by dendritic L-type channels exists in spinal MNs (Heckman et al., 2003) and in insect MNs (this study), we found no evidence for PIC in Drosophila MNs. Given that spinal MN dendrites extend spatially by approximately a factor 10–20 further than larval Drosophila MNs (Duch and Ryglewski, 2016), we speculate that PIC might be a specialization of the immensely far reaching vertebrate MN dendrites (Heckman et al., 2003, 2004). In addition, our data indicate that Ca2+ influx through axonal L-type channels is required in Drosophila MNs to support high-frequency firing through BK channel-mediated fAHP and Na+ channel de-inactivation. Therefore, L-type channels subserve distinctly different functions in different neuronal compartments, but both functions cooperate in increasing MN firing rates during locomotion. BK channel activation by axonal Ca2+ influx through VGCCs is also required for maximum firing rates of other central neurons, including CA1 pyramidal neurons (Gu et al., 2007), indicating that similar mechanisms are used by different types of central neurons. However, this axonal mechanism does not seem to exist in vertebrate MNs and might well be a unique feature of unmyelinated invertebrate motor axons. We suggest that facilitation of MN excitability by L-type channels is an old strategy, but interactions with other conductances and strategic placement to specific motoneuronal compartments seem well adapted to the distinctly different requirements in morphologically different nervous systems of different species.
Footnotes
This work was supported by the German Research Foundation to C.D. (Du-331/6-1) and S.R. (Ry 117/3-1). We thank Dr. R. B. Levine (University of Arizona, Tucson, Arizona) for helpful comments on the paper.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Carsten Duch, Institute of Developmental Biology and Neurobiology, Johannes Gutenberg University of Mainz, Colonel-Kleinmann Weg 2, 55099 Mainz, Germany. cduch{at}uni-mainz.de