Introduction

One of the most studied phenomena triggered by prolonged activity blockade is the increase in amplitudes of miniature excitatory postsynaptic currents (mEPSCs), first demonstrated in cultures of dissociated cortical neurons (Turrigiano et al., 1998) and spinal cord neurons (O’Brien et al., 1998). This compensatory response is now termed homeostatic synaptic plasticity (Turrigiano and Nelson, 2004; Pozo and Goda, 2010) and has been shown to be dysregulated in several neurodevelopment, psychiatric, and neurodegenerative disorders. One of the initial studies of homeostatic synaptic plasticity demonstrated that AMPAergic mEPSCs were uniformly increased after prolonged activity blockade (Turrigiano et al., 1998), and the authors named this process synaptic scaling. An accompanying increase in AMPA-type glutamate receptors (AMPARs) was observed in both of the early studies of homeostatic synaptic plasticity (O’Brien et al., 1998; Turrigiano et al., 1998), and has been confirmed many times (a non-exhaustive list includes (Ju et al., 2004; Thiagarajan et al., 2005; Shepherd et al., 2006; Stellwagen and Malenka, 2006; Hou et al., 2008; Gainey et al., 2009; Soden and Chen, 2010; Correa et al., 2012; Altimimi and Stellwagen, 2013; Letellier et al., 2014; Xu and Pozzo- Miller, 2017).

Homeostatic synaptic plasticity is becoming increasingly implicated in both pathological brain conditions, for example, epilepsy, neuropsychiatric disorders, Huntington’s, and alcohol use disorder (Trasande and Ramirez, 2007; Fernandes and Carvalho, 2016; Wang et al., 2017; Lovinger and Abrahao, 2018; Smith-Dijak et al., 2019; Lignani et al., 2020; Suzuki et al., 2021; Kavalali and Monteggia, 2023), and essential normal functions such as sleep (Tononi and Cirelli, 2014; Diering et al., 2017; Torrado Pacheco et al., 2021), so a molecular understanding of the process is an extremely important next step. The homeostatic synaptic plasticity field has identified several proteins required for synaptic scaling of mEPSC amplitudes, the majority of which are involved in the regulation of AMPAR levels (Shepherd et al., 2006; Seeburg and Sheng, 2008; Gao et al., 2010; Anggono et al., 2011; Beique et al., 2011; Diering et al., 2014; Gainey et al., 2015; Tan et al., 2015; Pastuzyn and Shepherd, 2017; Sanderson et al., 2018).

In our previous work studying the increase in miniature endplate currents (mEPC) following prolonged in vivo activity-block of synaptic transmission at the mouse neuromuscular junction (NMJ), we were surprised to find no evidence of changes in acetylcholine receptor (AChR) levels (Wang et al., 2005). This result led us to search for presynaptic molecules that might homeostatically regulate mEPC amplitude, perhaps via the presynaptic quantum. In previous studies in chromaffin cells, we identified the small GTPase Rab3A, a synaptic vesicle protein, as a regulator of synaptic vesicle fusion pore opening (Wang et al., 2008), so we examined whether deletion of Rab3A (Rab3A-/-) might prevent homeostatic upregulation of mEPC amplitude. The results were clear: in the Rab3A-/- mouse, the homeostatic increase in mEPC amplitude was strongly reduced, and was completely abolished in mice expressing a point mutation in Rab3A, the Earlybird mutant (Rab3AEbd/Ebd) (Wang et al., 2011).

The Rab3A-/- mouse has minimal phenotypic abnormalities, with evoked synaptic transmission and mEPSCs essentially normal in hippocampal slices (Geppert et al., 1994). At the Rab3A-/-NMJ, reductions in evoked transmission were detected, but only under conditions of reduced extracellular calcium (Coleman et al., 2007). The most dramatic effect of loss of Rab3A is the disruption of a presynaptic form of long-term potentiation (LTP) at the mossy fiber-CA3 synapse (Weisskopf et al., 1994; Castillo et al., 1997). To our knowledge, there is no evidence of Rab3A involvement in the expression or trafficking of postsynaptic receptors. Our results at the mammalian NMJ suggest that in addition to its importance in mossy fiber LTP, Rab3A may be required for the homeostatic plasticity of mEPSC amplitude via a presynaptic mechanism.

In the current work, we explored whether the in vivo findings at the NMJ might also apply to the more typically studied homeostatic synaptic plasticity of mEPSC amplitude in dissociated cortical neuron cultures. We report that our findings in vivo at the NMJ were almost exactly recapitulated in cultures of dissociated cortical neurons: 1. strong reduction of the homeostatic increase in mEPSC amplitude in the absence of Rab3A; 2. complete abolishment of the homeostatic increase in mEPSC amplitude in the presence of the Rab3A Earlybird mutant; and 3. increased mEPSC amplitude in the Rab3A Earlybird mutant prior to activity blockade.

However, in contrast to the unchanged AChR levels at the mammalian NMJ, there was a modest increase in levels of the GluA2-type AMPARs at cortical synapses after activity- blockade, which appeared to be disrupted in cultures prepared from Rab3A-/-mice.

Importantly, when compared within the same cultures, GluA2 receptor levels did not always parallel mEPSC amplitudes. We also determined that Rab3A must be present in neurons, but not astrocytes, for full expression of homeostatic plasticity. The NMJ and cortical culture data taken together strongly suggest neuronal Rab3A is important for both postsynaptic receptor upregulation and a second mechanism affecting mEPSC amplitude, likely the amount of transmitter released by a single vesicle.

Materials and Methods

Animals

Rab3A+/- heterozygous mice were bred and genotyped as previously described (Kapfhamer et al., 2002; Wang et al., 2008). Rab3AEbd/Ebdmice were identified in an EU- mutagenesis screen of C57BL/6J mice, and after a cross to C3H/HeJ, were backcrossed for 3 generations to C57BL/6J (Kapfhamer et al., 2002). Rab3A+/Ebd heterozygous mice were bred at Wright State University and genotyped in a two-step procedure: 1. a PCR reaction with RabF1 and Dcaps3R as primers; and 2. a digestion with enzyme Bsp1286I (New England Biolabs) that distinguishes the Earlybird mutant by its different base-pair products. Rab3A+/- mice were backcrossed with Rab3A+/+ mice from the Earlybird heterozygous colony for 11 generations in an attempt to establish a single wild type strain, but differences in mEPSC amplitude and adrenal chromaffin cell calcium currents persisted, likely due to genes that are close to the Rab3A site, resulting in two wild type strains: 1. Rab3A+/+ from the Rab3A+/- colony, and 2. Rab3A+/+ from the Rab3A+/Ebd colony.

Primary Culture of Mouse Cortical Neurons

Primary dissociated cultures of mixed neuronal and astrocyte populations were prepared as previously described (Hanes et al., 2020). Briefly, postnatal day 0-2 (P0-P2) Rab3A+/+, Rab3A-/- or Rab3AEbd/Ebd neonates were euthanized by rapid decapitation, as approved by the Wright State University Institutional Animal Care and Use Committee, and brains were quickly removed. Each culture was prepared from the cortices harvested from two animals; neonates were not sexed. Cortices were collected in chilled Neurobasal-A media (Gibco) with osmolarity adjusted to 270 mOsm and supplemented with 40 U/ml DNAse I (ThermoFisher Scientific). The tissues were digested with papain (Worthington Biochemical) at 20 U/ml at 37°C for 20 minutes followed by trituration with a sterile, fire- polished Pasteur pipette, then filtered through a 100 μm cell strainer, and centrifuged at 1100 rpm for 2 minutes. After discarding the supernatant, the pellet was resuspended in room temperature Neurobasal-A media (270 mOsm), supplemented with 5% fetal bovine serum for astrocyte growth, and 2% B-27 supplement to promote neuronal growth (Gibco), L-glutamine, and gentamicin (ThermoFisher Scientific). Neurons were counted and plated at 0.15 * 106 cells/coverslip onto 12 mm coverslips pre-coated with poly-L-lysine (BioCoat, Corning). The culture media for the first day (0 DIV) was the same as the above Neurobasal-A media supplemented with FBS, B-27, L-glutamine, and gentamicin, and was switched after 24 hours (1 DIV) to media consisting of Neurobasal-A (270 mOsm), 2% B-27, and L-glutamine without FBS to avoid its toxic effects on neuronal viability and health (Stellwagen and Malenka, 2006). Half of the media was changed twice weekly and experiments were performed at 13-14 DIV. Two days prior to experiments, tetrodotoxin (TTX) (500 nM; Tocris), a potent Na+ channel blocker, was added to some cultures to chronically silence all network activity and induce homeostatic synaptic plasticity mechanisms, while untreated sister cultures served as controls. Cultures prepared from mutant mice were compared with cultures from wild-type mice from their respective colonies. Note that the cultures comprising the Rab3A+/+ data here are a subset of the data previously published in Hanes et al., 2020, and therefore the plots in Figure 1 are not identical to those in the previously published work. This smaller data set was restricted to the time period over which cultures were prepared from Rab3A-/- mice.

Block of activity with 48 hours TTX treatment increased mEPSC amplitudes in cortical cultures from wild type mice but not Rab3A-/- mice. A) 10 second current traces recorded at -60 mV in pyramidal cortical neurons from an untreated (“CON”) and TTX-treated (“TTX”) neuron in cultures prepared from Rab3A+/+ mice from the Rab3A+/- colony. B) Average traces for the recordings shown in (A). C) Box and whisker plots for average mEPSC amplitudes from 30 CON and 23 TTX cells (Rab3A+/+, Rab3A+/- colony; mean mEPSC amplitudes, CON, 13.9 ± 0.7 pA; TTX, 18.2 ± 0.9 pA, p = 4.58 * 10-4, Kruskal-Wallis test). D) Current traces recorded in pyramidal cortical neurons from a CON and TTX-treated neuron in cultures prepared from Rab3A-/- mice. E) Average traces for the recordings shown in (D). F) Box and whisker plots for average mEPSC amplitudes from 25 CON and 26 TTX cells, Rab3A-/- data (Rab3A-/-, mean mEPSC amplitudes, CON, 13.6 ± 0.1; TTX, 14.3 ± 0.6, p = 0.318, Kruskal-Wallis test). G) Analysis of pooled mEPSC amplitudes for cultures from Rab3A+/+ mice of the Rab3A+/- colony (i) Cumulative distribution functions (CDFs) of quantiles of mEPSC amplitudes pooled from CON (black curve) and TTX- treated cultures (red curve) (Kolmogorov Smirnov (KS) test, test statistic (D) = 0.149, p = 8.53 * 10-13. (ii) Quantiles pooled from CON and TTX-treated cells were ranked from smallest to largest, plotted against each other, and fit with linear regression functions. Gray line, identity; red dashed line, linear regression fit. Inset, full range of the data. (iii) The same ranked mEPSC quantiles were paired and the ratio of CON/TTX computed for each pair, then plotted as a function of the CON mEPSC for the pair. Dashed line at ratio = 1.0 (no plasticity). H) Analysis of pooled mEPSCs amplitudes for cultures from Rab3A-/- mice. (i) CDFs of quantiles pooled from CON (black curve) and TTX-treated cultures (red curve) (KS test, D = 0.070, p = 0.042). (ii) Quantiles pooled from CON and TTX-treated cells were ranked from smallest to largest, plotted against each other, and fit with a linear regression function. Gray lines, identity; red dashed lines, linear regression fits. Inset: full range of the data. (iii) The same ranked mEPSC quantiles were paired and the ratio of CON/TTX computed for each pair, then plotted as a function of the CON mEPSC for the pair. Dashed line at ratio = 1 (no plasticity).

