Abstract
The spatiotemporal organization of neurotransmitter receptors in postsynaptic membranes is a fundamental determinant of synaptic transmission and information processing by the brain. Using four independent super-resolution light imaging methods and EM of genetically tagged and endogenous receptors, we show that, in rat hippocampal neurons, AMPARs are often highly concentrated inside synapses into a few clusters of ∼70 nm that contain ∼20 receptors. AMPARs are stabilized reversibly in these nanodomains and diffuse freely outside them. Nanodomains are dynamic in their shape and position within synapses and can form or disappear within minutes, although they are mostly stable for up to 1 h. AMPAR nanodomains are often, but not systematically, colocalized with clusters of the scaffold protein PSD95, which are generally of larger size than AMPAR nanoclusters. PSD95 expression level regulates AMPAR nanodomain size and compactness in parallel to miniature EPSC amplitude. Monte Carlo simulations further indicate the impact of AMPAR concentration in clusters on the efficacy of synaptic transmission. The observation that AMPARs are highly concentrated in nanodomains, instead of diffusively distributed in the PSD as generally thought, has important consequences on our understanding of excitatory neurotransmission. Furthermore, our results indicate that glutamatergic synaptic transmission is controlled by the nanometer-scale regulation of the size of these highly concentrated nanodomains.
Introduction
The molecular composition and organization of synapses is a fundamental determinant of synaptic transmission. Considerable efforts have been devoted to determine precisely the spatial organization of presynaptic and postsynaptic molecular elements and to correlate it with synaptic activity. Presynaptic vesicles contain a limited number of glutamate molecules, whose release creates a transient local glutamate gradient (Lisman et al., 2007). As AMPA receptors (AMPARs) have low affinity for glutamate, only a fraction of synaptic AMPARs are activated by this gradient (Liu et al., 1999). Thus, the precise localization of AMPARs with respect to the presynaptic terminal and glutamate release sites is crucial to tune synaptic transmission. Modeling and experimental data indicate that heterogeneity of AMPAR distribution may have a strong impact on synaptic transmission (Glavinovic, 1999; Franks et al., 2002; Franks et al., 2003; Xu-Friedman and Regehr, 2004; Wu et al., 2007; Tarusawa et al., 2009; Budisantoso et al., 2012). Three categories of approaches have been used to determine AMPAR organization: light or electron microscopy (EM) and electrophysiology. All have concluded that AMPARs are concentrated in the postsynaptic density (PSD), directly opposed to presynaptic terminals (Baude et al., 1995; Kharazia and Weinberg, 1997; Takumi et al., 1999b). AMPARs are also present in the extrasynaptic membrane and in intracellular pools, and constant exchange of AMPARs between these different pools by lateral diffusion and recycling tunes receptor numbers at synapses and regulates synaptic transmission (Luscher et al., 1999; Makino and Malinow, 2009; Petrini et al., 2009; Choquet, 2010).
The precise distribution of AMPARs in nanometer-sized subdomains has so far only been studied in fixed samples. Immunogold labeling of AMPARs for EM is the highest-resolution imaging technique (Kharazia et al., 1996; Petralia et al., 1997; Nusser et al., 1998; Takumi et al., 1999a). However, the exact distribution of AMPAR density is still controversial. Studies demonstrate either that synaptic AMPARs are localized at the PSD periphery (Kharazia et al., 1996; Matsubara et al., 1996; Bernard et al., 1997; Kharazia and Weinberg, 1997; Chen et al., 2008) or randomly distributed (Nusser et al., 1994; Masugi-Tokita and Shigemoto, 2007; Masugi-Tokita et al., 2007). Pre-embedding immunogold EM and replica-based labeling indicate the existence of AMPAR clusters on the membrane (Masugi-Tokita et al., 2007; Tarusawa et al., 2009; Tao-Cheng et al., 2011; Budisantoso et al., 2012). Complementary to EM, super-resolution light microscopy of immunofluorescently labeled samples, which allows breaking the diffraction barrier to ∼250 nm, has revealed the molecular organization of postsynaptic molecules at nanometric resolution (Dani et al., 2010). In many synapses, AMPARs were organized centrally, whereas in others they were organized laterally to the PSD. We recently performed mathematical analysis of live super-resolution tracking of AMPAR movements and found that AMPARs move to attracting interaction potential wells at the subdiffraction level in hippocampal dendrites (Hoze et al., 2012).
Here, we performed a detailed analysis of AMPAR and PSD95 localization and mobility at synapses using a combination of high-resolution imaging approaches. We demonstrate that AMPARs are strongly nonuniformly organized in synapses.
Materials and Methods
Constructs.
mEos2::GluA1 (noted Eos::GluA1), GluA2::tdEos (noted Eos::GluA2), and Homer1c::Cerulean were subcloned from SEP::GluA1, GluA2::GFP, Homer1C::Dsred, and Homer1c::GFP, which were described previously (Saglietti et al., 2007; Opazo et al., 2010). Eos::TEV::GluA1 was modified from mEos2::GluA1 by inserting a Tobacco Etch Virus Protease recognition site comprising a seven amino acid sequence (Glu-Asn-Leu-Tyr-Phe-Gln-Gly, cleaving between Gln and Gly), for removing Eos tag from surface-expressed Eos::GluA1 of live neurons. PSD95::mEos2 (noted Eos::PSD95) was subcloned from PSD95::EGFP (used for overexpression experiments), which was described previously (Mondin et al., 2011). The GFP-tagged shRNA construct designed to silence expression of PSD95 (SH::PSD95) and the PSD95 rescue construct (Rescue::PSD95), which expresses shRNA against PSD95 simultaneously with an shRNA-resistant form of GFP-tagged PSD95, were characterized previously (Schluter et al., 2006; Mondin et al., 2011). A generic scrambled RNA (SH::Scramble) vector containing a noneffective mammalian protein coexpressed either with GFP (SH::Scramble::GFP) or Homer1c::GFP (SH::Scramble::Homer1c) was used as a negative control to compare the effect of SH::PSD95::GFP. The compensatory mutant for stagazin and PSD95 were described and characterized previously (Schnell et al., 2002; Bats et al., 2007). The Stargazin/PSD95 compensatory mutants fused with GFP were obtained by modifying the −2 position threonine of Stargazin to phenylalanine (StgT231F) and the 225 position histidine of PSD-95 to valine (PSD95H225V).
Cell culture and transfection.
Preparation of cultured neurons was performed as previously described (Heine et al., 2008). Hippocampal neurons from 18-d-old rat embryos of either sex were cultured on glass coverslips following the Banker protocol (Kaech and Banker, 2006). Neurons were transfected using Effectene at 9–11 d in vitro (DIV) or by Nucleofection at the time of plating (Nucleofector II Device, Lonza Cologne) with mEos2::GluA1, Eos::TEV::GluA1, GluA2::tdEos, and PSD95::mEos2, alone or in combination with Homer1c::Cerulean, HA::GluA1, or Homer1c::GFP, and experiments were performed at 14–20 DIV.
Single particle tracking photoactivation localization microscopy (sptPALM).
Cells were imaged at 37°C in an open chamber (Ludin chamber, Life Imaging Services) mounted on an inverted motorized microscope (Nikon) equipped with a 100× 1.49 NA PL-APO objective and a perfect focus system, allowing long acquisition in oblique illumination mode. Imaging was performed for <20 min unless otherwise stated in an extracellular solution as in Petrini et al. (2009). Cells expressing Eos::GluA1 or Eos::GluA2 were photoactivated using a 405 nm laser (Omicron), and the resulting photoconverted single molecule fluorescence was excited with a 561 nm laser (Cobolt). Both lasers illuminated the sample simultaneously. Their respective power was adjusted to keep the number of the stochastically activated molecules constant and well separated during the acquisition. Laser intensities were tuned to leave the single molecule fluorescent during multiple frames before bleaching. The fluorescence was collected by the combination of a dichroic and emission filters (D101-R561 and F39–617, respectively, Chroma) and a sensitive EMCCD camera (Evolve, Photometric). The acquisition was steered by MetaMorph software (Molecular Devices) in streaming mode at 50 frames per second (20 ms exposure time) using a 200 × 200 pixels region of interest. The native nonactivated fluorescence form of Eos molecules was excited using a conventional GFP filter cube (ET470/40, T495LPXR, ET525/50, Chroma). Homer1c::Cerulean fluorescent protein was observed using a CFP filter (ET436/20, T455LP, ET480/40, Chroma). Multicolor fluorescent micro-beads (Tetraspeck, Invitrogen) were used as fiduciary markers to register long-term acquisitions and correct for lateral drifts.
Universal point accumulation in nanoscale topography (uPAINT).
uPAINT was applied as reported previously (Giannone et al., 2010). Dissociated neurons were nucleofected at the time of plating with Homer1c::GFP, allowing synapse identification. Experiments took place at 13–15 DIV. Coverslips were mounted on a Ludin chamber filled with 600 μl of extracellular solution at 37°C as in Petrini et al. (2009). A low concentration of ATTO-647 nm-coupled anti-GluA2 antibodies was then added to the chamber. Stochastic labeling of endogenous GluA2-containing AMPA receptors by dye-coupled antibodies allowed us to record thousands of trajectories lasting longer than 1 s. The amount of nonspecific binding of antibodies was evaluated using hippocampal cultures derived from GluA2 knock-out mice and was found to be ∼15%. Potential effects of the antibody binding on receptor dynamics was tested by comparing two different GluA2 antibodies and results were found to be similar (see Fig. 6A,B). Multicolor fluorescent micro-beads (Tetraspeck, Invitrogen) were used as fiduciary markers to register long-term acquisitions and correct for lateral drifts.
Single molecule localization and tracking.
A typical single-cell sptPALM or uPAINT experiment acquired with the microscope setup and protocol described above produced a set of 20,000 images that were analyzed to extract molecule localization and dynamics. Single molecule fluorescent spots were localized in each image frame and tracked over time using a combination of wavelet segmentation (Izeddin et al., 2012) and simulated annealing (Racine et al., 2006, 2007) algorithms. Under the experimental conditions described above, the image resolution of the experimental setup was quantified to 46.6 nm for sptPALM. The localization accuracy, which depends on the image signal-to-noise ratio and the segmentation algorithm (Kubitscheck et al., 2000; Cheezum et al., 2001; Izeddin et al., 2012), was determined experimentally using fixed samples expressing Eos-tagged proteins. We analyzed 133 2D distributions of single molecule positions belonging to long trajectories (>50 frames) by bidimensional Gaussian fitting, and the resolution was determined to be 2.3σ. The software package used to derive quantitative data on protein localization and dynamics is custom written as a plug-in running within the MetaMorph software environment. For the trajectory analysis, synapses were identified by wavelet image segmentation of the Homer1c postsynaptic marker. The corresponding binary mask was then used to sort single particle data analyses to specific synaptic regions. Dendrites were identified as subcellular regions lacking the synaptic marker.
Live-cell surface staining and stimulated emission depletion microscopy (STED).
The monocolonal antibody against an extracellular epitope of GluA2 was prepared as described previously (Giannone et al., 2010). Primary neuronal cultures transfected with Homer1c::GFP were incubated with the monoclonal anti-GluA2 antibody (3 μg/ml) in the neuronal growth medium for 4–6 min at 37°C, and then fixed in 4% paraformaldehyde and 4% sucrose in PBS. Fixed samples were rinsed in PBS and then blocked in PBS containing 1% BSA. The primary antibodies were then revealed by incubating ATTO647N (ATTO-tec) coupled anti-mouse IgG secondary antibodies for 30 min at room temperature. Finally, samples were rinsed in PBS and mounted in Vectashield mounting medium. We used a commercial STED inverted microscope (DMI6000 TCS SP5 AOBS, Leica) to obtain super-resolved images of GluA2-ATTO647N-labeled neurons expressing Homer1c::GFP as a postsynaptic marker. The Homer1c::GFP images were recorded using the regular confocal mode on the same microscope. A spatial resolution of 63.7 ± 1.2 nm was measured using 40 nm crimson beads.
Direct stochastic optical reconstruction microscopy (dSTORM).
Primary neuronal cultures were incubated with rabbit-anti-GluA1 antibody (targeted to the extracellular region RTSDSRDHTRVDWKRC, Agro-Bio) for 20–30 min. They were then fixed using 4% paraformaldehyde and 4% sucrose in PBS, and washed with PBS. They were incubated with NH4Cl 50 mm for 30 min before permeabilization. They were permeabilized using 0.1% Triton X-100 and incubated with PBS containing 1% BSA for 30 min. They were then incubated with mouse-anti-PSD95 antibody (MA1–046, Thermo Fisher Scientific) for 30 min and washed several times with PBS containing 1% BSA. The primary antibodies were then revealed by incubating Alexa-647-coupled anti-rabbit IgG secondary (A21245, Invitrogen) and RhodamineRed coupled anti-mouse secondary antibodies (715-295-151, Jackson ImmunoResearch Laboratories) for 30 min at room temperature. For single-color dSTORM experiments, primary neuronal cultures were incubated with mouse-anti-GluA2 antibody (Giannone et al., 2010) for 5–7 min. They were then fixed, as described previously. The primary antibodies were then revealed by incubating Alexa-647 coupled anti-mouse IgG secondary (Invitrogen) for 30 min at room temperature.
The stained coverslips were imaged the next day at room temperature in a closed chamber (Ludin Chamber, Life Imaging Services) mounted on an inverted motorized microscope equipped with a 100× 1.49 NA PL-APO objective and a Perfect Focus System (Nikon), allowing long acquisition in oblique illumination mode. Imaging was performed in an extracellular solution containing reducing agents and oxygen scavengers. For dSTORM, ensemble fluorescence of Rhodamine Red and Alexa-647 was first converted in to dark state using a 561 nm laser or 640 nm laser (Coherent) at 30–50 kw/cm2 intensity. Once the ensemble fluorescence was converted into the desired density of single molecules per frame, the laser power was reduced to 7–15 kw/cm2 and imaged continuously at 50 fps for 20,000 frames. The level of single molecules per frame was controlled by using a 405 nm laser (Omicron). The dyes were sequentially imaged (Alexa-647 first and followed by Rhodamine) to collect the desired single molecule frames. The laser powers were adjusted to keep an optimal level of stochastically activated molecules during the acquisition. Both the ensemble and single molecule fluorescence was collected by the combination of a dichroic and emission filter (D101-R561 and F39-617, respectively, Chroma; and quad-band dichroic filter (Di01-R405/488/561/635, Semrock). The fluorescence was collected using a sensitive EMCCD (Evolve, Photometric). Single molecule localization and reconstruction were performed online using automatic feedback control on the lasers, enabling optimal molecule density during the acquisition (Kechkar et al., 2013). The acquisition and localization sequences were driven by MetaMorph software (Molecular Devices) in streaming mode at 50 frames per second (20 ms exposure time) using an area equal to or <256 × 256 pixels region of interest. We used multicolor fluorescent microbeads (Tetraspeck, Invitrogen) as fiduciary markers to register long-term acquisitions and correct for lateral drifts and chromatic shifts. A spatial resolution of 14 nm was measured using centroid determination on 100 nm Tetraspeck beads acquired with similar signal-to-noise ratio than dSTORM single molecule images.
Synaptic surface area quantification.
Epifluorescence images of the synaptic marker protein Homer1c (Homer1c::GFP or Homer1c::Cerulean) were used to differentiate between synapses and the rest of the dendrite. Epifluorescence images of Homer1c were first thresholded and segmented using the morphometric image analysis module of MetaMorph software (Molecular Devices) for structures >0.02 μm2. Morphological features, such as area, length, and breadth of each segmented structure, were exported to calculate their respective distributions. The mean synaptic area of Homer1c labeled structures was measured to be 0.11 ± 0.04 μm2, which is consistent with the surface area of the spine head. The Homer1c segmented images were then used to define synapses area for single particle tracking experiments.
Super-resolution cluster analysis.
AMPAR nanodomains and PSD95 subclusters were identified from super-resolution images by custom software written as a plug-in running inside MetaMorph. Single-molecule-based super-resolution images were reconstructed from the 20,000 frames before being analyzed, whereas STED bulk images were analyzed raw. Nanodomains, which corresponded to clustered areas where the signal density was higher, were first identified by wavelet segmentation. Nanodomain number and dimensions were then computed by 2D anisotropic Gaussian fitting, from which the principal and the auxiliary axes were extracted as 2.3σlong and 2.3σshort, respectively. The shape factor was calculated as a ratio between the auxiliary and the principal axes.
For the quantification of PSD95 cluster and nanodomain surface areas, the single molecule-based super-resolution images were first segmented by the morphometric image analysis module of MetaMorph software for structures >0.02 μm2. Morphologic features, such as area, length, and breadth of each cluster, were exported to calculate their respective distribution. The dimensions of PSD95 subclusters were calculated the same way as AMPAR nanodomains.
Proteolytic cleavage assay of extracellular Eos.
The proportion Eos::TEV::GluA1 that was localized on neuronal surface was assessed using AcTEV Protease (Invitrogen) to cleave the Eos tag from surface expressed receptors. sptPALM was used to measure fluorescence intensity of Eos::TEV::GluA1-expressing cells before and after incubation with AcTEV protease in extracellular medium (1:60 dilution in presence of 1 mm DTT, for 4–5 min, 37°C). Control neurons were incubated in AcTEV protease, which was rendered inactive by boiling for 10 min.
Live cell nanogold immunolabeling and electron microscopy of cultured neurons.
Coverslips with attached neurons (DIV 18–21) were incubated for 3 min at 37°C in GluA2 monoclonal antibody diluted in culture medium (6 μg/ml, gift from E. Gouaux). Neurons were removed from antibody and fixed in freshly prepared, prewarmed EM-grade 4% paraformaldehyde (EMS) in 0.15 m Sorensen's phosphate buffer (EMS) at room temperature for 45 min. All subsequent steps were performed at room temperature. Neurons were rinsed 3 times in 0.15 m Sorensen's PB, once in 0.1 m Millonig's PBS, and then blocked in 0.1 m Millonig's PBS with 2% BSA and 0.1% cold water fish skin gelatin (Aurion, EMS) for 60 min. Next, neurons were incubated in FluoroNanogold anti mouse Fab′ AlexaFluor-488 (Nanoprobes) diluted 1:100 in 0.1 m Millonig's PBS blocking solution for 90 min, then rinsed once in blocking solution, once in Sorensen's PB, and placed in freshly prepared 2% glutaraldehyde (EMS) in Sorensen's PB for 30 min. Afterward, neurons were stored in Sorensen's PB until silver intensification.
Processing and embedding for electron microscopy.
FluoroNanogold label was enhanced for 6–8 min using HQ Silver Reagent (Nanoprobes) according to manufacturer's instructions and processed immediately for EM; all steps were performed at room temperature. After several rinses in Sorensen's PB, neurons were incubated in 0.2% OsO4 in Sorensen's PB for 30 min, and then rinsed 10 times in dH2O to remove all traces of PB before placing neurons in filtered 0.25% uranyl acetate dissolved in dH2O for 30 min. After several water rinses, neurons were dehydrated by 3 min incubations in a graded series of ethanol: 50%, 70%, 95%, and 100%. No propylene oxide was used to prevent loss of immunogold label. Samples were infiltrated during 1–2 h steps in 70% Epon812/ethanol mixture followed by 2 exchanges of 100% freshly prepared Epon812 (Taab Laboratory and Microscopy), and finally embedded in freshly prepared Epon812. To allow cutting of en face sections of neurons, coverslips were placed cell side facing up on a glass slide and gelatin capsules filled with Epon812 were inverted and placed on top of coverslip, and polymerized at 60°C for 48 h. Coverslips were removed from polymerized samples by gentle heating over a flame while pulling slightly on the glass slide and gelatin capsules. Ultrathin sections (60–70 nm) were cut using an Ultra 35°diamond knife (Diatome) and a Leica Ultracut UCT M26 (Leica Microsystems) and picked up on 2 mm slot grids with a 1% formvar support film (EMS). When cutting serial sections, a 50:50 mix of rubber cement/xylene was applied to the top and bottom sides of the trimmed block to prevent sections from separating. Sections were contrasted with 3% aqueous uranyl acetate for 5 min, and then Reynolds's lead citrate for 7 min before imaging using a Hitachi H7650 transmission electron microscope operated at 80 kV. Images were captured using an Orius CCD (Gatan).
Analysis of the pattern and distribution of immunogold labeling on 2D EM sections.
Micrographs were analyzed using MetaMorph software (Molecular Devices). To quantify the density of immunolabel on different membrane domains, ROIs were drawn along membranes that were categorized as follows: spine PSD membrane, extrasynaptic spine membrane, or dendritic shaft membrane. For each ROI, the surface was measured and the number of immunolabels that were within 20 nm of the membrane was counted. The size of the clusters was measured using the Integrated Morphometry Analysis module, after adjusting the threshold to select label from the background. Immunolabel particles that were overlapping or separated by <20 nm were identified as a cluster.
3D-reconstruction of synaptic spines by serial section EM.
Series of images taken at 40,000–50,000 magnification from 6 to 8 serial sections of 70 nm thickness were used to generate reconstructions of individual immunogold labeled spines. Stacks of images were aligned using the MetaMorph software auto-align module. Reconstructions were generated from the aligned stacks using Imaris software (Bitplane). The membrane contours of the spine were traced in each section to create the spine volume. Immunolabel was reconstructed using the Spots tool. The number and size of nanodomains per spine were measured on each image of the serial sections as previously described. Label that appeared in the same location in successive serial sections was considered to be the same immunolabel and counted only once.
Quantification of AMPAR content in nanodomains.
The number of AMPARs per nanodomain was estimated from dSTORM super-resolution reconstructions of endogenous GluA2 using a custom-made analysis module operating inside MetaMorph software. For each cell, single AMPARs were identified using wavelet segmentation and Gaussian fitting as isotropic and isolated objects of size <40 nm in diameter. Isolated single AMPARs could be differentiated from nanoclusters resulting from the high spatial resolution provided by the dSTORM imaging technique. Then, the histogram and median of the integrated intensity of each individual AMPAR per cell were computed. This corresponds to the intensity distribution of one tetrameric structure, in function of the ratio of immunolabeled GluA2 with the other nonlabeled AMPAR subunits. Finally, synaptic nanodomains were detected using wavelet segmentation similarly to the cluster analysis described in Materials and Methods, and the number of receptors per cluster was estimated by dividing the cluster's total intensity by the median intensity of the identified isolated AMPARs.
Electrophysiology.
The extracellular recording solution contained in mm concentration as follows: 130 NaCl, 2.5 KCl, 2 MgCl2, 2 CaCl2, 10 HEPES, and 10 d-glucose, pH 7.4. To block GABA-A receptors, 50 μm picrotoxin was added to solution, and mini-EPSC recordings were performed in the presence of 1 μm TTX. The bath temperature was kept at 34°C. Patch-clamp microelectrodes were pulled for a 4–6 MΩ resistance, from borosilicate on a P-97 model puller (Sutter Instruments). Mini recordings were performed with an EPC 10 double patch-clamp amplifier (HEKA Elektronik). Data analysis was performed with a homemade software (Penn et al., 2012). Cells are clamped at −70 mV.
Statistics.
Statistical values are given as medians ± 25%/75% interval or mean ± SEM, except for Figure 10 where mean ± SD is used. Statistical significances were performed using Sigmaplot software (Systat Software). Non-Gaussian distribution datasets were tested by Mann–Whitney U test, and Gaussian distributions were tested by t test.
Calculation of required sample size.
The experiments were designed to compare two conditions per sample, including a control condition. The experimental observations of interest were changes in the localization and mobility of the nanoscale organization of synaptic molecules on primary cultured hippocampal neurons. In these experimental conditions, the sample size was calculated with a power factor of 0.6–0.8 and α of 0.2–0.5, which required sample sizes of 5–12 cells per condition depending upon the SD of the sample. To account for variability between cultures, samples were chosen over 2–4 independent cultures. These values were obtained from the power and sample size calculator from statistical solutions. (Source: http://www.statisticalsolutions.net/pss_calc.php).
Modeling.
AMPAR kinetics were modeled using Monte Carlo simulation essentially as described previously (Glavinovic and Rabie, 1998; Franks et al., 2002; Raghavachari and Lisman, 2004; Heine et al., 2008) using a 7 state model (see Fig. 10) and published values for kinetics rate constants, adjusted to fit the experimentally recorded miniature EPSCs (mEPSCs) in our system. The basic model considered a 300 × 300 nm postsynaptic membrane with 400 receptors disposed at varying distances. Glutamate (3000 molecules) was released in the presynaptic cell lying 15 nm away from the postsynapse. Glutamate was allowed to diffuse at 0.1 μm2/ms (Nielsen et al., 2004; Budisantoso et al., 2012) as in the synaptic cleft (15 nm wide). Monte Carlo simulations were performed with a ΔT = 1 μs time step. For the kinetic constants, the values were modified from Heine et al. (2008): C1 to C2, K1 = 4,590,000 m−1 s−1; C2 to C1, K1r = 4260 s−1; C2 to C3, K2 = 28,400,000 m−1 s−1; C3 to C2, K2r = 3260 s−1; D7 to D6, K3 = 1,270,000 m−1 s−1; D6 to D7, K3r = 45.7 s−1; C3 to O4, α = 1000 s−1; C2 to D7, α1 = 2890 s−1; C3 to D6, α2 = 120 s−1; O4 to D5, α3 = 200 s−1; D6 to D5, α4 = 16.8 s−1; 04 to C3, β = 200 s−1; D7 to C2, β1 = 20 s−1; D6 to C3, β2 = 0.727 s−1; D5 to 04, β3 = 4 s−1; and D5 to D6, β4 = 190 s−1.
Results
Visualization of AMPAR subunit organization by super-resolution light microscopy and electron microscopy
We applied four super-resolution imaging techniques to identify AMPAR distribution and mobility at a resolution well below the diffraction limit of light microscopy. On living hippocampal neurons, we applied sptPALM (Manley et al., 2008) and uPAINT (Giannone et al., 2010) to obtain high-density maps of AMPAR localization and movement on neuronal surfaces. These techniques are based on the localization of individual fluorescent molecules, either by stochastic activation of isolated fluorescent molecules (PALM) or by stochastic labeling (uPAINT) of the molecule of interest to provide high-density localization of isolated single molecule events (∼1000 localizations/μm2/min) over neuronal surfaces at a spatial resolution ∼45 nm. Third, STED microscopy (Hell, 2009) was used as an additional bulk super-resolution imaging technique to investigate the nanoscale AMPAR localization in fixed hippocampal neurons. Fourth, dSTORM (Heilemann et al., 2008; van de Linde et al., 2011), also called Ground State Depletion followed by Individual Molecule return (Folling et al., 2008), was used as a dual-color single molecule-based super-resolution microscopy method to describe the respective nanoscale distribution of AMPARs and PSD95 at a resolution of ∼20 nm. Finally, we compared these super-resolution light microscopy data with AMPAR distribution observed by pre-embedding immunogold EM.
sptPALM revealed that AMPARs are confined in stable subsynaptic clusters
sptPALM experiments were performed on neurons transfected with GluA1 and GluA2 AMPAR subunits genetically fused with the variant of Eos photo-switchable protein, mEos2 (referred to as Eos::GluA1 and Eos::GluA2). Homer1c fused to the Cerulean fluorescent protein (Homer1c::Cerulean) was used as a synaptic marker because it labels postsynaptic sites and does not interfere with the mobility and localization of AMPARs (Bats et al., 2007).
Diffraction limited epifluorescence images of Eos::GluA1 displayed a relatively homogeneous distribution of GluA1 proteins along the dendrites and enrichment in the spine head (Fig. 1A; Epifluorescence). In contrast, sptPALM super-resolved intensity images obtained from living neurons presented a strongly nonuniform distribution of Eos::GluA1 in different subcellular regions (Fig. 1A; sptPALM-Intensity). Reconnecting the localization of individual Eos::GluA1 molecules between sequential images allowed us to compute the trajectories of individual GluA1 molecules along the neuronal surface (Fig. 1A; sptPALM-Trajectories) and to map the instantaneous diffusion coefficient distribution (Fig. 1A; sptPALM-Diffusion). We could observe several thousand trajectories in a single dendritic segment within 1–2 min. This coarse cartography revealed a high heterogeneity of receptor mobility between subcellular regions. To ensure the surface expression of Eos-tagged receptors, we incubated neurons expressing a GluA1 construct with a TEV cleavage site between Eos and the GluA1 N terminus (Eos::TEV::GluA1) with the active TEV protease. This resulted in a 90–95% reduction in single molecule fluorescence (Fig. 2A), indicating that most Eos::TEV::GluA1 fluorescence was on the cell surface.
We used sptPALM data to compute the localization and trajectories of AMPARs (Figs. 1B and 2E) as well as the distribution of receptor mobility inside synaptic sites labeled by Homer1c::Cerulean (Fig. 1C) and compared the differences in global mobility of Eos::GluA1, Eos::GluA2 (Fig. 2E), and surface native AMPARs (Fig. 2F). Overall, the mobility of overexpressed AMPAR and endogenous surface receptors was faster in the dendritic shafts than in synaptic areas (Fig. 2E), as previously shown (Bats et al., 2007). Close examination of super-resolved intensity images and trajectory maps of Eos::GluA1 at spines indicated that Eos::GluA1 is concentrated in several distinct clusters of sizes smaller than the diffraction limit (Fig. 1B; GluA1-Intensity map) and that over half of Eos::GluA1 moves <0.01 μm2/s inside synapses (Fig. 1B; GluA1-Trajectories; and Fig. 1C; Distribution). Within these clusters, the receptors undergo confined movements, as indicated by the plateau reached by the mean square displacement (MSD) curve (Fig. 1C; MSD plateau ∼0.006 μm2, median values of the diffusion coefficient D in μm2/s with interquartile range [IQR] for Eos::GluA1; dendrites 0.0153 IQR 0.00036–0.0687, synapses 0.0032 IQR 0.00001–0.0195). In contrast, expression of Eos::GluA2 yielded different images, the formation of synaptic clusters being largely reduced (Figure 1B; GluA2-Intensity map). Eos::GluA2 moved faster (Fig. 1C; GluA2 Distribution) and trajectories (Fig. 1B; GluA2-Trajectories) displayed higher diffusion and much less confinement (Fig. 1C; MSD, plateau ∼0.02 μm2, diffusion coefficient of Eos::GluA2 in μm2/s; dendrites 0.052 IQR 0.014–0.127, synapses 0.015 IQR 0.001–0.052). This indicates that Eos::GluA2 visits 4 times more synaptic surface than Eos::GluA1 in the same time. For Eos::GluA2-expressing cells, >90% of the spines did not exhibit synaptic clusters. Expression of individual subunits mostly forms homomers, whereas endogenous receptors are mostly formed of GluA1/GluA2 heteromers. Interestingly, the coexpression of Eos::GluA2 with nontagged GluA1 partially rescued the formation of subsynaptic clusters (Fig. 1B; GluA1/2-Intensity map). Cells coexpressing Eos::GluA2 and GluA1 displayed an increased confinement of Eos::GluA2 and a decreased diffusion (Fig. 1C; GluA1/2 MSD, plateau ∼0.013 μm2, diffusion coefficient of Eos::GluA2+GluA1 in μm2/s; dendrites 0.032 IQR 0.003–0.103, synapses 0.0138 IQR 0.0002–0.056). This was also consistent with changes in the global mobility of Eos::GluA2 where the fraction of Eos::GluA2 moved slower when coexpressed with untagged GluA1 (Fig. 2E).
Together, these data suggest that expressed GluA1 subunits are less mobile and are organized into subdiffraction sized clusters on hippocampal spines compared with ectopic GluA2 subunits. However, the subdiffraction limited organization and the slow moving fraction of receptors can be partially restored on ectopic GluA2 when coexpressing GluA2 with GluA1. This indicates the importance of GluA1 subunit in the subspine organization and retention of AMPARs on the synaptic membrane.
uPAINT revealed that AMPARs are confined in stable subsynaptic clusters
uPAINT was performed as a complementary high-density single particle tracking technique to monitor the localization and movement of endogenous AMPARs located on the neuronal surface plasma membrane. We used a monoclonal antibody directed against an extracellular region of the GluA2 AMPAR subunit conjugated to the ATTO647N fluorescent dye to randomly label surface AMPARs containing the GluA2 subunit on living hippocampal neurons (Fig. 1D; Surf-GluA2). This method specifically labeled endogenous GluA2 containing AMPARs as demonstrated by the loss of labeling observed in neurons cultured from GluA2 KO mice (Fig. 2B). Furthermore, it provided information on the general behavior of endogenous AMPARs because the majority of AMPARs are either GluA2/GluA1 or GluA2/GluA3 heterotetramers (Wenthold et al., 1996). In contrast to sptPALM, uPAINT is a technique based on specific stochastic labeling (Fig. 2C) that does not allow for the labeling of all available receptors on the surface. This drawback is counterbalanced by the advantage that no overexpression is needed and that each antibody is labeled with several organic fluorophores, providing a bright probe that can be localized with better accuracy and observed for longer durations than Eos molecules.
uPAINT of GluA2-containing AMPARs yielded a strongly nonuniform distribution of fluorescence intensity (Fig. 1D; Surf-GluA2-Intensity map) similar to that obtained during sptPALM images of Eos::GluA1-containing receptors. In each spine, AMPARs were concentrated in a few subdiffraction-sized clusters. Receptors were mainly slowly diffusing and highly confined in these domains (Figure 1D; Trajectories; and Fig. 1E; MSD; plateau ∼0.006 μm2, diffusion coefficient of Surf GluA2 in μm2/s; dendrites 0.0012 IQR 0.00001–0.007, synapses 0.0008 IQR 0.00001–0.0048). As previously reported (Bats et al., 2007; Opazo et al., 2010), endogenous GluA2-containing receptors moved slower than overexpressed subunits. Together, both sptPALM of expressed Eos::GluA1 homomers or Eos::GluA2/untagged GluA1 heteromers and uPAINT of endogenous GluA2-containing AMPARs yielded a comparable distribution of receptors confined in subsynaptic clusters.
STED confirmed that AMPARs are concentrated in subsynaptic clusters
The intensity localization maps of AMPARs reconstructed from sptPALM and uPAINT imaging experiments rely on the accumulation of multiple sequential single molecule localizations and may contain undersampling biases. We thus acquired images of AMPAR distribution using the alternative STED nanoscopic imaging approach (Hell, 2009). We live-stained neurons expressing the synaptic marker Homer1c::GFP with the same monoclonal antibody against an extracellular epitope of GluA2 used in uPAINT, followed by a secondary antibody coupled to ATTO647N. Postsynapses labeled with Homer1c::GFP were first identified by scanning confocal microscopy (Fig. 3A; Confocal Homer1c). STED images of the surface localized endogenous GluA2 were then acquired (Fig. 3A; STED GluA2). Consistent with sptPALM and uPAINT experiments, STED images confirmed that AMPARs are highly clustered into subdiffraction sized domains at synapses (Fig. 3B; sptPALM, uPAINT, STED).
Subsynaptic AMPAR clusters are nanometer-scale “nanodomains”
The super-resolution techniques used in this study enabled the subdiffraction size and shape of AMPAR clusters to be measured. Anisotropic bidimensional Gaussian fitting was applied (Fig. 3C) as described in Materials and Methods (clusters contained in n = 323 synapses by sptPALM, 126 by uPAINT, 118 by STED). The long (principal) and short (auxiliary) axes were determined as the full width at half-maximum (FWHM) of the fitted Gaussian. Most of the principal axes were measured to be smaller than the diffraction limit, with a median of 69.6 nm (IQR 54–93 nm) by sptPALM, 78.1 nm (IQR 64.9–87.1 nm) by uPAINT, and 132 nm (IQR 98–191 nm) by STED (Fig. 3D). These values are in accordance with the small size of confinement obtained from MSD curves (Fig. 1C). Because of their nanoscale dimensions, we henceforth refer to these AMPAR clusters detected by super-resolution microscopy as AMPAR nanodomains.
Most AMPAR nanodomains exhibited an ovoid shape with an average shape factor of 0.6–0.7 (Fig. 3E). As a control, we performed similar nanodomain quantifications on fixed samples. The principal axes of AMPAR nanodomains in fixed samples were distributed around a median of 78.1 nm (IQR 54–95.7 nm) and 62.8 nm (IQR 48.9–80.3 nm) using sptPALM and uPAINT, respectively, similar to that obtained in live-cell experiments (Fig. 3E). The confinement area of AMPAR movements measured by uPAINT experiments on fixed neurons (0.0028 μm2) is almost twofold smaller than the average size of the nanodomains (0.005 μm2) (Fig. 2D), indicating that the determination of the nanodomain size is not limited by the pointing accuracy. Similarly, the resolution provided by sptPALM under our experimental conditions was better than 50 nm (see Materials and Methods). The larger nanodomain size measured by STED likely reflects its poorer resolution in our acquisition conditions.
Finally, we quantified the fraction of AMPARs in nanodomains and the number of nanodomains per spine. All imaging techniques showed that ∼65% of the receptors were concentrated in nanodomains and that >80% of spines contained at least one nanodomain (88% by sptPALM, 82% by uPAINT, and 85% by STED) (Fig. 3F). The average number of nanodomains per spine was ∼2.5 (2.42 ± 0.18, n = 584 spines) when measured with sptPALM. STED yielded similar results (2.38 ± 0.17, n = 351 spines), whereas uPAINT led to a lower average of 1.4 ± 0.1 (n = 299 spines), probably because of the intrinsic underlabeling of the technique (Fig. 3F). Spines without nanodomains usually exhibited a low level of diffuse AMPAR staining that was not further quantified.
AMPAR nanodomains detected by live-labeling and pre-embedding immunogold electron microscopy
Despite the advent of super-resolution optical microscopy, immunogold labeling by EM still provides the highest resolution for AMPAR localization and remains the “gold standard” (Kharazia et al., 1996; Bernard et al., 1997; Kharazia and Weinberg, 1997; Nusser et al., 1998; Takumi et al., 1999b; Ganeshina et al., 2004), with the important drawback that it can only image fixed tissue. Most previous EM studies have used the postembedding immunolabeling method, which provides incomplete information about the distributions of receptors because of its relatively low sensitivity (Nusser et al., 1998; Masugi-Tokita and Shigemoto, 2007). Pre-embedding immunogold labeling (Bernard et al., 1997) allows mapping the distribution of surface AMPARs with increased sensitivity and has been recently used to observe the distribution of GluA2, the most prevalent subunit in AMPARs (Tao-Cheng et al., 2011). Here, we performed similar pre-embedding immunogold labeling, with two deviations. First, we immunolabeled surface GluA2-containing AMPARs in live neurons in contrast to the common postfixation protocols. Second, instead of the commercial 6C4 anti-GluA2 clone, we used a GluA2 antibody with excellent avidity, as confirmed by uPAINT experiments (Fig. 2B) (Giannone et al., 2010). After 3 min of live-labeling, neurons were fixed, labeled with a secondary fab′ fragment coupled to nanogold (1.4 nm size), and processed for EM including silver enhancement of nanogold as indicated in Materials and Methods.
A typical view of surface labeling is shown in Figure 4A. The density of immunogold labeling at the PSD membrane in spines was 7.23 particles/μm. This was ∼3.5 times more dense than labeling on extrasynaptic spine membrane (2.12 particles/μm) and 5 times more concentrated than along the dendritic shaft membrane (1.24 particles/ μm) (n = 538 label particles measured on 31.94 μm of PSD membrane, 74.59 μm of extrasynaptic spine membrane and 120.65 μm of dendritic shaft membrane). Interestingly, immunolabel was often organized in clusters located in the PSD as well as on extrasynaptic membranes (Fig. 4B). We measured the size of the immunolabel clusters (defined as 2 or more particles separated by <20 nm) (Fig. 4C) and found that the average cluster length was 77.48 ± 23.98 nm (n = 123 clusters). This size is in close agreement with the size of nanodomains measured by the super-resolution light microscopy techniques. The size of particle clusters did not vary between membrane regions (PSD clusters = 81.48 ± 27.12 nm; spine clusters = 72.72 ± 20.04 nm; dendritic shaft clusters = 77.04 ± 23.32 nm), and each cluster contained on average ∼3 label particles. Similar to findings at the super-resolution level, ∼60% of the label was contained in a cluster, a value that was consistent between membrane regions (% of particles in cluster: PSD 55%, spine 60%, dendritic shaft 66%). Based on these similarities, we concluded that the clustered immunolabel viewed at the EM level corresponds to nanodomains detected using super-resolution light microscopy.
To count the number of nanodomains present on individual spines and specify their location in the PSD or extrasynaptic spine membrane, seven spines were reconstructed by serial section EM (Fig. 4D). On average, there were 2.14 ± 0.90 nanodomains in the PSD and 4.14 ± 1.86 nanodomains on the spine extrasynaptic membrane.
Evolution and stability of AMPAR nanodomains in spines
We investigated the dynamics of AMPARs between nanodomains by uPAINT because this approach provides the longest observation time of individual fluorescent probes (median recording time for a single anti-GluA2-ATTO647N conjugate = 150 ms, IQR 100/300, n = 15000, with >10% of trajectories >0.5 s). Single receptor trajectories could be classified into three distinct patterns (Fig. 5A). In the first pattern, molecules diffused in and out of nanodomains with random variations in the instantaneous diffusion coefficient that were uncorrelated to a specific identified localization of the protein. This suggested that the receptors did not engage in stabilizing interactions (Fig. 5A,B; Type 1). In the second pattern, receptors alternated between nanodomains and the rest of the spine space, displaying clear periods of confinement when colocalized with the nanodomains (Fig. 5A,B; Type 2). This was particularly apparent on the plot of the instantaneous diffusion coefficient versus time of these molecules that displayed a drastically reduced diffusion coefficient once the molecule localized in the nanodomain. Finally, a third group of molecules were observed to remain confined in the nanodomains with very low diffusion coefficients for several minutes, indicating a strong association with molecular components of the nanodomain (Fig. 5A,B; Type 3). The much stronger confinement of receptors in nanodomains was also apparent when comparing the average MSD for the type 1 and type 3 mobility modes (Fig. 5A; MSD). These plots also indicate that, in between nanodomains, the receptors undergo rapid free diffusion. Based on this observation, synapses were divided into two functional zones discriminated by receptor mobility: zones inside nanodomains and zones outside nanodomains. Interestingly, the residence time of receptors outside nanodomains was short while they randomly moved rapidly (median residence in extra-nanodomain = 0.25 s, IQR 0.1/0.7, n = 161; median total trajectory length 0.75 s, IQR 0.55/1.4, n = 5000). In contrast, molecules were confined inside nanodomains for much longer durations, going up to few minutes (median residence inside nanodomains = 2.25 s, IQR 0.1/16.2, n = 155; Fig. 5C). Together, these data indicate that, inside synapses, AMPARs diffuse freely between nanodomains then rapidly hop inside nanodomains where they are strongly immobilized and reside for long periods. These stabilization domains could be related to the attracting molecular interaction potential wells we recently found by mathematical analysis performed on sptPALM data (Hoze et al., 2012).
Then, we investigated the temporal stability of nanodomains using time-lapse sptPALM. We performed 50 s periods of sptPALM acquisitions every 5 min over a time period of 45 min. We could identify and track individual nanodomains in the spine for several tens of minutes (Fig. 6A,B). As quantified in Figure 5D, >20% of the spines (n = 75) had a stable nanodomain during the entire acquisition time (45 min). Approximately 40% of the spines contained nanodomains, which could not be detected for longer than 5 min of observation, indicating that a fraction of the nanodomains in the spine are either dynamic or transient, whereas others have a very stable localization and lifetime. Nanodomains in the remaining 60% of the spines were stable for >5 min inside the synapses. More than 50% of spines had multiple nanodomains, bearing combinations of stable and nonstable nanodomains in the same synapse (Fig. 5D). Kymographs performed on highly stable nanodomains illustrate that nanodomains are visited by many molecules (Fig. 6A,B).
Finally, we analyzed the differences of mobility between in and out nanodomains, from single-molecule trajectories of Eos::GluA1 molecules obtained by sptPALM. Most of the immobilized molecules resided within the nanodomains, whereas molecules that were either outside or moving in and out of the nanodomains were either diffusive or weakly confined (Fig. 6C). We categorized the trajectories into three families, based on the nanoscale spatial heterogeneity as follows: (1) exclusively inside the nanodomain (In), (2) exclusively outside nanodomain (Out), and (3) those who visited the nanodomains for a minimum of 10% and maximum of 80% of the duration of their observed trajectories (In/Out). The cumulative distribution of the instantaneous diffusion of the trajectories displayed a shift toward lower mobility within, and higher mobility outside the nanodomains (Fig. 6D). The instantaneous diffusion distributions of In/Out molecules indicated an intermediate shift of mobility compared with In and Out ones (Fig. 6D). We also plotted the average of mean square diffusion curves for each of the three families of receptors (In, Out, and In/Out). The mean square diffusion of the molecules within nanodomains (In) displayed a strong confinement indicating immobilization (Fig. 6E). Receptors outside nanodomains (Out) exhibited unhindered Brownian diffusion. The mean square diffusion of receptors temporarily associated with nanodomains (In/Out) displayed weak confinement, indicating the transient changes in the mobility between strong confinement and unhindered diffusion. These data confirm that nanodomains are indeed zones of strong confinement where AMPA receptors are transiently trapped with the potential for receptor to exchange and diffuse freely into the membrane outside the nanodomain.
Nanoscale co-organization of PSD95 and AMPARs
In the previous sections, the synaptic localization of nanodomains by super-resolution light microscopy was achieved with the epifluorescence image of Homer1c::Cerulean. Even though Homer1c is a postsynaptic protein, recent results indicate that it might not be closely associated with the PSD (Dani et al., 2010). To monitor the distribution of AMPAR nanodomains within the PSD, we used PSD95, a key organizer for postsynaptic scaffolding. PSD95 is also involved in the positioning and synaptic retention of glutamate receptors within the PSD and is localized closer to the postsynaptic membrane than Homer1c (Dani et al., 2010). We investigated the organization of PSD95 in spines by performing sptPALM on neurons expressing PSD95 genetically fused with mEos2 (Eos::PSD95). Epifluorescence images showed that the majority of Eos::PSD95 was localized to spine heads (Fig. 7A; Epifluorescence). Corresponding images obtained using sptPALM displayed similarly strong Eos::PSD95 labeling in spine heads with >80% of the molecules observed by sptPALM located in spines and sparse labeling outside (Fig. 7A; sptPALM-Intensity). At the spine level, most of the Eos::PSD95 signal was concentrated into a single large zone or, more rarely, in multiple zones (Fig. 7B). The area of Eos::PSD95 clusters was computed using a combination of thresholding and morphometry analysis on sptPALM intensity images. The size of Eos::PSD95 clusters (n = 317 synapses) followed a Gaussian distribution centered on 0.0405 ± 0.025 μm2 (Fig. 7C), with the mean of the principal and auxiliary axis calculated to be ∼424.7 ± 177.2 nm (Fig. 7C, inset) and 282.8 ± 78.85 nm, respectively. These measurements are consistent with previous reports on size of the PSD (Petersen et al., 2003). Thus, we concluded that Eos::PSD95 is a reliable marker of the PSD.
A line-scan through the Eos::PSD95-labeled PSDs indicated that within PSDs, PSD95 is heterogeneously distributed (data not shown, but see Fig. 7J for similar results obtained with dSTORM on endogenous PSD95). At the subdiffraction limit scale, PSD95 intensity distribution fluctuated from a basal level to several higher intensity peaks, suggesting organization into subdomains. These subdomains were measured by anisotropic bidimensional Gaussian fitting (Fig. 7D) as described in Materials and Methods. The principal and auxiliary axes of PSD95 subdomains (n = 465) were, respectively, 158.2 ± 100 nm (Fig. 7E) and 102.2 ± 39 nm. For comparison, the distribution of AMPAR nanodomain sizes, which were smaller, is also plotted in Figure 7E (solid line). Finally, the number of PSD95 subdomains was counted for each of the Eos::PSD95-labeled PSDs, revealing that >90% of PSDs contained one or more subdomains (Fig. 7E, inset). We also took advantage of sptPALM experiments to quantify the mobility of Eos::PSD95. As shown previously (Sharma et al., 2006), PSD95 showed low mobility in synapses and the distribution of PSD95 mobility within synapses was lower than the synaptic mobility of AMPARs (Fig. 7F).
To compare the distribution of PSD95 with respect to that of AMPARs at the PSD, we performed two-color dSTORM imaging (Heilemann et al., 2005; van de Linde et al., 2011) on neurons stained with antibodies for endogenous surface GluA1 coupled to Alexa-647 and for endogenous PSD95 to Rhodamine Red. As with the other super-resolution light imaging techniques, AMPAR nanodomains of comparable size and frequency were detected by dSTORM (Fig. 7H; nanodomains size: 87 ± 6 nm and 1.2 ± 1.0 nanodomains per spine, n = 87 spines). Like Eos::PSD95, endogenous PSD95 exhibited a heterogeneous distribution inside the PSD with higher concentrations in subdomains (Fig. 7I,J). Although AMPARs were strictly localized to nanodomains and had low basal levels in the regions outside the PSD and outside nanodomains within the PSD, the basal level of PSD95 was higher between subdomains in the PSD than outside the PSD.
When viewed at the epifluorescence level, PSD95 and GluA1 displayed a high degree of overlap (Fig. 7G). However, dSTORM images revealed a varying degree of overlap between PSD95 and GluA1 (Fig. 7H,I). A detailed analysis of a series of individual PSDs (Fig. 7I,J) indicated that PSD95 and AMPARs exhibited only partial colocalization. In many instances, the AMPAR nanodomains overlapped with the PSD95 subdomains, indicating an organizational relationship between the two proteins at the synapse (Fig. 7J). However, this perfect colocalization at the nanoscale is not systematic; AMPAR nanodomains and PSD95 subdomains were often tens of nanometers away from each other. This illustrates differences in the synaptic organization of AMPARs and PSD95. This result is at variance with a recent study that reported similar nanodomains than in our work but found a high degree of colocalization between AMPAR and PSD95 (Macgillavry et al., 2013). This may arise from the fact that this study used overexpressed PSD95 or AMPARs, manipulations that may alter the normal reciprocal distribution of both elements.
Next, we correlated the endogenous PSD area labeled by PSD95 with the number of AMPAR nanodomains from the two-color dSTORM super-resolution images. PSDs were grouped into three different size categories: <0.05, 0.05–0.1, and >0.1 μm2. The number of nanodomains per group was 0.24 ± 0.04, 0.72 ± 0.06, and 1.08 ± 0.78, respectively. These data indicated that small PSDs were more likely to contain less or no nanodomains compared with larger ones (Fig. 8A). Interestingly, the number of nanodomains per PSD increased linearly with the PSD area (correlation coefficient = 0.997, SD = 0.107) (Fig. 8B). Finally, we measured the distance between centroids of nanodomains contained within the same spine or the same PSD. The average distance between the nanodomains contained within a spine was 512.88 ± 26.5 nm (n = 91 spines). In <30% of the cases, the distance between nanodomains within a spine was <250 nm (median 455 nm, IQR 309.11/732.7, n = 91). When looking specifically at nanodomains within the same PSD95 cluster, the average distance between nanodomains was 317.30 ± 168.34 nm, and the average PSD diameter was 339.14 ± 228.89 (18 PSDs). In ∼50% of the cases, the distance between nanodomains in the same PSD was <250 nm (median 247.6 nm, IQR 191.1–439.7 nm), confirming that nanodomains are more concentrated in PSDs than overall in the spine. The average size of the PSDs obtained by endogenous PSD95 staining was smaller than that obtained by overexpression of mEos::PSD95. This is consistent with the observation that PSD95 overexpression can increase the size of the dendritic spines (El-Husseini et al., 2000).
Modulation of endogenous PSD95 levels and its effect on AMPAR nanodomains
PSD95 is one of the core components of the PSD and could thus play a role in organizing AMPAR nanodomains. To test this hypothesis, we modulated the level of PSD95 in neurons using shRNA to reduce PSD95 levels or overexpression of PSD95. An shRNA-resistant PSD95 rescue condition was included to verify the specificity of shRNA effects. These conditions will be referred to as SH::PSD95, OverEx::PSD95, and Rescue::PSD95. Immunofluorescence staining for PSD95 was used to measure variations in PSD95 expression level. In cells expressing SH::PSD95, we observed a 20–25% reduction in the fluorescence intensity of endogenous PSD95 compared with nontransfected cells. In neurons expressing Rescue::PSD95 and OverEx::PSD95, we observed a 6- to 7-fold and 15- to 17-fold, respectively, increase in PSD95 levels compared with nontransfected cells. These results were consistent with previous reports using these DNA constructs (Schluter et al., 2006; Mondin et al., 2011).
To quantify the effect of PSD95 knockdown or overexpression on AMPAR nanodomains, we performed dSTORM imaging of immunostained surface GluA2 (endogenous) in neurons expressing the various constructs. The average length of the nanodomains significantly increased from control neurons compared with neurons expressing SH::PSD95, OverEx::PSD95, and Rescue::PSD95 (Table 1; Fig. 8C; Mann–Whitney test, p = 0.017, p = 0.039, and p = 0.015, respectively). Although the Rescue::PSD95 resulted in larger nanodomains (116 ± 8 nm), this increase was found to be not significant compared with SH::PSD95 and Overex::PSD95. The number of nanodomains per spine was also significantly affected by the modulation of PSD95 levels (Table 1; Fig. 8D). Nanodomain number decreased in cells expressing SH::PSD95 compared with control ones (Table 1; Fig. 8D; control number of nanodomain/spine 1.21 ± 0.13, SH-PSD95 0.98 ± 0.11, p = 0.046). In contrast, the number of nanodomains per spine increased dramatically in cells expressing OverEx::PSD95 (Table 1). Although expression of the Rescue::PSD95 displayed an increase in the nanodomain density of 1.34 ± 0.12 compared with the control, this increase was not found to be significant (Table 1).
The coexpression of a generic scrambled shRNA against a noneffective mammalian protein with GFP (SH::Scramble::GFP) or with another synaptic protein Homer 1c::GFP (SH::Scramble::Homer1c) displayed no significant changes in the nanodomain characteristics compared with the control untransfected cells or cells expressing Homer1c::GFP (Table 1). No significant difference was detected either in the nanodomain size or number between OverEx::PSD95 and Rescue::PSD95, indicating that the Rescue::PSD95 mimics mild overexpression of PSD95. Consequently, the dynamic experiments on live cells and further quantifications were performed using control, SH::PSD95, and OverEx::PSD95 unless otherwise stated.
We also estimated the surface covered by nanodomains in spines. In control conditions, nanodomains covered 28 ± 2.08% of the spine surface and concentrated 65.7 ± 3.11% of AMPAR signal (n = 8 neurons). For the whole neurons, the nanodomains accounted for 43 ± 5.14% of the total AMPAR signal but represented only 2.09% of the total observed neuronal surface. Interestingly, on cells overexpressing PSD95, AMPARs were more clustered as nanodomain content increased to 55.1 ± 4.62% of the total AMPAR signal in only 1.27% of the neuronal surface. SH::PSD95 expression induced a reduction in clustering as nanodomain content decreased to 36 ± 3.52% of AMPARs signal, present on 0.98% of the total neuronal surface. This indicated that overexpression of PSD95 concentrated AMPARs in nanodomains, whereas reduction of PSD95 reduced the number of AMPARs on the neuronal membrane as well as in nanodomains.
Relationship between PSD95 expression level and association with stargazin on the dynamic organization of AMPARs at the synapse
The intracellular scaffold protein PSD95 has been suggested to play a major role in stabilizing AMPARs at synaptic sites through binding to their auxiliary TARP (Schnell et al., 2002; Bats et al., 2007; Opazo et al., 2012). To investigate the role of PSD95 in the dynamic organization of nanodomains, we quantified AMPAR dynamics in neurons with uPAINT while modulating PSD95 levels. We computed the distribution of GluA2 mobility globally (on the whole neuron) and selectively inside synapses (Fig. 8E) and compared it with GluA2 mobility in control neurons expressing Homer-DsRed or neurons expressing OverEx::PSD95 and SH::PSD95. When OverEx::PSD95 was expressed, the global mobility of AMPARs was slower than in control neurons. No significant difference in AMPAR mobility was observed at synapses, probably because of the already saturated concentration of PSD95 (global median values Homer-DsRed, n = 18 cells, 0.0012 IQR 0.00001–0.007, synapses 0.0008 IQR 0.00001–0.0048; Homer-DsRed + OverEx::PSD95, n = 14 cells, 0.0009 IQR 0.00001–0.0052, synapses 0.00076 IQR 0.00001–0.005). Expression of SH::PSD95 induced the opposite effect. Both the global and synaptic mobility of AMPARs was significantly increased (Homer-DsRed + SH::PSD95, n = 15 cells, 0.0034 IQR 0.00006–0.021, synapses 0.0023 IQR 0.00001–0.011).
We next analyzed variations in AMPAR nanodomain size in spines as a function of PSD95 expression level using uPAINT (Fig. 8F). The average length of the nanodomain principal axis increased from 74 ± 1 nm for control spines to 164 ± 5 nm for OverEx::PSD95 (n = 311) and 155 ± 4 nm for SH::PSD95 (n = 229) (Fig. 8F, inset). These experiments confirm the quantifications obtained from dSTORM data and show that both increase and decrease in PSD95 expression levels trigger a significant increase in AMPAR nanodomain size.
We demonstrated previously that overexpression of a mutant form of stargazin (StgT231F) (Schnell et al., 2002; Bats et al., 2007), which prevents C-terminal association of stargazin to PSD95, increases the mobility of AMPARs globally and at synapses (0.04 IQR 0.0036–0.15 synapses 0.035 IQR 0.006–0.12), demonstrating that the interaction between stargazin and PSD95 is critical for AMPAR immobilization. Here, we examined the effect of StgT231F expression on the density and size of nanodomains (Table 2). StgT231F induced a major reduction in nanodomain density compared with the control (control 0.13 ± 0.03 nanodom · μm−2 and StgT231F 0.007 ± 0.002, p < 0.001). On the contrary, coexpression of StgT231F with a rescue variant of PSD95 (PSD95H225V), which can bind StgT231F (Schnell et al., 2002), partially rescued global nanodomain density and AMPAR immobilization (0.03 IQR 0.0032–0.14 synapses 0.028 IQR 0.0039–0.11). The coexpression of StgT231F and PSD95H225V also resulted in a large increase in the global density of nanodomains compared with the StgT231F expression alone (coexpression 0.03 ± 0.01 compared with Stg T231F 0.007 ± 0.002, p = 0.042). The expression of PSD95H225V alone did not significantly affect the mobility of AMPA receptors globally or at synapses (0.007 IQR 0.00003–0.072 synapses 0.00358 IQR 0.00002–0.032). It also only had a minor effect on the global nanodomain density of AMPARs (Table 2). This is expected as only one of the PDZ domains in PSD95 is mutated while the other two PDZ domains can still bind to endogenous stargazin. On the contrary, StgT231F cannot bind to any of the PSD95 family PDZ domains, explaining its drastic effect on AMPAR organization.
Expression of StgT231F alone, PSD95H22V alone, or coexpression of StgT231F and PSD95H225V also increased the size of nanodomains (Table 2; p < 0.001 for the three conditions compared with the control), similarly to the effect of modulating endogenous PSD95 levels, confirming that the stargazin-PSD95 interaction is vital for the formation of nanodomains as well as AMPAR immobilization.
Correlating AMPAR organization and synaptic transmission
Next, we sought to estimate the number of AMPARs per nanodomains from dSTORM super-resolution images obtained from surface GluA2 staining. As a reference for the signal provided by an individual labeled AMPAR, we measured the distribution and median of the fluorescence intensity arising from isolated single emitting species in the shaft membrane, as detailed in Materials and Methods (Fig. 9A, red circles). The signal variability includes fluctuations that depend on the antibody labeling efficacy and fluorophore blinking properties. The intensity distribution of these single emitters fit perfectly with a bimodal Gaussian distribution (Fig. 9B) with the second peak mean intensity located at twice the first peak mean intensity (239.0 ± 138.1 and 118.4 ± 62.11, respectively, 9 cells, 4405 isolated AMPARs). This is coherent with the hypothesis that the signal from the isolated single emitting species arises from individual AMPARs that usually contain two GluA2 subunits and in which either one or both the GluA2 subunits could be labeled. Interestingly, the cell-to-cell fluctuations of the median intensity was <10%. To estimate the AMPAR content per nanodomain, we used the median intensity value of this single emitter distribution for each cell, as a putative reference for single receptor fluorescence intensity. Figure 9A illustrates various single receptors (red circles) and nanodomains (green regions) with their respective number of AMPAR estimation. We extracted the distribution of AMPARs per nanodomain for the control and the SH::PSD95 conditions (Fig. 9C). We observed a decrease of the number of AMPARs per nanodomain for the SH::PSD95 (21.6 ± 0.58, n = 12 cells) compared with the control conditions (27.3 ± 0.61, n = 9 cells), corresponding to a reduction of 21% (Fig. 9D). Interestingly, no significant difference was observed in the number of AMPARs per nanodomain between control and OverEx::PSD95 (27.76 ± 0.76, n = 14 cells) (Fig. 9D).
To correlate these fluctuations in nanodomain AMPAR content to synaptic responses, we measured the mEPSC amplitudes from neurons in similar conditions, and compared measurements in control, upregulated, or downregulated PSD95 expression level. For correlation purposes, we recorded cells from the same primary hippocampal cultures that were used for uPAINT experiments to quantify AMPAR dynamics and nanodomain size (Fig. 8C,D). The mEPSC amplitude distribution displayed a shift toward higher values when PSD95 was overexpressed compared with the control condition (Fig. 9E). mEPSC amplitude average values (Fig. 9F) revealed that OverEx::PSD95 is not significantly different from control (control, n = 17 cells, 19.7 ± 1.3 pA; OverEx::PSD95, n = 21 cells, 21.7 ± 0.7 pA). In contrast, PSD95 knock-down induced a 20% decrease in mEPSC amplitude (SH::PSD95, n = 14 cells, 15.8 ± 1.5 pA). Together, these experiments demonstrate a positive correlation between nanodomain AMPAR content and mEPSC responses.
Modeling the impact of AMPAR organization on synaptic transmission
To estimate the potential impact of AMPAR nanoscale organization on synaptic transmission, we performed Monte Carlo simulations as in Heine et al. (2008), but using various AMPAR organizations (Fig. 10). We used previously established parameters of glutamate release in the synaptic cleft, AMPAR activation, and desensitization kinetics schemes (Fig. 10A) (Jonas et al., 1993; Franks et al., 2002; Lisman et al., 2007; Heine et al., 2008). The model perfectly mimicked the kinetics and amplitude of experimentally recorded mEPSCs (Fig. 10C). We then systematically varied the localization (Fig. 10D–F) and density (Fig. 10G–I) of an AMPAR nanodomain containing 25 receptors and also analyzed the crosstalk between multiple nanodomains located at varying distances (Fig. 10J,K). All parameters strongly impacted the amplitude of mEPSCs. However, our simulations indicate that there is a degree of tolerance on the location of the glutamate release site with respect to the AMPAR nanodomain because the amplitude of the mEPSC only starts to decrease when glutamate is released >100 nm away from the nanodomain center (Fig. 10D–F). In contrast, the mEPSC amplitude was very sensitive to AMPAR density because it started to decrease as soon as it was lower than the maximal packaging we used of one receptor every 20 nm (Fig. 10G–I).
We then modeled the presence of two similar nanodomains in the synapse. First, glutamate was released on top of one cluster and the second one positioned at varying distances. The second cluster had no more impact on the AMPAR response when it was distant of >300 nm (Fig. 10 J1,2), which is about the mean value that we measured for intercluster distances. Second, we fixed the intercluster distance at 300 nm and varied the location of glutamate release. Interestingly, the AMPAR response was markedly increased when glutamate was released in between the two clusters (maximum 19.2 ± 3.6 activated AMPARs) compared with the condition with a single cluster in the synapse (maximum 16.2 ± 3.3 activated AMPARs, p < 0.002) (Fig. 10K1,2).
Together, our experimental data and these simulations indicate that AMPAR nanoscale organization as well as location of the glutamate release strongly impact mEPSC amplitude.
Discussion
Using four independent super-resolution light imaging methods on both genetically tagged and endogenous receptors, as well as pre-embedding EM, we demonstrate that, in live and fixed hippocampal neurons, AMPARs are organized inside most synapses into a few highly concentrated domains of ∼70 nm, and not homogeneously distributed as previously described. AMPARs are stabilized in nanodomains and freely randomly diffuse between them. Nanodomains are dynamic in their composition and overall structure. AMPARs can enter and exit nanodomains within seconds by lateral diffusion. Inside synapses, nanodomains can have two different behaviors: they can be stable for up to an hour or change shape and move and/or disappear within minutes, albeit at rather slow rates. Nanodomain estimated size is robust as measures of AMPAR distribution by sptPALM and uPAINT using both localization maps and confinement of receptor movements from MSD curves yielded comparable values. This was confirmed on fixed cells using STED, dSTORM, and EM. The observation of similar clusters in fixed and living cells reinforces that they are functional receptor domains. Live-cell single molecule imaging demonstrates that nanodomains are dynamic as AMPARs could enter and exit them frequently. Previous low-density single particle tracking data showed the alternating trapping and free diffusion of AMPARs in the PSD (Groc et al., 2004), as well as confinement in subspine compartments (Ehlers et al., 2007). However, these data could not provide enough spatial information to map nanodomains inside synapses. A recent super-resolution light microscopy study using immunocytochemistry and STORM on fixed samples (Dani et al., 2010) has highlighted the molecular organization of many postsynaptic molecules, including AMPARs. Even though an inhomogeneous AMPAR distribution was observed, with clusters apparent on the micrographs, the heterogeneity of receptor staining within the synapses appeared lower compared with our study and the existence of nanodomains was not reported.
Studies of AMPAR distribution at the nanometer level using EM and either postembedding immunogold labeling (Kharazia et al., 1996; Petralia et al., 1997; Nusser et al., 1998; Takumi et al., 1999a) or pre-embedding labeling (Tao-Cheng et al., 2011) or SDS-digested freeze fracture replica labeling (Masugi-Tokita et al., 2007; Tarusawa et al., 2009; Budisantoso et al., 2012; Fukazawa and Shigemoto, 2012) have suggested a nonhomogeneous distribution of receptors within the PSD. SDS-digested freeze fracture replica labeling revealed a microcluster intrasynaptic arrangement of AMPARs in corticogeniculate synapses and a more diffuse arrangement in retinogeniculate synapses (Tarusawa et al., 2009). Similarly, AMPARs are organized in intrasynaptic microclusters in parallel fiber-Purkinje cell synapses, whereas they are homogeneously distributed in climbing fiber-Purkinje cell synapses (Masugi-Tokita et al., 2007). AMPAR density in the clusters reported in these preparations ranged from ∼500 to 2000 AMPAR/μm2. We found that AMPAR density inside nanodomains reached ∼2500 AMPAR/μm2, which is comparable. The number of nanodomains detected per synapse was slightly different in between the different techniques, although in the same range, with more domains detected by EM compared with dSTORM, sptPALM, and uPAINT. This can be explained by the more efficient detection efficacy reached by EM.
A novel approach to immunogold labeling in which neurons were live-labeled with a highly avid GluA2 primary antibody, followed by pre-embedding immunogold secondary application and silver intensification, resulted in high-density labeling of surface GluA2 and the presence of AMPAR nanodomains. The underestimation of the degree of AMPAR clustering in earlier EM studies (Kharazia et al., 1996; Petralia et al., 1997; Nusser et al., 1998; Takumi et al., 1999a) could arise from undersampling of the epitopes with postembedding methods and use of different primary antibodies. Importantly, our super-resolution experiments in live cells indicate that the AMPAR organization into nanodomains does not result from a fixation artifact.
The stabilization of AMPARs at postsynaptic sites is regulated by the interaction of AMPAR subunits and auxiliary proteins, such as TARPs, with intracellular and extracellular interactors (Nicoll et al., 2006; Bats et al., 2007; Frischknecht et al., 2009). AMPAR nanodomains are relatively stable, as they can last for at least an hour, and most likely originate from interactions of AMPAR complex components with intracellular scaffolds, such as PSD95. Stabilization of AMPARs in nanodomains is reversible in the time span of seconds, and diffusion outside nanodomains is as rapid as in the extrasynaptic membrane. This observation reaffirms that AMPARs are in a dynamic equilibrium between sites of retention, consistent with the notion that postsynaptic slots control the heterogeneity in receptor immobilization locations within the synapse (Lisman and Raghavachari, 2006).
The size of nanodomains is surprisingly homogeneous, suggesting that it may be set by the physicochemical properties of a limited number of interacting building blocks. This subsynaptic organization is linked to the homeostasis of PSD95, as an increase or decrease in PSD95 expression levels modifies the nanodomain characteristics. A decrease in endogenous PSD95 expression levels broadens the surface of nanodomains and decreases AMPAR content and mEPSC amplitude. This is likely the result of an alteration in the nanodomain structure and AMPAR destabilization. This is in agreement with the recent EM tomography observation of RNA interference knockdown of PSD95 leading to patchy loss of PSD material. This loss correlates with the loss of putative AMPAR-type structures in the PSD (Chen et al., 2011). In contrast, overexpression of PSD95 increases AMPAR nanodomain size, but not its total content (i.e., lowering AMPAR density in the nanodomain). PSD95 overexpression does not significantly change the mEPSC amplitude, in agreement with previous work (Ehrlich and Malinow, 2004). However, it should be noted that PSD95 is most likely not the only core organizer of nanodomains, in particular as AMPAR and PSD95 do not perfectly colocalize, and as knocking down PSD95 (this study) or disrupting the PSD95-Stargazin interaction (Sainlos et al., 2011) does not fully disrupt nanodomain organization or synaptic transmission.
Together, we find it particularly interesting that mEPSC amplitude is well correlated with AMPAR nanodomain total content and not their size. Indeed, the amplitude of unitary synaptic responses is determined largely by the number of receptors binding neurotransmitters after presynaptic terminal release. This parameter is directly linked to the q value, which is the size of a quantum of glutamate inside one vesicle. Previously, this quantum was thought to be sufficient to saturate all synaptic AMPARs. More recently, biophysical models of glutamatergic synaptic transmission have outlined that, because of the steep glutamate gradient in the synaptic cleft after vesicle release, AMPARs can only be activated in an area of ∼150 nm (Lisman et al., 2007; Wu et al., 2007; Budisantoso et al., 2012) around the release site. As a consequence, organization of AMPARs inside the PSD, together with their proximity with presynaptic release sites, tune the strength and the variability of synaptic responses (Franks et al., 2002; Franks et al., 2003). Existing models of synaptic transmission rely mostly on a classical uniform distribution of AMPARs at the PSD. It means that, even if the q value is not sufficient to saturate all the receptors of the synapse, the AMPAR response will be directly correlated with the size of the quantum: the bigger the quantum, the more receptors are opened. The demonstration that AMPARs are highly concentrated in nanodomains and that these domains are long lasting opens the new possibility that single vesicle release may activate a nanodomain, and this could represent a postsynaptic quantum. Synaptic transmission could then be powerfully regulated by the composition and the position of the nanodomains relative to the presynaptic release site.
To investigate this concept further, we measured the distance between nanodomains inside single spines by STORM and sptPALM and performed modeling. In PSDs bearing multiple nanodomains, only approximately half of them are closer one to the other by <250 nm. With respect to the area over which AMPARs are thought to be activated by glutamate (diameter of <200 nm) (Lisman et al., 2007), FWHM computed from our simulations 156 ± 15 nm (Fig. 10B2), this implies that a large fraction of AMPAR nanodomains are functionally isolated one from the other with respect to single vesicle release (Fig. 10J1). Using Monte Carlo simulations, we also found that the distance between glutamate release site and an AMPAR nanodomain strongly impacts the efficacy of transmission. The identity, location, and number of vesicle release sites at individual terminals are still a matter of debate. Nevertheless, it is tempting to postulate a correlation between release sites and AMPAR nanodomains, likely through trans-synaptic adhesion proteins, such as neurexin/neuroligins.
Finally, previous work demonstrated that fast AMPAR diffusion in and out release sites renews the pool of naive receptors that can be activated and removes desensitized ones, allowing the synapse to sustain high-frequency stimulation (Heine et al., 2008). We demonstrated that this effect of AMPAR diffusion is possible only if AMPAR exchange takes place on a small area (diameter of a couple of hundred nanometers) or if AMPARs are moving rapidly. The concentration of AMPARs in nanodomains, and their high mobility outside these, provides a framework for sufficiently fast exchange in AMPARs as well as for changes in receptor density in front of release sites by fast aggregation and dissociation of nanodomains aided by lateral diffusion.
Footnotes
This work was supported by the Centre National de la Recherche Scientifique, the Agence Nationale pour la Recherche (StimTrafPark), the Regional Council of Aquitaine, the Fondation pour la Recherche Médicale, the European Union's seventh framework program for research and development for the Marie-Curie-Intra European Fellowship, the ERC GrantNano-Dyn-Syn, and the ERA-NET project MODIFSYN. We thank C. Breillat, D. Bouchet, B. Tessier, A. Frouin, L. Royer, and N. Retailleau for help in molecular biology and cell cultures; A. Penn for providing custom-made mEPSC analysis software; M. Heilemann for helpful advice on dSTORM; members of the Bordeaux Imaging Center, especially P. Legros for the use of the STED microscope; S. Lacomme and E. Gontier of the Electron Microscopy Facility; S. Marais for help with EM reconstructions; L. Cognet and N. Rebola for helpful discussions; and E. Gouaux for the generous gift of the anti-GluA2 antibody.
The authors declare no competing financial interests.
- Correspondence should be addressed to either Dr. Daniel Choquet or Dr. Jean-Baptiste Sibarita, University of Bordeaux, Interdisciplinary Institute for Neuroscience, and Centre National de la Recherche Scientifique, Unité Mixte de Recherche 5297, 146 rue Léo Saignat, 33077, Bordeaux Cedex, France. dchoquet{at}u-bordeaux2.fr or jean-baptiste.sibarita{at}u-bordeaux2.fr