ARTICLES

Quantitative Investigation of Calcium Signals for Locomotor Pattern Generation in the Lamprey Spinal Cord

Published Online:https://doi.org/10.1152/jn.00138.2004

Abstract

Locomotor pattern generation requires the network coordination of spinal ventral horn neurons acting in concert with the oscillatory properties of individual neurons. In the spinal cord, N-methyl-d-aspartate (NMDA) activates neuronal oscillators that are believed to rely on Ca2+ entry to the cytosol through voltage-operated Ca2+ channels and synaptically activated NMDA receptors. Ca2+ signaling in lamprey ventral horn neurons thus plays a determinant role in the regulation of the intrinsic membrane properties and network synaptic interaction generating spinal locomotor neural pattern activity. We have characterized aspects of this signaling quantitatively for the first time. Resting Ca2+ concentrations were between 87 and 120 nM. Ca2+ concentration measured during fictive locomotion increased from soma to distal dendrites [from 208 ± 27 (SE) nM in the soma to 335 ± 41 nM in the proximal dendrites to 457 ± 68 nM in the distal dendrites]. We sought to determine the temporal and spatial properties of Ca2+ oscillations, imaged with Ca2+-sensitive dyes and correlated with fluctuations in membrane potential, during lamprey fictive locomotion. The Ca2+ signals recorded in the dendrites showed a great deal of spatial heterogeneity. Rapid changes in Ca2+-induced fluorescence coincided with action potentials, which initiated significant Ca2+ transients distributed throughout the neurons. Ca2+ entry to the cytosol coincided with the depolarizing phase of the locomotor rhythm. During fictive locomotion, larger Ca2+ oscillations were recorded in dendrites compared with somata in motoneurons and premotor interneurons. Ca2+ fluctuations were barely detected with dyes of lower affinity providing alternative empirical evidence that Ca2+ responses are limited to hundreds of nanomolars during fictive locomotion.

INTRODUCTION

Central pattern generators (CPGs) are networks of neurons the activity of which underlies many vital behaviors like digestion, circulation, respiration, and locomotion (Arshavsky et al. 1997; Calabrese 1995; Grillner et al. 2001; Marder 2000; Selverston 1993). The interplay of intrinsic properties of neurons and their synaptic interactions determines the operation of these networks (Arshavsky et al. 1997; Grillner and Wallén 1985, 2002). Both network connectivity and neuronal intrinsic properties depend on Ca2+ signaling (Bootman et al. 2002; Ghosh and Greenberg 1995; Zucker 1999). Ca2+ entry not only determines synaptic transmitter release, playing a fundamental role in communication between cells, but also is a critical factor that affects intrinsic membrane properties. A deeper understanding of CPGs requires knowledge of how molecular components determine the properties of individual neurons, and consequently, the operation of the network in which they are embedded (Dale and Kuenzi 1997).

The CPG that coordinates muscle contractions during lamprey swimming can be recorded in the spinal cord in vitro (Grillner et al. 1998) after the application of exogenous glutamate receptor agonists like N-methyl-d-aspartate (NMDA). This pattern, “fictive” locomotion, is defined as the presence of bursts of action potentials in the ventral roots that alternate between the left and right sides of the spinal cord at any given segment (Cohen and Wallén 1980). Several mechanisms underlie these oscillating bursts. A network of reciprocally connected excitatory and inhibitory interneurons sequentially depolarizes and hyperpolarizes the motoneuron pool (Buchanan and Grillner 1987; Ekeberg et al. 1991). Burst termination in these interneurons has been attributed to the secondary effects of Ca2+ entry through voltage-gated Ca2+ channels and NMDA receptors and to spike adaptation due to summation of spike afterhyperpolarizations by Ca2+ activation of a potassium current (IKCa) (El Manira et al. 1994). Similarly, exogenous NMDA application activates intrinsic neuronal oscillators that are believed to require Ca2+ entry through the NMDA receptor in lamprey and in other vertebrates (Hochman et al. 1994; Reith and Sillar 1998; Tell and Jean 1993; Wallén and Grillner 1987). Thus in tetrodotoxin (TTX), to inhibit intercellular coordination, NMDA receptor-mediated oscillations are recorded in lamprey ventral horn neurons (Sigvardt et al. 1985). Burst termination is also strongly modified by the influence of stretch-sensitive neurons (Viana Di Prisco et al. 1990) while several neuromodulators, serotonin, dopamine, and GABA, are thought to act on calcium channels modifying intracellular Ca2+ concentration or on the K+ currents that Ca2+ entry may activate (Kemnitz 1997; Matsushima and Grillner 1992; Matsushima et al. 1993).

Qualitative changes in intracellular Ca2+ have been very briefly described in lamprey motoneurons, dendrites, and axons during fictive swimming (Bacskai et al. 1995) but the sources, types of neurons, subcellular locations, and amplitudes of somatic and dendritic cytosolic-free Ca2+ fluctuations remain unresolved. Elevation in intracellular Ca2+ concentration may be due to Ca2+ influx through voltage-dependent Ca2+ channels, synaptic glutamate-activated channels, in particular NMDA receptors, and Ca2+ release from internal stores. Spatial heterogeneity of Ca2+ signals may have important implications for motor pattern generation. The relative locations of sources of Ca2+, and the K+ channels at which this Ca2+ acts, will impact the Ca2+ concentrations necessary for neuronal oscillation and the time course of the oscillations (Jahromi et al. 1999; Marrion and Tavalin 1998; Sah and Bekkers 1996). A better understanding of Ca2+ dynamics in CPG neurons is clearly needed (Ivanov and Calabrese 2000; Kloppenburg et al. 2000).