Preparation of Astrocyte Feeder Layers

Astrocyte feeder layers were prepared from the cortices of P0-P2 Rab3A+/+ or Rab3A-/- mouse pups as described previously (Stellwagen and Malenka, 2006). Briefly, cortices were dissected and cells were dissociated as described above. Cell suspensions of mixed neuronal and astrocyte populations were plated onto glass coverslips pre-coated with poly-L-lysine in Dulbecco’s Modified Eagle Media (ThermoFisher Scientific) supplemented with 5% FBS (to promote astrocyte proliferation and to kill neurons), L- glutamine, and gentamicin, and maintained in an incubator at 37°C, 5% CO2; cultures were maintained in this manner for up to 1 month to generate purely astrocytic cultures (all neurons typically died off by 7 DIV). Culture media was replaced after 24 hours, and subsequent media changes were made twice weekly, replacing half of the culture media with fresh media. Feeder layers were not used for neuronal seeding until all native neurons were gone and astrocytes approached 100% confluency (visually inspected).

Plating of Neurons on Glial Feeder Layers

Cortical neurons were obtained as described above. The cell pellet obtained was resuspended in Neurobasal-A (osmolarity adjusted to 270 mOsm) containing B27 (2%, to promote neuronal growth), L-glutamine, and 5-fluorodeoxyuridine (FdU, a mitotic inhibitor; Sigma). Addition of FdU was used to prevent astrocyte proliferation and contamination of the feeder layer with new astrocytes, promoting only neuronal growth on the feeder layers (FdU-containing media was used for the maintenance of these cultures and all subsequent media changes). Astrocyte culture media was removed from the feeder layer cultures, and the neuronal cell suspension was plated onto the astrocyte feeder cultures. The culture strategy used to distinguish the relative roles of neuronal and astrocytic Rab3A is outlined in Figure 9. At 1 DIV, all of the culture media was removed and replaced with fresh Neurobasal-A media containing FdU described above, and half of the media was replaced twice per week for all subsequent media changes. Cultures were maintained in a 37°C, 5% CO2 incubator for 13-14 DIV.

Whole-Cell Voltage Clamp to Record mEPSCs

At 13-14 DIV, mEPSCs from TTX-treated and untreated sister cultures of Rab3A+/+ or Rab3A-/- neurons from the Rab3A+/- colony, or Rab3A+/+ or Rab3AEbd/Ebd neurons from the Rab3A+/Ebd colony, were recorded via whole-cell voltage clamp to assess the role of Rab3A in homeostatic synaptic plasticity. Recordings were taken from pyramidal neurons, which were identified visually by a prominent apical dendrite; images were taken of all cells recorded from. Cells were continuously perfused with a solution consisting of (in mM): NaCl (115), KCl (5), CaCl2 (2.5), MgCl2 (1.3), dextrose (23), sucrose (26), HEPES (4.2), pH = 7.2 (Stellwagen and Malenka, 2006). On the day of recording, the osmolarity of the media from the cultures was measured (normally 285 – 295 mOsm) and the perfusate osmolarity was adjusted to match the culture osmolarity, to protect against osmotic shock to the neurons. To isolate glutamatergic mEPSCs, TTX (500 nM) and picrotoxin (50 μM) were included in the perfusion solution to block action potentials and GABAergic currents, respectively. The NMDA receptor antagonist, APV, was not included in the perfusion solution because all mEPSCs were blocked by CNQX and picrotoxin, demonstrating no APV-sensitive mEPSCs were present (data not shown). Patch electrodes (3.5 – 5 MΩ) were filled with an internal solution containing (in mM): K-gluconate (128), NaCl (10), EGTA (1), CaCl2 (0.132), MgCl2 (2), HEPES (10), pH = 7.2.

Osmolarity was adjusted to 10 mOsm less than the perfusion solution osmolarity. Neurons were clamped at a voltage of -60 mV using an Axopatch 200B patch-clamp (Axon Instruments), recorded from for 2 – 5 minutes, and data were collected with Clampex 10.0/10.2 (Axon Instruments). The antagonist of Ca2+-permeable AMPA receptors (including GluA1 but not GluA2), N-naphthyl acetylspermine (NASPM, 20 μM; Tocris), was applied during recordings in a subset of experiments. Because NASPM is an open channel blocker, it was applied with a depolarizing high K+ solution (25 mM KCl, 95 mM NaCl). Baseline recordings were performed for 2 minutes in our standard perfusate, then were suspended while NASPM + High K+ solution was applied for 45 seconds, followed by a NASPM only solution for 5 minutes, after which recording was recommenced for 5 minutes (because we found in pilot experiments that frequency was reduced following NASPM application).

Data Analysis and Statistics of mEPSCs

Miniature excitatory postsynaptic currents were manually selected using Mini Analysis software (Synaptosoft) to identify mEPSCs. The program threshold was set at 3 pA but the smallest mEPSC selected was 4.01 pA. Records were filtered at 2 kHz using a low-pass Butterworth filter prior to selection. For computing mEPSC mean in each experimental condition, individual cell means were pooled across multiple cultures and compared with the non-parametric Kruskal-Wallis test, with n = the number of cells and the overall means presented as ± SEM. For cumulative probability distribution functions (CDFs) of mEPSC amplitude, 30 quantiles were computed for each cell and pooled across cultures (see (Hanes et al., 2020)), and a Kolmogorov-Smirnov test (KS test) was used to test for significant differences, with n = the number of mEPSC quantiles. The rank order plots were created by computing a matched number of quantiles from the two experimental conditions (usually control (CON) and TTX). We used an algorithm to identify the product closest to, but above, 24 quantiles for a data set. For example, if there were 13 CON cells and 12 TTX cells, we computed 24 quantiles for each of the 13 cells and 26 quantiles for each of the 12 cells, for an equal total number of 312 quantiles for both data sets. The quantiles were sorted from smallest to largest, TTX amplitudes plotted vs. CON amplitudes, and the relationship fit with a linear regression function with the intercept term allowed to vary. The ratio plots were created by taking the ratio of TTX/CON and plotting as a function of the CON amplitude quantiles.

Immunocytochemistry, microscopy, and data analysis

Primary cultures of mouse cortical neurons were grown for 13-14 DIV. Antibodies to GluA2 (mouse ab against N-terminal, EMD Millipore) were added directly to live cultures at 1:40 dilution, and incubated at 37 °C in a CO2 incubator for 45 minutes. Cultures were rinsed 3 times with PBS/5% donkey serum before being fixed with 4% paraformaldehyde. After 3 rinses in PBS/5% donkey serum, cultures were incubated in CY3-labeled donkey-anti-mouse secondary antibodies for 1 hour at room temperature, rinsed in PBS/5% donkey serum, permeabilized with 0.2% saponin, and incubated in chick anti-MAP2 (1:2500, AbCAM) and rabbit anti-VGLUT1 (1:4000, Synaptic Systems) for 1 hour at room temperature in PBS/5% donkey serum. After rinsing with PBS/5% donkey serum, coverslips were incubated with 488-anti chick and CY5-anti rabbit secondary antibodies for 1 hour at room temperature, rinsed, blotted to remove excess liquid, and inverted on a drop of Vectashield (Vector Labs). Coverslips were sealed with nail polish and stored at 4 °C for < 1 week before imaging. All secondary antibodies were from Jackson Immunoresearch and were used at 1:225 dilution.

Coverslips were viewed on a Fluoview FV1000 laser scanning confocal microscope with a 60x oil immersion, 1.35 NA objective. Once a pyramidal neuron was identified, Fluoview 2.1 software was used to zoom in on the primary dendrite (5X) and confocal sections were taken every 0.5 μm. Images were analyzed offline with ImagePro 6 (Cybernetics). The composite image was used to locate synaptic sites containing both VGLUT1 and GluA2 immunoreactivity in close apposition to each other and to the primary dendrite or a secondary branch. An area of interest (AOI) was manually drawn around the receptor cluster in the confocal section in which it was the brightest. The AOIs for a dendrite were saved in a single file; the AOI number and the confocal section it was associated with were noted for later retrieval. For quantification, AOIs were loaded, an individual AOI was called up on the appropriate section, and the count/measurement tool used to apply a threshold (400 or 450; identical for CON and TTX coverslips but different for different cultures); pixels within the cluster that were above the threshold were automatically outlined, and size, average intensity, and integral of the outlined region reported. For quantile sampling, only dendrites with a minimum of 6 synaptic sites were included. For creating CDFs, 30 quantiles were computed for each dendrite. For the rank order and ratio plots, we matched the total quantiles for CON and TTX GluA2 receptor characteristics as was done above for mEPSC amplitudes.

Results

We previously reported that mixed cultures of cortical neurons and astrocytes prepared from postnatal day 0-2 mouse pups responded to a block of action potential-mediated activity by a 48 hour TTX treatment with an increase in mEPSC amplitude (Hanes et al., 2020). Here, we asked the question whether cortical cultures prepared from mice lacking the small GTPase Rab3A, or, expressing a point mutation of Rab3A, Rab3A Earlybird, have an altered homeostatic plasticity response to a loss of network activity. To obtain Rab3A-/- and Rab3AEbd/Ebdhomozygotes, we established two mouse colonies of heterozygous breeders with cultures prepared from pups derived from a final breeding pair of homozygotes. Although we backcrossed Rab3A+/- with Rab3A+/Ebd for 11 generations, clear differences in mEPSC amplitudes in untreated cultures (see below) and in calcium current amplitudes in adrenal chromaffin cells (unpublished obs.) remained. Therefore, throughout this study we keep the two strains separate and there are two Rab3A+/+ or ‘wild type’ phenotypes.

