Skip to main content

ORIGINAL RESEARCH article

Front. Synaptic Neurosci., 15 February 2022
This article is part of the Research Topic Molecular Nanomachines of the Presynaptic Terminal, Volume II View all 12 articles

Correlative Live-Cell and Super-Resolution Imaging to Link Presynaptic Molecular Organisation With Function

Updated
  • 1Centre for Developmental Neurobiology, Institute of Psychiatry, Psychology and Neuroscience, King’s College London, London, United Kingdom
  • 2MRC Centre for Neurodevelopmental Disorders, Institute of Psychiatry, Psychology and Neuroscience, King’s College London, London, United Kingdom

Information transfer at synapses occurs when vesicles fuse with the plasma membrane to release neurotransmitters, which then bind to receptors at the postsynaptic membrane. The process of neurotransmitter release varies dramatically between different synapses, but little is known about how this heterogeneity emerges. The development of super-resolution microscopy has revealed that synaptic proteins are precisely organised within and between the two parts of the synapse and that this precise spatiotemporal organisation fine-tunes neurotransmission. However, it remains unclear if variability in release probability could be attributed to the nanoscale organisation of one or several proteins of the release machinery. To begin to address this question, we have developed a pipeline for correlative functional and super-resolution microscopy, taking advantage of recent technological advancements enabling multicolour imaging. Here we demonstrate the combination of live imaging of SypHy-RGECO, a unique dual reporter that simultaneously measures presynaptic calcium influx and neurotransmitter release, with post hoc immunolabelling and multicolour single molecule localisation microscopy, to investigate the structure-function relationship at individual presynaptic boutons.

Introduction

Synapses are remarkably heterogenous subcellular compartments, formed of a presynaptic terminal and a postsynaptic element that vary in structure, function and molecular composition. At the presynaptic terminal, calcium influx through voltage gated calcium channels (VGCCs) triggers the fusion of a neurotransmitter-filled synaptic vesicle (SV) with the plasma membrane at the active zone (AZ), a specialised site containing the necessary machinery for vesicle tethering, docking and release (Südhof, 2012). Individual synapses display a high degree of heterogeneity in all stages of this process, differing in their expression of calcium channel subtypes (Luebke et al., 1993; Wheeler et al., 1994; Reid et al., 1997; Gasparini et al., 2001), the level of calcium influx (Ermolyuk et al., 2012) and the sensitivity of release to calcium (Ermolyuk et al., 2012; Jackson and Burrone, 2016; Rebola et al., 2019). Most importantly, neurotransmitter release is a probabilistic process, with SV fusion occurring with a release probability (Pr) that varies widely between synapses (Murthy et al., 1997). Presynaptic boutons are also heterogenous in their structure, such as volume and shape, as well as in their sub-synaptic features, such as AZ area and the number of neurotransmitter-filled vesicles, including their distribution in different vesicle pools (Schikorski and Stevens, 1997; Holderith et al., 2012). Although several of these aspects of presynaptic structure correlate with functional measures, our understanding of this important relationship remains incomplete, in part due to the technical challenges associated with studying both structure and function at the level of the individual synapse. The small size of presynaptic structures, below the spatial resolution limit of light microscopy, requires the use of electron or super-resolution microscopy, whilst the small and rapid events associated with neurotransmitter release require highly sensitive probes with good temporal resolution. It is therefore necessary to combine imaging modalities to uncover how structure informs function at the chemical synapse.

Methods to Study Synapse Function

Our fundamental understanding of presynaptic function comes from classical electrophysiological studies that identified the quantal nature of neurotransmitter release (del Castillo and Katz, 1954) and its relationship with calcium influx (Dodge and Rahamimoff, 1967). Whilst electrophysiological techniques have since been refined to enable recording from individual presynaptic terminals, including small central boutons (Novak et al., 2013), imaging techniques remain the simplest approach to study synapses in the large numbers required to understand their heterogeneity. For this, a range of chemical indicators and genetically encoded reporters have been developed to measure two key aspects of presynaptic function, calcium influx and neurotransmitter release (reviewed in Wang et al., 2019). More recently, the development of genetically encoded voltage indicators has also allowed imaging of the shape of the AP along axons and at individual boutons (Hoppa et al., 2014; Sabater et al., 2021).

Reporters of neurotransmitter release consist of a pH-sensitive form of GFP (pHluorin), typically fused to the luminal domain of a synaptic vesicle protein such as VAMP2 (Miesenböck et al., 1998), synaptophysin (Granseth et al., 2006), or VGLUT1 (Voglmaier et al., 2006). The pHluorin fluorescence is quenched at the acidic pH inside the vesicle and undergoes a 20-fold increase upon vesicle fusion and exposure to the extracellular medium at pH7.4 (Sankaranarayanan et al., 2000). Variants of these probes now include the addition of an invariant fluorophore (tdimer2) for normalisation of the pHluorin signal (Ratio-sypHy, Rose et al., 2013) as well as red-shifted pHluorin analogues for multicolour imaging (Li and Tsien, 2012; Liu et al., 2021). More recently, fluorescent sensors that bind glutamate or GABA have been constructed (Marvin et al., 2013, 2018, 2019) to directly sense neurotransmitter release.

The most widely used genetically encoded calcium indicators are the GCaMP family (Nakai et al., 2001), which have undergone numerous iterations to improve sensitivity and kinetics. Red-shifted calcium indicators have also been produced, such as RGECO and RCaMP1, from which jRCaMP1a and jRGECO1a have been developed (reviewed in Wang et al., 2019). Many of these probes have been targeted to specific subcellular compartments, including the presynaptic terminal (Dreosti et al., 2009; Walker et al., 2013). To monitor both presynaptic calcium influx and neurotransmitter release simultaneously we constructed SypHy-RGECO (Jackson and Burrone, 2016), a dual sensor probe which also provides an optical measure of the calcium dependence of release with single synapse resolution.

This range of tools has been essential in demonstrating the heterogeneity of synapse function, but to understand the underlying structural features that govern it, they must be combined with tools that provide information on synapse morphology and molecular organisation.

Methods to Study Synapse Structure

Electron microscopy (EM) provided the first insights into synapse structure, revealing the existence and arrangement of SVs within a presynaptic terminal, the electron dense AZ with tethered vesicles poised for release and the opposing postsynaptic density (PSD) (Harris and Weinberg, 2012), now known to precisely cluster neurotransmitter receptors (Nair et al., 2013). Although unparalleled in its spatial resolution, EM does not provide information on the distribution of molecules unless combined with immunogold labelling which remains challenging to achieve with high labelling efficiency. Cryo-electron tomography holds great promise to study the structure of the synapse and its proteins in their native state but remains a specialised and experimentally demanding technique.

Over the last two decades, super-resolution microscopy (SRM) has emerged as a tool that allows examination of protein organisation at the nanoscale using standard fluorescent labelling approaches (reviewed in Sigal et al., 2018). Whilst it does not reach the spatial resolution of EM, it is a vast improvement over conventional light microscopy, provides information on the nanoscale organisation of molecules and is more amenable to higher throughput imaging. SRM has shown that several key AZ proteins such as Bassoon (Dani et al., 2010; Glebov et al., 2017), RIM1/2 (Tang et al., 2016; Haas et al., 2018) and Munc13 (Sakamoto et al., 2018) are arranged into subsynaptic domains (SSDs), as are postsynaptic proteins including PSD95 (Fukata et al., 2013; MacGillavry et al., 2013; Nair et al., 2013) and AMPA receptors (Nair et al., 2013). VGCCS have been shown to be highly dynamic at the presynaptic membrane, but transiently stabilise into SSDs at the AZ through protein interactions (Schneider et al., 2015; Heck et al., 2019). Importantly, SRM has brought the pre- and post-synaptic compartments together at the nanoscale by uncovering a novel arrangement in which RIM1/2 and PSD-95 SSDs align across the synaptic cleft, creating a transsynaptic nanocolumn for the efficient transfer of information across the cleft (Tang et al., 2016; Haas et al., 2018). In an intriguing twist to this arrangement, the nanoscale organisation of synaptic proteins was also shown to be plastic, suggesting an important role of subsynaptic domains in regulating synaptic function (Tang et al., 2016; Glebov et al., 2017; Hruska et al., 2018; Heck et al., 2019; Compans et al., 2021). However, the link between nanoscale organisation of synaptic molecules and synapse function is only just beginning to be understood.

Studying the Structure-Function Relationship at the Synapse

Correlative light and electron microscopy studies have revealed that release probability is tightly coupled to the number of docked vesicles, which constitutes the readily releasable pool (RRP) (Murthy et al., 2001; Branco et al., 2010), and that both properties correlate with area of the AZ (Holderith et al., 2012). Whilst these correlations are maintained within synapses from a single type of neuron, synapses from different neuron types but with similar ultrastructural features, such as number of docked vesicles, display different release probabilities (Xu-Friedman et al., 2001).