Although Ca2+ entry to neurons of vertebrate central pattern generators is considered a vital component of the operation of all neurons involved in this activity, we know very little of this phenomenon. We have therefore sought to determine the temporal and spatial properties of Ca2+ entry during fictive locomotion in the lamprey spinal cord, a neural system the basic components of which are believed to be retained in higher vertebrates (Burke et al. 2001; Cohen 1992; Pribe et al. 1997; Rossignol and Dubuc 1994). We wished to monitor Ca2+ transients evoked during fictive locomotion in the somata and dendrites of motoneurons and neurons of the CPG. CPG neurons are the population of interneurons that similar to motoneurons show oscillating membrane potentials during fictive locomotion. To load the neurons with Ca2+-sensitive dye, we have used techniques that are both invasive and noninvasive to these neurons. Noninvasive loading ensured that the recording technique did not interfere with the final Ca2+ concentration to validate absolute Ca2+-transient amplitudes prior to and during fictive locomotion. Invasive recording, in which dye is loaded into the neurons through a recording microelectrode, enabled correlations to be made between electrophysiological activity and Ca2+ transients. Preliminary results have been published in abstract form (Viana Di Prisco et al. 2001).

METHODS

Experiments were performed on the isolated spinal cords of both adult and larval lampreys (Petromyzon marinus and Ichthyomyzon unicuspis). The animals were anesthetized with tricaine methanesulfonate (MS-222; 100 mg/l; Sigma Chemical, St. Louis, MO), decapitated, and dissected in a cold saline solution (Ringer) of the following composition (in mM): 100 NaCl, 2.1 KCl, 2.6 CaCl2, 1.8 MgCl2, 4 glucose, 5 HEPES, adjusted to a pH of 7.60 (modified from Wickelgren 1977). The skin and muscle were removed. The spinal cord was isolated and removed from the protective meninx primitiva and placed in a cooled small-volume chamber with a cover-slip floor that is inserted onto the stage of an upright microscope (Olympus BXW50) The recording chamber was continually superfused with cold oxygenated Ringer (8–10°C). Solutions of pharmacological agents were bath-applied at a perfusion rate of ∼1 ml/min. The tissue was held down after attaching it to a platinum staple long enough to record from 15 segments of spinal cord.

Electrophysiology

Microelectrode recordings were made conventionally with thin-walled glass sharp micropipettes (20–50 MΩ) filled with 1.5 M potassium acetate to which 2.5 mM Ca2+-sensitive dye (Oregon Green 488 BAPTA-1 dextran or Oregon Green 488 BAPTA-2) was added. Ventral root activity was monitored with suction electrodes.

Labeling of neurons with Ca2+-sensitive dyes

Retrograde labeling of neurons with dextran-amine-conjugates of Ca2+-sensitive dyes has proven to be very effective in the lamprey brain stem preparation (McClellan et al. 1994; Schwartz and Alford 2000). We have adapted this technique to allow dye loading of spinal motoneurons (Takahashi and Alford 2002). Two techniques were used for loading dye into neurons:

RETROGRADE LABELING OF MOTONEURONS FROM THE MUSCLE WALL.

The animals were anesthetized. A dextran-conjugated Ca2+-sensitive dye, either 10,000 MW Oregon Green 488 BAPTA-1 dextran; Fluo-4 dextran; Fura dextran, or 10,000 MW Oregon Green 488 dextran (Molecular Probes), was injected into the body musculature (5 mM, 5 μl). The animal was then allowed to recover in a cooled aquarium. After 48 h, the animal was killed under anesthesia, and the spinal cord adjacent to the injection site was removed. Imaging revealed Ca2+-sensitive dye labeling in motor and sensory neurons of corresponding hemisegments (Fig. 1A).

FIG. 1.