Example current traces of spontaneously occurring mEPSCs recorded from pyramidal neurons in untreated (CON) 13-14 DIV cortical cultures and sister cultures treated with 500 nM TTX for 48 hours prepared from wild type animals in the Rab3A+/- colony are shown in Figure 1A. Average mEPSC waveforms from the same recordings are shown in Figure 1B. The mean mEPSC amplitudes for 30 control and 23 TTX-treated neurons are displayed in the box and whisker plot in Figure 1C; after activity blockade the average mEPSC amplitude increased from 13.9 ± 0.7 pA to 18.2 ± 0.9 pA (p = 4.58 * 10-4, Kruskal-Wallis test). Example current traces of mEPSCs and average mEPSC waveforms are shown in Figures 1D and E, respectively, for cortical cultures prepared from Rab3A-/- mice. In contrast to the behavior of neurons in cultures prepared from the wild type strain, the increase in mEPSC amplitudes in Rab3A-/- neurons after activity blockade was dramatically reduced, and the average mEPSC amplitude was not significantly increased, for 25 untreated cells and 26 TTX-treated cells (Figure 1F, 13.6 ± 0.1 vs. 14.3 ± 0.6, p = 0.318, Kruskal-Wallis test). To further examine homeostatic plasticity in the presence and absence of Rab3A, we computed 30 quantiles for the mEPSC amplitude distribution of each neuron (Hanes et al., 2020), pooled the quantiles, and plotted the data as cumulative distribution functions (CDFs) for CON and TTX. For mEPSC amplitudes from cultures prepared from Rab3A+/+ mice, the difference between CDFs was highly significant (Figure 1Gi, test statistic D = 0.172, p = 1.62 * 10-10, Kolmogorov Smirnov (KS) test). In contrast, the increase in mEPSC amplitude after activity blockade was dramatically reduced in cortical cultures from Rab3A-/- mice (Figure 1Hi). A KS test for the CDFs was significant, but the test statistic much smaller (0.070) and the p value much larger (0.042) than for the wild type strain.

It was originally proposed that loss of activity produced a uniform multiplicative increase in mEPSC amplitude across the entire distribution of mEPSCs (Turrigiano et al., 1998). This uniform scaling would be expected to preserve original differences in synaptic weights resulting from other forms of plasticity (Turrigiano, 1999; Turrigiano and Nelson, 2004). In addition, uniform scaling suggested a cell-wide mechanism (in fact, a culture-wide mechanism) to identically modify all synapses. However, we recently showed that in mouse cortical neurons, rat cortical neurons, and mouse hippocampal neurons, scaling is non-uniform, a phenomenon we called “divergent scaling;” the multiplicative factor is smallest (close to 1) for small mEPSC amplitudes and increases to ∼1.4 for larger mEPSC amplitudes (Hanes et al., 2020). In subsequent reviewing of the current literature (Koesters et al., 2022), we concluded that the majority of previous studies of homeostatic plasticity of mEPSC amplitude following activity blockade also show divergent scaling, most notably the study that originally defined synaptic scaling (Turrigiano et al., 1998). Divergent scaling is not obvious on the standard plot used to demonstrate synaptic scaling, the rank-order plot (Turrigiano et al., 1998), but is only apparent once the ratio of TTX mEPSC amplitude/CON mEPSC amplitude is computed and plotted as a function of CON mEPSC amplitude (Hanes et al., 2020). To compare between uniform and divergent scaling in the presence and absence of Rab3A, we created both rank order plots and ratio plots for Rab3A+/+ data. The rank order plot for Rab3A+/+ neurons (Figure 1Gii) shows that treatment with TTX appears to cause a uniform multiplication of control mEPSC amplitudes based on the linear fit to the data, which had a slope of 1.43 (R2 = 0.982). However, in the ratio plot for mEPSCs in Rab3A+/+ cultures, the scaling was clearly divergent with the smallest mEPSC amplitudes having the smallest ratios (at 1) and the ratio increasing with mEPSC amplitude to a maximum of 1.51 (Figure 1Giii). In the absence of Rab3A, the rank ordered mEPSC data follow the line of identity and the slope of the linear fit was 0.97 (R2 = 0.985), suggesting complete absence of homeostatic plasticity (Figure 1Hii). However, the more sensitive ratio plot reveals a residual divergent scaling to a peak ratio of 1.13 (Figure 1Hiii). We previously observed a residual homeostatic effect of activity blockade on mEPCs at the mouse NMJ in the absence of Rab3A (Wang et al., 2011) suggesting that there may be compensatory or redundant mechanisms after deletion of Rab3A.

We next examined whether expression of a single point mutant of Rab3A (Rab3AEbd/Ebd) would abolish synaptic scaling in the mouse cortical cultures as we previously showed for mEPCs at the mouse NMJ (Wang et al., 2011). Example current traces of mEPSCs, average mEPSC waveforms, and box and whisker plots are shown in Figures 2A and B, respectively, for cortical cultures prepared from wild type mice in the Rab3A+/Ebd colony. Treatment with TTX for 48 hours leads to a significant increase in the average mEPSC amplitude of 23 TTX-treated cells when compared to 20 untreated cells (Figure 2C, CON, 11.0 ± 0.6 pA; TTX, 15.0 ± 1.3 pA, p = 0.02, Kruskal-Wallis test). We note here that while the two strains respond very similarly to activity blockade, the mean mEPSC amplitude in untreated cultures was significantly different in the two wild type strains, 13.9 ± 0.7 pA, wild type from Rab3A+/- colony; 11.0 ± 0.6 pA, wild type from Rab3A+/Ebd colony (p = 0.004, Kruskal-Wallis test).

Block of activity with 48 hours TTX treatment increased mEPSC amplitudes in cortical cultures from wild type mice but not Rab3AEbd/Ebd mice, which had an increased mEPSC amplitude at baseline. A) 10 second current traces recorded at -60 mV in pyramidal cortical neurons from an untreated (“CON”) and TTX-treated (“TTX”) neuron in cultures prepared from Rab3A+/+ mice from the Rab3A+/Ebd colony. B) Average traces for the recordings shown in (A). C) Box and whisker plots for average mEPSC amplitudes from 20 CON and 23 TTX cells (Rab3A+/+, Rab3A+/Ebd colony, mean mEPSCs amplitudes, CON, 11.0 ± 0.6 pA; TTX, 15.0 ± 1.3 pA, p = 0.02, Kruskal-Wallis test). D) Current traces recorded in pyramidal-shaped cortical neurons from a CON and TTX-treated neuron in cultures prepared from Rab3AEbd/Ebd mice. E) Average traces for the recordings shown in (D). F) Box and whisker plots of average mEPSC amplitudes from 21 CON and 22 TTX cells (Rab3AEbd/Ebd, mean mEPSC amplitudes, CON, 15.1 ± 1.0 pA; TTX, 14.6 ± 1.1 pA, p = 0.81, Kruskal-Wallis test). G) Analysis of pooled mEPSC amplitudes for cultures from Rab3A+/+ mice of the Rab3A+/Ebd colony. (i) CDFs of quantiles pooled from CON (black curve) and TTX-treated cultures (red curve) (KS test, D = 0.177, p value = 2.88 * 10-9). (ii) Quantiles pooled from CON and TTX-treated cells were ranked from smallest to largest, plotted against each other, and fit with linear regression function. Gray line, identity; red dashed line, linear regression fit. Inset: full range of the data. (iii) The same ranked mEPSC quantiles were paired and the ratio of CON/TTX computed for each pair, then plotted as a function of the CON mEPSC for the pair. Dashed line at ratio = 1 (no plasticity). H) Analysis of pooled mEPSC amplitudes for cultures from Rab3AEbd/Ebdmice. (i) CDFs of quantiles pooled from CON (black curve) and TTX- treated cultures (red curve). (KS test, D = 0.067, p value = 0.097 (TTX mEPSCs were slightly reduced)). (ii) Quantiles pooled from CON and TTX-treated cells were ranked from smallest to largest, plotted against each other, and fit with a linear regression function. Gray line, identity; red dashed line, linear regression fit. Inset: full range of the data. (iii) The same ranked mEPSC quantiles were paired and the ratio of CON/TTX computed for each pair, then plotted as a function of the CON mEPSC for the pair. Dashed line at ratio = 1 (no plasticity).

We found a complete disruption of homeostatic plasticity in cortical cultures prepared from Rab3AEbd/Ebd mice, as can been seen in viewing example mEPSC traces and average mEPSC waveforms (Figures 2D and E, respectively). The lack of TTX effect was confirmed in a comparison of mEPSC amplitude means for 21 untreated and 22 TTX-treated cells (Figure 2F, CON, 15.1 ± 1.0 pA vs. TTX, 14.6 ± 1.1 pA, p = 0.81, Kruskal-Wallis test). For cultures prepared from the wild type strain of the Earlybird heterozygote colony, the CDFs of the pooled quantiles from control and TTX-treated cells were significantly different (Figure 2Gi, D = 0.177, p value = 2.88 * 10-9, KS test), but for cultures prepared from Rab3AEbd/Ebd mice, the CDF of mEPSC amplitudes from TTX treated cells was shifted to the left, indicating a slight reduction after activity blockade, which did not reach statistical significance (Figure 2Hi; D = 0.067, p = 0.097, KS test). The rank ordered data indicate that the increased slope after TTX observed in the wildtype data (Figure 2Gii, slope 1.73, R2 = 0.993) was completely abolished in the Earlybird data (Figure 2Hii, slope 1.03, R2 = 0.989). Finally, in the wild type data, the ratio of TTX mEPSC amplitude/CON mEPSC amplitude increases from 1.08 for the smallest amplitudes to over 1.60 for the largest amplitudes (Figure 2Giii), whereas for mEPSC amplitudes from Rab3AEbd/Ebdcultures, the ratio of TTX/CON actually falls below 1 (minimum 0.9, Figure 2 Hiii).

Our results show that the homeostatic increase of mEPSC amplitude after activity blockade is disrupted in both Rab3A-/- and the Rab3AEbd/Ebd cortical neurons, strongly supporting a crucial role of functioning Rab3A in the synaptic scaling process. However, it is important to note that the disruption differs for Rab3A-/- and Rab3AEbd/Ebd. In the Rab3A-/- data set, mEPSCs from untreated cultures were indistinguishable from mEPSCs from Rab3A+/+ untreated cultures, demonstrating loss of Rab3A has no impact on basal activity mEPSC amplitudes, but the increase in mEPSC amplitudes after activity blockade was strongly diminished in Rab3A-/-neurons. In the Rab3AEbd/Ebd data set, mEPSC amplitudes from untreated cultures were significantly larger than those of untreated cultures from wild type mice, as can be seen in Figure 3 for the CDFs (Figure 3A, D = 0.245; p = 1.52 * 10-16, KS test) and the box and whisker plot (Figure 3A Inset, Rab3A+/+, 11.0 ± 0.7 pA, vs. Rab3AEbd/Ebd, 15.1 ± 1.0 pA p = 0.0027). The linear fit of the rank ordered data of Rab3AEbd/Ebdvs. wildtype mEPSCs had a slope of 1.63 (Figure 3B). In the ratio plot (Figure 3C), scaling was more uniform than in any other data set we have examined so far, with values starting at 1.4, declining to 1.2, then rising again to 1.4 for the majority of the data, instead of beginning near 1 and rising to 1.5, as happened for homeostatic plasticity (Figure 1Giii). Given this distinct behavior, it cannot be stated that the presence of mutated Rab3A causes the identical effect on mEPSC amplitudes as that of activity blockade. It may be that the much longer time in the presence of the Rab3A Earlybird mutant compared to the 48 hour TTX treatment leads to a stable increase in the smallest mEPSCs.

