Identification of a stereotypic molecular arrangement of endogenous glycine receptors at spinal cord synapses

  1. Stephanie A Maynard
  2. Philippe Rostaing
  3. Natascha Schaefer
  4. Olivier Gemin
  5. Adrien Candat
  6. Andréa Dumoulin
  7. Carmen Villmann
  8. Antoine Triller  Is a corresponding author
  9. Christian G Specht  Is a corresponding author
  1. Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, France
  2. Institute for Clinical Neurobiology, University Hospital, Julius-Maximilians-University, Germany
  3. Diseases and Hormones of the Nervous System (DHNS), Inserm U1195, Université Paris-Saclay, France

Abstract

Precise quantitative information about the molecular architecture of synapses is essential to understanding the functional specificity and downstream signaling processes at specific populations of synapses. Glycine receptors (GlyRs) are the primary fast inhibitory neurotransmitter receptors in the spinal cord and brainstem. These inhibitory glycinergic networks crucially regulate motor and sensory processes. Thus far, the nanoscale organization of GlyRs underlying the different network specificities has not been defined. Here, we have quantitatively characterized the molecular arrangement and ultra-structure of glycinergic synapses in spinal cord tissue using quantitative super-resolution correlative light and electron microscopy. We show that endogenous GlyRs exhibit equal receptor-scaffold occupancy and constant packing densities of about 2000 GlyRs µm-2 at synapses across the spinal cord and throughout adulthood, even though ventral horn synapses have twice the total copy numbers, larger postsynaptic domains, and more convoluted morphologies than dorsal horn synapses. We demonstrate that this stereotypic molecular arrangement is maintained at glycinergic synapses in the oscillator mouse model of the neuromotor disease hyperekplexia despite a decrease in synapse size, indicating that the molecular organization of GlyRs is preserved in this hypomorph. We thus conclude that the morphology and size of inhibitory postsynaptic specializations rather than differences in GlyR packing determine the postsynaptic strength of glycinergic neurotransmission in motor and sensory spinal cord networks.

Editor's evaluation

The manuscript presents a quantitative and comprehensive characterization of glycinergic synapses and spinal cord synapse organization that will be a highly valuable resource for the field. The use of cutting-edge imaging and quantitative techniques allows an understanding of these understudied synapses at a high level.

https://doi.org/10.7554/eLife.74441.sa0

Introduction

Synaptic transmission relies on the integration of spatially and temporally controlled signals by neurotransmitter receptors in the postsynaptic membrane. The molecular arrangement of postsynaptic receptors and scaffold proteins is therefore key to the synaptic function. However, the heterogeneity and complexity of postsynaptic sites have made it difficult to resolve its internal organization, ascertain whether distinct compositional states exist, and determine how the organization is affected in disease.

Glycine receptors (GlyRs) are the main inhibitory neurotransmitter receptors in the adult spinal cord and brainstem. Glycinergic neurons arise from different embryonic origins, with specific types of neurons residing in characteristic layers of the spinal cord (Lu et al., 2015). Depending on their location, glycinergic neurons mediate sensory and motor information in the dorsal and ventral spinal cord, respectively, which requires high reliability and fidelity of transmission (Alvarez, 2017). Consequently, deficits in glycinergic transmission are involved in pain mechanisms (Harvey et al., 2004) and motor-related neurological diseases (Schaefer et al., 2018). The electrophysiological properties of glycinergic currents indicate that only a limited number of receptors are activated by the release of a single synaptic vesicle (Oleskevich et al., 1999; Singer and Berger, 1999), suggesting that the nanoscale organization of the receptors determines signal amplitude.

It has been shown that neurotransmitter receptors at excitatory and inhibitory synapses are organized within sub-synaptic domains (SSDs) that are aligned with presynaptic elements of the active zone (AZ) (Crosby et al., 2019; MacGillavry et al., 2013; Pennacchietti et al., 2017; Tang et al., 2016; Yang et al., 2021). These so-called trans-synaptic nanocolumns are thought to increase the efficacy of synaptic transmission (Haas et al., 2018). At mixed inhibitory synapses, both glycine and GABAA receptors are immobilized opposite to presynaptic release sites through direct interactions with their common scaffold protein gephyrin (Maric et al., 2011; Specht et al., 2013; Yang et al., 2021). However, accurate quantification of receptor numbers and their precise arrangement within postsynaptic sites in native tissue is lacking. Further, the question is raised as to whether the structure of glycinergic synapses varies in functionally diverse circuits of the dorsal and ventral spinal cord, whether it changes over time, and whether it is disturbed in GlyR pathologies such as the neuromotor disease hyperekplexia in humans. In hyperekplexic patients, mutations in the receptor subunit genes GLRA1 and GLRB lead to decreased receptor availability and disturbances in glycinergic transmission, resulting in exaggerated startle reflexes, muscle hypertonia, and stiffness in infancy (Chung et al., 2013; Chung et al., 2010; Schaefer et al., 2013).

To investigate whether the molecular arrangement of GlyRs may account for differences in the functional specificity of sensory and motor circuits, we have quantitatively analyzed the ultra-structural organization of inhibitory synapses in spinal cord tissue. We have combined molecule counting of endogenous GlyRs using single-molecule localization microscopy (SMLM) with correlative light and electron microscopic analysis (CLEM) to obtain receptor numbers as well as detailed spatial information of the synapse at the nanometer scale. We have further examined to what extent the molecular organization is maintained throughout adult development and during GlyRα1 deficiency. We show that GlyRs are packed at a constant density of about 2000 receptor complexes per µm2 at mature synapses, suggesting that they are assembled in a stereotypic fashion. This GlyR molecular organization is maintained in the hyperekplexia model oscillator despite a decrease in ventral synapse size, indicating that GlyRα1 deficiency does not affect the integrity of the synaptic arrangement as such. Collectively, our results provide the structural basis for understanding the mechanisms underlying receptor availability and the integration of neurotransmitter-induced signals.

Results

Generation of a KI model expressing endogenous levels of mEos4b-GlyRβ

In order to quantify GlyR numbers and their precise distribution at synapses, we generated a knock-in (KI) mouse model expressing endogenous mEos4b-tagged GlyRβ subunits (Figure 1 and Figure 1—figure supplement 1). The β-subunit drives the synaptic localization of the receptor through direct interactions with the synaptic scaffold protein gephyrin at inhibitory synapses (Kim et al., 2006; Meyer et al., 1995). To date, labeling of GlyRβ in situ using immunocytochemistry has proven difficult due to a lack of reliable antibodies that recognize the native β-subunit (only antibodies for Western blotting recognizing the denatured protein are available), which has severely limited the study of the receptor. The coding sequence of mEos4b was inserted in exon 2 of the Glrb gene by homologous recombination supported by a guide RNA (CRISPR-Cas9) to increase recombination efficiency (ICS, Illkirch, France). Specifically, the fluorophore sequence was inserted after the signal peptide and before the N-terminus of the mature GlyRβ subunit, meaning that it does not interrupt the coding sequence of the receptor (Figure 1—figure supplement 1B). The correct insertion was confirmed by amplification and sequencing of genomic DNA. Semi-quantitative RT-PCR revealed that equal concentrations of GlrbEos and the wild-type (WT) transcript (GlrbWT) are expressed in heterozygous animals. When bred to homozygosity, KI animals follow Mendelian inheritance (Figure 1—figure supplement 2A), exhibit normal lifespans (Figure 1—figure supplement 2B), and display no overt phenotype, suggesting that the GlyR expression and/or function are not altered.

Figure 1 with 4 supplements see all
Quantitative confocal imaging of endogenous glycine receptors (GlyRs) in spinal cord tissue.

(A) Representative confocal images of 40 μm spinal cord tissue sections from homozygous (+/+) and heterozygous (+/-) mEos4b-GlyRβ mice (green). Both mice are homozygous for mRFP-gephyrin (magenta). Scale bar = 0.5 mm. (B) Quantification of mEos4b-GlyRβ and mRFP-gephyrin fluorescent intensity of homozygous and heterozygous 2-month-old animals measured in the area indicated by the white square in (A). Plots show median and quartiles. N = 5–6 images from 5 to 6 tissue slices per genotype from 2 mice per age group. (C) Normalized fluorescent decay traces of homozygous mEos4b measured in the area indicated by the white square in (A) over eight consecutive frames. Intensities were binned in the first image and tracked on an individual synapse basis across the eight frames. (D) Mean intensity ± 95% confidence interval of mEos4b-GlyRβ and mRFP-gephyrin measured from distal edge of spinal cord in 2-month-old homozygous mice. Intensities measured in region as indicated by white rectangle in (A). N = 3 images from three tissue slices from two mice per condition. (E) Representative images of homozygous mEos4b-GlyRβ and mRFP-gephyrin at dorsal and ventral synapses. Scale bar = 10 μm. (F) Quantification of numbers of gephyrin clusters across the spinal cord. Plot shows mean ± 95% confidence interval. N = 3 images from three tissue slices from two mice per condition.

To further confirm that GlyR function is not altered by the introduction of mEos4b, we carried out whole-cell recordings in cultured spinal cord neurons of GlrbWT/WT and homozygous GlrbEos/Eos animals (Figure 1—figure supplement 2C). The agonist glycine was applied in a concentration series from 1 µM to 300 μM. The maximal chloride currents at saturating glycine concentrations of 300 μM were not significantly different in GlrbEos/Eos animals, despite a minor increase in the EC50 (GlrbWT/WT 100 ± 5 μM, GlrbEos/Eos 130 ± 9 μM, p=0.0123 t-test). In view of the millimolar concentration of glycine present during synaptic transmission (Beato, 2008; Legendre, 1998), these data indicate that the presence of the N-terminal fluorophore does not affect GlyR function under physiological conditions. Hill coefficients for GlrbWT/WT and GlrbEos/Eos were in the range of 3.5–4, arguing for cooperativity of the subunits during glycine binding. The glycinergic origin of the chloride influx was confirmed by blocking the currents recorded in the presence of 100 μM glycine with 10 μM strychnine.

Quantitative confocal imaging of endogenous mEos4b-GlyRβ and mRFP-gephyrin at spinal cord synapses in tissue

To verify the expression and synaptic targeting of the mEos4b-GlyRβ protein, we carried out quantitative confocal imaging in 40 µm vibratome tissue sections. GlrbEos/Eos animals were crossed with a previously established KI mouse line expressing mRFP-tagged gephyrin to visualize inhibitory postsynaptic sites (Specht et al., 2013). Since the synaptic localization of the GlyR is strictly dependent on its interaction with gephyrin (Patrizio et al., 2017), we expected a high degree of co-localization of the two proteins in the brainstem and spinal cord (Zeilhofer et al., 2005). Indeed, mEos4b-GlyRβ was widely expressed at inhibitory synapses in the thalamus, midbrain, pons, and medulla (Figure 1—figure supplement 3). Very little fluorescence was detected in the forebrain despite the high reported expression of the Glrb transcript (Fujita et al., 1991; Malosio et al., 1991), suggesting that protein levels are controlled by post-transcriptional mechanisms in a region-specific manner, as previously proposed (Weltzien et al., 2012).

In the spinal cord, we observed bright punctate mEos4b-GlyRβ signals throughout the gray matter, with the exception of the superficial laminae of the dorsal horn, where the intensity of the green fluorescence was markedly lower (Figure 1A). The expression of mEos4b-GlyRβ and mRFP-gephyrin in homozygous and heterozygous animals was quantified in confocal images of thoracic and lumbar spinal cord slices at 2 months (Figure 1B) and 10 months of age (Figure 1—figure supplement 4). The same region of the ventral horn, indicated by the white square in Figure 1A, was taken for quantification of mEos4b-GlyRβ and mRFP-gephyrin expression in all conditions. The integrated mEos4b intensity at gephyrin-positive ventral horn synapses was exactly two times higher in GlrbEos/Eos mice than in GlrbEos/WT, demonstrating that both alleles are expressed with the same efficiency, and that the mEos4b fluorophore does not affect the synaptic localization of the receptor complexes. To confirm that the confocal image acquisition was in the linear dynamic range, we bleached the mEos4b fluorophores by repeatedly scanning the same tissue area at constant laser power (Figure 1A, white square), which resulted in a linear decay of pixel intensities over a range of more than 20-fold (Figure 1C).

