Abstract
Direct measurements of synaptic inhibition (I) and excitation (E) to spinal motoneurons can provide an important insight into the organization of premotor networks. Such measurements of flexor motoneurons participating in motor patterns in turtles have recently demonstrated strong concurrent E and I as well as stochastic membrane potentials and irregular spiking in the adult turtle spinal cord. These findings represent a departure from the widespread acceptance of feedforward reciprocal rate models for spinal motor function. The apparent discrepancy has been reviewed as an experimental artifact caused by the distortion of local networks in the transected turtle spinal cord. We tested this assumption in the current study by performing experiments to assess the integrity of motor functions in the intact spinal cord and the cord transected at segments D9/D10. Excitatory and inhibitory synaptic inputs to motoneurons were estimated during rhythmic motor activity and demonstrated primarily intense inputs that consisted of qualitatively similar mixed E/I before and after the transection. To understand this high functional resilience, we used mathematical modeling of networks with recurrent connectivity that could potentially explain the balanced E/I. Both experimental and modeling data support the concept of a locally balanced premotor network consisting of recurrent E/I connectivity, in addition to the well known reciprocal network activity. The multifaceted synaptic connections provide spinal networks with a remarkable ability to remain functional after structural divisions.
Introduction
Synaptic inputs that drive motoneurons during rhythmic motor patterns have been of great interest since the birth of modern neuroscience (Graham Brown, 1924). These inputs are essential for generating any type of movement and provide vital information about the underlying network architecture. The investigation of synaptic inputs has traditionally been performed using indirect measures such as by recording synaptic potentials during different phases of rhythmic motor activity (Shefchyk and Jordan, 1985; Cazalets et al., 1996; Gabriel et al., 2009), systemic pharmacology (Pratt and Jordan, 1987; Talpalar et al., 2011), lesioning (Cangiano and Grillner, 2005), and activity of identified premotor interneurons (Geertsen et al., 2011). Membrane depolarization/hyperpolarization has also been suggested as a measure of excitation(E)/inhibition (I) (Wallén et al., 1993; Hochman and Schmidt, 1998), especially when combined with intracellular chloride injections (Robertson and Stein, 1988).
A recent study in the adult turtle has quantified, for the first time, excitatory and inhibitory conductances during rhythmic motor patterns (Berg et al., 2007). It demonstrated concurrent waves of E/I during the on cycle, seemingly in contradiction with the half-center model. The cornerstone of the half-center model is reciprocal inhibition (Jankowska et al., 1965; Moult et al., 2013) that does not explicitly encompass the possibility of on-cycle inhibition (see Fig. 1A–E). Nevertheless, concurrent E/I during the on cycle does not contradict the concept of reciprocal inhibition, because inhibition is also detected in the off cycles. Concurrent E/I merely indicates the importance of other sources of inhibitory regulation in addition to the basic reciprocal circuitry. Interestingly, an architecture that could provide such input, known as a balanced network, is widespread across regions of the brain and could therefore represent a manifestation of general principles of network organization, which possibly has its evolutionary origin in the spinal cord.
In balanced networks, the membrane potential (Vm) is depolarizing as E increases with I, while keeping the E/I ratio constant (see Fig. 1E–G). However, the purpose of such a metabolically expensive mechanism is not yet understood. Balanced networks have been found in the brainstem (Parkis et al., 1999; Magnusson et al., 2008; de Almeida and Kirkwood, 2010), in superior colliculus (Populin, 2005), in invertebrates (Baca et al., 2008; Sasaki et al., 2009), and in the hippocampus (Csicsvari et al., 1999; Atallah and Scanziani, 2009). It is well documented in the neocortex during sensory processing (Borg-Graham et al., 1998; Haider et al., 2006), generating high-conductance states (Destexhe et al., 2003). Despite the characterization of such networks in several regions of the brain, the finding of balanced E/I was considered “heterodox” in the turtle spinal cord (Grillner and Jessell, 2009). A similar analysis of E/I drive was performed during neurochemically induced locomotor activity in the neonatal rat spinal cord, and it was reported that inhibition alternated with weak excitation (Endo and Kiehn, 2008). The discrepancy was explained as an artifact recorded in the transected turtle preparation in which the network integrity could have been structurally compromised.
The purpose of the present study was to determine whether the consequential finding of balanced E/I was indeed an experimental artifact. First, we investigated the sensitivity of the network activity to transections at segments D9/D10. If the balanced E/I observed in the original study was attributable to a structural division, then the intact spinal cord should have qualitatively different dynamics lacking balanced activity. Second, we used a network model consistent with the experiments to investigate the ability of networks to maintain their functional integrity when compromised structurally.
Materials and Methods
Synaptic drive during motor activity
It is well established that Vm has depolarizing waves during rhythmic motor activity caused by synaptic inputs and possibly aided by intrinsic properties (Stein et al., 1982; Grillner, 2006). What type of synaptic input can cause these slow waves? Traditionally, three categories were considered (Fig. 1): E driven, I driven, and reciprocal E/I. The latter possibility is appealing because of its direct relationship with the classical half-center model (Graham Brown, 1924; Jankowska et al., 1965). However, there is a fourth possibility (Berg et al., 2007), which encompasses the concurrence of E and I (Fig. 1E). In this scheme, E and I have simultaneous and opposite pull. Here, it is important to differentiate between the excitability of a neuron and the excitation/inhibition to a neuron. A neuron responding with more spikes to the same input is often misinterpreted as receiving less inhibition and/or more excitation. Nevertheless, this is not necessarily the case. A neuron can be most excitable when it receives most inhibition, as long as the excitation is sufficient to overcome this inhibition. The depolarization in Vm can be explained by the Millman equation, which is the weighted average of conductances with respect to their reversal potentials:
where 〈…〉 is the time average over a small window; Gtot, GL, Gexc, and Ginh, are the total, leak, and mean E and I conductances, respectively; and EL, Eexc, and Einh are their reversal potentials. The term balanced E and I refers to the situation where the ratio between E and I is kept constant while varying the synaptic intensity. Thus, if we define the ratio β of synaptic conductance 
When Gtot becomes large (Gtot ≫ GL), Vm is approximated by
Four possible scenarios of inhibition and excitation leading to the same dynamics in Vm. A–D, The dynamics in the membrane potential in A is driven mainly by rhythmic excitation (B), or mainly by rhythmic inhibition (C), or a combination of the two (D). E, Alternatively, Vm can be driven by balanced excitation/inhibition, where the fraction of E/I conductances (β = Gexc/Ginh) is kept constant. F, Sample intracellular current-clamp recording with negative injection to avoid spikes (Iinj = −2.3 nA) defining three states: quiescence (no synaptic input), on cycle (the spike region), and off cycle (midpoint between on cycles). G, When β is kept constant, Vm will follow the black line (in the direction of red arrows) as the total conductance is waxing and waning [Adapted from Kolind et al. (2012) with permission].
Integrated preparation
Red-eared turtles (Trachemys scripta elegans) of either sex were placed on crushed ice for 2 h to ensure hypothermic anesthesia. The turtles were killed by decapitation, and blood was substituted by perfusion with a Ringer's solution containing (in mm) 120 NaCl, 5 KCl, 15 NaHCO3, 2 MgCl2, 3 CaCl2, and 20 glucose, saturated with 98% O2 and 2% CO2 to obtain pH 7.6. The carapace containing the D4–Ca2 spinal cord segments was isolated by transverse cuts (Alaburda and Hounsgaard, 2003; Alaburda et al., 2005), and the cord was perfused with Ringer's solution through the vertebral foramen, via a steel tube and gasket pressing against the D4 vertebra. The spinal column was opened on the ventral side along D8–D10, and the dura mater was gently removed with a fine scalpel and forceps. The pia mater was opened with a longitudinal cut along the spinal cord with the bent tip of a needle. The cut was performed parallel to the ventral horn along the ventral roots. The surgical procedures comply with Danish legislation and were approved by the controlling body under the Ministry of Justice.
Nerve recording and spinal transection.
Electroneurogram recordings (ENGs) were performed with suction electrodes. The spinal cord was exposed at D9–S1, and the scratch behavior was recorded before and after transection. The transverse cut was performed at the caudal end of D10 of the spinal cord. The scratch behavior was measured by the activity of the four nerves: hip flexor, knee extensor, dD8, and hip extensor. The coherence between the nerve activities (see Fig. 3) was calculated as the trial-averaged cross-power spectrum between them, normalized by their root mean power spectra. The coherence values were selected at their dominant frequency in the power spectra. The ENGs were recorded with a differential amplifier, Iso-DAM8. The bandwidth was 100 Hz to 1 kHz.
Network activation.
Mechanical stimulation was performed with the fire polished tip of a bent glass rod. The rod was either mounted onto the membrane of a loudspeaker or controlled with a piezoelectric active stimulator. A function generator controlled the loudspeaker activation: the duration, frequency, and amplitude of the stimulus (Alaburda et al., 2005). A digital step pulse controlled the piezoelectric active stimulator: contraction duration, timing, and number of repetitions and hereby the frequency and amplitude.
Intracellular recordings.
Intracellular recordings in current-clamp mode were performed with a Multiclamp 700B and amplifiers (Molecular Devices). Glass pipettes were pulled (P-1000; Sutter Instruments) and filled with a mixture of 0.9 m potassium acetate and 0.1 m KCl. For some experiments, the electrodes also contained 4% w/v biocytin. Data were sampled at 20 kHz with a 12-bit analog-to-digital converter, displayed by means of Axoscope and Clampex software. The intracellular recordings were performed with two different methods in the transected and intact preparation. In the transected preparation, the electrodes were inserted into segments D9 and D10 from the caudal end of the transected surface, in parallel with the spinal cord. Motoneurons were directly accessible from the surface at a typical depth of 50–300 μm. In the intact preparation, the electrodes were inserted perpendicular to the spinal cord from the ventral side of D9–D10. Motoneurons were found at depths ranging from ∼300 to 800 μm. The electrodes that were inserted through the white matter had less stable resistance than those inserted through the gray matter. Therefore, the electrode resistance was more stable in the transected preparation, but in both preparations there were instances of recordings, which were stable enough for multiple trials, i.e., about 30 min.
Estimating inhibition and excitation.
All calculations were performed in Matlab. We used two methods of estimating synaptic conductances: The ohmic method and the auto-correlation function (ACF) method. Both methods rely on intracellular recordings, where Vm is considered stationary within a small window, so Gtot and 〈Vm〉 are constant. The methods differ in how Gtot is estimated. The ohmic method uses Vm paired with an injected current, Iinj (Fig. 2A), over at least two or more trials to estimate the slope of the I/V curve (Fig. 2B) (Anderson et al., 2000; Shu et al., 2003; Mariño et al., 2005). This requires that the same network state can be reproduced. In this sample data, the same motor behavior was activated by tactile stimulation (onset indicated by filled triangle). Spikes were removed by low-pass filtering Vm to give an I/V relationship as a function of time, where the slope was Gtot (indicated for two time points in Fig. 2B). Once Gtot was estimated, we could estimate excitation and inhibition for every point in time (Fig. 2C) via the following expressions:
where Ginh is given by the following:
GL, Eexc, Einh, and EL are assumed constant. The major contribution to membrane conductance is assumed synaptic, such that the total conductance is Gtot = GL + Gexc + Ginh. The resulting time-dependent Gtot (Fig. 2C, gray), Ginh (Fig. 2C, magenta), and Gexc (Fig. 2C, cyan) using Eqs. 3 and 4 are shown. Linearity should be assessed, for instance, by plotting the estimates for different depolarizations against each other (Fig. 2D). The above expressions are based on the assumption that Vm is stationary, such that the derivative is zero in a one-compartment scheme giving the following:
where the mean membrane potential, 〈Vm〉, is entered instead of Vm itself. Cm is the membrane capacitance.
The alternative method of estimating Gtot, the ACF method, is based on the time constant (τm) of the fluctuation in Vm, using the fact that Gtot = Cm/τm (Berg and Ditlevsen, 2013). The time constant is estimated via the decay in the autocorrelation function, hence the ACF method. This can conveniently be done on single-trial data for small windows. Two windows are highlighted with individual ACF decays in Figure 2E. This method is useful when experimental circumstances do not allow reproduction of same-state or different current injections. When applying both methods on the same data, there was a significant correlation between the estimates (r = 0.64, p ≪ 0.001). Part of the uncoupled variance is attributable to asymmetry in number of trials in the two methods; the ohmic estimate is based on three trials, whereas the ACF estimate is based on only one trial, making it more prone to statistical fluctuations. Thus, the ohmic method is preferable when the trial-to-trial variability is small, and the ACF method is preferable when it is large. We recommend using both measures when possible. For the intact preparation, stable recordings usable for multitrial analysis were rare. The majority of the estimation of conductance was, therefore, accomplished via the ACF method. In the few cases where the ohmic method was possible, the first five peaks of the traces were aligned around the peaks in 100 ms windows.
Estimating reversal potentials and leak
GL was estimated in each trial as the total conductance during quiescence (GL = 25–80 nS). The estimates of GL, Ginh, and Gexc rely on the simplifying assumption that leak, inhibition, and excitation are homogeneously distributed over the soma-dendritic surface and that each of these conductances can be represented by a single reversal potential, i.e., the one-compartment model. For the sample data shown (see Fig. 5), the reversal potentials were estimated to Eexc ≈ 0 mV, Einh ≈ −81 mV, and ELeak = −79 mV to minimize the presence of negative conductances. Einh was otherwise adjusted according to the action potential voltage threshold. The thresholds were determined by the first peak in the third derivative of Vm before an action potential (Henze and Buzsáki, 2001). Einh was then set to −65 mV and adjusted according to the difference between the minimum action potential threshold for the given recording and the population average (Ethreshold = −42 mV).
Synaptic fluctuations.
The fluctuations in Vm caused by synaptic inputs were estimated in both on and off cycles. In the on cycle, a spike-triggered value 18 ms before the spike occurred to diminish contamination from spikes. Multiple trials of off cycles were aligned at their center trough of Vm, and a 200 ms window was selected for processing. The SD of the high-pass-filtered window was used as a measure for the synaptic fluctuations.
Identification of motoneurons
Motoneurons were mainly identified by their location in the ventral horn, size (via Rm), and spiking relationship with nerve activity. A subset was filled with biocytin for histological processing. The tissue containing the motoneuron was carefully removed and left in PBS with 4% paraformaldehyde for 24–48 h. The tissue was then rinsed with and stored in PBS. The tissue section was mounted in an agar mount and sliced into several 100 μm slices using a DSK Microslicer (model DTK-3000W). The slices were incubated for 3–4 h at room temperature with cyanine-3 conjugated to streptavidin (1:500 or 1:250; Jackson ImmunoResearch Laboratories) in blocking buffer (PBS with 5% donkey serum and 0.3% Triton X-100). The slices were washed with PBS and incubated overnight at 5°C with primary choline acetyltransferase antibodies and goat anti-ChAT antibodies (1:500, AB144P; Millipore) diluted in blocking buffer. The slices were washed three times with PBS and incubated for 1 h at room temperature with the secondary antibody, Alexa488, conjugated to donkey anti-goat antibodies (1:1000; Jackson ImmunoResearch Laboratories) diluted in blocking buffer. After three washes with PBS, the slice was mounted and coverslipped using ProLong Gold antifade reagent (Invitrogen) and cured overnight at room temperature before microscopy. The slice was viewed using a confocal microscope (Zeiss LSM 700 with diode lasers) on an Axiolmager M2 using 10×/0.30 EC Plan-Neofluar, 40×/0.6 Corr LD Plan-Neofluar, and 63×/1.40 oil DIC Plan-Apochromat objectives (Zeiss). The fluorophore Cy3 is excited at 555 nm and detected in the range 560–1000 nm, whereas Alexa488 is excited at 488 nm and detected in the range 0–544 nm. Images were handled with ZEN 2009 software (Zeiss) in the LSM and 8-bit TIFF format.
Balanced network model
A simple model of a balanced network (see Fig. 7) was constructed, to assess the impact of structural divisions. The model network projects to motoneurons belonging to the same nucleus. Since it is known that there are important nonlocal descending inputs (Berkowitz and Stein, 1994) with persistent and excitatory activity (Currie and Stein, 1990), we assumed that there was an external drive to the local network, which was unaffected by the transection. The local network consisted of n inhibitory and n excitatory neurons with random connectivity. Each neuron had, on average, K inhibitory and K excitatory presynaptic connections. Each neuron was modeled as a leaky integrate-and-fire neuron with Vm given by the following:
where Vrest is the resting potential and Isyn,k is the synaptic current of the kth presynaptic neuron. An action potential was produced if Vm reached threshold Vthr. Vm was reset to Vrest immediately after an action potential. To simplify, the variables were rescaled so that the threshold became 1 and the resting potential was 0. The synaptic currents were modeled as double-exponential and with a strength of Jk/
Expected behavior.
The mean excitatory and inhibitory input currents to an average neuron are as follows:
where re, ri, and rext are the average firing rates for the excitatory, inhibitory and the external population, respectively. For the currents not to diverge for large K, the terms in Eqs. 7 and 8 have to balance as follows:
where 
Both rates are proportional to the external rate and independent on the number of neurons. When transecting the network and reducing the network size to a fraction α, the number of presynaptic inputs to an average neuron will decrease to K = αK′, where K′ is the original input to an average neuron. Since Kext remains the same, we arrive at the counterintuitive result that a structural division will cause higher spike rates of the remaining neurons. The rates are increased by approximately 1/
Results
Concurrence between excitation and inhibition of hip flexor motoneurons was found in a previous study where the spinal cord was transected at D10 (Berg et al., 2007). This study implemented the ohmic analysis to extract synaptic conductances from multiple trial data (Fig. 2A–D). The synaptic input was characterized by having a large rhythmic conductance, both for inhibition and excitation with high concurrence, in qualitative agreement with the balanced scheme (Fig. 1E). We repeated this analysis for cells in the intact cord. The analysis was augmented by an additional method of estimating synaptic conductances, the ACF method, which does not require multiple trial data (Fig. 2E,F). First, we performed a basic comparison between the activity in the intact and the transected spinal cord to determine the impact of the transection.
The ohmic and the ACF methods for estimating synaptic inputs illustrated on transected sample data. A, The ohmic method requires at least two trials connecting Vm and current injection (Iinj) to estimate Gtot. Shown here are three trials each with constant but different Iinj. The onset of tactile induced motor behavior is indicated (▴). B, A low-pass filtering of the traces gives an I/V relationship as a function of time (indicated for two time points) where the slope is Gtot. C, The resulting time-dependent Gtot (gray), Ginh (magenta), and Gexc (cyan) using Eqs. 3 and 4. D, Linearity should be confirmed, for instance, by plotting the estimates for different depolarizations against each other. E, The ACF method. The total conductance is estimated via the membrane time constant (GACF = C/τm), which is estimated via autocorrelation decay. F, The estimates of the total conductance using the two methods have a significant correlation (r = 0.64, p ≪ 0.001). [Adapted from Berg and Ditlevsen (2013) with permission.]
Transection
Impact on motor behavior
The activity of four motor nerves from the lumbar spinal cord enlargement (Fig. 3A) in turtles (n = 3) was recorded during tactile induced pocket scratch activity in the intact condition (Fig. 3B). Next, the spinal cord was carefully transected at the border of D10, and the motor activity was observed in response to the same tactile activated pocket scratch (Fig. 3C). The nerve activity of dD8 (Fig. 3C, top trace), the hip flexor (Fig. 3C, top middle), and the knee extensor (Fig. 3C, bottom middle) had little or no qualitative change in amplitude and cycle period after the transection. The hip extensor nerve (Fig. 3C, bottom trace) was severely affected by the transection, which is likely caused by the more caudal location (Fig. 3, compare B,C). The nerve activities were quantified as the normalized ENG spectral power, which had a gradual decrease from trial to trial before transection (Fig. 3D, cyan columns). After transection, the activity level declined, though both hip flexor and knee extensors had a transient increase in power (Fig. 3D, magenta columns). The average normalized power was not significantly affected (Fig. 3E). The respective phase relationships between the nerves were quantified via the coherence function, and these were not significantly affected by the transection (Fig. 3F).
A caudal transection has marginal effects on motor activity. A, Lumbar enlargement of the turtle spinal cord. Spinal segments and nerves are labeled. The transection is performed at the caudal end of segment D10 (dashed line), to allow access to the motoneuron in D10. B, C, Four simultaneous ENGs during scratch reflex measured before (cyan) and after (magenta) transection. One cycle from each recording is highlighted. D, Normalized single-trial ENG power of scratch behavior (frequency band, 0.3–0.9 Hz) in chronological order. The recordings in B and C are indicated by ★ and †, respectively. E, Average power across trials in D. Significant difference is indicated (★★★). F, Coherence between three nerve pairs (I, hip flexor vs extensor; II, dD8 vs hip flexor; III, hip flexor vs knee extensor) at the dominant frequency. Dots indicate single trials, and arrows indicate trial average.
Synaptic fluctuations
We compared the synaptic fluctuating potentials in neurons in the intact and transected spinal cord (Fig. 4) and separated the motor behavior into on and off cycles (Fig. 4A). Since the on cycles often have high spike density, we estimated the fluctuation via the SD of the spike-triggered distribution of Vm before the peak of the action potential (Fig. 4B). Fluctuations in the off cycles were determined as the SD of Vm over time (Fig. 4C). The distributions of synaptic fluctuations did not show a significant difference between the intact spinal cord and the transected spinal cord neither during the on-cycle periods (Fig. 4D) nor during the off-cycle periods (Fig. 4E).
Synaptic fluctuations in motoneurons are unaffected by transection. A, Membrane potential from a sample neuron during a scratch reflex (Iinj = 0). Spike times are indicated by cyan vertical lines. Selected off cycles are indicated by magenta bars. B, Synaptic fluctuations in on cycles are assessed as the trial variance of Vm, 18 ms before the triggered spikes (Berg et al., 2008), for intact (cyan) and transected (gray). C, Synaptic fluctuations in off cycles are assessed as SDs in 200 ms windows of Vm across time and cycles. D, E, Population distribution of synaptic fluctuations in intact (n = 10 cells) and transected (n = 8 cells) for the on cycle (cyan, D) and off cycle (magenta, E). Intact and transected populations are statistically indistinguishable in on cycles and in off cycles, using a Kolmogorov-Smirnov two-sample test (α = 5%).
Excitation and inhibition to motoneurons
We analyzed intracellular recordings from neurons in the transected (Fig. 2) and the intact (Figs. 5, 6) spinal cord to estimate E and I conductances during motor patterns. The synaptic inputs were estimated using both the ohmic and the ACF methods on data for which these methods were suitable (Berg and Ditlevsen, 2013). Motoneurons were identified according to their spiking activity with respect to the nerve activity (Fig. 5A) and their anatomical location (Fig. 5E). Sample estimates of Ginh and Gexc using the ACF method on single sweep data (Fig. 5B,C) had a significant correlation (Fig. 5D). The majority of cells in the population showed a similar positive correlation between Ginh and Gexc, although two cells showed a mild negative correlation (Fig. 6A). A sample experiment and estimates from one of these two cells is shown (Fig. 6B–D). Excitation is weak and in opposite phase with inhibition. The spike rate is low, and Gtot is strongest in the off cycles, which was atypical for the population.
Sample motoneuron in D10 from an intact spinal cord receiving balanced synaptic input. A, Hip flexor nerve recording during tactile induced pocket scratching. B, Vm during constant current injection (Iinj = 0). Action potentials clipped at −40 mV. C, The total conductance (gray), inhibition (magenta), and excitation (cyan) estimated with the ACF method using a 250 ms window. Shaded regions indicate the location for which the windows contained spikes. D, Inhibition and excitation estimated with the ACF method without spikes for three trials. The gray dashed line connects the x–y averages in six bins on the abscissa. Error bars indicate their respective SD. E, Transverse histological section to confirm the biocytin-filled neuron (red) used in A–D belongs to the group of motoneurons in the ventral horn (ChAT-positive neurons shown in green). Scale bar, 200 μm.
Balanced E/I in the majority of cells, with the exception of two cells. A, The correlation between excitation and inhibition for a population of neurons (n = 15), estimated via the ACF method (black) and ohmic method (gray). Most cells (n = 13) have a positive correlation between E and I, with the exception of two cells. A sample cell from Figure 5 is indicated by a star. B–D, Sample neuron receiving opposite timing of E and I. B, Flexor nerve recording. C, Membrane potential (Iinj = 0) of sample neuron 14 (†). Action potentials clipped at −30 mV. D, The total conductance (gray), inhibition (magenta), and excitation (cyan) estimated with the ACF method using a 200 ms window. Shaded regions indicate the omitted locations in which the windows contained spikes.
Structural division in a balanced network model
Since the experiments demonstrated mainly balanced E/I and resemblance between intact and transected spinal cord, it is relevant to consider the impact of structural divisions in a balanced network. We therefore constructed a simple model of the premotor network consisting of two populations of E and I neurons (Fig. 7A), with fixed random connectivity (van Vreeswijk and Sompolinsky, 1996; Brunel, 2000). What happens to the dynamics when the network is reduced to a fraction, α, such that the motoneurons receive inputs from fewer premotor neurons (Fig. 7B)? First, we noticed that the average firing rate of inhibitory neurons (Fig. 7C, red) was lower than the excitatory neurons (Fig. 7C, blue) and that their spiking was irregular (Fig. 7C, left). The irregular spiking was expected since it is a hallmark of balanced networks in the asynchronous state (Vogels et al., 2005). The lower firing rate of the inhibitory population was attributable to the choice of connection strengths and the balancing condition that recurrent inhibition has to dominate recurrent excitation to maintain stability, which lead to explicit expressions of rates (Eqs. 9, 10).
The model network remains balanced after structural division. A, Schematic of the architecture of the local network. Input is provided by an external population of excitatory neurons (CPG), which connects to both the excitatory (blue boxes) and inhibitory (red boxes) populations. Subsets of these populations project to each other resulting in recurrent inhibition and excitation. Other subsets of E/I project to motoneurons (MN; right gray box). B, Detailed schematic of a miniature premotor network consisting of inhibitory (red; Ninh = 100) and excitatory (blue; Nexc = 100) neurons. A subset of the population receives excitatory inputs from the external source (CPG; left gray cells), and another subset projects to local motoneurons (bottom right gray cells). A transection line is indicated with a dashed line. C, Rastergram of 10 sample neurons of both E and I populations before and after transection (50%). D, The mean spike rates of E (blue) and I (red) for the intact, 10 and 50% transected network as a function of external drive. E, The E (blue) and I (red) spike rate input to a sample motoneuron normalized by the total number of output projections (Nout,I and Nout,E) as a function of external input to the premotor network for intact, 10% and 50% transected network. F, The mean rates of E versus I neurons for intact and transected premotor network, and their slopes have marginal change after transection (inset).
A structural division will decrease the number of neurons, which should have little or no effect on the firing rates (Eqs. 9, 10). Nevertheless, there was a slight increase in the mean firing rate for both E and I neurons after transection (Fig. 7C). This was not caused by the reduction of the population directly, but rather by the elimination of some of the recurrent connections. The average number of presynaptic inputs to a neuron was diminished by α, such that K = αK′, where K′ is the original input to an average neuron. This will increase the mean firing rates after the division since they are proportional to 1/
Despite increasing firing rates for the remaining premotor neurons, the synaptic inputs to motoneurons will decrease because of the decreasing number of presynaptic cells after the division. Consequently, the higher average firing rate of the active neurons was counterbalanced by fewer active connections to motoneurons, and therefore we saw only a small effect on the motoneuron activity (Fig. 7E). Interestingly, the structural division of the network did not alter the balanced state. After transection of the spinal cord, the activity of E neurons still remained in proportion to I neurons (Fig. 7F). These observations are in qualitative accordance with the experimental observations: there was very little change in the ENG amplitude, i.e., the mean motoneuron input after transection (Fig. 3), the inputs to motoneurons were still balanced E/I (Figs. 5, 6), and the activity was still irregular (Fig. 4).
Discussion
In the present study, we investigated the sensitivity of the functional network integrity to a structural division at the D9/D10 segment in the turtle spinal cord. The study can be summarized into the following observations. First, there was no significant change in the ENG amplitudes after the transection (Fig. 3B–E). This suggests that the motoneurons had the same levels of activity and therefore received similar amounts of E/I input before and after the structural division. Second, the relative timing of synergists/antagonists seemed to be unaffected by the transection (Fig. 3F), which indicates the disconnected network could not have served a vital purpose for the network integrity. Third, the synaptic fluctuations were qualitatively similar before and after transection in both on and off cycles (Fig. 4B–E). This suggests that there was no substantial change in the amounts of synaptic inhibition and/or excitation in both phases. Fourth, waves of excitation and inhibition were qualitatively similar in intact versus transected spinal cord. The synaptic motoneuron inputs were balanced in most cells (n = 13 of 15), although two cells had stronger off-cycle inhibition (Figs. 5, 6). The balanced activity of the premotor network and its resilience to structural division is in accordance with the proposed model (Fig. 7), which also remained balanced after the transection (Fig. 7F). These results lend support to a premotor network with recurrent E/I circuitry, which has a remarkable ability to preserve functional integrity.
Resolving discrepancies between turtle and rat
Our data are at odds with the observation by Endo and Kiehn (2008) that locomotion in rats is mainly controlled by reciprocal inhibition with little or no excitation, which was confirmed in mice (Talpalar et al., 2011). There are a number of different possibilities for the discrepancy between the results. The premotor networks controlling scratching and locomotion could have different specialization as proposed for cats (Frigon and Gossard, 2010), or there could be a species difference. Alternatively, differences between rodents and turtles could be explained by differences in the precision of the movements that relate to developmental stage: neonatal for rodents and adult for turtles. For instance, in movements where the control of force is critical, such as the precision grip tasks in monkeys, there is coactivation between E/I premotoneurons (Takei and Seki, 2013; Wu and Perlmutter, 2013). Furthermore, inhibition has an important role in adjusting force during hindlimb control in adult cats (Johnson et al., 2012), suggesting that the degree of balance E/I depends on the particular task and force. Neonatal/embryonic rodents performing fictive locomotion are unable to carry themselves and therefore unable to adjust force appropriately. This motor inability is linked to immaturity of the spinal motor network and therefore could explain the discrepancy between the results obtained from adult turtles and neonatal rats. The importance of the developmental stage in forming a balanced network is seen in other parts of the nervous system. For instance, the auditory cortex is known to have balanced E/I in the adult rat (Wehr and Zador, 2003), but in early development, the network organization is qualitatively different. The balance only becomes significant between P15 and P25 and is fully developed at P30 (Dorrn et al., 2010). Since the rats used by Endo and Kiehn (2008) were newborn (P0–P04), they are substantially younger than the point at which they develop balanced E/I in the cortex.
Even if an E/I balance is present in the immature spinal circuit, resolving the synaptic conductances is difficult in these experiments, because (1) the estimation of synaptic conductance is distorted by the applied neurochemicals (NMDA/5-HT/dopamine) that chronically affect the receptors and (2) gap-junction couplings are present at this developmental stage. Gap junctions are important in circuitry formation and biochemical signaling, but since their abundance decreases with development (Chang and Balice-Gordon, 2000), they are probably less important in adult locomotion. Nevertheless, they have a tremendous impact on the electrical dynamics associated with neonatal fictive locomotion. For instance, rhythmic motor patterns can be communicated purely through electrical coupling without action potentials (Tresch and Kiehn, 2000). This electrical communication also distorts the estimation of synaptic conductance, which relies entirely on chemical neurotransmission (Berg and Ditlevsen, 2013).
Other work in turtle
The results from the transection experiments presented here (Fig. 3) are in accord with the early transection experiments by Mortin and Stein (1989). They demonstrated that the posterior 40–80% of the hindlimb enlargement is not essential for both rostral and pocket scratch generation. Currie and Gonsalves (1999) performed similar transection experiments, indicating that at least some of D10 is required to have on/off phase separation in the unilateral rhythm. Since our ENG recordings have clear on/off separation (Fig. 3C), the transection must have been more caudal than those performed by Currie and Gonsalves (1999). Nevertheless, the essential signals for motor pattern activation are carried by descending fibers (Berkowitz and Stein, 1994).
In their seminal work, Robertson and Stein (1988) performed a comprehensive assessment of the synaptic input during the off cycles in the key synergist and antagonist motoneurons. Though their work confirmed the presence of off-cycle inhibition, the authors did not consider on-cycle inhibition, perhaps because of their awareness of the absence of Renshaw inhibition (Ruigrok et al., 1984). Nevertheless, using the ohmic method on their trials with different current injections [Robertson and Stein (1988), their Fig. 3], we are able to corroborate the existence of an inhibitory on-cycle component (data not shown). The authors also manipulated the chloride reversal potential, which broadened the on cycle. This is likely attributable to a general broader inhibition and sharper excitation (Figs. 2C, 5C). Based on this meta-analysis, we suggest that there is a qualitative accord between past and present data.
Purpose of balanced input to motoneurons
We investigated inputs to hip flexor motoneurons, whereas the hip extensor motoneuron inputs are yet unresolved. A possible discrepancy between inputs to flexor and extensor motoneurons could be related to the particular motor behavior, which in this case is pocket scratching (Earhart and Stein, 2000). In this type of scratching, the knee is moved forward to rub against the location in the receptive field where the tactile stimulus occurred. Therefore, in this part of the cycle, referred to as the rub phase (Mortin et al., 1985), the hip flexors have to exert a force of high amplitude and precision, whereas the hip extensors are active in the postrub phase, which serves only to retract the limb without a well defined target. The movement in the postrub phase has high variability, whereas the rub phase has low variability (Mortin et al., 1985). This need for finer motor control makes it necessary for the premotor network to have accurate control of the hip flexor motoneurons, including their gain. Suitable gain modulation of neurons can be accomplished by balancing synaptic E/I (Silver, 2010), thus offering an explanation for the balanced input, and a prediction that hip extensors receive less balanced E/I, since their control is less critical. Similarly, balanced E/I has been observed among spinal premotor interneurons during precision grip tasks by monkey, where fine motor control is required (Takei and Seki, 2013; Wu and Perlmutter, 2013).
Premotor network architecture
The participating interneurons in spinal premotor networks are large in numbers and scattered across multiple segments (Stepien et al., 2010). The number of synaptic contacts onto a motoneuron can be up to 60,000 (Ulfhake and Cullheim, 1988), and motoneurons often exhibit irregular activity suggesting intense synaptic bombardment (Calvin and Stevens, 1967; Berg et al., 2008; Jahn et al., 2011; Dideriksen et al., 2012). There are eight interneurons for every motoneuron in the turtle spinal segment D9, and the number is similarly high in the lumbar enlargement in chickens (6:1) and mice (15:1) (Walløe et al., 2011). The relative high number of interneurons implies a more complicated architecture than the pure feedforward half-center model. If all neurons were feedforward connected, the signal propagation would either be slow, because of the vast number of intermediate neurons, or highly redundant through parallel connections. Alternatively, the large number of interneurons could constitute recurrent connectivity, which would provide higher processing capabilities. A balanced recurrent network has the potential of exhibiting a rich repertoire of dynamical states (Shu et al., 2003) and could therefore offer a novel direction in investigating the flexibility and multifunctionalism of motor networks (Berkowitz et al., 2010; Berkowitz and Hao, 2011). Nevertheless, recurrent excitatory connections will tend to make networks unstable (Vogels et al., 2005), and a dominant recurrent inhibition is therefore required (Brunel, 2000). Reciprocal connections will not stabilize the activity as they only give off-cycle inhibition.
Based on these observations, we propose a minimalistic model (Fig. 8). The model is composed of the generic reciprocal inhibition combined with a locally balanced premotor network (Stein, 2010), both receiving descending connections from a rhythm-generating network.
Sample data and proposed model for motor circuitry. A, B, Tactile induced motor behavior consisting of coordinated rhythmic muscle contractions (A) and depolarizations in the membrane potential of the corresponding motoneuron (B). C, The rhythmic behavior is caused by waves of inhibition (magenta) and excitation (cyan) to the motoneuron. In the on cycles, excitation and inhibition are concurrently peaking, hence the motoneuron is receiving balanced input from the premotor network (arrow), whereas in the off cycles, there is relatively more inhibition than excitation, which is likely caused by the activation of reciprocal inhibition from the antagonist (arrow and shaded region). Data are adapted from Berg and Ditlevsen (2013) with permission. D, The proposed model: descending inputs from a rhythm-generating network provides drive to both a premotor network, where recurrent excitation (blue population) is balanced with recurrent inhibition (red population), as well as generic reciprocal inhibitory interneurons. The connectivity and number of E/I neurons in the network is unknown. MN, Motoneuron; Ia, reciprocal inhibitory premotor interneuron.
Footnotes
This work was supported by the Novo Nordisk Foundation (R.W.B.) and The Danish Council for Independent Research Medical Sciences (R.W.B. and P.C.P.). We thank John Hertz and Jørn Hounsgaard for useful discussions.
- Correspondence should be addressed to Rune W. Berg, Faculty of Health and Medical Sciences, Department of Neuroscience and Pharmacology, University of Copenhagen, Blegdamsvej 3, DK-2200 KBH N, Denmark. runeb{at}sund.ku.dk