FIG. 1.Dye labeling of lamprey spinal ventral horn neurons. Ai: motoneurons from a single spinal hemisegment were retrogradely labeled 48 h after injection of the trunk musculature with the dextran conjugated dye Oregon Green 488 dextran. Note that the neurons are filled to the finest dendrites. The image is a reconstruction of 100 confocal sections through a spinal hemisegment (top is the lateral edge of the spinal cord, bottom is the midline. Aii: a single confocal section through the boxed region in Ai. The ellipse surrounds an image of a pipette filled with 1 μM Oregon Green placed through the imaged section. This fluorescence intensity was used as a calibration point for the dye concentration within the surrounding dendrites. Aiii: a comparison of fluorescence intensities of pipettes with a range of dye concentrations positioned as in Aii. The dendrites in this fill had a mean dye concentration of 1.0 μM. B: neuron filled with the Ca2+-sensitive dye, Oregon Green 488 BAPTA1 dextran, and imaged in situ on the confocal microscope. The dendritic tree is similarly resolved.


INJECTION FROM A SHARP MICROELECTRODE USED TO RECORD VENTRAL HORN NEURONS.

Conventional sharp microelectrode recordings were made from ventral horn neurons under visual observation through a ×40 objective lens on an upright microscope. The electrodes were loaded with 2.5 mM Ca2+-sensitive dye, and by pressure application (200 ms, 2 kPa), rapid dye labeling of the neuronal soma was achieved. The labeled neurons are then left without manipulation for 20 min to allow labeling of the dendrites (Fig. 1B). With this method, we labeled both motoneurons (identified by axon joining a ventral root) and interneurons (identified by axonal projections ipsi or contralateral within the spinal cord). Only interneurons that showed rhythmic activity during fictive locomotion and, hence most likely to be part of the spinal CPG, were included in the study. Because we did not find differences in Ca2+ signals between CPG interneurons and motoneurons their results were pooled.

Dye calibration

Spinal motoneurons are readily labeled by retrograde transport of dextran amine-conjugated dye from the trunk musculature (Fig. 1Ai). Using this technique, it was possible to label, noninvasively, many of the motoneurons in one segment of the spinal cord such that the complete dendritic arborization was visualized. This technique avoids structural damage from microelectrode recording. However, the use of a Ca2+-sensitive dye may alter Ca2+ responses by buffering free Ca2+ in the cytosol. Estimates of endogenous Ca2+-buffering capacities vary across neuron types, however, use of 1–10 μM concentrations of fura-2 (kd = 140 nM) interferes only minimally with endogenous Ca2+ signaling where previously investigated (Neher and Augustine 1992; Tank et al. 1995). We sought to determine that dextran conjugated Oregon Green 488 BAPTA1 (kd = 566 nM) concentrations used in the lamprey were <10 μM to ensure that the dye did not significantly alter Ca2+ dynamics.

Dye concentrations after loading were calibrated by comparing the intensity of neurons loaded with an inert fluorescent dye (Oregon Green 488 dextran) after injection into the trunk musculature (identical concentration and volume to injections with Ca2+-sensitive dyes, 5 mM, 5 μl) with sample concentrations of the same dye contained within electrodes inserted into the tissue. Microelectrodes containing Oregon Green 488 dextran, at different concentrations from 0.1 to 5 μM dissolved in Ringer solution were sequentially positioned in the tissue and imaged in the same plane as the labeled dendrites and somata with a confocal microscope (n = 3 preparations). This ensured that the image of the calibration pipette and the labeled neurons were at the same optical depth in the tissue (Fig. 1Aii). The fluorescence intensity of the dendrites and soma was compared with that of dye contained in the pipette over the range of concentrations (Fig. 1Aiii). Dye concentrations measured with this technique did not exceed 10 μM. These are concentrations sufficiently low that they should not significantly alter Ca2+ responses in the cells due to problems with buffering (Neher and Augustine 1992). It is also apparent from the confocal image that the motoneurons are labeled to the full extent of their dendrites, and this labeling is resolved in live tissue.

We have demonstrated how readily ventral horn neurons may be labeled with the same dye through the recording microelectrode. In Fig. 1B, a live neuron filled with the Ca2+-sensitive dye Oregon Green 488 BAPTA1 dextran was imaged on the confocal microscope. The cell was recorded by using a sharp microelectrode containing 1 M K acetate and 1 mM Ca2+ dye. Five pressure pulses (200 ms, 1.5 kPa) were applied to the interior of the pipette to label the neuron, and the image was taken after allowing 15 min for dye diffusion to occur. The injection caused no decrement in input impedance or membrane potential.

Imaging

Most imaging was accomplished with a CCD system (Hamamatsu ORCA) mounted onto a compound microscope (Olympus BX50WI) equipped with a rapidly switchable Xenon source (Sutter DG4). For analysis of nonratiometric dyes (e.g., Oregon Green BAPTA-1 or Fluo-4 dextran), fluorescence intensities after background subtraction were normalized to the baseline (prestimulus), giving a baseline value of ΔF/F = 0. Baseline was defined as the fluorescence at the location analyzed prior to initiation of fictive swimming for retrogradely filled neurons. For neurons filled through the microelectrode, this was defined as the minimum fluorescence during fictive locomotion. For ratiometric imaging, we used dextran-conjugated Fura (Fura dextran) and Ca2+ concentrations calculated using the method of Grynkiewicz et al. (1985). In some cases, dye-filled neurons were imaged on a Biorad MRC 600 confocal microscope (Cochilla and Alford 1998; Schwartz and Alford 2000).

Calibration of Ca2+-sensitive dyes

Ca2+-fluorescence curves were created from Ca2+ standards obtained form molecular probes. The calibration solutions contained 100 mM KCl and 30 mM 3-(N-morpholino)propanesulfonic acid (MOPS), pH 7.2. To achieve a range of Ca2+ concentrations from 0 to 40 μM, 11 buffer solutions were made up comprising 0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, and 10.0 mM CaEGTA and corresponding descending concentrations of K2EGTA to maintain constant osmolarity. Fluorescence was measured on the recording microscope after the addition of 5 μM dye (Oregon 488 BAPTA 1 dextran Oregon Green 488 BAPTA2, Fluo-4 dextran, or Fura dextran to the solution). All calibrations were performed at 10°C to conform with the experimental recording temperature.

RESULTS

Ca2+ oscillations recorded during fictive locomotion

Ventral horn interneurons and motoneurons were labeled by direct injection into the neurons as described in the preceding text (Figs. 1B and 2C; n = 9 cells). Kinetic data that resolved Ca2+ oscillations during fictive locomotion were obtained using the nonratiometric dye Oregon Green 488 BAPTA1 dextran. These data could then be used to relate dynamic changes in Ca2+ concentration to electrophysiological measures of depolarization and spiking.

FIG. 2.

FIG. 2.Ca2+ and membrane potential oscillations during fictive locomotion. A: a neuron, labeled by dye injection through the microelectrode (shown in C), was imaged during fictive locomotion. Regions defined as soma, proximal dendrite (Ai), and distal dendrites (Aii) were separately analyzed. With no current injection and recording in current-clamp conditions the membrane potential (mp) oscillated in phase with the ipsilateral ventral root (vr). Synchronous to this oscillation the Ca2+ dye fluorescence increased in phase with the membrane potential oscillations. The amplitude of these oscillations increased with distance from the soma. B: DC current injection (0.2–0.3 nA) caused the neuron to fire action potentials in the depolarizing phase of the locomotor-induced membrane potential oscillations. Associated with the action potentials was an increase in the amplitude of the fluorescence oscillations (Bi, soma and proximal dendrite; Bii, distal dendrite). Ci: image of the soma and proximal dendrite and (Cii) the distal dendrite from the image stack used to generate the data in (A and B). Fictive locomotion was induced by the application of 200 μM NMDA to the superfusate.


Fictive locomotion in the lamprey spinal cord comprises alternating ventral root bursting activity between opposite sides of the spinal cord. Motoneurons and premotor interneurons depolarize and may spike in phase with the bursting activity on the same side of the spinal cord. In Fig. 2, a neuron was labeled with Oregon Green 488 BAPTA 1 dextran and the fluorescence imaged at 20 Hz simultaneously to electrophysiological recording of membrane potential and of firing in the adjacent ipsilateral ventral root. Ca2+ oscillations are observed in both soma and dendrites. In all nine neurons examined, the Ca2+ concentration always rose during the depolarizing phase of activity of the recorded neuron, and the simultaneous ipsilateral ventral root burst and fell during the hyperpolarizing phase.

Lamprey ventral horn neurons show a great deal of variability in action potential firing frequency during fictive locomotion from zero to three or four action potentials per depolarizing cycle (Buchanan 2001). It is, however, possible to prevent or initiate spiking by injection of DC current through the recording microelectrode. The mean membrane potential of the recorded neuron was controlled by injection of DC current under current-clamp conditions. In Fig. 2A, the neuron was recorded from with no action potentials. Under these conditions, the amplitude of the Ca2+ signal varied considerably across different dendritic compartments. The cells were analyzed as three separate compartments as shown in Fig. 2C: the soma, the proximal dendrite (defined as the unbranched dendrite within 50 μm of the soma), and the distal dendrites (comprising branched dendritic compartments distal to the proximal dendrites ≥100 μm from the soma). The amplitude of Ca2+ oscillations measured from peak to trough increased from a low at the soma (ΔF/F = 0.035 ± 0.015) to the proximal dendrites (ΔF/F = 0.102 ± 0.038) to the distal dendrites (ΔF/F = 0.201 ± 0.066; where F was defined as the interburst fluorescence after background fluorescence was subtracted, ΔF as the difference between F and fluorescence at peak). All measurements were the means from nine neurons of the means of at least five locomotor cycles.

The neurons were depolarized with DC current injection to allow spiking in phase with locomotor-induced oscillations (Fig. 2B). The mean amplitude of the Ca2+ oscillations was significantly increased during this activity. A single action potential in each cycle increased the Ca2+ signal at the soma significantly more than at the distal dendrites (n = 4 cells recorded from 4 separate preparations; P < 0.05; soma increase of 200 ± 50%, proximal dendrite by 117 ± 34%, distal dendrite 86 ± 33%). This result was consistently seen in all cells examined. This increase in Ca2+ signal amplitude was caused by Ca2+ entry evoked during the action potential and not by a depolarization-induced alteration in synaptically driven Ca2+ entry or a reduction in Ca2+ sequestration due merely to the current injection necessary to cause action potentials. This is clear for two reasons: first, the Ca2+ signal was only increased in cycles where an action potential occurred; and second, the differential (which represents Ca2+ flux into and out of the cytosol) of the Ca2+ signal reveals rapid increases in Ca2+ concentration that coincide precisely with the timing of action potentials in the neurons (Fig. 3A). Because the fluorescence of the Ca2+-sensitive dye is a measure of cytosolic Ca2+ concentration, the differential represents Ca2+ flux to and from this cellular compartment. The precise timing of action potentials during fictive locomotion and the resulting Ca2+ transient was seen throughout the neuron, from soma to distal dendrites, and was seen in all cells examined in which action potentials were recorded.

FIG. 3.

FIG. 3.Ca2+ entry to the cytosol occurs synchronously to depolarizing events in the ventral horn neurons. A: action potentials evoke synchronized rapid Ca2+ entry events. Ai: the same trace recorded in Fig. 2Ai in which the Ca2+ signal is enhanced during action potentials. Aii: the differential of the trace in Ai. Data above the solid line indicates Ca2+ entry to the cytosolic compartment. Below the line indicates removal of Ca2+ from this compartment. Note that the large abrupt transients coincide precisely with the action potentials in Aiii. Aiii: microelectrode recording of the neuronal membrane potential. B: Ca2+ entry is also synchronized with subthreshold membrane potential oscillations. Bi: the same trace recorded in Fig. 2Bi in which the Ca2+ signal is recorded during subthreshold locomotor-induced membrane potential oscillations. Bii: the differential of the trace in Bi after 5-point box smoothing function (±60 ms) was applied to the data. Ca2+ entry is synchronized to membrane potential oscillations. Biii: microelectrode recording of the neuronal membrane potential.


The differential of the Ca2+ signal was also calculated during fictive locomotion when membrane potential oscillations were subthreshold for action potentials. To reduce noise, the graph in Fig. 3Bi was filtered with a three-point box filter spanning a sampling period of 60 ms. The differential clearly demonstrates the coincidence of cytosolic Ca2+ entry to with depolarization of the neuron. This result was similar throughout the dendritic tree of all nine neurons examined. We were unable to record significant differences in the timing relationship between depolarizations and Ca2+ signals in the neurons examined regardless of whether they were interneurons or motoneurons.

Locomotor-induced Ca2+ transients recorded in retrogradely labeled neurons

Microelectrode recording from the somata may have interfered with the amplitude of the somatic Ca2+ signal, particularly at the impalement site. This would lead to an artifactual apparent distribution of the Ca2+ transient amplitudes. To determine whether somatic recordings were altered by the recording paradigm, cells were imaged after retrograde labeling of the dye. These recordings confirmed the results using microelectrodes and reveal larger dendritic than somatic Ca2+ signals essentially the same as the microelectrode recordings (Fig. 4, A and B; and quantified in Fig. 5B). Additionally, it is difficult to normalize the Ca2+ transient amplitude to fluorescence intensities recorded prior to the application of NMDA and the initiation of fictive locomotion in neurons filled through the microelectrode. This is because slow leakage of dye from the pipette leads to a drift in the resting fluorescence intensity. This drift is not present in neurons labeled retrogradely and so the fluorescence measurements were normalized to resting conditions prior to the initiation of fictive locomotion.

FIG. 4.

FIG. 4.Calcium transients recorded in retrogradely labeled motoneurons. Ca2+ concentration changes were recorded in motoneurons retrogradely labeled with Oregon Green 488 BAPTA1. Ca2+ transients are expressed normalized to resting fluorescence recorded in the neurons prior to the initiation of fictive locomotion. Fictive locomotion was activated by the bath application of NMDA (200 μM). Ai: single time point profiles of normalized fluorescence along the region of the dendrite marked by the red box in (Aiii) and the oscillations recorded in (B). The graphs represent fluorescence along the length of the dendrite at the peak (a) and trough (b) of oscillations in the dendrite recorded during fictive locomotion. The color coding is of fluorescence intensity and is identical to that in (Aii) Aii: pseudocolor plot of Ca2+ concentration recorded along the length of the dendritic profile shown in red in (Aiii). The x axis records time with the same scale as the graphs in B (5 s in total). The y axis represents length along the dendrite from the proximal dendrite at bottom to the edge of the spinal cord at top. The length of dendrite measured is shown boxed in (Aiii). Aiii: image of the dendritic trees of 2 retrogradely neurons. The red boxed area represents the region from which Ca2+ transients were obtained during fictive locomotion. B: Ca2+ transients recorded at the locations marked by lettered arrows (w, x, y, z) in (Ai) and recorded with time. The time axis is the same as in Aii. C: Ca2+ transients recorded during fictive locomotion in the cell body of the dendrite profiled in (Aiii). Inset: the cell body from which the Ca2+ transients were recorded and a second labeled neuron. The record is from the neuron on the left outlined in red.


FIG. 5.

FIG. 5.Spatial distribution of Ca2+ fluctuations during fictive locomotion. Ai: Ca2+ oscillations (top) recorded from the distal dendrites of a motoneuron retrogradely labeled with Oregon Green 488 BAPTA 1. The neuron's Ca2+ concentration oscillated in antiphase to ventral root activity recorded from a contralateral ventral root. Locomotion was activated by application of NMDA (100 μM). Aii: ratio of dendritic fluorescence recorded as the mean of the fluorescence at the peak of the locomotor evoked Ca2+ signal divided by the mean of the minimal signal [(ΔF + F)/F; peak in distal dendrites was ∼2]. Aiii: fluorescence image of the Ca2+ dye-filled dendrites from which the data in Ai and Aii were taken. Bi: Ca2+ oscillations and simultaneous ventral root activity recorded from the somata during fictive locomotion as for A. Bii: ratio calculated over the soma as in Aii.Biii: fluorescence image of the Ca2+ dye-filled soma from which the data in Bi and Bii were taken.


Motoneurons retrogradely labeled with Oregon Green 488 BAPTA1 dextran were imaged before and after initiation of locomotion (by the application of NMDA, 100 or 200 μM, to the superfusate). The fluorescence was normalized to levels recorded in the resting neurons [resting fluorescence = 1; thus fluorescence is expressed as (ΔF + F)/F]. In the neuron labeled in Fig. 4 the normalized fluorescence recorded along the length of the dendrite outlined in red (Fig. 4Aiii) was displayed versus time in a pseudocolor plot. The amplitudes of the peak fluorescence and trough fluorescence was compared along the length of the dendrite and showed great variance in both these minima and maxima (Fig. 4Ai). In addition, the fluorescence measurements recorded at different locations along the dendritic tree were recorded with time (Fig. 4B). Soma Ca2+ fluxes were also recorded over a separate time course (Fig. 4C), demonstrating much smaller concentration changes than those recorded in the dendrites. As for data recorded following microelectrode injection with Ca2+-sensitive dye there is an increase in Ca2+ signal amplitude from soma to distal dendrite.

To characterize the spatial distribution of Ca2+ concentration changes during fictive locomotion, images of fluorescence changes in soma and dendrites were also examined by generating ratio images of data from the troughs and peaks of Ca2+ measurements during fictive locomotion. Motoneurons were labeled with dextran-conjugated Oregon Green BAPTA 1 and imaged during fictive locomotion initiated by bath application of NMDA (100 or 200 μM) and monitored by ventral root recording (Fig. 5, Ai and Bi). This is the same pool of neurons analyzed in Fig. 4; n = 9). False color images presented in Fig. 5, Aii and Bii, were calculated from the ratio of mean fluorescence at the peak and the trough of the oscillations shown in Fig. 5, Ai and Bi, respectively, after subtraction of background fluorescence. The background was masked after identification of fluorescent structures in the images shown in Fig. 5, Aiii and Biii, with a thresholding function. It is clear from this result that substantial variation is seen in the Ca2+ concentration changes evoked by fictive locomotion recorded across the dendritic tree while substantially lower Ca2+ concentration changes are seen in the soma and proximal dendrite. This variation in Ca2+ transient amplitude across the dendritic tree held true for all nine neurons examined and did not correlate to small shifts in focal plane.