Across the spinal cord slices, the intensity of synaptic mEos4b-GlyRβ puncta increased from dorsal to ventral both in homozygous (Figure 1D and E) and in heterozygous animals (Figure 1—figure supplement 4). Similarly, mRFP-gephyrin fluorescence was higher and more variable in the ventral horn, suggesting that synapses were on average about twice as big as those in the dorsal horn, despite being fewer in number (Figure 1F and Figure 1—figure supplement 4). The mEos4b/mRFP ratio was relatively equal across the spinal cord with the exception of the superficial layers of the dorsal horn, where GlyR levels were largely exceeded by gephyrin (Figure 1A and D). The lower GlyR-scaffold occupancy of synapses in laminae I–III can be explained by the predominant expression of GABAARs that compete for receptor binding sites at these mixed inhibitory synapses (Alvarez et al., 1996; Lorenzo et al., 2014; Todd et al., 1996).

Dual-color super-resolution imaging of glycinergic spinal cord synapses

To quantify the observed structural differences at super-resolution, we combined radial fluctuation (SRRF) analysis of mRFP-gephyrin and photo-activated localization microscopy (PALM, a form of SMLM) of mEos4b-GlyRβ in spinal cord tissue from double KI animals. Sucrose-impregnated cryosections of 2 µm thickness were prepared from dorsal and ventral tissue and placed on gridded coverslips (Figure 2A). SRRF and PALM images were acquired sequentially. First, mRFP signals were recorded with high-intensity 561 nm laser illumination until all mRFP fluorophores were bleached (10,000 frames). mEos4b was then photoconverted with increasing 405 nm laser intensity and imaged at 561 nm for 25,000 frames until all available fluorophores were exhausted. By acquiring both fluorophores, mRFP and photoconverted mEos4b, in the same emission band (607/36 nm), any chromatic misalignment between the two super-resolved images was eliminated. SRRF reconstruction was carried out on the raw mRFP image sequence and PALM images were generated from individual mEos4b detections using Gaussian peak fitting (Figure 2B). The spatial resolution was estimated using Fourier ring correlation (FRC), which measures the similarity of two images as a function of spatial frequency by comparing the odd and even frames of the raw image sequence. According to this analysis, the spatial resolution of SRRF was 46 nm and that of PALM 21 nm. It should be noted that the synaptic puncta in the SRRF images appear somewhat smaller and brighter due to differences in the reconstruction methods that result in differences in the dynamic intensity range.

Figure 2 with 4 supplements see all
Dual-color super-resolution imaging of mEos4b-GlyRβ and mRFP-gephyrin.

(A) 2 µm cryosections of spinal cord tissue were cut from dorsal and ventral tissue and placed on gridded glass coverslips. (B) Representative photo-activated localization microscopy (PALM) reconstruction of mEos4b-GlyRβ and super-resolution radial fluctuation (SRRF) reconstruction of mRFP-gephyrin at single dorsal and ventral synapses. Scale bar = 500 nm. (C) Intensity correlation quotient (ICQ) of mEos4b-GlyRβ and mRFP-gephyrin in heterozygous and homozygous 2-month-old mice. Plot shows median ± interquartile range. N = 357–604 synapses from 23 dorsal and 23 ventral images from 9 and 8 tissue slices, respectively, from 2 mice per condition. (D) Quantification of GlyR-gephyrin occupancy in 2-month-old homozygous mice (GlrbEos/Eos/GphnmRFP/mRFP). Nonparametric Spearman’s rank shows the same positive correlation at dorsal and ventral synapses. N = 1115 dorsal synapses and 1107 ventral synapses. (E) Pointillist reconstruction of mEos4b-GlyRβ detections. Inset shows a single synapse; red arrows indicate extrasynaptic receptor complexes. Scale bar = 500 nm. (F) Area of dorsal vs. ventral synapses in 2- and 10-month-old homozygous mice. Plot shows median ± interquartile range. N = 234–433 synapses from 20 to 23 images from 7 to 9 tissue slices from 2 mice per age per condition. Nonparametric Kruskal–Wallis ANOVA with Dunn’s multiple comparison test. (G) Number of synapses in dorsal and ventral tissue in 2- and 10-month-old homozygous mice. Plot shows median ± interquartile range. N = 20–23 images from 7 to 9 tissue slices from 2 mice per age per condition. Nonparametric Kruskal–Wallis ANOVA with Dunn’s multiple comparison test. (H) Histogram of the number of glycine receptors (GlyRs) per synapse calculated from the molecular conversion of detections (see Figure 2—figure supplements 3 and 4) in 2-month-old homozygous mice. N = 433 dorsal synapses and 304 ventral synapses from 23 dorsal and 23 ventral images from 9 and 8 tissue slices, respectively, from 2 mice. CV, coefficient of variation. (I) Scatter plot of GlyR density vs. synapse area in 2-month-old homozygous mice shows no difference between dorsal and ventral synapse densities. N = 433 dorsal synapses and 304 ventral synapses. *p<0.05, **p<0.01, ***p<0.001.

The majority of synaptic clusters in the dual super-resolution images were small and spherical or elongated. Larger clusters displayed a variety of morphologies including elongated shapes seen in side view (Figure 2B) as well as convoluted structures, and were more frequently observed in the ventral horn (Figure 2 and Figure 2—figure supplement 1). All mEos4b-GlyR clusters closely matched the mRFP-gephyrin clusters, confirming the localization of the receptors in the postsynaptic membrane. As expected of two directly interacting synaptic components, the degree of co-localization of mEos4b-GlyRβ and mRFP-gephyrin was very high, with mean intensity correlation quotients (ICQ) around 0.3 (Figure 2C and Figure 2—figure supplement 1), a value indicative of close spatial correlation (Li et al., 2004). Minor mismatches between the super-resolution images are explained by the fact that the majority of synapses are seen in cross-section and that the two fluorophores are located on opposite sides of the postsynaptic membrane (~30 nm distance; Specht et al., 2013). There were no obvious differences between the ICQ values of dorsal versus ventral synapses (0.28–0.3), indicating equivalent GlyRβ-gephyrin binding in the two regions. It should be noted that the ICQ reflects relative fluctuations between images and is not sensitive to absolute differences in signal intensities, resulting in similar ICQ values in animals that are heterozygous or homozygous for mRFP-gephyrin (Figure 2C). Quantitative comparison of mEos4b-GlyRβ and mRFP-gephyrin intensities confirmed that the amounts of receptor and scaffold proteins are closely correlated, and that the occupancy of receptor binding sites is the same in the dorsal and ventral horn, independent of synapse size (Figure 2D and Figure 2—figure supplement 1). To estimate the sizes of the synapses, we applied a density threshold to the PALM pointillist images and calculated the areas of the mEos4b-GlyRβ clusters (Figure 2E). The mean synapse area in the ventral horn was larger and more variable than in the dorsal region, both in animals of 2 and 10 months of age (Figure 2F). We also remarked that the overall number of synapses in ventral horn tissue was lower compared to the dorsal horn, significantly so by 10 months (Figure 2G; see also Figure 1F and Figure 1—figure supplement 4). The inverse relationship between synapse number and size is likely due to the presence of different cell types in the two regions, specifically Renshaw cells and large motor neurons in the ventral horn that require strong inhibitory control (e.g., Bhumbra et al., 2014).

Quantification of GlyR numbers and densities at native spinal cord synapses

Given that mEos4b-GlyRβ subunits are expressed at endogenous levels in GlrbEos/Eos animals and that all the fluorophores were exhausted during the PALM recordings, we were able to count the number of GlyRs at spinal cord synapses. The total number of mEos4b-GlyRβ detections at synapses (Figure 2 and Figure 2—figure supplement 2) was converted into molecule numbers taking into account the blinking properties of the fluorophore and the α32 stoichiometry of the pentameric GlyR complex (Durisic et al., 2014; Patrizio et al., 2017). To this aim, the average number of detections per fluorophore (detections/burst) and the fraction of functional fluorophores (probability of detection, Pdet; Figure 2—figure supplement 3) were determined in each set of experiments using extrasynaptic receptor complexes (Figure 2E, red arrowheads). We calculated a median copy number of 114 pentameric GlyR complexes at dorsal horn synapses and twice that number at ventral horn synapses in 2-month-old animals (Figure 2H). Copy numbers were almost identical at 10 months (Figure 2—figure supplement 4), indicating that the glycinergic network was mature at both time points. These numbers exceed estimates derived from electrophysiological recordings in newborn, juvenile, and adult rat spinal cord neurons that suggest the activation of as few as 7 and up to about 110 GlyRs during an average miniature inhibitory postsynaptic current (mIPSC) (Chéry and de Koninck, 1999; Oleskevich et al., 1999; Singer and Berger, 1999; Takahashi, 1992). The high numbers of GlyRs measured by fluorophore counting therefore imply that the available receptors are not saturated by quantal release, which is likely to increase the dynamic range of postsynaptic inhibition (Alvarez, 2017).

Our quantitative PALM data further demonstrate that differences in receptor numbers result from differences in synapse size (Figure 2F). By combining the two parameters, we derived mean GlyR densities of ~2000 µm–2 (Figure 2I and Figure 2—figure supplement 4). Similar receptor densities of 1250 µm–2 and ~2000 µm–2 have been measured at GABAergic synapses in cerebellar stellate cells and in cultured hippocampal neurons, respectively (Liu et al., 2020; Nusser et al., 1997). We saw no differences in the GlyR packing density at dorsal and ventral horn synapses, nor did we find a clear size dependence, as determined by linear regression of all synapses (Figure 2I) and the evolution of the coefficient of variation of GlyR density with respect to synapse area (Figure 2—figure supplement 4). These findings are significant because they indicate that GlyR density is constant and largely independent of neuron type, embryonic origin, or physiological function. Assuming that gephyrin molecules are clustered at densities of up to 9000 µm–2 (Specht et al., 2013), our data also suggest that close to 50% of the receptor binding sites are occupied by GlyRs at native spinal cord synapses, in line with earlier observations of GlyR subunits that were overexpressed in cultured neurons (Patrizio et al., 2017).

Quantitative SR-CLEM of GlyRβ molecular organization

To integrate the results of molecule counting with three-dimensional ultra-structural information and the exact synapse size, we further analyzed dorsal and ventral horn synapses by SR-CLEM. Previously imaged cryosections of GlrbEos/Eos tissue from 10-month-old animals were embedded in epoxy resin, and ultra-thin (70 nm) serial sections were collected on electron microscopy (EM) slot grids with an ultramicrotome (Figure 3A). After osmium tetroxide enhancement, electron micrographs of identified synapses were acquired in all serial sections and registered manually using the coverslip grids and cellular structures as reference (Figure 3B and C). All of the synapses that were both imaged by PALM and reconstructed with EM were functionally mature as judged by the apposition of a single presynaptic terminal containing synaptic vesicles. In line with our PALM data, we found that glycinergic synapses in the ventral horn were substantially larger and more variable in size than those in the dorsal horn (Figure 3C and D). There was good correspondence between the calculated synapse areas in the EM and PALM image reconstructions, even though PALM underestimated the sizes of some large ventral horn synapses (Figure 3—figure supplement 1). This is probably due to the fact that a majority of synapses are tilted, and that the axial component of the area is not captured in the SMLM projections. Whereas most synapses in the dorsal horn were macular, ventral synapses were frequently composed of subdomains (Figure 3E–G). In agreement with earlier studies (Alvarez et al., 1997; Lushnikova et al., 2011; Santuy et al., 2018), the degree of complexity scales with the size of the postsynaptic specialization (Figure 3F), and was taken into account for the calculation of the combined area in the EM serial sections.

Figure 3 with 1 supplement see all
Quantitative super-resolution correlative light and electron microscopy (SR-CLEM) of endogenous glycine receptor (GlyR) molecular organization.

