Original Research ARTICLE
The heterogeneity in GABAA receptor-mediated IPSC kinetics reflects heterogeneity of subunit composition among inhibitory and excitatory interneurons in spinal lamina II
- 1Unité de Neurosciences Cellulaires et Moléculaire, Institut Universitaire en Santé Mentale de Québec, Québec, QC, Canada
- 2Department of Biological Applications and Technology, University of Ioannina, Ioannina, Greece
- 3Laboratory of Genetic Neuropharmacology, McLean Hospital, and Department of Psychiatry, Harvard Medical School, Belmont, MA, USA
- 4Department of Psychiatry and Neuroscience, Université Laval, Québec, QC, Canada
GABAergic inhibition displays rich functional diversity throughout the CNS, which arises from variations in the nature of inputs, subunit composition, subcellular localization of receptors and synapse geometry, or reuptake mechanisms. In the spinal dorsal horn (SDH), GABAA and glycine receptors play a major role in the control of excitability and accuracy of nociceptive processing. Identifying which components shape the properties of the inhibitory synapses in different cell types is necessary to understand how nociceptive information is integrated. To address this, we used transgenic mice where inhibitory interneurons express GAD65-EGFP. We found that GABAA, but not glycine receptor-mediated evoked IPSCs displayed slower kinetics in EGFP+ vs. EGFP− interneurons. GABAA miniature IPSC decay kinetics showed a large variability in both populations, however the distribution of decays differed between EGFP+ and EGFP− interneurons. The range of mIPSC decay kinetics observed was replicated in experiments using rapid application of GABA on outside-out patches taken from SDH neurons in slices. Furthermore, GABAA decay kinetics were not affected by uptake blockers and were not different in mice lacking δ or α5 subunits, indicating that intrinsic channel properties likely underlie the heterogeneity. To identify whether other α subunits shape the various kinetic properties observed we took advantage of knock-in mice carrying point mutations in either the α1, α2, or α3 subunits rendering Ro 15-4513 a selective agonist at the benzodiazepine modulatory site. We found that α1 and α2 subunit underlie the fast decaying component of IPSCs while the slow component is determined by the α3 subunit. The differential distribution of GABAA subunits at inhibitory synapses thus sculpts the heterogeneity of the SDH inhibitory circuitry. This diversity of inhibitory elements can be harnessed to selectively modulate different components of the spinal nociceptive circuitry for therapeutic interventions.
The spinal dorsal horn (SDH) serves both as a relay and processing station for somatic sensory information entering from the peripheral sensors in transit to the higher brain centers. Primary afferents carrying somatic sensory input terminate in the spinal cord in a laminar distribution pattern. Aδ - and peptidergic C-fibers carry nociceptive, thermal, itch and innocuous tactile information and terminate in lamina I and outer II while non-peptidergic C-fibers terminate predominantly in inner lamina II (Ribeiro-da-Silva and De Koninck, 2008). Incoming nociceptive information converges onto projection neurons in lamina I following processing through a local network of interneurons in lamina II that control the final output. The inhibitory neurotransmitters GABA and glycine and their receptors in the dorsal horn play multiple roles in the control of information flow, the discrimination of sensory modalities and thus in securing the accuracy in the transmission of sensory information. Indeed, intrathecal GABAA receptor agonists increase the nociceptive threshold in rats (Hammond and Drower, 1984), while administration of GABAA or glycine receptor antagonists produce pain hypersensitivity and allodynia (Yaksh, 1989; Sivilotti and Woolf, 1994; Sorkin and Puig, 1996). Recent work suggests that reduced inhibition in the dorsal horn underlies neuropathic and inflammatory pain (Coull et al., 2003; Harvey et al., 2004). Such a disinhibition lowers the control of interconnected sensory networks allowing unorthodox information flow from normally innocuous inputs to nociceptive pathways (Baba et al., 2003; Torsney and MacDermott, 2006; Keller et al., 2007).
GABAergic transmission displays a high level of heterogeneity in the brain. Several types of interneurons that release GABA form synapses with different properties on various compartments of the target cells. An important potential substrate of heterogeneity at GABAA synapses is their molecular composition: 19 subunits are available to assemble the pentameric GABAA receptors (Olsen and Sieghart, 2009). Several subunits, including α1, α2, α3, and α5, are expressed in the SDH (Bohlhalter et al., 1996). Additionally, several distinct morphological and histochemical classes of interneurons that can realease GABA have been identified (Laing et al., 1994). However, despite the important functions fulfilled by inhibitory interneurons, little is known on the functional organization of inhibitory networks, their connections, mode of operation or properties. Such knowledge is imperative in light of recent work unveiling new opportunities for the development of selective treatments for chronic pain by targeting GABAA receptors (Knabl et al., 2008; Zeilhofer et al., 2009). However, as we asserted previously, inhibitory interneurons in the SDH can play a dual role in moderating spinal excitability, restraining it by feed forward inhibition of the polysynaptic excitatory relay circuit, but also promote it by reciprocal silencing of other inhibitory interneurons (Labrakakis et al., 2009). Hence, the complex functionality of inhibitory networks has to be taken into account (Cossart et al., 2005). In this context, we explored the properties of GABAA–mediated synaptic events in subclasses of interneurons in the SDH, using mice that express EGFP under the control of the inhibitory neuron promoter glutamate decarboxylase (GAD65). We found that the properties of GABAA IPSCs are different between inhibitory and presumably excitatory interneurons. In addition we show that the heterogeneity of kinetic properties we observed results from difference in subunit composition of the receptors. These findings have bearings on strategies for selective modulation of subcomponents of the inhibitory circuitry in the SDH.
Materials and Methods
All experiments were performed in accordance with regulations of the Canadian Council on Animal Care. Experimental procedures were approved by the Comité de protection des animaux de l'Université Laval.
Several genetically modified mouse lines were used in this study. Heterozygous GAD65-EGFP transgenic mice express the enhanced green fluorescent protein (EGFP) under the control of the GAD65 promoter (Lopez-Bendito et al., 2004; Labrakakis et al., 2009). Gabra5−/− and Gabrd−/− mice are deficient in the GABAA receptor α5 and δ subunits respectively (Mihalek et al., 1999; Collinson et al., 2002). Homozygous α1(H101R), α2(H101R), and α3(H126R) mice are carrying a mutation substituting a histidine for an arginine at the benzodiazepine binding site of the GABAA receptor α1, α2 and α3 subunits, respectively (Rudolph et al., 1999; Low et al., 2000). Wild type and genetically modified mice were all of the C57BL/6J, C57BL/6J x129/SvEv (α5−/−) or C57BL/6J × 129/SvJ (δ −/− and δ +/+) background. GAD65-EGFP mouse phenotype was determined by the presence of GFP. The genotype of all other mouse lines was ascertained by PCR.
Preparation of Spinal Cord Slices
Adult mice (3–6 months) of both sexes, were anesthetized with ketamine/xylazine and perfused intracardially with ice-cold oxygenated (95% O2, 5% CO2) sucrose substituted ACSF containing (mM) 252 sucrose, 2.5 KCl, 1.5 CaCl2, 6 MgCl2, 10 glucose, 26 NaHCO3, 1.25 NaH2PO4 and 5 kynurenic acid as previously described (Chery et al., 2000; Labrakakis et al., 2009). Mice were decapitated, the spinal cord was removed by hydraulic extrusion and 250 μm thick parasagittal slices were cut from the lumbar portion. Slices were transferred in normal oxygenated ACSF (126 NaCl, 2.5 KCl, 2 CaCl2, 2 MgCl2, 10 glucose, 26 NaHCO3, 1.25 NaH2PO4) and incubated at 33°C for 1 h and then kept at room temperature until recording.
Slices were transferred in the recording chamber and continuously superfused at 2–3 ml/min with oxygenated ACSF at room temperature (23–26°C). Dorsal horn neurons were visualized with a Zeiss Axioplan2 microscope equipped with infrared “gradient-contrast” optics, epifluorescence and a × 40 water immersion objective. Patch pipettes (borosilicate glass; 6–8 MΩ) were filled with (in mM) 135 CsCl, 10 HEPES, 2 MgCl2, 0.5 EGTA. For some recordings, CsCl was substituted by 130 CsSO3CH3/5 CsCl. Whole cell patch clamp recordings were made using a Multiclamp 700B amplifier (Molecular Devices, Sunnyvale, CA). Access resistance was monitored periodically throughout the experiment. Recordings were not analyzed if access resistance was unstable or exceeded 30 MΩ. Data were low-pass filtered at 3 kHz, digitized at 10 kHz and acquired with the Strathclyde electrophysiology software (WinWCP and WinEDR courtesy of Dr. J. Dempster, University of Strathclyde, Glasgow, UK). All recordings in this study are from neurons with their somata located in laminae II.
Monosynaptic IPSCs were evoked focally by electrical stimulation (30–70 μ A, 250 μs) via patch pipette filled with ACSF and placed 50–100 μm from the soma of the recorded cell as described previously (Labrakakis et al., 2009). Single stimuli were delivered every 10 s. The mIPSCs were detected and analyzed using Mini Analysis (Synaptosoft, Decatur, GA) and locally designed software (YDK). Detection thresholds for mIPSCs were set at three times the RMS of noise. Decay time constants were fitted using automated least square algorithms. The necessity to introduce additional exponential components to the fits was first judged on the basis of visual inspection. When the merit of additional components was not obvious, further statistical analysis was applied as previously described (Chery and De Koninck, 1999). The weighted decay time constant (τw) was calculated from dual-exponential fits using the following equation: τw = (τ1A1 + τ2A2)/(A1 + A2) where τ1 and τ2 are the fast and slow decay time constants and A1 and A2 are the equivalent amplitude weighting factors. The cumulative probability plots of mIPSC decay τ's were fitted by mixtures of Gaussian distributions:
Where R1, …, Rn, are the ratios of the n normal distributions (such that ), xi, …, xn are the means and σ1, …, σn, the standard deviations (Cordero-Erausquin et al., 2009). The necessity to introduce additional Gaussian components to the fits was judged first on the basis of visual inspection of the fitted curves superimposed onto the data. When the merit of additional components was not obvious, an F-test was used to assess how the additional component improved the value of the reduced chi-square: Fx = Δχ2/χ2ν (where χ2ν = χ2/ν and the factor ν = N − n is the number of degrees of freedom left after fitting N data points to the n parameters; df1 = 3 and df2 = ν). The critical value for the merit of additional components was set at a low level (p < 0.0001) to favor parsimony of the fitted function (De Koninck and Mody, 1994; Chery and De Koninck, 1999). For the analysis of the effects of Ro 15-4513 (Sigma) on different populations of mIPSCs we categorized them as fast if their decay τw was <100 ms or slow if their decay τw was >100 ms. The 100 ms cut off was chosen on the basis of the Gaussian distributions fitted from the cumulative probability plots; 100 ms is the interface between the slower and the faster Gaussian components (Figure 4A right).
Rapid Agonist Application on Excised Outside Out Patches
Stable outside out membrane patches were excised by pulling the pipette away from a whole-cell patched neuron. Excised patches were placed at the interface of a double-bore glass flow pipe with control ACSF and 1 mM GABA-containing solutions. Rapid exchange was achieved by fast displacement using a piezoelectric positioning system (Physik Instrumente, Germany) as previously described (Bowie et al., 1998). Solution exchange speed was determined at the end of each experiment by measuring open tip currents resulting from the liquid junction potentials between control and 0.5x ACSF (rise and decay typically ranged between 400 and 500 μs). Data were discarded from patches in which the liquid junction currents exhibited slow rise times.
Simulation, Analysis, and Statistics
To simulate the effect of dendritic filtering on mIPSC decay kinetics, a simple ball and stick model was used in NEURON (Hines and Carnevale, 1997) software. Dendrite diameter was 1.5 μm, axial resistivity 300 Ω cm, membrane capacitance 1 μF/cm2. A passive leak conductance of 0.2 mS/cm2 was distributed throughout the cell. Voltage clamp series resistance was set to 30 MΩ. The synapses were modeled as a conductance change at different distances along the dendrite. For this a recorded mIPSC with fast decay kinetics (decay 27.3 ms, 10–90% rise: 1.2 ms) was used. This model is adequate to compare the relative effects of dendritic filtering on rise time vs. decay time but cannot accurately describe the effect of neuron morphology on dendritic filtering.
To model the effect of receptor binding affinity on the decay kinetics, simulated synaptic currents were generated with Channelab (Synaptosoft, Decatur, GA). The 5th-order Runga-Kutta numerical integrator was used for simulated macroscopic currents. For simulated mIPSCs, the Monte Carlo simulator was used. Gaussian noise (3 kHz) was added to the simulated mIPSCs. The latter were analyzed similarly to real mIPSCs, using locally developed software.
Normality of the data was tested with the Shapiro-Wilk test. Parametric data are expressed as mean ± s.e.m. and Student's independent t-test was used to analyze differences. One-Way ANOVA was used to analyze IPSC differences between wt and point mutated knock-in mice and post-hoc tests were obtained with Bonferroni or Tukey corrections to compare means. The non-parametric, distribution-free Kolmogorov–Smirnov test was used to compare cumulative probability distributions. Unless otherwise noted, n's represent number of neurons.
Differential Decay Kinetics of Evoked GABAA IPSCs between Lamina II Interneuron Subpopulations
To investigate heterogeneity of GABAergic transmission in the dorsal horn we used parasagittal slices from GAD65-EGFP transgenic mice. Neurons were visually identified under epifluorescence as EGFP+ and thus inhibitory or as EGFP− and mainly excitatory interneurons (Lopez-Bendito et al., 2004; Labrakakis et al., 2009). Focally evoked monosynaptic GABAA IPSCs were recorded in the presence of glutamate and glycine receptor blockers, CNQX (10 μM), APV (40 μM), and strychnine (0.5 μM). At a holding potential of −70 mV with high Cl− pipettes, eIPSCs were inward and showed complex deactivation kinetics consisting of two exponentials both in EGFP− and EGFP+ neurons (Figure 1A). eIPSC amplitudes were similar in the two populations (46.6 ± 5.5 pA in EGFP− and 38.6 ± 5.2 pA EGFP+ neurons, p > 0.05, n = 21 and 28 respectively). Comparison of the decay kinetics, revealed differences in the weighted decay time constant (τw; 69 ± 6.3 ms, n = 21 in EGFP− and 132.5 ± 11.1, n = 28 in EGFP+ neurons; p < 0.001). Decay kinetics in EGFP− or EGFP+ cells did not differ between male and female mice (65.8 ± 6.4 and 75.4 ± 12.3 for EGFP−, 150.7 ± 17.8 and 116.6 ± 12.9 for EGFP+, male and female respectively; ANOVA, p > 0.05). Additional experiments were performed using low Cl− pipettes and a holding potential of 0 mV. Under these conditions the reversal potential for Cl− and GABAA currents is more negative than the holding potential. As expected, focal stimulation produced outward eIPSCs that displayed double exponential decay kinetics. In agreement with the previous experiments eIPSCs decay kinetics were significantly slower in EGFP+ neurons (68.7 ± 7.4 vs. 190.3 ± 28.2 ms for EGFP− and EGFP+ respectively, n = 12 and 14, p < 0.05) while amplitudes were not significantly different (61.8 ± 13.8 pA, n = 12 for EGFP− and 36.9 ± 6.3 pA, n = 14 for EGFP+, p > 0.05; Figure 1B). Finally, comparison of the biexponential decay components separately showed that the main difference between the two interneuron subpopulations was in the slow decay component in both sets of experiments (Table 1).
Figure 1. GABAA receptor-mediated eIPSCs differ in their decay kinetics between lamina II dorsal horn neuron subpopulations. (A) IPSCs evoked by focal stimulation recorded from a EGFP negative (GFP−) and a EGFP positive (GFP+) neuron. Recordings were made at a −70 mV holding potential with CsCl containing pipettes. Histograms show the mean amplitude and weighted decay time constant (τw) in the two cell populations. The asterisk denotes significant difference (n = 21 for GFP− and n = 28 for GFP+). (B) Example traces of eIPSCs recorded from GFP− and GFP+ cells at a holding potential of 0 mV and using CsSO3CH3-containing patch pipettes. Histograms show the mean amplitude and decay τw. Asterisk denotes significant difference (n = 12 for GFP− and n = 14 for GFP+).
Table 1. Summary of decay kinetic values for GABAA eIPSCs recorded at holding potentials of −70 and 0 mV.
Decay Kinetics of Glycinergic eIPSCs were Similar Among Interneuron Subpopulations
In addition to GABA, fast inhibitory neurotransmission in the spinal cord is also carried by glycine. A large proportion of the GABAergic inhibitory interneurons in the lamina II of the dorsal horn also co-release glycine, while all the glycinergic neurons also contain GABA (Todd and Sullivan, 1990; Chery and De Koninck, 1999). We asked the question if glycine receptor mediated synaptic transmission also displays different decay properties between the EGFP+ and EGFP− neuron populations. Glycinergic eIPSCs were pharmacologically isolated in the presence of APV, CNQX and the GABAA receptor antagonist SR95531 (10 μM) and recorded at 0 mV. Evoked eIPSCs under these conditions displayed monoexponential decay kinetics (Figure 2). Comparison of amplitudes (128.4 ± 23.9 pA vs. 130 ± 15.34, p > 0.05) and decay time constants (τ = 12.35 ± 2.1 vs. 12.2 ± 1.5) showed that they were not different in EGFP− (n = 8) and EGFP+ (n = 10) neurons. This indicates that decay kinetic differences are restricted to the GABAA transmission and is not due to a difference in electrotonic properties of the cell populations.
Figure 2. Glycine receptor-mediated eIPSCs have similar decay kinetics in GFP+ and GFP− neurons. Example traces of eIPSCs were recorded from GFP− and GFP+ neurons at a holding of 0 mV. Neither amplitude nor decay time constant were significantly different in the two populations (n = 6 for GFP− and n = 6 for GFP+).
Differences in GABAA mIPSCs between Interneuron Subpopulations
Evoked stimulation results in the synchronous recruitment of multiple release sites which may lead to accumulation of the neurotransmitter in the synaptic cleft that may spillover out of the synapse to extrasynaptic receptors (Isaacson et al., 1993). Alteration of the transmitter time course in the cleft or activation of distant receptors by long range diffusion could be responsible for the slow decay component in the eIPSCs and for the observed differences. To minimize the effect that diffusion or massive accumulation of neurotransmitter might have in shaping the decay kinetics of GABAA IPSCs we looked at quantal release events in EGFP− and EGFP+ interneurons. Miniature IPSCs (mIPSCs) were recorded in the presence of TTX (1 μM) and strychnine (0.5 μM; Figures 3A,B) at a holding potential of −70 mV with high Cl− pipettes. Mean mIPSC frequencies were higher in EGFP− interneurons (0.15 ± 0.02 Hz, n = 6) compared to EGFP+ interneurons (0.04 ± 0.01 Hz, n = 6, p < 0.05). In addition mIPSCs in EGFP− cells were found to have larger amplitudes (−19.54 ± 1.8 pA, n = 6) than EGFP+ cells (−11.57 ± 1.2 pA, n = 6, p < 0.05; Figure 3C). The mIPSCs showed a large variability in their decay kinetics (Figure 3C) in both cell populations. Figure 3D shows the cumulative probability plot of the τw showing differential distributions in the two interneuron populations. To reveal if the large variability of decay kinetics could be the result of electrotonic filtering we compared the decay τw to the rise slope (10–90 rise time/amplitude) of the mIPSCs. We found no correlation between rise and decay times (Figure 3D), indicating that the decay kinetics variability was not a result of filtering of mIPSC from remote synaptic locations. In addition, to limit the influence of presumably strongly filtered mIPSCs, we confined our analysis to a subset of mIPSCs with faster rise times possibly arising closer to the soma. Cumulative probability plots of the decay τw (Figure 3E) confirm the differential distribution of decays in the two interneuron populations, while no correlation could be found between rise time and decay τw for this subset of our data.
Figure 3. GABAA mIPSCs show both fast- and slow-decay kinetics. (A) Example of a recording of mIPSCs from a GFP+ (top) and a GFP− (bottom) neuron displaying variable decay kinetics. (B) Different examples of single mIPSCs with variable decay kinetics: fast (left), mixed (middle) and slow decay (right) events. (C) Histograms of mIPSC frequency and amplitude showing significant differences (asterisk) between GFP− (n = 6) and GFP+ (n = 6) cell populations. The population distribution histogram of decay time constants (left, bin width: 20 ms) shows that GFP− and GFP+ display different proportions of fast and slow mIPSCs. (D) The cumulative probability plot (left) shows differential distribution of mIPSC decay τw between GFP− and GFP+ cell populations (left, p < 0.05 Kolmorgorov-Smirnov test). The plot consists of pooled mIPSC decay time constants from 6 GFP− and 6 GFP+ cells. The same number of consecutively occurring mIPSCs from each cell was used. The relationship between 10 and 90% rise slope and the decay τw in EGFP− and EGFP+ mIPSCs is plotted on the right. There was no correlation between mIPSC rise and decay neither in GFP− (r = 0.22) nor in GFP+ (r = 0.05) neurons. (E) Cumulative probability plot (left) and rise-slope/decay τw relationship plot (right) of the subset of data with rise times <4 ms. The distribution of decay τw between GFP− and GFP+ cells was different (p < 0.05), while there was no correlation between rise time and decay τw in GFP− (r = 0.11) and GFP+ (r = 0.14) neurons.
The above results indicate that electrotonic filtering is not the main factor shaping decay kinetics. We thus considered the possibility that the large decay variability was due to different mIPSC populations with variable properties. Indeed, the cumulative probability plots indicate multiple possible subpopulations of mIPSC decays. We were able to fit the data with a mixture of three Gaussian distributions (Figure 4A). The means and standard deviations (SD) of the three distributions were similar in both interneuron populations (Table 2). This indicates that the same subpopulations of mIPSCs make up the synaptic events of both EGFP− and EGFP+, however the relative contribution of each mIPSC subpopulation was different. To investigate the rise time and amplitude distributions within these mIPSCs subpopulations we divided the mIPSCs in three groups based on the means and SD derived from the Gaussian fits (mean ± 1 × SD). Rise time distribution for each group of mIPSCs was similar within each interneuron type and followed the distribution of the total population; however, mIPSCs in EGFP+ interneurons showed slower rise times than EGFP− mIPSCs (Figure 4B). This shows that rise time is cell type dependent and confirms that decay τw does not depend of the rise time. Amplitude distribution for the two faster decaying mIPSC groups was similar within each interneuron population with those in EGFP+ interneurons showing smaller size (Figure 4C). Interestingly, the amplitudes of the slower decaying mIPSCs in EGFP− interneurons show smaller size distributions akin to those in EGFP+ cells.
Figure 4. Analysis of mIPSC populations. (A) The cumulative probability plots for GFP− and GFP+ neurons (left and middle) were fitted with the sum of three normal distribution functions (red line). The resulting three distributions are plotted (right) separately (black lines for GFP−, red lines for GFP+. The green line is showing the 100 ms cut-off point used to separate slower from faster events). (B) Cumulative probability plots of the 10–90% rise times for GFP− (black) and GFP+ (red) mIPSCs. Thick lines show the total mIPSC population while the thin lines show subpopulations of mIPSCs based on the three Gaussian distributions distinguished in (A). (C) Cumulative probability plots of the mIPSC amplitudes in GFP− (black) and GFP+ (red) interneurons. Thick lines show the total mIPSC population while the thin lines show subpopulations of mIPSCs based on the three Gaussian distributions distinguished in (A). The arrow denotes the cumulative distribution of the amplitudes in the slow decaying subpopulation of GFP− interneurons. All plots consist of pooled data from 6 GFP+ and 6 GFP− neurons. The same number of consecutively occurring mIPSCs from each cell was used.
Table 2. Decay τ means (μ), standard deviations (SD) and weighting factor results from the mIPSC population analysis.
These results show that GABAergic synapses in SDH cells express a variety of decay kinetics and that the relative distribution of synapses with different kinetics is cell type dependent.
Variations in eIPSCs Decay Kinetics are not Due to Differences in GABA Reuptake
Although transmitter released by single vesicles is not always sufficient to diffuse to activate extrasynaptic receptors or spillover (Isaacson et al., 1993; Overstreet and Westbrook, 2003), it might occur in certain circumstances because of the presence of extrasynaptic high affinity receptors, as was shown for certain glutamatergic synapses (Diamond, 2001). To directly examine if diffusion of GABA and activation of extrasynaptic receptors could be responsible for the slow decay component of GABAA IPSCs, we looked at the effect of GABA uptake block. Evoked IPSCs were recorded in EGFP+ interneurons before and during the application of the GABA transporter blocker SKF89976A. Bath application of 100 μM SKF89976A did not significantly affect the decay time constant (τw, 143 ± 6.9 vs. 137.3 ± 4.9; Figure 5A) in EGFP+ interneurons. Similarly, in EGFP− interneurons, SKF89976A did not have a significant effect on the decay time constant (69.6 ± 5.3 vs.74.5 ± 5.1; Figure 5B). On the other hand SKF89976A caused significant prolongation of the currents recorded after a 20 Hz train stimulus (Figure 5A). These results suggest that under our experimental conditions, diffusion and spillover of GABA does not contribute in the shaping of the slow decay kinetics observed in the eIPSCs of the dorsal horn.
Figure 5. Slow GABAA IPSC decay kinetics in lamina II dorsal horn neurons are not explained by spillover or extrasynaptic receptor activation. (A) Application of the GABA uptake blocker SKF89976A (SKF) does not affect the decay time of single GABAA eIPSCs in GFP+ neurons. Traces on the left show average eIPSCs (0 mV holding potential) before (light gray) and during SKF application (black) recorded from a GFP+ neuron. On the right, the cumulative probability plots of eIPSC decay τw's from 5 GFP+ cells are shown (p > 0.05). The inset illustrates IPSCs evoked by a train of 12 stimuli at 20 Hz in the same cell (as in left) and shows a prolonged response after SKF application, indicating that block of uptake, in this case, revealed accumulation of GABA (the traces shown are normalized to the peak value). SKF prolonged the IPSC decay half-width after a train stimulus from 0.54 ± 0.08 to 0.65 ± 0.07 s (n = 5 GFP+ cells, p < 0.05). (B) Application of SKF does not affect the decay time of GABAA eIPSCs in GFP− neurons. Traces on the left show average eIPSCs before (light gray) and during SKF application (black) recorded from a GFP− neuron. On the right, the cumulative probability plots of eIPSC decay τw's from 4 GFP− cells are shown (p > 0.05). (C) Traces of mIPSC recordings (scale bars: 10 pA, 0.5 s) and examples of single mIPSCs (scale bars: 10 pA, 50 ms) recorded from α5−/− (upper) and wild type (lower) mice. The cumulative probability plot shows identical distribution of decay τw between knock-out (6 cells) and control mice (6 cells). (D) Traces of mIPSC recordings (scale bars: 10 pA, 0.5 s) and examples of single mIPSCs (scale bars: 10 pA, 50 ms) recorded from δ −/− (upper) and δ +/+ (lower) mice. The distribution of decay time constants is similar in δ −/− (6 cells) and δ+/+ (5 cells) as shown by the cumulative probability plot.
mIPSc Decay Kinetics are not Altered in α5−/− and δ−/− Mice
In additional experiments, we investigated the existence of slow decaying GABAA mIPSCs in mice lacking the α5 or δ subunits. GABAA receptors containing the α5 or δ subunits are high affinity receptors and are thought to be located extrasynaptically, mediating tonic inhibition (Mody and Pearce, 2004; Farrant and Nusser, 2005). We postulated that if the slowly decaying component of the synaptic current is due to activation of extrasynaptic receptors it might be mediated by either α5 or δ containing subunits. Figure 5C shows mIPSC recordings from α5−/− mice. Both mIPSCs with faster and prolonged decaying kinetics were observed. As shown in the cumulative probability plot the distribution of decay time constants in α5−/− mice did not differ from that of age-matched wild type mice. Similarly, slowly decaying mIPSCs were also observed in δ −/− mice (Figure 5D). As with the α5−/−, the decay time constant distribution in δ −/− was not significantly different from that of wild type littermates (δ +/+), confirming that GABA spillover and extrasynaptic receptor activation are not responsible for the slow decays.
Heterogeneity of Decay Kinetics is Replicated in Experiments with Rapid GABA Application to Spinal Lamina II Outside-Out Patches
The above results indicate that intrinsic receptor properties and not extrinsic factors, like neurotransmitter time course and diffusion, shape the decay kinetics of GABAA IPSCs. We sought to confirm this using conditions of controlled GABA application. Outside out membrane patches were isolated from lamina II neurons in mouse slices and exposed to brief pulses (1 ms) of GABA (1 mM) using a piezoelectric-controlled rapid application system. GABA induced the opening of several channels on the membrane patches (Figure 6A). Average traces for each patch were obtained from multiple brief pulse applications and the decay time constant was calculated by fitting a monoexponential decay function. Outside-out membrane patches displayed variable (Figure 6B) decay time constants ranging from 22.3 to 281.1 ms (n = 6), which corresponds to the range of decays we observed for mIPSCs. These data indicate that the intrinsic properties of GABAA channels in lamina II interneurons are sufficient to explain the diversity of IPSC decay kinetics, consistent with the above findings that the range of decays observed are not due to differences in transmitter release or GABA reuptake mechanisms.
Figure 6. Heterogeneity of decay kinetics of transient GABAA currents in outside-out patches. (A) Current traces showing responses to rapid application (1 ms pulse) of 1 mM GABA of varying decay times in different membrane patches. On top, the average of the liquid junction current recording is shown. In the middle, traces from a single application are shown. Average traces of multiple applications on the same membrane patches are shown at the bottom. (B) Graph showing the decay τ of individual excised patches. In black are shown outside-out patches recorded with high Cl− pipette solution and a holding of −70 mV while in red are recordings from excised patches made with low Cl− pipettes at 0 mV holding potential.
Simulations Reveal that Dendritic Filtering is not Sufficient to Explain the Variability in Decay Kinetics, however Diverse Subunit Composition could Yield the Observed Variability
The SDH consists of a variety of interneurons with different morphologies. Soma size and dendritic length can differ considerably between excitatory and inhibitory inteneurons (Yasaka et al., 2007). These morphological differences could influence the dendritic filtering of synaptic events. To identify how much of the decay variability could be attributed to dendritic filtering we used simulation experiments on a simple model neuron consisting of a soma and a single dendrite. We used a recorded event and fed it as a synaptic conductance at different distances on the dendrite. Using this approach we were able to yield filtered rise times that compared to those in our recorded mIPSCs (Figure 7A). Filtering also slowed the decay kinetics, however the resulting variability is only a fraction of the recorded variability. This indicates that dendritic filtering would influence decay kinetics to a limited extent and cannot explain the full variability of mIPSC kinetics in our experiments (Figure 7B).
Figure 7. Simulation of mIPSC decay kinetics. (A) The graph illustrates the changes in 10–90% rise time and decay τ in relation to the distance of the synapse from the soma in a model cell. A recorded fast mIPSC was used as a synaptic conductance. (B) Top, superimposed synaptic currents measured at the soma for synapses located at 0 μm (black) and 1000 μm (red) away from the soma. Traces were aligned at the peak. Bottom, superimposed traces from a synaptic current originating 1000 μm away from the soma of the model cell (red) and a slow decaying recorded mIPSC (green). (C) Schematic representation of the six-state Markov model of GABA receptor activation used (left). U, unbound; B, bound; O, open,; and D, desensitized states. A list of rate constants used to simulate fast and slow decay GABA currents is shown on the right. on1 and on2 binding rate constants are in M−1s−1 and the rest rate constant in s−1. (D) Top, computationally generated macroscopic current simulations for fast and slow decay parameters. Bottom, examples of simulated GABAA mIPSCs resulting from the summation of stochastic opening of 20 single channels. (E) Population histogram showing the distribution of decay τ for fast (black) and slow (red) simulated mIPSCs.
Subunit composition of GABAA receptors can be a major factor in determining IPSC decay kinetics. A prominent example of cell type-specific and subunit-dependent differences in kinetic properties has been shown in the thalamus (Browne et al., 2001; Schofield and Huguenard, 2007). Inhibitory neurons from the nucleus reticularis display IPSCs with characteristic slow decay kinetics, while relay neurons in the ventrobasal area display fast decaying IPSCs, a difference attributable to GABAA receptor affinity (Browne et al., 2001; Schofield and Huguenard, 2007). This raises the question of whether differences in the binding characteristics that result from distinct subunit composition are sufficient to explain the wide range of the mIPSC decay kinetics observed in spinal lamina II neurons. To test this, we used a six-state Markov model of channel gating to simulate mIPSC decays. This model (Figure 7C) was structurally similar to those previously described (Jones and Westbrook, 1995; Schofield and Huguenard, 2007). It included two binding states (mono-liganded B1 and bi-liganded B2), two open and two desensitized states. We first used rate constants similar to ones previously reported (Schofield and Huguenard, 2007) and modified them to obtain monoexponential macroscopic currents with a decay time constant equivalent to the fast decay component of our eIPSCs (average τ1 from EGFP+ and EGFP− eIPSCs; Table 1, high Cl− recordings). We then reused the same parameters, but modified only GABA binding and unbinding rate constants to generate macroscopic currents with slow decay kinetics comparable to that of our eIPSCs (average τ2 from the eIPSCs in the two interneuron populations). Using these parameters, 100 “fast” and 100 “slow” mIPSCs were simulated each from stochastic opening of 20 single channels (Figure 7D). Figure 7E shows the distribution of decay time constants of the simulated mIPSCs, yielding a similar diversity of kinetics to that of mIPSCs recorded in lamina II interneurons. These results indicate that a difference in intrinsic parameters, such as binding and unbinding rates, is sufficient to replicate the heterogeneity of mIPSC observed among lamina II interneurons as was observed in the thalamus (Schofield and Huguenard, 2007). While not definitive proof, it is consistent with a differential subunit composition at different GABAA synapses.
Distinct α Subunit Composition Explains the Differential Decay Kinetics of GABAA IPSCs
To test whether differences in α subunit composition of GABAA receptors could explain the heterogeneity of IPSC decay kinetics, we exploited three lines of knock-in mice with altered α-subunit pharmacology. In these knock-in mice a histidine residue is substituted by an arginine at the benzodiazepine binding site of the α1 [α1(H101R)], α2 [α2(H101R)] or α3 [α3(H126R)] subunits (Rudolph et al., 1999; Low et al., 2000). The mutated subunits form functional GABAA receptors (Benson et al., 1998) but are insensitive to diazepam modulation. On the other hand, Ro 15-4513 which acts as a partial inverse agonist on wild type GABAA receptors, acts as an agonist on mutated benzodiazepine binding site. Thus, it selectively potentiates responses of GABAA receptors that contain the point mutated α-subunits (Benson et al., 1998). This provided the missing pharmacological tool to test for the contribution of each of these subunits to the different kinetic components of GABAA IPSCs.
In spinal slices taken form wild type and each of these knock-in mice, we recorded mIPSCs before (CTRL) and during the application of 1 μM Ro 15-4513 (Figure 8A). Amplitude, decay τw and frequency of the mIPSCs in CTRL did not differ between wild type and knock-in mice (not shown). To investigate decay-kinetic specific differences in the effect of Ro 15-4513 we categorized mIPSCs as fast, if their decay τw was <100 ms or as slow for mIPSCs with a τw > 100 ms. Benzodiazepines can modulate the amplitude, the decay time kinetics as well as the frequency of mIPSCs (in cases where it unmasks perisynaptic receptors for example, Chery and De Koninck, 1999). Hence, to measure the overall effect of Ro 15-4513, for each category of mIPSCs (fast and slow), we calculated the charge transfer for each mIPSC and then the sum of charge transfer for all mIPSCs in each category occurring over a period 10 min of CTRL and Ro 15-4513 conditions. We expressed the result as the ratio of Ro15-4513 to CTRL (Δ charge).
Figure 8. Differential role of α-subunits in shaping GABAA IPSC decay kinetics. Effect of the partial benzodiazepine inverse agonist Ro15-4513 on fast and slow mIPSCs in control (wt), α1(H101R), α2(H101R), and α3(H126R) knock-in mice. (A) Example traces of mIPSC recordings before (ctrl) and during Ro 15-4513 in the different mouse lines. (B) The relative change in total mIPSC charge transfer (Δ charge-transfer) mediated by fast decay mIPSCs after the application of Ro15-4513. (C) The relative change in the Δ charge-transfer mediated by slow decay mIPSCs after the application of Ro15-4513. (D) The relative changes in mIPSC decay τ, amplitude and frequency for fast decaying mIPSCs after the application of Ro15-4513. (E) The relative changes in mIPSC decay τ, amplitude and frequency for slow decaying mIPSCs after the application of Ro15-4513. Asterisks denote significant difference between cells taken from wild type and knock-in mice (ANOVA with post-hoc Tukey test, p < 0.05, n = 9, 6, 6, and 6 respectively). (F) The cumulative probability plot of mIPSC decay τw for α3(H126R) knock-in mice (left) is shown before (ctrl) and during Ro 15-4513. Application of Ro 15-4513, shifted the decay τ distribution to the right, as expected by the increase in frequency of the slowly decaying mIPSCs subpopulation. Indeed, the cumulative probability plots for the fast (middle) and slow (right) subpopulations of the α3(H126R) mIPSC decay τw show no differences between ctrl and Ro15 in their distribution. All plots consist of pooled data from six α3(H126R) neurons. The same number of consecutively occurring mIPSCs from each cell was used.
Comparison of Ro 15-4513 effect in the four groups of mice revealed a significant potentiation (p < 0.05) of the relative Δ charge contributed by fast mIPSCs in α1(H101R) and α2(H101R) mice, but not in α3(H126R) mice (p > 0.05) when compared with wt mice Figure 8B). In contrast, Ro 15-4513 significantly potentiated (p < 0.05) the relative Δ charge contributed by slow mIPSCs only in α3(H126R) mice (Figure 8C). The Ro 15-4513 effects on Δ charge for both slow and fast mIPSCs was mostly due to changes in event frequency (Figures 8D–F). These results indicate that GABAA receptors containing α1 and α2 subunits are majorly responsible for shaping fast decaying mIPSCs while α3 underlies slow decay mIPSCs.
Here we demonstrate that GABAA transmission is distinct in lamina II interneuron subpopulations as shown by the slower kinetics of evoked GABAA IPSCs in inhibitory interneurons. Detailed analysis of mIPSC decay times showed that although both EGFP+ and EGFP− interneurons consist of mIPSC populations with similar decay kinetics, the relative distribution of faster and slower decay kinetics is distinct in the two populations. While in EGFP− interneurons are dominated by faster decaying mIPSCs, in EGFP+ interneurons slow decaying mIPSCs were more frequent on a relative scale. Furthermore, our data indicate that these differences reflect differential distribution of GABAA receptor subunits.
Several lines of evidence suggest that the slower decaying mIPSCs are not the result of electrotonic filtering. The 10–90% rise time and decay time constants of mIPSCs do not correlate. When mIPSCs are subdivided to populations of faster, medium and slower decaying events each subpopulation shows similarly distributed rise times for each interneuron subtype. In addition, the existence of mIPSCs with mixed decay kinetics suggests that both components arise from the release of a single vesicle and thus the same electrotonic distance. Finally, in experiments from isolated membrane patches and rapid drug application a wide range in decay kinetics is observed. Thus, it is improbable that the observed variability in IPSC decay kinetics is due to dendritic filtering. However, we also observed slower mIPSC rise times in EGFP+ compared to EGFP− interneurons. This could be due to differences in electrotonic filtering in the two interneuron populations as a consequence of differences in morphology (Yasaka et al., 2007). However, other factors could also influence rise time, including differences in synaptic geometry or distinct subunit composition. Recent work in the rat dorsal horn using GABA and glycine uncaging suggests that the distribution of inhibitory neurotransmitter receptors in islet cells is mainly on the soma and proximal dendrites (Kato et al., 2007). In the same cells, electrotonic filtering of transients caused by glutamate uncaging only modestly affected decay time constants, providing evidence of electrically compact cells (Kato et al., 2007). Analogous results were provided by Chery and De Koninck (1999) for lamina I neurons where similar rise and decay kinetics were observed for both proximally and distally evoked IPSCs. In agreement with this the glycinergic eIPSC decay time constants in our experiments, as well as amplitudes, were equivalent in the two populations. Nonetheless, filtering could alter mIPSCs properties differently in EGFP− and EGFP+ interneurons. On the other hand, the similar rise time distribution for different decay kinetic populations, indicate that such filtering would affect the whole range of mIPSCs uniformly, thus it is unlikely it contributes to the differences in the mIPSC variability we observed between the two interneuron populations.
Slowly decaying GABAA IPSCs have been observed before in several CNS areas. In the spinal ventral horn, Renshaw cells have significantly slower GABAA IPSC decay kinetics than non-Renshaw cells which correlates with expression α3/α5-containing GABAA receptors in these neurons (Geiman et al., 2002; Gonzalez-Forero and Alvarez, 2005). In the thalamus, GABAergic neurons in the reticular nucleus display slow decay IPSCs that are thought to arise from the expression of α3 subunits that confer lower receptor affinity for GABA (Browne et al., 2001; Schofield and Huguenard, 2007). Slow GABAA IPSCs in hippocampal pyramidal neurons (Pearce, 1993), on the other hand, arise from low concentration and long lasting GABA transients that are evoked by neurogliaform cells (Karayannis et al., 2010) which also involve the activation of α5 containing GABAA receptors (Zarnowska et al., 2009). Moreover, in the cerebellum, granule cells receive both fast and slow decay GABAA IPSCs from the same type interneuron, the Golgi cell. While fast transients arise from activation of synaptic α1 containing receptors, slow transients involve spillover and the activation of high affinity α6 containing receptors (Rossi and Hamann, 1998). Neocortical low threshold-spiking (LTS) interneurons display slower IPSC decay kinetics than fast spiking (FS) interneurons, which correlates with the presence of α1 subunit in FS but lack of this subunit in LTS interneurons (Bacci et al., 2003). Finally, in pyramidal cortical neurons, similar to their hippocampal counterparts, low-concentration, slow GABAA transients from neuroglia yield IPSCs with slow decay kinetics (Szabadics et al., 2007). These paradigms show that a multitude of mechanisms are involved in creating GABAA IPSC kinetic diversity, including subunit composition but also variable synaptic arrangements yielding a wide variety of flavors of inhibition.
GABA clearance from the cleft by uptake transporters determines the duration of the neurotransmitter availability, its concentration and diffusion out of the cleft. All of these parameters also shape the decay kinetics of GABAA IPSCs (Isaacson et al., 1993; Nusser et al., 2001; Overstreet and Westbrook, 2003). Our experiments indicate that the neurotransmitter time course in the synapse plays little role in shaping the decay time course of IPSCs in spinal lamina II, since pharmacological blockade of GABA uptake did not affect the kinetics of single eIPSCs. The presence of slow decay mIPSCs also argues against spillover. The occurrence of slow mIPSCs in knockouts of the high affinity α5 and δ subunits is consistent with these conclusions. Finally, we confirmed the role of intrinsic receptor properties as sufficient to explain the heterogeneity of decay kinetics with outside out membrane patches under controlled GABA applications.
The subunit composition of GABAA receptors is a major determinant of their pharmacological and biophysical properties (Macdonald and Olsen, 1994; Gingrich et al., 1995; Olsen and Sieghart, 2009). Immunocytochemical and in situ hybridization studies in the spinal cord have shown the expression of several subunits, including α1, α2, α3, and α5, the β 2 and β 3, and the γ2 subunit (Persohn et al., 1991; Wisden et al., 1991; Ma et al., 1993; Bohlhalter et al., 1996; Todd et al., 1996). Furthermore, in recent functional studies, the spinal α2 subunit and, to a lesser extent, α3 and α5 have been implicated in the analgesic effect of diazepam in neuropathic and inflammatory pain (Knabl et al., 2008). The distinct effects of Ro 15-4513 on the fast and slow decay mIPSCs in the three knock-in mice establishes the GABAA receptor composition as the main factor contributing to the heterogeneity in IPSC decay kinetics in spinal lamina II. Consistent with our findings, fast decay IPSCs have been previously associated with the α1 subunit in other brain areas (Browne et al., 2001; Vicini et al., 2001; Bacci et al., 2003). Also consistent with our findings, α3 subunits have been shown to underlie slow decay IPSCs in the ventrobasal thalamus (Browne et al., 2001; Schofield and Huguenard, 2007). In addition to these observations, further variability in GABAA receptor kinetics may result from receptor phosphorylation (Nusser et al., 1999) or from modulation by endogenous neurosteroids (Keller et al., 2004; Poisbeau et al., 2005). In the SDH a regional variability of GABAA decay kinetics was shown to be a result of local differences in steroid synthesis (Inquimbert et al., 2008). Variation in menstrual cycle hormones might also contribute to an additional variation in kinetics (Poisbeau et al., 2014).
The heterogeneity in kinetics between different cell types and even within the same neuron at different synapses, implies multiple functionality of the GABAA synapse. Indeed, distinct GABAA receptor kinetics can differentially affect information processing. While fast inhibitory currents can modulate the threshold of input-output transfer (Crowley et al., 2009), slow and tonic inhibition can influence both the slope (gain) and threshold of input-output transfer (Mitchell and Silver, 2003; Prescott and De Koninck, 2003; Crowley et al., 2009). Consequently, activation of a subset of synapses with a specific subunit composition can result in distinct processing and information flow. This raises the question of whether synapses with certain subunit composition and kinetics (fast or slow) receive specialized input. At present our understanding of the information processing in the SDH is too poor to conclude on this issue because several key pieces of information on the structural and functional organization in this area are still missing and further investigations are warranted.
Taken together the present results provide evidence that distinct subunit composition underlies the heterogeneity in time course of inhibition in the dorsal horn of the spinal cord. This finding may have important implications for pain control and drug design, as they may contribute to our understanding of the net impact drugs may have at the network level. In recent studies, molecular dissection of the analgesic effect of diazepam has shown that it is mainly mediated by the α2 and, to a lesser extent, α3 subunits at the spinal level (Knabl et al., 2008). These results open new avenues for specific pharmacological targeting of GABAA receptors subtypes for pain treatment. Yet, the α3 subunit did not appear to contribute to benzodiazepine-mediated analgesia to the same degree in the different pain tests performed (Knabl et al., 2008). In this context, it is interesting to note that, from our results, the α3 subunit appears to dominate inhibition of inhibitory interneurons. Potentiation of this subunit on these interneurons may thus cause significant disinhibition. Because distinct elements of the dorsal horn circuitry may be affected differentially in each chronic pain condition and in each sensory response, the net effect of modulating the α3 subunit may differ. In certain cases, for example, the “pro-algesic” effect of disinhibition may counteract any concomitant analgesic action resulting from inhibition of dorsal horn excitatory circuits. Thus, our results suggest that selective targeting of the α2 subunit may yield a more broadly effective analgesic.
Charalampos Labrakakis and Yves De Koninck conceived the study. Charalampos Labrakakis designed, performed and analyzed the experiments. All authors contributed to the interpretation of data and preparation of the final manuscript.
Conflict of Interest Statement
The authors would like to thank Dr. D. Bowie for advice on rapid agonist application, Drs. F. Erdelyi and G. Szabo for the generous gift of GAD-GFP mice and K. Vandal for excellent technical assistance. This work was supported by the Canadian Institute of Health Research (CIHR) grant MOP-12942 to Yves De Koninck. Charalampos Labrakakis was supported by a CIHR postdoctoral fellowship. Yves De Koninck was a Chercheur National of the Fonds de la recherche en santé du Quebec (FRSQ).
Baba, H., Ji, R. R., Kohno, T., Moore, K. A., Ataka, T., Wakai, A., et al. (2003). Removal of GABAergic inhibition facilitates polysynaptic A fiber-mediated excitatory transmission to the superficial spinal dorsal horn. Mol. Cell Neurosci. 24, 818–830. doi: 10.1016/S1044-7431(03)00236-7
Benson, J. A., Low, K., Keist, R., Mohler, H., and Rudolph, U. (1998). Pharmacology of recombinant gamma-aminobutyric acida receptors rendered diazepam-insensitive by point-mutated alpha-subunits. Febs. Lett. 431, 400–404. doi: 10.1016/S0014-5793(98)00803-5
Bohlhalter, S., Weinmann, O., Mohler, H., and Fritschy, J. M. (1996). Laminar compartmentalization of GABA(A)-receptor subtypes in the spinal cord: an immunohistochemical study. J. Neurosci. 16, 283–297.
Browne, S. H., Kang, J., Akk, G., Chiang, L. W., Schulman, H., Huguenard, J. R., et al. (2001). Kinetic and pharmacological properties of GABA(A) receptors in single thalamic neurons and GABA(A) subunit expression. J. Neurophysiol. 86, 2312–2322.
Chery, N., and De Koninck, Y. (1999). Junctional versus extrajunctional glycine and GABA(A) receptor-mediated IPSCs in identified lamina I neurons of the adult rat spinal cord. J. Neurosci. 19, 7342–7355.
Chery, N., Yu, X. H., and De Koninck, Y. (2000). Visualization of lamina I of the dorsal horn in live adult rat spinal cord slices. J. Neurosci. Methods 96, 133–142. doi: 10.1016/S0165-0270(99)00195-8
Collinson, N., Kuenzi, F. M., Jarolimek, W., Maubach, K. A., Cothliff, R., Sur, C., et al. (2002). Enhanced learning and memory and altered GABAergic synaptic transmission in mice lacking the alpha 5 subunit of the GABA(A) receptor. J. Neurosci. 22, 5572–5580.
Cordero-Erausquin, M., Allard, S., Dolique, T., Bachand, K., Ribeiro-da-Silva, A., and De, K. Y. (2009). Dorsal horn neurons presynaptic to lamina I spinoparabrachial neurons revealed by transynaptic labeling. J. Comp. Neurol. 517, 601–615. doi: 10.1002/cne.22179
Cossart, R., Bernard, C., and Ben Ari, Y. (2005). Multiple facets of GABAergic neurons and synapses: multiple fates of GABA signalling in epilepsies. Trends Neurosci. 28, 108–115. doi: 10.1016/j.tins.2004.11.011
Coull, J. A., Boudreau, D., Bachand, K., Prescott, S. A., Nault, F., Sik, A., et al. (2003). Trans-synaptic shift in anion gradient in spinal lamina I neurons as a mechanism of neuropathic pain. Nature 424, 938–942. doi: 10.1038/nature01868
Crowley, J. J., Fioravante, D., and Regehr, W. G. (2009). Dynamics of fast and slow inhibition from cerebellar golgi cells allow flexible control of synaptic integration. Neuron 63, 843–853. doi: 10.1016/j.neuron.2009.09.004
Geiman, E. J., Zheng, W., Fritschy, J. M., and Alvarez, F. J. (2002). Glycine and GABA(A) receptor subunits on Renshaw cells: relationship with presynaptic neurotransmitters and postsynaptic gephyrin clusters. J. Comp. Neurol. 444, 275–289. doi: 10.1002/cne.10148
Gingrich, K. J., Roberts, W. A., and Kass, R. S. (1995). Dependence of the GABA(A) receptor gating kinetics on the alpha-subunit isoform: implications for structure-function relations and synaptic transmission. J. Physiol. (Lond). 489, 529–543.
Gonzalez-Forero, D., and Alvarez, F. J. (2005). Differential postnatal maturation of GABAA, glycine receptor, and mixed synaptic currents in Renshaw cells and ventral spinal interneurons. J. Neurosci. 25, 2010–2023. doi: 10.1523/JNEUROSCI.2383-04.2005
Hammond, D. L., and Drower, E. J. (1984). Effects of intrathecally administered THIP, baclofen and muscimol on nociceptive threshold. Eur. J. Pharmacol. 103, 121–125. doi: 10.1016/0014-2999(84)90197-3
Harvey, R. J., Depner, U. B., Wassle, H., Ahmadi, S., Heindl, C., Reinold, H., et al. (2004). GlyR alpha 3: an essential target for spinal PGE(2)-mediated inflammatory pain sensitization. Science 304, 884–887. doi: 10.1126/science.1094925
Inquimbert, P., Rodeau, J. L., and Schlichter, R. (2008). Regional differences in the decay kinetics of GABA(A) receptor-mediated miniature IPSCs in the dorsal horn of the rat spinal cord are determined by mitochondrial transport of cholesterol. J. Neurosci. 28, 3427–3437. doi: 10.1523/JNEUROSCI.5076-07.2008
Karayannis, T., Elfant, D., Huerta-Ocampo, I., Teki, S., Scott, R. S., Rusakov, D. A., et al. (2010). Slow GABA transient and receptor desensitization shape synaptic responses evoked by hippocampal neurogliaform cells. J. Neurosci. 30, 9898–9909. doi: 10.1523/JNEUROSCI.5883-09.2010
Kato, G., Kawasaki, Y., Ji, R. R., and Strassman, A. M. (2007). Differential wiring of local excitatory and inhibitory synaptic inputs to islet cells in rat spinal lamina II demonstrated by laser scanning photostimulation. J. Physiol. (Lond). 580, 815–833. doi: 10.1113/jphysiol.2007.128314
Keller, A. F., Beggs, S., Salter, M. W., and De Koninck, Y. (2007). Transformation of the output of spinal lamina I neurons after nerve injury and microglia stimulation underlying neuropathic pain. Mol. Pain 3:27. doi: 10.1186/1744-8069-3-27
Keller, A. F., Breton, J. D., Schlichter, R., and Poisbeau, P. (2004). Production of 5 alpha-reduced neurosteroids is developmentally regulated and shapes GABA(A) miniature IPSCs in lamina II of the spinal cord. J. Neurosci. 24, 907–915. doi: 10.1523/JNEUROSCI.4642-03.2004
Knabl, J., Witschi, R., Hosl, K., Reinold, H., Zeilhofer, U. B., Ahmadi, S., et al. (2008). Reversal of pathological pain through specific spinal GABA(A) receptor subtypes. Nature 451, 330–334. doi: 10.1038/nature06493
Labrakakis, C., Lorenzo, L. E., Bories, C., Ribeiro-da-Silva, A., and De Koninck, Y. (2009). Inhibitory coupling between inhibitory interneurons in the spinal cord dorsal horn. Mol. Pain 5:24. doi: 10.1186/1744-8069-5-24
Laing, I., Todd, A. J., Heizmann, C. W., and Schmidt, H. H. (1994). Subpopulations of GABAergic neurons in laminae I-III of rat spinal dorsal horn defined by coexistence with classical transmitters, peptides, nitric oxide synthase or parvalbumin. Neuroscience 61, 123–132. doi: 10.1016/0306-4522(94)90065-5
Lopez-Bendito, G., Sturgess, K., Erdelyi, F., Szabo, G., Molnar, Z., and Paulsen, O. (2004). Preferential origin and layer destination of GAD65-GFP cortical interneurons. Cereb. Cortex 14, 1122–1133. doi: 10.1093/cercor/bhh072
Low, K., Crestani, F., Keist, R., Benke, D., Brunig, I., Benson, J. A., et al. (2000). Molecular and neuronal substrate for the selective attenuation of anxiety. Science 290, 131–134. doi: 10.1126/science.290.5489.131
Ma, W., Saunders, P. A., Somogyi, R., Poulter, M. O., and Barker, J. L. (1993). Ontogeny of GABA(A) receptor subunit mRNAs in rat spinal cord and dorsal root ganglia. J. Comp. Neurol. 338, 337–359. doi: 10.1002/cne.903380303
Mihalek, R. M., Banerjee, P. K., Korpi, E. R., Quinlan, J. J., Firestone, L. L., Mi, Z. P., et al. (1999). Attenuated sensitivity to neuroactive steroids in gamma-aminobutyrate type A receptor delta subunit knockout mice. Proc. Natl. Acad. Sci. U.S.A. 96, 12905–12910. doi: 10.1073/pnas.96.22.12905
Nusser, Z., Naylor, D., and Mody, I. (2001). Synapse-specific contribution of the variation of transmitter concentration to the decay of inhibitory postsynaptic currents. Biophys. J. 80, 1251–1261. doi: 10.1016/S0006-3495(01)76101-2
Nusser, Z., Sieghart, W., and Mody, I. (1999). Differential regulation of synaptic GABA(A) receptors by cAMP-dependent protein kinase in mouse cerebellar and olfactory bulb neurones. J. Physiol. (Lond). 521, 421–435. doi: 10.1111/j.1469-7793.1999.00421.x
Persohn, E., Malherbe, P., and Richards, J. G. (1991). In situ hybridization histochemistry reveals a diversity of GABAA receptor subunit mRNAs in neurons of the rat spinal cord and dorsal root ganglia. Neuroscience 42, 497–507. doi: 10.1016/0306-4522(91)90392-2
Poisbeau, P., Keller, A. F., Aouad, M., Kamoun, N., Groyer, G., and Schumacher, M. (2014). Analgesic strategies aimed at stimulating the endogenous production of allopregnanolone. Front. Cell Neurosci. 8:174. doi: 10.3389/fncel.2014.00174
Poisbeau, P., Patte-Mensah, C., Keller, A. F., Barrot, M., Breton, J. D., Luis-Delgado, O. E., et al. (2005). Inflammatory pain upregulates spinal inhibition via endogenous neurosteroid production. J. Neurosci. 25, 11768–11776. doi: 10.1523/JNEUROSCI.3841-05.2005
Prescott, S. A., and De Koninck, Y. (2003). Gain control of firing rate by shunting inhibition: roles of synaptic noise and dendritic saturation. Proc. Natl. Acad. Sci. U.S.A. 100, 2076–2081. doi: 10.1073/pnas.0337591100
Ribeiro-da-Silva, A., and De Koninck, Y. (2008). “Morphological and neurochemical organization of the spinal dorsal horn,” in The Science of Pain, eds A. I. Basbaum and M. C. Bushnell (New York, NY: Academic Press), 279–310.
Rossi, D. J., and Hamann, M. (1998). Spillover-mediated transmission at inhibitory synapses promoted by high affinity alpha(6) subunit GABA(A) receptors and glomerular geometry. Neuron 20, 783–795. doi: 10.1016/S0896-6273(00)81016-8
Rudolph, U., Crestani, F., Benke, D., Brunig, I., Benson, J. A., Fritschy, J. M., et al. (1999). Benzodiazepine actions mediated by specific gamma-aminobutyric acid(A) receptor subtypes. Nature 401, 796–800. doi: 10.1038/44579
Sivilotti, L., and Woolf, C. J. (1994). The contribution of GABAA and glycine receptors to central sensitization: disinhibition and touch-evoked allodynia in the spinal cord. J. Neurophysiol. 72, 169–179.
Sorkin, L. S., and Puig, S. (1996). Neuronal model of tactile allodynia produced by spinal strychnine: effects of excitatory amino acid receptor antagonists and a mu-opiate receptor agonist. Pain 68, 283–292. doi: 10.1016/S0304-3959(96)03130-2
Szabadics, J., Tamas, G., and Soltesz, I. (2007). Different transmitter transients underlie presynaptic cell type specificity of GABA(A,slow) and GABA(A,fast). Proc. Natl. Acad. Sci. U.S.A. 104, 14831–14836. doi: 10.1073/pnas.0707204104
Todd, A. J., and Sullivan, A. C. (1990). Light microscope study of the coexistence of GABA-like and glycine-like immunoreactivities in the spinal cord of the rat. J. Comp. Neurol. 296, 496–505. doi: 10.1002/cne.902960312
Torsney, C., and MacDermott, A. B. (2006). Disinhibition opens the gate to pathological pain signaling in superficial neurokinin 1 receptor-expressing neurons in rat spinal cord. J. Neurosci. 26, 1833–1843. doi: 10.1523/JNEUROSCI.4584-05.2006
Vicini, S., Ferguson, C., Prybylowski, K., Kralic, J., Morrow, A. L., and Homanics, G. E. (2001). GABA(A) receptor alpha 1 subunit deletion prevents developmental changes of inhibitory synaptic currents in cerebellar neurons. J. Neurosci. 21, 3009–3016.
Wisden, W., Gundlach, A. L., Barnard, E. A., Seeburg, P. H., and Hunt, S. P. (1991). Distribution of Gaba(A) receptor subunit messenger-rnas in rat lumbar spinal-cord. Brain. Res. Mol. Brain Res. 10, 179–183. doi: 10.1016/0169-328X(91)90109-B
Yaksh, T. L. (1989). Behavioral and autonomic correlates of the tactile evoked allodynia produced by spinal glycine inhibition: effects of modulatory receptor systems and excitatory amino acid antagonists. Pain 37, 111–123. doi: 10.1016/0304-3959(89)90160-7
Yasaka, T., Kato, G., Furue, H., Rashid, M. H., Sonohata, M., Tamae, A., et al. (2007). Cell-type-specific excitatory and inhibitory circuits involving primary afferents in the substantia gelatinosa of the rat spinal dorsal horn in vitro. J. Physiol. 581, 603–618. doi: 10.1113/jphysiol.2006.123919
Zarnowska, E. D., Keist, R., Rudolph, U., and Pearce, R. A. (2009). GABA(A) receptor alpha 5 subunits contribute to GABA(A, slow) synaptic inhibition in mouse hippocampus. J. Neurophysiol. 101, 1179–1191. doi: 10.1152/jn.91203.2008
Keywords: GABAA receptors, IPSCs, decay kinetics, spinal dorsal horn, subunit composition
Citation: Labrakakis C, Rudolph U and De Koninck Y (2014) The heterogeneity in GABAA receptor-mediated IPSC kinetics reflects heterogeneity of subunit composition among inhibitory and excitatory interneurons in spinal lamina II. Front. Cell. Neurosci. 8:424. doi: 10.3389/fncel.2014.00424
Received: 12 September 2014; Accepted: 24 November 2014;
Published online: 11 December 2014.
Edited by:Arianna Maffei, Suny Stony Brook, USA
Copyright © 2014 Labrakakis, Rudolph and De Koninck. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yves De Koninck, Unité de Neurosciences Cellulaires et Moléculaire, Institut Universitaire en Santé Mentale de Québec, 2601, Chemin de la Canardière, Québec, QC G1J 2G3, Canada e-mail: firstname.lastname@example.org