Calibration of Ca2+ concentrations

To obtain quantitative data about changes in Ca2+ concentration during the oscillations, we required a measurement of the Ca2+ concentration at rest and an estimate of the concentration ranges seen during fictive locomotion. We have used the fura-2 ratiometric method described by Grynkiewicz et al. (1985; see also Neher 1995). Both techniques described in the preceding text (see methods) were used to load ratiometric dye into spinal neurons. Fura dextran was retrogradely loaded into motoneurons by in vivo injections (5–10 μl, 5 mM) into the muscle wall. Ca2+ signals were recorded with excitation wavelengths at 340 and 380 nm in a total of seven neurons in five preparations. A further two neurons were similarly recorded after microinjection with Fura dextran. Resting Ca2+ concentrations were calculated from images obtained with 340- and 380-nm excitation on a standard Ca2+ calibration buffer kit (Molecular Probes) using the same batch of Fura dye and at the same temperature as the physiological measurements (10°C). Resting Ca2+ concentrations (Fig. 6A, unfilled bars) were fairly uniform across the recorded neurons with a mean Ca2+ concentration value of 100 ± 15 nM.

FIG. 6.

FIG. 6.Measurements of Ca2+ concentration in various neuronal compartments during fictive locomotion. A: histogram of Ca2+ concentration obtained from fura dextran ratios measured with excitation at 340 and 380 nm. □, concentrations calculated from neurons under resting conditions; ▦, mean concentrations calculated during fictive locomotion activated by the application of 200 μM NMDA to the superfusate. B: histogram of fluorescence recorded with Oregon green 488 BAPTA1 dextran during fictive locomotion. The fluorescence is recorded as a ratio of fluorescence during fictive locomotion to that during resting conditions. The Ca2+ concentration values on the right-hand axis are calculated from Ca2+ calibration curves obtained for the nonratiometric dyes and the resting Ca2+ concentrations obtained with fura dextran shown in (A).