(A) Samples used for photo-activated localization microscopy (PALM) were re-sectioned as serial ultra-thin 70 nm sections for CLEM. (B) Representative dorsal and ventral synapses imaged with PALM (mEos4b-GlyRβ; green) and super-imposed with their corresponding electron micrographs. Scale bar = 500 nm. (C) Representative 3D reconstructions of dorsal and ventral synapses. Green, postsynaptic site; gray, presynaptic bouton. Scale bar = 500 nm. (D) Area of dorsal and ventral synapses measured by EM. Plot shows median ± interquartile range. Nonparametric unpaired two-tailed t-test, Mann–Whitney post hoc. (E) Percentage of total synapses measured in EM with segmented shapes. (F) Comparison of segmentation index with synaptic area in dorsal and ventral synapses. (G) Juxtaposition of a raw electron micrograph and reconstructed PALM/super-resolution radial fluctuation (SRRF) images of sub-synaptic domains (SSDs) in the same ventral synapse. Scale bar = 500 nm. (H) Analysis of GlyR density following correction for EM area measurements. Plot shows mean ± SD. (I) Combined scatter plot of dorsal and ventral synapse densities shows density is independent of synapse size. White, all synapses; black, SSDs. N = 13 dorsal and 23 ventral synapses from 2 mice. *p<0.05, ns, not significant.

The ratio of GlyR copy numbers and the area of the inhibitory postsynaptic specialization obtained by EM resulted in average receptor densities of approximately 2000 µm–2 (Figure 3H). Consistent with our PALM estimates, we did not observe significant differences between synapses in the dorsal and the ventral horn (Figure 2I and Figure 2—figure supplement 4). Furthermore, the GlyR packing density was not dependent on synapse size (Figure 3I), supporting an earlier proposal (Lim et al., 1999). This suggests that GlyRs are assembled in a systematic manner, where receptor numbers increase linearly with synapse size. Since the morphological complexity of synapses increases with size, it can also be concluded that GlyR occupancy at individual subdomains of the postsynaptic sites is uniform. GlyR densities were indeed not significantly different within subclusters of reconstructed synapses (Figure 3G and I, black data points). Together, these findings point to a tight regulation of the architecture of glycinergic synapses across different molecular length scales, where GlyRs are arranged in sub-synaptic signaling units.

GlyR packing density is unaltered in the hyperekplexia mouse model oscillator

Having identified that GlyRs have a stereotypic molecular organization that is maintained throughout adulthood and across synapses in different neuronal circuits, we questioned whether this arrangement is altered in a mouse model of hyperekplexia, a motor-related neurological disease that significantly impacts motor processing in the ventral horn of the spinal cord. The mouse mutation oscillator (Glra1spd-ot/spd-ot) is recessively inherited and causes a microdeletion and frameshift in the TM3-4 intracellular loop of the GlyRα1 subunit, leading to subunit truncation and subsequent loss of functional GlyRs at synapses (Kling et al., 1997). Homozygous oscillator mice do not live past 3 weeks of age (Buckwalter et al., 1994). In contrast, heterozygous animals have a normal lifespan and exhibit a more subtle phenotype. Glra1spd-ot/WT mice display an increased startle reflex and lower GlyRα1 levels (Kling et al., 1997), making them a suitable model for human hyperekplexia. We generated mutant mice that were homozygous for mEos4b-GlyRβ (as described above) and heterozygous for oscillator (GlrbEos/Eos / Glra1spd-ot/WT) as well as WT littermates (GlrbEos/Eos/Glra1WT/WT). In these experiments, inhibitory synapses were detected by immunolabeling of endogenous gephyrin using the mAb7a antibody.

Quantitative confocal imaging in 40 µm vibratome tissue sections showed bright punctate mEos4b-GlyRβ signals localized at synapses (Figure 4A). Across the spinal cord slices, the intensity of synaptic mEos4b-GlyRβ puncta increased from dorsal to ventral in WT animals (Figure 4B), replicating the intensity profile observed previously (Figure 1D). The intensity of synaptic mEos4b-GlyRβ puncta in heterozygous oscillator animals was substantially lower than those of WT littermates, which explains the reductions in membrane levels of GlyR and gephyrin previously observed by Western blotting (Kling et al., 1997). The mEos4b/gephyrin-7a ratio was relatively equal across the spinal cord with greater variation seen in oscillator. No ectopic GlyRβ clusters were detected, meaning that GlyRs and gephyrin always colocalized (Figure 4—figure supplement 1). The number of gephyrin-positive synapses across the spinal cord remained unchanged between WT and oscillator.

Figure 4 with 1 supplement see all
Quantitative confocal imaging of endogenous glycine receptors (GlyRs) in the oscillator mouse model.

(A) Representative confocal images of dorsal and ventral synapses from heterozygous oscillator mice (+/-) compared to homozygous WT (+/+) littermates. All mice are homozygous for mEos4b-GlyRβ (green), with gephyrin-7a immunolabeling (magenta). Scale bar = 10 μm. (B) Mean intensity ± 95% confidence interval of mEos4b-GlyRβ and gephyrin-7a at gephyrin-positive puncta, and numbers of gephyrin clusters measured from distal edge of spinal cord in 2-month-old mice. N = 9–11 images from 9 to 11 tissue slices from 2 mice per genotype.

Using our quantitative PALM approach, we determined the number and size of glycinergic synapses in tissue slices of dorsal and ventral spinal cord, as well as the mEos4b detection density (Figure 5 and Figure 5—figure supplement 1) in order to understand the alterations in glycinergic synapse architecture in this mutant mouse model. The detections per synapse were converted into molecule numbers as described before (Figure 2—figure supplement 3 and Figure 5—figure supplement 2). In WT animals, we observed small and spherical dorsal synapses and larger, elongated ventral synapses, while oscillator synapses appeared small both in dorsal and ventral horn tissue (Figure 5A). This was confirmed by quantitative analysis. The mean synapse area in the ventral horn was significantly larger than in the dorsal region in WT animals (Figure 5B), in agreement with our earlier data (Figure 2F). However, this difference was lost in Glra1spd-ot/WT littermates. The overall number of synapses was lower in ventral horn tissue compared to the dorsal horn in both WT and oscillator animals (Figure 5C). These data suggest that the decrease in functional receptors in heterozygous oscillator mice manifests itself as a reduction in the size of ventral synapses.

Figure 5 with 2 supplements see all
Photo-activated localization microscopy (PALM) imaging of the glycine receptor (GlyR) in the oscillator mouse model.

(A) Representative PALM reconstructions of mEos4b-GlyRβ at single dorsal and ventral synapses in heterozygous (+/-) oscillator and homozygous (+/+) wild-type (WT) littermates. All mice are homozygous for mEos4b-GlyRβ. Scale bar = 500 nm. (B) Area of dorsal vs. ventral synapses in heterozygous oscillator vs. WT littermates. Plot shows median ± interquartile range. N = 120–282 synapses from 24 images from 9 to 11 tissue slices from 2 mice per genotype. Nonparametric Kruskal–Wallis ANOVA with Dunn’s multiple comparison test. (C) Number of synapses in dorsal and ventral tissue in heterozygous oscillator vs. WT mice. Plot shows median ± interquartile range. N = 12–13 images from 2 mice per genotype. Parametric one-way ANOVA with Tukey’s multiple comparison test. (D) Histogram of the number of GlyRs per synapse calculated from the molecular conversion of detections (see Figure 2—figure supplement 3, and Figure 5—figure supplements 1 and 2). N = 282 WT dorsal and 120 ventral synapses, 273 oscillator dorsal and 156 ventral synapses from 24 images from 9 to 11 tissue slices from 2 mice per genotype. CV, coefficient of variation. (E) Scatter plots of GlyR density vs. synapse area show no difference between dorsal and ventral synapse densities in WT and oscillator. N = same as in (D). *p<0.05, ***p<0.001, ns, not significant.

We further quantified the total number of endogenous GlyRs at synapses (Figure 5D) and found that the number and distribution of receptors in WT mice matched that of 2- and 10-month-old WT animals analyzed previously (Figure 2H and Figure 2—figure supplement 4). This is remarkable considering that these experiments were performed independently and subsequent to the previous dataset, attesting to the stability of the measurements. It should also be pointed out that the GlyR copy numbers were indistinguishable from those in mRFP-gephyrin double KI animals, confirming that the presence of the fluorophores does not affect their expression and synaptic clustering. We found similar copy numbers in dorsal synapses of heterozygous oscillator mice compared to WT littermates. However, in ventral tissue, we found a shift towards lower receptor numbers per synapse, suggestive of smaller synapses. By combining the measurements of GlyR copy numbers with synapse area, we could derive the receptor density. Consistent with our previous PALM data, we found a constant receptor density independent of synapse size in WT mice, as judged by the shallow slope of the linear regression (Figure 5E). We observed a similar and constant GlyR packing density in synapses of oscillator mice. This suggests that despite a decrease in the total number of functional GlyRs in the heterozygous oscillator mouse model the molecular organization underlying receptor clustering within the synapse is maintained, in line with the lack of an overt neuromotor phenotype in these animals. In other words, the receptors are assembled stereotypically in the disease model, as in WT, with synapse size consistently scaling with receptor number. Together, our findings describe a highly regulated architecture of glycinergic synapses in both WT animals as well as in a model of synaptic pathology, providing a structural basis of glycinergic signaling.

Discussion

Combining single-molecule PALM imaging, molecular counting, and 3D EM, we have shown that glycinergic synapses in different regions of the spinal cord follow the same structural principle insofar as their receptor-scaffold occupancy and packing densities are the same. This uniformity extends to the sub-synaptic level. The presence of so-called sub-synaptic domains at inhibitory synapses has been shown by super-resolution microscopy (Crosby et al., 2019; Dzyubenko et al., 2016; Pennacchietti et al., 2017; Specht et al., 2013; Yang et al., 2021). However, it remains controversial whether the identified patterns represent the overall structure of the postsynaptic specialization itself or whether they reflect intra-synaptic variations in molecule clustering. Our quantitative SR-CLEM data lend support to the first model, whereby inhibitory postsynaptic sites in the spinal cord are composed of subdomains that determine the distribution of the GlyRs. This organization is achieved through direct interactions between GlyRs and gephyrin, as shown by the close correspondence between the receptors and scaffold proteins. As such, the GlyR subclusters at spinal cord synapses do not constitute SSDs within the postsynaptic membrane in the strict sense since they exhibit uniform binding to the synaptic scaffold (discussed in Yang and Specht, 2019). The stereotypic GlyR density within SSDs observed in our study supports the idea that these structures can instead be equated with the convolutions of the synaptic junction observed by EM (Alvarez et al., 1997; Lushnikova et al., 2011; Peters and Palay, 1996; Santuy et al., 2018). The formation of these convolutions is probably a consequence of gephyrin oligomerization that appears to introduce an asymmetry in the synaptic scaffold. The situation may be different at GABAergic synapses, where the coexistence of gephyrin-dependent and gephyrin-independent clustering mechanisms could lead to the formation of spatially more restricted SSDs containing different GABAAR subtypes (Pennacchietti et al., 2017; Specht, 2020).

The nanoscale organization of inhibitory synapses is the same in glycinergic neurons in the dorsal and the ventral spinal cord, despite their different embryonic origins. Ventral horn synapses are generally larger, more complex, and contain more GlyRs, suggesting that the size of the synapse is differentially regulated in a regional and cell-type-specific manner to adjust the level of glycinergic inhibition. As such, these synapses may be particularly well adapted to motor circuits, assuring strong and reliable inhibition of the postsynaptic neuron (Alvarez, 2017). At the same time, the long and tortuous perimeter of the postsynaptic specialization is likely to accelerate the dynamic exchange of GlyRs and other synaptic components (Chow et al., 2017; Santuy et al., 2018), thereby promoting the molecular plasticity at complex inhibitory synapses (Specht, 2020). Whether glycinergic plasticity results from transient changes in GlyR occupancy and/or from the recruitment of extrasynaptic GlyR-gephyrin complexes (Chapdelaine et al., 2021) has not been proven thus far. However, our data suggest that ultimately it is the size and complexity rather than the GlyR packing density that is dynamically regulated. It can further be argued that large and morphologically complex synapses may be particularly well adapted to integrate fast repetitive, or indeed multivesicular release arising from one or more presynaptic sites, thus providing a strong and reliable inhibition of the postsynaptic neuron while maintaining fast neurotransmitter clearance (discussed in Alvarez, 2017; Rudolph et al., 2015).

