Abstract
The structural plasticity of dendritic spines is considered to be an important basis of synaptic plasticity, learning, and memory. Here, we induced input-specific structural LTP (sLTP) in single dendritic spines in organotypic hippocampal slices from mice of either sex and performed ultrastructural analyses of the spines using efficient correlative light and electron microscopy. We observed reorganization of the PSD nanostructure, such as perforation and segmentation, at 2–3, 20, and 120 min after sLTP induction. In addition, PSD and nonsynaptic axon–spine interface (nsASI) membrane expanded unevenly during sLTP. Specifically, the PSD area showed a transient increase at 2–3 min after sLTP induction. The PSD growth was to a degree less than spine volume growth at 2–3 min and 20 min after sLTP induction but became similar at 120 min. On the other hand, the nsASI area showed a profound and lasting expansion, to a degree similar to spine volume growth throughout the process. These rapid ultrastructural changes in PSD and surrounding membrane may contribute to rapid electrophysiological plasticity during sLTP.
SIGNIFICANCE STATEMENT To understand the ultrastructural changes during synaptic plasticity, it is desired to efficiently image single dendritic spines that underwent structural plasticity in electron microscopy. We induced structural long-term potentiation (sLTP) in single dendritic spines by two-photon glutamate uncaging. We then identified the same spines at different phases of sLTP and performed ultrastructural analysis by using an efficient correlative light and electron microscopy method. We found that postsynaptic density undergoes dramatic modification in its structural complexity immediately after sLTP induction. Meanwhile, the nonsynaptic axon–spine interface area shows a rapid and sustained increase throughout sLTP. Our results indicate that the uneven modification of synaptic and nonsynaptic postsynaptic membrane might contribute to rapid electrophysiological plasticity during sLTP.
Introduction
Dendritic spines are the major sites for receiving excitatory synaptic inputs and play important roles in neuronal signal transduction, memory storage, and neuronal circuit organization. Plasticity of spine structure plays a critical role in electrophysiological plasticity, learning, and memory (Sala and Segal, 2014). In particular, it has been reported that activity-dependent spine enlargement [i.e., structural LTP (sLTP)] is correlated with LTP (Lang et al., 2004; Matsuzaki et al., 2004; Kopec et al., 2006) and memory formation (Hayashi-Takagi et al., 2015). Ultrastructural correlates of LTP have also been studied using electron micrographs of LTP-induced brain slices under a condition where a large fraction of spines was electrically stimulated. These studies revealed that LTP is correlated with, on average, a larger spine volume and PSD area, a higher number of concave spine profiles, an increase in perforated PSD, and a higher fraction of polyribosome-containing spines (Van Harreveld and Fifkova, 1975; Desmond and Levy, 1986; Harris et al., 1992; Toni et al., 2001; Ostroff et al., 2002). Furthermore, LTP was found to be associated with the generation of multi-innervating spines and multiple spine synapses (Toni et al., 1999; Giese et al., 2015).
However, strong and synchronous LTP stimulation induced in a large population of dendritic spines makes it challenging to differentiate spines that receive direct LTP stimulation from spines that are influenced by surrounding ones. Recent studies have addressed this issue by inducing sLTP and LTP with two-photon (2P) glutamate uncaging in single dendritic spines, and the spines have been visualized under transmission electron microscopy (TEM) by correlating light and electron microscopy (Bosch et al., 2014; Meyer et al., 2014). It has been found that, while the spine volume and synaptic transmission increase rapidly after sLTP induction, the PSD area increases slowly over several hours (Matsuzaki et al., 2004; Bosch et al., 2014; Meyer et al., 2014). The difference in timing implies the reorganization of substructures of synapses during sLTP. However, the throughput of these methods has been low, and thus detailed analyses of PSD morphology have not been done.
In addition to these EM studies, nanometer-scale PSD structure and plasticity have also been studied using optical super-resolution microscopy (MacGillavry et al., 2013; Nair et al., 2013; Hruska et al., 2018; Masch et al., 2018; Wegner et al., 2018). These studies suggest that, in dissociated culture, the synaptic scaffold is made of isolated nanoclusters (MacGillavry et al., 2013; Nair et al., 2013; Hruska et al., 2018), which increase in number during chemically induced LTP (MacGillavry et al., 2013; Hruska et al., 2018). However, a recent in vivo study using stimulated emission depletion microscopy suggested that PSD-95 does not form segregated clusters, but nearly half of the postsynaptic scaffolds show complex spatial organization (Masch et al., 2018). Thus, overall, how the PSD and surrounding membrane change their nanometer-level structure during sLTP still remains elusive.
Here, we demonstrate that single spines stimulated and imaged under a 2P microscope can be reliably and efficiently reimaged with automated tape-collecting ultramicrotome (ATUMtome) sectioning and scanning EM (SEM). This technique allows for clear visualization of subspine ultrastructures, including the PSD and synaptic membrane (Kamasawa et al., 2015; Sun et al., 2020). 3D reconstruction of the stimulated spines showed that the PSD underwent rapid and sustained structural changes, including perforation and segmentation, while its size slowly developed. Associated with changes in the PSD structure, we also observed an immediate and sustained expansion of the nonsynaptic axon–spine interface (nsASI), which has a narrow extracellular space bordered by axon and spine membrane without PSD. The nsASI expansion, coordinated with the increase in structural complexity of PSD, may support a rapid increase in synaptic transmission during sLTP.
Materials and Methods
Hippocampal organotypic slice culture and gene gun transfection.
All of the animals were handled according to approved Institutional Animal Care and Use Committee protocols of the Max Planck Florida Institute for Neuroscience. Hippocampal organotypic slice cultures were prepared from postnatal day 4–6 mice of either sex, as previously described (Gogolla et al., 2006). After 10–12 d in culture, CAG-EGFP plasmid (catalog #16664, Addgene) transfection was performed by a Helios gene gun (Woods and Zito, 2008). Two days after transfection, slices with sparsely labeled CA1 pyramidal neurons were selected for glutamate uncaging.
2p glutamate uncaging.
4-Methoxy-7-nitroindolinyl-caged-L-glutamate (MNI)-glutamate uncaging was performed with 2P laser microscopy (720 nm) with 4 mW power, 6 ms dwell time, 0.5 Hz × 30 stimulation for LTP induction at room temperature. The buffer for 2P imaging and uncaging was Mg2+-free ACSF containing 127 mm NaCl, 2.5 mm KCl, 25 mm NaHCO3, 1.25 mm NaH2PO4, 4 mm CaCl2, and 25 mm d-glucose plus 1 μm TTX and 4 mm MNI-glutamate aerated with 95% O2 and 5% CO2 (Hedrick et al., 2016). The slices were kept in ACSF for at least 10 min before uncaging and during imaging until fixation. Target spines were selected on secondary dendrites of a CA1 pyramidal neuron (Fig. 1A, Extended Data Fig. 1-1A1,A2). Uncaging was performed at ∼0.5 µm from the tip of the selected spine (Fig. 1B, Extended Data Fig. 1-1A3). Zoom 1 (z-stack, 5 µm/step, 200 × 200 µm2), zoom 5 (z-stack, 3 µm/step, 40 × 40 µm2), and zoom 25 (z-stack, 1 µm/step, 8 × 8 µm2) images of the whole neuron and specific uncaged region were acquired within 5 min before uncaging with 2P microscopy, and zoom 25 images of the same region were taken immediately after uncaging to confirm the success of stimulation (Fig. 1C, Extended Data Fig. 1-1A4). Hippocampal slices were transferred into a fixative buffer [2% paraformaldehyde and 2% glutaraldehyde in 0.1 m phosphate buffer (PB) at pH 7.4] at different time points (2–3 min, 20 min, 2 h) after uncaging and kept for 1 h. Slices were then washed in 0.1 m PB three times for a total of 30 min. Zoom 25 images of target spines were taken again to confirm the maintenance of sLTP (Extended Data Fig. 1-1B1).
Laser-burning marks.
Some slices were chosen to test additional fiducial marks for correlation. For these slices, burning marks were introduced with a 2P laser near the target neurons and target secondary dendrites after slice fixation, positioned > 5 µm from the target spines. For easy operation, an uncaging laser at 720 nm (30 mW) was used in line scan mode to introduce horizontal burning lines while slices kept perfused in 0.1 m PB, pH 7.4. The number of line scans needs to be tested for each sample. Successful burning marks show autofluorescence under 2P microscopy (Extended Data Fig. 1-1B2) and can be visualized in the slice embedded in resin under both light microscopy and EM.
Laser burning significantly improved the efficiency of relocating stimulated spines, especially in estimating the z-position for target neurons and spines. However, laser burning is not a must for our workflow, and should be avoided if a long-distance tracing of axons is desired since the structures of interest might be disrupted.
Confocal scanning microscopy imaging.
Confocal images were taken with a Zeiss LSM780 confocal scanning microscope. GFP+ neurons were viewed with a laser excitation wavelength of 488 nm. Images were taken in tile mode to cover the whole slice and in z-stack mode (5 µm interval) to cover the z-axis range of the target neuron (Extended Data Fig. 1-1B3).
Immuno-EM processing and ultra-thin sectioning.
Slices were incubated in 50 mm glycine in 0.1 m PB for 10 min to block excessive aldehyde sites. Slices were then cryoprotected with 15−30% sucrose, followed by two freeze–thaw cycles of 1 min liquid nitrogen permeation. After blocking with normal goat serum and fish skin gelatin, the slices were incubated with anti-GFP primary antibody (0.1 µg/ml; catalog #ab6556, Abcam) for 40 h, followed by nanogold-conjugated secondary antibody (1:100; catalog #2003, Nanoprobes) for 16 h. Silver enhancement was then performed with the HQ Silver Enhancement Kit (catalog #2012, Nanoprobes) to increase the visibility of gold particles. Slices were then postfixed with 0.5% aqueous OsO4 for 40 min at 4°C; stained with 1% aqueous uranyl acetate for 35 min; dehydrated with gradually increased concentrations of ethanol, acetone, and propylene oxidase; and infiltrated by Durcupan ACE (Sigma-Aldrich).
Slices were trimmed down to the region containing the target neuron according to confocal image localization and light microscopy visualization. The final sample size was ∼1 × 1.5 mm. The samples were then sectioned with a 4 mm Diatome 35° knife by ATUMtome (RMC/Boeckeler Instruments) at a thickness of 50 nm and were collected onto rolling Kapton tape (Extended Data Fig. 1-1C1). The Kapton tape was then cut and aligned on 4 inch silicon wafers covered with carbon tape. Three 150 mesh copper grids were put onto blank regions of carbon tape as fiducial marks (Extended Data Fig. 1-1C2). A 5-nm-thick layer of carbon was then coated onto the surface of the wafer using a high-vacuum sputter coater (model ACE600, Leica).
SEM imaging.
SEM imaging was performed with Zeiss Merlin VP Compact using Zeiss Atlas 5 AT software (version 5.0.49.4). An overview image of the wafer was taken with a digital camera for the alignment of the physical position of the wafer in the SEM chamber (Extended Data Fig. 1-1C2). The alignment was performed by using the grids on the wafer as fiducial marks. Approximate locations of sections could be visualized in the digital camera overview image. Then low-resolution (2 µm/pixel) images of whole sections were first taken to accurately locate the position of the sections by secondary electron (SE2) detector (2.5–3 kV, 0.8 µs dwell time; Extended Data Fig. 1-1D1). The relative location of the neuron labeled with immunogold in the serial sections was then estimated according to the overview confocal image. A backscattered electron (BSE) detector (BSD4A, 8–10 kV, 10 µs dwell time) was used at this step to search for profiles labeled with the silver-enhanced nanogold particles. This was followed by the acquisition of a series of medium-resolution (30 nm/pixel) images of the target neuron (Extended Data Fig. 1-1D2). A rough 3D reconstruction of the neuron containing its soma and target secondary dendrite regions was made by Amira software (version 2019.3; Thermo Fisher Scientific; Extended Data Fig. 1-1E1). By comparing the shape of the 3D reconstruction of a labeled neuron and the originally captured 2P microscopy image, the target neuron was confirmed, and the target secondary dendrite was determined. For slices with laser-burning marks, the step to reconstruct the whole neuron for the correlation was omitted. Burning marks were easily found in EM as empty straight lines (Extended Data Fig. 1-1D3) and could be used directly for the identification of targeted dendritic segments.
Final high-resolution images (4 nm/pixel) were then taken using an In-lens Duo Detector in BSE mode (2.5 kV, 75 µs dwell time) or an aBSD detector (8–10 kV, 25–40 µs dwell time) to visualize ultrastructure of the target spines and their surrounding regions (Fig. 1D, Extended Data Fig. 1-1D4). Target spines were identified by comparing the spine shape and location between 3D reconstructed dendritic segments and 2P images. The ultrastructure required for analyzing spines and synapses, such as presynaptic and postsynaptic membranes, synaptic cleft, and PSD, is clearly visualized. 3D reconstruction of dendritic spines was performed in Amira (Fig. 1E, Extended Data Fig. 1-1E2).
In this study, the success rate for collecting final high-resolution images for the target spines was 44% (15 of 34 slices). The reasons for failed experiments include the following: (1) poor ultrastructure preservation; (2) weak immunoreaction to locate target neuron under a dissecting microscope; (3) heavy silver enhancement causing large deposition of nanogold particles disturbed the ultrastructure observation; (4) nonoptimal conditions of the Kapton tape (having scratches or debris on its surface) and/or ultra-thin sections (having wrinkles, chatters, or damages) that interrupted target spine imaging; and (5) accident loss of the region of interest in the trimming process. The success rate significantly increased with experience and improved quality control on slice health and Kapton tape.
Estimation of spine volume before 2P glutamate uncaging.
sLTP spines were imaged in 2P microscopy within 5 min before glutamate uncaging and after chemical fixation. The spine volume increase was quantified as a change in GFP fluorescence intensity. The spine volume after fixation was calculated from 3D reconstruction of the spines with EM images, and the spine volume before uncaging was calculated as the spine volume measured in EM divided by the volume change measured optically.
Image analysis.
ASI was traced as the postsynaptic part of the interface between the spine head and the opposing presynaptic bouton. PSD was identified as the electron-dense region in ASI, and nsASI was identified as the ASI membrane without PSD.
Spine shape complexity and PSD complexity.
Dendritic spines and PSDs were reconstructed in 3D in Amira, then exported into the MATLAB format to create a shadowed PSD region on the spine surface. Spine head volume and surface area were measured after reconstruction and were used to calculate the spine shape complexity index with the following formula, with which a sphere-shaped spine head generates an index equal to 1, as follows:
The smoothed PSD area and perimeter were measured and used to calculate the complexity index with the following formula, with which a circle-shaped PSD generates an index equal to 1, as follows:
Experimental design and statistical analyses.
We introduced single-spine sLTP using 2P glutamate uncaging in hippocampal organotypic slice culture because of their easy operation and stability to produce good ultrastructure preservation. We imaged ∼100 sections around the sLTP spines. All spines located on the same secondary dendrite with sLTP spines were used as control spines, except those with PSD parallel to the section plane. Based on the dendrite orientation, 1–31 control spines were taken for each sLTP spine, and control spines were all located within 20 µm from the sLTP spine. Since we compare reconstructed dendrites with 2P images to ensure correct correlation, it is not feasible to make the spine ultrastructure measurement process blind.
Statistical tests were performed in GraphPad Prism 6. Measured ultrastructural data were first tested for normal distribution and log-normal distribution with a D'Agostino–Pearson test. Data without skewness were compared with Welch's t test and shown with the mean and SEM for each control versus sLTP group. Data without skewness after log-transformation (consistent with log-normal distribution) were log transformed and compared with the Welch's test. For data that did not follow normal distribution or log-normal distribution, a Mann–Whitney test was performed for comparison. To compare the growth of PSD, ASI, and nsASI with spine volume growth, their ratio was analyzed with two-way ANOVA. Data consistent with log-normal distribution were log transformed before two-way ANOVA. A post hoc Tukey's multiple-comparison test was performed to compare among different sLTP phases, and a post hoc Bonferroni's multiple-comparison test was performed to compare between control and sLTP conditions. An extra sum-of-squares F test was used to compare regression lines. All statistics are summarized in Tables 1 and 2. Raw data are available in Extended Data Figure 1-2.
Comparison of the ultrastructures of dendritic spines in control and sLTP conditions at early, intermediate and late phases of sLTP
Comparison of the relative changes of dendritic spine ultrastructures in control and sLTP conditions at early, intermediate, and late phases of sLTP
Results
Correlating dendritic spines imaged by 2P microscopy and SEM
We established an efficient method to analyze the 3D ultrastructure of spines undergoing spine-specific sLTP in organotypic hippocampal slices of mice. Spine structural plasticity was first imaged in 2P microscopy and then identified in EM by combining pre-embedding immunolabeling of GFP with ATUMtome serial sectioning and SEM imaging (Fig. 1A–E, Extended Data Fig. 1-1, whole flow). We induced sLTP at single spines with glutamate uncaging by delivering 2P laser pulses located ∼0.5 µm from the tip of the spine head in GFP-expressing CA1 pyramidal neurons in organotypic hippocampal slices (Fig. 1A–C, Extended Data Fig. 1-1A; Matsuzaki et al., 2001, 2004; Harvey et al., 2008). These slices were then fixed at different time points after uncaging stimulation to capture the structural changes at the following different sLTP phases: early phase (2–3 min), intermediate phase (20 min), and late phase (120 min; Extended Data Fig. 1-1B1). GFP expression in entire slices was imaged in confocal microscopy, and optional burning marks were introduced near the target dendritic segment with 2P lasers to facilitate EM correlation (Extended Data Fig. 1-1B2,B3). After immunostaining of GFP with nanogold particles, we performed silver enhancement, heavy metal staining, dehydration, and resin embedding. Using ATUMtome, we cut and collected ultra-thin sections (50 nm), which were imaged in SEM at low resolution (2 µm/pixel) and intermediate resolution (30 nm/pixel; Extended Data Fig. 1-1C,D1,D2). Burning marks facilitated the correlation as they can be easily visualized in EM as empty space (Extended Data Fig. 1-1D3). A rough 3D reconstruction of the nanogold-positive profile was used for comparison with 2P neuronal images to identify target neurons (Extended Data Fig. 1-1E1). Following the identification of target neurons and secondary dendrites, we imaged spines that underwent sLTP (sLTP spine) and multiple spines located within the imaging range (∼100 sections) on the same dendrite (control spine) with high resolution (4 nm/pixel; Fig. 1D, Extended Data Fig. 1-1D4). 3D morphology of dendritic spines and PSDs were reconstructed manually from the sections and subjected to further analyses (Fig. 1E, Extended Data Fig. 1-1E2).
2P glutamate uncaging at single spine induces spine expansion. A, 2P image of an EGFP+ CA1 pyramidal neuron in mice organotypic hippocampal slice selected for glutamate uncaging. Red arrow points to the location of target spine. Scale bar, 20 µm. B, C, 2P images of a dendritic spine before (B) and after (C) 2P glutamate uncaging stimulation. Red dot in B shows the uncaging position, and red arrow in C points to the enlarged dendritic spine. Scale bar, 1 µm. D, A 4 nm/pixel high-resolution SEM image taken with an InlensDuo detector shows detailed ultrastructure of the modulated dendritic spine in C. Spine head is shown in yellow. Scale bar, 0.5 µm. E, 3D reconstruction of the same dendritic spine shown in D. Spine head is shown in yellow, and PSD is shown in red. Scale bar, 0.5 µm. The full workflow for correlative light and electron microscopy is shown in Extended Data Figure 1-1. F, The volume of control spines was similar to the volume of sLTP spines before stimulation. Welch's t test was used on log-transformed data. G, Histograms with log-normal distribution fit curves of raw spine head volume for control spines in early, intermediate, and late phases of sLTP. H, Spine head volume for control spines and random spines located >10 µm from sLTP spine at late phase of sLTP. Welch's t test was used on log-transformed data. I, Cumulative distribution of raw spine head volume for control and sLTP spines in early, intermediate, and late phases of sLTP. J, Spine head volume for control and sLTP spines in early, intermediate, and late phases of sLTP. Welch's t test was used on log-transformed data. K, Spine shape complexity index for control (Ctrl) and sLTP spine at early, intermediate, and late phases of sLTP. Mann–Whitney t test was used. Raw measurements for spine parameters and 3D reconstruction for each spine is shown in Extended Data Figure 1-2. *p < 0.05, **p < 0.01, ****p < 0.0001.
Figure 1-1
Correlative light and electron microscopy workflow of a 2P uncaging modulated dendritic spine. Left, Workflow and estimated time cost. A1, A2, Zoom 1 (200 × 200 µm2) and zoom 5 (40 × 40 µm2) 2P images of an EGFP+ CA1 pyramidal neuron in organotypic hippocampal slice selected for glutamate uncaging. Target dendrite pointed out by arrow. Scale bars: A1, 20 µm; A2, 5 µm. A3, A4, Zoom 25 (8 × 8 µm2) 2P images of a dendritic spine before (A3) and after (A4) 2P glutamate uncaging stimulation. Red dot in A3 shows the uncaging position, and red arrow in A4 points to the enlarged dendritic spine. Scale bar, 1 µm. B1, Zoom 25 2P image of the uncaged spine after fixation. Arrow points at the target spine. Scale bar, 1 µm. B2, Laser-burning fiducial mark is introduced by 2P laser near the target spine after tissue fixation. Scale bar, 1 µm. B3, A confocal image shows the overview of GFP expression pattern in the whole hippocampal organotypic slice. Target neuron is labeled in red square. Scale bar, 200 µm. C1, ATUMtome setup for sectioning and collecting ultra-thin sections. C2, Silicon wafer holding aligned Kapton tape containing serial sections. Three copper grids on the edge of the wafer are used as fiducial marks for alignment. D1, A 2 µm/pixel image taken with SE2 detector shows the outline of two ultra-thin sections. Scale bar, 200 µm. D2, A 30 nm/pixel image taken with BSE detector shows tissue structure and the silver enhanced nanogold labeling pattern. Nanogolds are seen as black dots in the image. Scale bar, 50 µm. D3, A 4 nm/pixel high-resolution image taken with an InlensDuo detector shows the laser-burning mark (white empty region) and nearby nanogold-labeled structure (boxed in red and zoomed to show details). Scale bar, 2 µm. D4, A 4 nm/pixel high-resolution image taken with an InlensDuo detector shows the ultrastructure of the modulated dendritic spine in A4. Scale bar, 0.5 µm. E1, Partial 3D reconstruction of the GFP+ neuron shown in A1. Targeted dendrite pointed out by arrow. Scale bar, 20 µm. E2, 3D reconstruction of the target dendritic spine. Spine head is shown in yellow, and PSD is shown in red. Scale bar, 200 nm. Download Figure 1-1, TIF file.
Figure 1-2
Raw measurements for spine parameters and 3D reconstruction for each spine. Download Figure 1-2, XLSX file.
The mean ± SEM volume of control spines was 0.060 ± 0.003 µm3 (n = 194), which is consistent with previous studies (Harris et al., 1992). From fluorescence changes following sLTP and the volume of sLTP spines in EM, we estimated the average volume of sLTP spines before the stimulation. The estimated volume (mean ± SEM, 0.063 ± 0.006 µm3; n = 36) was similar to the volume of control spines, suggesting that our sLTP spines represent average spines in terms of their volume before sLTP (Fig. 1F, Extended Data Fig. 1-2). We found that control spine volume distribution was skewed from the normal distribution and can be expressed by log-normal function (Fig. 1G; Loewenstein et al., 2011). Thus, for the rest of the study, we performed parametric statistics after log transformation for datasets consistent with a log-normal distribution. Spine volume in control spines in the late phase was significantly smaller than that in the early phase (ANOVA on log-transformed data, p =0.002). However, since the same trend was also observed in spines on different neurons located away from the stimulated spines (>10 µm; Fig. 1H), this phenomenon is likely because of the prolonged incubation of slices in ACSF containing MNI-glutamate and TTX, instead of the heterosynaptic LTD (Harvey and Svoboda, 2007; Oh et al., 2015; Hedrick et al., 2016). For such reasons, we compared only sLTP spines and surrounding spines fixed at the same time point. As expected, sLTP spines showed a significantly larger size during the early and intermediate phases compared with control spines. At the late phase, the volume of sLTP spines was larger than that before stimulation (mean ± SEM fold change, 1.33 ± 0.06; Extended Data Fig. 1-2). However, because of the large variability in control spine volume and the relatively small increases in sLTP spine volume at the late phase, the spine volume was not statistically different between the control and sLTP groups (Fig. 1I,J, Table 1). Similarly, we found that the spine surface area increased after sLTP induction during early and intermediate phases, but not in the late phase (Table 1). Furthermore, the shape of the spine head became dramatically more complicated after sLTP induction. To quantify the spine shape complexity, we defined the spine shape complexity index as the ratio of the square root of the spine surface to the cube root of the spine volume, normalized so that a sphere spine has an index equal to 1. The index was significantly higher in sLTP spines at all phases (Fig. 1K).
PSD reorganization after sLTP induction
We next reconstructed the PSD in 3D for all of the spines serially imaged to further characterize their structural features (Fig. 2A,B). The PSD area size in sLTP spines was significantly larger than that in control spines at the early phase of sLTP (Fig. 2C, Table 1), albeit to a lower degree than the spine volume change (Fig. 2D,E). In control spines, the PSD area correlated tightly with the spine volume (Fig. 2D), which is consistent with previous studies (Harris and Stevens, 1989; Bosch et al., 2014). sLTP induction broke the correlation between PSD size and spine volume. The ratio of the PSD size to spine head volume dramatically decreased at the early and intermediate phases, suggesting that the PSD growth was much smaller than the spine volume changes during these phases. The slope of regression lines of PSD size on spine volume at the late phase for sLTP spines was higher than that for control spines (Fig. 2D3). To compare PSD area and spine volume growth, we analyzed the ratio of PSD area to spine volume at different time points (Fig. 2E). Two-way ANOVA for stimulation conditions (control and sLTP) and time points (early, intermediate, and late) showed that the datasets were significantly different between control and sLTP spines (p < 0.0001), and also over time points (p = 0.0043), with significant interaction between stimulation conditions and time points (p = 0.0059). Post hoc multiple-comparison analyses showed that sLTP spines had significantly smaller PSD-to-spine volume ratios in the early phase (p < 0.0001) and intermediate phase (p = 0.0036), but not in the late phase (p = 1.0). Furthermore, sLTP spines had a significant increase in the ratio in the late phase, compared with the early phase (p = 0.0108) and intermediate phase (p = 0.0096; Table 2). Thus, these results suggest that PSD growth is slower than spine volume growth by approximately hours, which is consistent with previous studies (Bosch et al., 2014; Meyer et al., 2014).
sLTP induction disrupts the correlation between PSD and spine head volume. A, B, Serial images of an example spine (A) and 3D reconstruction of the example spine (B), with spine head shown in yellow and PSD shown in red. Scale bar, 200 nm. C, PSD size at early, intermediate, and late phases of sLTP. Welch's t test on log-transformed data. D, Correlation between spine head volume and PSD size for control spines (Ctrl) and sLTP spines at early, intermediate, and late phases of sLTP. Regression line forced through origin. Black line and open circles represent control spines, and red line and solid dots represent sLTP spines at each phase of sLTP. The p value was calculated by comparing the two regression lines with an extra sum-of-squares F test. E, Ratio of PSD area to spine volume at early, intermediate, and late phases of sLTP. Two-way ANOVA followed by Bonferroni's post hoc test was used for comparison between control and sLTP spines, and Tukey's post hoc test was used for comparison among different sLTP phases. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
We further analyzed the detailed nanometer-scale structure of the PSD (Fig. 3A–H). We observed PSDs with simple, fenestrated, horseshoe, and segmented shapes, as previously reported (Fig. 3A–D). In addition, we also observed PSDs with irregular shapes that did not fit previous categories (Fig. 3E,F) and with shapes including multiple features (Fig. 3G, horseshoe and fenestration, H, segmentation and fenestration). At all of the time points, a higher percentage of sLTP spines contained PSDs with complex structural features, including fenestration, horseshoe shape, and segmentation, compared with control spines (Fig. 3I). In the early phase, nearly half of the stimulated spines contained segmented PSD, composed of either two or three fragments (Fig. 3J). To quantitatively evaluate the complexity of PSD shape, we calculated a complexity index, defined as the ratio of the square of the perimeter to the area, and normalized so that the index for a circular disk shape equals 1. We found that a higher PSD complexity index was associated with larger spine volume and PSD area (Fig. 3K,L). Because of this correlation, the complexity index was highly variable among different spines. We observed that the PSD shape was significantly more complex in sLTP spines than in control spines at the early phase and showed a trend toward higher complexity at the intermediate and late phases of sLTP (Fig. 3M). Overall, our results indicate that a dramatic PSD shape transformation happened immediately after sLTP induction. We also observed that two sLTP spines with segmented PSDs were multiply innervated (intermediate phase, 1 of 13; late phase, 1 of 13), while three control spines were multiply innervated (early phase, 1 of 80; late phase, 2 of 73; Fig. 4).
sLTP induces increased PSD complexity. A–H, 3D reconstruction of example spines (gray) containing PSD (red) with different complexity indexes (CIs) and different shape categories: simple (A), fenestration (B), horseshoe (C), segmentation (D), irregular horseshoe (E, F), and combined shape categories (G, H). Scale bar, 200 nm. PSD CI is shown for individual spine. I, Percentage of spines contains PSD in simple, fenestrated, horseshoe, and segmented categories at early, intermediate, and late phases of sLTP in sLTP and control (Ctrl) spines. J, Percentage of spines containing one to four PSD segments at early, intermediate, and late phases of sLTP. K, L, PSD shape complexity is positively correlated with spine head volume (K) and PSD size (L) in control spines at all phases of sLTP. r, Pearson's correlation coefficient. M, PSD complexity index at early, intermediate, and late phases of sLTP. Mann–Whitney t test was used. **p < 0.01.
Observation of multi-innervated spines. A, Percentage of single- and multi-innervated spines at early, intermediate, and late phases of sLTP for sLTP and control (Ctrl) spines. B, C, 3D reconstruction of example single-innervated spine (B) and multi-innervated spine (C) with opposing presynaptic boutons and synaptic vesicles. Spines are shown in gray, PSD is shown in red, presynaptic bouton is shown in transparent blue, and synaptic vesicles are shown as golden balls. Scale bar, 500 nm.
Nonsynaptic axon–spine interface expansion after sLTP induction
ASI is the area with thin extracellular space bordered by axon and spine membrane (Fig. 5A,B). A rapid and maintained expansion of ASI was observed after sLTP induction (Fig. 5C). To quantify the time course of ASI growth compared with spine volume increase, we calculated the ratio of ASI area to spine volume using two-way ANOVA and a post hoc multiple comparison. We found that the ratio was significantly smaller for sLTP in the early phase (p = 0.003), but not in the intermediate phase (p = 0.50) and late phase (p = 0.20). Furthermore, in sLTP spines, the ratio is significantly larger in the late phase compared with early phase (p = 0.0002) and intermediate phase (p = 0.031; Fig. 5D, Table 2). Thus, the results suggest that ASI growth is slower than spine volume growth, and the ratio changes progressively toward the control level.
sLTP induces expansion of total ASI postsynaptic membrane. A, B, Serial images of an example spine (A) and 3D reconstruction of the example spine (B) with spine head shown in yellow and ASI shown in green. Scale bar, 200 nm. C, ASI area size at early, intermediate, and late phases of sLTP for sLTP and control (Ctrl) spines. Welch's t test was used on log-transformed data. D, Ratio of ASI area to spine volume at early, intermediate, and late phases of sLTP. Two-way ANOVA followed by Bonferroni's post hoc test was used for comparison between control and sLTP spines, and Tukey's post hoc test was used for comparison among different phases. *p < 0.05, **p < 0.01, ***p < 0.001.
ASI consists of synaptic area, which contains PSD, and nonsynaptic area (nsASI), which lacks PSD. We found that nsASI expanded immediately after sLTP induction and maintained at a larger size until the late phase (Fig. 6A–C, Table 1). Since nsASI may also be able to support synaptic transmission (Ventura and Harris, 1999), this expanded nsASI may contribute to synaptic strengthening during LTP.
sLTP induces sustained expansion of nsASI. A, B, Serial images of an example spine (A) and 3D reconstruction of the example spine (B) with spine head shown in yellow and nsASI area shown in cyan, pointed out by red arrows. Scale bar, 200 nm. C, nsASI size at early, intermediate, and late phases of sLTP. Welch's t test on log-transformed data. D, Ratio of nsASI area to spine volume at early, intermediate, and late phases of sLTP for sLTP and control (Ctrl) spines. Two-way ANOVA followed by Bonferroni's post hoc test was used for comparison between control and sLTP spines, and Tukey's post hoc test was used for comparison among different phases. E, nsASI size is positively correlated with PSD size in control spines at all phases of sLTP. r, Pearson's correlation coefficient. F, Ratio of nsASI area to PSD area at early, intermediate, and late phases of sLTP. Welch's t test on log-transformed data. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
Similar to analyses of PSD and ASI area, we compared the growth of nsASI and spine volume by calculating the ratio of nsASI area to spine volume at different time points (Fig. 6D). Two-away ANOVA for stimulation conditions (control and sLTP) and time points (early, intermediate, and late) showed that the datasets were significantly different over time points (p = 0.0048). Unlike PSD, post hoc multiple-comparison analyses showed that sLTP spines had trended to a larger ratio in the late phase (p = 0.0506), but not in the early phase (p = 0.8004) and intermediate phase (p = 1.0). The sLTP spines showed a significant increase in the ratio in the late phase, compared with the early phase (p = 0.011; Table 2). Thus, these results suggest that nsASI area increases immediately, similar to spine volume growth, and trends to exceed the spine volume growth at the late phase. Furthermore, the size of the nsASI area was proportional to that of the PSD area in control spines (Fig. 6E), which is consistent with a previous report (Ventura and Harris, 1999). However, the ratio of nsASI area to PSD area was higher in sLTP spines than in control spines in all sLTP phases, further proving an uneven expansion of synaptic and nonsynaptic spine membrane after sLTP induction (Fig. 6F).
Discussion
In this study, we developed an efficient method for analyzing ultrastructural changes of single dendritic spines after input-specific sLTP. Using this method, we found that PSD complexity increased within 2–3 min by nanometer-scale structural changes, including perforation and segmentation. In particular, a large fraction of PSDs showed segmentation into two to three segments during the early, intermediate, and late phases of sLTP. Associated with these changes, PSD size increased during the early phase of sLTP, although it was not as pronounced as changes in spine volume changes. On the other hand, our results suggested that the nsASI membrane underwent rapid, pronounced, and sustained expansion, which was proportional to the spine volume growth in the early and intermediate phases but exceeded the spine volume growth at the late phase. Our data indicated that the total increase in ASI area (= nsASI + PSD areas) during sLTP was majorly contributed by the nsASI area expansion.
Previously ultrastructural analyses following sLTP induced by uncaging have been performed using serial sectioned TEM (Bosch et al., 2014; Meyer et al., 2014). However, perhaps because of low yield, these studies restricted their analyses on PSD area size and spine volume. Our methods using ATUMtome serial sectioning and SEM imaging allowed us to achieve correlation more efficiently and to perform detailed analyses of the nanostructure of PSDs and spines. Although SEM imaging generally has a slight disadvantage in terms of resolution, our usage of the advanced detector system visualized sufficient ultrastructural details for the current study.
Organotypic slices have been used in numerous studies to study structural synaptic plasticity and functional synaptic plasticity (Maletic-Savatic et al., 1999; Matsuzaki et al., 2004; Kopec et al., 2007; De Roo et al., 2008; Lee et al., 2009; Makino and Malinow, 2009; Bosch et al., 2014; Meyer et al., 2014). Though the ultrastructure might be different in organotypic slices from those in acute slices and living animals, we chose to use organotypic slices because of their easy operation and stability to produce good ultrastructure preservation. Spine volume in control spines in the late phase was significantly smaller than that in the early phase (ANOVA on log-transformed data, p = 0.002). This is likely because of the incubation of slices in ACSF containing MNI-glutamate and TTX for a long time. This trend was also observed in spines on different neurons located away from the stimulated spines (>10 μm; Fig. 1H). Previous studies also observed changes in ultrastructure after long-term incubation of slices with ACSF, including reduction in PSD area (Bell et al., 2014; Meyer et al., 2014) and active zone size (Bell et al., 2014).
It is known that LTP induction, either by glutamate uncaging or electrical stimulation, is associated with a rapid expansion of spine head volume coupled with an increase in postsynaptic sensitivity (Lang et al., 2004; Matsuzaki et al., 2004; Harvey et al., 2008; Yang et al., 2008b; Lee et al., 2009; Bosch et al., 2014) as well as presynaptic bouton volume change (Meyer et al., 2014) and presynaptic structure reorganization (Bell et al., 2014). The change of spine head volume is similar to the onset and time course of electrophysiological LTP indicated by AMPAR-mediated current (Matsuzaki et al., 2004; Kopec et al., 2006). However, it has been reported that the potentiation of PSD size is much slower than that of volume change, apparently inconsistent with rapid potentiation of AMPAR-mediated current (Bosch et al., 2014).
Our finding of the rapid complication of PSD nanostructure in the early phase, as well as in the intermediate and late phases, associated with a rapid expansion of nsASI area in these phases may be a mechanism underlying rapid recruitment of AMPARs to the synapse. AMPARs exist in extrasynaptic membrane near the PSD (Tao-Cheng et al., 2011; Nair et al., 2013) and likely were able to support synaptic transmission (Franks et al., 2003; Tardin et al., 2003), although their contribution to synaptic transmission is expected to be smaller (Franks et al., 2003). In particular, extrasynaptic AMPARs may play even more significant roles in LTP since LTP is known to be associated with a greater amount of glutamate release (Dolphin et al., 1982; Bekkers and Stevens, 1990; Malinow and Tsien, 1990; Zakharenko et al., 2001). This idea is consistent with previous electrophysiological experiments showing rapid recruitment of AMPARs to perisynaptic sites before the full expression of LTP (Yang et al., 2008a).
We induced sLTP at single spines with glutamate uncaging by delivering 2P laser pulses located ∼0.5 µm from the tip of the spine head in GFP-expressing CA1 pyramidal neurons in organotypic hippocampal slices (Fig. 1A–C, Extended Data Fig. 1-1A; Matsuzaki et al., 2001, 2004; Harvey et al., 2008). This protocol is known to increase AMPAR-evoked current on the spine (Matsuzaki et al., 2004; Harvey et al., 2008). While this protocol is likely to activate extrasynaptic receptors (Franks et al., 2003), the plasticity depends on NMDAR and CaMKII, but not on mGluR, in a manner similar to that of electrophysiological LTP (Matsuzaki et al., 2004; Lee et al., 2009; Chang et al., 2017; Murakoshi et al., 2017). However, since previous studies suggest that metabotropic receptors (mGluRs) exist in perisynaptic areas within ∼60 nm of synapses (Luján et al., 1996), the expansion of nsASI area observed in all phases of LTP in our study may cause enhanced mGluR function after LTP induction.
We observed that two sLTP spines containing segmented PSDs are innervated by multiple presynaptic boutons. An increase of multi-innervated spines has been reported to be associated with LTP (Nikonenko et al., 2003; Radwanska et al., 2011; Giese et al., 2015). However, probably because of the rare existence of multi-innervated spines (∼1% of total spines), it has rarely been reported in previous population EM studies that analyzed the ultrastructural correlates for LTP. Since LTP is known to be correlated with a larger number of synapses in CA1 neurons (Tominaga-Yoshino et al., 2008; Watson et al., 2016), the rapid multi-innervated spines formation may serve as the basis for new synapse formation during sLTP. The multi-innervated spines observed in our experiments that synapsed with different presynaptic boutons likely support the theory of new bouton recruitment instead of spine splitting during LTP (Fiala et al., 2002).
One potential caveat of this study is that the uncaging laser irradiation may have caused some changes in the ultrastructure of dendritic spines. While previous studies have suggested that there is minimum impact on the electrophysiological properties from an uncaging laser (Matsuzaki et al., 2001, 2004; Lee et al., 2009; Murakoshi et al., 2011), this point remains to be resolved in future work.
Overall, our ultrastructural study revealed a time course of coordinated synaptic and extrasynaptic ultrastructural changes during sLTP induced at a single spine. Our results suggest that the combination of rapid and sustained increases in PSD complexity and nsASI membrane, as well as slow development of PSD size, are all likely substrates of synaptic potentiation at different temporal phases.
Footnotes
This work was supported by National Institutes of Health Grants R35-NS-116804, DP1-NS-096787, and R01-MH-080047, and the Brain Research Foundation. We thank Jessica Martin for performing some of the initial 3D reconstruction, Jaime Richards for preparing organotypic slice cultures, Takayasu Minuki and Debbie Guerrero-Given for technical assistance, and Lesley Colgan for critical reading. We also thank the Max Planck Florida Institute for Neuroscience Imaging Center for instrument and technical support.
The authors declare no competing financial interests.
- Correspondence should be addressed to Ryohei Yasuda at ryohei.yasuda{at}mpfi.org