Ideally Ca2+ concentrations should be calibrated in vitro in the neurons investigated. This is typically achieved by the application of an ionophore to the preparation to permeabilize the preparation to allow maximal and minimal Ca2+ binding to be determined. This technique is not feasible in the intact spinal cord because ionophores (e.g., ionomycin) do not penetrate the tissue well. Instead we buffered intracellular Ca2+ throughout the spinal cord by the application of the membrane-permeable Ca2+ buffer BAPTA-AM. Baseline fluorescence ratios were determined, under these circumstances, for comparison with those obtained with the Ca2+ calibration kit. Ratios of emission intensities after excitation at 340 and 380 nm were 0.192 ± 0.03 in neurons in which Ca2+ was buffered with BAPTA-AM and 0.196 in the Ca2+ calibration kit at 0 Ca2+.

Fura dextran recording does not give sufficient temporal resolution in this preparation to investigate dynamic fluctuations in Ca2+ concentrations during locomotor activity. At dye concentrations sufficiently low to ensure minimal buffering of endogenous Ca2+, insufficient fluorescent intensity is resolved to image significant detail with a frame rate more rapid than 5 Hz. However, mean Ca2+ concentrations were recorded in various neuronal compartments during resting conditions and then during fictive locomotion evoked by the application of NMDA (100 μM) to the superfusate. The nine retrogradely labeled neurons were analyzed by three compartments for this study, the soma, the proximal dendrite (defined as the unbranched dendrite within 50 μm of the soma), and the distal dendrites (distal to the primary branches ≥100 μm from the soma). Mean Ca2+ concentration measured during fictive locomotion increased from soma to distal dendrites (Fig. 6A) similar to the results obtained with the nonratiometric dye. The Ca2+ concentration ranged from 208 ± 27 nM in the soma to 335 ± 41 nM in the proximal dendrites to 457 ± 68 nM in the distal dendrites. Mean concentration was measured by calculating the fluorescence ratio obtained by excitation of Fura dextran at 340 and 380 nm over ≥10 frames recorded at 5 Hz during fictive locomotion. Ca2+ concentrations were calculated from standard curves obtained as described in methods.