The increase in mEPSC amplitudes in neurons from cultures prepared from Rab3AEbd/Ebd mice, relative to those from Rab3A+/+ mice in the same strain, showed less divergent scaling than that of the activity-block induced increase in mEPSC amplitudes typically observed in cultures from Rab3A+/+ mice. A) Cumulative distribution function of quantiles pooled from Rab3A+/+ (black curve) and Rab3AEbd/Ebd cultures (red curve), D = 0.245; p = 1.52 * 10-16, KS test. Inset: box and whisker plot of average mEPSC amplitudes from 20 Rab3A+/+ and 21 Rab3AEbd/Ebdcells, mean mEPSC amplitudes, Rab3A+/+, 11.0 ± 0.7 pA; Rab3AEbd/Ebd, 15.1 ± 1.0 pA, p = 0.0027, Kruskal-Wallis test. B) Quantiles pooled from Rab3A+/+ and Rab3AEbd/Ebd cells were ranked from smallest to largest, plotted against each other, and fit with linear regression functions. Gray line, identity; red dashed line, linear regression fit. C) The same ranked mEPSC quantiles were paired and the ratio of Rab3A+/+/Rab3AEbd/Ebd computed for each pair, then plotted as a function of the Rab3A+/+ mEPSC for the pair. Dashed line at ratio = 1.0 (no difference).

Alternatively, these are two distinct mechanisms of mEPSC amplitude augmentation that occlude each other. In any case, the increase in mEPSC amplitude in cultures from Rab3AEbd/Ebdmice is consistent with the increase in mEPC amplitude we observed at the Rab3AEbd/Ebd NMJ (Wang et al., 2011).

The small GTPase Rab3A is generally thought to function presynaptically to regulate synaptic vesicle trafficking, possibly in an activity-dependent manner (Castillo et al., 1997; Lonart et al., 1998; Leenders et al., 2001; Schluter et al., 2006; Coleman and Bykhovskaia, 2009; Tian et al., 2012). In contrast, homeostatic plasticity of mEPSC amplitude has been attributed to an increase in postsynaptic receptors on the surface of the dendrite (O’Brien et al., 1998; Turrigiano et al., 1998). Is Rab3A required for the increase in surface AMPA-type glutamate receptors that has been confirmed by multiple studies (see Introduction)? It should be noted that at the NMJ in vivo we could find no evidence for an increase in AChRs after TTX block of the sciatic nerve in vivo (Wang et al., 2005). The type of AMPA receptor that has been shown to be increased, calcium-permeable GluA1-containing AMPA receptor or calcium-impermeable GluA2-containing AMPA receptor, appears to depend on the experimental manipulation. For example, block of activity by APV (NMDA antagonist) combined with TTX (inhibitor of Na channels) induces an increase in mEPSC amplitude in dissociated hippocampal cultures that is completely reversed by acute application of the Ca-permeable receptor-specific inhibitor NASPM (Sutton et al., 2006), indicating the entirety of the homeostatic effect is due to those receptors, likely GluA1. In contrast, the increase in mEPSC amplitude in mouse hippocampal slice cultures induced by TTX alone was not affected by another GluA1-specific inhibitor, philanthotoxin, suggesting mediation by either Ca-impermeable receptors or a presynaptic effect on quantal size (Soden and Chen, 2010); see also, (Dubes et al., 2022). Because we used TTX alone to block network activity, we expected that NASPM would not reverse the TTX- induced increase in mEPSC amplitude in our mouse cortical cultures, and we found that this was indeed the case. Figure 4A shows that the TTX-induced increase in mean mEPSC amplitude was nearly identical in a set of 11 CON and 11 TTX-treated cells before and after NASPM treatment (before NASPM, CON 12.9 ± 3.5 pA; TTX, 17.5 ± 3.1 pA, p = 0.009; after NASPM, CON 11.9 ± 2.6 pA; TTX 16.1 ± 3.5 pA, p = 0.006, Kruskal-Wallis test). Application of NASPM caused a modest decrease in mEPSC amplitude in both untreated and TTX-treated cultures (Figure 4B, CON, before NASPM, 12.9 ± 3.5 pA; after, 11.9 ± 2.6 pA, p = 0.08; TTX, before NASPM, 17.5 ± 3.1 pA; after, 16.1 ± 3.5 pA, p = 0.08, paired t test), and also decreased mEPSC frequency (Figure 4C, CON, before NASPM, 1.84 sec-1; after NASPM, 1.56 sec-1; p = 0.003, paired t-test; TTX, before NASPM, 4.40 ± 3.51 mEPSCs sec-1; after NASPM, 2.68 ± 2.25 mEPSCs sec-1, p = 0.02, paired t-test), indicating that effective concentrations of NASPM were reached and mEPSCs due to AMPA receptors composed entirely of Ca-permeable subunits were abolished.

Ca2+-permeable, NASPM-sensitive glutamate receptors did not mediate the TTX- induced increase in mEPSC amplitude in cultures of dissociated mouse cortical neurons. A) mEPSC amplitudes for 11 mouse cortical neurons from untreated cultures (“CON”) and 11 neurons from TTX treated cultures (“TTX”) were recorded before and after application of 20 μM NASPM. mEPSC amplitudes were still significantly larger in TTX-treated cultures compared to untreated control cultures after NASPM treatment (TTX effect pre-NASPM, p = 0.009; TTX effect post-NASPM, p = 0.006, Kruskal-Wallis non-parametric test). B) Acute application of NASPM had a modest effect on mEPSC amplitude in cells from CON or TTX cultures that did not reach statistical significance (CON, pre-NASPM, 12.9 ± 1.1 pA, post-NASPM, 11.9 ± 0.8 pA, p = 0.08; TTX, pre-NASPM, 17.5 ± 0.9 pA, post-NASPM, 16.1 ± 1.0 pA, p = 0.08; student’s paired t-test). C) Acute application of NASPM consistently reduced mEPSC frequency in cells from untreated and cells from TTX-treated cultures (CON, pre-NASPM, 1.84 ± 1.82 mEPSCs sec-1, post-NASPM, 1.56 ± 1.74 mEPSCs sec-1; p = 0.003; TTX, pre-NASPM, 4.40 ± 3.51 mEPSCs sec-1, post-NASPM, 2.68 ± 2.25 mEPSCs sec-1, p = 0.02; student’s paired t-test).

Having established that GluA1 receptors are not contributing to the homeostatic increase in mEPSC amplitude, we turned to immunohistochemistry and confocal imaging to assess whether GluA2 receptor expression was increased in our wild type mouse cortical cultures following 48 hour treatment with TTX. Since mEPSCs necessarily report synaptic levels of receptors, we used VGLUT1-positivity to identify synapses on pyramidal primary apical dendrites labeled with MAP-2 immunofluorescence. Figure 5 shows 6 pairs of VGLUT1- and GluA2-immunofluorescent clusters (white frames) along each of two primary dendrites, one in an untreated cortical culture (CON, left), the other in a culture treated with TTX for 48 hours (TTX, right), both from Rab3A+/+ mice. A dendrite typically required ∼10 confocal sections to be fully captured, and the total number of synaptic pairs for all the sections imaged on a dendrite was usually < 20, so this is an atypically high number of pairs within a single section; these particular dendrites and sections were selected for illustration purposes. In addition to the synaptic pairs, we observed many GluA2-immunoreactive clusters not associated with VGLUT1 immunoreactivity; those along the dendrites may be extrasynaptic receptors, and those outside the dendrites may be on astrocytes (Fan et al., 1999). There are also GluA2 immunoreactive clusters that are close to VGLUT1 immunoreactivity but are not located along any apparent MAP2-positive neurite, suggesting axon-axonal contacts, although VGLUT1 has also been detected in astrocytes (Ormel et al., 2012). Only sites that contained both VGLUT1 and GluA2 immunoreactivity close to the primary MAP2-postive dendrite or a secondary branch were selected for analyses.

Confocal images of GluA2 and VGLUT1 immunofluorescence on MAP2-positive primary dendrites of pyramidal-shaped neurons in high density dissociated cortical cultures prepared from Rab3A+/+ mice. Shown are 5X confocal sections of two primary dendrites, one from an untreated culture (CON) and the other from a culture treated with TTX (TTX). Postsynaptic GluA2-positive receptor clusters were identified by their proximity to VGLUT1 immunoreactive presynaptic terminals within the same confocal section. GluA2-VGLUT1 synaptic pairs, which were of maximal brightness in these confocal slices, are indicated by white trapezoids. It can be seen that substantial numbers of GluA2-positive clusters located on the dendrites are unopposed by VGLUT1-positive terminals, and that some GluA2-positive and VGLUT1-positive clusters are found away from any MAP2 positivity, suggesting localization to axons or astrocytes. Scale bar = 10 μm.

Variability in the magnitude of the homeostatic response from culture to culture is averaged out in physiological experiments by the pooling of data from many cells across many cultures. To reduce the necessity for many cultures, we chose to pair experiments in the same cultures by recording mEPSCs from one set of coverslips, and processing another set of coverslips from the same culture for immunohistochemistry. We completed this matched paradigm of physiology and immunohistochemistry, which to our knowledge has never been done before, for 3 cultures prepared from Rab3A+/+ mice and 3 cultures prepared from Rab3A-/- mice. We present the results for mEPSC data and imaging data pooled across the 3 experiments in Figure 6 (Rab3A+/+), Figure 7 (Rab3A-/-), and Table 1. Levels of GluA2 immunoreactivity at synaptic sites were quantified by the size of the GluA2-positive receptor cluster and the average intensity value of the receptor cluster. We analyzed the imaging data the same way we analyzed mEPSC data, sampling 30 quantiles from each dendrite’s data set, and calculating the means, pooling quantiles across dendrites to create CDFs, and sorting quantiles from smallest to largest to produce rank order and ratio plots.