Neurotransmitter release is driven by calcium influx such that the average release probability of a vesicle in the RRP is correlated with the size of the calcium transient (Ermolyuk et al., 2012). The number of VGCCs also scales linearly with AZ area and release probability in some synapse types (Holderith et al., 2012; Miki et al., 2017). However, weak cerebellar granule cell synapses express greater numbers of Cav2.1 channels than the stronger stellate cell synapses, but these channels are located at a greater distance from SV release sites and are therefore less effective at initiating neurotransmitter release (Rebola et al., 2019). This relationship can be further complicated by the expression of different calcium channel subtypes (Wheeler et al., 1994; Reid et al., 1997; Gasparini et al., 2001) and splice isoforms (reviewed in Lipscombe et al., 2013) that display different dynamics at the plasma membrane as well as distinct biophysical properties, which affect Pr (Heck et al., 2019). In addition, molecular crowding may affect the recruitment and retention of calcium channels in the plasma membrane. For example, blockade of neuronal activity causes a decrease in the density of Bassoon clusters at the active zone and the preferential recruitment of Cav2.1 to synapses with lower density (Glebov et al., 2017). Taken together it is clear that both the number of synaptic proteins and their spatial relationships are important in determining synaptic function.

A key spatial relationship at the synapse is the alignment of pre and postsynaptic protein clusters across the synaptic cleft in transsynaptic nanocolumns, which define the SSDs at which neurotransmitter release will maximally activate postsynaptic receptors (Nair et al., 2013; Tang et al., 2016; Haas et al., 2018). Until recently it was thought that the AZ constituted a single release site, however precise localisation of vesicle fusion using vGlut-pHluorin has revealed that there are multiple hotspots for release at each AZ, ranging from 5 to 14 (Maschi and Klyachko, 2017). Combining this with live SRM has shown that these hotspots are marked by dense nanoclusters of RIM1/2 (Tang et al., 2016), a protein whose levels at the AZ have been intimately linked to Pr (Holderith et al., 2012), and constitute distinct release sites. Munc13-1 and syntaxin-1 are similarly organised into SSDs, and their numbers also correlate with the number of release sites as well as the number of vesicles in the RRP (Sakamoto et al., 2018). However, in contrast to previous studies (Branco et al., 2010; Holderith et al., 2012), the number of Munc13-1 molecules, and consequently the number of release sites and RRP size, did not correlate with Pr at individual synapses (Sakamoto et al., 2018), suggesting that our understanding of the relationship between the molecular composition of release sites and release probability remains incomplete.

Many studies relating synaptic nanoscale structure with function have done so using indirect methods (Glebov et al., 2017) and only a few have attempted to combine SRM techniques and functional imaging at the same synapses (Tang et al., 2016; Sakamoto et al., 2018; Holderith et al., 2020), despite the insight this type of correlation can provide. Here we present a correlative approach combining live imaging of SypHy-RGECO, a genetically encoded reporter that simultaneously measures calcium influx and neurotransmitter release (Jackson and Burrone, 2016), with multicolour 3D-dSTORM using spectral demixing (Testa et al., 2010; Platonova et al., 2015). Using the well-characterised presynaptic proteins Bassoon and Cav2.1, we provide proof of principle that both nanoscale organisation and functional measures can be investigated at the level of individual boutons.

Materials and Methods

An overview of the pipeline for functional and 3D-dSTORM correlative imaging is shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Functional and 3D-dSTORM correlative imaging pipeline. (1–3) Preparation of neurons for functional imaging of a presynaptic calcium and vesicle release sensor (SypHy-RGECO). (4) Fixation and immunolabelling of neurons. (5–6) Relocation and super-resolution (dSTORM) imaging of SypHy-RGECO + ve neurons. (7–8) Registration between imaging modalities, analysis of calcium influx (RGECO) and neurotransmitter release (SypHy) and nanoscale organisation of synaptic proteins (dSTORM), followed by single synapse correlation between modalities.

Step 1: Preparation of Hippocampal Neurons on Gridded Dishes

Hippocampi were dissected from E18.5 Wistar rat embryos (Charles River Laboratories, United Kingdom), treated with trypsin (Worthington, United Kingdom) at 0.5 mg/ml and triturated with fire-polished Pasteur pipettes. Dissociated cells (80–100k) were plated onto 35 mm low profile Grid500 μ-dishes (Ibidi, Germany) pre-treated with 100 μg/mL poly-D-lysine (Sigma, United Kingdom) and 10 μg/mL laminin (Thermo Fisher Scientific, United Kingdom). Neurons were maintained in Neurobasal media supplemented with 2% B27, 2% foetal bovine serum, 1% glutamax, and 1% penicillin/streptomycin (all Thermo Fisher Scientific, United Kingdom) in a 37°C humidified incubator with 5% CO2. After 3 days in vitro (DIV) the media was exchanged for serum and antibiotic-free media.

Step 2: Transfection With SypHy-RGECO