Peak and trough Ca2+ concentrations recorded during fictive locomotion

Estimation of absolute Ca2+ concentrations with nonratiometric dyes is not straightforward. However, approximate Ca2+ concentrations may be estimated using resting Ca2+ concentrations measured with Fura dextran (mean resting concentration ranged from 87 to 120 nM across the three identified neuronal compartments). The affinity of the Ca2+ dye (Oregon Green 488 BAPTA1 dextran to Ca2+ was obtained from standard buffered Ca2+ solutions (Molecular Probes) to which dye was added and Ca2+ concentration/fluorescence curves measured at physiological temperature (10°C). The affinity obtained in this way was somewhat lower than published data from the supplier (we calculated the affinity as 566 nM, as compared with 464 nM indicated by Molecular Probes at room temperature. This value varies considerably between batches and with temperature) and a sixfold increase in fluorescence between 0 μM Ca2+ bound to a saturating concentration of 40 μM (Fmax/Fmin = 6) was observed. Using the resting Ca2+ concentrations calculated from Fura dextran measurements, the kd of Oregon Green 488 BAPTA1 dextran at 10°C, and the fluorescence measured in control conditions and during fictive locomotion, estimates were made of Ca2+ transient concentrations during fictive locomotion (Fig. 5B) using the following equation

(1)

Given the ratios of Fmax/Fmin = 6, and Fr/Fmin = 1.9 (calculated from the resting Ca2+ concentration obtained with Fura dextran data), and the calibration curve obtained for Oregon Green 488 BAPTA1 dextran, we can calculate [Ca2+] (in nM) during locomotion as

(2)

Where Ft/Fr = ratio of transient fluorescence during fictive locomotion (Ft) to resting fluorescence (Fr) prior to the addition of NMDA to the superfusate. The Ca2+ concentrations estimated using this method during fictive locomotion are similar to those obtained without measuring kinetic information with Fura dextran. All neuronal compartments show Ca2+ oscillations during swimming with a simultaneous increase in baseline Ca2+. The peak values of the measure (ΔF + F)/F (where F = baseline fluorescence prior to the addition of NMDA, ΔF + F = fluorescence measured during fictive locomotion) were at peak 1.19 ± 0.07, 1.34 ± 0.17 and 1.86 ± 0.13 in soma, proximal dendrites and distal dendrites respectively and at minimum in the interburst period 1.15 ± 0.06, 1.16 ± 0.06 and 1.47 ± 0.05 in the same neuronal compartments. These data were calculated from responses during locomotion in 5 neurons from 5 preparations and are shown graphically in Fig. 6B. The left hand vertical axis of that graph shows the equivalent calculated value of Ca2+ concentration given the above equations.

Imaging Ca2+ with low affinity dyes

To provide an alternative estimate of the magnitude of Ca2+ oscillations during fictive locomotion, two dyes with significantly lower affinities were used. Low affinity Fluo 4 dextran has an affinity of approximately 3 μM at room temperature. At 10°C its affinity exceeded 5 μM. Neurons (n = 3) were filled with this dye both by retrograde labeling and by direct injection through the recording microelectrode. No Ca2+ oscillations or increases in baseline fluorescence were recorded with this dye. A further 5 neurons were labeled with an intermediate affinity dye (Oregon Green 488 BAPTA2). Calibration curves obtained at 10°C indicated an affinity of 1.3 μM and Fmin/Fmax =18. From this data and the calibration curve the Ca2+ oscillation amplitude during fictive locomotion was predicted to lead to a change in fluorescence of 0.1 ΔF/F. In two of the four recorded neurons a noticeable Ca2+ oscillation was recorded during locomotion (Fig. 7) The mean oscillation amplitude recorded from peak (amplitude during the burst) divided by the trough (amplitude during the interburst) from all five neurons was 1.1 ± 0.1 (ΔF + F)/F. This data were not significantly different from the estimate given by higher affinity dyes supporting the results of Ca2+ oscillation amplitude data obtained with those higher affinity dyes and indicating that Ca2+ did not reach a sufficiently high concentration to saturate those dyes.

FIG. 7.

FIG. 7.Ca2+ transients measured with a low-affinity Ca2+-sensitive dye. A: neuron injected with dye (Oregon Green 488 BAPTA2) through the recording microelectrode. The contrast is enhanced over the distal dendritic region to allow visualization. Bi: Ca2+ oscillations normalized to the inter-burst minimum in the distal dendrites during fictive locomotion. Bii: Ca2+ oscillations recorded in the proximal dendrites. C: Membrane potential of the recorded neuron.


DISCUSSION

Baseline and transient Ca2+ levels

Ca2+ oscillations recorded in neurons represent two phases of activity. On the one hand, rising Ca2+ concentrations represent Ca2+ entry into the cytosolic compartment (in which the Ca2+-sensitive dye is localized) from the extracellular fluid, and from intracellular organelles that function as Ca2+ releasable stores. On the other hand, falling Ca2+ concentrations represent removal of Ca2+ from the cytosolic compartment. Surprisingly little information is available on the kinetics of Ca2+ signaling in neurons that comprise vertebrate central pattern generators (Backsai et al. 1995; Lev-Tov and O'Donovan 1995).

Here we have loaded lamprey spinal motoneurons and interneurons with Ca2+-sensitive dyes that enable high-speed imaging of physiological Ca2+ signaling during locomotor activity. We have measured the magnitude of Ca2+ responses and their baseline concentration. Resting Ca2+ concentrations were found to be between 87 and 120 nM. Mean Ca2+ concentrations recorded during fictive locomotion increased to a high of 450 nM in the distal dendrites. Our absolute estimated value for Ca2+ oscillations during rhythmic activity, are in basic agreement with other reports. For instance Ca2+ concentration increases of 50–200 nM have been reported in mouse hypoglossal motoneurons (Ladewig and Keller 2000). It has been estimated that during locomotor activity recovery of Ca2+ transients must occur on a time scale of tens of milliseconds to avoid a potentially excitotoxic accumulation of basal Ca2+ levels (Palecek et al. 1999). At least in the short term, the endogenous buffering allows motoneurons to rapidly recover from Ca2+ transients at a low energy cost (Palecek et al. 1999). Recent evidence shows that lamprey spinal cord putative CPG neurons and motoneurons are indeed rich in Ca2+ buffering molecules like calbindin and calretinin (Megias et al. 2003).

Spatial heterogeneity of Ca2+ signals due to entry through NMDA channels and other sources

Traditionally much attention has been focused on temporal aspects on membrane potential and spiking activity in the study of rhythmic neural networks, but recently several investigations have aimed at the problem of how synaptic integration from dendritically distributed inputs is conveyed to the soma, and at the functional dynamics of the spatial dendritic branching. For example, modeling studies have shown that bistable oscillatory changes in the dendritic arbor generate different charge transfer to the soma under tonic activation of NMDA conductances (Korogod et al. 2002). Thus here we have investigated the spatial distribution of Ca2+ responses in lamprey ventral horn spinal neurons. We have found that larger Ca2+ concentration changes take place in distal and proximal dendrites, whereas smaller Ca2+ oscillations are present at the soma. This is not at all unexpected because there is evidence that excitatory synaptic contacts and thus glutamate receptors are predominantly located in the dendrites rather than the soma (Moore et al. 1995, 1999).

The spatial distribution of Ca2+ currents likely play an important role in determining synaptic integration and activity control in CPG neurons and motoneurons. Not only presynaptic inputs are distributed along different dendritic compartments but also intrinsic conductances can display a variable spatial distribution over the somato-dendritic membrane. The spatial segregation of Ca2+ entry in dendritic regions and its coupling to Ca2+ -dependent ionic conductances is known to play a key role in determining patterns of neural activity (Hallworth et al. 2003). The intrinsic conductances, in particular Ca2+ -dependent ones, shape the dendritic response to receptor-gated conductances, and provide the dendrites with a dynamic way of regulating integration and spike timing critical for the operation of CPG neurons. Although our amplitude measurements indicated that Ca2+ transients do not exceed 1 μM, we cannot exclude the possibility that mechanisms important for locomotion are activated by very local changes in Ca2+ not resolved by visible light imaging.

Timing of Ca2+ signals and the locomotor rhythm.

In all of the neurons examined Ca2+ entry to the cytosol coincided with the depolarizing phase of the locomotor rhythm. This was true regardless of the location examined within the neuron or whether the depolarization was an action potential or a subthreshold oscillation in membrane potential.

Implications for plasticity in pattern generation

Finally, it is also important to consider that Ca2+ entry oscillations during induced locomotion may be important to other functions beyond the immediate control of the central pattern generation. Various forms of synaptic plasticity such as long-term depression (LTD) and potentiation (LTP) that affect the long-term function of CPGs (Parker and Grillner 1999, 2000; Soto-Trevino et al. 2001; Wolpaw and Tennissen 2001) may rely on Ca2+ signaling. The spatial segregation of Ca2+ signals in dendritic regions could affect the interplay of signaling mechanisms in plasticity, as shown in rat hippocampal pyramidal neurons (Nakamura et al. 2002). Furthermore, the pathway and temporal course of Ca2+ entry may be critical for the activation of different intracellular signal transduction processes (Gallin and Greenberg 1995). For instance, in the stomatogastric ganglion, changes in the network operation rely on the kinetics of Ca2+ currents (Turrigiano et al. 1995) and it has been postulated that intracellular Ca2+ levels affect gene expression and membrane conductance in an activity-dependant fashion (Liu et al. 1998). Perhaps long term effects of Ca2+ on the CPG for locomotion are altered by a similar mechanism.

In summary

Vertebrate locomotor pattern generation requires the network coordination of spinal ventral horn neurons acting in concert with their intrinsic oscillatory properties (Bianchi et al. 1995; Butt et al. 2002; Del Negro et al. 2002; Grillner et al. 2001; Marder and Thirumalai 2002). Ca2+ oscillations detected during NMDA-induced fictive swimming in lamprey spinal cord neurons by imaging an high affinity Ca2+-sensitive dye (Oregon green 488 BAPTA 1) were correlated with fluctuations in membrane potential, and instantaneous changes in Ca2+ induced fluorescence coincided with spike events. Using nonratiometric dye and ratiometric measurements with Fura 2 dextran we estimated Ca2+ fluctuations do not exceed 1 μM. Ca2+ fluctuations were barely detected with dyes of lower affinity, Oregon Green 488 BAPTA 2 (kd approximately 1 μM) and not at all with Fluo-4 dextran (kd approximately 3 μM), providing alternative empirical evidence that Ca2+ responses are limited to hundreds of nm during fictive locomotion.

GRANTS

This work was supported by National Institute of Neurological Disorders and Stroke Grant NS-31713 and National Science Foundation Grant IBN 0094444 to S. Alford.

FOOTNOTES

  • The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

We thank E. Schwartz, H. Photowala, T. Gerachshenko, and T. Blackmer for their helpful comments on the manuscript and for many helpful discussions.

REFERENCES

AUTHOR NOTES

  • Address for reprint requests and other correspondence: S. Alford, Dept. of Biological Sciences, University of Illinois, 840 West Taylor Street, Chicago, IL 60607 (E-mail ).