Abstract
The corelease of several neurotransmitters from a single synaptic vesicle has been observed at many central synapses. Nevertheless, the signaling synergy offered by cotransmission and the mechanisms that maintain the optimal release and detection of neurotransmitters at mixed synapses remain poorly understood, thus limiting our ability to interpret changes in synaptic signaling and identify molecules important for plasticity. In the brainstem and spinal cord, GABA and glycine cotransmission is facilitated by a shared vesicular transporter VIAAT (also named VGAT), and occurs at many immature inhibitory synapses. As sensory and motor networks mature, GABA/glycine cotransmission is generally replaced by either pure glycinergic or GABAergic transmission, and the functional role for the continued corelease of GABA and glycine is unclear. Whether or not, and how, the GABA/glycine content is balanced in VIAAT-expressing vesicles from the same terminal, and how loading variability effects the strength of inhibitory transmission is not known. Here, we use a combination of loose-patch (LP) and whole-cell (WC) electrophysiology in cultured spinal neurons of GlyT2:eGFP mice to sample miniature inhibitory post synaptic currents (mIPSCs) that originate from individual GABA/glycine co-releasing synapses and develop a modeling approach to illustrate the gradual change in mIPSC phenotypes as glycine replaces GABA in vesicles. As a consistent GABA/glycine balance is predicted if VIAAT has access to both amino-acids, we test whether vesicle exocytosis from a single terminal evokes a homogeneous population of mixed mIPSCs. We recorded mIPSCs from 18 individual synapses and detected glycine-only mIPSCs in 4/18 synapses sampled. The rest (14/18) were co-releasing synapses that had a significant proportion of mixed GABA/glycine mIPSCs with a characteristic biphasic decay. The majority (9/14) of co-releasing synapses did not have a homogenous phenotype, but instead signaled with a combination of mixed and pure mIPSCs, suggesting that there is variability in the loading and/or storage of GABA and glycine at the level of individual vesicles. Our modeling predicts that when glycine replaces GABA in synaptic vesicles, the redistribution between the peak amplitude and charge transfer of mIPSCs acts to maintain the strength of inhibition while increasing the temporal precision of signaling.
Introduction
A well-established and long accepted postulate known as “Dale’s principle” proposes that neurons have chemical unity for signaling (Strata and Harvey, ; Tritsch et al., ) —meaning that they communicate by releasing the same neurotransmitter at all their synapses. When the quantal nature of synaptic transmission was later discovered, application of Dale’s principle implied that each presynaptic terminal shares the particular set of enzymes and transporters that are required for the synthesis, recapture and accumulation of a specific neurotransmitter(s) into the synaptic vesicles, which are recycled and refilled locally before being reused (Hnasko and Edwards, ). Now it is clear that many central synapses corelease multiple neurotransmitters, expanding the repertoire of chemical signaling (Vaaga et al., ; Granger et al., ) and the extent of presynaptic coordination require to uphold Dale’s principal.
In the hindbrain, GABA and glycine are coreleased from single vesicles (Jonas et al., ; Lu et al., ). Although GABA and glycine have separate presynaptic supply mechanisms, postsynaptic receptors, and modulators, they share the same vesicular transporter, VIAAT (Wojcik et al., ; Aubrey et al., ) and when GABA/glycine are coreleased from the same single vesicle, mixed miniature inhibitory postsynaptic currents (mIPSCs) have a characteristic biphasic time course (Jonas et al., ; Keller et al., ; Russier et al., ). The regulation of GABA/glycine corelease is known to play a significant role during the maturation of spinal inhibitory circuits (Chéry and de Koninck, ; Keller et al., ; Russier et al., ; Coull et al., ; Rahman et al., ; Medelin et al., ), as well as in other hind-brain regions including the auditory nucleus (Kotak et al., ; Lu et al., ; Fischl and Burger, ; Nerlich et al., ) and cerebellum (Dugué et al., ; Rousseau et al., ). Early biochemical examinations of GABA and glycine uptake into synaptosomes (an isolated synaptic terminal preparation) indicated that both amino acids compete for the same vesicular transporter (Fykse and Fonnum, ; Burger et al., ; Christensen and Fonnum, ). This idea was endorsed when VIAAT was found in both GABAergic and glycinergic terminals (McIntire et al., ; Sagné et al., ) and taken together, these data suggested that the vesicular content is determined by the relative occupancy of GABA/glycine at VIAAT (Gasnier, ) and therefore, will be primarily determined by the relative presynaptic cytosolic concentrations of the two inhibitory neurotransmitters.
In active terminals, GABA is synthesized from glutamate by two isoforms of the glutamate decarboxylase enzyme (GAD65/67; Martin and Rimvall, ; Kakizaki et al., ), whereas the concentrative power of GlyT2, a 3 Na+-coupled glycine transporter (Roux and Supplisson, ), is required to increase basal glycine concentration to the appropriate levels for VIAAT uptake (Gomeza et al., ; Rousseau et al., ; Apostolides and Trussell, ). Experimental evidence confirms that the cytosolic concentration of glycine and GABA does indeed influence synaptic vesicle content (Rousseau et al., ; Apostolides and Trussell, ; Ishibashi et al., ). Furthermore, two invitro studies have succeeded in sampling the inhibitory transmitter content of VIAAT-expressing vesicles at single terminals and found evidence that the mIPSCs phenotypes were not homogenous. Instead, GABA-, glycine- and mixed mIPSCs were all detected at some individual synapses, suggesting that the signaling phenotype of a single terminal may considerably vary from vesicle to vesicle (Katsurabayashi et al., ; Aubrey et al., ; Supplementary Figure S1).
Here, we record mIPSCs in networks of culture spinal cord neurons and extract the subpopulation of mIPSCs that originate from an individual terminal by simultaneously recording mIPSCs with a loose-cell patch-clamp electrode placed over a single pre-synaptic varicosity. Then, we examine whether vesicles originating from the same co-releasing terminal evoke homogenous postsynaptic currents, as expected if all of a presynaptic terminal’s vesicles contain a similar concentration of GABA and glycine. In addition, we develop a simulation model of cotransmission to examine how mIPSC peak amplitude and charge transfer change as function of glycine vs. GABA release. Using activation kinetic schemes of GlyR and GABAAR, our simulations show that alterations in vesicular content would not compromised the strength, but rather would shape the time-course of postsynaptic inhibition.
Materials and Methods
Embryonic Mouse Spinal Cord Neurons
Primary cultures of spinal cord neurons were prepared as described in Hanus et al. () from embryonic day 13 or 14 (E13–14) C57BL/6J wild-type or heterozygous GlyT2-EGFP mouse pups (Zeilhofer et al., ). Embryos were obtained by cesarean section from pregnant mice anesthetized by intraperitoneal injection of ketamine-xylazine (100 and 10 mg/kg) and killed by cervical dislocation. Spinal cords were dissected under sterile conditions into PBS with 33 mM glucose at pH = 7.4 and then incubated in trypsin/EDTA solution (0.05% v/v, Sigma, St. Louis, MO, USA) for 10 min at 37°C. Cells were dissociated mechanically in a modified L15 Leibowitz’s medium (Life Technologies, Cergy Pontoise, France) and plated at a density of 1.0 × 105 cells/cm2 on sterilized glass coverslips coated with 60 μg/ml poly D-L ornithine and with medium containing 5% inactivated fetal calf serum (Sigma, St. Louis, MO, USA). To insure easy visualization of a few eGFP-GlyT2 positive axons and their boutons, each coverslip was composed of 5%–10% neurons from GlyT2-eGFP mice; the remaining 90%–95% neurons were from unlabeled WT littermates. Cells were maintained at 37°C in 5% CO2 in serum-free Neurobasal™ medium containing supplement B27 (Invitrogen, Carlsbad, CA, USA; Brewer et al., ) for up to 3 weeks. Medium was changed every 4–5 days.
Electrophysiology
Whole-cell (WC) patch clamp recordings of spinal cord neurons (14–22 DIV) were performed at 30°C. Voltage-clamp was imposed by a Multiclamp 700B amplifier controlled by pCLAMP 9 or 10 acquisition software (Molecular Devices). Currents were filtered at 4 kHz and sampled at 20 kHz using a Digidata 1440A (Molecular Devices). Neurons were continuously bathed with an external solution containing (mM): NaCl 140, KCl 2.4, CaCl2 2, MgCl2 2, Glucose 10, HEPES 10, pH = 7.4. WC patch clamped mIPSCs were recorded in the presence of 0.2 μM of the sodium channel blocker tetrodotoxin (TTX) and 2 μM of the benzodiazepine flunitrazepam. AMPA and NMDA receptors were blocked with 2 μM NBQX (6-nitro-7-sulfamoylbenzo[f]quinoxaline-2, 3-dione) and 5 μM MK-801, respectively. When indicated GABAA receptors (GABAARs) were blocked with 5 μM gabazine (SR 95531), and glycine receptors (GlyRs) were blocked with 1 μM strychnine. The majority of whole cell mIPSCs have fast rise times, and a few cells were found to have a correlation between their glycinergic mIPCS rise and decay time constants, suggesting that dendritic filtering of mIPSCs from more distant synapses influences the kinetic of some, but not all, events included in the WC mIPSC data (Supplementary Figure S2C; Gardner et al., ). There was no evidence of dendritic filtering at any of the single synapse mIPSCs, as expected for events that originated from synaptic terminals that are on, or very close to, the neuronal soma (Supplementary Figure S2C).
Patch-pipettes were pulled from borosilicate glass capillaries (Hilgenberg, Maisfeld, Germany) and pipettes for WC recording had typical resistances of 4–6 MΩ. The mIPSCs were recorded at a holding potential (VH) of −70 mV (taking into account the junction potential) using pipettes filled with a standard internal solution containing (in mM): CsCl 140, CaCl2 1, EGTA 10, BAPTA 1, MgCl2 1, Mg-ATP 4, QX314-Cl 5, Hepes 10, adjusted to pH 7.4 with CsOH. Loose-patch (LP) recording pipettes had resistances of 1–2 MΩ and LP mIPSCs were recorded at 0 mV. One out of every 10 attempts at obtaining a LP recording were successful. LP pipettes were filled with the high calcium external solution (4 mM Ca2+; 0 mM Mg2+), 0.2 μM TTX, 2 μM flunitrazepam, 2 μM NBQX, 5 μM MK-801. Baclofen, a GABABR agonist, was added to the extracellular solution to reduce the frequency of mIPSC when required (O’Brien et al., ). Unless indicated, drugs were purchased from Sigma (St. Louis, MO, USA) or Tocris Bioscience (Bristol, UK).
Data Analysis
Miniature synaptic currents were detected semi-automatically and analyzed using SpAcAn, a custom-made integral detection Igor package developed by Guillaume Dugué and Charly Rousseau and available at www.spacan.net. or in AxoGraph. LP and WC mIPSCs whose peak currents occurred simultaneously (<0.6 ms) and were verified to have a smooth rise and decay phase, as well as a strongly correlated peak amplitudes, were included in the analysis. Electrophysiological results are reported as mean ± SEM. All statistical tests were nonparametric and performed using Prism software unless indicated. The Mann–Whitney and Wilcoxon Matched-Pairs test was used to assess differences between two independent and two related samples respectively. For all tests, the number of asterisks (*) in the figures corresponds to level of significance: ***p < 0.001, **p < 0.01 and *p < 0.05.
Populations of glycine and GABA mIPSCs recorded in the presence of SR 95331 or strychinine, respectively, were classified as glycinergic, GABAergic or mixed on the basis of their peak amplitude to charge transfer relationship and then fitted with a linear regression and 95% prediction intervals. We also attempted to classify mIPSC events by fitting their currents with biexponential and monoexponential fits (Jonas et al., 1993; Rahman et al., ) but found that the peak amplitude to charge relationship robust and fast (Supplementary Figure S2). We set a conservative limit of ≥30% to considering groups of single synapse mIPSCs as a separate population. This value was based on our pharmacological characterization of glycine and GABA mIPSCs.
Numerical Simulation
The set of linear differential equations that defined the GlyR and GABAAR kinetic models was numerically solved in Mathematica 11 using the Q-matrix approach (Colqhoun and Hawkes, ). The current amplitude corresponds to the summation of 120 GlyRs and 60 GABAARs, with elementary currents of 1.05 and 3 pA for the GABAAR mono- and di-liganded receptors, respectively, and 5 pA for GlyR. The rate constants are indicated in Figure 3A, that correspond to the values determined by Burzomato et al. () for GlyR and Labrakakis et al. () for GABAAR except for the dissociation rate constants that were slightly increased to match the experimental time course of the mIPSC and peak-charge amplitude relationship.
Figure 1
Figure 2
Figure 3
We used the simplified assumption of a monoexponential time course for the transient of GABA and glycine (AAN: amino-acid neurotransmitter) with a piecewise function:
The time constant of the transient in Figure 3 was τ = 0.7 ms and the [AAN]peak = [Gly]peak + [GABA]peak = 3 mM (Beato, ). The time offset (t0) was 2 ms and the time of integration was 200 ms.
The simulated mixed mIPSC corresponds to the algebraic sum of the glycine and GABA mIPSCs. We used the FindMinimum and NIntegrate functions of Mathematica 11 (Wolfram Research) to determine the peak amplitude and charge transfer of simulated mIPSCs.
The parameter α set the relative glycine concentration in the transient:
In Figure 4 and Supplementary Figure S3, we examined the distribution of mIPSC phenotypes that were simulated with a broad variation in [Gly]peak and [GABA]peak and/or GABAA and GlyR. Normal distributions of [Gly]peak, [GABA]peak and/or GABAA and GlyR were generated using the Random Variate and NormalDistribution functions of Mathematica with [Gly]peak = N[3α, (3α CV)2] and [GABA]peak = N[3(1−α), (3(1−α)CV)2]. Accordingly, the average [AAN]peak was constant (3 mM) in all simulations.
Figure 4
Ethics Statements
This study was carried out in accordance with the recommendations of guidelines of the Centre National de la Recherche Scientifique. The protocol was approved by the CNRS under number 02235.02 of the general agreement C750520.
Results
Identification of mIPSCs From Individual Terminals Using a Combination of Loose-Patch and Whole-Cell Electrophysiology in Cultured Spinal Neurons From GlyT2-eGFP Mice
We combined the spinal cords of wild-type and GlyT2-eGFP (5–10%) embryos (Zeilhofer et al., ) and made low-density cultures of eGFP-expressing neurons that were used for electrophysiological experiments at DIV14–21. When viewed under fluorescent light, we observed thin, parallel neurites with beaded varicosities (1.05 ± 0.2 μm, n = 30), suggestive of axons (Figure 1A). The GlyT2-eGFP+ beaded varicosities rapidly accumulated the styryl dye FM 4–64, a marker of cycling presynaptic vesicles (Figure 1B), confirming that they correspond to presynaptic terminals.
Next, we voltage-clamped (VH = −70 mV) the cell body of an eGFP− neuron that appeared to have direct contact with one or a few eGFP+ varicosities at or near to its soma, and recorded WC mIPSCs in the presence of TTX (Figure 1C). WC mIPSCs are heterogeneous because they correspond to the postsynaptic neurons response to the spontaneous release of single synaptic vesicle at one of its many inhibitory terminals (Hubbard et al., ). To identify the specific subset of mIPSCs that originated from a single GlyT2-eGFP+ varicosity, we simultaneously recorded extracellular mIPSCs from a second LP recording electrode (VH = 0 mV) placed directly over a single synapse, identified by its eGFP+ varicosity, located at or near to the soma of the recorded neurons (Figure 1C; Forti et al., ). Single-synapse mIPSCs were identified from whole cell mIPSCs (inward) that were time- and amplitude-locked to extracellular LP mIPSCs (outward, Figures 1C,D). LP-mIPSCs were never detected if the LP-pipette was not centered precisely on the GlyT2-eGFP+ varicosity, consistent with the prediction that the LP-pipette will detect less than 1% of currents that originate outside the LP pipette (Forti et al., ) and suggestive that detection of current from two independent synases is unlikely.
GABA- and Glycine-Only mIPSCs Events Were Defined Based on Their Peak Amplitude to Charge Relationship
Although GABA and glycine activate ligand-gated chloride channels with distinct biophysical and pharmacological properties, GABAergic and glycinergic mIPSCs are difficult to separate because of their overlapping receptor kinetics (Burzomato et al., ; Labrakakis et al., ). We used flunitrazepam, a positive allosteric modulator of GABAA receptors, to selectively slow the decay of GABAergic mIPSCs and allow a crude separation of GABA vs. glycine mIPSCs based on their peak amplitude and charge transfer (Figures 2A,B). Nevertheless, the two mIPSC populations were highly variable both within and between neurons (Figure 2C; n = 15 cells, 600 mIPSC events each). The variability is attributable to the presence of multiple independent inhibitory connections with potential differences in: receptor numbers, receptor subunit composition, synaptic architecture, as well as synaptic vesicle content (Ropert et al., ; Frerking et al., ; Auger and Marty, ; Nusser et al., ) and dendritic filtering (Gardner et al., ). This variability was not reduced by analyzing data from neurons prepared from the same animals, or on the same number of days in vitro (data not shown). Linear regressions of the peak amplitude vs. charge data and 95% prediction intervals frame the GABA and glycine mIPSC populations and were used to classify pure and mixed phenotypes (Figure 2C). When we recorded whole cell mIPSCs in the absence of strychnine and SR 95531, and selected mIPSC events with clear biphasic decay currents characteristic of mixed GABA/glycine cotransmission (Jonas et al., ; Lu et al., ), these events mapped in-between the GABA and glycine mIPSC prediction intervals on graphs of peak amplitude vs. charge (Figure 2D) or within the GABA mIPSC prediction intervals (Figure 2E). Thus, the linear regression and 95% prediction intervals of the glycine peak to charge data provides functional limits to define glycinergic and mixed mIPSCs. In contrast, the GABAergic mIPSC data was not able to differentiating between mixed and GABAergic mIPSCs, especially when mixed mIPSC currents had a large GABA component (Figure 2E).
To explore all possible phenotype outcomes in response to variable GABA and glycine vesicular release, we simulate glycinergic, GABAergic and mixed mIPSCs with activation kinetic models for GABAAR and GlyR that have been previously established (Burzomato et al., ; Labrakakis et al., ; Figures 3A,B). We limited our simulations to the simplest, linear substitution of GABA by glycine, keeping a constant [neurotransmitter release] (AANpeak = 3 mM, τr= 0.7 ms (Beato, ), Figure 3C). Therefore, increasing α potentiates the glycine component and decreases the GABA component of mixed mIPSCs (Figures 3D–F) in such a way that the increase in peak amplitude may compensate for the reduction in charge transfer, thus preserving the strength of synaptic inhibition (Figure 3G).
To map a more realistic distribution of inhibitory phenotypes, we plot the peak-charge amplitude distribution of simulated mIPSCs for a randomly generated number of postsynaptic receptors and peak concentration transients (see “Materials and Methods” section). First, we plot the smooth density histograms for 200 simulated mIPSCs when there is a unique, pre- or post-synaptic source of variability with a coefficient of variation of 0.2 (CVR = 0.2 or CVAAN = 0.2, Supplementary Figure S3). Supplementary Figure S3 illustrates that a lower variability of the GABA mIPSCs is expected when α values are low, because GABAAR receptors are almost saturated by 3 mM GABA transient (EC50 = 0.81 mM for GABAAR and 2.68 mM for GlyR, with comparable apparent hill-coefficients of 2.17 and 2.06, respectively). Then we map the distribution of 500 randomly simulated mIPSCs with average α values ranging from 0 (pure GABA) to 1 (pure Glycine) by 0.05 step increments (Figure 4A), with high pre- (CVR = 0.4, Figure 4B) and post- (CVAAN = 0.4, Figure 4C) synaptic variability. By definition, the simulated mIPSCs are mixed and have two components when 0 < α < 1 and would therefore be arbitrarily qualified as mixed, even if one component is small and the mIPSC resembles pure GABAergic or pure glycinergic events. Therefore, we delimit three areas in the charge-peak amplitude plot (Figure 4A) for the attribution of GABAergic (blue area), glycinergic (orange area) and mixed phenotypes (green area), based on the 95% prediction from the linear regression of the experimental values for charge-peak amplitude. The average peak-charge amplitudes for each α values follow the model predictions (Figure 4D) and the repartition of mIPSC phenotypes as function of α (Figure 4E) shows that 10% GABA release (α = 0.9) reduces the glycinergic component of the mixed mIPSCs by 50%, whereas in contrast a similar reduction in GABAergic phenotypes is not achieved until α = 0.3. This finding is consistent with our experimental data (Figures 2D,E) and indicates that mixed mIPSCs with small glycine components are likely to be characterized as a GABAergic phenotype. Overall, the model confirms that neurotransmitter release by vesicles containing higher proportions of glycine than GABA will evoke mixed mIPSCs that have intermediate peak amplitude-charge relationships, falling in-between the GABA and glycine prediction intervals defined by our experimental data. When α = 0.45–0.85 the model predicts that mixed phenotypes will be the dominant phenotype, with <20% of mixed events being mischaracterized.
A Combination of Mixed and Glycinergic mIPSCs Were Detected at Some Individual Terminals
Individual-terminal mIPSCs were recorded in the absence of strychnine and SR 95331, and analyzed when the peak amplitude and charge values were measurable in ≥25 LP-linked mIPSCs. We classified the mIPSC phenotype from 18 individual terminals. In 4 of the 18 synapses only glycine mIPSC were detected (Figure 5A), with a reduced co-efficients of variation (CV) compared to glycine mIPSCs from whole cell recording (0.42 ± 0.5, range 0.30–0.51; and 0.65 ± 0.05, range 0.25–0.86, respectively). The rest (14/18) were cotransmission terminals, with a substantial proportion of mixed mIPSCs. According to our classification system, we found that 36% (5/14) of the cotransmission terminals we sampled had a homogenous mixed mIPSC phenotype as hypothesized (Figure 5B). The majority (64%, 9/14) had a heterogeneous phenotype made up of a combination of mixed and glycine-only mIPSCs (Figures 5C,D). In these combination synapses, the mixed eIPSCs were always larger than glycine-only mIPSCs, as predicted for an additive current, and the mixed current was carried primarily by glycine (Figures 5E,F). Examination of the raw mIPSCs currents grouped into their designated phenotypes (averaged from 10 consecutive events), demonstrates that phenotypically glycinergic mIPSCs have characteristic fast monophasic decay kinetic (Figures 5A,C,D). In contrast, mIPSCs events that fell into the mixed or GABA phenotype regions all had a discernable biphasic kinetics (Figures 5B–D) indicating this population corresponded to mixed events. Together, these data supports the identification of distinct mixed and pure mIPSC phenotypes in combination synapses.
Figure 5
Discussion
Heterogeneity of GABA/Glycine Vesicular Content as a Presynaptic Source of mIPSC Variability
Because identifying the glycine and GABA components in individual mixed IPSCs is not straightforward, we previously used a giant-synapse, made from a secretory BON cell and a sniffer HEK cell, to sample all vesicular phenotypes when a VIAAT-expressing BON cell had access to both neurotransmitters (Aubrey et al., ). In this cellular model, sniffer HEK cells expressed GABA and glycine receptors that had cationic and anionic permeability respectively, thus GABA- and glycine-evoked currents had opposite polarities and could be identified unambiguously (Aubrey et al., , Supplementary Figures S1A–C). Surprisingly, the distribution of pure and mixed events from individual co-releasing BON cells suggested heterologous vesicular GABA/glycine content (Aubrey et al., , Supplementary Figure S1C), as reported in acutely dissociated spinal cord neurons (Katsurabayashi et al., ). We can excluded a significant contribution of the sniffer cell in the BON/HEK model because scanning across the entire sniffer HEK cell surface with brief iontophoretic applications of a fixed [GABA] and [glycine] evoked stable biphasic currents (Supplementary Figures S1D,E), again supporting heterogeneous vesicular content that favors glycine- or GABA-like mIPSCs.
In the present work, we combined LP extracellular recordings with WC recordings to sample mIPSCs from single inhibitory varicosity in spinal cord cultures from GlyT2:EGFP mice. We have shown previously that evoked IPSCs from GlyT2-expressing neurons are predominately mixed, with a larger glycinergic component (Rousseau et al., ). As expected, we detected mixed mIPSCs in the majority (78%) of GlyT2-GFP+ synapses sampled and found that most mixed GABA/glycine synapses signaled with a combination of pure and mixed mIPSCs. This phenotypic heterogeneity may have a presynaptic origin that results from variability in the accumulation and packaging of GABA and glycine in VIAAT-containing synaptic vesicles, similar to the BON/HEK model system (Aubrey et al., ; Supplementary Figure S1), however we cannot rule out a postsynaptic origin because inhibitory synapses are dynamic structures (Choquet and Triller, ), with glycine and GABAA receptors being constantly trapped or exchanged between gephyrin containing nano-domains (Calamai et al., ; Maric et al., ; Tyagarajan and Fritschy, ; Alvarez, ; Pennacchietti et al., ). In addition, varicosities are larger structures than the small en passant boutons of typical central synapses that may contain multiple active zones facing postsynaptic clusters with different receptor compositions (Nusser et al., ; Biró et al., ; Lévi et al., ). When we generated simulated mIPSCs including a pre- and a post-synaptic source of variability (Supplementary Figure S3), we found that both were able to recapitulate the mostly glycinergic mIPSC distributions observed experimentally. We did not attempt to differentiate the source of variability any further in the present work.
mIPSC Detection Bias
The distribution of mIPSC phenotypes, which is based on their peak to charge relationship, indicates that single varicosities predominantly signal with glycine in these cultures (α ≥ 0.7). This dominant glycinergic phenotype was anticipated because: 1. all the synapses sampled had GlyT2-eGFP+ presynaptic varicosities. 2. Glycine mIPSCs have faster kinetics and larger peak amplitude than GABA mIPSCs. As a result, 3. time-integration of mIPSCs was more consistently resolved for glycine events than for GABA events because often GABA mIPSCs that did not return to baseline before a subsequent mIPSC occurred and thus could not be reliably integrated. 4. Mixed mIPSCs with a high or low α-values are likely to be mis-classified, although mixed mIPSCs with small GABA components (α > 0.9) were more easily detected than those with small glycine components (α < 0.3). 5. Finally, mIPSCs with larger peak amplitudes, typically carried by glycine currents in these cultures, were more likely to be detected by the low resistance, LP recordings.
A Model to Study Cotransmission
We modeled mixed mIPSCs as the summation of the GABAARs and GlyRs mediated currents activated by a fast synaptic transient of GABA and glycine. Each randomly generated transient aims to represent the possible content of one VIAAT expressing vesicle. As the time course of glycine and GABA in the synaptic cleft have been thoroughly modeled (Overstreet and Westbrook, ; Beato, ), we used a fixed average peak concentration of 3 mM and a synaptic transient time-constant of 0.7 ms for the simulation. We fixed the peak-concentration transient because the amount of neurotransmitter stored in a vesicle may be limited by VIAAT thermodynamics, VIAAT kinetics or set by other vesicular factors (Edwards, ).
We used kinetic models and rate constants well-established for both receptors (Burzomato et al., ; Labrakakis et al., ), and assumed for simplicity that glycine and GABA act independently, although it has been shown that GABA/glycine corelease accelerates the kinetics of glycinergic IPSCs in pure glycinergic synapses (Lu et al., ). The distribution of randomly generated simulated mIPSCs on the charge-peak amplitude plot provided a graphical overview of all transitions between pure GABAergic and pure glycinergic phenotypes when the proportion of glycine in the neurotransmitter transient is sequentially increased. We operationally distinguished pure and mixed phenotypes for individual mIPSCs based on their location on three experimentally determined areas in the charge-peak amplitude plot. The distribution of “pure” GABAergic, mixed and “pure” glycinergic phenotypes shows that “pure” glycinergic phenotype is only compatible with pure glycine release since α = 0.95 (about 150 μM GABA and 2.85 mM Glycine) is enough to evoked mixed phenotypes. In contrast “pure” GABAergic phenotypes are reported up to α = 0.2 (about 2.4 mM GABA and 0.6 mM glycine), indicating that minor glycine corelease may more easily go unnoticed than minor GABA corelease. This asymmetry reflects the apparent saturation of GABAARs with a 3 mM GABA transient.
Functional Roles for GABA/Glycine Corelease
GABA/glycine co-signaling is habitually thought to be a property of immature inhibitory synapses in the spinal cord, brainstem and granular layer of the cerebellum that diminishes during postnatal development, with few exceptions (Inquimbert et al., ; Bhumbra et al., ). Indeed mixed GABA/glycine cotransmission was initially detected in juvenile rat spinal and brainstem motoneurons (Jonas et al., ; O’Brien and Berger, ; Russier et al., ) and lamina I and II of the spinal cord dorsal horn (Keller et al., ). Furthermore, a critical shift toward pure glycinergic transmission takes place during the postnatal maturation in spinal cord locomotor centers (Gao et al., ) and brainstem auditory nuclei where the large amplitude and rapid decay of glycinergic IPSCs provide the microsecond precision needed for the processing and the coding of acoustic information (Brand et al., ; Nabekura et al., ).
However, corelease with or without cotransmission persists in adults at some synapses, suggesting a continuing specific need for each neurotransmitter. In the superficial lamina of the spinal cord, GABA and glycine cotransmission is replaced by pure glycinergic transmission after P23 (Keller et al., ). Vesicular corelease persists in into adulthood however, and coreleased GABA influences the release probability by acting at presynaptic GABABR (Chéry and de Koninck, 2000). Furthermore, extrasynaptic GABAA receptor signals can be detected at higher stimulation intensities (Chéry and de Koninck, ). At physiological temperatures evidence suggests that GABA and glycine corelease may not contribute to spillover and tonic inhibition in the dorsal horn of the spinal cord, as the diffusion of GABA and glycine out of the synaptic cleft is readily controlled by transporters in this region (Mitchell et al., ). Nevertheless, spillover of synaptically released glycine has been shown to potentiate the NMDA responses in the superficial dorsal horn of the spinal cord, by increasing the occupancy of the NMDA receptor coagonist site (Ahmadi et al., ). The same developmental pattern of inhibitory transmitter use has been reported in the rat MNTB, where GABA/glycine cotransmission is detected until P9–P12 after which time pure glycinergic transmission dominates (Awatramani et al., ). Again vesicular corelease persists in older animals, and coreleased GABA has been shown to shorten glycinergic mIPSC decay times here (Lu et al., ). In the central nucleus of the inferior colliculus of the auditory midbrain, GABA and glycine cotransmission persists in brain slices from P19 to P35 mice. Coreleased GABA and glycine evoke IPSCs with nearly identical amplitude and time course, suggesting that GABA and glycine are operationally fully interchangeable in this nucleus, thus securing inhibitory signaling by redundancy (Moore and Trussell, ). Finally, in the cerebellum of juvenile rats, inhibitory corelease occurs between Golgi cells and vestibulocerebellar Unipolar Brush Cells (Dugué et al., ; Rousseau et al., ), and once again evidence suggests that the transmission phenotype is under postsynaptic control. Interestingly, the phenotype of transmission in this region is not fixed developmentally, but seems to be coordinated with the phenotype of the glutamatergic input onto the Unipolar Brush Cell (Rousseau et al., ).
Our data suggests that a dynamic vesicular GABA/glycine balance at individual synapses may contribute to the tuning of phasic synaptic inhibition in spinal neurons. If the origin of the variability is indeed the neurotransmitter content in presynaptic vesicles, then glycine IPSC kinetics (Lu et al., ) and the probability of transmitter corelease (where GABAB receptors are expressed (Chéry and de Koninck, , 2000) would also contribute to this fine tuning.
Possible Mechanisms for Vesicular GABA/Glycine Variations at Single Terminals
The mechanisms that specify or regulate the vesicular GABA/glycine content at mixed synapses are not well understood. As previously discussed, changes to GABA (Mathews and Diamond, ; Wang et al., ) and glycine supply/resupply (Rousseau et al., ; Apostolides and Trussell, ; Ishibashi et al., ) can strongly shift the GABA/glycine balance and alterations in IPSC characteristics can be observed within minutes. Indeed, metabolic alteration such as these have been implicated in disease states including chronic pain (Coull et al., ; Imlach et al., 2016; Takazawa et al., 2017) and amyotrophic lateral sclerosis (Medelin et al., ).
With an apparent lower affinity for glycine than for GABA (about 25 and 6 mM, respectively), VIAAT uptake is likely to be rate limiting and a major source of variability for vesicular loading (Burger et al., ; Gasnier, ; Edwards, ; Farsi et al., ; Takamori, ). In addition, VIAAT has presumably the lowest driving force for uphill transport among the vesicular transporters (Edwards, ), being coupled to the exchange of a single H+ (Hell et al., ; McIntire et al., ; Farsi et al., ) and with smaller protonmotive force (Egashira et al., ). If VIAAT kinetic is limiting for filling vesicles, as suggested by a slow recovery from synaptic depression (Katsurabayashi et al., ; Rousseau et al., ; Wang et al., ; Yamashita et al., ), then the initial vesicular content may continue to change as VIAAT exchanges a fast filling neurotransmitter (presumably GABA) with cytosolic glycine, providing opportunity for alterations in the transient neurotransmitter composition.
Other potential mechanisms include protein-protein interactions between the transmitter supply machinery and different populations of vesicles (Jin et al., ), differential modulation of the pool of cycling vesicle in response to presynaptic GABA or glycine supply (Mathews and Diamond, ; Wang et al., ; Truckenbrodt et al., ) or differential modes of GABA and glycine uptake by VIAAT (Aubrey et al., ). Given the diverse and highly regulated roles of the presynaptic terminal, this list is in not exhaustive.
In summary, this work adds to a growing body of evidence that suggest that the presynaptic terminal and neurotransmitter corelease can significantly regulate synaptic transmission. We present evidence that the GABA/glycine balance at individual co-releasing terminals is variable, and predict these variations would serve to fine tune the timing of inhibition and the integration of sensory inputs. Given the complexity and high level of control the presynapse has over synaptic vesicle release (Rizzoli and Betz, ; Edwards, ; Südhof, ), small, dynamic alternations in the balance of GABA/glycine signaling are possible and likely to contribute to synaptic plasticity.
Cotransmission is known to shape synaptic plasticity of other brain regions (Vaaga et al., ; Tritsch et al., ). For example, in the lateral habenula the GABA component of GABA/glutamate co-releasing inputs from the basal ganglia are selectively reduced in an animal model of depression. This deficiency is restored following sustained treatment with the antidepressant citalopram (Shabel et al., ). A similar reduction in the GABA component of mixed GABA/glutamate release has been observed in the same region in animal models of cocaine withdrawal, where it is attributed to a selective decrease in VIAAT expression in vesicles (Meye et al., ).
Statements
Author contributions
SS and KA conceived and planned the experiments. KA carried out the experiments and their analysis. SS carried out the modeling and simulations. SS and KA contributed to the interpretation of the results and wrote the manuscript.
Funding
This project received support from the Centre National de la Recherche Scientifique (CNRS), the ENS, the Institut National de la Santé et de la Recherche Médicale (INSERM), the program «Investissements d’Avenir» (ANR-10-LABX-54 MEMO LIFE, ANR-11-IDEX-0001-02 PSL), the FRM (DEQ 20140329498) and grants from the Association Française contre les Myopathies. KA was supported by a fellowship from the Association Française contre les Myopathies and received support from PMRI Ltd.
Acknowledgments
We thank Hanns U. Zeilhofer for providing the GlyT2-EGFP mice and Bruno Gasnier for providing BON:VIAAT cells. We thank Guillaume Dugué and Charly Rousseau for SpAcAn (http://www.spacan.net/) adaptations; and Dan Costa and Rachel O’Connell for statistical advice.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnsyn.2018.00040/full#supplementary-material
References
1
AhmadiS.Muth-SelbachU.LauterbachA.LipfertP.NeuhuberW. L.ZeilhoferH. U. (2003). Facilitation of spinal NMDA receptor currents by spillover of synaptically released glycine. Science300, 2094–2097. 10.1126/science.1083970
2
AlvarezF. J. (2017). Gephyrin and the regulation of synaptic strength and dynamics at glycinergic inhibitory synapses. Brain Res. Bull.129, 50–65. 10.1016/j.brainresbull.2016.09.003
3
ApostolidesP. F.TrussellL. O. (2013). Rapid, activity-independent turnover of vesicular transmitter content at a mixed glycine/GABA synapse. J. Neurosci.33, 4768–4781. 10.1523/JNEUROSCI.5555-12.2013
4
AubreyK. R.RossiF. M.RuivoR.AlboniS.BellenchiG.Le GoffA.et al. (2007). The transporters GlyT2 and VIAAT cooperate to determine the vesicular glycinergic phenotype. J. Neurosci.27, 6273–6281. 10.1523/JNEUROSCI.1024-07.2007
5
AugerC.MartyA. (1997). Heterogeneity of functional synaptic parameters among single release sites. Neuron19, 139–150. 10.1016/s0896-6273(00)80354-2
6
AwatramaniG. B.TurecekR.TrussellL. O. (2005). Staggered development of GABAergic and glycinergic transmission in the MNTB. J. Neurophysiol.93, 819–828. 10.1152/jn.00798.2004
7
BeatoM. (2008). The time course of transmitter at glycinergic synapses onto motoneurons. J. Neurosci.28, 7412–7425. 10.1523/JNEUROSCI.0581-08.2008
8
BhumbraG. S.MooreN. J.MoroniM.BeatoM. (2012). Co-release of GABA does not occur at glycinergic synapses onto lumbar motoneurons in juvenile mice. Front. Cell. Neurosci.6:8. 10.3389/fncel.2012.00008
9
BiróA. A.HolderithN. B.NusserZ. (2006). Release probability-dependent scaling of the postsynaptic responses at single hippocampal GABAergic synapses. J. Neurosci.26, 12487–12496. 10.1523/JNEUROSCI.3106-06.2006
10
BrandA.BehrendO.MarquardtT.McAlpineD.GrotheB. (2002). Precise inhibition is essential for microsecond interaural time difference coding. Nature417, 543–547. 10.1038/417543a
11
BrewerG. J.TorricelliJ. R.EvegeE. K.PriceP. J. (1993). Optimized survival of hippocampal neurons in B27-supplemented neurobasal™, a new serum-free medium combination. J. Neurosci. Res.35, 567–576. 10.1002/jnr.490350513
12
BurgerP. M.HellJ.MehlE.KraselC.LottspeichF.JahnR. (1991). GABA and glycine in synaptic vesicles: storage and transport characteristics. Neuron7, 287–293. 10.1016/0896-6273(91)90267-4
13
BurzomatoV.BeatoM.Groot-KormelinkP. J.ColquhounD.SivilottiL. G. (2004). Single-channel behavior of heteromeric α1β glycine receptors: an attempt to detect a conformational change before the channel opens. J. Neurosci.24, 10924–10940. 10.1523/JNEUROSCI.3424-04.2004
14
CalamaiM.SpechtC. G.HellerJ.AlcorD.MachadoP.VannierC.et al. (2009). Gephyrin oligomerization controls GlyR mobility and synaptic clustering. J. Neurosci.29, 7639–7648. 10.1523/JNEUROSCI.5711-08.2009
15
ChéryN.de KoninckY. (1999). Junctional versus extrajunctional glycine and GABAA receptor- mediated IPSCs in identified lamina I neurons of the adult rat spinal cord. J. Neurosci.19, 7342–7355. 10.1523/JNEUROSCI.19-17-07342.1999
16
ChéryN.de KoninckY. (2000). GABAB receptors are the first target of released GABA at lamina I inhibitory synapses in the adult rat spinal cord. J. Neurophysiol.84, 1006–1011. 10.1152/jn.2000.84.2.1006
17
ChoquetD.TrillerA. (2013). The dynamic synapse. Neuron80, 691–703. 10.1016/j.neuron.2013.10.013
18
ChristensenH.FonnumF. (1991). The ontogeny of the uptake systems for glycine, GABA and glutamate in synaptic vesicles isolated from rat spinal cord-medulla. Dev. Brain Res.64, 155–159. 10.1016/0165-3806(91)90219-9
19
ColqhounD.HawkesA. (1995). A Q-Matric Cookbook.New York, NY: Plenum Press.
20
CoullJ. A.BoudreauD.BachandK.PrescottS. A.NaultF.SíkA.et al. (2003). Trans-synaptic shift in anion gradient in spinal lamina I neurons as a mechanism of neuropathic pain. Nature424, 938–942. 10.1038/nature01868
21
DuguéG. P.DumoulinA.TrillerA.DieudonnéS. (2005). Target-dependent use of co-released inhibitory transmitters at central synapses. J. Neurosci.25, 6490–6498. 10.1523/JNEUROSCI.1500-05.2005
22
EdwardsR. H. (2007). The neurotransmitter cycle and quantal size. Neuron55, 835–858. 10.1016/j.neuron.2007.09.001
23
EgashiraY.TakaseM.WatanabeS.IshidaJ.FukamizuA.KanekoR.et al. (2016). Unique pH dynamics in GABAergic synaptic vesicles illuminates the mechanism and kinetics of GABA loading. Proc. Natl. Acad. Sci. U S A113, 10702–10707. 10.1073/pnas.1604527113
24
FarsiZ.PreobraschenskiJ.van den BogaartG.RiedelD.JahnR.WoehlerA. (2016). Single-vesicle imaging reveals different transport mechanisms between glutamatergic and GABAergic vesicles. Science351, 981–984. 10.1126/science.aad8142
25
FischlM. J.BurgerM. R. (2014). Glycinergic transmission modulates GABAergic inhibition in the avian auditory pathway. Front. Neural Circuits8:19. 10.3389/fncir.2014.00019
26
FortiL.BossiM.BergamaschiA.VillaA.MalgaroliA. (1997). Loose-patch recordings of single quanta at individual hippocampal synapses. Nature388, 874–878. 10.1038/42251
27
FrerkingM.BorgesS.WilsonM. (1995). Variation in GABA mini amplitude is the consequence of variation in transmitter concentration. Neuron15, 885–895. 10.1016/0896-6273(95)90179-5
28
FykseE. M.FonnumF. (1988). Uptake of gamma-aminobutyric acid by a synaptic vesicle fraction isolated from rat brain. J. Neurochem.50, 1237–1242. 10.1111/j.1471-4159.1988.tb10599.x
29
GaoB.-X.StrickerC.Ziskind-ConhaimL. (2001). Transition from GABAergic to glycinergic synaptic transmission in newly formed spinal networks. J. Neurophysiol.86, 492–502. 10.1152/jn.2001.86.1.492
30
GardnerS. M.TrussellL. O.OertelD. (1999). Time course and permeation of synaptic AMPA receptors in cochlear nuclear neurons correlate with input. J. Neurosci.19, 8721–8729. 10.1523/JNEUROSCI.19-20-08721.1999
31
GasnierB. (2000). The loading of neurotransmitters into synaptic vesicles. Biochimie82, 327–337. 10.1016/s0300-9084(00)00221-2
32
GomezaJ.OhnoK.HülsmannS.ArmsenW.EulenburgV.RichterD. W.et al. (2003). Deletion of the mouse glycine transporter 2 results in a hyperekplexia phenotype and postnatal lethality. Neuron40, 797–806. 10.1016/s0896-6273(03)00673-1
33
GrangerA. J.WallaceM. L.SabatiniB. L. (2017). Multi-transmitter neurons in the mammalian central nervous system. Curr. Opin. Neurobiol.45, 85–91. 10.1016/j.conb.2017.04.007
34
HanusC.VannierC.TrillerA. (2004). Intracellular association of glycine receptor with gephyrin increases its plasma membrane accumulation rate. J. Neurosci.24, 1119–1128. 10.1523/JNEUROSCI.4380-03.2004
35
HellJ. W.EdelmannL.HartingerJ.JahnR. (1991). Functional reconstitution of the gamma-aminobutyric acid transporter from synaptic vesicles using artificial ion gradients. Biochemistry30, 11795–11800. 10.1021/bi00115a009
36
HnaskoT. S.EdwardsR. H. (2012). Neurotransmitter corelease: mechanism and physiological role. Annu. Rev. Physiol.74, 225–243. 10.1146/annurev-physiol-020911-153315
37
HubbardJ. I.StenhouseD.EcclesR. M. (1967). Origin of synaptic noise. Science157, 330–331. 10.1126/science.157.3786.330
38
InquimbertP.RodeauJ. L.SchlichterR. (2007). Differential contribution of GABAergic and glycinergic components to inhibitory synaptic transmission in lamina II and laminae III-IV of the young rat spinal cord. Eur. J. Neurosci.26, 2940–2949. 10.1111/j.1460-9568.2007.05919.x
39
ImlachW. L.BholaR. F.MohammadiS. A.ChristieM.-J. (2016). Glycinergic dysfunction in a subpopulation of dorsal horn interneurons in a rat model of neuropathic pain. Sci. Rep.6:e37104. 10.1038/srep37104
40
IshibashiH.YamaguchiJ.NakahataY.NabekuraJ. (2013). Dynamic regulation of glycine-GABA co-transmission at spinal inhibitory synapses by neuronal glutamate transporter. J. Physiol.591, 3821–3832. 10.1113/jphysiol.2012.250647
41
JinH.WuH.OsterhausG.WeiJ.DavisK.ShaD.et al. (2003). Demonstration of functional coupling between gamma–aminobutyric acid (GABA) synthesis and vesicular GABA transport into synaptic vesicles. Proc. Natl. Acad. Sci. U S A100, 4293–4298. 10.1073/pnas.0730698100
42
JonasP.MajorG.SakmannB. (1993). Quantal components of unitary EPSCs at the mossy fibre synapse on CA3 pyramidal cells of rat hippocampus. J. Physiol.472, 615–663.
43
JonasP.BischofbergerJ.SandkühlerJ. (1998). Corelease of two fast neurotransmitters at a central synapse. Science281, 419–424. 10.1126/science.281.5375.419
44
JonesM. V.WestbrookG. L. (1995). Desensitized states prolong GABAA channel responses to brief agonist pulses. Neuron15, 181–191. 10.1016/0896-6273(95)90075-6
45
KakizakiT.OriuchiN.YanagawaY. (2015). GAD65/GAD67 double knockout mice exhibit intermediate severity in both cleft palate and omphalocele compared with GAD67 knockout and VGAT knockout mice. Neuroscience288, 86–93. 10.1016/j.neuroscience.2014.12.030
46
KatsurabayashiS.KubotaH.HigashiH.AkaikeN.ItoY. (2004). Distinct profiles of refilling of inhibitory neurotransmitters into presynaptic terminals projecting to spinal neurones in immature rats. J. Physiol.560, 469–478. 10.1113/jphysiol.2004.067017
47
KellerA. F.CoullJ. A.CheryN.PoisbeauP.De KoninckY. (2001). Region-specific developmental specialization of GABA-glycine cosynapses in laminas I-II of the rat spinal dorsal horn. J. Neurosci.21, 7871–7880. 10.1523/jneurosci.21-20-07871.2001
48
KotakV. C.KoradaS.SchwartzI. R.SanesD. H. (1998). A developmental shift from GABAergic to glycinergic transmission in the central auditory system. J. Neurosci.18, 4646–4655. 10.1523/jneurosci.18-12-04646.1998
49
LabrakakisC.RudolphU.De KoninckY. (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. 10.3389/fncel.2014.00424
50
LéviS.SchweizerC.BannaiH.PascualO.CharrierC.TrillerA. (2008). Homeostatic regulation of synaptic glyr numbers driven by lateral diffusion. Neuron59, 261–273. 10.1016/j.neuron.2008.05.030
51
LuT.RubioM. E.TrussellL. O. (2008). Glycinergic transmission shaped by the corelease of GABA in a mammalian auditory synapse. Neuron57, 524–535. 10.1016/j.neuron.2007.12.010
52
MaricH. M.KasaragodV. B.HausratT. J.KneusselM.TretterV.StrømgaardK.et al. (2014). Molecular basis of the alternative recruitment of GABA(A) versus glycine receptors through gephyrin. Nat. Commun.5:5767. 10.1038/ncomms6767
53
MartinD. L.RimvallK. (1993). Regulation of gamma-aminobutyric acid synthesis in the brain. J. Neurochem.60, 395–407. 10.1111/j.1471-4159.1993.tb03165.x
54
MathewsG. C.DiamondJ. S. (2003). Neuronal glutamate uptake contributes to GABA synthesis and inhibitory synaptic strength. J. Neurosci.23, 2040–2048. 10.1523/jneurosci.23-06-02040.2003
55
McIntireS. L.ReimerR. J.SchuskeK.EdwardsR. H.JorgensenE. M. (1997). Identification and characterization of the vesicular GABA transporter. Nature389, 870–876. 10.1038/39908
56
MedelinM.RancicV.CellotG.LaishramJ.VeeraraghavanP.RossiC.et al. (2016). Altered development in GABA co-release shapes glycinergic synaptic currents in cultured spinal slices of the SOD1(G93A) mouse model of amyotrophic lateral sclerosis. J. Physiol.594, 3827–3840. 10.1113/jp272382
57
MeyeF. J.Soiza-ReillyM.SmitT.DianaM. A.SchwarzM. K.MameliM. (2016). Shifted pallidal co-release of GABA and glutamate in habenula drives cocaine withdrawal and relapse. Nat. Neurosci.19, 1019–1024. 10.1038/nn.4334
58
MitchellE. A.GentetL. J.DempsterJ.BelelliD. (2007). GABAA and glycine receptor-mediated transmission in rat lamina II neurones: relevance to the analgesic actions of neuroactive steroids. J. Physiol.583, 1021–1040. 10.1113/jphysiol.2007.134445
59
MooreL. A.TrussellL. O. (2017). Corelease of inhibitory neurotransmitters in the mouse auditory midbrain. J. Neurosci.37, 9453–9464. 10.1523/jneurosci.1125-17.2017
60
NabekuraJ.KatsurabayashiS.KakazuY.ShibataS.MatsubaraA.JinnoS.et al. (2004). Developmental switch from GABA to glycine release in single central synaptic terminals. Nat. Neurosci.7, 17–23. 10.1038/nn1170
61
NerlichJ.KeineC.RübsamenR.BurgerR. M.MilenkovicI. (2014). Activity-dependent modulation of inhibitory synaptic kinetics in the cochlear nucleus. Front. Neural Circuits8:145. 10.3389/fncir.2014.00145
62
NusserZ.Cull-CandyS.FarrantM. (1997). Differences in synaptic GABA(A) receptor number underlie variation in GABA mini amplitude. Neuron19, 697–709. 10.1016/s0896-6273(00)80382-7
63
O’BrienJ. A.BergerA. J. (1999). Cotransmission of GABA and glycine to brain stem motoneurons. J. Neurophysiol.82, 1638–1641. 10.1152/jn.1999.82.3.1638
64
O’BrienJ. A.SebeJ. Y.BergerA. J. (2004). GABAB modulation of GABAA and glycine receptor-mediated synaptic currents in hypoglossal motoneurons. Respir. Physiol. Neurobiol.141, 35–45. 10.1016/j.resp.2004.03.009
65
OverstreetL. S.WestbrookG. L. (2003). Synapse density regulates independence at unitary inhibitory synapses. J. Neurosci.23, 2618–2626. 10.1523/jneurosci.23-07-02618.2003
66
PennacchiettiF.VasconS.NieusT.RosilloC.DasS.TyagarajanS. K.et al. (2017). Nanoscale molecular reorganization of the inhibitory postsynaptic density is a determinant of GABAergic synaptic potentiation. J. Neurosci.37, 1747–1756. 10.1523/jneurosci.0514-16.2016
67
RahmanJ.LatalA. T.BesserS.HirrlingerJ.HülsmannS. (2013). Mixed miniature postsynaptic currents resulting from co-release of glycine and GABA recorded from glycinergic neurons in the neonatal respiratory network. Eur. J. Neurosci.37, 1229–1241. 10.1111/ejn.12136
68
RizzoliS. O.BetzW. J. (2004). The structural organization of the readily releasable pool of synaptic vesicles. Science303, 2037–2039. 10.1126/science.1094682
69
RopertN.MilesR.KornH. (1990). Characteristics of miniature inhibitory postsynaptic currents in CA1 pyramidal neurones of rat hippocampus. J. Physiol.428, 707–722. 10.1113/jphysiol.1990.sp018236
70
RousseauF.AubreyK. R.SupplissonS. (2008). The glycine transporter GlyT2 controls the dynamics of synaptic vesicle refilling in inhibitory spinal cord neurons. J. Neurosci.28, 9755–9768. 10.1523/jneurosci.0509-08.2008
71
RousseauC. V.DuguéG. P.DumoulinA.MugnainiE.DieudonnéS.DianaM. A. (2012). Mixed inhibitory synaptic balance correlates with glutamatergic synaptic phenotype in cerebellar unipolar brush cells. J. Neurosci.32, 4632–4644. 10.1523/JNEUROSCI.5122-11.2012
72
RouxM. J.SupplissonS. (2000). Neuronal and glial glycine transporters have different stoichiometries. Neuron25, 373–383. 10.1016/s0896-6273(00)80901-0
73
RussierM.KopysovaI. L.AnkriN.FerrandN.DebanneD. (2002). GABA and glycine co-release optimizes functional inhibition in rat brainstem motoneurons in vitro. J. Physiol.541, 123–137. 10.1113/jphysiol.2001.016063
74
SagnéC.El MestikawyS.IsambertM. F.HamonM.HenryJ. P.GirosB.et al. (1997). Cloning of a functional vesicular GABA and glycine transporter by screening of genome databases. FEBS Lett.417, 177–183. 10.1016/s0014-5793(97)01279-9
75
ShabelS. J.ProulxC. D.PirizJ.MalinowR. (2014). Mood regulation. GABA/glutamate co-release controls habenula output and is modified by antidepressant treatment. Science345, 1494–1498. 10.1126/science.1250469
76
StrataP.HarveyR. (1999). Dale’s principle. Brain Res. Bull.50, 349–350. 10.1016/S0361-9230(99)00100-8
77
SüdhofT. C. (2012). The presynaptic active zone. Neuron75, 11–25. 10.1016/j.neuron.2012.06.012
78
TakamoriS. (2016). Presynaptic molecular determinants of quantal size. Front. Synaptic Neurosci.8:2. 10.3389/fnsyn.2016.00002
79
TakazawaT.ChoudhuryP.TongC.-K. K.ConwayC. M.ScherrerG.FloodP. D.et al. (2017). Inhibition mediated by glycinergic and GABAergic receptors on excitatory neurons in mouse superficial dorsal horn is location-specific but modified by inflammation. J. Neurosci.37, 2336–2348. 10.1523/jneurosci.2354-16.2017
80
TritschN. X.GrangerA. J.SabatiniB. L. (2016). Mechanisms and functions of GABA co-release. Nat. Rev. Neurosci.17, 139–145. 10.1038/nrn.2015.21
81
TruckenbrodtS.ViplavA.JähneS.VogtsA.DenkerA.WildhagenH.et al. (2018). Newly produced synaptic vesicle proteins are preferentially used in synaptic transmission. EMBO J.37:e98044. 10.15252/embj.201798044
82
TyagarajanS. K.FritschyJ. M. (2014). Gephyrin: a master regulator of neuronal function?Nat. Rev. Neurosci.15, 141–156. 10.1038/nrn3670
83
VaagaC. E.BorisovskaM.WestbrookG. L. (2014). Dual-transmitter neurons: functional implications of co-release and co-transmission. Curr. Opin. Neurobiol.29, 25–32. 10.1016/j.conb.2014.04.010
84
WangL.TuP.BonetL.AubreyK. R.SupplissonS. (2013). Cytosolic transmitter concentration regulates vesicle cycling at hippocampal GABAergic terminals. Neuron80, 143–158. 10.1016/j.neuron.2013.07.021
85
WojcikS. M.KatsurabayashiS.GuilleminI.FriaufE.RosenmundC.BroseN.et al. (2006). A shared vesicular carrier allows synaptic corelease of GABA and glycine. Neuron50, 575–587. 10.1016/j.neuron.2006.04.016
86
YamashitaM.KawaguchiS. Y.HoriT.TakahashiT. (2018). Vesicular GABA uptake can be rate limiting for recovery of ipscs from synaptic depression. Cell Rep.22, 3134–3141. 10.1016/j.celrep.2018.02.080
87
ZeilhoferH. U.StudlerB.ArabadziszD.SchweizerC.AhmadiS.LayhB.et al. (2005). Glycinergic neurons expressing enhanced green fluorescent protein in bacterial artificial chromosome transgenic mice. J. Comp. Neurol.482, 123–141. 10.1002/cne.20349
Summary
Keywords
GABA, glycine, inhibitory neurotransmission, VIAAT, VGAT, GlyT2-eGFP mouse, cotransmission, quantal release
Citation
Aubrey KR and Supplisson S (2018) Heterogeneous Signaling at GABA and Glycine Co-releasing Terminals. Front. Synaptic Neurosci. 10:40. doi: 10.3389/fnsyn.2018.00040
Received
02 August 2018
Accepted
12 October 2018
Published
06 November 2018
Volume
10 - 2018
Edited by
Alberto A. Rasia-Filho, Federal University of Health Sciences of Porto Alegre, Brazil
Reviewed by
Dominique Debanne, INSERM U1072 Neurobiologie des canaux Ioniques et de la Synapse, France; Werner Kilb, Johannes Gutenberg-Universität Mainz, Germany
Updates
Copyright
© 2018 Aubrey and Supplisson.
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) and the copyright owner(s) 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: Karin R. Aubrey karin.aubrey@sydney.edu.auorcid.org/0000-0002-1808-4041 Stéphane Supplisson stephane.supplisson@ens.frorcid.org/0000-0002-0062-9752
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.