GluA2 receptor cluster size and average intensity increased following activity blockade in parallel with mEPSC amplitudes recorded in the same cortical cultures prepared from Rab3A+/+ mice, but the effects on receptor levels were less robust and differed from mEPSC amplitudes in characteristics of scaling. A) Cumulative distribution functions (CDFs) of mEPSC amplitudes recorded from untreated (“CON”) and TTX-treated cells (“TTX”) pooled from 3 cultures were significantly different (KS test, D = 0.162, p = 1.42 * 10-8, CON, n = 690, TTX, n = 720). Inset, mean mEPSC amplitudes for 23 CON cells and 24 TTX cells were significantly different (CON, 13.7 ± 4.5 pA; TTX, 16.4 ± 4.3 pA, p = 0.016, Kruskal-Wallis test). B) 24 mEPSC amplitude quantiles from each of 23 untreated cells were sorted from smallest to largest and plotted vs. 23 mEPSC amplitude quantiles from each of 24 TTX-treated cells, and the resulting relationship fit with a linear regression equation, Y = 1.106x + 1.22 (R2 = 0.984). Inset, full range of data used for fit. C) Sorted data from (B) were used to calculate the ratio of TTX mEPSC amplitude/CON mEPSC amplitude and the ratios plotted vs. CON mEPSC amplitude. D) CDFs of size of GluA2 receptor clusters imaged from primary dendrites in the same untreated cultures and TTX-treated cultures as in (A) – (C) were significantly different (KS test, D = 0.089, p = 0.002, CON, n = 870, TTX, n = 870). Inset, mean GluA2 cluster size for 29 CON dendrites and 29 TTX dendrites were not significantly different (CON, 0.97 ± 0.38 μm2; TTX, 1.15 ± 0.59 μm2, p = 0.44, Kruskal-Wallis test). E) 29 GluA2 cluster size quantiles from each of 29 untreated dendrites were sorted from smallest to largest and plotted vs. 29 cluster size quantiles from each of 29 TTX-treated dendrites, and the resulting relationship fit with a linear regression equation, Y = 1.328x -0.136 (R2 = 0.98). Inset, full range of data used for fit. F) Sorted data from (E) were used to calculate the ratio of TTX GluA2 cluster size/CON GluA2 cluster size and the ratios plotted vs. CON GluA2 cluster size. G) CDFs of intensity of GluA2 receptor clusters imaged from primary dendrites in untreated cultures and TTX-treated cultures were significantly different (KS test, D = 0.120, p = 6.73 * 10-6, CON, n = 870, TTX, n = 870). Inset, mean GluA2 cluster intensity for 29 CON dendrites and 29 TTX dendrites were not significantly different (CON, 673 ± 90 a.u.; TTX, 687 ± 72 a.u., p = 0.25, Kruskal-Wallis test). H) 29 GluA2 cluster intensity quantiles from 29 untreated dendrites were sorted from smallest to largest and plotted vs. 29 cluster intensity quantiles from 29 TTX-treated dendrites, and the resulting relationship fit with a linear regression equation, Y = 0.899x + 82 (R2 = 0.987). Inset, full range of data used for fit. I) Sorted data from (H) were used to calculate the ratio of TTX cluster intensity/CON cluster intensity and plotted vs. CON cluster intensity.

GluA2 receptor cluster size, average intensity, and mEPSC amplitudes recorded in the same cortical cultures prepared from Rab3A-/- mice were not increased following activity blockade. A) Cumulative distribution functions (CDFs) of mEPSC amplitudes recorded from untreated (“CON”) and TTX-treated cells (“TTX”) pooled from 3 Rab3A-/- cultures were not significantly different (KS test, D = 0.072, p = 0.085, CON, n = 630, TTX, n = 570). Inset, mean mEPSC amplitudes for 21 CON cells and 19 TTX cells were not significantly different (CON, 14.9 ± 3.8 pA; TTX, 14.0 ± 4.0 pA, p = 0.34, Kruskal-Wallis test). B) 19 mEPSC quantiles from each of 21 untreated cells were sorted from smallest to largest and plotted vs. 21 mEPSC quantiles from each of 19 TTX-treated cells, and the resulting relationship fit with a linear regression equation, Y = 0.93x + 0.13 (R2 = 0.997). Inset, full range of data used for fit. C) Sorted data from (B) were used to calculate the ratio of TTX mEPSC amplitude/CON mEPSC amplitude and plotted vs. CON mEPSC amplitude. D) CDFs of GluA2 receptor cluster size imaged from primary dendrites in the same untreated cultures and TTX-treated cultures as in (A) – (C) were not significantly different (KS test, D = 0.045, p = 0.31, CON, n = 900, TTX, n = 870). Inset, mean GluA2 cluster size for 30 CON dendrites and 29 TTX dendrites were not significantly different (CON, 0.93 ± 0.27 μm2; TTX, 0.91 ± 0.28 μm2, p = 0.74, Kruskal-Wallis test). E) 29 GluA2 cluster size quantiles from each of 30 untreated dendrites were sorted from smallest to largest and plotted vs. 30 cluster size quantiles from each of 29 TTX-treated dendrites, and the resulting relationship fit with a linear regression equation, Y = 1.01x -0.03 (R2 = 0.993). Inset, full range of data used for fit. F) Sorted data from (E) were used to calculate the ratio of TTX GluA2 cluster size/CON GluA2 cluster size and the ratios plotted vs. CON GluA2 cluster size. G) CDFs of GluA2 receptor cluster average intensity imaged from primary dendrites in untreated cultures and TTX-treated cultures were not significantly different (KS test, D = 0.0.060, p = 0.08, CON, n = 900, TTX, n = 870). Inset, mean GluA2 cluster intensity for 30 CON dendrites and 29 TTX dendrites were not significantly different (CON, 766 ± 68 a.u.; TTX, 776 ± 79 a.u., p = 0.47, Kruskal-Wallis test). H) 29 GluA2 cluster intensity quantiles from 30 untreated dendrites were sorted from smallest to largest and plotted vs. 30 cluster intensity quantiles from 29 TTX-treated dendrites, and the resulting relationship fit with a linear regression equation, Y = 1.09x - 60 (R2 = 0.993). Inset, full range of data used for fit. I) Sorted data from (H) were used to calculate the ratio of TTX cluster intensity/CON cluster intensity and the ratios plotted vs. CON cluster intensity.

Comparison of effect of activity blockade with TTX on mEPSC amplitude and GluA2 receptor cluster characteristics. Voltage clamp recordings of mEPSCs and confocal imaging of GluA2 receptor immunofluorescence at VGLU1-positive synapses were performed on coverslips from the same 3 Rab3A+/+ cultures (left) and 3 Rab3A-/- cultures (right). p values are not corrected for multiple comparisons.

In the data pooled from the 3 matched experiments, neurons from Rab3A+/+ cultures showed a significant increase in mean mEPSC amplitudes following activity blockade (Figure 6A inset, CON, 13.7 ± 4.5 pA, n = 23; TTX, 16.4 ± 4.3 pA, n = 24; p = 0.016, Kruskal-Wallis test); a significant shift of the TTX CDF to larger mEPSC amplitude values (Figure 6A, D = 0.162, p = 1.42 * 10-8, KS test); a slope of 1.11 in the rank order plot (Figure 6B, R2 = 0.98); and a mean ratio of the 50th to 75th percentile of 1.24 in the ratio plot (Figure 6C). These data indicate that the homeostatic response averaged across the 3 cultures was very similar to, but slightly smaller than, that of the previous data set presented in Figure 1. The means for size and intensity of GluA2 receptor clusters, while showing trends to higher values after activity blockade, were not significantly different (Figure 6D inset, size, CON, 0.97 ± 0.38 μm2; TTX, 1.15 ± 0.59 μm2, p = 0.44; Figure 6G inset, intensity, CON, 673 ± 90, TTX, 687 ± 72, p = 0.25, Kruskal-Wallis test; Table 1). However, when comparing the CDFs, where the sample size is much greater (equal to usually close to 30 quantiles * number of cells), the shifts to larger values for the TTX CDF did reach statistical significance for GluA2 receptor cluster size (Figure 6D, D = 0.089, p = 0.002, KS test) and intensity (Figure 6G, D = 0.120, p = 6.73 * 10-6, KS test). In the rank order plot for GluA2 receptor cluster sizes (Figure 6E), the slope was actually larger than that for mEPSC amplitude (1.33), but for intensity (Figure 6H), the slope value was below 1 (0.90). The weaker effects on synaptic GluA2 receptor levels relative to mEPSC amplitude suggest the possibility that GluA2 receptors are not the sole determining factor in the increased mEPSC amplitude following activity blockade. In addition, the plot of the ratio of TTX GluA2 receptor cluster size/control size as a function of ranked control values showed a surprising deviation from the divergent scaling we have observed for mEPSC amplitude (compare Figure 6C to Figure 6F). The smallest GluA2 receptor clusters showed a TTX/CON ratio as high as 3, but the mean ratio across the 50th to 75th percentile settled down to 1.11 within < 20 samples. This finding of larger ratios for the smallest GluA2 receptor clusters is supported by a previous study that followed the same fluorescently labeled postsynaptic sites over time following activity blockade. Wang and colleagues found that the biggest increases after activity blockade occurred at the smallest/dimmest synaptic sites (Wang et al., 2019). A similar inverse relationship was reported for changes in individual fluorescently labeled postsynaptic sites followed over time under normal activity conditions (Minerbi et al., 2009; Statman et al., 2014). The mismatch in ratio values between the smallest mEPSC amplitudes and smallest GluA2 receptor cluster sizes may be due to the inability of physiological assays to detect the mEPSCs coming from the smallest synapses.

We proceeded to determine the homeostatic responses of mEPSC amplitude and GluA2 receptor levels in 3 cultures prepared from Rab3A-/- mice. As shown for the data set in Figure 1, mean mEPSC amplitude was not increased following activity blockade in the data pooled from this new set of 3 Rab3A-/- cultures (Figure 7A inset, CON 14.9 ± 3.8 pA, n = 21; TTX, 14.0 ± 4.0 pA, n = 19; p = 0.34, Kruskal-Wallis test). Notably, the TTX CDF was shifted to smaller mEPSC amplitude values (Figure 7A, D = 0.072, p = 0.085, KS test), the slope of the rank order plot was below 1 (Figure 7B, slope 0.93, R2 = 0.997), and the mean ratio across the 50th to 75th percentile was below 1 (Figure 7C, 0.96). Thus, the disruption of homeostatic plasticity of mEPSC amplitude shown with the previous data sets was recapitulated in this data set. For immunohistochemistry results from the same Rab3A-/- cultures, we found that mean GluA2 receptor cluster size was unchanged following activity blockade (Figure 7D inset, CON, 0.93 ± 0.27 μm2, TTX, 0.91 ± 0.28 μm2, p = 0.74, Kruskal-Wallis test), the TTX CDF was not significantly shifted (Figure 7D, D = 0.045, p = 0.31, KS test), the slope on the rank order plot was 1.01 (Figure 7E), and the ratio hovered at or below 1 except for a small group of ratios at the high end of the control values (Figure 7F). For intensity, mean values were not increased after activity blockade (Figure 7G inset, CON, 766 ± 68, TTX, 776 ± 79, p = 0.47, Kruskal-Wallis test), the CDFs were not significantly different (Figure 7G, D = 0.060, p = 0.080), the slope on the rank order plot was 1.09 (Figure 7H), and the mean ratio value across the 50th to 75th percentile was 1.02 (Figure 7I). Taken together, these results indicate that the modest increases in GluA2 receptor cluster size and intensity following activity blockade observed in wild-type cultures do not occur in the absence of Rab3A.

