Dendritic spines are post-synaptic compartments of most excitatory synapses. Alterations to dendritic spine morphology occur during neurodevelopment, aging, learning, and many neurological and psychiatric disorders, underscoring the importance of reliable dendritic spine analysis. This protocol describes quantifying dendritic spine morphology accurately and reproducibly using automatic three-dimensional neuron reconstruction software.
Synaptic connections allow for the exchange and processing of information between neurons. The post-synaptic site of excitatory synapses is often formed on dendritic spines. Dendritic spines are structures of great interest in research centered around synaptic plasticity, neurodevelopment, and neurological and psychiatric disorders. Dendritic spines undergo structural modifications during their lifespan, with properties such as total spine number, dendritic spine size, and morphologically defined subtype altering in response to different processes. Delineating the molecular mechanisms regulating these structural alterations of dendritic spines relies on morphological measurement. This mandates accurate and reproducible dendritic spine analysis to provide experimental evidence. The present study outlines a detailed protocol for dendritic spine quantification and classification using Neurolucida 360 (automatic three-dimensional neuron reconstruction software). This protocol allows for the determination of key dendritic spine properties such as total spine density, spine head volume, and classification into spine subtypes thus enabling effective analysis of dendritic spine structural phenotypes.
Dendritic spines are protrusions of dendrites often comprising the post-synaptic site of glutamatergic synapses1,2. Dendritic spines are of particular interest in the field of synaptic plasticity. Spines are often altered when synaptic strength changes, becoming larger and stronger in long-term synaptic potentiation or smaller and weaker in long-term synaptic depression3,4,5,6,7. Beyond synaptic plasticity, the profile of dendritic spines changes throughout the lifespan. In early development, there is a period of dendritic spine formation and growth, followed by dendritic spine pruning until reaching a steady state8,9,10. In the aging brain, spine loss accompanies brain shrinkage and cognitive decline11. Additionally, many neurological, neurodegenerative, and psychiatric disorders are characterized by aberrant dendritic spines. Multiple brain regions in individuals affected with schizophrenia have fewer dendritic spines, likely resulting from altered synaptic pruning12. Autism spectrum disorders are also characterized by dendritic spine pathologies13. Dendritic spine loss is a hallmark of both Alzheimer's and Parkinson's disease14,15. Given the wide array of research topics encompassing investigations into dendritic spine properties, techniques for accurate spine quantification are of paramount importance.
Staining, i.e., the Golgi method, or labeling neurons via dye filling or expressing fluorescent proteins are common methods for dendritic spine visualization16,17,18. Once visualized, spines can be analyzed with a variety of free and commercially available software clients. The desired output of the analysis is an important factor in determining which software will be of the most use. Fiji is a viable software option for questions centered around dendritic spine density. However, this technique largely relies on time consuming manual counting that can introduce the potential for bias. New plugins such as SpineJ allow for automatic quantification, additionally allowing for more accurate spine neck analysis19. A drawback of these approaches is the loss of a three-dimensional analysis for determining spine volume, as SpineJ is limited to two-dimensional image stacks. Additionally, obtaining spine subtype information becomes challenging via these processes. The four predominant spine subtypes, thin, mushroom, stubby, and filopodia, all connote individual functions and are largely classified via morphology20. Thin spines are characterized by an elongated neck and defined head21. Mushroom spines have a much larger and pronounced spine head22. Stubby spines are short and have little variance between head and neck23. Filopodia are immature spines with a long, thin neck and no obviously observable head24. While classification provides valuable information, spines exist on a continuum of dimensions. Classification into categories is based on ranges of morphological measurements25,26. Manually measuring spines for classification compounds the logistical burden for researchers in this approach.
Other software options focusing specifically on three-dimensional dendritic spine analysis are better suited for investigations into spine volume and subtype properties27,28,29,30,31. Despite the difficulty presented by three-dimensional analysis, such as poor z-plane resolution and smear, these software options allow for reliable three-dimensional reconstruction of dendrites and dendritic spines in a user-guided semi-automated fashion. Automatic classification of identified spines into their subtypes is also a feature present in some of these spine analysis software packages. This can ameliorate concerns of potential workload and experimental bias. Neurolucida 360 is one commercially available software allowing for reliable and reproducible three-dimensional dendritic spine identification and classification32. Here, we present a comprehensive protocol to effectively prepare fixed tissue, acquire images, and ultimately quantify and classify dendritic spines using this software.
All animal procedures followed the US National Institutes of Health Guidelines Using Animals in Intramural Research and were approved by the National Institute of Mental Health Animal Care and Use Committee.
1. Preparation of fixed hippocampal slices
2. High-resolution confocal imaging
3. Dendritic spine quantification
Effectively utilizing this analysis method begins with the selection of dendritic segments for tracing. As described in Figure 1, the ideal dendrites for tracing are not in close proximity to other dendrites. Dendrites running in parallel can result in improperly identifying spines from a neighboring dendrite. Dendrites directly intersecting or running perpendicular in a different z-plane add significant difficulty to accurate dendritic tracing as well. It is also important to note the differences in dendrite thickness. As previously reported, there are key differences in spine density with dendrites of varying thickness36. There can also be differences in the same dendrite with increased distance from the branch point37. Tracing dendrites of the same order and thickness, ideally with similar branch point origins, can control for the existing heterogeneity of dendritic spine density. Identifying the branch point in some preparations may prove unfeasible, but the thickness of the dendrite should always be a controllable factor in dendrite tracing. The accurate tracing of dendritic segments is vital to obtaining accurate results from this analysis. It is necessary to ensure that all points of the traced dendrite are truly within the dendrite. Viewing the three-dimensional dendrite from different directions can assist with this process. As demonstrated in Figure 2A,B, the top-down view shows what appears to be a properly traced dendrite. In the side view; however, numerous points are not located on the dendrite itself. These issues are not present in the side view of Figure 2C. It is also vital to ensure dendrites are properly filled during tracing. A dendrite that is underfilled can result in pieces of dendrites being inappropriately identified as spines. A dendrite that is overfilled can prevent true spines from being identified due to the minimum height threshold. This manual assessment of the user-guided tracing is critical to allow for accurate dendritic spine analysis.
The identification of dendritic spines also requires a user-guided approach. Using the "Detect All" function to set the uniform detector sensitivity threshold is inadequate for numerous reasons. Using the "Detect All" feature is useful for identifying the most blatantly obvious spines, but the filling of these spines must be checked to verify. The identified spines with the initial "Detect All" may be underfilled. To correct this, the identified spine must be individually deleted and then reidentified manually at a higher detector sensitivity (Figure 3A–C). This ensures that the spine is adequately filled. There is substantial heterogeneity in the required detector sensitivity for spines that must be accounted for manually. Increasing the detector sensitivity to detect all may result in overly filled spines, which require manual correction as well (Figure 3D). An additional issue with improper detector sensitivity is the inappropriate creation of a conglomerate spine, one filled dendritic spine that encompasses multiple spines. Two spines in close proximity to each other can be improperly merged into one conglomerate spine (Figure 4A,B). The spine detection software has a "Split" feature, which can be used to separate spines that have been merged by overfilling. The "Split" feature allows for the individual spines to be readily generated from the conglomerate spine (Figure 4C). Accurate dendrite tracing and dendritic spine filling allow for accurate classification into spine subtypes. Spine classification relies on morphology from the filled spines and distance from dendrites, so every step in the process plays a role in the morphological classification (Figure 5).
Due to the necessity for manual selection and thresholding, it is crucial to follow a uniform standard for all analyses. This is especially pertinent if multiple users contribute to data analysis. To ensure that all investigators performing analysis are following the same standard, investigators should compare data from the same traced dendrites. This can reduce the potential for experimenter bias by ensuring that each researcher is identifying spines based on shared, uniform criteria in a blinded fashion. There is also the possibility for bias from a single researcher between days or even on the same day due to fatigue. This should be monitored throughout the process of data analysis. To further ensure the validity of the analysis, comparing initial results to those published in the literature ensures that the protocol is being effectively followed. It is critical to note that this comparison will only be effective if the preparation and parameters are shared. Differences in staining, acquisition of fluorescent signals, order and thickness of dendrites, or brain region can contribute to different results8,36. In the case of missing published results, using multiple researchers to validate spine identification allows for increased confidence in the reliability and reproducibility of the analysis. A supplemental analysis folder has been included in this manuscript. This folder contains files of sample images of dendritic segments, traced dendrites, traced dendrites with identified and classified spines, and data output (Supplementary Table 1, Supplementary File 1, Supplementary File 2, Supplementary File 3, and Supplementary File 4). New users can train on this data set to practice the procedures described in this paper. User-generated results within 10% of the provided sample dataset are considered acceptable for reproducing the standard of analysis. Due to the potentially subjective criteria of a fully filled spine and the need for manual examination of automatically detected spines, variance between and within researchers is a normal part of the analysis. Should the generated results exceed that threshold; however, a side-by-side comparison should be conducted to determine instances of different spine volumes as well as improperly included or excluded spines. The sample dataset can then be reanalyzed until the acceptable threshold is reached.
Figure 1: Selecting dendrites for dendritic spine analysis. (A) 3D-volume display of z-stack confocal images taken from CA1 proximal dendrites in the THY1-YFP transgenic mouse line. Note the heterogeneity of dendrite order with thicker primary dendrites in blue ovals and thinner, secondary and tertiary dendrites in pink ovals. (B) Ideal candidates for dendrite tracing are denoted by green ovals. Note the thickness and limited intersections, overlaps, and proximity to other dendrites. The red oval denotes dendritic segments to be avoided for dendritic tracing due to high intersections, overlaps, and proximity to other dendrites. Thicker, primary dendrites are also not candidates suitable for tracing. Scale bar = 25 µm. Please click here to view a larger version of this figure.
Figure 2: Accurately tracing dendritic segments. (A) 3D-volume display of z-stack confocal images taken from CA1 proximal dendrites in the THY1-YFP transgenic mouse line to be traced via the user-guided directional kernel method. Scale bar = 10 µm. (B) Example of poor dendrite tracing. The dendrite appears to be properly traced in the top-down view. The side view shows the dendrite is improperly filled with points deviating from the dendrite. (C) Example of a proper dendrite tracing. The top-down view appears similar to B, but the side view differs substantially. The dendrite in C is properly traced as indicated by being fully filled with no deviations from the dendrite. Please click here to view a larger version of this figure.
Figure 3: Accurately filling dendritic spines using manual selection. (A) 3D-volume display of z-stack confocal images taken from CA1 proximal dendrites in the THY1-YFP transgenic mouse line of a spine awaiting manual detection. Scale bar = 0.5 µm. (B) Example of an underfilled dendritic spine. There is a substantial fluorescent signal still visible due to incomplete filling. (C) Example of a properly filled dendritic spine. The presence of a "corona" of signal just barely visible around the exterior of the filling is the standard for accurately filling dendritic spines. (D) Example of an overfilled dendritic spine. The detector sensitivity is too high, resulting in an overfilled spine. The filling has gone beyond the borders of the fluorescence and has an almost imperceptible corona. Please click here to view a larger version of this figure.
Figure 4: Splitting conglomerate dendritic spines. (A) 3D-volume display of z-stack confocal images taken from CA1 proximal dendrites in the THY1-YFP transgenic mouse line with two spines in close proximity. Scale bar = 0.15 µm. (B) An example of two independent spines improperly filled as one conglomerate dendritic spine. (C) Following the use of the "Split" feature, the conglomerate spine is split into two distinct properly filled dendritic spines. Please click here to view a larger version of this figure.
Figure 5: Dendritic spine identification and classification into subtypes. (A) 3D-volume display of z-stack confocal images taken from CA1 proximal dendrites in the THY1-YFP transgenic mouse line of a traced dendritic segment isolated for dendritic spine quantification and classification. Scale bar = 5 µm. (B) Traced dendritic segment with all dendritic spines identified and examined to ensure proper filling and splitting. The software arbitrarily assigns colors to identified spines in this step. (C) Classification of all identified dendritic spines into subtypes using defined parameters in the software. Blue = mushroom, yellow = thin, and green = stubby. Filopodia are not present due to the age of this tissue. (D) Representative images of mushroom, thin, and stubby spines unfilled (top) and filled (bottom). Scale bar = 0.3 µm. Please click here to view a larger version of this figure.
Supplementary Figure 1: Accessing the 3D Environment. Z-stack of confocal images viewed in the software interface. 3D Environment navigation from the Trace tab in the main viewer has been highlighted in yellow. Please click here to download this File.
Supplementary Figure 2: Image parameters and orientation settings for 3D Environment. 3D Environment viewer for confocal z-stack images. Parameters in the highlighted Change Image Display tab denoted by yellow arrows are set to Display Image As: 3D Volume and Show Surface As: Max Projection. Move Pivot Point and Reset Orientation are identified by yellow arrows. Please click here to download this File.
Supplementary Figure 3: Dendrite segment tracing. (A) 3D-volume of z-stack confocal images for dendrite tracing. With the tree tab, user-guided, and directional kernels all selected, tracing begins by placing the initial kernel on the dendrite with a left click. (B) Propagation of directional kernels down the dendrite following cursor movement. (C) Left-clicking further down the dendrite fills the directional kernels. (D) Example of directional kernels not populating down dendrite. Instead, a lone kernel is present further down the segment. (E) Left-clicking at the lone kernel fills the dendrite between the two points. Right-clicking ends the tracing. Please click here to download this File.
Supplementary Figure 4: Adjusting points in traced dendrites. (A) Traced dendrite segment pending point adjustment. Dendrite editing requires the "Tree" tab and "Edit" tab to be selected. Both are highlighted in yellow. Dendrite has been selected for editing with a left click. (B) Selecting the points tab, highlighted in yellow, allows for the selection of individual points on the dendrite segment. The green point has a thickness of 1.2 µm. (C) Adjusted point to fill the dendrite more accurately. The new thickness value of the green point is 0.6 µm. Please click here to download this File.
Supplementary Table 1: Sample image analysis results. Please click here to download this File.
Supplementary File 1: Sample image tracings with dendrites and spines.dat Please click here to download this File.
Supplementary File 2: Sample tracings with dendrites.dat Please click here to download this File.
Supplementary File 3: Sample dendrite image file.czi Please click here to download this File.
Supplementary File 4: Sample dendrite image file.jpx Please click here to download this File.
This protocol details the specific steps of sample preparation, imaging, and the process of dendritic spine quantification and classification using three-dimensional reconstruction software. This software is a powerful tool capable of producing robust structural data that contributes to a diverse array of investigations. Throughout the process, there are some critical steps that make this protocol less of a methodological burden and enhance the overall output of the data. The method for labeling dendritic spines is one of the first things researchers should consider before embarking on this protocol. Issues with spine quantification can arise from insufficient labeling methods. Staining for certain proteins expressed at low levels in spines can result in signals that are too low for the software. It must also be noted that bias can be introduced by tracing the brightest fluorescing dendrites. While it is unclear if different fluorescing dendrites have different physiological properties, it is still a limitation of the protocol to consider. Additionally, in some transgenic lines, such as the THY1-YFP-H line, fluorescence in dendritic spines does not appear until around P21. This makes this line unsuitable for investigations into younger developmental time points. Consideration of the method utilized to label spines for the avenue of investigation is not a trivial aspect, as without sufficient fluorescence, the software has diminished usability. Similarly, the image acquisition hardware requires consideration. There are some file types that prove to be less compatible with the analysis software than others. Specifically, ND2 files have been identified as problematic file types for effective use of the software. The software providers recommend conversion to file types such as JPEG2000 should issues occur.
Tissue preparation and image acquisition are also important steps for high-quality spine quantification. Proper fixation, slicing, and mounting of tissues ensures a long-lasting sample with minimal artifacts that can interfere with the data analysis. Imaging the tissue is also not simply a matter of taking z-stack images of the entire brain slice. During imaging, the intent should always be to acquire stacks containing dendrites for spine quantification. An emphasis should be put on acquiring z-stacks incorporating dendrites that will be easy to trace. A thicker z-stack typically results in more background dendrites. This makes it more difficult to effectively trace dendrites with the software. Taking extra time during imaging to find better candidates for tracing will save more time during the spine quantification analysis. Additionally, ensure that the z-stack images incorporate the entire dendrites to be traced. If the dendrites are only partially viewable in the z-stack, dendrite tracing, and spine identification will prove to be difficult and inaccurate due to incomplete 3D rendering.
The process of tracing dendrites and obtaining an accurate profile of identified dendritic spines can be an arduous process. There is a degree of nuance to it. During dendrite tracing the user-guided function can occasionally not function as intended. On occasions, the directional kernels will not populate over a certain segment or begin populating over an undesired segment. One way to circumvent this is to start with a smaller typical process width. This makes the dendrite more detectable in the software, allowing for easier tracing. Should the directional kernel fail to populate entirely, left-clicking on the dendrite will place one manually. It will come to a very small point and not fill the dendrite, but that can be corrected with the thickness adjustment as described in step 3.8 of the protocol. While dendrites can be manually traced should the software prove inadequate, manually tracing spines is not a capability of the software. There will be instances where a clear spine appears to be visible, but no matter how high the detector sensitivity goes, the software will fail to detect it. One thing to check is if the suspected spine is out of range. If it is in range but still not being identified by the software, then this spine will be excluded from the analysis. While this may rarely occur, it is a limitation to consider. As with any analysis requiring thresholding and a manual classification component, there is the possibility of introducing bias. This issue can be further compounded when comparing data generated by multiple users. The semi-automated nature of this analysis seeks to minimize the introduction of this bias, but it is not wholly eliminated. In our lab, a 10% variance between researchers on a standard dataset was reasonable with sufficient practice and training. While efforts have been made to minimize bias, it is still important to consider inter-researcher bias when evaluating data generated through this protocol.
Taking into consideration the minor drawbacks of the software, the output of dendritic spine analysis using this technique is very robust. As previously described, there are a myriad of metrics that can be extrapolated from accurate dendritic segment tracing and spine identification. The ability to obtain spine subtype information provides valuable insight at a deeper level than the basic metrics. This data is important due to the interconnectivity between structure and function. Each spine subtype connotes a function. Thin spines are the predominant subtype undergoing abundant turnover21. Thin spines also have the potential to develop into mushroom spines38. This is consistent with mushroom spines being strongly associated with learning and memory38,39. Stubby spines are additionally believed to be a component of learning, potentially as remnants of mushroom spines40. Filopodia, while not prevalent in many adult tissues, are spine precursors that are of key interest in development41,42. 3D electron microscopy remains the gold standard for the most accurate classification of spine subtypes. While valuable, this technique is limited by arduous manual sorting and classification that are prone to human errors. The semi-automated design of this analysis reduces the instances where subjective bias could be introduced. While there may be drawbacks in terms of the absolute resolution and fluorescence intensity required to make perfect classifications in this protocol, it still provides a methodologically less taxing alternative to 3D electron microscopy and manual classification. Furthermore, full dendrite analysis of multiple dendritic branches from a wider swatch of brain regions is possible using the analysis outlined in this work. This is not the case with electron microscopy. Through the use of this protocol, it is possible to address structure-centered questions in multiple disciplines, including but not limited to synaptic plasticity, development, and neurological and psychiatric disorders, in a reliable and reproducible manner.
The authors have nothing to disclose.
We would like to acknowledge Carolyn Smith, Sarah Williams Avram, Ted Usdin, and the NIMH SNIR for technical assistance. We would additionally like to acknowledge the Colgate University Bethesda Biomedical Research Study Group. This work is supported by the NIMH Intramural Program (1ZIAMH002881 to Z.L.).
518F Immersion Oil | Zeiss | 444960-0000-000 | |
Cryostat | Leica | CM3050S | For slice preparation |
Fine Forceps | FST | 11150-10 | |
Hemostat Forceps | FST | 13020-12 | |
Large Surgical Scissors | FST | 14002-16 | |
LSM 880 Confocal Microscope | Zeiss | LSM 880 | |
Microscope Cover Glass | Fisherbrand | 12-541-035 | |
Mini-Peristaltic Pump II | Harvard Apparatus | 70-2027 | For perfusions |
Neurolucida 360 | MBF Bioscience | v2022.1.1 | Spine Analysis Software |
Neurolucida Explorer | MBF Bioscience | v2022.1.1 | Spine Analysis Software |
OCT Compound | Sakura Finetek | 4583 | For cryostat sectioning |
Paraformaldehyde (37%) | Fisherbrand | F79-1 | |
Plan-Apochromat 63x/1.40 Oil DIC | Zeiss | 440762-9904-000 | |
Scalpel Blade | FST | 10022-00 | |
Small Surgical Scissors | FST | 14060-09 | |
Spatula | FST | 10091-12 | |
Sucrose | FIsherbrand | S5-500 | |
Superfrost Plus Microslides | Diagger | ES4951+ | |
Vectashield HardSet Mounting Medium | Vector Laboratories | H-1400-10 |