Abstract
Following prolonged activity blockade, amplitudes of miniature excitatory postsynaptic currents (mEPSCs) increase, a form of plasticity termed “homeostatic synaptic plasticity.” We previously showed that a presynaptic protein, the small GTPase Rab3A, is required for full expression of the increase in miniature endplate current amplitudes following prolonged blockade of action potential activity at the mouse neuromuscular junction in vivo (Wang et al., 2011), but it is unknown whether this form of Rab3A-dependent homeostatic plasticity shares any characteristics with central synapses. We show here that homeostatic synaptic plasticity of mEPSCs is impaired in mouse cortical neuron cultures prepared from Rab3A−/− and mutant mice expressing a single point mutation of Rab3A, Rab3A Earlybird mice. To determine if Rab3A is involved in the well-established homeostatic increase in postsynaptic AMPA-type receptors (AMPARs), we performed a series of experiments in which electrophysiological recordings of mEPSCs and confocal imaging of synaptic AMPAR immunofluorescence were assessed within the same cultures. We found that Rab3A was required for the increase in synaptic AMPARs following prolonged activity blockade, but the increase in mEPSC amplitudes was not always accompanied by an increase in postsynaptic AMPAR levels, suggesting other factors may contribute. Finally, we demonstrate that Rab3A is acting in neurons because only selective loss of Rab3A in neurons, not glia, disrupted the homeostatic increase in mEPSC amplitudes. This is the first demonstration that neuronal Rab3A is required for homeostatic synaptic plasticity and that it does so partially through regulation of the surface expression of AMPA receptors.
Introduction
One of the most studied phenomena triggered by prolonged activity blockade is the increase in amplitudes of miniature excitatory postsynaptic currents (mEPSCs) in neurons. First demonstrated in cultures of dissociated cortical neurons (Turrigiano et al., 1998) and spinal cord neurons (O’Brien et al., 1998), the compensatory response to prolonged loss of activity has been dubbed “homeostatic synaptic plasticity” (Turrigiano and Nelson, 2004; Pozo and Goda, 2010). These first two studies provided evidence that glutamate receptors were also increased after inactivity, via an increased response to exogenously applied glutamate (Turrigiano et al., 1998) and increased immunofluorescent labeling of GluA1 receptors at synaptic sites (O’Brien et al., 1998). The accompanying increase in synaptic AMPA receptors has been confirmed with immunofluorescence in multiple studies of hippocampal and cortical cultures treated with TTX (Wierenga et al., 2005; Hou et al., 2008; Ibata et al., 2008; Jakawich et al., 2010a; Tatavarty et al., 2013; Xu and Pozzo-Miller, 2017; Dubes et al., 2022). It is now well accepted that the homeostatic increase in mEPSC amplitudes in neurons is due to an increase in postsynaptic AMPA receptors.
In our previous work, we found that a TTX cuff applied for 48 hr around the sciatic nerve in mice led to an increase in the amplitude of the miniature endplate currents (mEPCs) recorded in tibialis anterior muscles. Surprisingly, and unlike central synapses, we could find no evidence that the acetylcholine receptors levels at the neuromuscular junction (NMJ) were increased in the TTX-cuff treated mice (Wang et al., 2005). This result led us to search for presynaptic molecules that might homeostatically regulate the presynaptic quantum. In previous studies in chromaffin cells, we identified the small GTPase Rab3A, a presynaptic 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 at the NMJ. The results were clear: the homeostatic increase in mEPC amplitude induced by the TTX cuff was strongly reduced in the Rab3A−/− mouse, and was completely abolished in the Earlybird mutant (Rab3AEbd/Ebd), which has a single point mutation in Rab3A that causes a shift towards early awakening that is more dramatic than that of the Rab3A deletion mouse (Kapfhamer et al., 2002; Wang et al., 2011). These results led us to conclude that the homeostatic increase in mEPC amplitude at the NMJ is a presynaptic phenomenon.
It has remained somewhat of a mystery what the role of Rab3A is in synaptic transmission. 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). There are some modest changes in short term plasticity during repetitive stimulus trains, with increased depression observed in response to moderate frequencies in hippocampal slices (Geppert et al., 1994) and at the mammalian NMJ (Coleman and Bykhovskaia, 2009), but increased facilitation in response to high frequency stimulation at the mammalian NMJ (Coleman and Bykhovskaia, 2010). The most dramatic effect of loss of Rab3A is the abolishment of a presynaptic form of long-term potentiation (LTP) at the mossy fiber-CA3 synapse that is induced by 25 Hz stimulation in the presence of NMDA blockers (Weisskopf et al., 1994; Castillo et al., 1997). Our result that prolonged inactivity-induced plasticity of mEPC amplitude at the mammalian NMJ is disrupted with loss of function of Rab3A further cements a role for Rab3A in activity-dependent plasticity of synaptic strength. However, a crucial question that remains is whether the TTX-cuff induced effect on the amplitude of the spontaneous synaptic current at the NMJ shares underlying mechanisms with central synapses. We set out to answer this question for the more well-studied phenomenon of homeostatic plasticity of mEPSCs in dissociated mouse cortical neurons prepared from wild type, Rab3A−/− mice and Rab3aEbd/Ebd mice treated with TTX for 48 hr. We show here that homeostatic synaptic plasticity of mEPSC amplitude was strongly reduced in cultures from Rab3A−/− mice, and completely abolished in cultures from Rab3AEbd/Ebd mice, supporting a shared mechanism with the neuromuscular junction. Unexpectedly, the modest increase in GluA2 immunofluorescence we observed in wild type cultures was not present in Rab3A−/− cultures. Loss of Rab3A in neurons, but not in glia, disrupted the increase in mEPSC amplitudes after TTX treatment. These results support the idea that the increase in mEPSC amplitude is due to more postsynaptic receptors, and that neuronal Rab3A is necessary for this regulation. However, we found that GluA2 receptors and mEPSC amplitudes did not always change in parallel. This mismatch suggests there may be contributors other than postsynaptic receptor numbers to the homeostatic plasticity of mEPSC amplitude.
Results
We previously reported that mixed cultures of cortical neurons and glia prepared from postnatal day 0-2 mouse pups responded to a block of action potential-mediated activity by a 48 hr 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. The Rab3A Earlybird mutation was discovered in a mutagenesis screen in mice for genes involved in the rest-activity cycle (Kapfhamer et al., 2002). The Earlybird mutant was named for its phenotype, a significantly shorter circadian period which resulted in early awakening. In that study, the authors concluded that the Earlybird mutation was likely dominant-negative, as the shift in circadian period was substantially more dramatic than that observed in the Rab3A−/− mice. As described in the Introduction, the effects of loss of Rab3A have been modest, so we included experiments using the Rab3A Earlybird mutant mice to see if there would be a more robust phenotype, and in doing so, strengthen conclusions based on the Rab3A deletion mice. To obtain Rab3A−/− and Rab3AEbd/Ebd homozygotes, 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.
We found that loss of Rab3A resulted in a greatly impaired homeostatic increase after activity blockade. 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 hr prepared from wild type and Rab3A−/− mice in the Rab3A+/− colony are shown in Figures 1A and B, respectively. Average mEPSC waveforms from the same recordings are shown in Figures 1C and D. The mean mEPSC amplitudes for 30 control and 23 TTX-treated neurons from Rab3A+/+ cultures are displayed in the box and whisker plot in Figure 1E; after activity blockade the average mEPSC amplitude increased from 13.9 ± 0.7 pA to 18.2 ± 0.9 pA. In contrast, in cultures prepared from Rab3A−/− mice, the average mEPSC amplitude was not significantly increased, for 25 untreated cells and 26 TTX-treated cells (Figure 1F, 13.6 ± 0.1 pA vs. 14.3 ± 0.6 pA). We found that TTX treatment also resulted in an increase in mEPSC frequency in cultures prepared from Rab3A+/+ mice, as shown in Figure 1G (CON, 2.26 ± 0.37 sec−1; TTX, 4.62 ± 0.74 sec−1). This TTX-induced frequency increase was strongly reduced in cultures from Rab3A−/− mice, but there were a few outliers with very high frequencies still present after TTX treatment (Figure 1H, CON, 2.74 ± 0.49 Hz; TTX, 3.23 ± 0.93 Hz).
We next examined homeostatic plasticity in cultures prepared from wild type mice in the Rab3A+/Ebd colony. Example current traces of mEPSCs, and average mEPSC waveforms, are shown in Figures 2A and C, respectively, for cortical cultures prepared from wild type mice in the Rab3A+/Ebd colony. Treatment with TTX for 48 hr led to a significant increase in the average mEPSC amplitude of 23 TTX-treated cells when compared to 20 untreated cells (Figure 2E, CON, 11.0 ± 0.6 pA; TTX, 15.0 ± 1.3 pA). We note here that while the two wild type 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).
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 2B and D, 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). There was a trend towards increased mEPSC frequency in TTX-treated cultures from this wild type strain, but the difference was not as robust as for the wild type cultures from Rab3A+/− colony, and did not reach statistical significance (Figure 2G, CON, 1.15 ± 0.19 sec−1; TTX, 2.54 ± 0.55 sec−1). This trend was not observed in neurons from cultures prepared from Rab3AEbd/Ebd mice, likely due to an increase in frequency in untreated cells—frequencies remained high in TTX-treated cells (Figure 2H, CON, 1.71 ± 0.41 sec−1; TTX, 3.05 ± 0.80 sec−1).
Our results show that the homeostatic increase in 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 this 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/Ebddata set, mEPSC amplitudes from untreated cultures were significantly larger than those of untreated cultures from wild type mice, as can be seen comparing the CON values in Figures 2E and F (CON, Rab3A+/+, 11.0 ± 0.7 pA, vs. CON, Rab3AEbd/Ebd, 15.1 ± 1.0 pA, p = 0.0027, Kruskal-Wallis test). The increase in mEPSC amplitude in cultures from Rab3AEbd/Ebd mice is consistent with the increase in mEPC amplitude we observed at the Rab3AEbd/EbdNMJ (Wang et al., 2011). Similar increases in mEPSC amplitude at baseline, combined with a failure to increase further after activity blockade, were noted in cultures prepared from mECP2, AKAP, Homer1a, and Arc/Arg3.1 deletion mice (Shepherd et al., 2006; Hu et al., 2010; Xu and Pozzo-Miller, 2017; Sanderson et al., 2018) as well as in cultures treated with ubiquitin proteasome inhibitors (Jakawich et al., 2010a) or microRNA 186-5p (Silva et al., 2019), suggesting the possibility of a generalized phenomenon in which increased mEPSC amplitude is induced when homeostatic regulation is disrupted. Since we did not observe such an increase in mEPSC amplitude at baseline in cultures from Rab3A−/− mice, it remains a possibility that the point mutant may bind to novel partners, causing activities that would not be either facilitated or inhibited by Rab3A. Still, it is a strong coincidence that the novel activity of the mutant Earlybird Rab3A would affect mEPSC amplitude, the same characteristic that is modulated by activity blockade in a Rab3A dependent manner. Finally, although we cannot rule out that once mEPSC amplitude is increased at baseline, it is not possible for homeostatic plasticity to increase mEPSC amplitude further due to a non-specific ceiling effect, the simplest interpretation is that the presence of the mutant protein mimics the condition of activity blockade.
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)? 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 Ca2+-permeable receptor-specific inhibitor NASPM (Sutton et al., 2006), indicating the entirety of the homeostatic effect is due to an increase in homomeric Ca2+-permeable AMPA receptors, likely GluA1, although GluA3 and GluA4 are also Ca2+-permeable (Hollmann et al., 1991; Burnashev et al., 1992). Heteromeric receptors of any of these subunits in combination with GluA2 would have rendered the receptors Ca2+-impermeable and no longer sensitive to NASPM (Blaschke et al., 1993; Bochet et al., 1994; Jonas and Burnashev, 1995). In contrast, the increase in mEPSC amplitude in mouse hippocampal slice cultures induced by TTX alone was not affected by another Ca2+-permeable AMPA receptor inhibitor, philanthotoxin, suggesting mediation by either GluA2-containing AMPA 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 3A 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; after NASPM, CON 11.9 ± 2.6 pA; TTX 16.1 ± 3.5 pA). This result indicates that the TTX-induced increase in mEPSC amplitude does not depend on Ca2+-permeable receptors, since the effect of TTX clearly remained when their presence was removed by acute NASPM application. We do not think that there was any technical issue with the NASPM application, because overall, mEPSC amplitudes were reduced modestly in both untreated and TTX-treated cultures (Figure 3B, CON, before NASPM, 12.9 ± 3.5 pA; after, 11.9 ± 2.6 pA; TTX, before NASPM, 17.5 ± 3.1 pA; after, 16.1 ± 3.5 pA). Furthermore, NASPM consistently reduced mEPSC frequency (Figure 3C, CON, before NASPM, 1.84 ± 0.55 sec−1; after NASPM, 1.56 ± 0.53 sec−1; TTX, before NASPM, 4.40 ± 3.51 sec−1; after NASPM, 2.68 ± 2.25 sec−1). Similar to the data in Figure 1, frequency is significantly increased after TTX treatment (p = 0.04), but this difference is no longer significant after NASPM application (p = 0.20). A loss in frequency after acute application of NASPM suggests there are synaptic sites that express only homomers of Ca2+-permeable receptors; mEPSCs from these sites would be expected to be completely blocked by NASPM (see cartoon description in Figure 3D, Left). The magnitude of the frequency reduction following acute NASPM appears to be greater after TTX treatment, suggesting that loss of activity promotes the establishment of synaptic sites that contain only Ca2+-permeable receptors. However, these new sites do not appear to contribute to the increase in mEPSC amplitude (Figure 3A). The lack of effect of NASPM also rules out the contribution of Ca2+-permeable forms of Kainate receptors GluA5 and GluA6 (Koike et al., 1997; Sun et al., 2009).
Having established that homomeric Ca2+-permeable receptors are not contributing to the homeostatic increase in mEPSC amplitude, we turned to immunohistochemistry and confocal imaging to assess whether GluA2 receptor expression, which will identify either GluA2 homomers or GluA1-GluA2 heteromers, was increased in our wild type mouse cortical cultures following 48 hr treatment with TTX. Since mEPSCs necessarily report synaptic levels of receptors, we used VGLUT1-immunofluorescence to identify synapses on pyramidal primary apical dendrites labeled with MAP-2 immunofluorescence. We focused on the primary dendrite of pyramidal neurons as a way to reduce variability that might arise from being at widely ranging distances from the cell body, or, from inadvertently sampling dendritic regions arising from other neurons, possibly inhibitory neurons. Figure 4 (Top) shows two pyramidal neurons, one from an untreated culture on the left, and one from a TTX-treated culture on the right, prepared from Rab3A+/+ mice. To measure characteristics of synaptically located GluA2 receptor clusters, we zoomed in on the primary dendrite (region in white rectangle). The zoomed regions for single confocal sections of the selected area are shown below, and contain pairs of VGLUT1– and GluA2-immunofluorescent clusters indicated by white trapezoidal frames. 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-immunofluorescent clusters not associated with VGLUT1 immunofluorescence (a few are indicated with white arrows), most probably extrasynaptic receptors. There were also GluA2-positive clusters present outside of MAP2-positive dendrites, which may be located on astrocytes (Fan et al., 1999). Finally, GluA2 immunofluorescent clusters close to VGLUT1 immunofluorescence but not located along any apparent MAP2-positive neurites suggests the presence of axon-axonal contacts, although VGLUT1 has also been detected in astrocytes (Ormel et al., 2012). Only sites that contained both VGLUT1 and GluA2 immunofluorescence close to the primary MAP2-postive dendrite or a secondary branch along this primary dendrite were selected for analyses. In the experiment from which these images were collected, we analyzed the distance from the cell body for each synapse. The average distance from the cell body, for dendrites from the untreated cell, was 38.5 ± 2.8 μm; for the TTX-treated cell, it was 42.4 ± 3.2 μm (p = 0.35, Kruskal-Wallis test). Since the greatest distance from the cell body is set by the outer limit of the zoom window, and that window was placed adjacent to the cell body, distances in other experiments would be in this same range.
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 immunofluorescence. We completed this matched paradigm of physiology and immunofluorescence, which to our knowledge has never before been done, 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 5, with the data summarized in Table 1. Levels of GluA2 immunofluorescence at synaptic sites were quantified by the size of the GluA2-immunofluorescent receptor cluster, the average intensity value, and the integral of intensity values across the area of the cluster. 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 5A; Table 1). This result indicates that the homeostatic response averaged across the 3 cultures was very similar to, but slightly smaller than, that of the data set presented in Figure 1. The means for size, intensity and integral of GluA2 immunofluorescent receptor clusters showed trends to higher values after activity blockade, but the differences did not reach statistical significance (Figure 5B size; Figure 5C intensity; Figure 5D integral; Table 1). These results support previous studies demonstrating an increase in AMPA receptors following activity blockade, with the magnitude of the changes in GluA2 receptor cluster size and integral matching that of mEPSC amplitude, at around 20% (Table 1), but not achieving the same level of statistical significance.
As expected, mean mEPSC amplitude was not increased following activity blockade in the data pooled from the new set of 3 Rab3A−/− cultures (Figure 5E). For immunofluorescence results from the same Rab3A−/− cultures, no trend towards higher values was apparent for GluA2 receptor cluster characteristics (Figure 5F size; Figure 5G intensity; Figure 5H integral; Table 1). Taken together, these results indicate that the modest increases in GluA2 receptor cluster size and integral following activity blockade observed in wild-type cultures do not occur in the absence of Rab3A. These results suggest the surprising finding that the presynaptic vesicle protein Rab3A is involved in the regulation of postsynaptic AMPA receptors during homeostatic plasticity.
As noted above, the homeostatic increase in mEPSC amplitude in cultures from Rab3A+/+ mice reached significance whereas that of receptor levels did not. To determine where this disparity was coming from, we compared mEPSC amplitude and GluA2 receptor data within individual experiments. 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 usable recordings in one day (6 CON and 6 TTX; i.e. low holding current, stable baseline, low noise, etc.). Ultimately, our sample size ranged from 6 to 10 cells per condition per culture. Figures 6A-C show the box plots for mEPSC amplitudes, and GluA2 receptor cluster size, intensity and integrals, for 3 individual Rab3A+/+ cultures. With these small sample sizes, the differences between untreated and TTX-treated cells were not statistically significant for mEPSC amplitudes in any of the 3 cultures, but an upward trend occurred in every case, ranging from 12 to 28% (Table 1). Second, in two of the three Rab3A+/+ experiments (Cultures #1 and #2) GluA2 receptor cluster size increased 22 and 56%, respectively, significantly so for Culture #2 (Table 1). However, in one experiment (Culture #3), GluA2 receptor cluster size showed a downward trend of 9.5%. Because we have 10 CON mean mEPSC amplitudes and 9 TTX mean mEPSC amplitudes that indicate this particular culture very likely has a normal homeostatic increase in mEPSC amplitude induced by activity blockade, the lack of increase in GluA2 receptor cluster characteristics suggests that mEPSC amplitudes can increase in the absence of an increase in receptors. Furthermore, Culture #3’s opposite results from Cultures #1 and #2, explain why in the pooled data set, the increase in GluA2 receptor characteristics after TTX treatment do not reach statistical significance.
Figures 7A-C show the box plots for mEPSC amplitudes, and GluA2 receptor cluster size, intensity and integrals, for 3 individual Rab3A−/− cultures. Here, mEPSC amplitudes showed a downward trend in all 3 cultures, as did GluA2 receptor characteristics for two of the three cultures. But we found another obvious mismatch between the mEPSC amplitude data and the GluA2 receptor cluster data for one of the Rab3A−/− cultures. For Culture #3 (Figure 7C), the box plot of mEPSC mean amplitudes showed a downward trend after TTX treatment, decreasing by 7% (Table 1), whereas GluA2 receptor cluster size, intensity and integral all showed upward trends, the increase ranging from 6 to 23% (Table 1). This result may suggest that the disruption of the homeostatic increase in receptors in the absence of Rab3A is not absolute. It also indicates that synaptic GluA2 receptor levels can increase without an accompanying increase in mEPSC amplitudes.
In summary, our matched mEPSC and receptor cluster results indicate that: 1. activity blockade resulted in an increase in mEPSC amplitudes usually accompanied by an increase in synaptic GluA2 receptors; 2. loss of Rab3A disrupted both the increase in mEPSC amplitudes and the increase in GluA2 receptor levels; and 3. The differences between homeostatic responses of mEPSC amplitudes and the GluA2 cluster sizes within the same cultures suggest there may be other factors that contribute to mEPSC amplitude.
Having used VGLUT1 immunofluorescence to mark synaptic sites, we could also ask whether there was a presynaptic effect of activity blockade on the size, intensity or integral of the VGLUT1 signal at sites apposed to the previously examined GluA2 postsynaptic sites. In the pooled data from cultures prepared from Rab3A+/+ mice, we found no change in the size of the VGLUT1-positive regions, and trends towards reduced intensity and integral (Table 2). These data do not support the idea that activity blockade caused an increase in the amount of VGLUT1 per vesicle, although we cannot rule out that a simultaneous reduction in the number of vesicles obscured the effect, since we did not have an independent label for synaptic vesicles. In cultures prepared from Rab3A−/− mice, the size, intensity and integral of the VGLUT1-positive regions showed upward trends (Table 2), suggesting the possibility that the modest decreases observed for cultures prepared from wild type mice depend on Rab3A, but it is difficult to make any conclusions when the p values are far from reaching significance. Similarly to the GluA2 receptor cluster data, the likely reason that the pooled data effects did not reach significance is because the effects were not consistent across all three cultures, with one Rab3A+/+ culture showing an upward trend instead of a downward trend, and one Rab3A−/− culture showing a downward trend instead of an upward trend (Table 2). It should be noted that there was no consistent inverse pattern for GluA2 vs. VGLUT1. While in cultures prepared from Rab3A−/− mice, when GluA2 characteristics trended downward, VGLUT1 characteristics trended upward, this pattern did not hold for cultures prepared from Rab3A+/+ mice. In Rab3A+/+ Culture #2, GluA2 showed significant increases, and VGLUT1 showed essentially no change. In Rab3A+/+ Culture #3, both GluA2 and VGLUT1 trended downwards. Taken together, the GluA2 and VGLUT1 immunofluorescence data indicate that GluA2 receptors likely contribute to the increase in mEPSC amplitudes after activity blockade in Rab3A+/+ cultures, and the disruption of that increase in Rab3A−/− cultures, but VGLUT1 levels change in the opposite direction from mEPSC amplitudes. Interestingly, VGLUT1 levels are strongly affected by loss of Rab3A, with a highly significant decrease in integral in untreated cultures from Rab3A−/− mice, compared to untreated cultures from Rab3A+/+ mice (size, p = 0.019; intensity, p = 0.004; integral, p = 4.6 * 10−4; Table 2).
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 glial 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 glial cells (Stellwagen and Malenka, 2006). It has recently been shown that the source of TNFα is astrocytes rather than microglia (Heir et al., 2024). 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 8, Left). It should be noted that under our culturing conditions, the glial feeder layers would be expected to be composed mainly of astrocytes (Heir et al., 2024). We compared the effect of activity blockade on mEPSC amplitudes recorded from cortical neurons from Rab3A+/+ mice plated on Rab3A+/+ glial feeder layers (Figure 8A); neurons from Rab3A+/+ mice plated on Rab3A−/− glia (Figure 8B), and neurons from Rab3A−/− mice plated on Rab3A+/+ glia (Figure 8C). If Rab3A is required for TNFα release from glia, then Rab3A+/+ neurons plated on Rab3A−/− glia should not show a homeostatic increase in mEPSC amplitude after treatment with TTX, and any cultures with Rab3A+/+ glia should have a normal homeostatic response. We found the opposite result: mEPSC amplitudes increased dramatically in cultures where Rab3A was present in neurons (WT on WT, CON, 13.3 ± 0.5 pA, TTX, 16.7 ± 1.2 pA; WT on KO, CON, 13.3 ± 1.0 pA, TTX, 18.8 ± 1.4 pA), but the increase was greatly diminished in cultures where Rab3A was present only in glia (KO on WT, CON, 15.2 ± 1.1 pA, TTX, 16.9 ± 0.7 pA). It is interesting to note that in the glial feeder layer culture, the loss of Rab3A in the neurons may have caused an increase in mEPSC amplitude at baseline, similar to the Earlybird mutant in dissociated cultures, but it is also possible this can be attributed to culture to culture variability. Regarding mEPSC frequency, in cultures with Rab3A+/+ neurons, frequency trended upward after activity blockade, whereas in the culture with Rab3A−/− neurons, frequency was already increased, and trended downward after activity blockade.
In summary, in the absence of Rab3A from glia, the homeostatic plasticity of mEPSC amplitude following activity blockade was completely normal, whereas in the absence of Rab3A in neurons, homeostatic plasticity was greatly diminished. This result makes it highly unlikely that Rab3A is required for the release of TNFα, or another factor from glia, 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 glia behaved remarkably similar to the mouse NMJ in response to loss of Rab3A function: in cultures from Rab3A−/− mice, the increase in mEPSC amplitudes following prolonged network silencing by TTX was strongly diminished, and in cultures from Rab3AEbd/Ebd mice, basal mEPSC amplitude was increased compared to that of wild-type cultures and was not further modified by TTX treatment. These results suggest that normal function of the synaptic vesicle protein Rab3A is required for the homeostatic increase of mEPSC amplitude in cortical cultures. The effect of TTX on mEPSC amplitude was not due to increases in Ca2+-permeable homomers because acute NASPM application was unable to reverse the increase in amplitude after TTX treatment. We used imaging of immunofluorescence to demonstrate a strong trend towards increased GluA2 levels at synapses in activity-blocked cultures from Rab3A+/+ mice that was no longer apparent in the absence of Rab3A, suggesting the surprising idea that this protein, previously implicated primarily in the regulation of calcium-dependent vesicle fusion, is also important for the regulation of surface expression of postsynaptic receptors. We found a significant increase in mEPSC frequency after activity blockade that was no longer significant in the absence of Rab3A, which was less surprising as frequency is known to depend on presynaptic factors. At the neuromuscular junction, higher mEPC frequencies that occur in 0.5 mM extracellular calcium and above were significantly reduced in the Rab3A deletion mouse (Coleman et al., 2007). On the other hand, the increase in frequency after TTX treatment is also strongly diminished after acute application of NASPM, suggesting there may be an increase in the number of synapses with only Ca2+-permeable homomers. It remains to be determined if this effect is Rab3A-dependent or if there are two distinct mechanisms contributing to homeostatic plasticity of frequency.
The levels of the glutamate transporter, VGLUT1, were not increased after activity blockade under our experimental conditions. A similar lack of increase in VGLUT1 levels, relative to synapsin in imaging of immunofluorescence, was previously observed in hippocampal cultures treated with NBQX, an AMPA receptor blocker (see supplemental data, (Wilson et al., 2005)). On the other hand, treatment with AP5, an NMDA receptor blocker, was able to induce an increase in VGLUT1 levels relative to synapsin (Wilson et al., 2005), suggesting that this presynaptic response depends on the type of activity manipulation. Another study showed increased VGLUT1 mRNA and protein in western blots in cortical cultures treated with TTX, but only examined VGLUT1 at synapses using imaging following treatment with bicuculline (De Gois et al., 2005). Taken together with our findings, it appears that there is not strong evidence that the amount of VGLUT1 at synapses is a contributing factor to the increase in mEPSC amplitude after activity blockade with TTX.
Homeostatic plasticity of postsynaptic receptors was less consistent than that of mEPSC amplitude
While we found that the effect of TTX treatment on mEPSC amplitude and GluA2 receptor levels were of very similar magnitude, the effect on GluA2 receptor levels was not statistically significant for the pooled data sets acquired in the same cultures. We determined that this discrepancy was due to one of the three cultures showing a downward trend in GluA2 receptor levels after activity blockade. In light of the extensive literature documenting increases in GluA1 or GluA2 receptors, or both, after activity blockade, the less than convincing increase in GluA2 receptor levels in the current work came as a surprise. Therefore we reviewed what has previously been shown, and were able to identify only 5 studies that examined the effect of activity blockade using TTX alone, reported levels of AMPA receptors specifically localized to synapses, and recorded mEPSCs in the same study, although not in the same cultures (Hou et al., 2008; Ibata et al., 2008; Jakawich et al., 2010a; Xu and Pozzo-Miller, 2017; Dubes et al., 2022). Since our issue is not with the magnitude of the changes being different for receptors and mEPSC amplitudes, but with statistical significance (or consistency) being different between the two methods, we wanted to determine whether the previous studies showed any evidence of greater variability in receptor measurements. Only a single study had comparable sample sizes for receptors and mEPSC amplitudes (either because number of receptor sites is used instead of cells, or because the number of cells differs for the two types of experiments), and the p value was < 0.01 for mEPSC amplitudes (n = 9 and 10 cells for CON and TTX, respectively), but < 0.05 for receptor data (n = 13 and 14 cells for CON and TTX) (Xu and Pozzo-Miller, 2017). This single example is consistent with the idea that the receptor data are more variable than mEPSC data. Even if a difference in consistency could be established in multiple studies, it remains a possibility that the disparity is attributable to the unique limitations of the two types of experiments. Nevertheless, the finding that mEPSC amplitude can increase despite a decrease in receptors suggests receptors may not be the only contributor to the TTX effect on mEPSC amplitude.
Neuronal, not glial, Rab3A may be required for the homeostatic release of a signaling molecule
Our results support the surprising possibility that the presynaptic vesicle protein Rab3A regulates postsynaptic AMPARs. Astrocytic release of TNFα was shown to mediate the increase in mEPSC amplitudes in prolonged TTX-treated cultures of dissociated hippocampal neurons and cultures of hippocampal slices (Stellwagen and Malenka, 2006; Heir et al., 2024), but we found that loss of Rab3A in glia did not disrupt the increase in mEPSC amplitude 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). This leaves neuronal Rab3A as the potential controller of a soluble factor regulating surface expression of receptors.
One possible signaling molecule being controlled by Rab3A is brain-derived neurotrophic factor (BDNF). 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). Reduction of Rab3A in astrocytes decreases astrocytic BDNF release (Hong et al., 2016). Taken together, these studies suggest that activity blockade could cause a Rab3A-dependent inhibition of BDNF release, leading to an increase in postsynaptic receptors. However, we have no evidence that Rab3A is required for the regulation of BDNF release from neurons, so it is just as likely that after activity blockade, Rab3A acts in a positive way to promote release of a trophic molecule that stimulates the insertion of postsynaptic receptors. It is unlikely that Rab3A is required when activity is normal because receptor levels do not change in the same way with loss of Rab3A as they do with loss of activity. Rather, our data fit with a model in which Rab3A mainly becomes engaged in response to activity blockade (Figure 9).
Is Rab3A acting in the postsynaptic cell to cause homeostatic increase in mEPSC amplitude?
To our knowledge, there is no evidence localizing Rab3A to postsynaptic sites. Other 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 Rab3A localization was not examined in those studies. On the other hand, a recent report implicates Rab3A in the localization of plasma membrane proteins in rafts in HEK cells and Jurkat-T cells, indicating Rab3A can be found in locations other than the presynaptic nerve terminal (Diaz-Rohrer et al., 2023). In this study a reduction in all Rab3 isoforms led to an impairment of plasma membrane localization of epidermal growth factor receptor (EGFR), so it is possible that Rab3A is acting similarly in postsynaptic dendrites to promote transport of GluA2 receptors to the surface after activity blockade (Figure 9). Another possibility, not depicted in Figure 9, is that Rab3A acts postsynaptically to regulate the release of a trophic factor that in turn regulates surface expression of GluA2 receptors. This idea is consistent with evidence that BDNF is synthesized and released locally from dendrites (Jakawich et al., 2010b; Thapliyal et al., 2022).
Alternate contributors to homeostatic increase in mEPSC amplitude
We do not think Rab3A impacts homeostatic synaptic plasticity solely through pre-or 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 any increase in AChR levels, yet it was dependent on Rab3A (Wang et al., 2011). Second, mEPSC amplitudes increased in three cortical cultures treated with TTX, but GluA2 receptor levels in one of the cultures actually went down. An increase in conductance could explain the increase in mEPSC amplitudes in the absence of receptor level increases, but all evidence to date suggests that the only way to increase conductance is to either switch from GluA2-containing heteromeric receptors to GluA1 homomers, which have much greater conductance (Oh and Derkach, 2005; Soares et al., 2013; Hou et al., 2015; Xu and Pozzo-Miller, 2017; Benke and Traynelis, 2019; Silva et al., 2019; Dubes et al., 2022), or, to post-translationally modify the receptor. Such conductance-altering modifications have been demonstrated only for GluA1 (Derkach et al., 1999; Kristensen et al., 2011; Jenkins and Traynelis, 2012). Importantly, any change in conductance due to an increase in GluA1 homomers or modification of GluA1 should have been detected by NASPM being able to disrupt the homeostatic increase in mEPSC amplitude, and this was not what we observed.
The other contributor to mEPSC amplitude besides receptors is the size of the presynaptic quantum, or the amount of transmitter released during a single fusion event. As stated above, we found a downward trend in VGLUT1 levels after activity blockade. This downward trend was not apparent in cultures prepared from Rab3A−/− mice, suggesting a possible Rab3A-dependent effect. Interestingly, Rab3A levels are reduced 50% in hippocampal extracts from mice lacking VGLUT1 (Fremeau et al., 2004). Coupled with our finding that loss of Rab3A was associated with a highly significant decrease in VGLUT1 levels in untreated cultures, these results suggest an interdependence between Rab3A and VGLUT1. More importantly, our results rule out VGLUT1 levels as the source for increased mEPSC amplitude.
The amount of transmitter released during vesicle fusion can be affected by the kinetics of the fusion pore opening (Chang et al., 2017). We previously used amperometric measurements in mouse adrenal chromaffin cells to show loss of Rab3A increased the occurrence of very small amplitude fusion pore feet (Wang et al., 2008). Paired with our finding that there is an increase in the occurrence of very slow-rising, abnormally shaped mEPCs at the neuromuscular junction of Rab3A−/− mice, our data suggest that small synaptic vesicles may have a fusion pore step under certain circumstances (see also (Chiang et al., 2018)). If activity blockade causes an increase in mEPSC amplitude due to a more rapid opening of a fusion pore, or a larger fusion pore conductance, it is possible that Rab3A is required for this modulation.
Finally, a simple way to increase the amount of transmitter released by a vesicle is to enlarge vesicle size. Larger vesicle diameter has been observed following activity disruption in multiple paradigms, including ocular deprivation, cortical ablation, and immobilization (Kawana et al., 1971; Lloret and Saavedra, 1975; Cheresharov et al., 1978). Loss of Rab3 family members is associated with increased vesicle size (Riedel et al., 2002; Schonn et al., 2010). Finally, an increase in volume has been shown to increase miniature endplate junctional current amplitude at the Drosophila NMJ (Karunanithi et al., 2002). Rab3A could be required for an activity-dependent increase in vesicle size, but the increase in diameter needed for a 25% increase in transmitter, and therefore mEPSC amplitude, would be difficult to detect.
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. Homeostatic synaptic plasticity is becoming increasingly implicated in essential normal functions such as sleep (Tononi and Cirelli, 2014; Diering et al., 2017; Torrado Pacheco et al., 2021), and pathological conditions such as 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). Therefore our identification of a novel regulatory pathway involving Rab3A may have important therapeutic implications. In future work it will be important to determine whether Rab3A is acting presynaptically or postsynaptically to regulate GluA2 receptors, and what other contributor to mEPSC amplitude depends on Rab3A, possibly fusion pore characteristics or vesicle volume.
Materials and methods
Animals. To determine the role of Rab3A in homeostatic plasticity in mouse cortical cultures, we employed two distinct genetic mouse strains with altered Rab3A function. Rab3A deletion mice were generated as follows: Rab3A+/− heterozygous mice were bred and genotyped as previously described (Kapfhamer et al., 2002; Wang et al., 2008), and maintained as a colony using multiple heterozygous breeding pairs. Homozygous Rab3A+/+ or Rab3A−/− mouse pups were obtained by breeding Rab3A+/+ pairs and Rab3A−/− pairs, respectively, so that the pups were the first generation progeny of a homozygous mating. This protocol minimizes the tendency of multiple generations of Rab3A+/+ and Rab3A−/− homozygous breedings to produce progeny that are more and more genetically distinct. Rab3AEbd/Ebd mice 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 (forward primer: TGA CTC CTT CAC TCC AGC CT; reverse primer: TGC ACT GCA TTA AAT GAC TCC T) followed by 2. a restriction digest 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 glia 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 glia 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). This plating density produces a “high density” culture characterized by a complex mesh of neuronal processes criss-crossing the field of view, completely filling the space between cell bodies (see Figure 4). We found in preliminary experiments that an increase in mEPSC amplitude in mouse cortical cultures was inconsistent when cell density was low, likely due to lower levels of baseline activity, although it is a limitation of this study that we did not directly measure activity levels in cultures prepared from Rab3A+/+, Rab3A−/− or Rab3AEbd/Ebd mice, and therefore do not know if differences in results are due to differences in activity levels. While unlikely for Rab3A+/+ and Rab3A−/− studies, given the mEPSC amplitude in untreated cultures that was relatively small, it remains possible that in Rab3AEbd/Ebd studies, the larger mEPSC amplitude in untreated cultures was due to a loss of activity in these cultures. 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. This smaller data set was restricted to the time period over which cultures were also prepared from Rab3A−/− mice.
Preparation of Glial Feeder Layers. Glial 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 glial 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 glial 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 glial 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 glia 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 glial proliferation and contamination of the feeder layer with new glia, promoting only neuronal growth on the feeder layers (FdU-containing media was used for the maintenance of these cultures and all subsequent media changes). Glial culture media was removed from the feeder layer cultures, and the neuronal cell suspension was plated onto the glial feeder cultures. The culture strategy used to distinguish the relative roles of neuronal and glial Rab3A is outlined in Figure 8, Left. 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 in pilot studies, all mEPSCs were blocked by 10 µM CNQX and 50 μm picrotoxin, demonstrating no APV-sensitive mEPSCs were present (data not shown). Previous studies have described NMDA mEPSCs in neuronal cultures, but recordings were performed in 0 Mg2+ (Watt et al., 2000; Sutton et al., 2006). In our recordings, the presence of extracellular Mg2+ (1.3 mM) and TTX would mean that the majority of NMDA receptors are blocked by extracellular Mg2+ at resting Vm. Patch electrodes (3.5 – 5 MΩ) were filled with an internal solution containing (in mM): K-gluconate (130), 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, which we have found improves success rate of achieving low access resistance recordings. 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).
NASPM application. The antagonist of Ca2+-permeable AMPA receptors, N-naphthyl acetylspermine (NASPM, 20 μM; Tocris), was applied during recordings in a subset of experiments. NASPM is a synthetic analog of Joro Spider Toxin (JSTX) (Koike et al., 1997), and chemicals in this family block Ca2+ permeable glutamate receptors (Blaschke et al., 1993; Herlitze et al., 1993; Iino et al., 1996). The presence of the edited GluA2 subunit disrupts both the Ca2+ permeability and the sensitivity to JSTX and related substances (Blaschke et al., 1993; Bochet et al., 1994; Jonas and Burnashev, 1995). Therefore, 20 μM NASPM is expected to completely block AMPA receptors containing GluA1, GluA3 or GluA4 subunits (Hollmann et al., 1991; Burnashev et al., 1992), or the Kainate receptors containing GluA5 or GluA6 (Koike et al., 1997; Sun et al., 2009), but be ineffective against any heteromer or homomer containing GluA2 subunits. Because the effect of the spider toxins and their analogues are use-dependent (Herlitze et al., 1993; Iino et al., 1996; Koike et al., 1997), NASPM was applied with a mildly depolarizing high K+ solution (25 mM KCl, 95 mM NaCl) which we expected to trigger release of glutamate and opening of glutamate-activated receptors, allowing entry of the pore blocker. 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). Analysis of mEPSCs. Miniature excitatory postsynaptic currents (mEPSCs) were manually selected using Mini Analysis software (Synaptosoft), a standard program used for mEPSC event detection and analysis. Records were filtered at 2 kHz using a low-pass Butterworth filter prior to selection. The program threshold was set at 3 pA but the smallest mEPSC selected was 4.01 pA. MiniAnalysis selects many false positives with the automated feature when a small threshold amplitude value is employed, due to random fluctuations in noise, so manual re-evaluation of the automated process is necessary to eliminate false positives. If the threshold value is set high, there are few false positives but small amplitude events that visually are clearly mEPSCs are missed, and manual re-evaluation is necessary to add back false negatives or the population ends up biased towards large mEPSC amplitudes. As soon as there is a manual step, bias is introduced. Interestingly, a manual reevaluation step was applied in a recent study that describes their process as ‘unbiased’ (Wu et al., 2020). In sum, we do not believe it is currently possible to perform a completely unbiased detection process. A fully manual detection process means that the same criterion (“does this look like an mEPSC?”) is applied to all events, not just the false positives, or the false negatives, which prevents the bias from being primarily at one end or the other of the range of mEPSC amplitudes. It is important to note that when performing the MiniAnalysis process, the researcher did not know whether a record was from an untreated cell or a TTX-treated cell.
Immunofluorescence, 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 and a single confocal section of the cell acquired (see Top, Figure 4), Fluoview 2.1 software was used to zoom in on the primary dendrite (5X) and a series of 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 identified with MAP2 immunoreactivity. An area of interest (AOI) was manually drawn around the GluA2 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. Multiple thresholds were used, 400-550, but we found that within a given experiment, the magnitude of the TTX effect was not altered by different thresholds, therefore, all data presented here use a threshold of 450. Pixels within the cluster that were above the threshold were automatically outlined, and size, average intensity, and integral of the outlined region saved to Origin2020 (OriginLab) file for calculating mean values. For analyzing VGLUT1 sites, previous AOIs were examined, and if not directly over the VGLUT1 site, moved to encompass it (this happened when GluA2 and VGLUT1 were side by side rather than overlapping). The threshold was 200 for all VGLUT1 data except Rab3A+/+ Culture #1, where threshold was 400. There was no threshold that worked for both that culture and the other cultures, with too much of the background regions being included for Culture #1 with a threshold of 200, and too little of the synaptic site area being included for the other Cultures with a threshold of 400.
Statistics. In our previous publication (Hanes et al., 2020), we described the characteristics of homeostatic plasticity of mEPSC amplitudes in mouse cortical cultures using multiple mEPSC amplitude quantiles from each cell to create data sets with thousands of samples. We compared the cumulative distribution functions for mEPSC amplitudes in untreated vs. TTX-treated cultures, as well as sorting the values from smallest to largest, and plotted TTX values vs. control values in the rank order plot first described in (Turrigiano et al., 1998) to demonstrate synaptic scaling. Finally, we plotted the ratio of TTX values/control values vs. control values to show that the ratio, which represents the effect of TTX, increased with increasing control values, which we termed “divergent scaling.” Here, we wanted to avoid the inflation of sample size caused by pooling multiple mEPSCs per cell in order to statistically compare the effects of TTX in cultures of wild type mice, Rab3A−/− mice and Rab3AEbd/Ebd mutant mice. Therefore, we computed the mean mEPSC amplitude for each cell and compared means from untreated and TTX-treated cultures with the non-parametric Kruskal-Wallis test, with N = the number of cells. The overall means are presented ± SEM in the text, and p values are displayed above each plot. The same process was carried out for the characteristics of synaptic GluA2 receptor clusters and VGLUT1 sites in imaging experiments, with N = the number of dendrites (one dendrite was sampled for each cell). Although we quantified three characteristics of the GluA2 clusters and VGLUT1 sites, size, average intensity, and integral, we have not corrected p values for multiple comparisons. Means are presented as box plots + data, with the box at 25 and 75 percentiles, whiskers at 10 and 90 percentiles, a dot indicating the mean, and a line indicating the median. In Figure 3, mEPSC amplitudes were measured within the same cell before and after NASPM application, so a paired T test was used to compare the pre– and post-NASPM mEPSC amplitudes. Statistics were computed and plots were created using Origin2020 (OriginLab).
Acknowledgements
The authors thank E. Whitlock and H. Ghouse for their contributions to confocal image analyses. This work was partially supported by NINDS P01 NS057228, Tim Cope, PI, KE, PI project 4.
References
- 1.Plasma membrane insertion of the AMPA receptor GluA2 subunit is regulated by NSF binding and Q/R editing of the ion poreProc Natl Acad Sci U S A 107:11080–11085
- 2.Histone Methylation by the Kleefstra Syndrome Protein EHMT1 Mediates Homeostatic Synaptic ScalingNeuron 91:341–355
- 3.AMPA-Type Glutamate Receptor Conductance Changes and Plasticity: Still a Lot of NoiseNeurochem Res 44:539–548
- 4.A single amino acid determines the subunit-specific spider toxin block of alpha-amino-3-hydroxy-5-methylisoxazole-4-propionate/kainate receptor channelsProc Natl Acad Sci U S A 90:6528–6532
- 5.Subunit composition at the single-cell level explains functional properties of a glutamate-gated channelNeuron 12:383–388
- 6.Divalent ion permeability of AMPA receptor channels is dominated by the edited form of a single subunitNeuron 8:189–198
- 7.Rab3A is essential for mossy fibre long-term potentiation in the hippocampusNature 388:590–593
- 8.Light regulates expression of brain-derived neurotrophic factor mRNA in rat visual cortexProc Natl Acad Sci U S A 89:9444–9448
- 9.Fusion pores and their control of neurotransmitter and hormone releaseJ Gen Physiol 149:301–322
- 10.Ultrastructural changes in the anterior horn synapses of rat spinal cord under different locomotor conditionsJ Neural Transm 42:9–21
- 11.The Transmembrane Domain of Synaptobrevin Influences Neurotransmitter Flux through Synaptic Fusion PoresJ Neurosci 38:7179–7191
- 12.Rab3a-mediated vesicle recruitment regulates short-term plasticity at the mouse diaphragm synapseMol Cell Neurosci 41:286–296
- 13.Cooperative regulation of neurotransmitter release by Rab3a and synapsin IIMol Cell Neurosci 44:190–200
- 14.Rab3a deletion reduces vesicle docking and transmitter release at the mouse diaphragm synapseNeuroscience 148:1–6
- 15.Homeostatic scaling of vesicular glutamate and GABA transporter expression in rat neocortical circuitsJ Neurosci 25:7121–7133
- 16.Ca2+/calmodulin-kinase II enhances channel conductance of alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionate type glutamate receptorsProc Natl Acad Sci U S A 96:3269–3274
- 17.Rab3 mediates a pathway for endocytic sorting and plasma membrane recycling of ordered microdomainsProc Natl Acad Sci U S A 120
- 18.Homer1a drives homeostatic scaling-down of excitatory synapses during sleepScience 355:511–515
- 19.miR-124-dependent tagging of synapses by synaptopodin enables input-specific homeostatic plasticityEMBO J 41
- 20.AMPA receptor protein expression and function in astrocytes cultured from hippocampusJ Neurosci Res 57:557–571
- 21.Mechanisms of homeostatic plasticity in the excitatory synapseJ Neurochem 139:973–996
- 22.Vesicular glutamate transporters 1 and 2 target to functionally distinct synaptic release sitesScience 304:1815–1819
- 23.The role of Rab3A in neurotransmitter releaseNature 369:493–497
- 24.Differential vesicular sorting of AMPA and GABAA receptorsProc Natl Acad Sci U S A 113:E922–931
- 25.Divergent Synaptic Scaling of Miniature EPSCs following Activity Blockade in Dissociated Neuronal CulturesJ Neurosci 40:4090–4102
- 26.Astrocytes Are the Source of TNF Mediating Homeostatic Synaptic PlasticityJ Neurosci 44
- 27.Argiotoxin detects molecular differences in AMPA receptor channelsNeuron 10:1131–1140
- 28.Ca2+ permeability of KA-AMPA--gated glutamate receptor channels depends on subunit compositionScience 252:851–853
- 29.Mutant Huntingtin Impairs BDNF Release from Astrocytes by Disrupting Conversion of Rab3a-GTP into Rab3a-GDPJ Neurosci 36:8790–8801
- 30.Homeostatic regulation of AMPA receptor expression at single hippocampal synapsesProc Natl Acad Sci U S A 105:775–780
- 31.MicroRNA miR124 is required for the expression of homeostatic synaptic plasticityNat Commun 6
- 32.Homeostatic scaling requires group I mGluR activation mediated by Homer1aNeuron 68:1128–1142
- 33.Postsynaptic VAMP/Synaptobrevin Facilitates Differential Vesicle Trafficking of GluA1 and GluA2 AMPA Receptor SubunitsPLoS One 10
- 34.Rapid synaptic scaling induced by changes in postsynaptic firingNeuron 57:819–826
- 35.Voltage-dependent blockage of Ca(2+)-permeable AMPA receptors by joro spider toxin in cultured rat hippocampal neuronesJ Physiol 496:431–437
- 36.An essential postsynaptic role for the ubiquitin proteasome system in slow homeostatic synaptic plasticity in cultured hippocampal neuronsNeuroscience 171:1016–1031
- 37.Local presynaptic activity gates homeostatic changes in presynaptic function driven by dendritic BDNF synthesisNeuron 68:1143–1158
- 38.PKC phosphorylates GluA1-Ser831 to enhance AMPA receptor conductanceChannels (Austin 6:60–64
- 39.Molecular mechanisms controlling calcium entry through AMPA-type glutamate receptor channelsNeuron 15:987–990
- 40.LTP requires a unique postsynaptic SNARE fusion machineryNeuron 77:542–558
- 41.Mutations in Rab3a alter circadian period and homeostatic response to sleep loss in the mouseNat Genet 32:290–295
- 42.Quantal size and variation determined by vesicle size in normal and mutant Drosophila glutamatergic synapsesJ Neurosci 22:10267–10276
- 43.Rapid homeostatic plasticity and neuropsychiatric therapeuticsNeuropsychopharmacology 48:54–60
- 44.Enlargement of synaptic vesicles as an early sign of terminal degneration in the rat caudate nucleusJ Comp Neurol 142:297–307
- 45.Syntaxin-4 defines a domain for activity-dependent exocytosis in dendritic spinesCell 141:524–535
- 46.Blocking effect of 1-naphthyl acetyl spermine on Ca(2+)-permeable AMPA receptors in cultured rat hippocampal neuronsNeurosci Res 29:27–36
- 47.Mechanism of Ca2+/calmodulin-dependent kinase II regulation of AMPA receptor gatingNat Neurosci 14:727–735
- 48.Rab3a is involved in transport of synaptic vesicles to the active zone in mouse brain nerve terminalsMol Biol Cell 12:3095–3102
- 49.Homeostatic Plasticity in EpilepsyFront Cell Neurosci 14
- 50.Postsynaptic membrane fusion and long-term potentiationScience 279:399–403
- 51.Enlargement of synaptic vesicles in degenerating nerve endings: a comparison between cat and monkeyJ Neurocytol 4:1–6
- 52.Mechanism of action of rab3A in mossy fiber LTPNeuron 21:1141–1150
- 53.Synaptic plasticity mechanisms common to learning and alcohol use disorderLearn Mem 25:425–434
- 54.A plethora of presynaptic proteins associated with ATP-storing organelles in cultured astrocytesGlia 26:233–244
- 55.Neuronal Activity Patterns Regulate Brain-Derived Neurotrophic Factor Expression in Cortical Cells via Neuronal CircuitsFront Neurosci 15
- 56.NSF binding to GluR2 regulates synaptic transmissionNeuron 21:87–97
- 57.Activity-dependent modulation of synaptic AMPA receptor accumulationNeuron 21:1067–1078
- 58.Dominant role of the GluR2 subunit in regulation of AMPA receptors by CaMKIINat Neurosci 8:853–854
- 59.VGLUT1 is localized in astrocytic processes in several brain regionsGlia 60:229–238
- 60.The AMPA receptor GluR2 C terminus can mediate a reversible, ATP-dependent interaction with NSF and alpha– and beta-SNAPsNeuron 21:99–110
- 61.Unraveling mechanisms of homeostatic synaptic plasticityNeuron 66:337–351
- 62.Rab3D is not required for exocrine exocytosis but for maintenance of normally sized secretory granulesMol Cell Biol 22:6487–6497
- 63.BDNF has opposite effects on the quantal amplitude of pyramidal neuron and interneuron excitatory synapsesNeuron 21:521–530
- 64.Control of Homeostatic Synaptic Plasticity by AKAP-Anchored Kinase and Phosphatase Regulation of Ca(2+)-Permeable AMPA ReceptorsJ Neurosci 38:2863–2876
- 65.Rab3 superprimes synaptic vesicles for release: implications for short-term synaptic plasticityJ Neurosci 26:1239–1246
- 66.Rab3 proteins involved in vesicle biogenesis and priming in embryonic mouse chromaffin cellsTraffic 11:1415–1428
- 67.Arc/Arg3.1 mediates homeostatic synaptic scaling of AMPA receptorsNeuron 52:475–484
- 68.MicroRNA-186-5p controls GluA2 surface expression and synaptic scaling in hippocampal neuronsProc Natl Acad Sci U S A 116:5727–5736
- 69.Impairment and Restoration of Homeostatic Plasticity in Cultured Cortical Neurons From a Mouse Model of Huntington DiseaseFront Cell Neurosci 13
- 70.Differential subcellular targeting of glutamate receptor subtypes during homeostatic synaptic plasticityJ Neurosci 33:13547–13559
- 71.Fragile X protein FMRP is required for homeostatic plasticity and regulation of synaptic strength by retinoic acidJ Neurosci 30:16910–16921
- 72.Interaction of the N-ethylmaleimide-sensitive factor with AMPA receptorsNeuron 21:393–400
- 73.Synaptic scaling mediated by glial TNF-alphaNature 440:1054–1059
- 74.Differential regulation of AMPA receptor and GABA receptor trafficking by tumor necrosis factor-alphaJ Neurosci 25:3219–3228
- 75.A neuronal role for SNAP-23 in postsynaptic glutamate receptor traffickingNat Neurosci 13:338–343
- 76.Calcium-permeable presynaptic kainate receptors involved in excitatory short-term facilitation onto somatostatin interneurons during natural stimulus patternsJ Neurophysiol 101:1043–1055
- 77.Miniature neurotransmission stabilizes synaptic function via tonic suppression of local dendritic protein synthesisCell 125:785–799
- 78.Convergence of distinct signaling pathways on synaptic scaling to trigger rapid antidepressant actionCell Rep 37
- 79.How to scale down postsynaptic strengthJ Neurosci 33:13179–13189
- 80.Retinoic acid-gated BDNF synthesis in neuronal dendrites drives presynaptic homeostatic plasticityElife 11
- 81.Rab3A mediates vesicle delivery at photoreceptor ribbon synapsesJ Neurosci 32:6931–6936
- 82.Sleep and the price of plasticity: from synaptic and cellular homeostasis to memory consolidation and integrationNeuron 81:12–34
- 83.Sleep Promotes Downward Firing Rate HomeostasisNeuron 109:530–544
- 84.Activity deprivation leads to seizures in hippocampal slice cultures: is epilepsy the consequence of homeostatic plasticity?J Clin Neurophysiol 24:154–164
- 85.Homeostatic plasticity in the developing nervous systemNat Rev Neurosci 5:97–107
- 86.Activity-dependent scaling of quantal amplitude in neocortical neuronsNature 391:892–896
- 87.Exosomes and Homeostatic Synaptic Plasticity Are Linked to Each other and to Huntington’s, Parkinson’s, and Other Neurodegenerative Diseases by Database-Enabled Analyses of Comprehensively Curated DatasetsFront Neurosci 11
- 88.Regulation of quantal shape by Rab3A: evidence for a fusion pore-dependent mechanismJ Physiol 586:3949–3962
- 89.Impaired activity-dependent plasticity of quantal amplitude at the neuromuscular junction of Rab3A deletion and Rab3A earlybird mutant miceJ Neurosci 31:3580–3588
- 90.Activity-dependent presynaptic regulation of quantal size at the mammalian neuromuscular junction in vivoJ Neurosci 25:343–351
- 91.Activity coregulates quantal AMPA and NMDA currents at neocortical synapsesNeuron 26:659–670
- 92.Mediation of hippocampal mossy fiber long-term potentiation by cyclic AMPScience 265:1878–1882
- 93.Postsynaptic expression of homeostatic plasticity at neocortical synapsesJ Neurosci 25:2895–2905
- 94.Presynaptic regulation of quantal size by the vesicular glutamate transporter VGLUT1J Neurosci 25:6221–6234
- 95.Postsynaptic synaptotagmins mediate AMPA receptor exocytosis during LTPNature 544:316–321
- 96.Homeostatic mechanisms regulate distinct aspects of cortical circuit dynamicsProc Natl Acad Sci U S A 117:24514–24525
- 97.EEA1 restores homeostatic synaptic plasticity in hippocampal neurons from Rett syndrome miceJ Physiol 595:5699–5712
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
Copyright
© 2023, Koesters et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 523
- downloads
- 6
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.