As noted above, the magnitude of the homeostatic effect on mEPSC amplitude appeared to be more robust than that on receptor levels in cultures from Rab3A+/+ mice. Further evidence that there is not a one-to-one correspondence in the homeostatic response of mEPSC amplitudes and GluA2 receptor levels was apparent when the individual experiments’ data were compared. The electrophysiology experiments were technically difficult because in order to evaluate whether a particular culture displayed homeostatic plasticity of mEPSC amplitude, we set a minimum requirement of recording from 6 cells per condition, meaning a successful experiment required at least 12 (6 CON and 6 TTX) usable recordings (i.e. low holding current, stable baseline, low noise, etc.) in one day. Ultimately, our sample size ranged from 6 to 10 cells per condition per culture. Figures 8A-C show the CDFs and ratio plots for 3 individual Rab3A+/+ cultures. Rank order plots are not included here because the ratio plot is more sensitive than the rank order plot when determining effect magnitude. We found that the mEPSC amplitude effects differed from the GluA2 cluster size effects at both the broad level (increases vs decreases), and in the specific ways in which the effect was non-uniform. Broadly, although two of the three Rab3A+/+ experiments show a homeostatic increase in the TTX CDF and a ratio of TTX/CON > 1 for both mEPSC amplitudes (Figures 8Ai,iii and 8Bi,iii) and synaptic GluA2 receptor cluster sizes (Figures 8Aii,iv and 8Bii,iv), Culture #3 shows an increase in mEPSC amplitude but a decrease in the GluA2 receptor cluster size, based on both the shifts in CDFs (mEPSC, Figure 8Ci, vs GluA2 receptor cluster size, Figure 8Cii) and the majority of ratio values (mEPSC, Figure 8Ciii, vs. GluA2 receptor cluster size, Figure 8Civ). Similarly, we found an obvious mismatch between the mEPSC amplitude data and the GluA2 cluster size data for one of the three Rab3A-/- experiments. For Culture #3, the mEPSC TTX CDF shows a shift to the left (Figure 8Fi), and the ratios were below 1 for control mEPSC amplitudes < 20 pA and above 1 for control mEPSC amplitudes > 20 pA (Figure 8Fiii). In the same culture, GluA2 receptor cluster size CDF shows a shift to the right after activity blockade (Figure 8Fii), and in the ratio plot, values were at or above 1 throughout the range of control sizes (Figure 8Fiv).

mEPSC amplitudes and GluA2 receptor cluster sizes did not always change in parallel after treatment with TTX in cultures from Rab3A+/+ and Rab3A-/- mice. A) Rab3A+/+ Culture #1: mEPSC amplitudes and GluA2 receptor cluster sizes increased in parallel after treatment in TTX. i) CDFs of mEPSC amplitudes of cells from untreated and TTX treated dishes were significantly different (KS test, D = 0.16, p = 0.010; CON, n = 180, TTX, n = 210). ii) CDFs of GluA2 receptor cluster sizes of dendrites from untreated and TTX treated dishes were significantly different (KS test, D = 0.17, p = 3.79 * 10-4; CON, n = 300, TTX, n = 300). iii) 35 mEPSC amplitude quantiles from each of 6 untreated cells and 30 quantiles from each of 7 TTX-treated cells were sorted from smallest to largest and used to calculate the ratio of TTX mEPSC amplitude/CON mEPSC amplitude, which was plotted vs. CON mEPSC amplitude. iv) 30 GluA2 receptor cluster size quantiles from each of 10 untreated dendrites and 30 GluA2 receptor cluster size quantiles from each of 10 TTX-treated dendrites were used to calculate the ratio of TTX GluA2 cluster size/CON GluA2 cluster size, and plotted vs. CON GluA2 cluster size. Inset, full range of data. B) Rab3A+/+ Culture #2: mEPSC amplitudes and GluA2 receptor cluster sizes increased in parallel after treatment in TTX. i) CDFs of mEPSC amplitudes of cells from untreated and TTX treated dishes were significantly different (KS test, D = 0.21, p = 1.06 * 10-4; CON, n = 210, TTX, n = 240). ii) CDFs of GluA2 receptor cluster sizes of dendrites from untreated and TTX treated dishes were significantly different (KS test, D = 0.21, p = 5.39 * 10-6; CON, n = 270, TTX, n = 300). iii) 32 mEPSC amplitude quantiles from each of 7 untreated cells, and 28 mEPSC amplitude quantiles from each of 8 TTX treated cells were sorted from smallest to largest and were used to calculate the ratio of TTX mEPSC amplitude/CON mEPSC amplitude, and plotted vs. CON mEPSC amplitude. iv) 30 GluA2 receptor cluster size quantiles from each of 9 untreated dendrites and 27 GluA2 receptor cluster size quantiles from each of 10 TTX-treated dendrites were sorted from smallest to largest and used to calculate the ratio of TTX GluA2 cluster/CON GluA2 cluster size, and plotted vs. CON GluA2 cluster size. Inset, full range of data. C) Rab3A+/+ Culture #3: mEPSC amplitudes increased, but GluA2 receptor cluster sizes decreased, after treatment in TTX. i) CDFs of mEPSC amplitudes of cells from untreated and TTX treated dishes were significantly different (KS test, D = 0.15, p = 0.002; CON, n = 300, TTX, n = 270). ii) CDFs of GluA2 receptor cluster sizes of dendrites from untreated and TTX treated dishes were significantly different, but the TTX CDF is shifted to the left (KS test, D = 0.16, p = 0.001; CON, n = 300, TTX, n = 270) iii) 27 mEPSC amplitude quantiles from each of 10 untreated cells, and 30 mEPSC amplitude quantiles from each of 9 TTX-treated cells, were sorted from smallest to largest, used to calculate the ratio of TTX mEPSC amplitude/ CON mEPSC amplitude, and plotted vs. CON mEPSC amplitude. iv) 27 GluA2 receptor cluster size quantiles from each of 10 untreated dendrites and 30 GluA2 receptor cluster sizes from each of 9 TTX-treated dendrites were sorted from smallest to largest, used to calculate the ratio of TTX GluA2 receptor cluster size/CON GluA2 receptor cluster size, and plotted vs. CON GluA2 receptor cluster size. D) Rab3A-/- Culture #1 (note, this culture is distinct from that of Rab3A+/+ Culture #1). Neither mEPSC amplitudes, nor GluA2 receptor cluster size, increased following TTX treatment. i) CDFs of mEPSC amplitudes of cells from untreated and TTX-treated cultures were not significantly different (KS test, D = 0.05, p = 0.90; CON, n = 210, TTX, n = 180). ii) CDFs of GluA2 receptor cluster sizes of dendrites from untreated and TTX (Figure 8 continued) treated dishes were not significantly different, although the TTX CDF is shifted slightly to the left (KS test, D = 0.10, p = 0.08; CON, n = 300, TTX, n = 300). iii) 30 mEPSC amplitude quantiles from each of 7 untreated cells and 35 mEPSC quantiles from each of 6 TTX-treated cells were sorted from smallest to largest, used to calculate the ratio TTX mEPSC amplitude/CON mEPSC amplitude, and plotted vs. CON mEPSC amplitude. iv) 30 GluA2 receptor cluster size quantiles from each of 10 untreated dendrites and 10 TTX-treated dendrites were sorted from smallest to largest, used to calculate the ratio TTX cluster size/CON cluster size, and plotted vs. CON cluster size. E) Rab3A-/- Culture #2 (distinct from Rab3A+/+ Culture #2) mEPSC amplitudes and GluA2 receptor cluster sizes decreased after TTX treatment. i) TTX mEPSC amplitudes CDF is shifted slightly to the left compared to the CON mEPSC amplitudes CDF but the difference is not significant (KS test, D = 0.08, p = 0.28, CON, n = 240, TTX, n = 240). ii) TTX GluA2 receptor cluster sizes CDF is significantly shifted to the left compared to CON GluA2 receptor cluster sizes CDF (KS test, D = 0.13, p = 0.01; CON, n = 300, TTX, n = 300). iii) 24 mEPSC amplitude quantiles from each of 8 untreated cells and 8 TTX- treated cells were sorted from smallest to largest, used to calculate the ratio TTX mEPSC amplitude/CON mEPSC amplitude, and plotted vs. CON mEPSC amplitude. iv) 27 GluA2 receptor cluster size quantiles from each of 10 untreated dendrites, and 30 GluA2 receptor cluster size quantiles from each of 9 TTX treated dendrites, were sorted from smallest to largest, used to calculate the ratio TTX GluA2 cluster size/CON cluster size, and plotted vs. CON cluster size. Both mEPSC amplitude and GluA2 receptor cluster size ratios were below 1 across the majority of data. F) Rab3A-/- Culture #3 (distinct from Rab3A+/+ Culture #3). mEPSC amplitudes were decreased following TTX treatment, but GluA2 receptor cluster sizes were increased. i) TTX CDF of mEPSC amplitudes was significantly shifted to the left of CON CDF (KS test, D = 0.21, p = 4.85 * 10-4, CON, n = 180, TTX, n = 180). ii) TTX CDF of GluA2 receptor cluster sizes was significantly shifted to the right of CON CDF (KS test, D = 0.15, p = 0.003, CON, n = 300, TTX, n = 300). iii) 30 mEPSC amplitude quantiles from each of 6 untreated and 6 TTX-treated cells were sorted from smallest to largest, used to calculate the ratio TTX mEPSC amplitude/CON mEPSC amplitude, and plotted vs. CON mEPSC amplitude. iv) 30 GluA2 receptor cluster size quantiles from each of 10 untreated dendrites and 10 TTX-treated dendrites were sorted from smallest to largest, used to calculate the ratio TTX GluA2 cluster size/CON GluA2 cluster size, and plotted vs. CON CluA2 cluster size.