In addition, our data show that GlyR density and occupancy do not change between 2 and 10 months of age, indicating that receptor clustering is fully mature by the earlier time point. Studies of normal aging of spinal cord synapses are scarce, and its effect on receptor organization has not been studied. Broadhead and colleagues (Broadhead et al., 2020) report no difference between the number of excitatory synapses in the ventral horn and only a slight increase in dorsal synapses between 2- and 9-month-old mice. Broadly in line with these findings, we found no difference in the number of synapses, synapse area, and GlyR packing density in dorsal and ventral tissue between 2 and 10 months. Thus, glycinergic postsynaptic sites show considerable control over their molecular composition throughout adulthood, further emphasizing the functional significance of their synaptic architecture in both sensory and motor signaling. Our data therefore suggest that a constant GlyR density potentially provides the most efficient organization of the glycinergic postsynaptic site, while enabling the refinement of the size and complexity of the synapse due to ongoing neural activity.

Their molecular organization sets glycinergic synapses apart from excitatory synapses that do not exhibit systematic receptor clustering. Different glutamate receptors are highly variable and occupy separate SSDs within the overall postsynaptic density (PSD) (Goncalves et al., 2020). The number of AMPARs at excitatory synapses can range from essentially 0 (at silent synapses) to more than 100 (Nusser et al., 1998). Within SSDs, the average AMPAR density has been estimated at around 5000 µm–2 (Goncalves et al., 2020), although the absolute values will have to be confirmed using quantitative labeling strategies such as the one described here. Interestingly, AMPAR content of SSDs was shown to vary in response to synaptic plasticity without an apparent change in SSD size (Compans et al., 2021). This suggests that the packing density of AMPARs is not constant even at the nanoscale, which may be the basis for AMPAR plasticity at excitatory synapses.

The neuromotor disease hyperekplexia results from defects in glycinergic inhibition in humans. Several mouse models with analogous mutations in the Glra1 gene recapitulate the phenotype of exaggerated startle reflexes and muscle stiffness (Schaefer et al., 2018). In general, the mutations in the various mouse models are less well tolerated than in humans and often have lethal phenotypes in homozygotes. Hyperekplexia can be recessively or dominantly inherited in humans, but neither has been shown to cause lethality. The oscillator mouse model, whilst lethal in the homozygous form, displays a relatively mild phenotype in heterozygous animals, with a measurable startle reflex and normal lifespan (Kling et al., 1997). Heterozygous oscillator therefore represents a relevant model for the subtler phenotype in humans and the long-term stability of glycinergic synapses beyond the developmental stage at which lethality occurs in homozygous animals. However, it is not known how the reduced GlyR levels in the Glra1spd-ot/WT hypomorph can affect (and sustain) functional motor networks. Our characterization of the molecular organization of heterozygous oscillator synapses shows that GlyR packing follows the same principle as in WT synapses, even though the total number of available functional receptors is reduced, resulting in smaller synapses in the ventral spinal cord. This further emphasizes that the stereotypic arrangement of GlyRs dictates the size of the postsynaptic domain. Most of the synapses that are formed in heterozygous oscillator likely achieve a size threshold capable of sustaining glycinergic signaling without serious motor defects. The lack of fundamental structural changes at glycinergic synapses further suggests that no or only limited compensatory effects take place in oscillator, in agreement with earlier findings that α1β heteropentameric GlyR complexes cannot be compensated for by other subunit configurations, α1-homopentamers, or GABAARs (Schaefer et al., 2012). Our findings thus provide a new perspective into the molecular basis of GlyRα1 deficiency in an animal model of human hyperekplexia.

Taken together, our data show that dorsal and ventral synapses are distinct populations. Ventral horn synapses have much higher GlyR copy numbers, even though receptor density is not different. In contrast to the relatively compact, macular synapses in the dorsal horn, ventral horn synapses achieve a greater receptor number by enlarging the synaptic surface, thus multiplying the sites of signal transmission. These region-specific glycinergic synapse morphologies are likely to underlie functional differences at sensory (dorsal) versus motor (ventral) circuits.

Materials and methods

KI mouse model generation

Request a detailed protocol

The KI mouse line C57BL/6N-Glrbtm1Ics (MGI:6331106) carrying the mutant allele Glrbtm1(Eos4)Ics (MGI:6331065) was created by homologous recombination at the Institut Clinique de la Souris (ICS). Flanked by 5′ and 3′ homology arms of 1.23 kb and 3.49 kb, respectively, the targeting vector encompassed exon 2 of the Glrb gene with an insertion of the coding sequence of mEos4b, as well as a floxed neomycin selection cassette containing the Cre recombinase under control of protamine promoter in introns 2–3. The selection cassette was excised in the F1 generation by germline expression of Cre, resulting in a single loxP site in introns 2–3 of the Glrb locus (Figure 1—figure supplement 1A). The correct insertion of the mEos4b coding sequence was confirmed by sequencing of genomic GlrbEos/Eos tail DNA. Genotyping was done using three primers (primer 1: TACCTTCTTGTTTTCTCTCC; primer 2: GTCTGTTTTCCCTCATAAGG; primer 3: TCGCTTTTGTAAATGATATGG) for the amplification of the mutant GlrbEos (243 bp product) and/or the WT alleles (404 bp).

Purified spinal cord mRNA of GlrbEos/Eos, GlrbEos/WT, and GlrbWT/WT animals was reverse transcribed (primer 6: GGAGTCTAACAGTAATCTGG), and amplified (primer 4: AGGCGCGTCAAACTCGG; primer 5: CCATACCAACCAATGAAAGG). The correct splicing of the mutant transcript was confirmed by sequencing of amplified cDNA. For semi-quantitative RT-PCR, the GlrbEos/Eos mRNA was spiked with WT mRNA at a ratio of 1:2, 1:1, and 2:1 and amplified (Figure 1—figure supplement 1B).

All experiments (with the exception of the data in Figures 4 and 5, Figure 1—figure supplements 1 and 2, Figure 4—figure supplement 1, and Figure 5—figure supplements 1 and 2) were carried out with F2 animals resulting from a cross between the KI line C57BL/6N-Glrbtm1Ics (see above) and a KI mouse line expressing mRFP-tagged gephyrin (GphnmRFP) (Machado et al., 2011) in the C57BL/6J strain, meaning that the mice had a mixed C57BL/6N × C57BL/6J genetic background.

