Abstract
Renshaw cells (RCs) receive excitatory inputs from motoneurons to which then they inhibit. The gain of this spinal recurrent inhibitory circuit is modulated by inhibitory synapses on RCs. Inhibitory synapses on RCs mature postnatally, developing unusually large postsynaptic gephyrin clusters that colocalize glycine and GABAA receptors. We hypothesized that these features potentiate inhibitory currents in RCs. Thus, we analyzed glycinergic and GABAergic “inhibitory” miniature postsynaptic currents (mPSCs) in neonatal [postnatal day 1 (P1) to P5] and mature (P9-P15) RCs and compared them to other ventral interneurons (non-RCs). Recorded neurons were Neurobiotin filled and identified as RCs or non-RCs using post hoc immunohistochemical criteria. Glycinergic, GABAergic, and mixed glycine/GABA mPSCs matured differently in RCs and non-RCs. In RCs, glycinergic and GABAA mPSC peak amplitudes increased 230 and 45%, respectively, from P1-P5 to P9-P15, whereas in non-RCs, glycinergic peak amplitudes changed little and GABAA amplitudes decreased. GABAA mPSCs were slower in RCs (P1-P5, τ = 58 ms; P9-P15, τ = 43 ms) compared with non-RCs (P1-P5, τ = 27 ms; P9-P15, τ = 14 ms). Thus, fast glycinergic currents dominated “mixed” mPSC peak amplitudes in mature RCs, and GABAA currents dominated their long decays. In non-RCs, GABAergic and mixed events had shorter durations, and their frequencies decreased with development. Functional maturation of inhibitory synapses on RCs correlates well with increased glycine receptor recruitment to large gephyrin patches, colocalization with α3/α5-containing GABAA receptors, and maintenance of GABA/glycine corelease. As a result, charge transfer in GABA, glycine, or mixed mPSCs was larger in mature RCs than in non-RCs, suggesting RCs receive potent inhibitory synapses.
Introduction
Inhibitory synapses develop a large variety of functional, structural, and molecular properties to counteract synaptic excitation and shape-firing patterns. These characteristics adjust inhibitory current amplitudes and timing to the integration properties of neurons and synaptic networks (Browne et al., 2001; Pouille and Scanziani, 2001; Semyanov et al., 2004). They also homeostatically regulate the balance between synaptic excitation and inhibition (Turrigiano and Nelson, 2004). The number and position of inhibitory synapses along dendritic arbors, the amount of GABA or glycine released or coreleased, and the characteristics of postsynaptic neurotransmitter actions influence postsynaptic currents, and all are developmentally regulated. Postsynaptic actions are in turn dependent on the number and properties of GABAA and glycine receptors. Postsynaptic receptor selection results in synapses that are mainly glycinergic, others GABAergic, and some mixed (Jonas et al., 1998; Kotak et al., 1998; O'Brien and Berger, 1999; Smith et al., 2000; Gao et al., 2001; Keller et al., 2001; Nabekura et al., 2004). These neurotransmitter/receptor choices greatly influence inhibition strength and timing (Russier et al., 2002). GABAergic transmission predominates over glycinergic transmission during early development, when both neurotransmitters depolarize neurons. Thereafter, the balance switches to glycine-dominated transmission in many spinal and brainstem synapses. The present consensus is that GABAergic currents are more efficient triggering long depolarizations and calcium entry related to developmental mechanisms, including synaptic formation and stabilization. Thereafter, inhibitory synapses shift their function to synaptic integration, and a bias toward “faster” glycinergic mechanisms seems advantageous for many ventral spinal cord synapses. GABAA receptor recruitment influences GABAergic current amplitudes (Nusser et al., 1997, 1998), and their time course is strongly related to subunit composition (Browne et al., 2001; Vicini et al., 2001). Glycinergic currents show an interesting property in that postsynaptic current amplitude is correlated with the amount of postsynaptic gephyrin (Lim et al., 1999; Van Zundert et al., 2004), a glycine receptor-anchoring protein (Kirsch and Betz, 1993; Feng et al., 1998; Meier et al., 2000, 2001).
One type of spinal interneuron, the Renshaw cell (RC), stands out because it displays a high density of proximal inhibitory synapses with large postsynaptic densities intensely immunoreactive to gephyrin and containing colocalized GABAA and glycine receptors (Alvarez et al., 1997; Geiman et al., 2002). These receptors are postsynaptic to boutons expressing markers, suggesting significant GABA/glycine corelease. We hypothesized that these structural and neurochemical properties increase synaptic inhibition by simultaneously enhancing peak amplitudes and time courses of postsynaptic currents (Geiman et al., 2002). These characteristics are likely adapted to RC function because, unlike other interneurons, they receive motor axon excitatory inputs of large amplitude and long duration that evoke high-frequency burst firing (Renshaw, 1946; Eccles et al., 1954). Synaptic structure on RCs matures postnatally in conjunction with ventral horn motor output and locomotion (Geiman et al., 2000). To investigate functional features characteristic of inhibitory synapses on RCs, we analyzed the amplitude and time courses of spontaneous miniature inhibitory currents and their GABAergic and glycinergic components during postnatal development and compared them to other ventral interneurons of similar ages.
Materials and Methods
Tissue preparation and immunohistochemical processing. The distribution, density, and immunocytochemical properties of RCs were analyzed in spinal cord sections obtained from early [postnatal day 1 (P1), P5] and late (P10, P15) postnatal rat spinal cords. This information was then used to target our recordings to a specific area of the spinal cord containing high densities of RCs and to classify recorded cells as “Renshaw” (RC) or “non-Renshaw” (non-RC) using post hoc immunohistochemistry.
For histological processing, a total of 12 pups (three rats in each age group) were used. Rats were anesthetized with an overdose of sodium pentobarbital (>90 mg/kg, i.p.) and perfused transcardially with cold vascular rinse (0.01 m phosphate buffer with 137 mm NaCl, 3.4 mm KCl, and 6 mm NaHCO3, pH 7.3), followed by fixative (4% paraformaldehyde in 0.1 m phosphate buffer). Spinal cords were removed, postfixed for 1 h at room temperature, and cryoprotected at 4°C in 0.1 m phosphate buffer with 15% sucrose. Twenty-micrometer-thick transverse sections were obtained on a cryostat from L4-L5 spinal cord segments and collected on slides. Some sections were dual immunolabeled for neuronal nuclear protein (NeuN) [mouse monoclonal antibody (mAb); Chemicon, Temecula, CA], a generalized neuronal marker, and for calbindin [rabbit polyclonal antibody (pAb); Swant, Bellinzona, Switzerland], an RC marker in the ventral horn of the developing spinal cord (Geiman et al., 2000). In other sections, these markers were combined with immunostaining for the vesicular acetylcholine transporter (VAChT) (goat pAb; PharMingen, San Diego, CA). Different dual or triple immunolabelings were used to map the distribution of RCs to estimate their density in ventral horn regions targeted for recording (see Fig. 1 A, B) and to analyze the competence of criteria for distinguishing RCs from neighboring non-RC interneurons in neonatal rat spinal cord based on the presence of abundant “cholinergic contacts” (see Fig. 1C-F).
The sections were washed in 0.01 m phosphate buffer containing 0.9% saline and 0.1% Triton X-100 (PBS/TX) and blocked with normal horse serum (NHS) (1:10 in PBS/TX). For dual (calbindin/NeuN or calbindin/VAChT) or triple (calbindin/NeuN/VAChT) color immunofluorescence, the sections were incubated overnight at 4°C with different combinations of primary antibodies against calbindin D28K (rabbit pAb; dilution, 1:2000), NeuN (mouse mAb; 1:1000), or VAChT (goat pAb; 1:1000). Immunofluorescence was visualized with donkey anti-rabbit, anti-mouse, or anti-goat antibodies (1:50 dilutions) coupled with α methylcoumarin (AMCA), cyanine 3 (Cy3), Cy5, or FITC (Jackson ImmunoResearch, West Grove, PA). All primary and secondary antibodies were diluted in PBS/TX. Finally, after washing the excess of antibodies with PBS, the sections were coverslipped with Vectashield (Vector Laboratories, Burlingame, CA).
Confocal microscopy and image analysis. Dual Cy3 and FITC immunofluorescence and triple Cy3, FITC, and AMCA immunofluorescence were visualized with a Fluoview FX laser-scanning confocal microscope (Olympus, Tokyo, Japan), except that AMCA fluorescence was recorded with normal epifluorescence and a digital camera (Spot2; Diagnostic Instruments, Sterling Heights, MI). Triple color immunolabelings including FITC, Cy3, and Cy5 fluorochromes were imaged in a TCS confocal system (Leica, Nussloch, Germany). The proportions of calbindin-immunoreactive (IR) profiles to all small-sized NeuN-IR profiles in ventral LVII and LIX were estimated in low-magnification confocal images (10×) using ImagePro-Plus software (version 3.0.01; Media Cybernetics, Silver Spring, MD). The border between LIX and LVII was delineated by the distribution of motoneuron cell bodies (large NeuN-IR profiles). To accommodate for the progressive spinal cord growth and to minimize sampling variations, counting in ventral LVII was performed using a rectangular window, the linear dimensions of which were increased in steps of 25 μm for each successive older age group (200-275 μm width and 100-175 μm height from P1 to P15) (see Fig. 1 A, B). The box was placed in ventral LVII with the lateral border aligned to the medial margin of the LIX lateral motoneuron pool and the ventral border with the ventral funiculus gray-white matter boundary. The box in LIX was placed as a lateral extension in LIX of boxes drawn in LVII at each age. Sections from 12 animals (3 per age group), with at least three sections per animal (i.e., six ventral horns), were included in the analysis. To reduce errors because of assumptions about similarities in cell body size, shape, or orientation of ventral interneurons, only nucleolus-containing IR profiles were counted. Errors caused by differential splitting of cells at the surfaces were avoided by counting only nucleoli at least 3 μm within the tissue sections. These precautions were taken to assure unbiased sampling of profiles from RCs versus other neuronal subtypes. All counts were performed in the lumbar 4 and 5 segments.
To reconstruct VAChT coverage on calbindin-IR or NeuN-IR cells, optical section stacks (step size, 0.5 μm) were obtained throughout the cell body and proximal dendrites of the neurons. Dual immunofluorescence was captured in two-channel mode with a 60× oil immersion objective (numerical aperture, 1.4) digitally zoomed 1.5×. Cross talk between channels was avoided as in previous reports (Geiman et al., 2002). Images used for figure composition were filtered (high-Gauss filter, Image Pro-Plus 4.0; Media Cybernetics) and adjusted for contrast, brightness, and dynamic resolution for best quality presentation without changing or altering the information content in the images.
Slice preparation. We recorded ventral horn interneurons from spinal cord slices obtained from 1- to 5- or 9- to 15-d-old postnatal rat pups. All animals were anesthetized (50 mg/kg pentobarbital) and decapitated, and their spinal cords were removed quickly. Dissection was in ice-cold (<4°C) sucrose artificial CSF (S-aCSF) bubbled with 95% O2 and 5% CO2. The S-aCSF composition was as follows (in mm): 26 NaHCO3, 10 glucose, 3 KCl, 1.25 NaH2PO4, 2 MgCl2, 2.4 CaCl2, and 218 sucrose. Transverse slices from the L4-L5 segments (300-400 μm thick) were cut using an OTS-4000 tissue slicer (Electron Microscopy Sciences, Fort Washington, PA). Slices were transferred to normal oxygenated aCSF (in which sucrose was removed and 130 mm NaCl and 2 mm CaCl2 were added) and incubated for 1 h at 36°C and then allowed to stabilize at room temperature (22-25°C) for at least 30 min before electrophysiological recordings. Slices were transferred into the recording chamber and perfused continuously (at a rate of 4 ml/min) with oxygenated normal aCSF.
Whole-cell recordings. Whole-cell voltage-clamp recordings were obtained from somata of visually identified ventral horn interneurons (<25 μm diameter) using an Olympus BX 50WI microscope equipped with infrared-differential interference contrast (DIC), a 40× water immersion objective, and a CCD camera (CCD-Iris; Sony, Tokyo, Japan). Preferential recording areas were ventral LVII and LIX regions containing the highest densities of immunocytochemically identified RCs (see Fig. 2 A, B). Patch electrodes (2-5 MΩ resistance) contained the following (in mm): 120 CsCl, 4 NaCl, 4 MgCl2, 1 Cl2Ca, 10 HEPES, 0.2 EGTA, 3 Mg-ATP, and 0.3 GTP-Tris. In all of the experiments, 0.4% Neurobiotin (Vector Laboratories) was added to the internal solution. Only recordings with access resistance between 5 and 20 MΩ were considered acceptable for analysis. The access resistance was checked throughout the experiments, and recording was abandoned if it changed >15%. Junction potentials (usually <5 mV) were corrected after inserting the pipette into the bath. Cells were voltage clamped at -75 mV. In our recording conditions, the reversal potential for chloride currents was close to 0 mV; thus, GABAergic and glycinergic currents were readily detected as inward currents well separated from baseline noise. Synaptic currents were recorded and low-pass bessel filtered at 5 kHz with an Axopatch 200B amplifier (Axon Instruments, Union City, CA). Data were digitized at 10 kHz and acquired using Axograph (version 4.6; Axon Instruments). For each cell, we obtained two to three segments of 5 min continuous recording of spontaneous activity under drug combinations that pharmacologically isolated the synaptic currents of interest.
To isolate miniature spontaneous synaptic currents of GABAergic and/or glycinergic origin [miniature postsynaptic currents (mPSCs)], recordings were performed in the presence of tetrodotoxin (TTX) (1 μm; Alomone Labs, Jerusalem, Israel), 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX) (10 μm; Tocris Cookson, Ballwin, MO), and d-tubocurarine chloride (10 or 30 μm; Sigma, St. Louis, MO) applied by bath perfusion. Glycinergic or GABAergic mPSCs were subsequently isolated by adding either bicuculline methiodide (10 μm; Sigma) or strychnine hydrochloride (0.25 μm; Sigma) to the aCSF solution.
Data analysis. Off-line data analysis and curve fitting was performed with both pClamp 9.0 and MiniAnalysis (Synaptosoft, Decatur, GA). Events were detected by setting the threshold value for detection at three times the level of the root mean square noise (∼3 pA). Therefore, the detection threshold was 8-10 pA. We routinely scrolled through detected events and visually rejected any superimposed or spuriously detected events. Peak mPSC amplitudes were measured at the absolute maximum of the currents, taking into account the noise at baseline and around the peak. Time to decay was calculated as the time from peak to 33% decay of peak amplitude. Rise times were measured between 10 and 90% of the peak. The charge transferred (in picocoulombs) was measured as the integrated area under averaged mPSCs.
Curve fitting to mPSC decay (from peak to end) was performed on averaged traces with simplex algorithm least squares exponential fitting provided by MiniAnalysis software and using single- or double-exponential equations of the form: y = A0 + A e(-x/τ) and y = A0 + Afe(-x/τf) + Ase(-x/τs), respectively. Decay phases of individual mPSCs recorded under TTX, CNQX, and d-tubocurarine chloride were fitted either by one or two exponential functions that usually resulted in events falling into one of three classes: “fast decay,” “slow decay,” or mixed events (see Results). The averaged glycinergic and GABAergic mPSCs in individual neurons were always best fitted by a biexponential decay. In contrast, individual and averaged mPSC decays of pharmacologically isolated glycinergic or GABAergic events were best fitted by single monoexponential curves. Although this procedure yielded good correlation fits (usually r > 0.95) for averaged events and many individual mPSCs, it failed in a proportion of individual events in which curve fits clearly departed from trace points (r < 0.7). For this reason, classification of events falling in fast-decaying, slow-decaying, and mixed categories in recordings containing mixed GABAA and glycinergic currents was aided by visual inspection of individual events in which the program failed to provide good curve fitting. These events showed fluctuations in their decay phases usually because of superimposition of other events or random baseline oscillations. The proportion of events added by visual identification of inflections in their decays were, on average (of all cells analyzed), as follows: P1-P5 RCs, 14.8% of all events; P9-P15 RCs, 4.0%; P1-P5 non-RCs, 13.4%; P9-P15 non-RCs, 1.8%. Not surprisingly, more events were added in younger cells displaying the longer decays (Tables 1, 2, 3). In a few young neurons with high frequencies of spontaneous mPSCs they represented 20-25% of all events. Using this method, we avoided underestimation of the total numbers of mPSCs for the purpose of calculating frequencies. Decay time constants of visually selected fast- and slow-only events were not significantly different from the decay constants of, respectively, glycinergic or GABAergic currents pharmacologically isolated. Mixed events were clearly identified by the presence of an inflection separating slow and fast kinetic components. We expect that any biases or errors introduced by the experimenter when estimating the proportion of different classes of events should be minimal and, in any case, smaller than the errors incurred by relying only in curve fitting of individual events by the analysis program. Two-way ANOVA and post hoc Duncan's multiple range tests with p < 0.05 were used for statistical comparisons among different ages and classes of neurons. Data were expressed as mean ± SE.
Immunohistochemical processing of recorded slices. Neurons were filled with Neurobiotin during recording. Spinal cord slices containing Neurobiotin-filled neurons were fixed in 4% paraformaldehyde for 30 min and then stored at 4°C in PBS. Immunohistochemistry was performed in the thick slices used for recording. To facilitate penetration of antibodies, tissue sections were dehydrated and rehydrated through a graded series of alcohols and xylene. Slices were then washed in PBS, blocked with NHS (1:10 in PBS/0.2% TX), and immersed free-floating in primary antisera against VAChT (diluted 1:1000 in PBS/TX; the concentration of Triton X-100 was adjusted in these thick slices to 0.3% to favor penetration). Incubations were performed overnight at room temperature and under agitation. Thereafter, slices were washed and incubated for 2 h in streptavidin-Cy3- and FITC-coupled donkey anti-goat antibodies diluted in PBS/TX (1:50 dilutions; Jackson ImmunoResearch). Finally, slices were mounted on glass slides and coverslipped with Vectashield (Vector Laboratories).
Neuronal localization (Cy3-Neurobiotin labeling) was first visualized under epi-fluorescence with an Olympus BX60 microscope. Low-magnification (10×) and high-magnification (40×) images were captured with a digital color camera (Spot2 camera; Diagnostic Instruments) to plot the location of the neuron somata and record major morphological properties of dendrites and axon trajectories. VAChT immunoreactivity in low-magnification images pointed to the location of motoneurons and LIX. The cell body locations of all neurons included in the study (n = 41 at P1-P5 and n = 51 at P9-P14) were plotted in relation to the section outline, the ventral horn gray matter, and the location of LIX. The average representation of these locations is shown in Figure 2 B. Finally, the density of VAChT-IR contacts was analyzed with confocal microscopy. Labeled neurons (Cy3) and VAChT immunofluorescence (FITC) were scanned in two-channel mode in an Olympus Fluoview FX confocal microscope. Cell morphology and VAChT-IR bouton contacts were reconstructed from stacks of optical sections (1 μm step size) obtained at 60×2 magnification. The recorded neurons were classified into RC and non-RC groups according to the density of VAChT-IR terminals contacting the somatodendritic surfaces, as well as the location of the cell body, extension of the dendrites, and the presence of axon collaterals in the ipsilateral motoneuron pool (frequently containing synaptic-like varicosities). The characteristics of GABAergic and glycinergic synaptic currents were compared throughout postnatal development between both groups of neurons.
Results
Targeting and identification of recorded neurons as RCs
RCs are characterized by monosynaptic cholinergic inputs from motor axon collaterals, strong calbindin expression, large gephyrin clusters on their proximal somatodendritic membranes, and preferential location of lower lumbar segments in ventral LVII and ventral LIX (Eccles et al., 1954; Jankowska and Lindstrom, 1971; Fyffe, 1990; Ardvisson et al., 1992; Alvarez et al., 1997, 1999; Carr et al., 1998). Two previous studies in spinal cord slices used stimulation of motor axons as search stimulus to identify and record RCs (Oleskevich et al., 1999; Dourado and Sargent, 2002). In both studies, ventral root cholinergic monosynaptic inputs were identified in <10% of recorded neurons. This low yield could have resulted from targeting ventral horn regions with low densities of RCs and perhaps also from motor axon damage during slicing. To overcome these problems, we determined anatomically the optimal sampling sites in the neonatal spinal cord and tested several properties for RC post hoc histochemical identification.
Intense calbindin immunoreactivity is a distinctive feature of RCs (Ardvisson et al., 1992; Carr et al., 1998; Geiman et al., 2000). Ventral horn neurons with intense calbindin immunoreactivity are mostly located in ventral LIX and LVII (Fig. 1A,B), in a region corresponding to the adult RC area (Thomas and Wilson, 1965; Jankowska and Lindstrom, 1971; Fyffe, 1990). The proportion of calbindin-IR neurons to all NeuN-IR small-sized neurons increased during postnatal development from 30.1 ± 2.6% (P1-P5; mean ± SE) to 44.3 ± 3.6% (P10-P15) in ventral LVII and from 31.6 ± 3.8% to 57.3 ± 5.4% in ventral LIX (Fig. 1F) (see Materials and Methods for a description of the boxed regions of interest that were analyzed). At the same time, expansion of the neuropil reduced neuronal density from P1 to P15, particularly in ventral LVII (Fig. 1, compare C, D). In conclusion, sampling in the low cell density region shown in Figure 2A should result in a 30-60% chance of recording from an RC. This region appears darker with DIC optics because of the increased density of myelinated axons, likely motor axons, crossing this area.
This distinctive regional distribution is important for targeting recordings to RCs, but it is insufficient to identify this neuronal phenotype. Therefore, we investigated other features that could permit RC identification. Adult calbindin-IR RCs display high densities of cholinergic contacts on their dendrites (Alvarez et al., 1999). Similarly, neonatal calbindin-IR RCs received higher densities of VAChT-IR contacts than other NeuN-IR neurons (Fig. 1C-E). Two other features of adult RCs, large gephyrin clusters and calbindin immunoreactivity, were inadequate for RC characterization in the present study. Gephyrin-IR clusters on RCs undergo postnatal maturation (Geiman et al., 2000), and in neonates they are not different from gephyrin-IR clusters on other neurons (data not shown). On the other hand, calbindin immunoreactivity in recorded neurons was weaker and more difficult to detect than in non-recorded cells. Possibly, calbindin dilution into recording pipettes significantly lowers calbindin immunoreactivity in neurons recorded in whole-cell mode.
Neurons selected for recording were visualized with DIC optics in the target region. Thereafter, post hoc immunohistochemical confirmation of a high density of VAChT-IR contacts on Neurobiotin-filled neurons was a necessary condition for RC classification (Fig. 2C,D); 27.5% (60 of 218) of intracellularly recorded and stained cells were classified as RCs in our target population. Only neurons with good quality recordings (>25 min of stable recordings at constant low series resistance) and intracellular staining and in sections with successful VAChT immunolabeling were analyzed. VAChT immunolabeling of C terminals on neighboring motoneurons (Alvarez et al., 1999) was used as internal control of immunolabeling quality. The locations of neurons analyzed are shown in Figure 2B. We analyzed 20 neurons identified as RCs at each age (P1-P5 and P9-P15). All RCs were preferentially located within the ventralmost 200 μm of LVII and LIX. RC axons were directed toward ipsilateral motor pools, and sometimes varicosities were observed in LIX. The dendrites of RCs formed small arbors, mostly restricted to the RC area. Non-RC interneurons were more disperse in the ventral horn, occasionally distant from the motoneuron pools, and their morphological features were very variable. Axons recovered from non-RCs projected toward ipsilateral LIX motoneurons pools, crossed the midline into the contralateral spinal cord, or coursed toward the lateral funiculus. Some dendritic arborizations in this cell group were large and widespread, whereas others had restricted distributions. Thus, non-RCs likely comprise several neuronal subtypes. We analyzed 21 non-RCs at P1-P5 and 31 at P9-P15.
Differential maturation patterns of GABAergic/glycinergic mPSCs in ventral horn interneurons
We first compared all “inhibitory” mPSCs isolated with TTX, CNQX, and d-tubocarine in RCs (n = 7 at P1-P5 and n = 7 at P9-P15) and non-RCs (n = 9 at P1-P5 and n = 8 at P9-P15). These mPSC events are mixtures of GABAergic and glycinergic events, hereafter called mPSCsGABA/Gly, and were blocked by 10 μm bicuculline and 0.2-1 μm strychnine (Fig. 2 E, F).
In both cell types, some events decayed rapidly (decay time, ≤15 ms), others slowly (decay time, >15 ms), and some were mixed (events with fast and slow components, frequently separated by inflections in their decays). RCs always showed a significant population of mPSCs with very slow decays. These were infrequent in non-RCs. This difference was evident at early stages (Fig. 2E,F) and increased with age (Fig. 3A,B). At late postnatal ages, few mPSCGABA/Gly events (3.76%) were longer than 20 ms in non-RCs (Fig. 3A,C,E); these were more common in RCs (16.6%) (Fig. 3B,D,F). Rise time versus decay time plots showed a compression of the decay time range in non-RCs during postnatal development (Fig. 3C). On the contrary, mPSCs-GABA/Gly in RCs of ages P1-P5 and P9-P15 showed similar wide ranges of decay times (Fig. 3D). Thus, decay time cumulative distributions shifted to the left in non-RCs, indicating faster decay kinetics in older animals (Fig. 3E), whereas the cumulative curve from P9-P15 RCs was slightly displaced to the left for faster events and to the right for slower events (Fig. 3F). It is unlikely that electrotonic filtering contributed to differences in decay between both cell types. First, mPSCGABA/Gly mean rise times in P9-P10 RCs and non-RCs were similar (1.5 ± 0.2 vs 1.3 ± 0.1 ms, respectively; p = 0.42; two-way ANOVA; Duncan's test). Second, correlations between rise and decay times were never found (r < 0.2) (Fig. 3C,D). Finally, the dendritic arbors of non-RCs (with shorter decay times) were usually larger than in RCs, and RCs have less distal inhibitory synapses than non-RCs (Alvarez et al., 1997).
Peak amplitudes of mPSCsGABA/Gly were also differentially regulated during postnatal development in each neuronal type. The mean peak mPSCGABA/Gly amplitude increased by 121.3% [100 × value(P9-P15) - value(P1-P5)/value(P1-P5)] in the RC group (p < 0.05; two-way ANOVA; Duncan's test). In contrast, the average peak amplitude increase in non-RCs was moderate and did not reach statistical significance (61.2%; p = 0.07). Cumulative probability plots of peak mPSCGABA/Gly amplitudes shifted to the right with development in both RCs and non-RCs. The magnitude of change was much greater for RCs (Fig. 3G,H), and peak amplitude distribution histograms in P9-P15 RCs showed strong skews toward high-amplitude events.
In addition, we observed developmental changes in frequency, but these are best reported below, following descriptions of individual and mixed GABA/glycine unitary events.
Glycinergic mPSC development
To analyze the specific neurotransmitter components of these inhibitory synaptic currents, we first isolated pharmacologically glycinergic mPSCs (called mPSCsGly) by blocking glutamatergic (CNQX), nicotinic (d-tubocurarine), and GABAA (bicuculline) receptors. At early postnatal times, mPSCsGly were similar in both cell types (Fig. 4A,B, top traces). Mean peak amplitude, 10-90% rise time, decay time, and time constant (τgly) were not significantly different between both cell types (p ≥ 0.05; two-way ANOVA; Duncan's test) (Table 1). During development, mPSCGly decays became faster in both cell types (Fig. 4A-F; C,D, insets; Table 1). Rise time versus decay time plots (Fig. 4C,D) and distribution histograms (Fig. 4E,F) revealed a shortening in the range of decays, the mean values of which for RCs and non-RCs were, respectively, 1.69 and 1.82 times faster to those in P1-P5 neurons (Table 1). Common kinetic changes in pharmacologically isolated mPSCsGly probably reflect similar maturation of postsynaptic glycine receptor functional properties in both cell types. In contrast, mPSCGly amplitudes matured differently in RCs and non-RCs (Fig. 4G,H; Table 1). P9-P15 RCs exhibited a large increase (230.4%) in mean mPSCGly peak amplitude, relative to P1-P5 RCs. In contrast, the smaller increase (31%) in mean peak amplitude in non-RCs did not reach statistical significance (p = 0.2; two-way ANOVA; Duncan's test). Consequently, non-RC mPSCGly amplitude distribution histograms of different ages essentially overlapped, and only a very small shift was noted in their cumulative probability functions (Fig. 4G). In contrast, amplitude distributions from P9-P15 RC mPSCsGly were mostly skewed toward higher values, and their probability functions significantly shifted to the right (Fig. 4H). Thus, whereas mPSCsGly with amplitudes >100 pA were infrequent in RCs at early ages (4.5% of all events) and in non-RCs at early and late postnatal ages (11.0 and 14.8%, respectively), events >100 pA comprised nearly 50% of all mPSCsGly in P9-P15 RCs.
The frequency of pharmacologically isolated glycinergic events increased more than twofold from early to late postnatal ages in both RCs and non-RCs (see Fig. 6A3, Table 1). Event frequencies were always higher in non-RCs of all ages. This observation parallels morphological findings indicating the presence of 5-15 small independent active zones per synapse on non-RCs compared with fewer numbers of larger active zones in inhibitory synapses contacting adult RCs (Alvarez et al., 1997).
GABAergic mPSC development
In the ventral spinal cord, spontaneous miniature GABAergic synaptic events (mPSCsGABA) are generally downregulated in frequency and amplitude, and their duration is shortened during postnatal development (Gao et al., 2001). However, significant numbers of bicuculline-sensitive mPSCs with remarkably slow kinetics persist in mature RCs (Fig. 3). To analyze their postnatal development, mP-SCsGABA were pharmacologically isolated in the presence of 0.25 μm strychnine added to the CNQX, d-tubocurarine, and TTX mixture (Fig. 5A,B). Pharmacologically isolated GABAergic mPSCs recorded from non-RCs (n = 8) and RCs (n = 6) at P1-P5 were similar in 10-90% rise time and peak amplitude (p > 0.05; two-way ANOVA; Duncan's test) (Fig. 5G,H; Table 2). However, mean decay time and τGABA of mPSCsGABA in neonatal RCs was 2.1 times longer than in non-RCs (p < 0.05; two-way ANOVA; Duncan's test) (Fig. 5A-F, Table 2).
Differences in pharmacologically isolated mPSCsGABA between both cell types increased with development (Fig. 5A-H). The average decay times and τGABA of mPSCsGABA were shortened in P9-P15 non-RCs, becoming nearly threefold faster than those measured in RCs of similar age (Table 2). Rise time versus decay time plots, distribution histograms, or cumulative probability curves were compressed or shifted toward shorter values in P9-P15 non-RCs. In contrast, their ranges essentially overlapped in RCs from both ages (Fig. 5C-F). Nevertheless, the averaged τGABA suggested an acceleration of mPSCGABA decay kinetics in RCs with development (p < 0.05; two-way ANOVA; Duncan's test) (Table 2). Mean peak amplitudes significantly increased by 45.1% in RCs (p < 0.05; two-way ANOVA; Duncan's test), while showing a small and nonsignificant decline in non-RCs (p = 0.27) (Table 2).
Together, these results suggest significant differences in the functional development of inhibitory synaptic currents in RCs compared with non-RC interneurons (Fig. 6).
Postnatal development of glycine-GABA mixed mPSCs
The characteristics of inhibitory currents in RCs and non-RCs are also influenced by the amount of cotransmission between GABA and glycine. We analyzed cotransmission in both cell types by studying the characteristics of dual-component mPSCs from recordings obtained in the absence of bicuculline and strychnine, as those shown in Figures 2 and 3. According to the characteristics of pharmacologically isolated glycinergic and GABAergic mPSCs, we characterized three types of events based on their decay kinetics and developmental maturation (Fig. 7A): fast (≤15 ms at P1-P5 and ≤10 ms at P9-P15), slow (>15 ms at P1-P5 and >10 ms at P9-P15), and mixed (with an inflection separating fast and slow components). Decay parameters were estimated from decay phases best fitted to monoexponential or biexponential curves (r > 0.9 in ∼80-95% fits in different cell groups). However, a population of mPSCs with long decays sometimes displayed superimposed events or other noise fluctuations that prevented good fits (r < 0.7). To avoid underestimation of the number of events in each category, clear mPSCs with random fluctuations in their decays were classified as fast, slow, or mixed by visual inspection and were added to the populations of mPSCs classified by curve fitting. Visually classified events are a small proportion of all of the mPSCs analyzed in each cell group (P1-P5 RCs, 14.8% of all events; P9-P15 RCs, 4.0%; P1-P5 non-RCs, 13.4%; P9-P15, non-RCs, 1.8%). Decay time constants of averaged single fast and slow mPSCs selected in this way were very similar to pharmacologically isolated glycinergic and GABAergic currents in each neuronal subgroup (p ≥ 0.1; two-way ANOVA). Moreover, the decay constants of fast (τf) and slow (τS) components of biexponentially decaying events coincided with those obtained from single “fast” and “slow” mPSCs (Fig. 7A, Table 3).
The relative percentage of fast events increased significantly from early to late postnatal ages (Fig. 7B,C) in non-RCs (from 41.9 to 77.5%) and RCs (from 26.7 to 56.3%), in parallel to the developmental increase in the frequency of glycinergic mPSCs in both cell types (Fig. 6A3). Higher percentage values of fast mPSCs in non-RCs relative to RCs result from greater numbers of mPSCsGly and a lower frequency of mPSCsGABA in non-RCs. The fraction of slow-decaying events was significantly reduced with age in non-RCs, whereas it remained unchanged in RCs. Finally, the percentage of mixed mPSCs decreased in both neuronal types, although to a larger extent in non-RCs (Fig. 7B,C). The reduction in the number of mixed events in RCs probably derives from the increased frequency of pure fast glycinergic events. Finally, the downregulation of GABAergic currents in non-RCs contributes to reduce further the incidence of dual mPSCs in these neurons compared with RCs.
To further analyze the contribution of fast (likely glycinergic) and slow (likely GABAergic) components to mixed events in RCs and non-RCs, we estimated the relative contribution to absolute peak amplitudes in mixed mPSCs of the fast and slow components (Af, As) (Fig. 7D, Table 3). Af and As were calculated from the parameters of decay curve fittings using the biexponential formula y = A0 + Afe(-x/τf) + Ase(-x/τs). As indicated by the Af/As ratios in Figure 7D, the relative contribution of the fast component was similar in RCs and non-RCs at early postnatal ages, accounting for ∼60% of the total peak amplitude (Table 3). The relative contribution of fast versus slow components was similar in non-RCs of different ages, whereas Af/As ratios were significantly increased in RCs at later postnatal ages (Fig. 7D) (p < 0.05; two-way ANOVA; Duncan's test). In conclusion, GABAergic mPSCs and GABAergic components in mixed mPSCs were maintained during development to a larger extent in RCs compared with non-RCs. Thus, GABAergic currents are responsible for the retention in RCs of mixed events with slow decays during development, but their contribution to peak amplitudes is reduced in older ages because the large increase in the fast glycinergic currents.
Charge transfer of GABA/glycine mPSCs in mature ventral horn interneurons
Although GABAergic components might contribute to peak amplitude less than glycinergic components, their long time courses suggest strong contributions to the total inhibitory synaptic current. We calculated the total charge transferred by glycinergic and GABAergic components in single and dual mPSCs in RCs and non-RCs recorded from P9-P15 animals. Absolute charge transfer associated with glycinergic, GABAergic, or mixed currents recorded from RCs was always larger than for non-RCs [Fig. 8, filled circles (RC), open circles (non-RCs); Tables 1, 2, 3). Both peak amplitudes and time courses contributed to increase charge transfer, but the relative contributions of each parameter was different for GABA or glycine mPSCs. Correlations between charge and peak amplitude or decay constant of glycinergic, GABAergic, or dual mPSCs suggest that absolute charge transfer is better related with peak amplitude, but not with decay time, in fast glycinergic mPSCs and with both decay time and peak amplitude in slow GABAergic mPSCs (Fig. 8A1-B2). As a result, in mixed mPSCs, both peak amplitude (dominated by the fast glycinergic current) and time course (dominated by the second decay from the slower GABAergic currents) are both well correlated with peak amplitude (Fig. 8C1-C3).
Discussion
Inhibitory synapses on mature RCs are characterized by larger peak amplitudes, slower decays, and higher percentages of GABAergic currents compared with other ventral interneurons. In combination, these characteristics increase inhibitory charge transfer in RCs compared with non-RCs. Given the high density of inhibitory synapses on RCs (Harrington et al., 1994), we conclude that during postnatal development, inhibitory synaptic inputs on RCs become greatly potentiated compared with other ventral interneurons. These functional features can be explained by the development of specific structural and molecular properties of inhibitory synapses on RCs and appear well adapted to RC function.
Maturation of postsynaptic inhibitory current peak amplitudes
The most striking feature of inhibitory synapses on RCs was the large increase in glycinergic mPSC peak amplitudes with development. The magnitude of this postnatal increase is larger than in other neurons studied previously (Gao et al., 1998, 2001; Kotak et al., 1998; Singer et al., 1998; Smith et al., 2000; Nabekura et al., 2004; Van Zundert et al., 2004) and correlates with the increase, specifically on RCs, of gephyrin clustering, a protein involved in trapping inhibitory receptors at postsynaptic densities. From P2 to P15, average gephyrin cluster size on RCs increased from 0.19 to 0.41 μm2, and the proportion of very large clusters (>0.5 μm2) increased from 5 to 35% (Geiman et al., 2000). These large gephyrin clusters set RCs apart from other neurons (Alvarez et al., 1997).
Larger gephyrin clusters present an opportunity for trapping and anchoring more postsynaptic receptors and increase postsynaptic currents. Indeed, average gephyrin cluster sizes were correlated with average glycinergic mPSC peak amplitudes (Lim et al., 1999; Van Zundert et al., 2004). Modification of postsynaptic receptor properties could also contribute to increased peak amplitude. However, glycine receptor properties mature similarly in RC and non-RCs. Similar studies on cranial motoneurons also concluded that receptor accumulation is the major contributor to postnatal increases in average glycinergic mPSC peak amplitude (Singer and Berger, 1999). Thus, the best explanation for the large increase in glycinergic current amplitudes in RCs is higher recruitment of glycine receptors to large postsynaptic gephyrin clusters.
Gephyrin has been also related to postsynaptic GABAA receptor clustering (Craig et al., 1996; Essrich et al., 1998; Sassoe-Pognetto and Fritschy, 2000); however, GABAergic mPSC peak amplitudes were enhanced to a smaller degree in developing RCs. GABAergic mPSC amplitudes can also increase through accumulation of postsynaptic GABAA receptors (Nusser et al., 1997, 1998), but a recent report found only weak correlation with gephyrin clustering (Levi et al., 2004). Gephyrin interactions with glycine and GABAA receptors suggest a tighter relationship between gephyrin clustering and glycine receptors than with GABAA receptors (Meyer et al., 1995; Meier et al., 2000, 2001; Dahan et al., 2003; Hanus et al., 2004). Gephyrin disruption causes complete declustering of glycine receptors (Kirsch et al., 1993; Feng et al., 1998) and loss of GABAA receptors (Essrich et al., 1998; Kneussel et al., 1999), but significant gephyrin-independent GABAA receptor clusters remain (Fischer et al., 2000; Kneussel et al., 2001; Levi et al., 2004). Moreover, during early synapse formation, postsynaptic GABAA receptors cluster postsynaptically before gephyrin (Dumoulin et al., 2000; Danglot et al., 2003; Levi et al., 2004). Thus, although gephyrin contributes to stabilize GABAA receptors at postsynaptic densities (Danglot et al., 2003), it is not always necessary for GABAA receptor postsynaptic recruitment and clustering. In addition, GABAA receptors are best stabilized by a gephyrin splice isoform containing the C5 cassette (Meier and Grantyn, 2004); therefore, postsynaptic GABAA receptor numbers and currents might be best related to the relative abundance of specific gephyrin isoforms than to overall gephyrin content.
Maturation of postsynaptic inhibitory current time course
The long time courses of GABAA postsynaptic currents are major contributors to total charge transfer underlying inhibitory currents in RCs and are well explained by the expression of α3/α5-containing GABAA receptors (Geiman et al., 2002). Recombinant α3-containing GABAA receptors are the least sensitive to GABA (Sigel et al., 1990; Verdoorn, 1994), resulting in slower opening, deactivation, inactivation, and postsynaptic current lengthening (Gingrich et al., 1995; Serafini et al., 1998). Decay time constant and total charge transfer in RC GABAergic mPSCs was strikingly similar to those reported in reticular thalamic (Rt) neurons (Browne et al., 2001), another neuron that expresses α3/α5 GABAA receptors. Interestingly, Rt neurons and RCs have common functional features: both provide recurrent negative feedback and respond to excitatory synaptic inputs with high-frequency burst firing. The properties of α3/α5-containing GABAA receptors might be well adapted to the common integrative properties exhibited by these cells.
Slow-decaying GABAA currents were good predictors of whether or not recorded neurons would display anatomical features of RCs. In addition, GABAergic currents decayed 25% faster in mature RCs compared with immature RCs. A parallel speeding of GABAA currents was found for the faster decays of non-RCs. These changes could suggest developmentally regulated alterations in GABAA receptor subunit structure affecting channel kinetics. Alternatively, they could be explained by alterations in receptor localization (Chery and De Koninck, 1999), β subunit phosphorylation (Hinkle and MacDonald, 2003), or GABA reuptake maturation (Draguhn and Heinemann, 1996) (but see Chery and De Koninck, 1999; Nusser and Mody, 2002; Overstreet and Westbrook, 2003). Our immunocytochemical data suggest that GABAA receptors in RCs are highly clustered at synapses in adults and neonates (Geiman et al., 2002; T. Culberston and F. J. Alvarez, unpublished observations), indicating that major developmental changes in receptor localization in RCs are unlikely. Our data also suggests that GABAA receptors in RCs contain β3 subunits (Geiman et al., 2002). Protein kinase A-dependent phosphorylation of β3 subunits slows down GABAA receptor deactivation and prolongs postsynaptic currents (Hinkle and MacDonald, 2003). In addition, neurosteroid prolongation of GABAergic currents is potent on β2/3-containing GABAA receptors (Wingrove et al., 1994). In the spinal cord, neurosteroid modulation weakens with age (Keller et al., 2004), and this decline is perhaps partly responsible for the acceleration of GABAergic currents.
A reduction in glycinergic current decay time was also observed in both RCs and non-RCs during development. A switch from neonatal α2- to adult α1-containing glycine receptors occurs during the first two postnatal weeks in the rat spinal cord (Akagi and Miledi, 1988; Becker et al., 1988; Malosio et al., 1991) and correlates with a decrease in glycine receptor mean open time and PSC decays (Takahashi et al., 1992; Singer et al., 1998). Decay parameters measured at P1-P5 were faster than those reported in neurons with predominant expression of α2 glycine receptors and similar to neurons at intermediate developmental stages expressing both α1 and α2 glycine receptors (Singer et al., 1998; Van Zundert et al., 2004). Previously, we detected significant amounts of α1 glycine receptors in RCs just after birth (Geiman et al., 2000). Therefore, glycinergic decay acceleration in RCs is best explained by progressive downregulation of neonatal α2 subunits.
Maturation of cotransmission
GABA and glycine corelease is generally demonstrated by spontaneous mPSCs displaying biphasic decays with partial time constants similar to those of GABAergic and glycinergic mPSCs. Because mPSCs are interpreted as the postsynaptic action of neurotransmitter released from single vesicles, the presence of mixed mPSCs suggests GABA and glycine corelease from individual vesicles (Jonas et al., 1998), a possibility also supported by the use of a common vesicular transporter (Sagne et al., 1997; Chaudry et al., 1998) (but see Katsurabayashi et al., 2004). Biphasic mPSCs indicative of corelease have been shown at synapses in the spinal cord (Jonas et al., 1998; Chery and De Koninck, 1999; Gao et al., 2001; Keller et al., 2001), cerebellum (Dumuolin et al., 2001), and brainstem (O'Brien and Berger, 1999; Russier et al., 2002; Nabekura et al., 2004). A postnatal decrease in the level of cotransmission was reported for many neurons, and it is usually associated with downregulation of GABAergic currents and a switch to mainly glycine-mediated inhibitory synaptic mechanisms (Kotak et al., 1998; Gao et al., 2001; Keller et al., 2001; Nabekura et al., 2004). Most ventral interneurons show a similar developmental pattern; however, RCs were different because of the remarkable maintenance of GABAergic and mixed currents. Correspondingly, morphological data indicate that 60-75% of inhibitory terminals on adult RCs coexpress presynaptic and postsynaptic markers suggesting mixed neurotransmission, a percentage higher than in other ventral horn synapses (Geiman et al., 2002). The mechanisms that regulate the amount of cotransmission are unknown (van den Pol, 2004), but the choice has important consequences on the timing and strength of synaptic inhibition. High levels of cotransmission over RCs enhance inhibitory current strength by permitting large peak amplitudes while maintaining very slow time courses. Russier et al. (2002) showed that inhibitory synapses decreased motoneuron firing more efficiently when GABAergic and glycinergic transmission occurred simultaneously compared with experimental situations in which one or other component was blocked. Cotransmission also allowed inhibition with complementary time courses, transient by glycine and sustained by GABA.
Inhibitory modulation of RC function
The functional properties of inhibitory synaptic currents in RCs increase synaptic inhibition in amplitude and duration and develop in parallel to the spinal circuitry underlying locomotion. In the adult spinal cord, RC-mediated recurrent inhibition modulates the dynamic behavior of motoneuron firing by acting as a gain regulator of their discharge and recruitment (Windhorst, 1996; Uchiyama et al., 2003; Hultborn et al., 2004). The fine organization of recurrent inhibition on different motoneuron types and pools allows RCs to influence motoneuron firing variability, force gain, and motor discharge synchronization in homonymous and synergistic motoneuron pools and muscles. RCs also facilitate activity in antagonists through regulation of reciprocal Ia inhibitory interneurons (Hultborn et al., 1979). Thus, despite the relative simplicity of the recurrent inhibitory loop, RCs exert a variety of complex functions during motor activity. Not surprisingly recurrent inhibition is finely modulated during different motor behaviors. For example, the level of recurrent inhibition changes during voluntary movements according to the strength of muscle contraction and motor task (Hultborn and Pierrot-Deseilligny, 1979; Katz and Pierrot-Deseilligny, 1999; Iles et al., 2000). This modulation implies the existence of effective mechanisms to decouple RC firing from motor axon inputs. However, RCs receive long-lasting cholinergic EPSPs from motor axons and display unusual low thresholds for action potential generation. These peculiarities result in the characteristic tendency of RCs to discharge bursts of spikes in response to incoming motor axon volleys, even when only one action potential travels down a single motor axon (for review, see Windhorst, 1990). Modulation of RC activity therefore requires complementary well-matched inhibition. Inhibitory synapses on RCs are strategically located more proximal to the cell body than cholinergic synapses (Alvarez et al., 1999) and their synaptic structure (Alvarez et al., 1997) and molecular/neurochemical organization (Geiman et al., 2002) seem to correlate with the large amplitude and relatively long time courses of inhibitory currents described here.
In conclusion, the functional, structural, and molecular characteristics of inhibitory synapses on RCs appear well adapted to counteract robust and long-lasting motor axon-driven excitatory inputs in RCs and allow effective modulation of RC activity and recurrent inhibition of motoneurons.
Footnotes
This work was supported by National Science Foundation Grant 9984441 to F.J.A. D.G.-F. was supported by a postdoctoral fellowship from the Ministerio de Educación, Cultura y Deporte (Spain). We thank Drs. Javier Stern and Bruce Walsmley for technical advice and helpful comments and suggestions on a previous version of this manuscript.
Correspondence should be addressed to Dr. Francisco J. Alvarez, Department of Anatomy and Physiology, Wright State University, 3640 Colonel Glenn Highway, Dayton, OH 45435. E-mail:francisco.alvarez{at}wright.edu.
Copyright © 2005 Society for Neuroscience 0270-6474/05/252010-14$15.00/0