Perhaps more striking than the broad mismatches, which could be attributed to the small sample sizes, is the complete lack of correspondence of TTX/CON ratios in the mEPSC amplitudes compared to GluA2 cluster sizes in the same cultures. In cultures from Rab3A+/+ mice, mEPSC amplitudes show divergent scaling similar to what we have previously reported, with the smallest CON amplitudes showing the smallest ratios, and the ratio increasing monotonically to a plateau, or, a peak followed by a gradual decline, as CON mEPSC amplitude becomes very large. In contrast, the very smallest GluA2 cluster sizes in Rab3A+/+ cultures have the largest increase, with ratios as high as 3 in the two cultures that showed a broad overall increase. The ratios decline dramatically to values around 1.5, but then remain substantially above 1 throughout the entire range of data. Interestingly, the dramatic increase in the size of the smallest GluA2 clusters is absent in the two cultures from Rab3A-/- mice that broadly showed no increase, suggesting this aspect is also dependent on Rab3A. It is not even possible to cut off the ratios at a certain point, to address a threshold effect such as GluA2 cluster sizes below a certain value not being measurable electrophysiologically, and get the ratios for the remaining data to match up. Taken together, these data indicate that variation in homeostatic effects on GluA2 cluster size is not responsible for the variation in homeostatic effects on mEPSC amplitude.

In summary, our matched mEPSC and receptor cluster results indicate that: 1. activity blockade resulted in an increase in mEPSC amplitudes and synaptic GluA2 receptors; 2. loss of Rab3A disrupted both the increase in mEPSC amplitudes and the increase in GluA2 receptor levels; and 3. The homeostatic effects on mEPSC amplitudes and the GluA2 cluster sizes within the same cultures differed at both the broad level, and in the specific way the ratios were non- uniform, making it difficult to conclude modulation of GluA2 cluster size is the sole determinant of homeostatic increases in mEPSC amplitude. Although it remains a possibility that the disparities in homeostatic effects on mEPSC amplitudes and GluA2 cluster size TTX/CON ratios might arise due to a sampling error, given the small number of cells sampled in each individual experiment (≤ 10), it seems unlikely that differences arising from sampling error alone would be so consistent. However, only when we are able to more easily sample larger numbers of cells, or it becomes possible to measure mEPSCs and receptor levels in the same cells, will we be able to resolve this question.

We have demonstrated here for the first time that the synaptic vesicle protein Rab3A influences the homeostatic regulation of postsynaptic GluA2 receptors. What are possible ways that Rab3A could exert a postsynaptic action? It has previously been shown that exogenous addition of TNFα to hippocampal cultures causes an increase in surface expression of GluA1 receptors (although not GluA2 receptors), and that the homeostatic increase in mEPSC amplitudes is abolished in cultures prepared from the TNFα deletion mouse (Stellwagen et al., 2005; Stellwagen and Malenka, 2006). Furthermore, neurons from TNFα+/+ mice plated on astrocytic feeder layers derived from TNFα-/- mice fail to show the increase in mEPSC amplitude after TTX treatment, indicating that the TNFα inducing the receptor increases following activity blockade comes from the astrocytes (Stellwagen and Malenka, 2006). Rab3A has been detected in astrocytes (Maienschein et al., 1999; Hong et al., 2016), so to determine whether Rab3A is acting via regulating TNFα release from astrocytes, we performed experiments similar to those of Stellwagen and Malenka (2006) (schema illustrated in Figure 9, left side). We compared the effect of activity blockade on mEPSC amplitudes recorded from cortical neurons from Rab3A+/+ mice plated on Rab3A+/+ astrocytic feeder layers (Figure 9A); neurons from Rab3A+/+ mice plated on Rab3A-/- astrocytes (Figure 9B), and neurons from Rab3A-/- mice plated on Rab3A+/+ astrocytes (Figure 9C). If Rab3A is required for TNFα release from astrocytes, then Rab3A+/+ neurons plated on Rab3A-/- astrocytes should not show a homeostatic increase in mEPSC amplitude after treatment with TTX, and any cultures with Rab3A+/+ astrocytes should have a normal homeostatic response. We found the opposite result: mEPSC amplitudes increased dramatically in cultures where Rab3A was present in neurons (Figures 9A and B), but much more modestly increased in cultures where Rab3A was present only in astrocytes (Figure 9C).

Rab3A in neurons, not astrocytes, was required for full homeostatic plasticity after activity blockade. A-C, Data from dissociated cortical neurons plated on an astrocyte feeder layer, each prepared separately from the type of mice depicted in the schema at left: A) Neurons from Rab3A+/+ mice plated on astrocytes from Rab3A+/+ mice; i) CDFs for mEPSC amplitudes from untreated (CON, black curve) and TTX-treated (TTX, red curve) cultures (KS test, test statistic, D = 0.141, p = 2.74 * 10-5, n = 510 CON, 600 TTX samples; inset, means, CON, 13.3 ± 0.5 pA, n = 17 cells, TTX, 16.7 ± 1.2 pA, n = 20 cells, p = 0.03, Kruskal Wallis test; ii) rank order plot and linear regression fit, slope = 1.89, intercept = -8.51, R2 = 0.97; iii) ratio plot of TTX mEPSC amplitude/CON amplitude vs. CON amplitude; B) Neurons from Rab3A+/+ mice plated on astrocytes from Rab3A-/- mice; i) CDFs of mEPSC amplitudes, CON, black curve; TTX, red curve (KS test, test statistic 0.273, p = 2.98 * 10-11), n = 330 CON, 330 TTX samples; inset, means, CON, 13.3 ± 1.0 pA, n = 11 cells, TTX, 18.8 ± 1.4 pA, n = 11 cells, p = 0.005, Kruskal-Wallis test; ii) rank order plot and linear regression fit, slope = 1.76, intercept = -4.64, R2 = 0.98; iii) ratio plot of TTX mEPSC amplitude/TTX mEPSC amplitude vs. CON mEPSC amplitude; C) Neurons from Rab3A-/- neurons plated on astrocytes from Rab3A+/+ mice; i) CDFs of mEPSC amplitudes, CON, black curve, TTX, red curve (KS test, test statistic 0.126, p = 0.005), n = 420 CON, 330 TTX samples; inset, means, CON, 15.2 ± 1.1 pA, n = 14 cells, TTX, 16.9 ± 0.7 pA, n = 11 cells, p = 0.125, Kruskal-Wallis test;; ii) rank order plot and linear regression fit; slope = 1.01, intercept = 1.53, R2 = 0.97; ratio plot, TTX mEPSC amplitude/CON mEPSC amplitude vs. CON mEPSC amplitude. Gray line, identity; red dashed line, linear regression fit. Insets in rank order plots depict full range of data.

Interestingly, the homeostatic effects appear to be larger in the feeder layer cultures, compared to the effects observed in our previous neuronal/astrocyte mixed cultures prepared directly onto poly-L-lysine coated coverslips, but the ratio plots clearly show scaling was divergent under both plating conditions (compare Figures 1Giii with Figures 9Aiii and Biii). We noticed that in the astrocyte feeder layer cultures, and in the matched mEPSC amplitude and GluA2 receptor measurement experiments, the mean mEPSC amplitude in the untreated cultures prepared from Rab3A-/- mice was slightly, but not significantly, larger, compared to the mean mEPSC amplitude for untreated cultures prepared from Rab3A+/+ mice in the same experiments (astrocyte feeder layer, Rab3A+/+ neurons on Rab2A+/+ astrocytes, 13.3 ± 0.5 pA, Rab3A-/- on Rab3A+/+ astrocytes, 15.2 ± 1.1 pA; matched mEPSC and receptor experiments, Rab3A+/+ 13.7 ± 4.5 pA, Rab3A-/- 14.9 ± 3.8 pA. It is therefore possible that loss of Rab3A, like expression of the Rab3A Earlybird mutant, affects basal mEPSC amplitude, albeit to a lesser extent. However, it could also be that these differences reflect random variation in mEPSC amplitude from culture to culture.

In summary, in the absence of Rab3A from astrocytes, the divergent scaling of mEPSC amplitude following activity blockade was completely normal, whereas in the absence of Rab3A in neurons, scaling was greatly diminished. This result makes it highly unlikely that Rab3A is required for the release of TNFα, or another factor from astrocytes, that induces a homeostatic upregulation of postsynaptic receptors and thereby increases mEPSC amplitude following TTX treatment. Neuronal Rab3A appears to mediate the homeostatic increase in mEPSC amplitude following activity blockade.

Discussion

We found that homeostatic synaptic plasticity of mEPSC amplitude in dissociated mixed cultures of mouse cortical neurons and astrocytes behaved remarkably similar to the mouse NMJ in response to loss of Rab3A function: scaling up of mEPSC amplitude following prolonged network silencing by TTX was strongly diminished in cultures from Rab3A-/-, and in cultures from Rab3AEbd/Ebd mice, basal mEPSC amplitude was increased compared to that of wild-type cultures and was not further modified following 48-hour treatment with TTX. These results suggest that normal function of the presynaptic vesicle protein Rab3A is required for the homeostatic scaling up of mEPSC amplitude in cortical cultures and at the NMJ in vivo.

Rab3A not likely to regulate receptor trafficking

We demonstrated that an increase in synaptic GluA2 receptors accompanies the increase in mEPSC amplitude in cultures of dissociated mouse cortical neurons after 48-hour TTX treatment. By examining the same cultures with both mEPSC and immunofluorescence measurements, we found: 1. the effect on GluA2 receptor levels, whether measured as the size of a synaptic cluster or the density of receptors in a cluster, was less robust than the effect on mEPSC amplitudes; 2. the homeostatic effect on mEPSC amplitudes and GluA2 receptors substantially differed from each other both broadly (mEPSC amplitudes could increase while GluA2 cluster sizes decreased, and vice versa), and in the details of non-uniform scaling, one example being that the smallest mEPSC amplitudes showed the least effects but the smallest GluA2 clusters showed the greatest effects. Interestingly, there are other examples where the effect of activity blockade by either TTX alone (Hu et al., 2010), TTX and AP5 (Letellier et al., 2014), or DNQX (Blackman et al., 2012) on AMPAR levels did not match that of mEPSC amplitudes, although in these cases, the increases in receptors were larger than the increases of mEPSC amplitudes. The striking differences between TTX effects on AMPARs and mEPSC amplitudes indicates that there may be another factor contributing to mEPSC amplitude besides levels of AMPAR expression. One possibility is alterations in the type of AMPA receptor expressed, since GluA1 homomers have much greater conductance than GluA2-containing receptors (Oh and Derkach, 2005; Benke and Traynelis, 2019). Increases in mEPSC amplitudes not accompanied by increases in receptor numbers have been attributed to such a switch (Hou et al., 2015; Silva et al., 2019; Dubes et al., 2022).

In addition to finding a disparity between synaptic GluA2 receptor and mEPSC responses to prolonged activity blockade, there are two other reasons we do not think Rab3A impacts homeostatic synaptic plasticity solely through postsynaptic regulation of GluA2 receptor trafficking to the plasma membrane. First, prolonged activity blockade with TTX at the NMJ led to an increase in mEPC amplitude that was not accompanied by an increase in AChR levels, yet it was dependent on Rab3A (Wang et al., 2011). Second, for Rab3A to modulate AMPAR trafficking, it must be located in the postsynaptic dendrite. Multiple presynaptic molecules, such as SNARE proteins, synaptotagmins, and NSF have been identified in the postsynaptic compartment (Lledo et al., 1998; Nishimune et al., 1998; Osten et al., 1998; Song et al., 1998; Araki et al., 2010; Kennedy et al., 2010; Suh et al., 2010; Jurado et al., 2013; Hussain and Davanger, 2015; Gu et al., 2016; Wu et al., 2017), but to our knowledge, Rab3A has not.

Rab3A may regulate quantal size presynaptically

Based on our findings, the lack of postsynaptic Rab3A expression and the known presynaptic expression of Rab3A, Rab3A may regulate mEPSC amplitude presynaptically by modulating: a) the amount of transmitter transported into the vesicle, b) the size of the vesicle, or c) fusion pore characteristics. mRNA for the glutamate transporter, VGLUT1, was increased in rat cortical cultures after 48-hr treatment with TTX (De Gois et al., 2005), and VGLUT1 relative to synapsin I in rat hippocampal cultures was increased after 48-hr treatment with AP-5 (Wilson et al., 2005). This activity-dependent regulation of VGLUT1 was recently corroborated in a proteome study (Dorrbaum et al., 2020), and it has also been shown that mEPSC amplitudes are increased after overexpression of VGLUT at the Drosophila NMJ (Daniels et al., 2004) and in rat hippocampal cultures (Wilson et al., 2005). These data suggest that increases in VGLUT1 after activity-blockade could contribute to the homeostatic increase in mEPSC amplitude. Although there is no direct evidence for an interaction between VGLUT1 and Rab3A, Rab3A levels are reduced 50% in hippocampal extracts from mice lacking VGLUT1 (Fremeau et al., 2004). A Rab3A-dependent modulation of presynaptic quantal size would explain the homeostatic effects on mEPSC amplitudes that are not matched by changes in receptor levels.