Adult Glra1spd-ot mice (oscillator, JAX stock #000536, RRID:IMSR_JAX:000536) from Jackson Laboratories (Bar Harbor, ME) were transferred to the animal facility of the Institute for Clinical Neurobiology (Würzburg, Germany). Genotyping was done using primer 7: GCCTCCGTGCTTTCTCCCTGC and primer 8: CCAGCCACGCCCCAAAG for the amplification of the mutant Glra1spt-ot (187 bp product) and/or the WT alleles (194 bp). Oscillator mice were backcrossed into the C57BL/6J background for at least 15 generations. Heterozygous GlrbEos/WT animals were crossed with heterozygous oscillator mice for two generations, giving rise to F2 heterozygous oscillator animals that are homozygous for the Glrbtm1(Eos4)Ics allele (Glra1+/spd-ot/GlrbEos/Eos). These animals had a mixed C57BL/6N × C57BL/6J genetic background and were used for the experiments shown in Figures 4 and 5, Figure 4—figure supplement 1, and Figure 5—figure supplements 1 and 2.

All experiments were in accordance with the European Union guidelines and approved by the local veterinary authorities. Animals at IBENS were treated in accordance with the guidelines of the French Ministry of Agriculture and Direction Départementale des Services Vétérinaires de Paris (École Normale Supérieure, Animalerie des Rongeurs, license B 75-05-20). Procedures carried out at the Institute for Clinical Neurobiology were approved by the Veterinäramt der Stadt Würzburg and the Committee on the Ethics of Animal Experiments (Regierung von Unterfranken, Würzburg) and authorized under reference numbers 55.2-2531.01-09/14; 55.2.2-2532.2-949-31.

Primary cultures of spinal cord neurons

Request a detailed protocol

Cultures of mixed spinal cord neurons were prepared at embryonic day 13 (E13) from C57BL/6J WT GlrbWT/WT and GlrbEos/Eos littermates. Spinal cord tissue was trypsinized in trypsin/EDTA (1 mg/ml) and DNase I (final concentration 0.1 mg/ml) for 20 min at 37°C. Trypsinization was stopped with 10% FCS. After trituration, spinal cord neurons were centrifuged at 800 rpm for 10 min. Cells were plated in a 3 cm dish on polylysine-coated coverslips at a density of 2–2.5 × 105 cells/dish. Neurons were kept at 37°C and 5% CO2 in neurobasal medium containing 2 mM L-glutamine and B27 supplement (Thermo Fisher Scientific) with an exchange of half the medium after 4 days in culture.

Electrophysiological recordings

Request a detailed protocol

Spinal cord neuronal cultures at day in vitro 13 (DIV13) were used for patch-clamp recordings in whole-cell configuration. Currents were amplified with a EPC-10 amplifier (HEKA). A laminar flow of increasing agonist concentration (1, 10, 30, 60, 100, and 300 µM glycine, and 100 µM glycine/10 µM strychnine) was applied to the suspended cell using an Octaflow II system (ALA Scientific Instruments), allowing 10–30 ms for equilibration. The external buffer for spinal cord neurons was (in mM) 142 NaCl, 8.1 KCl, 1 CaCl2, 6 MgCl2, 10 glucose, 10 HEPES, pH adjusted to 7.4 with NaOH. To block neuronal excitability and ligand-gated ion channels, the external buffer was complemented with 1 µM TTX, 10 µM bicuculline, 10 µM CNQX, and 25 µM AP-5. The internal buffer was (in mM) 153 KCl, 1 MgCl2, 5 EGTA, 10 HEPES, pH adjusted to 7.4 with CsOH. Recording pipettes were fabricated from borosilicate capillaries with an open resistance of 4–6 MΩ. Currents were measured at a holding potential of –70 mV. All experiments were performed at 22°C. The mean current at each glycine concentration was determined from the peak current amplitudes measured in N = 10–11 cells per genotype from three independent preparations (biological replicates).

Spinal cord and brain tissue preparation and vibratome slices

Request a detailed protocol

Mice were sacrificed at 2 and 10 months of age by perfusion with 4% w/v paraformaldehyde (PFA; Polysciences, EM grade) and 0.1% v/v glutaraldehyde (GA; CliniSciences) in phosphate buffered saline (PBS, pH 7.4). Perfused animals were kept on ice for 30 min, followed by the dissection of the brain and spinal cord in PBS. Tissue was post-fixed in 4% w/v PFA in PBS overnight at 4°C. Brain and spinal cord tissue was rinsed in PBS, cut into smaller segments of thoracic and lumbar regions of the spinal cord and sliced on a vibratome (Leica) at a thickness of 40 μm (for confocal imaging) and 300 μm (for Tokuyasu preparation), and stored in PBS at 4°C.

Confocal imaging and analysis

Request a detailed protocol

In order to label neuronal cells in brain slices (Figure 1—figure supplement 3) and inhibitory synapses in oscillator and WT littermates (Figure 4), free-floating vibratome slices (40 µm thickness) were blocked and permeabilized in PBS containing 0.25% Triton X100 (Sigma) and 0.1% fish gelatin (Sigma) for 1 hr, and immunolabeled with either a primary antibody against NeuN (guinea pig polyclonal, 1:500 dilution, Millipore, #ABN90, RRID:AB_11205592) or gephyrin (mouse monoclonal, mAb7a, 1:500 dilution, Synaptic Systems, #147011, RRID:AB_887717) in PBS containing 0.1% Triton X100 and 0.1% fish gelatin overnight, followed by 3 hr incubation with Alexa Fluor 647-conjugated secondary antibody (donkey anti-guinea pig, 1:1000 dilution, Jackson ImmunoResearch, #706-605-148, RRID:AB_2340476) or Cy3-conjugated secondary antibody (goat anti mouse, 1:1000 dilution, Jackson ImmunoResearch, #115-165-003, RRID:AB_2338680) respectively.

Glass slides (Vector Laboratories) were cleaned with 70% v/v ethanol (Sigma), and vibratome sections were rinsed three times in PBS and mounted onto the glass slides. The glass slides were then briefly rinsed in distilled water and dried. A drop of VectaShield (Vector Laboratories) was added to each spinal cord section and covered with a #1.5 glass coverslip, which was sealed with PicoDent Twinsil Speed (equal weights of catalyst and base). Slides were stored at 4°C for confocal imaging.

Confocal imaging was carried out on a Leica SP8 TCX microscope using a Leica HC PL APO ×40/1.30 NA oil-immersion objective (Leica) and captured in 8-Bit using the Leica LAS-X software (RRID:SCR_013673) with setting HyD3. Images were captured sequentially, with laser illumination at wavelength 570 nm (mRFP, Cy3) imaged first, followed by laser illumination at 491 nm wavelength (mEos4). A cross-section from the dorsal horn to the ventral horn was imaged at a zoom of 5, speed of 25, 512 × 512 pixel (px) format. For decay analysis, eight consecutive frames were captured at a zoom of 5, speed 25, 512 × 512 px format. To tile the whole spinal cord, images were captured in at a zoom of 1, speed 100, 256 × 256 px format.

To ensure alignment of the clusters for the decay traces, images were opened in the image analysis software ICY (RRID:SCR_010587), and the rigid registration plug-in used, taking the first frame of mRFP-gephyrin as reference. The mRFP-gephyrin/ and mEos4b-GlyRβ channels were then separated and the Spot Detector plug-in de Chaumont et al., 2012 used to identify the clusters in each frame in the mRFP-gephyrin channel, with the identified clusters saved as a region of interest (ROI) set. Using the image analysis software Fiji (RRID:SCR_002285), the identified mRFP-gephyrin-positive cluster ROI-Set was used to identify inhibitory synapses in the first frame of the mEos4b-GlyRβ channel. These inhibitory synapses were binned based on mEos4b intensity gray levels (5–12, 13–24, 25–49, 50–74, 75–99, 100–124, 125–255) in frame 1 and a new ROI-Set generated for each bin. Using the frame 1 intensity ROI-Sets, the integrated intensity of mEos4b was then measured at individual clusters across the eight frames. This enabled decay analysis of mEos4b intensity at synapses relative to their starting intensity (see Figure 1C).

In order to analyze the intensity of mRFP-gephyrin and mEos4b-GlyRβ clusters within the spinal cords from mice of different genotypes, the identified mRFP-gephyrin clusters from the first frame of the decay traces (as measured by the ICY Spot Detector plug-in, see above) were used to measure the relative intensity of mRFP-gephyrin and mEos4b-GlyRβ clusters at those locations. The ROI-Set of all mRFP-gephyrin-positive clusters was used in Fiji to identify inhibitory synapses, where the integrated intensity of mRFP and mEos4b was measured for each synapse (Figure 1B).

For the cross-sectional analysis, the mRFP-gephyrin/gephyrin-7a clusters were identified across the imaged tissue using the ICY Spot Detector plug-in, as described above, and saved as an ROI-Set. In Fiji, the integrated intensity of these identified clusters was measured in the mRFP-gephyrin/gephyrin-7a channel and the mEos4b-GlyRβ channel (Figure 1D and F).

Cryosectioning of sucrose impregnated spinal cord tissue

Request a detailed protocol

Sucrose impregnated cryosections were prepared using an ultracryotomy protocol adapted from Tokuyasu, 1973. The 300 μm spinal cord vibratome slices were transferred into a 2.3 M sucrose solution in PBS overnight at 4°C and micro-dissected to isolate gray matter of the dorsal and the ventral horn region. These fragments were placed individually on top of drops of sucrose solution on aluminum EM pins (Leica) and immediately frozen in liquid nitrogen. Sections of 2 µm thickness were sliced on an ultramicrotome (Leica EM UC6) at –80°C and placed onto gridded coverslips (type 1.5H, ibidi GmbH), covered in PBS, and stored at 4°C for a maximum of 5 days before imaging.

Single-molecule localization microscopy

Request a detailed protocol

Sucrose cryosections on gridded coverslips were rinsed once in PBS and imaged in PBS. Dual-color super-resolution images were acquired on an inverted Nikon Eclipse Ti microscope with a ×100/1.49 NA oil-immersion objective, with an additional 1.5× lens in the emission path, using an Andor iXon EMCCD camera (16-Bit, 107 nm pixel size), and NIS-Elements software (Nikon, RRID:SCR_014329). An emission filter 607/36 was chosen for imaging both mRFP-gephyrin and mEos4b-GlyRβ. Brightfield images were taken of the whole grid square identifying tissue structures. Lamp images were taken of the unconverted mEos4b-GlyRβ and mRFP-gephyrin (10 frames of 100 ms, ND8). mRFP-gephyrin movies of 10,000 frames were recorded with HiLo 561 nm continuous laser illumination (output power 50%, 400 mW, 50 ms frames). This was followed by 2 min of 100% 561 nm laser illumination to ensure all mRFP-gephyrin was bleached. Movies of 25,000 frames were recorded with photoconversion of mEos4b-GlyRβ by 0.5 ms pulsed 405 nm laser illumination (gradually increased to 100% by frame 22,000) with continuous 561 nm laser illumination (output power 50%, 400 mW, 50 ms frames). The focal plane was maintained using a Nikon perfect focus system.

SMLM image analysis (SRRF and PALM)

Request a detailed protocol

Frames 100–6000 of the mRFP-gephyrin movies were taken for analysis (to remove saturated frames at the beginning and bleached frames at the end) and were drift corrected and reconstructed using NanoJ-SRRF plug-in for Fiji (Gustafsson et al., 2016).

Quantification of mEos4b-GlyRβ was carried out using a lab script for MATLAB (MathWorks, RRID:SCR_001622). The mEos4b single fluorophores were detected by Gaussian fitting. The resulting pointillist images were drift corrected in the x/y plane using five dense clusters of detections over a sliding window of 2000 frames. Rendered images were produced with a pixel size of 10 nm, sigma 0.01.

FRC was used to estimate the spatial resolution by dividing the odd and even frames of a raw mRFP-gephyrin movie and a raw mEos4b-GlyRβ movie and analyzing the resulting image stacks by SRRF and PALM, respectively. The images reconstructed from the odd and even frames were then compared using the FRC tool of the NanoJ-SQUIRREL plug-in for Fiji (Culley et al., 2018).

The mRFP-gephyrin and mEos4b-GlyRβ-rendered images were aligned by rigid registration using the Fiji plug-in TurboReg. The co-localization of mRFP-gephyrin and mEos4b-GlyRβ was carried out by individually cropping each synapse as separate images. The Fiji plug-in ICQ was then applied to each synapse (Li et al., 2004). The occupancy analysis was analyzed by thresholding the synapses in the mEos4b-GlyRβ images and measuring the intensity of each synapse in both channels. To analyze the PALM mEos4b-GlyRβ clusters, a lab written script for MATLAB (CountMol; Patrizio et al., 2017) was used to identify synapses (minimum number of detections 250, minimum cluster size 200 nm, maximum cluster size 3000 nm) and an intensity threshold of 0.1. For molecule conversion, CountMol was used to identify extrasynaptic receptor complexes (minimum number of detections 5, minimum cluster size 10 nm, maximum cluster size 120 nm). The number of detections per burst (identified as a minimum of two detections, with one burst per 1,000 frames) and the probability of detection Pdet=2N1/N2+2 were calculated and used to convert the detections to mEos4b-GlyRβ molecules (Durisic et al., 2014; Patrizio et al., 2017); see Figure 2—figure supplement 3, and Figure 5—figure supplement 2.

Electron microscopy

Request a detailed protocol

Cryosections used for SMLM imaging on gridded coverslips were postfixed by incubation in 1% osmium tetroxide for 1 hr at 4°C, dehydrated in graded ethanol concentrations, and embedded in araldite epoxy resin. Grid squares imaged in SMLM were identified using the grid pattern imprinted in the resin. Serial ultra-thin 70 nm sections of these regions were cut and transferred onto formvar-coated EM grids (0.432 mm slot grids) using a UC6 ultramicrotome (Leica). Sections were counterstained with 5% uranyl acetate in 70% methanol for 10 min, then washed in distilled water and air-dried before observation on a Philips TECNAI 12 microscope (Thermo Fisher Scientific).

For 3D synapse reconstruction, synapses were manually outlined in each serial section image using Fiji, followed by manual rotation and coarse alignment using the software GIMP (RRID:SCR_003182), then fine alignment of the synaptic area with the Microscopy Image Browser (MIB) software (RRID:SCR_016560). The aligned images were then opened in IMOD software (RRID:SCR_003297) to generate the 3D reconstruction.

The length of the synaptic junction of identified synapses was measured in high-magnification EM images with Fiji software. The total area of the postsynaptic surface was calculated as the cumulative length of the postsynaptic domain in the entire stack of serial sections multiplied by the thickness of each section (70 nm). The segmentation index corresponds to the number of gaps in the postsynaptic site that were detected in the x/y plane of the images or along the z-axis (i.e., an interruption of the postsynaptic site in one or several continuous sections in the stack) and represents an estimate of the morphological complexity of the synapse.

Graphing and statistical analysis

Request a detailed protocol

All graphing and statistical analysis were carried out using the software GraphPad Prism (RRID:SCR_002798) v.8 for all except the electrophysiology experiments that were carried out in v.9. Data were tested for normality of distribution using D’Agostino–Pearson and Kolmogorov–Smirnov tests. Data are represented as dot plots with median ± interquartile range (IQR), or histograms, unless otherwise stated. *p<0.05, **p<0.01, ***p<0.001, ns, not significant.

Data availability

All data generated and analyzed during this study are included in the manuscript and supporting files. The raw data are available for download from the online Figshare repository at https://doi.org/10.6084/m9.figshare.17816747.v1 or via the online repository at https://www.opendata.bio.ens.psl.eu/Super-CLEM-GlyR/.

The following data sets were generated
    1. Maynard SA
    2. Rostaing P
    3. Schaefer N
    4. Gemin O
    5. Candat A
    6. Dumoulin A
    7. Villmann C
    8. Triller A
    9. Specht CG
    (2022) figshare
    Identification of a stereotypic molecular arrangement of endogenous glycine receptors at spinal cord synapses.
    https://doi.org/10.6084/m9.figshare.17816747.v1

References

Decision letter

  1. Julie A Kauer
    Reviewing Editor; Stanford University, United States
  2. John R Huguenard
    Senior Editor; Stanford University School of Medicine, United States

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

[Editors' note: this paper was reviewed by Review Commons.]

https://doi.org/10.7554/eLife.74441.sa1

Author response

Reviewer 1:

In this manuscript Maynard et al. describe a newly generated knockin mouse to study the endogenous distribution of Gly receptors in the spinal cord. Using quantitative confocal imaging and SMLM the distribution and levels of GlyRs at spinal cord synapses is compared between dorsal and ventral horn. They found that levels of synaptic GlyR are higher in dorsal than ventral spinal cord synapses. Nevertheless, the ratio to gephyrin seems constant, except for synapses in superficial layers of the dorsal horn, where gephyrin levels exceeded the levels of GlyRs. There are also fewer, but larger synapses in the ventral horn than in the dorsal horn. These findings are further corroborated by an SR-CLEM approach. Furthermore, it is shown that in a mouse model for hyperekplexia GlyR levels are lower, but still enriched at synapses, and the dorsal-ventral gradient in GlyR expression was maintained. The difference in size of ventral and dorsal synapses observed in WT animals was also lost in the oscillator mouse, suggesting that particularly the ventral synapses are affected. Despite these differences, the density of GlyRs per synapse remained similar.

Major comments:

1. Line 113: "labeling the β-subunit has proven difficult". This statement is unclear and it would be informative for readers to grasp what exactly has been difficult, and why the approach described here overcomes that? Related to that, the authors state "KI animals reach adulthood and display no overt phenotype, suggesting that the presence of the N-terminal fluorophore does not affect receptor expression and function". That is indeed reassuring, but it does not exclude that receptor numbers, function and distribution are altered. As it seems there is no prior literature on tagging the β subunit, additional evidence that the tag does not interfere with receptor trafficking or functioning would be desirable.

We have clarified why it has been difficult to label the GlyR β subunit until now, lines 113-115 “To date, labeling of GlyRβ in situ using immunocytochemistry has proven difficult due to a lack of reliable antibodies that recognize the native β-subunit (only antibodies for Western blotting recognizing the denatured protein are available), which has severely limited the study of the receptor.” Hence it was important to us to generate this knock-in mouse in order to study the endogenous GlyR at synapses, which is the least well studied receptor mediating fast synaptic transmission.

The reviewer makes an important point regarding the labeling of the GlyRβ-subunit with a fluorescent protein that has also been raised by the other reviewers. We have now verified receptor function by patch clamp recordings of glycine currents in whole-cell configuration in spinal cord neuron cultures from the mEos4b KI mouse (new Supplementary Figure S2C). At saturating glycine concentrations of 300 μM we found no difference in chloride influx between mEos4 KI and WT mice. Since glycine concentrations in the synaptic cleft are in the millimolar range during synaptic transmission, these data strongly suggest that glycinergic transmission is not affected by the presence of the mEos4b under physiological conditions, despite a minor shift in the EC50.

There are several other strong arguments that suggest that mEos4b-GlyRb expression, subcellular localization and function are the same as those of the native subunit. Firstly, the mEos4b sequence was inserted after the signal peptide and before the beginning of the coding sequence of the mature β-subunit (Figure S1). Since the mEos4b sequence does not interrupt the coding sequence it is less likely to affect the receptor conformation. Secondly, we did not notice any behavioural phenotypes in animals carrying the GlrbEos allele. At the time of weaning, the genotypes of the pups corresponded to the expected Mendelian frequency (new Figure S2A). Moreover, we did not observe a reduction in live expectancy of GlrbEos/Eos animals (new Figure S2B), demonstrating that the mEos4b-GlyRb does not cause pathology in older animals.

Most importantly, our imaging data (Figure 1-3) provide exhaustive evidence that mEos4b-GlyRb assembles with GlyR α subunits as heteropentameric receptor complexes that are trafficked to the plasma membrane and inserted into the synaptic membrane due to their interaction with the gephyrin scaffold at functional synapses. Using quantitative imaging, we have also shown that homozygous GlrbEos/Eos KI mice have exactly twice the number of receptors at synapses as heterozygous animals, strongly suggesting no interference in receptor trafficking to the plasma membrane and gephyrin binding. As the mEos4b mice were also bred with the oscillator mouse model of hyperekplexia, which is lethal when homozygous, we could further test the combined effect of GlrbEos and GlyRa1spt-ot. The presence of both alleles did not lead to any noticeable phenotypes in heterozygous oscillator mice. On the contrary, both synaptic targeting and the packing density of the receptors were not altered in this model, despite a region-specific reduction in synapse size due to the reduced availability of the intact GlyRa1 subunit.

We believe that these data overwhelmingly support our conclusion that the presence of the mEos4b tag does not alter the structure and function of the receptor, making this mouse model uniquely suited to study the dynamics and regulation of glycinergic synapses in a quantitative manner and at the molecular level.

2. In the Discussion the authors conclude that “Our quantitative SR-CLEM data lend support to the first model, whereby inhibitory PSDs in the spinal cord are composed of sub-domains that shape the distribution of the GlyRs”. This conclusion seems however based on one example image in Figure 3G that is not very convincing. The EM image seems to show two clearly separated PSDs opposed by two distinct active zones. So, although this conclusion is of high interest, more support should be given to substantiate this conclusion. More general, these subsynaptic domains (SSDs) are hardly further explored, but seem relevant for transmission, particularly given that the synaptic pool of GlyRs at these synapses is not saturated by single release events. How general are these SSDs at these synapses?

The representative image in Figure 3G shows two SSDs within the same postsynaptic site with a continuous presynaptic active zone. It should be noted that the PALM/SRRF images were taken of the entire 2 µm thick slice, whereas the electron micrograph shows only a single 70 nm section. We verified throughout the full 3D stack of serial sections that the presynaptic site remains continuous, which it does. We would also like to point out the scale of the image showing that the two SSDs are only around 170 nm apart, i.e. spatially very close. Our conclusions are however not based on this single image but the whole dataset. The graph in Figure 3I shows 3 synapses (out of N = 36), in which the GlyR density at separate SSDs could be quantified, demonstrating that the receptor density is not different between SSDs. The reviewer is correct that we do not further analyse the SSDs beyond their density and the analysis of the segmentation of the postsynaptic sites (Figure 3E-G). Further work on the functional role of SSDs in synaptic transmission is outside the scope of this manuscript and would indeed merit future study.

3. The approach for counting molecules based on the PALM acquisition has been developed in prior publications and seems robust. It would however be worth to present the reader with a bit more background and explain the assumptions of this approach in more detail. Particularly, since counting of mEos4b can be problematic, as there are multiple dark and fluorescent states of this fluorophore that could be influenced by the illumination scheme, see for instance De Zitter et al., Nat Methods 2019. Since the preceding SRRF acquisition already exposes the fluorophore to high and continuous 561-nm laser power this could skew the counting due to unaccounted conversion and perhaps bleaching of mEos4b. In line with this, although throughout the manuscript the term 'absolute copy numbers' is used the reported numbers are at best an estimate based on a number of assumptions. I think the wording 'absolute numbers' is therefore deceiving and should be nuanced.

We have clarified how the molecule conversion is calculated (Figure S7 legend), to provide a more complete description of the way in which the values were obtained. Further we have explained how we calculated the probability of detection. Since the probability of detection accounts for any unconverted or non-functional mEos4b molecules, our molecule counting approach is relatively resistant to potential pre-bleaching of fluorophores. It should be noted, that 561 nm illumination had no obvious effect on the non-converted (green) mEos4b fluorophores, as judged by the fact that the intensity of receptor puncta was unaffected by the SRRF recordings. We appreciate the reviewers point regarding the term ‘absolute copy number’ and we have adjusted our wording throughout the manuscript accordingly.

4. Related, most of the quantifications are in estimating the number of receptors, and not so much the distribution with the PSD. The term "molecular arrangement" – also used in the title – might therefore be misleading, there is in fact little characterization of how GlyRs are placed within the PSD. More focused analysis quantifying the distribution of receptors within the PSD and/or SSDs would strengthen the manuscript.

By estimating the number of receptors and the exact size of synapses, the main conclusion of our study is that receptor density at dorsal and ventral synapses is identical, independent of synapse size, subdomains, or in fact loss of GlyRs in a mouse model of hyperekplexia. This observation clearly relates to how receptors are packed within synapses, and thus describes their molecular arrangement.

5. The reported N is confusing and makes it hard to judge the reproducibility of the data. Sometimes it refers to number of images, sometimes number of synapses, but it is unclear from how many experiments these are drawn. This should be reported more completely (number of animals should be reported at least) and consistently. In figure 1, the N numbers (N=3-5 images) are particularly low and question how consistent these findings are across multiple animals.

We have clarified the N in the figure legends, to reflect the full size of the datasets that have been analysed.

6. The levels of mRFP-Gephyrin seem to differ between the different mouse lines, is this a significant difference?

No significant differences in mRFP-gephyrin levels were found in animals with different mEos4b-GlyRb genotype (Figure 1B). However, expression of mRFP-gephyrin in heterozygous animals is 50% of that in homozygous mRFP-gephyrin KI animals (not shown).

7. The ICQ analysis for co-localization is hardly explained. How do we interpret this parameter? What does an average value of ~0.3 mean? A comparison with sets of proteins that do not overlap as a negative control would strengthen the conclusion.

We have clarified that an ICQ value of 0.3 is indicative of a very high spatial correlation between pixels, and provided a corresponding reference for ICQ analysis (lines 209-210). We would like to point out that the scale of the ICQ is between -0.5 to 0.5, meaning that a value of 0.3 comes close to complete correlation.

Minor comments:

1. "Very little fluorescence was detected in the forebrain, despite the high reported expression of the Glrb transcript". Can the authors expand on this? What would explain this discrepancy?

We have clarified the text to include “suggesting that protein levels are controlled by post-translational mechanisms in a region-specific manner, as previously proposed (Weltzien et al., 2012)” (Lines 152-153). The reason for this discrepancy is not known. However, the distribution of mEos4b expression throughout the brain is as expected, based on the literature.

2. "What region is quantified in Figure 1B? is the same region in all conditions? This should be specified more clearly as the manuscripts presents a clear gradient in expression levels in the spinal cord and thus the location will influence the intensity measurements.

We have explained in the text that this is the region at the centre of the ventral horn identified by the white square in Figure 1A, and that the same region was analysed for all images across all animals. Page 5, lines 160-161 “The same region of the ventral horn, indicated by the white square in Figure 1A was taken for quantification of mEos4b-GlyRβ and mRFP-gephyrin expression in all conditions.”

3. The labeling approach does not differentiate between surface and internal receptors, this should be made more explicit in the text.

Whilst this is correct, we have only analysed mEos4b-positive synapses that had corresponding gephyrin clusters, meaning synapses where receptors are located in the postsynaptic membrane. Indeed we found that all mEos4b clusters imaged colocalised with mRFP-gephyrin clusters. We have adjusted the text accordingly, page 6, line 205-206 “All mEos4b-GlyR clusters closely matched the mRFP-gephyrin clusters, confirming the localization of the receptors in the postsynaptic membrane.”

Significance:

The presented data are interesting and the experiments are technically advanced and carefully performed. Particularly the SR-CLEM approach is technically advanced. The datasets present a quantitatively detailed characterization of spinal cord synapses and will be of interest for researchers working in the field of spinal cord circuitry, as well as super-resolution imaging. The conceptual advance for the field is however somewhat limited. It seems that the presented data confirm the general notion that receptor numbers and synapse size are highly correlated. So, although this manuscript describes very interesting observations, in its present form the manuscript does not provide any new mechanistic insight or significant advance in our understanding of how these synapses operate.

We thank the reviewer for his/her comments relating to the technicality of our manuscript. However we think that the statement “The conceptual advance for the field is however somewhat limited” is unfair, as this level of organisation of inhibitory synapses at the molecular scale has never been achieved before, as pointed out by the other reviewers, and especially not as regards different ages of animals and a disease model that directly affects receptor numbers in a region-specific manner. We therefore believe that our study will have a substantial impact within the fields of synaptic neuroscience as well as quantitative neurobiology.

Referee cross-commenting:

I agree with the other reviewers that this study is technically advanced, but I remain critical towards the extent of conceptual advancement this study brings and there are some important concerns with the presented data that need to be addressed. Nevertheless, indeed many of these concerns can be addressed without additional experiments. As pointed out also by other reviewers additional validation that the fusion proteins are not disrupting their function or organization would be important.

Reviewer 2:

Maynard et al. investigate (inhibitory) glycinergic synapses in mouse spinal cord, which regulate motor and sensory processes. The authors analyse the molecular architecture and ultra-structure of these synapses in native spinal cord tissue using quantitative super-resolution correlative light and electron microscopy. The major finding is that GlyRs exhibit equal receptor-scaffold occupancy and constant absolute packing densities across the spinal cord and throughout adulthood, although ventral and dorsal inhibitory synapses differ in size. Moreover, what the authors call a „stereotypic arrangement" is even maintained in a hypomorphic mutant (oscillator), which is deficient in the adult GlyR a1 subunit.

Specific comments:

1. To reach their conclusions the authors generate two knock-in mouse lines, one with mEOS-labelled GlyR ß-subunit and one with mRFP-labelled gephyrin, a subsynaptic scaffolding protein of inhibitory synapses, which are subsequently crossed. Both changes are not unproblematic, as mutations in the N-terminal end of the GlyR ß subunit polypeptide chain might interfere with the assembly of functional GlyR (consisting of a und ß subunits) and mutations at the N-terminal end of gephyrin interfere with its homo-oligomerization into higher molecular assemblies.

We have demonstrated that the function of mEos4b-GlyRb does not differ significantly from WT GlyRs, by carrying out electrophysiological experiments (new Figure S2C). For a detailed response, please see the response to the first comment of reviewer 1. The mRFP-gephyrin KI strain has been validated and published previously (see Machado et al., 2011, J Neurosci; Specht et al. 2013 , Neuron) and was not specifically generated for this study. The experiments with the oscillator mutant did not include the mRFP-gephyrin allele. In these experiments, the wildtype GlrbEos/Eos (Figure 4, 5) behaves exactly as the GlrbEos/Eos in the double knock-in (Figure 1, 2), further validating the mouse models used.

2. However, in this experimental design both labelled proteins reach postsynaptic membrane specialisations. In case of the ß-subunit quantitative evaluation confirms that heterozygous animals contain only half of the labelled protein as homozygous, which is an indication but not a proof that the correct stoichometry of adult GlyR is maintained. Likewise, mRFP-labelled gephyrin assembles with WT-gephyrin in subsynaptic domains, but it is not clear, if the size and density of the synapses is changed by the knock-in procedure as compared to WT-synapses.

An effect of the mRFP tag on gephyrin clustering can be ruled out, since we observed no difference in synapse size and receptor density in GlrbEos/Eos animals with (Figure 1, 2) and without the GphnmRFP allele (Figure 4, 5, oscillator wild-type controls). Similarly, the synaptic mEos4b-GlyRb levels in heterozygous animals were precisely half those of the homozygous animals, strongly suggesting that the expression and trafficking of the tagged receptor subunit is unchanged, as the reviewer acknowledges. In the absence of any obvious behavioural and/or functional phenotypes (Figure S2) this KI model is in our view is an exceptional tool to study GlyRs expressed at endogenous levels in a cell-type specific manner.

3. Accepting these constraints, which to the knowledge of this reviewer have never been addressed to satisfaction, the authors provide a technically excellent, comprehensive analysis of glycinergic synapses in the spinal cord of double knock-in mice. Therefore, it should be stated in the title, that the investigations were performed with double knock-in instead of „native" spinal cord. Text and figures are clear and accurate and represent the state of the art.

We thank the reviewer for the positive comments regarding the techniques used in the study, and the clarity of the text and figures. We have adjusted the title as requested.

4. Finally, the reviewer would like to raise a minor point: the term postsynaptic density is derived from electron microscopical studies of synapses, where asymmetrical synapses display a „postsynaptic density" but symmetrical synapses do not. The latter were identified as inhibitory synapses and therefore, by definition, inhibitory synapses do not have a postsynaptic density, but rather a postsynaptic membrane specialisation. The use of the term „postsynaptic density" should, therefore, be restricted to excitatory synapses.

We are conscious of the importance of correct definitions and have revised the terminology, referring to “postsynaptic sites”, “postsynaptic domains”, and “postsynaptic specializations” as appropriate throughout the manuscript.

Significance:

The authors provide a state of the art advanced light and electron microscopical analysis of glycinergic synapses in the mouse spinal cord. They suggest a robust “stereotypical” mechanism in place, which guarantees a fixed stoichiometry of relevant components, which is even maintained in a hypomorphic mutant, which is believed to represent a mouse model of human hyperekplexia (startle disease).

Referee cross-commenting:

I would like to corroborate the arguments of the previous reviewer: it is not clear to which extent the fusion proteins influence the measurements, which are technically very advanced and well done, however. The authors do definitely not investigate "native spinal cord" as stated in the title.

The argument concerning fusion proteins must be taken especially serious as the fusions were induced in regions known to be responsible for assembly of glycine receptors and oligomerization of gephyrin.

We have verified the receptor function with electrophysiological recordings and clarified exactly where the fluorescent protein was inserted (see reviewer 1 response). Given the similarity in synapse size, fluorescence intensities and molecule densities observed in neurons expressing different combinations of tagged and native receptors and scaffold proteins, we strongly believe that all animal models used are well suited to the experimental aims of our study.

Reviewer 3:

Glycinergic synapses are the least well understood of synapses that mediate fast synaptic transmission. The manuscript by Maynard et al. adds new information about the structural aspects of these synapses, using PALM and EM imaging of spinal cord synapses from mice at 2 and 10 months. The authors created a knock-in mouse that expresses a tagged GlyRβ subunit, allowing synaptic localization of glycine receptors; all synaptically localized glycine receptors are thought to require the β subunit to be tethered by gephyrin. The authors compare synaptic profiles from: 2 month old vs. 10 month old mice; dorsal vs. ventral horn; and GlyR1-reduced vs. wild type mice. Strikingly, they find a tight relationship across all of these variables between glycine receptor puncta and gephyrin puncta, as well as an apparently constant "packing density" of glycine receptors. They conclude that synaptic extent is likely to be the most important determinant of synaptic strength, as the density of receptors within the postsynaptic density is constant. These results use cutting-edge imaging and are analyzed with care, and add new information to our understanding of these relatively less well characterized synapses.

Major comments:

1. The key conclusions are convincing and the claims appear solid. Additional experiments are not needed to support these claims. The data and the methods are largely presented in such a way that they can be reproduced, although there are minor suggestions for improvement below.

We thank the reviewer for his/her positive comments.

Minor comments:

1. Do the authors have any comment on the requirement during, e.g. LTP, for insertion of a gephyrin-GlyR unit? The lead author has speculated that gephyrin creates "slots" for GlyRs; yet apparently each slot is already filled in the snapshots taken here. How might postsynaptic LTP occur (Kandler group, Kauer group papers)?

Given the reciprocity of GlyR and gephyrin clustering at synapses, the occupancy of binding sites (and in turn the number of available ‘slots’) is dependent on the strength of receptor-scaffold interactions, as discussed previously (Specht 2020, Neuropharmacol). In this study we demonstrate that the density of GlyRs at synapses is constant, which implies that the receptor occupancy is also the same, with the possible exception of mixed inhibitory synapses in the superficial dorsal horn that contain a majority of GABAARs. The PALM/SRRF data are represented as rendered image reconstructions and not as pointillist representations, and the detection of unoccupied binding sites is below the spatial resolution of our approach. However, the high spatial correlation of the signal intensities (ICQ ≈ 0.3) suggests that receptor occupancy is equal between and within synapses. It has previously been established that there are more scaffold proteins than receptors at synapses (Specht et al. 2013, Neuron; Patrizio et al. 2017, Sci Rep). Based on these studies we report that approximately half the gephyrin binding sites are occupied by receptors (lines 262-655). We have also expanded the discussion, describing how shape and size of synapses may affect synaptic transmission, as well as the possible role of receptor-gephyrin interactions in synaptic plasticity at glycinergic synapses.

2. It would be very interesting in the discussion to contrast the present observations with what is known about excitatory synapses (NMDA and AMPAR distributions) and GABAergic synapses. Are the authors at all surprised that receptor packing is constant across conditions? Can the authors speculate on how non-gephyrin binding receptors (homomeric α receptors, which are found in recordings) may function and be tethered to the membrane.

We have included additional information about receptor numbers and distributions at excitatory (lines 428-438) and GABAergic (lines 389-393) synapses in the discussion. So far, homomeric GlyRs composed of α subunits have been found to be exclusively extrasynaptic. As stated on page 4, lines 111-112 the β subunit is required for binding of the GlyR to gephyrin and subsequent anchoring at the synapse. Previous studies have shown exocytosis of receptors to occur at extrasynaptic sites followed by lateral diffusion to synapses. Homomeric GlyRs are therefore most likely targeted to the extrasynaptic plasma membrane where they remain due to the lack of the β subunit.

3. Figure S1. It would be most helpful to quantify this; at the least to include an atlas-like drawing to allow identification of the structures illustrated and containing Glrb; better yet would be quantification of staining in regions where this is strongest.

We have added an atlas indicating the different brain regions expressing mEos4b-GlyRb protein as a new Supplementary Figure S3. The regional expression pattern agrees with the available literature about protein expression of the GlyRb subunit in different brain regions and hence provides further evidence that mEos4b-GlyRb is expressed like the native receptor. Due to the relatively low resolution of the tiled image no accurate quantification was possible. We have however added higher magnification confocal images of representative brain regions expressing varying amounts of GlyRb.

4. The fact that the lower panel in B is labeled as +/+ across all groups is initially confusing; perhaps relabel as mEos4 -/-, +/- and +/+?

We assume that the reviewer is referring to Figure 1B. The genotype of both the GlrbEos and the GphnmRFP allele is now indicated on the x-axes, and the legend has been modified to clarify that all these animals were homozygous for GphnmRFP/mRFP. We have strived to remain consistent throughout the manuscript when referring to genotypes and protein levels.

5. Do gephyrin levels drop in WT mice as well as in the mEosr-GlyRb mouse between 2 and 10 months? Do the authors have any thoughts on this (Supp figure S2)?

We found no differences in gephyrin levels between 2 and 10 months. Figure S2 (now Figure S4C) shows the number of synaptic gephyrin clusters, which was the same at different ages and genotypes.

Significance:

Glycinergic synapses are the least well understood of synapses that mediate fast synaptic transmission. The manuscript by Maynard et al. adds new information about the structural aspects of these synapses, using PALM and EM imaging of spinal cord synapses from mice at 2 and 10 months. The authors created a knock-in mouse that expresses a tagged GlyRβ subunit, allowing synaptic localization of glycine receptors.

This will be of interest to those studying inhibitory synapses, and more broadly to synaptic morphologists, physiologists and imagers for comparison with other synapse types.

My own expertise is NOT in these techniques, but I am a synaptic physiologist with a standing interest in glycinergic synapses; thus I am not providing serious technical critiques.

Referee cross-commenting:

Hi all, I agree with the other two reviewers, and do not have anything else to add.

Reviewer 4:

The authors used a correlative approach and combined photo-activated localization microscopy with electron microscopy to characterise Glycinergic synapses in spinal cord tissue. Some of the major findings are:

– The receptor-scaffold occupancy and packing densities of glycinergic synapses in different regions of the spinal cord are the same.

– Gephyrin clusters in the spinal cord are composed of sub-domains that shape the GlyR clusters.

– Ventral horn synapses are generally larger, more complex (containing a number of gaps) and contain more GlyRs.

–In a mouse model of Hyperekplexia, the number of GlyRs is reduced resulting in smaller synapses in the ventral spinal cord.

Major comments:

1. Are the key conclusions convincing?

Yes.

2. Should the authors qualify some of their claims as preliminary or speculative, or remove them altogether?

No.

3. Would additional experiments be essential to support the claims of the paper? Request additional experiments only where necessary for the paper as it is, and do not ask authors to open new lines of experimentation.

No.

4. Are the suggested experiments realistic in terms of time and resources? It would help if you could add an estimated cost and time investment for substantial experiments.

N/A.

5. Are the data and the methods presented in such a way that they can be reproduced?

Yes.

6. Are the experiments adequately replicated and statistical analysis adequate?

Yes.

Minor comments:

1. Specific experimental issues that are easily addressable.

Please see below.

2. Are prior studies referenced appropriately?

Yes.

3. Are the text and figures clear and accurate?

Yes.

4. Do you have suggestions that would help the authors improve the presentation of their data and conclusions?

Please see below.

5. As the authors pointed out, fusing mEos to the extrasynaptic terminal of GlyRb has been difficult and therefore this construct would benefit the larger scientific community. Figure 1C is a nice imaging control for expression efficiency, however, it is in stark contrast with the lack of functional control. Do authors have any electrophysiological evidence showing that the insertion of mEos4b doesn't modulate channel function? I would assume that the construct would be tested in cell lines before the KI mouse line was created. Was any functional analysis done? If yes, it would be very useful to show it. I do appreciate that the authors used a standard insertion between the 4th and 5th AA in the extracellular domain, which in most cases does not abolish channel function. Given the lack of an obvious phenotype in the KI mouse model, I believe that this is also the case here. However, I disagree with the statement in lines 120-121: "the presence of the N-terminal fluorophore does not affect receptor expression and function." I believe that if there are no electrophysiological measurements of GlyR function, this statement remains speculative. As the authors pointed out in their previous publication: "receptor function and gephyrin binding are not independent properties. Instead, we think that conformational changes triggered at extracellular or intracellular protein domains have downstream consequences on channel opening as well as receptor clustering." In line with this, my concern is that the modulation of channel function by mEos4b could result in an altered cluster size at synapses. There is a large body of literature showing that just one missense mutation in the extracellular domain of ion channel subunits can lead to synaptopathies because the channel function gets modulated, and there is an abundance of similar examples involving mutations of GlyR and GABAAR subunits. In my view, comparing the function of GlyRs incorporating wt-GlyRb and mEos4b-GlyRb subunits is important for the correct interpretation of the main findings of this work and would strengthen the publications.

As the reviewer points out, the insertion of the mEos4b sequence was considered carefully in order to have the least impact on receptor function. GlyR channelopathies are often caused by point mutations within the coding sequence, which is not the case in the GlrbEos allele. Instead, the mEos4b sequence was inserted after the single peptide of GlyRb, duplicating several amino acid residues in order to maintain the correct cleavage site and N-terminus of the mature receptor, and to not interrupt the GlyRb coding sequence (Figure S1B). In order to verify that the mEos4b-tag does not affect GlyR function, we have now carried out electrophysiological experiments (new Figure 2C). For a detailed description please see the response to the first comment of reviewer 1.

6. Line 189: Are the authors making conclusions based on intensity comparison of red mEos4b and mRFP? The title of this section implies that the red form of mEos was compared to mRFP(?) But mEos converts from green to red only partially. Was the probability for conversion taken into account at this point? Please clarify which version of mEos was compared to mRFP.

In line 189 (now 218) we compared the intensities of mRFP-gephyrin with those of converted (red) mEos4b in SRRF / PALM super-resolution images of the synapses (Figure 2D). Since the absolute intensities are altered by the process of image reconstruction, the probability that mEos4b is photoconverted does not have to be taken into account. The constant ratio of the SRRF and PALM image intensities confirms the data in Figure 1D showing that GlyR and gephyrin amounts are highly correlated throughout the spinal cord (with the exception of the superficial layers of the dorsal horn). We have clarified in the text that this analysis was carried out on reconstructed SRRF images of mRFP-gephyrin and PALM images of mEos4, line 202.

7. Line 192: Please clarify how the density threshold was calculated/determined? This is important for the replication of the experiments, and it also has implications for the calculated probability of detection of mEos4b. I am not aware that this probability was calculated before for mEos4b and therefore other researchers may decide to rely on the value calculated here.

We have now clarified in more detail how the probability of detection was calculated (new Supplementary Figure S7 legend).

8. In Figure 2 Gephyrin clusters look consistently smaller than GlyR clusters, which is inconsistent with the published work. I assume that the difference in size is a consequence of different image reconstruction methods(?) However, I would assume that SRRF would have lower resolution than your PALM measurements and that would result in wider Gephyrin clusters. Could you please explain this discrepancy? Also, could you provide an estimate for the image resolution in SRRF and PALM techniques? For SMLM, localization precision would suffice.

We have provided an estimate of the resolution of the two techniques using Fourier ring correlation, which gave 46 nm for SRRF and 21 nm for PALM. Additionally we have precised the discrepancy between reconstruction methods, page 6, lines 194-200 “The spatial resolution was estimated using Fourier ring correlation (FRC), which measures the similarity of two images as a function of spatial frequency by comparing the odd and even frames of the raw image sequence. According to this analysis, the spatial resolution of SRRF was 46 nm and that of PALM 21 nm. It should be noted that the synaptic puncta in the SRRF images appear somewhat smaller and brighter due to differences in the reconstruction methods that result in differences in the dynamic intensity range.”

9. Why is the data in Figure 5D and E represented as Detections/Synapse instead of GlyRs/Synapse? Could you please re-plot this so that a comparison with Figure 2H and I is straightforward?

We have converted the detections to receptor copy numbers as requested (Figure 5D,E).

10. Figure S5C: for P=0.5, 2=0.25. Please correct. Also, I assume that the second graph is what would be observed experimentally for dimers and P=0.5. Please clarify in the figure caption.

This was a mistake and has been corrected. We have also clarified which parts of the calculations are theoretical and which values were derived from our experimental data. We have provided a more detailed description in the figure legend of Supplementary Figure S7.

11. Line 606: Please provide a complete derivation of this formula.

We have provided a full derivation of this formula (new Figure S7C).

Significance:

The work described here seem to be a natural progression of a publication by Patrizio et al., 2017 that came out from the same laboratory. This study uses advanced methodologies in the imaging space to visualise and characterise Glycinergic synapses in spinal cord tissue. The experiments described here are technically demanding as evidenced by the relatively small number of publications describing super-resolution measurements in tissue samples. Even more rare are studies that attempt to do single protein counting in neuronal culture and tissue sections. Therefore, I believe that this work brings significant technical advancement in the field of super-resolution and corelative microscopy. The findings are also highly significant for all fields of neuroscience in which the structure of inhibitory Glycinergic synapse is relevant, ranging from the fundamental understanding of inhibitory synapse function to pathologies involving Glycinergic signalling.

I have substantial experience in different microscopy methods, including quantitative super-resolution microscopy based on single molecule counting. My background also covers the structure and function of GABAA and Glycine receptors using electrophysiology. I am familiar with the methods used in electron microscopy and the process of creating KI mouse lines, however I don't have hands-on experience in these fields.

https://doi.org/10.7554/eLife.74441.sa2

Article and author information

Author details

  1. Stephanie A Maynard

    Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, Paris, France
    Contribution
    Formal analysis, Writing – original draft, Writing – review and editing, Conceptualization, Funding acquisition, Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7838-3676
  2. Philippe Rostaing

    Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, Paris, France
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
  3. Natascha Schaefer

    Institute for Clinical Neurobiology, University Hospital, Julius-Maximilians-University, Wuerzburg, Germany
    Contribution
    Formal analysis, Investigation, Funding acquisition
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9743-1963
  4. Olivier Gemin

    Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3210-7876
  5. Adrien Candat

    Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  6. Andréa Dumoulin

    Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  7. Carmen Villmann

    Institute for Clinical Neurobiology, University Hospital, Julius-Maximilians-University, Wuerzburg, Germany
    Contribution
    Formal analysis, Funding acquisition, Investigation, Resources
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1498-6950
  8. Antoine Triller

    Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, Paris, France
    Contribution
    Funding acquisition
    For correspondence
    triller@biologie.ens.fr
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7530-1233
  9. Christian G Specht

    1. Institut de Biologie de l'ENS (IBENS), École Normale Supérieure, CNRS, Inserm, PSL University, Paris, France
    2. Diseases and Hormones of the Nervous System (DHNS), Inserm U1195, Université Paris-Saclay, Paris, France
    Contribution
    Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing
    For correspondence
    christian.specht@inserm.fr
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6038-7735

Funding

H2020 European Research Council (Plastinhib)

  • Antoine Triller

Agence Nationale de la Recherche (Synaptune)

  • Antoine Triller

Agence Nationale de la Recherche (Syntrack)

  • Antoine Triller

Labex (Memolife)

  • Antoine Triller

France Bio-Imaging

  • Antoine Triller

Deutsche Forschungsgemeinschaft (VI586)

  • Carmen Villmann

Fondation pour la Recherche Médicale (SPF201809007132)

  • Stephanie Maynard

Bavarian State Ministry of Science and the Arts and the University of Würzburg (Graduate School of Life Sciences (GSLS))

  • Natascha Schaefer

European Research Council (MICRO-COPS)

  • Antoine Triller

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

SAM was supported by a Fondation pour la Recherche Médicale (FRM) postdoctoral fellowship (SPF201809007132). NS was supported by funds of the Bavarian State Ministry of Science and the Arts and the University of Würzburg to the Graduate School of Life Sciences (GSLS), University of Würzburg. Research in our laboratory at IBENS was funded by the European Research Council (ERC, Plastinhib and Microcops), Agence Nationale de la Recherche (ANR, Synaptune and Syntrack), Labex (Memolife) and France-BioImaging (FBI). CV was supported by the Deutsche Forschungsgemeinschaft (DFG, VI586). We acknowledge the use of the EM platform of IBENS. We thank Pascal Legendre (ENP, Inserm) for insightful discussions and Constant Morez and Nadine Schibille for helpful comments on the manuscript. We also thank Marie-Christine Birling and Eve Geronimus from the Institut Clinique de la Souris (ICS, Illkirch, France) for the generation of the C57BL/6N-Glrbtm1Ics mouse line.

Ethics

All experiments were in accordance with European Union guidelines and approved by the local veterinary authorities. Animals at IBENS were treated in accordance with the guidelines of the French Ministry of Agriculture and Direction Départementale des Services Vétérinaires de Paris (École Normale Supérieure, Animalerie des Rongeurs, license B 75-05-20). Procedures carried out at the Institute for Clinical Neurobiology were approved by the Veterinäramt der Stadt Würzburg and the Committee on the Ethics of Animal Experiments (Regierung von Unterfranken, Würzburg) and authorized under reference numbers 55.2-2531.01-09/14; 55.2.2-2532.2-949-31.

Senior Editor

  1. John R Huguenard, Stanford University School of Medicine, United States

Reviewing Editor

  1. Julie A Kauer, Stanford University, United States

Version history

  1. Preprint posted: September 10, 2021 (view preprint)
  2. Received: October 5, 2021
  3. Accepted: November 15, 2021
  4. Accepted Manuscript published: December 8, 2021 (version 1)
  5. Version of Record published: January 11, 2022 (version 2)
  6. Version of Record updated: January 13, 2022 (version 3)

Copyright

© 2021, Maynard et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,176
    Page views
  • 179
    Downloads
  • 12
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Stephanie A Maynard
  2. Philippe Rostaing
  3. Natascha Schaefer
  4. Olivier Gemin
  5. Adrien Candat
  6. Andréa Dumoulin
  7. Carmen Villmann
  8. Antoine Triller
  9. Christian G Specht
(2021)
Identification of a stereotypic molecular arrangement of endogenous glycine receptors at spinal cord synapses
eLife 10:e74441.
https://doi.org/10.7554/eLife.74441

Further reading

    1. Neuroscience
    Stijn A Nuiten, Jan Willem de Gee ... Simon van Gaal
    Research Article

    Perceptual decisions about sensory input are influenced by fluctuations in ongoing neural activity, most prominently driven by attention and neuromodulator systems. It is currently unknown if neuromodulator activity and attention differentially modulate perceptual decision-making and/or whether neuromodulatory systems in fact control attentional processes. To investigate the effects of two distinct neuromodulatory systems and spatial attention on perceptual decisions, we pharmacologically elevated cholinergic (through donepezil) and catecholaminergic (through atomoxetine) levels in humans performing a visuo-spatial attention task, while we measured electroencephalography (EEG). Both attention and catecholaminergic enhancement improved decision-making at the behavioral and algorithmic level, as reflected in increased perceptual sensitivity and the modulation of the drift rate parameter derived from drift diffusion modeling. Univariate analyses of EEG data time-locked to the attentional cue, the target stimulus, and the motor response further revealed that attention and catecholaminergic enhancement both modulated pre-stimulus cortical excitability, cue- and stimulus-evoked sensory activity, as well as parietal evidence accumulation signals. Interestingly, we observed both similar, unique, and interactive effects of attention and catecholaminergic neuromodulation on these behavioral, algorithmic, and neural markers of the decision-making process. Thereby, this study reveals an intricate relationship between attentional and catecholaminergic systems and advances our understanding about how these systems jointly shape various stages of perceptual decision-making.

    1. Neuroscience
    Manfred G Kitzbichler, Daniel Martins ... Neil A Harrison
    Research Article Updated

    The relationship between obesity and human brain structure is incompletely understood. Using diffusion-weighted MRI from ∼30,000 UK Biobank participants, we test the hypothesis that obesity (waist-to-hip ratio, WHR) is associated with regional differences in two micro-structural MRI metrics: isotropic volume fraction (ISOVF), an index of free water, and intra-cellular volume fraction (ICVF), an index of neurite density. We observed significant associations with obesity in two coupled but distinct brain systems: a prefrontal/temporal/striatal system associated with ISOVF and a medial temporal/occipital/striatal system associated with ICVF. The ISOVF~WHR system colocated with expression of genes enriched for innate immune functions, decreased glial density, and high mu opioid (MOR) and other neurotransmitter receptor density. Conversely, the ICVF~WHR system co-located with expression of genes enriched for G-protein coupled receptors and decreased density of MOR and other receptors. To test whether these distinct brain phenotypes might differ in terms of their underlying shared genetics or relationship to maps of the inflammatory marker C-reactive Protein (CRP), we estimated the genetic correlations between WHR and ISOVF (rg = 0.026, P = 0.36) and ICVF (rg = 0.112, P < 9×10−4) as well as comparing correlations between WHR maps and equivalent CRP maps for ISOVF and ICVF (P<0.05). These correlational results are consistent with a two-way mechanistic model whereby genetically determined differences in neurite density in the medial temporal system may contribute to obesity, whereas water content in the prefrontal system could reflect a consequence of obesity mediated by innate immune system activation.