At 7DIV, neurons were transfected with CAMKII:SypHy-RGECO [Jackson and Burrone, 2016, Addgene plasmid #84078, subcloned under a minimal CAMKII promoter in the vector backbone of pAAV-CW3SL-EGFP, a gift from Bong-Kiun Kaang (Choi et al., 2014), Addgene plasmid #61463] using Effectene, following the manufacturer’s protocol (QIAGEN, United Kingdom). Expression of SypHy-RGECO under the minimal CAMKII promoter should be restricted to excitatory neurons.

Step 3: Live Imaging of SypHy-RGECO

Neurons were imaged live at DIV17-21 on an inverted Olympus IX71 microscope equipped with a 60x/1.42 NA oil objective, with a heated stage and objective to maintain the bath temperature at 32 ± 2°C. The growth media was replaced with HEPES buffered saline (HBS; 139 mM NaCl, 2.5 mM KCl, 10 mM HEPES, 10 mM D-Glucose, 2 mM CaCl2, 1.3 mM MgCl2; pH7.3 at 32°C, 290 mOsmol) supplemented with 10 μM NBQX, 25 μM APV and 10 μM Gabazine (Tocris, United Kingdom) to prevent recurrent activity. Images were acquired with an ORCA-Flash4.0 V2 C11440-22CU scientific CMOS camera (Hamamatsu Photonics, Japan) using HCImageLive Software (Hamamatsu Photonics) and the camera was cooled to approximately −20°C with an Exos2 water cooling system (Koolance, United States). The triggering of LEDs, camera acquisition and stimulation were controlled externally through Clampex software (pCLAMP10, Molecular Devices, United States) and a Master-8 programmable stimulator (A.M.P.I, Israel).

A GFP-positive cell was identified, and an extracellular parallel bipolar electrode (FHC, United States) was positioned with the tips on either side of the soma using a PatchStar Motorised Micromanipulator (Scientifica, United Kingdom). A 10-action potential (AP) stimulation at 20 Hz was delivered by an SD9 stimulator (Grass Technologies, United States) using 10 pulses of 9.5–10 V for 1 ms, a stimulus previously shown to elicit a single AP per pulse (Sabater et al., 2021). Imaging was performed with a dual band-pass filter set optimised for EGFP and mCherry (#59022, Chroma Technology, United States) and with LED excitation light sources of 470 and 585 nm (CoolLED, United Kingdom) set at 70 and 100% power, for pHluorin and RGECO fluorophores, respectively. The full frame was acquired with 4 × 4 binning (pixel size = 433 nm), alternating 20 ms exposures of the 470 and 585 nm LEDs with a 20 ms interval between frames, resulting in a final acquisition rate of 12.5 Hz for each channel.

Step 4: Immunofluorescence

Immediately after live imaging, neurons were fixed with 4 or 1% PFA (for Cav2.1 labelling) for 10 min at room temperature (RT). Neurons were washed in PBS 3 times and 50 mM NH4Cl in PBS was applied as a quenching agent, followed by a further 3 washes in PBS. Dishes were stored in PBS at 4°C in the dark until staining. Neurons were permeabilised with 0.1% TritonX-100 in PBS for 10 min, washed three times in PBS and placed in blocking solution (2% BSA in PBS) for 60 min at room temperature (RT). Cells were incubated with primary antibodies diluted in blocking solution for 60 min at RT then washed three times in blocking solution for 5 min each. Secondary antibodies and nanobodies in blocking solution were applied for 60 min at RT and cells were washed three times in block and 3 times in PBS for 5 min each. Cells were kept in PBS at 4°C in the dark until STORM imaging.

Step 5–6: Relocate Neurons for Multicolour 3D-dSTORM Imaging

dSTORM imaging was performed using a spectral-demixing SAFeRedSTORM module (Abbelight, France) mounted on an Olympus IX3 equipped with an oil-immersion objective (100 × 1.49NA oil immersion, Olympus, United Kingdom) and fibre-coupled 642 nm laser (450 mW, Errol). Fluorescent signal was collected with two ORCA-fusion sCMOS cameras (Hamamatsu). A longpass dichroic beam splitter (700 nm; Chroma Technology) was used to split the emission light on the two cameras. 3D-dSTORM imaging was performed by using cylindrical lenses placed before each camera.

The imaging buffer contained 1 mL Glucose Base (10% Glucose, 0.1% Glycerol), 125 μL of Enzyme Buffer (42 μg/mL Catalase (C100, Sigma), 1 mg/mL Glucose oxidase (G2133, Sigma), 50% Glycerol, 20 mM Tris-HCl pH7.5, 4 mM TCEP (C4706, Sigma), 25 mM KCl) and 125 μL of 1 M MEA-HCl (pH∼8 with NaOH, M6500, Sigma). This allowed a final concentration of 4 μg/mL of Catalase and 100 μg/mL of Glucose oxidase, and a final concentration of 100 mM of reducing agent MEA. The final pH was adjusted to fall between pH7.6 and pH8. To limit the contact of the imaging buffer with oxygen, the ibidi dish well was sealed with a 25 mm glass coverslip.

Image acquisition was driven by NEO software (Abbelight). The image stack contained 60,000 frames of a selected ROI of 512 × 512 pixels (pixel size = 97 nm). Cross-correlation was used to correct for lateral drifts. Super-resolution images (.tiff) with a pixel size of 10 nm and localisation files (.csv) were obtained using NEO_analysis software (Abbelight).

Step 7: Image Registration and Synapse Correlation

A maximum z-projection of the live SypHy channel and the reconstructed dSTORM image of the amplified GFP signal were registered using landmark selection and affine transformation in Matlab (Mathworks, script available at https://github.com/jburrone/Jackson-Compans-and-Burrone). Individual synapses that were clearly defined in both images and well aligned in the registration were manually selected and included in correlative analysis.

Step 8: Analysis of Structure and Function

8a: Analysis of Live Images

Images saved in CXD format were converted to TIFF files and split into two separate channels using FIJI (ImageJ). Images were analysed using custom written Matlab scripts (Mathworks, script available upon request). Regions of interest (ROIs, 5 × 5 pixels) were selected for each punctum of SypHy fluorescence and a 15 × 15 pixel ROI was used to calculate background fluorescence. Mean background-subtracted fluorescence intensity values were calculated for each ROI in both SypHy(G) and RGECO(R) channels. Traces were smoothed by averaging over a sliding window of 4 frames. Baseline fluorescence (G0 and R0) was measured as the mean of 15 frames prior to the stimulus. ΔG and ΔR values were calculated by the change in signal intensity from the baseline, with the peak responses defined as the maximum ΔG and ΔR within 15 frames of the stimulus. Puncta in which the ΔR response to the 10AP at 20 Hz stimulus was greater than three times the standard deviation of the baseline were considered responding synapses and were included in further analysis, regardless of whether there was a measurable ΔG response. ΔG and ΔR responses were both normalised to the baseline R0 fluorescence, which represents the overall levels of reporter at the synapse.

8b: Multicolour 3D-dSTORM Analysis

Point Clouds Analyst software (PoCA, Levet et al., 2015, 2019) was used to identify SypHy-RGECO, Cav2.1 and Bassoon clusters from localised molecule coordinates. First a Voronoi diagram was applied to draw 3D polygons of various sizes centred on the localised molecules for each colour separately. Automatic segmentation was performed to detect clusters for each protein (colour1: SypHy-RGECO to identify presynaptic boutons and colour2: Bassoon or Cav2.1). Clusters for each colour were thresholded based on their density (δi1 for colour1 and δi2 for colour2) with respect to the average density of the dataset for each colour, such that δi1 ≥ 1δd and δi2 ≥ 1δd. All selected neighbouring molecules within a maximum distance between neighbours of 200 and 100 nm were considered as forming a cluster when having a minimum number of 200 and 20 molecules for colour1 and colour2, respectively. The colocalisation between SypHy-RGECO and Bassoon or Cav2.1 was computed from the overlapping clusters. We found that some colour2 clusters (Bassoon or Cav2.1) were not automatically detected as colocalising with colour 1 (SypHy-RGECO) clusters even though they clearly formed part of the same synapse. There are multiple reasons for this, including the segmentation thresholds used, which limit the extent of the overlap between colours, and the fact that SypHy-RGECO labels synaptic vesicles and does not cover the entirety of the presynaptic volume. We therefore manually verified and corrected for any missed colocalisations, then extracted cluster volume and density of molecules for colour1 and colour2 from all synapses. Synapses with no colour2 cluster or ambiguous clusters that could not be manually verified were excluded.

8c: Correlation of Live and dSTORM Analysis

Functional responses (ΔG/R, ΔR/R and the ratio ΔG/ΔR) and structural parameters (number of detections and cluster volume and density) were analysed for correlation using Matlab (Mathworks, United States) and Prism 9 (GraphPad, United States), across all synapses selected during image registration (step 7) that passed the inclusion criteria in both imaging modalities (steps 8a and b).

Results

The first imaging step in our pipeline (see section “Materials and Methods” and Figure 1) is to simultaneously measure calcium influx and neurotransmitter release by performing live imaging of SypHy-RGECO, a probe that localises to presynaptic boutons, seen as clear puncta along the axons of hippocampal neurons (Figure 2A). Axonal regions were imaged whilst the soma of the transfected cell was stimulated with a parallel bipolar electrode using 10 pulses that induce single APs (Sabater et al., 2021), a stimulus which lies within the linear range for both reporters (Jackson and Burrone, 2016). Puncta in which the change in RGECO fluorescence was greater than 3S.D. of the baseline were considered functional synapses and included in later analysis (magenta traces, Figure 2A). Using the same criterion for sypHy responses, almost a quarter of synapses (23.27%) did not display measurable neurotransmitter release to a 10AP stimulus despite experiencing a calcium influx, suggesting they had a low release probability (white vs. green traces, Figure 2A). When neurotransmitter release could be quantified, individual synaptic ΔG/R and ΔR/R responses exhibited a wide range of amplitudes, which were positively correlated (Figure 2B). However, due to the spread of points within this correlation, individual synapses with the same level of calcium influx could display different levels of neurotransmitter release (Figure 2C). This can be quantified using the ratio ΔG/ΔR, a measure of the calcium-dependence of release (Figure 2D).

FIGURE 2
www.frontiersin.org

Figure 2. SypHy-RGECO responses to 10AP stimuli. (A) Image of an axon transfected with SypHy-RGECO (green channel (SypHy) shown). Representative traces of SypHy (ΔG/R, green/white) and RGECO (ΔR/R, magenta) responses to 10AP stimuli are shown for boutons with a detectable calcium influx (peak response greater than 3SD of the baseline), associated with SV release (green traces) or failure to release (bottom synapse, white trace) assessed using the same threshold. White boxes indicate selected boutons, all traces shown on the same scale, white arrows indicate stimulation. (B) Cells were stimulated with a 10AP 20 Hz stimulus. The change in fluorescence from baseline (ΔG and ΔR) was normalised to the basal RGECO fluorescence (R) for each synapse. ΔG/R and ΔR/R responses from all synapses show a positive correlation (Spearman’s rank correlation, r = 0.633, p = < 0.0001, n = 211 synapses). (C) SypHy (left) and RGECO (right) traces from 3 synapses (coloured points in B) with a similar amount of calcium influx (RGECO) in response to a 10AP stimulus but a different level of neurotransmitter release (SypHy). Black arrows indicate stimulation. (D) Distribution of the ratio ΔG/ΔR indicating the calcium dependence of release (for the synapses shown in B).

To relocate the same synapses after fixation and staining, we used both their grid position and axonal morphology. We first amplified the SypHy signal using an anti-GFP antibody and an AF488-coupled secondary antibody and labelled endogenous Bassoon with a primary and CF680-coupled secondary antibody suitable for dSTORM imaging. SypHy-RGECO + ve boutons were identified and a widefield image was acquired in the green channel. 3D-dSTORM was then performed for CF680 and the super-resolved image of Bassoon was reconstructed and overlaid on the diffracted GFP image (Figure 3A). Due to the low resolution of the GFP image, chromatic aberrations and the high density of axons in the culture, identifying colocalisation between Bassoon subsynaptic domains (SSDs) and SypHy-RGECO + ve boutons was difficult. To overcome these limits and avoid misidentification, we took advantage of spectral demixing dSTORM (SD-dSTORM) to perform multicolour super-resolution imaging. SD-dSTORM uses fluorophores excitable with a single wavelength (here 642 nm) but with a small shift in their emission spectrum (Figure 3B). Using a long-pass dichroic to split the emission of fluorescence from single emitters onto two cameras, a photon ratio can be calculated to assign the detection to one of the fluorophores (Testa et al., 2010; Platonova et al., 2015; Figures 3C,D and Supplementary Movie 1). As SD-dSTORM uses fluorophores with similar emission wavelengths, chromatic aberrations are negligible compared to more classical multicolour SRM that requires further image processing. Using this technique increased the resolution of SypHy-RGECO + boutons and improved identification of their associated Bassoon SSDs (Figure 3E).

FIGURE 3
www.frontiersin.org

Figure 3. Multicolour 3D-dSTORM to investigate synaptic protein nanoscale organisation. (A) Left panel: schematic of two presynaptic bouton labelled for active zone protein Bassoon resolved with dSTORM. Only one of the two boutons expresses SypHy-RGECO (zoomed synaptic vesicle). The green area represents the diffracted signal from GFP staining imaged with conventional widefield fluorescence microscopy. Right panel: representative example of a widefield image of a SypHy-RGECO + ve axon immunolabelled with AF488 (green) overlaid with a dSTORM image of endogenous Bassoon (magenta) immunolabelled with CF680. Zoomed examples of 3 boutons are shown (square box). (B) Emission spectra of AF647 (light red) and CF680 (dark red), and transmission profile of D700 dichroic (black). (C) Example images of single emitters obtained from the 2 cameras on our SD-dSTORM system and zoomed examples of 2 single emitters. AF647 intensity is stronger on camera 1 than camera 2. CF680 intensity is stronger on camera 2 than camera 1. (D) Distribution of single emitters according to their photon ratio [Icam1/(Icam2 + Icam1)] allows separation of AF647 and CF680 detections. (E) Left panel: schematic of two presynaptic bouton labelled for Bassoon resolved with dSTORM. Only one of the two boutons expresses SypHy-RGECO (zoomed synaptic vesicle). The green area represents the super-resolved signal from GFP staining imaged with SD-dSTORM. Right panel: representative example of multicolour 3D-dSTORM of a SypHy-RGECO + axon immunolabelled with AF647 (green) and endogenous Bassoon (magenta) immunolabelled with CF680. Zoomed examples of 3 boutons are shown (square box).

To further enhance the localisation precision of the SypHy-RGECO + vesicular pool, we used GFP nanobodies directly coupled with AF647. As expected, due to the reduced distance between SypHy and AF647 (Früh et al., 2021), we observed a net reduction in SypHy cluster volume with nanobody labelling compared to classical primary/secondary staining (Supplementary Figures 1A–C). However, after 1% PFA fixation required for labelling Cav2.1, nanobody staining appeared weaker and did not easily allow relocation of SypHy-RGECO + ve neurons for dSTORM imaging, correct registration of image modalities or confident identification of presynaptic boutons (Supplementary Figure 1A, top panel). When compared to nanobody staining following 4% PFA fixation, used for Bassoon labelling, the staining following 1% PFA fixation was less reliable, as shown by the decreased number of detections of AF647 per SypHy-RGECO cluster (Supplementary Figures 1A,D). Since under these conditions the GFP nanobody was likely undersampling the labelling of SypHy, we continued to use the classical primary and secondary antibody labelling approach, which appeared to be less sensitive to fixing conditions and improved the correlation between the number of SypHy-RGECO detections in dSTORM and the RGECO baseline (R0) in live imaging (Supplementary Figures 1E–G). Whilst the primary aim of super-resolving SypHy was to precisely identify live-imaged boutons and their associated proteins of interest, we were also able to investigate the relationship between presynaptic functional properties and SypHy-RGECO detections in neurons labelled with the polyclonal antibody. As this reporter is mostly expressed on SVs, the number and volume of the SypHy-RGECO cluster will be proportional to the total pool of vesicles. We did not observe any correlation between SV release (ΔG/R) or calcium influx (ΔR/R) and either the number or volume of SypHy detections (Supplementary Figure 2).

Two colour super-resolution images of either Cav2.1 or Bassoon together with SypHy-RGECO showed that both AZ proteins are organised in sub-synaptic domains (SSDs), such that a single SypHy-RGECO cluster, representing an individual bouton, could contain anywhere from 1 to 6 Cav2.1 (mean = 1.6, Figure 4A) or Bassoon SSDs (mean = 1.7, Figure 4B). Point Clouds Analyst, a method for 3D two-colour segmentation based on single molecule detection densities (PoCA, Levet et al., 2015, 2019), was used to identify SypHy-RGECO clusters and their associated Bassoon or Cav2.1 SSDs. In most cases, the identification of Bassoon or Cav2.1 SSDs within a SypHy-RGECO cluster was straight forward (Supplementary Figure 3C). In a few cases, segmentation limits meant that associated SSDs were not automatically identified (Supplementary Figures 3D,E). Where possible, these were manually corrected, otherwise synapses were excluded from further analysis. This analysis showed that synapses with more than one Cav2.1 SSD had bigger SypHy cluster volumes (Figure 4C), a difference that was even more noticeable for synapses with multiple Bassoon SSDs (Figure 4D). Interestingly, in synapses with multiple SSDs, the individual clusters of both Cav2.1 and Bassoon display a slight but significant decrease in density compared to those from synapses with single SSDs (Supplementary Figure 4). Together these observations suggest that synapses can use multiple strategies for protein organisation such as forming one or multiple SSDs with different densities of crowding. These types of observation would not be possible using conventional fluorescence microscopy, highlighting the advantage of using SRM to investigate synaptic protein organisation.

FIGURE 4
www.frontiersin.org

Figure 4. Nanoscale organisation of Cav2.1 and Bassoon in SypHy-RGECO + ve boutons. (A,B) Examples of SypHy-RGECO + ve boutons labelled for SypHy-AF647 (left panels) and Cav2.1-CF680 (A: middle panels) or Bassoon-CF680 (B, middle panels) imaged with SD-dSTORM. Merged images of SypHy (green) and Cav2.1 (A, magenta) or Bassoon (B, magenta) are shown in the right panels. White arrows indicate Cav2.1 or Bassoon subsynaptic domains (SSD). Synapses can contain a single (top panels) or multiple SSDs (lower panels). (C) Cumulative distribution of SypHy cluster volume labelled with polyclonal anti-GFP and AF647-coupled secondary antibodies after 1%PFA fixation, for synapses with single or multiple Cav2.1 SSDs (Kolmogorov-Smirnov test, p = 0.025, nsingle = 59, nmulti = 38). (D) Cumulative distribution of SypHy cluster volume labelled with AF647-coupled anti-GFP nanobody after 4%PFA fixation, for synapses with single or multiple Bassoon SSDs (Kolmogorov-Smirnov test, p = 0.004, nsingle = 50, nmulti = 42). *p < 0.05 and **p < 0.01.

As VGCCs are directly involved in neurotransmitter release, we next investigated the potential link between Cav2.1 nanoscale organisation, one major subtype of VGCC found at synaptic boutons, and SV release at individual synapses. The number of Cav2.1 SSDs varied considerably across boutons, ranging from 1 to 6 SSDs, with around 60% of analysed synapses containing only a single Cav2.1 SSD (Figure 5A). We separated synapses into two groups according to the presence of a single or multiple Cav2.1 SSDs. For both groups, we observed a similar positive correlation between ΔG/R and ΔR/R (Figure 5B). Within this correlation, the spread of the responses indicates a variation in the calcium-dependence of release between synapses, quantified by the ratio ΔG/ΔR. We compared the distribution of ΔG/ΔR between synapses with single or multiple Cav2.1 SSDs and found no significant difference (Figure 5C). As we did not observe any correlation between Cav2.1 SSD properties and the level of calcium influx or neurotransmitter release (Supplementary Figure 5), we next split synapses based on their release properties—those with detectable levels of neurotransmitter release in response to a 10AP stimulus and those that do not exhibit any detectable release despite showing robust Ca2+ influx. No difference in the number of Cav2.1 SSDs, nor in SSD properties were observed between synapses which released and those which did not (Figures 5D–F). Our findings indicate that Cav2.1 abundance or spatial distribution are not a good predictor for measures of presynaptic function.

FIGURE 5
www.frontiersin.org

Figure 5. Effect of Cav2.1 nanoscale organisation on neurotransmitter release. (A) Distribution of the number of Cav2.1 SSDs per synapse. (B) ΔG/R and ΔR/R responses from synapses post hoc labelled for Cav2.1 and identified as containing a single (blue) or multiple Cav2.1 SSDs (black) show a positive correlation (Spearman’s rank correlation, rsingle = 0.439, p < 0.0001, nsingle = 74 synapses; rmulti = 0.303, p < 0.0001, nmulti = 47 synapses). (C) Distribution of the ratio ΔG/ΔR for the two group of synapses. Distributions are not significantly different (Kolmogorov-Smirnov test, p = 0.657, nsingle = 74, nmulti = 47). (D) Percentage of synapses with single or multiple Cav2.1 SSDs in synapses without detectable release (left, n = 38) or with release (right, n = 121). White lines indicate separation of multiple SSDs by number. (E) Cumulative frequency distribution of the mean density of detections in Cav2.1 SSDs per synapse (number of detections/nm3) in synapses with (green) or without (grey) detectable release to a 10AP stimulus (Kolmogorov-Smirnov test, p = 0.122, nrelease = 121, nno release = 38). (F) Cumulative frequency distribution of the total number of detections in Cav2.1 SSDs per synapse in synapses with (green) or without (grey) detectable release to a 10AP stimulus (Kolmogorov-Smirnov test, p = 0.820, nrelease = 121, nno release = 38).

The probability of neurotransmitter release was previously shown to correlate well with active zone area (Holderith et al., 2012). We therefore investigated the relationship between Bassoon, one of the main presynaptic scaffolding proteins, and presynaptic function. Around 55% of synapses contained a single Bassoon SSD, with the rest of the synapses containing from 2 to 5 SSDs, much like our observations for Cav2.1 (Figure 6A). A similar positive correlation was found between ΔG/R and ΔR/R for synapses with single and multiple Bassoon SSDs (Figure 6B) and again, there was no difference in the relationship between Ca2+ and release between the two group of synapses when comparing the distributions of the ratio ΔG/ΔR (Figure 6C). We also found no correlation between Bassoon SSD properties and the levels of release or calcium influx, with the exception that in synapses with multiple Bassoon SSDs, the number of Bassoon detections was negatively correlated with the calcium dependence of release (ΔG/ΔR) (Supplementary Figure 6). However, when comparing synapses with detectable levels of release to those where release is not detectable, we observed differences in the nanoscale organisation of Bassoon. Releasing synapses tended to have fewer Bassoon SSDs (Figure 6D), as well as less Bassoon detections per synapse (Figure 6F). There was, however, no difference in the mean density of Bassoon distribution (Figure 6E).

FIGURE 6
www.frontiersin.org

Figure 6. Effect of Bassoon nanoscale organisation on neurotransmitter release. (A) Distribution of the number of Bassoon SSDs per synapse. (B) ΔG/R and ΔR/R responses from synapses post hoc labelled for Bassoon and identified as containing a single (blue) or multiple (black) Bassoon SSDs show a positive correlation (Spearman’s rank correlation, rsingle = 0.828, p < 0.001, nsingle = 49 synapses, rmulti = 0.794, p < 0.001, nmulti = 37 synapses). (C) Distribution of the ratio ΔG/ΔR for the two group of synapses. Distributions are not significantly different (Kolmogorov-Smirnov test, p = 0.575, nsingle = 49, nmultiple = 37). (D) Percentage of synapses with single or multiple Bassoon SSDs in synapses without detectable release (left, n = 16) or with release (right, n = 86). White lines indicate separation of multiple SSDs by number. (E) Cumulative frequency distribution of the mean density of detections in Bassoon SSDs per synapse (number of detections/nm3) in synapses with (green) or without (grey) detectable release to a 10AP stimulus (Kolmogorov-Smirnov test, p = 0.355, nrelease = 86, nno release = 16). (F) Cumulative frequency distribution of the total number of detections in Bassoon SSDs per synapse in synapses with (green) or without (grey) detectable release to a 10AP stimulus (Kolmogorov-Smirnov test, p = 0.041, nrelease = 86, nno release = 16). *p < 0.05.

Discussion

In this study we have described the use of a correlative approach combining functional imaging and multicolour 3D-dSTORM to investigate presynaptic structure and function at individual synapses. Using this method, we examined the nanoscale organisation of Bassoon and the VGCC Cav2.1. We found that these proteins are organised into SSDs, of which there can 1–6 per synapse (mean = 1.7 and 1.6, respectively), in agreement with previous studies using SRM (dSTORM Bassoon = 1.5, Tang et al., 2016; STED 1–4 Bassoon SSD, Hruska et al., 2018) and immunogold-EM (1–5 Cav2.1 SSD; Miki et al., 2017). In synapses with a single SSD, the density of the protein cluster was increased in comparison to the individual densities of the multiple SSDs, suggesting different strategies can be employed to arrange the same number of total proteins. At present it is not clear what determines the strategy at different synapses and further examination of the co-organisation of protein SSDs will be needed.

We found no correlation between the mean density or total number of Cav2.1 localisations and the levels of either calcium influx or neurotransmitter release. However, immunogold labelling in CA3-CA3 recurrent synapses has previously shown that the number of Cav2.1 channels correlates with AZ area, which in turn correlates with Pr (Holderith et al., 2012). In our study, the population of synapses we sampled from was more heterogeneous, since the promoter driving the expression of SypHy-RGECO could be expressed in any excitatory hippocampal neuron, which may account for any difference between the studies. Importantly, the calcium channel subtypes responsible for calcium influx have been shown to vary between synapses, adding another level of heterogeneity. Cav2.1, Cav2.2, and Cav2.3 have all been shown to support presynaptic calcium influx in glutamatergic terminals (Luebke et al., 1993; Wheeler et al., 1994; Gasparini et al., 2001), and the proportions of VGCC subtypes varied not only across hippocampal regions (Parajuli et al., 2012) but, more intriguingly, across synapses themselves, even when these synapses belong to the same neuron (Reid et al., 1997). In addition, it has been shown that postsynaptic cell identity can also influence presynaptic calcium influx and neurotransmitter release properties (Éltes et al., 2017). This heterogeneity in the complement of VGCC types, where some synapses are almost entirely reliant on one subtype, whilst others utilise a mixed population (Reid et al., 1997; Scheuber et al., 2004), needs to be better understood at the single synapse level. The complexity increases even further at the level of neurotransmitter release. In addition to VGCC number and composition, neurotransmitter release is known to be regulated by the spatial coupling of VGCCs and docked SVs (Eggermann et al., 2012; Rebola et al., 2019). The precise positioning of calcium channels next to docked vesicles is thought to be regulated by an ensemble of proteins such as RIM, ELKS and Munc13 (Wang et al., 2000; Kaeser et al., 2011; Dong et al., 2018; Quade et al., 2019) as well as the diffusive behaviour of VGCCs at the plasma membrane (Schneider et al., 2015; Heck et al., 2019). Any or all of these factors could contribute to the lack of correlation observed between Cav2.1 and functional parameters. It would be particularly interesting to use the same correlative approach to explore how the diverse complement of VGCC subtypes modulates synapse function, as well as the link between vesicle position and specific calcium channel subtypes.

Similarly, we found no correlation between the levels of Bassoon and either calcium influx or neurotransmitter release, in agreement with studies using other correlative methods (Holderith et al., 2020). We previously found that Bassoon clustering density was negatively correlated with neurotransmitter release and hypothesised this was due to a molecular crowding effect that limited the recruitment of other AZ proteins such as RIM and Cav2.1 (Glebov et al., 2017). Here, although we did not observe a negative correlation, we did find a significant increase in Bassoon detections and a larger proportion of multiple Bassoon SSDs in the small subset of synapses without detectable neurotransmitter release, consistent with increased Bassoon enrichment limiting presynaptic function. Furthermore, in synapses with multiple Bassoon SSDs in which release could be quantified, both the number and volume of Bassoon detections negatively correlated with the ratio ΔG/ΔR, a measure of the calcium dependence of release, suggesting a complex interplay between the spatial relationships of the AZ matrix and VGCCs and neurotransmitter release, which will require further investigation.

Using the number of SypHy detections obtained with dSTORM as a proportional measure of the total vesicle pool, we also saw no correlation with release or calcium influx, in agreement with EM studies showing the size of the total vesicle pool does not correlate with Pr, and that bouton volume shows only a weak relationship (Branco et al., 2010; Holderith et al., 2012).

Whilst the proteins we examined in this study showed few significant correlations with presynaptic functional measures, we have provided proof of principle of a correlative approach that can be used to examine functional properties and nanoscale protein organisation at individual presynaptic boutons. Nevertheless, there are some limitations to this technique which must be considered. First, as with all types of image analysis, detection thresholds must be applied. Here, using Voronoi segmentation for dSTORM imaging, we excluded low density detections that may arise from non-specific background or from isolated molecules that are not part of a cluster, perhaps because they are diffusing. By their nature these detections would not significantly alter the total number of counts and they are unlikely to substantially contribute to synapse function since they represent a minority of detections when compared to the larger clusters of protein. For SypHy, the highest density localisations arise from the total pool of SVs, which can be localised away from the AZ (Marra et al., 2012), hence this segmentation method may exclude localisations from surface stranded protein or individual SVs docked at the AZ. As the SSDs of Bassoon and Cav2.1 are located at the AZ, in some synapses they did not overlap with the segmented cluster of SypHy and had to be manually selected as belonging to the same bouton by visually inspecting individual SypHy localisations. If there was any ambiguity as to which SSD was associated with a SypHy cluster the synapse was excluded (Supplementary Figure 3). Second, despite applying manual correction, we still found a significant number of SypHy-RGECO positive synapses without any Cav2.1 or Bassoon SSD, for which there are several possible explanations. The protein may not have been present at the bouton, which is more likely in the case of Cav2.1 as it is known that different synapses can preferentially express other Cav2 subtypes (Reid et al., 1997). This is unlikely the case for Bassoon as it is thought to be ubiquitously expressed at synapses, although boutons lacking Bassoon maintain some enrichment of VGCCs and capacity to release (Davydova et al., 2014). Another explanation could be that the protein was expressed, but at a level below our detection threshold. Whilst we used standard immunolabelling techniques with commercially conjugated antibodies that offer high efficiency labelling, this may also vary between proteins and samples. Finally, the z-range in which we can localise fluorophores is limited to ∼1 um and depending on the size and orientation of the bouton in relation to the focal plane, the SSD of Bassoon or Cav2.1 may lie outside of this range. This limitation could possibly mean that some synapses were incorrectly classified as containing a single rather than multiple SSDs. This could occur if the synapse was sufficiently large to extend beyond the 1um range or the AZ SSD was detected only at the edge of the focal range and extended beyond it, which is only likely to be the case in a small proportion of synapses.

For functional imaging we used SypHy-RGECO, due to its ability to measure two aspects of presynaptic function simultaneously. Again, thresholding must be applied for both SypHy and RGECO channels, and some synapses display a calcium response but not a measurable release. We interpret these as having low release probability, as there may be a response that cannot be distinguished from baseline noise. There is a possibility that these are silent presynaptic terminals (Crawford and Mennerick, 2012), but to determine this would require using a stronger stimulus or more sensitive indicators to improve response detection, such as iGlusNFR for neurotransmitter release (Marvin et al., 2013, 2018, 2019) and GCaMP8 (Zhang et al., 2020) or jRCaMP1a/jRGECO1a (Dana et al., 2016) for calcium influx. Our correlative approach would work equally well with any of these reporters. In addition, whilst we used a 10AP stimulation protocol, which provides a robust response within the linear range for both reporters (Jackson and Burrone, 2016), different stimulation protocols could be employed to explore synaptic function in more depth than presented here.

The correlative approach would also be applicable to link function with the organisation of other synaptic proteins of interest and could be adapted to use different SRM methods. The main drawback of dSTORM imaging is the limited number of proteins that can be imaged at once. Here we have used two-colour imaging, with one channel occupied by GFP labelling for SypHy-RGECO, which is necessary to precisely relocate the same boutons. However, SD-dSTORM could allow the addition of a third fluorophore to examine two proteins of interest alongside GFP. In the dense protein environment of the presynapse, where proteins of interest are in very close proximity, care must be taken to ensure that single molecule localisation without crosstalk is achieved. Alternatively, Exchange-PAINT, SIM or the more recently developed MINFLUX could be used to improve the multiplexing possibilities of this technique, despite their own limitations (Jungmann et al., 2014; Gwosch et al., 2020; Butkevich et al., 2021). Live-SRM or single particle tracking could also add further important information about the dynamics of the nanoscale organisation of synaptic proteins and its role in regulating the calcium-dependence of release (Heine and Holcman, 2020).

Overall, we have demonstrated a flexible method that combines two imaging modalities at the same synapse to investigate functional properties and nanoscale organisation. Large numbers of synapses can be examined using this approach, enabling the wide range of synaptic heterogeneity to be explored, which will be necessary to fully understand the synaptic structure-function link.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

Ethical review and approval was not required for the animal study because only schedule 1 procedures performed by a competent individual were used in these studies, which are exempt under the Animals (Scientific Procedures) Act 1986.

Author Contributions

RJ, BC, and JB designed experiments, edited, and approved the manuscript. RJ performed and analysed functional experiments. BC performed and analysed dSTORM experiments. RJ and BC drafted the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was funded in whole, or in part, by the Wellcome Trust (Grant No. 215508/Z/19/Z to JB). For the purpose of open access, the authors have applied a CC BY public copyright license to any Author-Accepted Manuscript version arising from this submission. This work was additionally funded by an IBIN grant to RJ and BC and this project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 894933.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank Nicolas Bourg and Abbelight for the SAFeRedSTORM module and software, Olympus for the IX3 microscope, Florian Levet and Jean-Baptiste Sibarita (IINS, Bordeaux University) for PoCA software and Marcio Guiomar de Oliveira for preparation of dissociated hippocampal neurons.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnsyn.2022.830583/full#supplementary-material

References

Branco, T., Marra, V., and Staras, K. (2010). Examining size-strength relationships at hippocampal synapses using an ultrastructural measurement of synaptic release probability. J. Struct. Biol. 172, 203–210. doi: 10.1016/j.jsb.2009.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Butkevich, A. N., Weber, M., Cereceda Delgado, A. R., Ostersehlt, L. M., D’Este, E., and Hell, S. W. (2021). Photoactivatable fluorescent dyes with hydrophilic caging groups and their use in multicolor nanoscopy. J. Am. Chem. Soc. 143, 18388–18393. doi: 10.1021/jacs.1c09999

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, J.-H., Yu, N.-K., Baek, G.-C., Bakes, J., Seo, D., Nam, H. J., et al. (2014). Optimization of AAV expression cassettes to improve packaging capacity and transgene expression in neurons. Mol. Brain 7:17. doi: 10.1186/1756-6606-7-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Compans, B., Camus, C., Kallergi, E., Sposini, S., Martineau, M., Butler, C., et al. (2021). NMDAR-dependent long-term depression is associated with increased short term plasticity through autophagy mediated loss of PSD-95. Nat. Commun. 12:2849. doi: 10.1038/s41467-021-23133-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Crawford, D. C., and Mennerick, S. (2012). Presynaptically silent synapses: dormancy and awakening of presynaptic vesicle release. Neuroscientist 18, 216–223. doi: 10.1177/1073858411418525

PubMed Abstract | CrossRef Full Text | Google Scholar

Dana, H., Mohar, B., Sun, Y., Narayan, S., Gordus, A., Hasseman, J. P., et al. (2016). Sensitive red protein calcium indicators for imaging neural activity. eLife 5:e12727. doi: 10.7554/eLife.12727

PubMed Abstract | CrossRef Full Text | Google Scholar

Dani, A., Huang, B., Bergan, J., Dulac, C., and Zhuang, X. (2010). Superresolution imaging of chemical synapses in the brain. Neuron 68, 843–856. doi: 10.1016/j.neuron.2010.11.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Davydova, D., Marini, C., King, C., Klueva, J., Bischof, F., Romorini, S., et al. (2014). Bassoon specifically controls presynaptic P/Q-type Ca2+ channels via RIM-binding protein. Neuron 82, 181–194. doi: 10.1016/j.neuron.2014.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

del Castillo, J., and Katz, B. (1954). Quantal components of the end-plate potential. J. Physiol. 124, 560–573. doi: 10.1113/jphysiol.1954.sp005129

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodge, F. A., and Rahamimoff, R. (1967). Co-operative action of calcium ions in transmitter release at the neuromuscular junction. J. Physiol. 193, 419–432. doi: 10.1113/jphysiol.1967.sp008367

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, W., Radulovic, T., Goral, R. O., Thomas, C., Suarez Montesinos, M., Guerrero-Given, D., et al. (2018). CAST/ELKS proteins control voltage-gated Ca 2+ channel density and synaptic release probability at a mammalian central synapse. Cell Rep. 24, 284.e6–293.e6. doi: 10.1016/j.celrep.2018.06.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Dreosti, E., Odermatt, B., Dorostkar, M. M., and Lagnado, L. (2009). A genetically encoded reporter of synaptic activity in vivo. Nat. Methods 6, 883–889. doi: 10.1038/nmeth.1399

PubMed Abstract | CrossRef Full Text | Google Scholar

Eggermann, E., Bucurenciu, I., Goswami, S. P., and Jonas, P. (2012). Nanodomain coupling between Ca 2+ channels and sensors of exocytosis at fast mammalian synapses. Nat. Rev. Neurosci. 13, 7–21. doi: 10.1038/nrn3125

PubMed Abstract | CrossRef Full Text | Google Scholar

Éltes, T., Kirizs, T., Nusser, Z., and Holderith, N. (2017). Target cell type-dependent differences in Ca2+ channel function underlie distinct release probabilities at hippocampal glutamatergic terminals. J. Neurosci. 37, 1910–1924. doi: 10.1523/JNEUROSCI.2024-16.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Ermolyuk, Y. S., Alder, F. G., Henneberger, C., Rusakov, D. A., Kullmann, D. M., and Volynski, K. E. (2012). Independent regulation of basal neurotransmitter release efficacy by variable Ca2+ influx and bouton size at small central synapses. PLoS Biol. 10:e1001396. doi: 10.1371/journal.pbio.1001396

PubMed Abstract | CrossRef Full Text | Google Scholar

Früh, S. M., Matti, U., Spycher, P. R., Rubini, M., Lickert, S., Schlichthaerle, T., et al. (2021). Site-specifically-labeled antibodies for super-resolution microscopy reveal In Situ linkage errors. ACS Nano 15, 12161–12170. doi: 10.1021/acsnano.1c03677

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukata, Y., Dimitrov, A., Boncompain, G., Vielemeyer, O., Perez, F., and Fukata, M. (2013). Local palmitoylation cycles define activity-regulated postsynaptic subdomains. J. Cell Biol. 202, 145–161. doi: 10.1083/jcb.201302071

PubMed Abstract | CrossRef Full Text | Google Scholar

Gasparini, S., Kasyanov, A. M., Pietrobon, D., Voronin, L. L., and Cherubini, E. (2001). Presynaptic R-Type calcium channels contribute to fast excitatory synaptic transmission in the rat hippocampus. J. Neurosci. 21, 8715–8721. doi: 10.1523/JNEUROSCI.21-22-08715.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Glebov, O. O., Jackson, R. E., Winterflood, C. M., Owen, D. M., Barker, E. A., Doherty, P., et al. (2017). Nanoscale structural plasticity of the active zone matrix modulates presynaptic function. Cell Rep. 18, 2715–2728. doi: 10.1016/j.celrep.2017.02.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Granseth, B., Odermatt, B., Royle, S. J. J., and Lagnado, L. (2006). Clathrin-mediated endocytosis is the dominant mechanism of vesicle retrieval at hippocampal synapses. Neuron 51, 773–786. doi: 10.1016/j.neuron.2006.08.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Gwosch, K. C., Pape, J. K., Balzarotti, F., Hoess, P., Ellenberg, J., Ries, J., et al. (2020). MINFLUX nanoscopy delivers 3D multicolor nanometer resolution in cells. Nat. Methods 17, 217–224. doi: 10.1038/s41592-019-0688-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Haas, K. T., Compans, B., Letellier, M., Bartol, T. M., Grillo-Bosch, D., Sejnowski, T. J., et al. (2018). Pre-post synaptic alignment through neuroligin-1 tunes synaptic transmission efficiency. eLife 7:e31755. doi: 10.7554/eLife.31755

PubMed Abstract | CrossRef Full Text | Google Scholar

Harris, K. M., and Weinberg, R. J. (2012). Ultrastructure of synapses in the mammalian brain. Cold Spring Harb. Perspect. Biol. 4:a005587. doi: 10.1101/cshperspect.a005587

PubMed Abstract | CrossRef Full Text | Google Scholar

Heck, J., Parutto, P., Ciuraszkiewicz, A., Bikbaev, A., Freund, R., Mitlöhner, J., et al. (2019). Transient confinement of CaV2.1 Ca2+-Channel splice variants shapes synaptic short-term plasticity. Neuron 103, 66.e12–79.e12. doi: 10.1016/j.neuron.2019.04.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Heine, M., and Holcman, D. (2020). Asymmetry between pre- and postsynaptic transient nanodomains shapes neuronal communication. Trends Neurosci. 43, 182–196. doi: 10.1016/j.tins.2020.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Holderith, N., Heredi, J., Kis, V., and Nusser, Z. (2020). A high-resolution method for quantitative molecular analysis of functionally characterized individual synapses. Cell Rep. 32:107968. doi: 10.1016/j.celrep.2020.107968

PubMed Abstract | CrossRef Full Text | Google Scholar

Holderith, N., Lorincz, A., Katona, G., Rózsa, B., Kulik, A., Watanabe, M., et al. (2012). Release probability of hippocampal glutamatergic terminals scales with the size of the active zone. Nat. Neurosci. 15, 988–997. doi: 10.1038/nn.3137

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoppa, M. B., Gouzer, G., Armbruster, M., and Ryan, T. A. (2014). Control and plasticity of the presynaptic action potential waveform at small CNS nerve terminals. Neuron 84, 778–789. doi: 10.1016/j.neuron.2014.09.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Hruska, M., Henderson, N., Le Marchand, S. J., Jafri, H., and Dalva, M. B. (2018). Synaptic nanomodules underlie the organization and plasticity of spine synapses. Nat. Neurosci. 21, 671–682. doi: 10.1038/s41593-018-0138-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, R. E., and Burrone, J. (2016). Visualizing presynaptic calcium dynamics and vesicle fusion with a single genetically encoded reporter at individual synapses. Front. Synaptic Neurosci. 8:21. doi: 10.3389/fnsyn.2016.00021

PubMed Abstract | CrossRef Full Text | Google Scholar

Jungmann, R., Avendaño, M. S., Woehrstein, J. B., Dai, M., Shih, W. M., and Yin, P. (2014). Multiplexed 3D cellular super-resolution imaging with DNA-PAINT and Exchange-PAINT. Nat. Methods 11, 313–318. doi: 10.1038/nmeth.2835

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaeser, P. S., Deng, L., Wang, Y., Dulubova, I., Liu, X., Rizo, J., et al. (2011). RIM proteins tether Ca2+ Channels to presynaptic active zones via a direct PDZ-Domain interaction. Cell 144, 282–295. doi: 10.1016/j.cell.2010.12.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Levet, F., Hosy, E., Kechkar, A., Butler, C., Beghin, A., Choquet, D., et al. (2015). SR-Tesseler: a method to segment and quantify localization-based super-resolution microscopy data. Nat. Methods 12, 1065–1071. doi: 10.1038/nmeth.3579

PubMed Abstract | CrossRef Full Text | Google Scholar

Levet, F., Julien, G., Galland, R., Butler, C., Beghin, A., Chazeau, A., et al. (2019). A tessellation-based colocalization analysis approach for single-molecule localization microscopy. Nat. Commun. 10:2379. doi: 10.1038/s41467-019-10007-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., and Tsien, R. W. (2012). pHTomato, a red, genetically encoded indicator that enables multiplex interrogation of synaptic activity. Nat. Neurosci. 15, 1047–1053. doi: 10.1038/nn.3126

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipscombe, D., Allen, S. E., and Toro, C. P. (2013). Control of neuronal voltage-gated calcium ion channels from RNA to protein. Trends Neurosci. 36, 598–609. doi: 10.1016/j.tins.2013.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, A., Huang, X., He, W., Xue, F., Yang, Y., Liu, J., et al. (2021). pHmScarlet is a pH-sensitive red fluorescent protein to monitor exocytosis docking and fusion steps. Nat. Commun. 12:1413. doi: 10.1038/s41467-021-21666-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Luebke, I., Dunlap, K., and Turner, T. J. (1993). Multiple calcium channel types control glutamatergic synaptic transmission in the hippocampus. Neuron 11, 895–902. doi: 10.1016/0896-6273(93)90119-c

CrossRef Full Text | Google Scholar

MacGillavry, H. D., Song, Y., Raghavachari, S., and Blanpied, T. A. (2013). Nanoscale scaffolding domains within the postsynaptic density concentrate synaptic ampa receptors. Neuron 78, 615–622. doi: 10.1016/j.neuron.2013.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Marra, V., Burden, J. J., Thorpe, J. R., Smith, I. T., Smith, S. L., Häusser, M., et al. (2012). A preferentially segregated recycling vesicle pool of limited size supports neurotransmission in native central synapses. Neuron 76, 579–589. doi: 10.1016/j.neuron.2012.08.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Marvin, J. S., Borghuis, B. G., Tian, L., Cichon, J., Harnett, M. T., Akerboom, J., et al. (2013). An optimized fluorescent probe for visualizing glutamate neurotransmission. Nat. Methods 10, 162–170. doi: 10.1038/nmeth.2333

PubMed Abstract | CrossRef Full Text | Google Scholar

Marvin, J. S., Scholl, B., Wilson, D. E., Podgorski, K., Kazemipour, A., Müller, J. A., et al. (2018). Stability, affinity, and chromatic variants of the glutamate sensor iGluSnFR. Nat. Methods 15, 936–939. doi: 10.1038/s41592-018-0171-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Marvin, J. S., Shimoda, Y., Magloire, V., Leite, M., Kawashima, T., Jensen, T. P., et al. (2019). A genetically encoded fluorescent sensor for in vivo imaging of GABA. Nat. Methods 16, 763–770. doi: 10.1038/s41592-019-0471-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Maschi, D., and Klyachko, V. A. (2017). Spatiotemporal regulation of synaptic vesicle fusion sites in central synapses. Neuron 94, 65.e3–73.e3. doi: 10.1016/j.neuron.2017.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Miesenböck, G., de Angelis, D. A., and Rothman, J. E. (1998). Visualizing secretion and synaptic transmission with pH-sensitive green fluorescent proteins. Nature 394, 192–195. doi: 10.1038/28190

PubMed Abstract | CrossRef Full Text | Google Scholar

Miki, T., Kaufmann, W. A., Malagon, G., Gomez, L., Tabuchi, K., Watanabe, M., et al. (2017). Numbers of presynaptic Ca2+ channel clusters match those of functionally defined vesicular docking sites in single central synapses. Proc. Natl. Acad. Sci. U.S.A. 114, E5246–E5255. doi: 10.1073/pnas.1704470114

PubMed Abstract | CrossRef Full Text | Google Scholar

Murthy, V. N., Schikorski, T., Stevens, C. F., and Zhu, Y. (2001). Inactivity produces increases in neurotransmitter release and synapse size. Neuron 32, 673–682. doi: 10.1016/S0896-6273(01)00500-1

CrossRef Full Text | Google Scholar

Murthy, V. N., Sejnowski, T. J., and Stevens, C. F. (1997). Heterogeneous release properties of visualized individual hippocampal synapses. Neuron 18, 599–612. doi: 10.1016/S0896-6273(00)80301-3

CrossRef Full Text | Google Scholar

Nair, D., Hosy, E., Petersen, J. D., Constals, A., Giannone, G., Choquet, D., et al. (2013). Super-resolution imaging reveals that AMPA receptors inside synapses are dynamically organized in nanodomains regulated by PSD95. J. Neurosci. 33, 13204–13224. doi: 10.1523/JNEUROSCI.2381-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakai, J., Ohkura, M., and Imoto, K. (2001). A high signal-to-noise Ca2+ probe composed of a single green fluorescent protein. Nat. Biotechnol. 19, 137–141. doi: 10.1038/84397

PubMed Abstract | CrossRef Full Text | Google Scholar

Novak, P., Gorelik, J., Vivekananda, U., Shevchuk, A. I., Ermolyuk, Y. S., Bailey, R. J., et al. (2013). Nanoscale-targeted patch-clamp recordings of functional presynaptic ion channels. Neuron 79, 1067–1077. doi: 10.1016/j.neuron.2013.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Parajuli, L. K., Nakajima, C., Kulik, A., Matsui, K., Schneider, T., Shigemoto, R., et al. (2012). Quantitative regional and ultra structural localization of the Ca v2.3 subunit of R-type calcium channel in mouse brain. J. Neurosci. 32, 13555–13567. doi: 10.1523/JNEUROSCI.1142-12.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Platonova, E., Winterflood, C. M., and Ewers, H. (2015). A simple method for GFP- and RFP-based dual color single-molecule localization microscopy. ACS Chem. Biol. 10, 1411–1416. doi: 10.1021/acschembio.5b00046

PubMed Abstract | CrossRef Full Text | Google Scholar

Quade, B., Camacho, M., Zhao, X., Orlando, M., Trimbuch, T., Xu, J., et al. (2019). Membrane bridging by Munc13-1 is crucial for neurotransmitter release. eLife 8:e42806. doi: 10.7554/eLife.42806

PubMed Abstract | CrossRef Full Text | Google Scholar

Rebola, N., Reva, M., Kirizs, T., Szoboszlay, M., Lőrincz, A., Moneron, G., et al. (2019). Distinct nanoscale calcium channel and synaptic vesicle topographies contribute to the diversity of synaptic function. Neuron 104, 693.e9–710.e9. doi: 10.1016/j.neuron.2019.08.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Reid, C. A., Clements, J. D., and Bekkers, J. M. (1997). Nonuniform distribution of Ca 2 channel subtypes on presynaptic terminals of excitatory synapses in hippocampal cultures. J. Neurosci. 17, 2738–2745. doi: 10.1523/JNEUROSCI.17-08-02738.1997

PubMed Abstract | CrossRef Full Text | Google Scholar

Rose, T., Schoenenberger, P., Jezek, K., and Oertner, T. G. (2013). Developmental refinement of vesicle cycling at schaffer collateral synapses. Neuron 77, 1109–1121. doi: 10.1016/j.neuron.2013.01.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabater, V. G., Rigby, M., and Burrone, J. (2021). Voltage-gated potassium channels ensure action potential shape fidelity in distal axons. J. Neurosci. 41, 5372–5385. doi: 10.1523/JNEUROSCI.2765-20.2021

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakamoto, H., Ariyoshi, T., Kimpara, N., Sugao, K., Taiko, I., Takikawa, K., et al. (2018). Synaptic weight set by Munc13-1 supramolecular assemblies. Nat. Neurosci. 21, 41–55. doi: 10.1038/s41593-017-0041-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Sankaranarayanan, S., de Angelis, D., Rothman, J. E., and Ryan, T. A. (2000). The use of pHluorins for optical measurements of presynaptic activity. Biophys. J. 79, 2199–2208. doi: 10.1016/S0006-3495(00)76468-X

CrossRef Full Text | Google Scholar

Scheuber, A., Miles, R., and Poncer, J. C. (2004). Presynaptic Cav2.1 and Cav2.2 differentially influence release dynamics at hippocampal excitatory synapses. J. Neurosci. 24, 10402–10409. doi: 10.1523/JNEUROSCI.1664-04.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Schikorski, T., and Stevens, C. F. (1997). Quantitative ultrastructural analysis of hippocampal excitatory synapses. J. Neurosci. 17, 5858–5867. doi: 10.1523/JNEUROSCI.17-15-05858.1997

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, R., Hosy, E., Kohl, J., Klueva, J., Choquet, D., Thomas, U., et al. (2015). Mobility of calcium channels in the presynaptic membrane. Neuron 86, 672–679. doi: 10.1016/j.neuron.2015.03.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Sigal, Y. M., Zhou, R., and Zhuang, X. (2018). Visualizing and discovering cellular structures with super-resolution microscopy. Science 361, 880–887. doi: 10.1126/science.aau1044

PubMed Abstract | CrossRef Full Text | Google Scholar

Südhof, T. C. (2012). The presynaptic active zone. Neuron 75, 11–25. doi: 10.1016/j.neuron.2012.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, A. H., Chen, H., Li, T. P., Metzbower, S. R., MacGillavry, H. D., and Blanpied, T. A. (2016). A trans-synaptic nanocolumn aligns neurotransmitter release to receptors. Nature 536, 210–214. doi: 10.1038/nature19058

PubMed Abstract | CrossRef Full Text | Google Scholar

Testa, I., Wurm, C. A., Medda, R., Rothermel, E., von Middendorf, C., Fölling, J., et al. (2010). Multicolor fluorescence nanoscopy in fixed and living cells by exciting conventional fluorophores with a single wavelength. Biophys. J. 99, 2686–2694. doi: 10.1016/j.bpj.2010.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Voglmaier, S. M., Kam, K., Yang, H., Fortin, D. L., Hua, Z., Nicoll, R. A., et al. (2006). Distinct endocytic pathways control the rate and extent of synaptic vesicle protein recycling. Neuron 51, 71–84. doi: 10.1016/j.neuron.2006.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, A. S., Burrone, J., and Meyer, M. P. (2013). Functional imaging in the zebrafish retinotectal system using RGECO. Front. Neural Circuits 7:34. doi: 10.3389/fncir.2013.00034

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Kim, C. K., and Ting, A. Y. (2019). Molecular tools for imaging and recording neuronal activity. Nat. Chem. Biol. 15, 101–110. doi: 10.1038/s41589-018-0207-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Sugita, S., and Südhof, T. C. (2000). The RIM/NIM family of neuronal C2 domain proteins: interactions with Rab3 and a new class of Src homology 3 domain proteins. J. Biol. Chem. 275, 20033–20044. doi: 10.1074/jbc.M909008199

PubMed Abstract | CrossRef Full Text | Google Scholar

Wheeler, D. B., Randall, A., and Tsien, R. W. (1994). Roles of N-type and Q-type Ca2+ channels in supporting hippocampal synaptic transmission. Science 264, 107–111. doi: 10.1126/science.7832825

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu-Friedman, M. A., Harris, K. M., and Regehr, W. G. (2001). Three-dimensional comparison of ultrastructural characteristics at depressing and facilitating synapses onto cerebellar purkinje cells. J. Neurosci. 21, 6666–6672. doi: 10.1523/JNEUROSCI.21-17-06666.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Rózsa, M., Bushey, D., Zheng, J., Reep, D., Broussard, G. J., et al. (2020). jGCaMP8: a new suite of fast and sensitive calcium indicators. doi: 10.25378/janelia.13148243

CrossRef Full Text | Google Scholar

Keywords: synapse, super-resolution imaging, neurotransmitter release, calcium, active zone (AZ), correlative imaging

Citation: Jackson RE, Compans B and Burrone J (2022) Correlative Live-Cell and Super-Resolution Imaging to Link Presynaptic Molecular Organisation With Function. Front. Synaptic Neurosci. 14:830583. doi: 10.3389/fnsyn.2022.830583

Received: 07 December 2021; Accepted: 12 January 2022;
Published: 15 February 2022.

Edited by:

Lucia Tabares, Seville University, Spain

Reviewed by:

Matteo Fossati, Institute of Neuroscience, National Research Council (CNR), Italy
Vincenzo Marra, University of Leicester, United Kingdom

Copyright © 2022 Jackson, Compans and Burrone. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Rachel E. Jackson, cmFjaGVsLmphY2tzb25Aa2NsLmFjLnVr; Benjamin Compans, YmVuamFtaW4uY29tcGFuc0BrY2wuYWMudWs=; Juan Burrone, anVhbi5idXJyb25lQGtjbC5hYy51aw==

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.