Interestingly, Rab3A plays a role in regulating vesicle size. In adrenal chromaffin cells from mice lacking all 4 Rab3 isoforms (A, B, C and D) or mice heterozygous for Rab3A and lacking the B, C, and D isoforms (Rab3A+/-BCD-/-) with only 1 copy of Rab3A, there was a 20% increase in large dense core vesicle (DCV) diameter that did not reach significance (p = 0.11) (Schonn et al., 2010). Similarly, in the exocrine pancreas and parotid gland of the Rab3D-/-, large DCV diameter was increased by 22% (Riedel et al., 2002). This magnitude of change in diameter would produce a greater than 80% increase in volume. Also, an increase in volume has been shown to increase miniature endplate junctional current amplitude at the Drosophila NMJ (Karunanithi et al., 2002). There have not been any reports that small synaptic vesicle size is increased in Rab3-/-mice, but small synaptic vesicle size is regulated by activity: vesicle size increased in retinotectal terminals within the optic tectum following enucleation in pigeon (Cuenod et al., 1970), with similar results in cat and monkey (Lloret and Saavedra, 1975), in rat caudate nucleus following cortical ablation (Kawana et al., 1971), and in rat ventral spinal cord following immobilization (Cheresharov et al., 1978).

Finally, we previously showed that at the mouse NMJ, both loss of Rab3A and activity blockade by TTX caused an increase in the frequency of abnormally slow rising and/or prolonged duration mEPCs (Wang et al., 2011). We further demonstrated that the loss of Rab3A increased the frequency of very small amplitude fusion pore feet in mouse chromaffin cells (Wang et al., 2008). Taken together, these data suggested that Rab3A acts to prevent abnormal fusion pore release events. However, at the NMJ and in chromaffin cells, the kinetics (rise time, decay, half width) of the vast majority of release events were unchanged in the absence of Rab3A, making it unlikely that the shift in the distribution of mEPSC amplitudes is mediated by a Rab3A-dependent change in fusion pore kinetics.

Neuronal, not astrocytic, Rab3A may be required for the homeostatic release of a signaling molecule

Our results support the surprising possibility that presynaptic Rab3A regulates postsynaptic AMPARs, rather than occurring through a direct effect on the trafficking of AMPARs by Rab3A in the dendrite. Astrocytic release of TNFα was shown to mediate the increase in AMPARs in prolonged activity-blocked hippocampal cultures (Stellwagen and Malenka, 2006), but we found that loss of Rab3A in astrocytes does not disrupt the increase in mEPSC after activity blockade. Further, we previously showed that the homeostatic increase in NMJ mEPC amplitude was completely normal in the absence of TNFα (Wang et al., 2011).

Another possible way Rab3A might mediate regulation of postsynaptic AMPARs is through the presynaptic release of a signaling molecule, such as brain-derived neurotrophic factor (BDNF), that acts anterogradely to alter postsynaptic AMPAR levels. Addition of exogenous BDNF to neuronal cultures prevents the increase in mEPSC amplitude following activity blockade (Rutherford et al., 1998; Benevento et al., 2016) but see (Smith-Dijak et al., 2019), and BDNF mRNA is reduced after activity blockade in vitro (Benevento et al., 2016; Miyasaka and Yamamoto, 2021) and in vivo (Castren et al., 1992). Also, after reduction of secreted BDNF in culture media via the BDNF scavengers TrkB-FC or TrkB-IgG, mEPSC amplitudes are increased (Rutherford et al., 1998; Benevento et al., 2016; Smith-Dijak et al., 2019). Rab3B, a Rab3 family member with strong homology to Rab3A but expressed more highly in inhibitory nerve terminals (Tsetsenis et al., 2011), is modified epigenetically following TTX treatment by the same DNA methylation pathway as that for BDNF (Benevento et al., 2016). In addition, reduction of Rab3A in astrocytes decreases astrocytic BDNF release (Hong et al., 2016), and treatment of cultures with BDNF can increase Rab3A levels (Takei et al., 1997; Thakker-Varia et al., 2001) but see (Shinoda et al., 2014). Taken together, these results suggest an interplay between BDNF and Rab3A in the activity-dependent regulation of AMPA receptors.

Rab3A is required for divergent scaling

We have established that homeostatic synaptic plasticity in mouse cortical cultures, rat cortical cultures, and mouse hippocampal cultures demonstrates divergent, rather than uniform, scaling (Hanes et al., 2020; Koesters et al., 2022). The deviation from uniformity suggests that individual synapses may be independently regulated, and do not require a global, cell-wide signal. Our conclusion is supported by multiple studies demonstrating local activity- dependent regulation of mEPSC amplitudes at subsets of synapses (Ju et al., 2004; Sutton et al., 2006; Hou et al., 2008; Beique et al., 2011; Hou et al., 2011; Letellier et al., 2014). It is further supported by the recent finding that in hippocampal cultures treated with TTX, AMPA receptors increase selectively at synaptopodin-positive sites, which are already larger and contain higher amounts of AMPA receptors (Dubes et al., 2022). The authors provide evidence that the process depends on the microRNA miR-124: miR-124 is reduced after activity blockade; and overexpression of miR-124 blocks the TTX-induced increase in AMPA receptors and mEPSC amplitudes. Increases in miR-124 have also been linked to decreases in mEPC amplitudes at the NMJs of the slow channel syndrome mouse (mSCS) in which prolongation of ACh currents causes excess excitation (Zhu et al., 2013). Zhu and colleagues provide evidence for a postsynaptic calpain-Cdk5-nNOS pathway that is activated by increased calcium levels, with the NO feeding back retrogradely to regulate miR-124, which in turn directly interacts with Rab3A to affect quantal content (Zhu et al., 2013), but it is also possible that Rab3A is contributing to the mEPSC amplitude changes in mSCS. In sum, these studies suggest miR-124 could tie together presynaptic changes (Rab3A, transmitter released from a single vesicle) and postsynaptic changes (synaptopodin, AMPA receptors) in divergent homeostatic synaptic plasticity.

A model for the role of Rab3A in homeostatic plasticity of mEPSC amplitude

We propose a model for Rab3A function in homeostatic synaptic plasticity of mEPSC amplitude in which Rab3A contributes to presynaptic quantal size through the amount of transmitter released by fusion of a single vesicle, and to postsynaptic quantal size through the number of postsynaptic receptors. Normally, Rab3A cycles between GTP and GDP-bound forms (“State 1,” indicated by red tag in Figure 10A) and promotes synaptic vesicle mobilization and fusion, at the same time maintaining vesicle size/transmitter content within a narrow range. In this normally cycling Rab3A state, an anterograde signal dependent on Rab3A maintains GluA2 receptors at their normal, restricted level. Activity blockade leads to a buildup of an alternate form of Rab3A, possibly one that is unable to cycle through its GTP and GDP-bound forms (“State 2,” indicated by gray tag in Figure 10B). Without the quality control mechanism conferred by State 1 Rab3A, vesicle size/transmitter content skews to higher values, and the anterograde signal is abnormal, leading to GluA2 receptors levels rising outside normal limits. Not shown here, the Earlybird point mutation may permanently be in “State 2,” as the mutation is in the guanine nucleotide binding region (Kapfhamer), which would lead to increased transmitter released from single vesicles and increased receptor levels.

Rab3A is required for increased mEPSC amplitude after activity blockade via 2 mechanisms: 1. increased amount of transmitter released from a single vesicle and 2. increased postsynaptic receptors. Left: in the presence of normal activity, Rab3A is in its control state (“State 1,” red), maintaining mEPSC amplitude within a limited range by regulating either vesicle size or transmitter content, or both, and by regulating the number of postsynaptic receptors, likely through release of an anterograde signal. Right: after activity is blocked, Rab3A accumulates in an alternate functional state (“State 2,” gray) resulting in quantal size expanding to larger ranges because of increased vesicle size or transmitter content and because of increased numbers of receptors no longer being under the normal control of the anterograde signal.

We have demonstrated an essential role of the presynaptic vesicle protein Rab3A in the homeostatic increase in mEPSC amplitude and GluA2 AMPA receptor levels at synaptic sites in mouse dissociated cortical cultures. To our knowledge, this is the first evidence of a presynaptic protein being implicated in the homeostatic regulation of mEPSC amplitude and AMPARs in neuronal cultures. The current results extend our previous findings that showed the homeostatic increase in mEPC amplitudes at the mouse NMJ in vivo depends on Rab3A. Under conditions where postsynaptic receptors are close to being saturated by the amount of transmitter in a single vesicle (possibly in cortical cultures but not at the NMJ), it may be necessary to increase both the amount of transmitter released and the number or function of postsynaptic receptors in order to observe a homeostatic increase in the physiological signal, mEPSC amplitude. Our results suggest that in cortical cultures, neuronal Rab3A is a key player regulating the homeostatic increase in synaptic strength on both sides of the synapse.