Rapid active zone remodeling consolidates presynaptic potentiation

Neuronal communication across synapses relies on neurotransmitter release from presynaptic active zones (AZs) followed by postsynaptic transmitter detection. Synaptic plasticity homeostatically maintains functionality during perturbations and enables memory formation. Postsynaptic plasticity targets neurotransmitter receptors, but presynaptic mechanisms regulating the neurotransmitter release apparatus remain largely enigmatic. By studying Drosophila neuromuscular junctions (NMJs) we show that AZs consist of nano-modular release sites and identify a molecular sequence that adds modules within minutes of inducing homeostatic plasticity. This requires cognate transport machinery and specific AZ-scaffolding proteins. Structural remodeling is not required for immediate potentiation of neurotransmitter release, but necessary to sustain potentiation over longer timescales. Finally, mutations in Unc13 disrupting homeostatic plasticity at the NMJ also impair short-term memory when central neurons are targeted, suggesting that both plasticity mechanisms utilize Unc13. Together, while immediate synaptic potentiation capitalizes on available material, it triggers the coincident incorporation of modular release sites to consolidate synaptic potentiation.

N eurotransmitter-laden synaptic vesicles (SVs) release their content at presynaptic active zones (AZs) in response to Ca 2+ influx through voltage gated channels that respond to action-potential (AP) depolarization. Neurotransmitter binding to postsynaptic receptors subsequently leads to their activation for synaptic transmission. Modulation of transmission strength is called synaptic plasticity. Long-term forms of synaptic plasticity are major cellular substrates for learning, memory, and behavioral adaptation 1,2 . Mechanisms of long-term synaptic plasticity modify the structure and function of the presynaptic terminal and/or the postsynaptic apparatus. AZs are covered by complex scaffolds composed of a conserved set of extended structural proteins. ELKS/Bruchpilot (BRP), RIM, and RIM-binding protein (RBP) functionally organize the coupling between Ca 2+ -channels and release machinery by immobilizing the critical (M)Unc13 release factors in clusters close to presynaptic Ca 2+ -channels and thus generate SV release sites, at both mammalian and Drosophila synapses [3][4][5][6][7][8][9][10][11][12] . Whether and how discrete AZ release sites and the associated release machinery are reorganized during plastic changes remains unknown.
One crucial form of presynaptic plasticity is the homeostatic control of neurotransmitter release. This process, referred to as presynaptic homeostatic potentiation (PHP), is observed in organisms ranging from invertebrates to humans, but is perhaps best illustrated at the larval neuromuscular junction (NMJ) of Drosophila melanogaster 13,14 . Here, PHP requires the core AZscaffolding proteins RIM, RBP and Fife [15][16][17] and physiologically coincides with the upregulation of SV release sites 17,18 . Yet it is unknown how these AZ-scaffolds mediate release site addition, which downstream molecules are needed for PHP, whether AZscaffold independent reactions occur, and whether these mechanisms extend to other forms of plasticity, e.g., during learning in the central nervous system.
Here, we combine genetic and electrophysiological analysis to reveal a molecular sequence that triggers structural remodeling of AZ scaffolding proteins (BRP/RBP), and which ultimately leads to Unc13 addition within minutes. Using super-resolution light microscopy, we identify a modular AZ nano-architecture built by these proteins (which correspond to SV release sites) that rapidly extends by incorporating additional modules for plasticity. This rapid remodeling critically depends on the core AZ scaffolding proteins RBP/BRP, but neither on the early AZ assembly factors Liprin-α/Syd-1, nor on RIM or Fife. Additionally, AZ-remodeling was abolished in transport mutants previously shown to promote BRP/RBP transport. Strikingly, rapid addition of AZ nano-clusters was not required for the immediate expression of PHP on a minutes' timescale, but was essential to sustain potentiation thereafter. We identify Unc13A as a direct molecular target for PHP in experiments where Unc13A was delocalized from the AZ scaffolds. This mutant displayed sizable synaptic transmission but completely lacked PHP and AZremodeling. The same interference in mushroom body Kenyon cells of the Drosophila brain eliminated short-term memory, indicating that Unc13A is also a plasticity target in the central nervous system. In summary, we show that synapses capitalize on the available AZ material for immediate potentiation, but coincidently undergo release site addition via modular building blocks to consolidate stable synaptic potentiation. Thus, our work lays a foundation that will help to understand the mechanisms of a likely conserved presynaptic plasticity process that is important for dynamically adjusting and stabilizing neurotransmission across multiple timescales.

Results
Homeostatic plasticity regulates AZ protein levels. As a robust paradigm for assessing presynaptic plasticity over different time scales, we focused on PHP, which is well characterized at Drosophila NMJs 13 . To induce plasticity on a timescale of minutes, postsynaptic ionotropic glutamate receptors were partially blocked using the non-competitive open-channel blocker Philanthotoxin-433 (PhTx) 19 (Fig. 1a-d). This reduces postsynaptic sensitivity to neurotransmitter release from single SVs (reflected in a reduction of the amplitude of spontaneously occurring "minis", single SV fusion events, Fig. 1b). Initially, this also leads to a proportional decrease in AP-evoked transmission, but in under 10 min, PHP increases the number of SVs released per AP (quantal content). This compensates for the postsynaptic interference, resulting in AP-evoked transmission comparable to baseline levels (Fig. 1b) 19 . To identify molecular adaptations during plasticity, we investigated whether the levels of any of the evolutionarily conserved AZ proteins were altered. We accordingly immunostained against BRP, RBP, Unc13A (we focused on Unc13A, the Unc13 isoform dominating evoked SV release at Drosophila NMJ synapses 4 ; flybase: unc-13-RA), Syx-1A, Unc18, and Syd-1 (as motoneuronally expressed Syd-1-GFP) ( Fig. 1c; Supplementary Fig. 1a). In agreement with previous observations 18,20 , we found that 10 min of PhTx treatment increased AZ BRP-levels by about 50% (Fig. 1c, d). In addition, we found that RBP, Unc13A, and Syx-1A increased by about 30%, 60%, and 65%, respectively (Fig. 1c, d). The AZ levels of RBP/BRP, Unc13A/BRP, and Syx-1A/BRP scaled proportionally over all AZ sizes (Supplementary Fig. 1b; Ctrl, black lines). This proportionality was preserved upon PhTx treatment (Supplementary Fig. 1b; PhTx, blue lines). Notably, the AZ-levels and distribution of the essential Sec1/(M)Unc18 family protein Unc18-which was recently found to function in PHP 21 -were unaffected (Fig. 1c, d), demonstrating specific up-regulation of a subset of AZ proteins. Another AZ protein, the assembly factor Syd-1, even displayed a slight reduction upon PhTx ( Supplementary Fig. 1a), further underscoring a high degree of specificity.
To verify that these AZ-adaptations were specific to functional glutamate receptor interference, and to address their relevance over a longer time window, we investigated larvae bearing mutations in a glutamate receptor subunit (Fig. 1e-h). Deletion of the high-conductance receptor subunit IIA (GluRIIA) resulted in a similarly reduced postsynaptic sensitivity to single SV fusion events (Fig. 1f) as PhTx treatment (Fig. 1b). Because under these circumstances PHP also increases quantal content to achieve similar AP-evoked transmission (Fig. 1f), gluRIIA mutants have extensively been used to investigate long-term PHP (over the 3-4 days of larval development) 22,23 . Immunostainings against BRP, RBP, Unc13A, and Syx-1A confirmed their (in this case larger) elevation on this longer timescale (compare Fig. 1g, h with 1c, d) (100%, 70%, 400% and 200%, respectively, compared to 50%, 30%, 60%, and 65% upon PhTx treatment). Unlike the stoichiometric increase observed for BRP/Unc13A within minutes, this long-term PHP revealed enhanced Unc13A AZincorporation ( Supplementary Fig. 1c). Another distinction was a remarkable reorganization and eight-fold increase of Unc18 (Fig. 1g, h; Supplementary Fig. 1d, e). Our data hence imply that considerable AZ restructuring occurs within minutes of PHP induction, which is further enhanced across longer-lasting timescales.
Rapid addition of release site modules during PHP. We next sought to investigate how the altered levels of AZ proteins during PHP affected their nanoscopic topology by super-resolution STED microscopy (x-y resolution~30-40 nm). As noted before, planar AZs revealed clearly distinguishable individual Unc13A/ BRP/RBP spots arranged in a ring-like geometry [4][5][6][7] (Fig. 2a). It was recently shown by single molecule imaging that these individual clusters likely contain several (probably few tens of) molecules in the case of BRP 24 . We detected the number of clusters per AZ in single AZ-images, which was largely heterogeneous for all three proteins, with a simple peak detection algorithm ( Fig. 2c-e (black bars)). However, in all cases the cluster number per AZ increased upon PhTx-treatment ( Fig. 2c-e (blue bars)), and slightly increased further in gluRIIA Null mutants (Supplementary Fig. 2a; brown bars). With increasing AZ-cluster number the AZ diameter (measured from the AZ center to the center of the clusters) also increased in both conditions (Supplementary Fig. 2b, c), consistent with previous STED analysis performed on BRP-rings 18 . Notably, the observed remodeling only affected cluster numbers, and did not alter their intensities ( Supplementary Fig. 2b, c). Consequently, the first conclusion of our analysis is that PHP increases the number of Unc13A/BRP/ RBP nano-clusters per AZ within minutes.
We also wanted to investigate whether the overall AZ-topology changed upon cluster incorporation. Notably, averaging of STEDimages was recently used to generate a three-dimensional model of an average synapse, displaying the mean protein localization at high resolution 25 . Thus, to compare the overall single AZtopology, AZ images were centered, sorted by the number of clusters, aligned by rotation, and averaged. Two different alignment methods were used. In the first procedure, images were simply rotated such that the cluster with the highest intensity was positioned at the top (Supplementary Fig. 3a and Methods for details). Even though this procedure only targeted a single pixel per image (the position of the brightest cluster), the remaining (lower intensity) clusters were often found in similar relative positions, such that averaging revealed a simple polygonal geometrical series ( Supplementary Fig. 3a regularity. In a more refined analysis, we simultaneously considered the position of all clusters per AZ (Supplementary Fig. 3c and Methods for details) which also revealed a simple geometrical pattern for Unc13A, BRP, and RBP ( Fig. 2f-h and Supplementary  Fig. 3c). This stereotypical arrangement was best seen for AZs containing two to six clusters but less clear for AZs containing more than that (which could mean that these are less regular; Supplementary Fig. 2d   Thus, these findings imply that complexes of the core AZscaffold form discrete nano-modular structures, which correspond to SV release sites, and that rapid presynaptic plasticity triggers their fast AZ-incorporation which is even further enhanced over longer-timescales. Impairing BRP/RBP transport disrupts rapid AZ-remodeling. The remarkable remodeling of AZ material within the short timeframe of PhTx-treatment (minutes) raised the question of how this is mechanistically achieved. We first considered whether local presynaptic protein translation could be required 26 . However, treatment of larvae with 50 μg/ml of the translation blocker cycloheximide (prior and during PhTx-treatment) did not disrupt structural remodeling of BRP or vGlut (vesicular glutamate transporter) at AZs ( Supplementary Fig. 5a), consistent with remodeling being translation-independent. Moreover, the functional increase in quantal content remained expressed in the presence of the blocker ( Supplementary Fig. 5b), consistent with previous reports 19,20 .
Because active kinesin-dependent protein transport is required for long-term homeostatic plasticity in gluRIIA Null mutants 27 , we asked whether BRP/RBP transport mechanisms might be employed for AZ remodeling. For this, we investigated proteins involved in BRP/RBP transport by their mutation which causes abnormal BRP/RBP accumulation in the moto-neuronal axons, such as Atg1 (Unc-51) 28 , serine-arginine (SR) protein kinase at location 79D (Srpk79D 29,30 ), and App-like interacting protein (Aplip-1, Jip1 or JNK interacting protein in mammals), a selective RBP transport-adaptor 31 (Fig. 3a, b). While we observed clear PhTx-induced BRP-/Unc13A-upscaling in wild-type controls as well as in atg1 mutants (Fig. 3a, b), remodeling was fully absent upon null-mutation of srpk79D, aplip-1 32 or in animals bearing an Aplip-1 point mutation that selectively prevents kinesin light chain interaction (aplip-1 ek4 ) 33 (Fig. 3a, b). Additionally, STED microscopy revealed that aplip-1 ek4 and srpk79D ATC mutants appeared to contain fewer BRP/Unc13A clusters per AZ on average than wild-type (compare Fig. 2c, d with Fig. 3c) and fully lost the capacity to increase cluster numbers upon PhTxtreatment (even a decrease was observed, Fig. 3c). We also discovered that upon motoneuronal Aplip-1 or Srpk79D knockdown, Unc13A-GFP co-accumulated with aberrant axonal BRP aggregates ( Supplementary Fig. 5c, d). Interestingly, a partial coaccumulation of BRP/Unc13A-GFP (but to a lower extent in comparison to the knock-downs) was also present in the control situation, in line with an at least partial co-transport that we occasionally observed in live-imaging experiments (Supplementary Fig. 5c-g; Movie 1).
We further wanted to elaborate on the involvement of active protein transport during PhTx-induced AZ-remodeling by interfering with the cytoskeletal tracks used for short-range   NMJs of segment A2-A5 from third instar larvae and quantification of synaptic levels of wild-type-1, aplip-1 Null , aplip-1 ek4 , wild-type-2, atg1 and srpk79D ATC labelled with the indicated antibodies without (Ctrl; black) and with 10 min of PhTx (+PhTx; blue) treatment. Two independent experiments were performed, wild-type-1 was used as control for aplip-1 Null and aplip-1 ek4 while wild-type-2 was used for atg1 Null and srpk79D ATC . c Average BRP and Unc13A cluster number per AZ either without (Ctrl, −PhTx; black) or with PhTx While the AZ levels of BRP/Unc13A were already slightly enhanced by Latrunculin B treatment ( Supplementary Fig. 6a, b), PhTx-treatment failed to induce the typical increase of the AZ levels of these proteins (in fact a reduction was observed, Supplementary Fig. 6a, b). Thus, the actin cytoskeleton and active BRP/RBP transport are required for rapid AZ remodeling.
Rapid homeostatic AZ-remodeling depends on BRP and RBP. We next investigated which of the evolutionarily conserved core AZ scaffolding proteins are required for rapid AZ-remodeling. Loss of RBP fully blocked the rapid, PhTx-induced increase of BRP/Unc13A (Fig. 4a, b). BRP was also essential, because the typical increase in Syx-1A and Unc13A observed upon PhTx treatment was abolished (Compare Fig. 4a, b with Fig. 1c, d).
Notably, BRP-amounts appear to be rate-limiting because PhTxinduced AZ-remodeling was blocked in larvae heterozygous for a brp null allele (brp Null /+) (Fig. 4a, b). In contrast, null mutation of RIM, which abolished PHP 17 , did not interfere with AZremodeling ( Fig. 4a, b). Furthermore, the simultaneous deletion of RIM and Fife (a possible RIM homologue which is required for PHP 15,34 ) did not interfere with AZ-remodeling. Thus, RIM and Fife appear to act downstream of BRP/RBP and are non-essential for AZ-remodeling (Fig. 4a, b). AZ assembly is initiated by the conserved scaffolding proteins Liprin-α and Syd-1, which both regulate AZ size [35][36][37][38][39] . We reasoned that AZ growth -as observed here during plasticitymay capitalize on the same molecular machinery as de novo AZ formation, and therefore tested whether BRP/Unc13A-scaling depended on those proteins. However, liprin-α Null and syd-1 Null mutants revealed normal PhTx-induced BRP/Unc13A-scaling ( Fig. 4a, b), indicating that these factors are dispensable. Thus, systematic investigation of evolutionarily conserved AZ scaffolding proteins reveals a selective dependence on the core AZscaffolds BRP and RBP for structural remodeling during plasticity.
AZ-remodeling is required for chronic -but not rapid-PHP. Several studies have shown that presynaptic release positively correlates with AZ-size 7,40-43 . Therefore, we expected that the increase of AZ-BRP/Unc13A observed upon PhTx treatment would functionally increase presynaptic release ( Fig. 1a-d). However, it is not entirely obvious whether the AZ-remodeling (which continues beyond the minutes' timescale during long-term PHP ( Fig. 1e-h)) would be essential for rapid PHP. For instance, loss of RBP was shown to occlude both AZ remodeling (Fig. 4) and the functional increase in quantal content 16 Unc13A Fig. 4 Rapid AZ-remodeling during PHP requires BRP and RBP. a, b Confocal images and quantification of synaptic intensities of muscle 4 NMJs of abdominal segment 2-4 from third instar larvae at wild-type, rbp Null , brp Null , brp Null /+, rim Null , rim Null /fife Null , liprin-α Null , and syd1 Null NMJs labelled with the indicated antibodies without (control; black) and with 10 min PhTx (+PhTx; blue) treatment. For wild-type, images and data were adapted and replotted from Fig. 1. Please note, although wild-type controls were performed in parallel to every mutant genotype to check for functional AZ-remodeling upon PhTx-treatment in each set of experiments, we do not show all WT-control here due to space limitations. Therefore, AZ-protein levels should not be compared between genotypes. Source data as exact normalized and raw values (also of additional Wild-type controls for genotypes where PhTx-treatment failed to induce AZ-remodeling), detailed statistics including sample sizes and P values are provided in the Source Data file. Scale bars: 5 µm. Statistics: Mann-Whitney U test. *P ≤ 0.05; **P ≤ 0.01; ***P ≤ 0.001; n.s., not significant, P > 0.05. All panels show mean ± s.e.m.
not go hand-in-hand in the case of RIM (and Fife) mutants, whose AZs remodel, but which cannot express PHP (increased quantal content 15,17 ). These observations prompted us to systematically investigate the relevance of AZ-remodeling for the rapid induction and sustained expression of PHP.
We first investigated the dependence of rapid PHP on the rapid remodeling of Unc13A/Syx-1A in brp Null larvae. Strikingly, AZremodeling was blocked (Fig. 4a, b), but functional PHP expression (an increased quantal content) persisted at levels comparable to the wild-type/control situation (Fig. 5a, b). In addition, mutants of Srpk79D (whose AZs also do not undergo PhTx-induced AZ-remodeling; Fig. 3a, b and replotted in Fig. 5d) were likewise able to increase their quantal content (Fig. 5c), showing that AZ remodeling can be uncoupled from rapid PHP expression. Thus, even though AZ-remodeling occurs on a similar time-scale, it is not required to rapidly enhance the quantal content in these cases.
We next investigated whether elevation of the AZ protein levels was required to consolidate the increased quantal content over chronic time-scales in gluRIIA Null mutants ( Fig. 5e-h). Indeed, PHP was severely impaired in brp Null , gluRIIA Null double mutants (Fig. 5f). In an otherwise wild-type background, the increase in quantal content (upon gluRIIA-null mutation) was much larger than when BRP was additionally deleted (Fig. 5f). We could ensure that the impairment was not due to the overall reduced release in brp Null mutants, as a loss to compensate for the gluRIIA ablation was also seen in srpk79D ATC mutants (which had synaptic transmission comparable to wild-type cells (Fig. 5g), no AZremodeling upon PhTx-treatment (Fig. 3a, b), intact PHP upon PhTx-treatment (Fig. 5c) and severely impaired PHP expression in gluRIIA Null (Fig. 5g)). Importantly, AZ-remodeling was also fully blocked in gluRIIA Null ; srpk79D ATC double mutants (Fig. 5h). An intermediate behavior was seen in the case of the aplip-1 ek4 mutant, ( Supplementary Fig. 7), possibly because other transport adapters might compensate in this situation. Together, this suggests that PHP rapidly increases neurotransmitter release through modulation of the available AZ components, but in addition immediately induces AZ-remodeling to ensure its consolidation.
Presynaptic potentiation requires Unc13A. We next sought to identify the molecular substrate of PHP. Previous experiments established a requirement for the α1 subunit of the voltage gated Ca 2+ -channel Cacophony (Cac) 19,44 . In line with this, the levels of Cac as well as the Ca 2+ -influx increase upon PhTx treatment (Fig. 6a, b) 44,45 . We furthermore investigated Unc13A. A slight PhTx-induced BRP-/RBP-scaling persisted upon Unc13A loss, but was weaker than in the wild-type situation (Fig. 6a, b), possibly due to slightly elevated BRP/RBP-AZ-levels already in the non-PhTx-treated unc13A Null situation 4 . Notably, Cac-levels were still increased, even to a slightly larger extent than in the wild-type situation (Fig. 6a, b). However, functional PHP, the increase in quantal content, was completely lost (Fig. 6c, d). This indicates that Unc13A-like RIM and RBP 16,17 -plays an essential role in the plastic enhancement of NT release during PHP.
Rapid PHP and learning rely on the Unc13A N-terminus. The observation that Unc13A is essential for PHP is fully consistent with the previous findings that RIM and RBP are required 16,17 , because these proteins likely function in Unc13A AZ recruitment and activation (see discussion). As in other species (M. musculus/C. elegans), this interaction depends on the (M)Unc13 N-terminus 4,46-50 . To investigate the functional relevance of the Unc13A N-terminus for rapid PHP and AZ-remodeling, we used an Unc13A mutant lacking the N-terminal AZ-localization sequence (named C-term-GFP; Fig. 6e), which uncouples Unc13A from the central BRP/RBP scaffold 7 (and therefore supposedly also uncouples SV fusion from a possible regulatory function of RIM -see discussion). Importantly, the magnitude of AP-evoked synaptic transmission in these mutants was largely restored compared to the detrimental effect of unc13 Null mutation (compare Fig. 6h (C-term-GFP; grey traces) with Fig. 6h (full-length Unc13A-GFP; black traces) and ref. 7 ). However, in contrast to control larvae (unc13 Null with full-length Unc13A-GFP rescue; Fig. 6f-i), C-term-GFP mutants (unc13 Null with C-term-GFP rescue; Fig. 6f-i) completely lacked AZ-remodeling (Fig. 6f, g) (note that BRP levels were already enhanced in the non-PhTxtreated group, Supplementary Fig. 8a, b), Furthermore, unlike in control larvae, no rescue of evoked transmission and no increase in quantal content was seen upon PhTx treatment, indicating that C-term mutants were deficient of functional PHP (Fig. 6h-i). This demonstrates a dependence of PHP on the Unc13A N-terminus.
Although synapses vary tremendously in their excitability, input/output relationship, and transmitter type, the presynaptic release machinery is remarkably conserved in most systems and across species 9 . We wondered whether the principles of rapid presynaptic adaptation of the peripheral nervous system might be utilized in other synapse types and in other forms of plasticity, like the ones involved in learning and memory formation. Because short-term memory functions on timescales comparable to the PhTx-induced rapid homeostatic plasticity 51 , we investigated whether the Unc13A C-term mutant also exhibited learning deficits.
We specifically expressed C-term-GFP either alone, or while simultaneously knocking down endogenous full-length Unc13A (Unc13A-RNAi), in mushroom body Kenyon cells (KCs). KC output synapses undergo learning-induced plasticity and are required for the formation of short-term memories 52 . Knockdown of endogenous Unc13A and expression of C-term-GFP were confirmed using antibodies against GFP (labelling C-term-GFP but not endogenous protein) and the Unc13A N-terminus (labelling the endogenous protein but not C-term-GFP where the N-terminal epitope is deleted) (Fig. 7a). Additionally, we confirmed the strong efficacy of the Unc13A-RNAi via Western blot (knock-down in the entire brain using the pan-neuronal elav-Gal4 driver; Supplementary Fig. 9a). Because the C-term-GFP construct largely rescued the detrimental effect of unc13 null mutation at the NMJ 7 , we predicted that general transmission from KCs should be functional for both the C-term-GFP and the C-term-GFP/ Unc13A-RNAi conditions. To verify this, we performed in vivo two photon Ca 2+ -imaging (GCaMP6f) experiments in adult flies to assess odor-evoked responses at the M4/6 (MBON-β′2 mp/ MBON-γ5β′2a/MBON-β2β′2a) mushroom body output neurons, a postsynaptic circuit element directly downstream of KCs 53,54 . Indeed, robust Ca 2+ transients in response to odor stimulation were observed in all genetic constellations tested (Fig. 7b). We conclude as a result that KC output synapses expressing C-term-GFP or C-term-GFP/Unc13A-RNAi are functional under naive conditions. Together with the finding that naive odor avoidance was not statistically different between all these groups (Supplementary Fig. 9b), this allowed us to test whether either condition would interfere with learning and memory.
We assessed short-term memory in adult Drosophila using classical aversive olfactory conditioning. Flies were trained by pairing an odor with an electric shock and learning was scored by subsequently assaying avoidance of that odor 55 . All control groups showed similarly robust memory performance, while the relevant C-term-GFP/Unc13A-RNAi mutant (where endogenous Unc13A is knocked down and replaced by the C-term mutant incapable of PHP (Fig. 6e-i)) showed severe short-term memory impairments (Fig. 7c). Whether these impairments are indeed a consequence of the loss of a similar plasticity mechanism as the one observed at the NMJ or whether they are also related to differences in the synaptic transmission profile (e.g., short-term plasticity which is also affected in this mutant 7 ) remains to be established. Nevertheless, our data clearly indicate that Unc13A is a target for similar forms of plasticity (Unc13A-RNAi knockdown also impaired short-term memory, Supplementary Fig. 9b, c). Thus, our data imply that structural, functional, and behavioral adaptations are linked and that different forms of presynaptic plasticity may converge on Unc13A.

Discussion
Synapses are able to modify their transmission strength by undergoing plastic changes. This synaptic plasticity is crucial for neuronal circuit adaptation including learning and memory processes 1,56 . Molecular mechanisms for postsynaptic plasticity have been defined in considerable detail 2 . However, presynaptic mechanisms also modulate transmission strength in many synapse types and species 13,51,57 . Homeostatic plasticity is a wellstudied form of presynaptic plasticity at the Drosophila NMJ where an enhancement of AP-evoked neurotransmitter release counterbalances decreased postsynaptic receptor sensitivity. A number of relevant signaling molecules and pathways, including BMP signaling, CaMKII signaling, TOR signaling, proteasomal degradation, and trans-synaptic signaling are required for this 13,19,[58][59][60] . These factors appear to converge on two principal avenues to enhance presynaptic transmitter release, via increased Ca 2+ channel amounts and AP-induced Ca 2+ influx 44,45,61 and secondly via an increase in the number of releasable SVs and their associated release sites 17,18,61 . Nevertheless, some conditions were observed where Ca 2+ influx or Ca 2+ -channel levels were increased but the quantal content was not (Fig. 6a-d and ref. 17 ), suggesting that release site addition or activation is a required contributor. On the minutes' time-scale, structural AZremodeling was observed, yet whether and how this contributes to the enhancement of NT-release remained unclear 18 .
In the present study, we uncover a sequence of presynaptic molecular events that mediate AZ-remodeling (Fig. 8). We identify the presynaptic cytomatrix as a highly dynamic structure that can add discrete nano-modules of core proteins within minutes. In the initial phase of this structural remodeling, RBP Shaded areas represent the SEM. c Short-term memory scores after mushroom body-specific C-term-GFP rescue after Unc13A downregulation via locally driven RNAi expression (Ok107::UAS-C-term-GFP::Unc13A RNAi ) compared to controls expressing the driver, but not the RNAi (Ok107::+), the RNAi without driver (UAS-Unc13A RNAi ::+) or mushroom body-specific overexpression of the C-term-GFP construct (Ok107::UAS-C-term-GFP). See also Supplementary  Fig. 9. Source data as exact raw values, detailed statistics including sample sizes and P values are provided in the Source Data file. Scale bars: a 20 µm; b 10 µm. Statistics: nonparametric one-way analysis of variance (ANOVA) test, followed by a Tukey's multiple comparison test. ***P ≤ 0.001; n.s., not significant, P > 0.05. All panels show mean ± s.e.m. For representative images experiments were repeated twice with at least 6-7 brains per genotype and BRP are needed as their loss occludes an increase in Syx-1Aand Unc13A levels (Fig. 4). Somewhat unexpectedly, although AZ-remodeling occurs on the same minutes' time-scale, it is dispensable for the rapid potentiation of NT-release during PhTxinduced PHP (e.g., brp Null , aplip-1 ek4 or srpk79D ATC ) (Fig. 6). However, the remodeling is essential for the long-term consolidation of the potentiation (the capacity to restore the APevoked response in gluRIIA Null mutants was severely impaired when combined with brp Null or srpk79D ATC ) (Fig. 6). Thus, the rapid release enhancement appears to capitalize on AZ-material that is already present, for example by increasing Ca 2+ influx 44 and by the activation of already present but dormant release sites consistent with PHP depending on RIM, RBP, and Unc13A (refs. 16,17 and this study). Interestingly, the rapid potentiation coincides with the accumulation of BRP, Cac, RBP, Unc13A, Syx-1A, and (later) Unc18 in the AZ, which is required to consolidate and possibly extend the release enhancement. Therefore, the synapse utilizes two coincident programs which together ensure immediate rescue and also supply the synapse with backup-material in the form of BRP/RBP/Unc13A nanomodules in case the disturbance should persist.
Notably, recent work using STED microscopy to characterize hippocampal synapses also identified AZ nano-modules by clusters of Bassoon, vGlut, and Synaptophysin, and observed a scaling of vGlut and Synaptophysin upon chemically induced LTP 62 . This nicely aligns with the structural AZ-remodeling described here, suggesting an evolutionarily conserved process that tunes synaptic transmission by adding nano-modular structures to both sides of the synapse. In addition, activity dependent alterations in Syx-1A nano-clusters were also recently described 63 , further pointing to the AZ being a highly dynamic structure which adapts to different environmental demands.
How does the AZ scaffold remodel within minutes? While local protein translation is required for some forms of plasticity 64 , acute translation block did not interfere with PhTx induced AZ-remodeling ( Supplementary Fig. 5a). However, we found evidence that effective AZ protein transport and a functional cytoskeleton is a precondition here: The BRP/RBP transport adaptor/regulator proteins Aplip-1 and Srpk79D were required for rapid enhancement of BRP/RBP AZ levels (and loss of Aplip-1-mediated BRP/ RBP transport impaired short-term memory ( Supplementary  Fig. 9d, e)). Moreover, acute actin-depolymerization prevented the PhTx-induced BRP/Unc13A addition into AZs, further supporting a crucial role for their transport and in line with a recent study where Drosophila Mical, a highly conserved, multi-domain cytoplasmic protein that mediates actin depolymerization, was shown to be necessary for PHP 59 .
Considering the short timeframe of this adaptation, long-range transport appears unlikely. Instead, we favor the idea that Aplip-1 and Srpk79D function to engage an AZ-proximal reserve pool of components for rapid integration. This pool could originate from a local reservoir in the distal axon or terminal, or even between AZs, from which plasticity may trigger integration into established AZs 18 . Transport processes may fill or empty the reservoir. Between AZs, the reservoir could be composed of diffusely distributed proteins falling below the detection limit 65 , or could reflect a local rearrangement of material (note the reduction of AZs containing few AZ-protein modules after PhTx treatment (Fig. 2c-e) or upon gluRIIA ablation (Supplementary Fig. 2a)). Reducing the amount of BRP (by removing one gene copy) blocked the rapid structural adaptation (Fig. 4a, b), possibly because all available material was required to build AZs of proper functionality leaving no material for the reservoir. Regardless of the specific molecular mechanism, guided active transport along the cellular cytoskeleton or rearrangements of the cytoskeleton itself appear to serve a general function in synaptic plasticity in multiple species 59,[66][67][68] .
RIM and RBP are established targets for multiple forms of presynaptic plasticity in several synapse types and species 10 (Fig. 6). In fact, with our data we can infer the interrelation between these factors in presynaptic plasticity. RIM proteins are known to activate (M)Unc13s in several species [46][47][48]69,70 . While C. elegans and mouse (M)Unc13 proteins interact via an N-terminal C2A domain with RIM, no such domain is known for Drosophila, but the principal functional interaction may well be conserved (via another region of the N-term). This was directly tested by an Unc13A mutant, whose N-term was deleted 7 . This mutant largely rescued the severe loss of synaptic transmission in the unc13 Null condition with APevoked transmission comparable to the full-length Unc13A rescue (compare black with grey traces in Fig. 6h). However, this mutant completely lacked the capacity to undergo PHP, consistent with a required RIM/RBP/Unc13A interplay for plasticity. Furthermore, expressing the same mutant in the Drosophila memory center, the mushroom body, severely impaired shortterm memory formation (Fig. 7), pointing to a relevance of the RIM/RBP/Unc13 plasticity module in the Drosophila central nervous system. Thus, the morphological and molecular similarities between long-term sensitization in the famous sea slug Aplysia, presynaptic LTP in the mammalian brain, and homeostatic plasticity or learning in Drosophila indicate that the sequence of molecular events we describe here might be highly conserved.
For translation block experiments semi-intact preparations were pre-incubated in HL3 containing 50 μg/ml Cycloheximide (Chx; Sigmal-Aldrich, MO, USA) dissolved in DMSO for 10 min, a concentration previously shown to block protein synthesis in Drosophila 80  For Western blots 81 , adult flies were dissected in cold Ringer's solution and homogenized in Lysis buffer (1× PBS, 0.5%Triton, 2%SDS, 1× Protease inhibitor, 1× Sample buffer) followed by full-speed centrifugation at 18 ℃. One brain's supernatant for each group was subjected to SDS-PAGE and immunoblotted according to standard procedures. The following antibodies were used: guinea-pig Unc13A (1:2000 4 ) and mouse Tubulin (1:100,000; Sigmal-Aldrich, MO, USA; Cat# T9026, AB Registry ID: AB_477593). Antibodies obtained from the Developmental Studies Hybridoma Bank were created by the NICHD of the NIH and maintained at The University of Iowa, Department of Biology, Iowa City, IA 52242.
Image acquisition, processing, and analysis. Confocal microscopy was performed with a Leica SP8 microscope (Leica Microsystems, Germany). Images of fixed and live samples were acquired at room temperature. Confocal imaging of NMJs was done using a z-step of 0.25 μm. The following objective was used: 63 × 1.4 NA oil immersion for NMJ confocal imaging. All confocal images were acquired using the LAS X software (Leica Microsystems, Germany). Images from fixed samples were taken from muscle 4 of third instar larval 1b NMJs (segments A2-A5) or nerve bundles (segments A1-A3). Images for figures were processed with ImageJ software to enhance brightness using the brightness/contrast function. If necessary, images were smoothened (0.5 pixel Sigma radius) using the Gaussian blur function. Confocal stacks were processed with ImageJ software (http://rsbweb. nih.gov/ij/). Quantifications of AZs (scored via BRP) were performed following an adjusted manual 82 , briefly as follows. The signal of a HRP-Cy5 antibody was used as template for a mask, restricting the quantified area to the shape of the NMJ. The original confocal stacks were converted to maximal projections, and after background subtraction, a mask of the synaptic area was created by applying a certain threshold to remove the irrelevant lower intensity pixels. The segmentation of single spots was done semi-automatically via the command "Find Maxima" embedded in the ImageJ software and by hand with the pencil tool and a line thickness of 1 pixel. To remove high frequency noise a Gaussian blur filter (0.5 pixel Sigma radius) was applied. The processed picture was then transformed into a binary mask using the same lower threshold value as in the first step. This binary mask was then projected onto the original unmodified image using the "min" operation from the ImageJ image calculator. For spots/μm² the number of spots was divided by the size of the mask of the synaptic area.
For colocalization analysis (Pearson correlation coefficient) the ImageJ plugin "JACOP" (http://rsb.info.nih.gov/ij/plugins/track/jacop2.html) was used. To determine the synaptic protein levels, a custom-written ImageJ script was used that detects the locations with highest local maxima in pixel values to generate regions of interests (ROIs) and sets a point selection at each. The intensities were then measured and all selections were deleted, leaving intensity values and (x, y) locations in the results list. The results list was then used to create a circle of size = 5 pixels (pixel size 100 nm) centered around each (x, y) location and the integrated density within these ROIs was measured and taken for further calculations. The same ROIs were then used in the channel containing signals of the co-stained protein. For scatter plots, co-stainings of BRP either with RBP, Unc13A or Syx-1A with or without PhTx-treatment or wild-type and gluRIIA Null were used. The AZ numbers were counted (number of BRP spots) and the local synaptic levels of both co-stained proteins were measured. AZs were then sorted into five bins (AZ number divided by five) depending on their synaptic BRP levels and then the respective second channel intensities distributed to the appropriate bin. Binned BRP levels were then plotted against binned levels of the second channel.
STED microscopy. Two-color STED images were recorded on custom-built STEDmicroscopes 83,84 , which either combine two pairs of excitation laser beams of 595 nm and 635 nm or 595 nm and 640 nm wavelength with one STED fiber laser beam at 775 nm. All STED images were acquired using Imspector Software (Max Planck Innovation GmbH, Germany). STED images were processed using a linear deconvolution function integrated into Imspector Software (Max Planck Innovation GmbH, Germany). Regularization parameters ranged from 1e −09 to 1e −10 . The point spread function (PSF) for deconvolution was generated by using a 2D Lorentz function with its half-width and half-length fitted to the half-width and half-length of each individual image. For Fig. 3c, dual-color STED imaging with time-gated detection was performed using a commercial Leica SP8 TCS STED microscope (Leica Microsystems, Germany) equipped with a 100× NA 1.4 objective (HC PL Apo CS2; Leica Microsystems, Germany). Briefly, the system includes an inverted DMi8 CS microscope equipped with a 100× pulsed white light laser (WLL; ∼80-ps pulse width, 80-MHz repetition rate; NKT Photonics, Denmark) with a STED lasers for depletion (pulsed) at 775 nm. Detection of Alexa594 after excitation at 594 nm and emission detection of 604-650 nm and Atto647 after excitation at 640 nm at emission of 656-751 nm was performed in frame sequential mode. Time-gated detection with Hybrid detectors was set from 0.3-6 ns for both dyes. Raw data were deconvolved with Huygens Professional software (Scientific Volume Imaging) using a theoretical PSF automatically computed based on pulsedwave STED optimized function and the specific microscope parameters. Default deconvolution settings were applied. Images for figures and for finding highintensity clusters (see below) were processed with ImageJ software to remove obvious background or neighboring AZs (if required), enhance brightness/contrast and smoothened (0.5 pixel Sigma radius) using the Gauss blur function.
Classification and alignment of single AZs. All analysis described below was done using MATLAB R2016b (Mathworks Inc., Natick MA, USA) with optional toolboxes, which are indicated in the respective sections. Classification of individual AZs was achieved by using a custom script to detect the position of cluster centers (local intensity maxima) in the images where speckles of neighboring AZs were removed (only for Unc13A, BRP, RBP), while the averaging procedure described hereafter was performed on the corresponding raw images. This code retains only pixels above a defined grey value threshold. For the analysis of Unc13A stainings shown in Fig. 2f, a threshold value of 25 was used (18 for BRP, RBP, Syx-1A, and Unc-18). All pixel values in the image below this threshold value were set to zero to remove background noise. We then identified the positions of highintensity pixel clusters as local intensity maxima in the images. This was achieved by finding local maxima in the vertical and horizontal pixel lines. First, the function searched the first derivative (using the function diff) of all pixel columns for zero values and the changes in the surrounding slopes were considered to identify local maxima. The same procedure was then applied in pixel rows, but only for those pixel column values associated with a local maximum in the previous step. All single pixels that were associated with a maximum in both a row and a column were detected using the function intersect. To prevent detection of the same cluster more than once, a defined minimum distance of clusters was used (50 nm for BRP, Unc13A, Syx-1A and 20 nm for RBP) and only the local maximum with the highest intensity value was considered. All subsequent translation and averaging procedures were performed on corresponding raw images of the same AZs using the determined classification and positions of clusters. The non-cleaned AZs were sorted by the number of protein clusters detected this way. To calculate the center of mass of all coordinates, the means of all x-and y-coordinates were taken according to equations (1) and (2).
Sy ¼ n À1 Á X n 1 y obs ðnÞ ð 2Þ where (Sx, Sy) is the (x, y)-coordinate center of mass in the initial image, n is the number of identified clusters, and x obs (n) and y obs (n) are the positions of the n-th cluster in the present image. To align the center of mass to the center of the image, the necessary shift (Δx and Δy) of the original coordinates was calculated according to Eqs. (3) and (4) and used subsequently in Eqs. (5) and (6).
x centered ðnÞ ¼ x obs ðnÞ þ Δx ð5Þ In Eqs. (3) and (4), "imsize" refers to the size of the image in x or y dimension. The resulting coordinates x centered and y centered represent cluster coordinates after shifting the original coordinates x obs and y obs . The same translation was applied to the corresponding AZ image (using the function imtranslate, part of the 'Image Processing' toolbox). Clusters were ranked in a counter-clockwise sequence in all images (see illustration in Fig. 2b; Supplementary Fig. 3c) by sorting them for increasing angle between the image center and the cluster location in relation to the vertical midline (x = 26).
To align protein cluster coordinates between all investigated AZs, central rotation (around the image center, x middle = y middle = 26), was done using the operation in Eq. (7).
x rotated ðnÞ To find the optimal angle to overlay all AZs, a cost reflecting the sum of distances between cluster positions of the same rank in all images was minimized. The cost function was defined as described in Eq. (8).
ððx rotated ðn; mÞ À x rotated ðn; lÞÞ 2 þðy rotated ðn; mÞ À y rotated ðn; lÞÞ 2 Þ In Eq. (8), n is the cluster number, totClusters is the total cluster number of the respective category of images, m and l are particular AZ images in the stack of images from one category, and totImgs is the total number of AZ images in that category. The squared Euclidean distances were calculated using the function pdist, which performs the operation shown in Eq. (8).
The optimal rotation angle was found using a genetic algorithm function (ga, part of the 'Global Optimization' toolbox). The rotation angles evaluated were constrained in a range from −180 to 180 degrees. For faster optimization, parallelization (setting the option 'UseParallel' to 'true'; this requires the 'Parallel Computing' toolbox) was employed to evaluate 500 individual cost functions per generation. Cluster coordinates from all images in one category were rotated simultaneously. A single individual in the genetic algorithm represented a set of rotation angles for each image. The convergence criterion (TolFun) was left at the default value (a relative cost value change of <10 6 over 50 generations). The output of this optimization was a vector containing all rotation angles that led to the best overlap of cluster coordinates. Finally, these rotations were then applied to the centered original AZ images. All images aligned this way were then combined in a stack and an average image was generated by calculating the mean intensity of all image pixels. For better illustration of the AZ structure, pixel intensities were linearly scaled such that the highest intensity pixel had a value of 255. The procedure was only performed if more than two images existed in the same cluster number class for at least 5 consecutive cluster number classes. The histograms shown in Fig. 2c-e and Supplementary Fig. 2a were generated by counting the number of images in each cluster number class, and dividing each value by the total amount of images detected in all classes. The mean cluster intensity shown in Supplementary Fig. 2b, c was calculated by taking the average image intensities per category and subtracting the mean intensity of each category from the next higher one. These differences in mean intensity between subsequent categories were then further averaged.
To investigate the AZ structure in an approach independent of the clusterdistance minimization procedure described above, we repeated the averaging in a different way as follows. We developed a MATLAB code for AZ centering and alignment by rotation of the highest intensity pixel to identical angles and therefore similar positions, which yielded qualitatively similar results ( Supplementary  Fig. 3a). Again, high-intensity clusters were detected in AZ images cleaned from clusters belonging to bordering AZs, and all translations and rotations were then performed on unretouched images. The center of mass of the found cluster coordinates was calculated and the image shifted so that the center of mass of the coordinates was in the center of the 51 by 51 pixel space (x = y = 26; see Eqs. (1) to (6)). Only the position of the brightest pixel was then considered for rotation. To determine the angle by which to rotate each image to place this highest intensity pixel to the same fixed position (on the vertical midline between the two top image quadrants-the "twelve o'clock" position), the x and y distances to the center of rotation (equal to the center pixel of the image x = y = 26) were calculated to find the length (l) of the hypotenuse and the opposite side of the right triangle using pdist. The angle α in degrees was then calculated by taking the inverse sine in degrees of this value (MATLAB function asind), as shown in Eq. (9).
In cases where the brightest peak was located above the horizontal midline, the adjacent side of the triangle was the vertical midline. In cases where the brightest peak was located below the horizontal midline, the angle was calculated with the horizontal midline being the adjacent side, and 90°were added to the final angle value. Additionally, in cases where the brightest peak was located to the right of the vertical midline, the angle was multiplied with −1. To generate the results shown in Supplementary Fig. 3b, which shows the averaging of AZ images within randomly assigned categories, we generated a number of random category values and then proceeded with the averaging procedure described above. For this, we reproduced the distribution of category values from the AZ dataset according to the histogram values of the category vector as follows. A histogram of the category vector was generated (MATLAB function histogram) with a bin width of 1, yielding the absolute amount of AZs per category. A cumulative sum vector was calculated from these histogram values (MATLAB function cumsum). For each position in the category vector, we then chose a random number between 0 and 1 and multiplied it by the number of images. We then found the first position in the cumulative sum vector that was larger than this random value. The position found was equal to the assigned category. This resulted in a randomly assigned category vector with a similar distribution of categories as the original vector.
In vivo live imaging and analysis. In vivo imaging of intact Drosophila larvae was performed as previously described 85 . Briefly, third instar larvae were put into a drop of Voltalef H10S oil (Arkema, Inc., France) within an airtight imaging chamber. Before imaging, the larvae were anaesthetized with 20 short pulses of a desflurane (Baxter, IL, USA) air mixture until the heartbeat completely stopped. For assessing axonal transport, axons immediately after exiting the ventral nerve cord were imaged for 5 min using timelapse confocal microscopy. Kymographs were plotted using a custom-written ImageJ script.
Induction of homeostatic plasticity and electrophysiology. Third-instar larvae were selected and placed individually on a Sylgard block. Using a very sharp pin, the tail of the larva was pinned between the posterior spiracles, in the absence of solution. The head was pinned, making sure not to stretch the larva, so that the animal was relatively loose between the two pins. A small horizontal incision was made in the dorsal cuticle at the tail with a sharp scissors. The larva was cut vertically from the tail incision in an anterior direction (towards the head), continuing beyond the head pin. Great care was taken not to stretch the cuticle or animal during this process. Forty microliters of a 20 μM PhTx in modified hemolymph-like solution (HL3 86 ; composition (in mM): NaCl 70, KCl 5, MgCl 2 10, NaHCO 3 10, trehalose 5, sucrose 115, HEPES 5, CaCl 2 0, pH adjusted to 7.2) was pipetted into the abdominal cavity with minimal force, making sure to fill the abdomen. After 10 min incubation, the preparation was completed without rinsing. The cuticle was gently pinned twice on each side (without stretching). The connection of the intestines and trachea to the body at the posterior were cut. Holding the now free ends of the intestines and trachea with a fine forceps, remaining connections were cut moving in an anterior direction (towards the head). The intestines and trachea could then be gently removed without stretching the larva. Finally, the brain was held firmly and slightly raised above the body so that the scissors could be placed underneath to cut the segmental nerves. Care was taken not to touch the underlying muscle and to avoid excessive pulling of the nerves before they were cut. The completed preparation was rinsed 3 times with PhTx-free HL3 (0 mM CaCl 2 , 10 mM MgCl 2 ). PhTx-free HL3 solution was used in control treatments. Sylgard blocks were kept separate for PhTx and control treatments and all implements were rinsed after each recording. The Sylgard block and completed larval preparation was placed in the recording chamber which was filled with 2 ml HL3 (0.4 mM CaCl 2 , 10 mM MgCl 2 ). Recordings were performed at room temperature (~22°C) in current clamp mode at muscle 6 in segments A2/A3 as previously described 86 using an Axon Digidata 1550 A digitizer, Axoclamp 900 A amplifier with HS-9A x0.1 headstage (Molecular Devices, CA, USA) and on a BX51WI Olympus microscope with a 40X LUMPlanFL/IR water immersion objective. Sharp intracellular recording electrodes were pulled using a Flaming Brown Model P-97 micropipette puller (Sutter Instrument, CA, USA) with a resistance of 20-35 MΩ, back-filled with 3 M KCl. Cells were only considered with a membrane potential less than -60 mV and membrane resistances greater than 4 MΩ. All recordings were acquired using Clampex software (v10.5) and sampled at 10-50 kHz, filtering with a 5 kHz lowpass filter. mEPSPs were recorded for 1 min. eEPSPs were recorded by stimulating the appropriate nerve at 0.1 Hz, five times (8 V, 300 μs pulse) using an ISO-STIM 01D stimulator (NPI Electronic, Germany). Stimulating suction electrodes were pulled on a DMZ-Universal Puller (Zeitz-Instruments GmbH, Germany) and fire polished using a CPM-2 microforge (ALA Scientific, NY, USA). A maximum of two cells were recorded per animal.
Analysis was performed with Clampfit 10.5 and Graphpad Prism 6 software. mEPSPs were further filtered with a 500 Hz Gaussian low-pass filter. Using a single template for all cells, mEPSPs were identified and analyzed, noting the mean mEPSP amplitude per cell. For Fig. 5 and Supplementary Fig. 7, templates were generated for each cell and the first 30 mEPSPs were identified and taken into account for further analysis. An average trace was generated from the 5 eEPSP traces per cell. The amplitude of the average eEPSP trace was divided by the mean mEPSP amplitude, for each respective cell, to determine the quantal content.
Dissection and current clamp recordings of w1118 vs gluRIIA Null were performed as above in male third-instar larvae. Cells with an initial membrane potential greater than −55 mV, resistances less than 5 MΩ or multiple responses to a single stimulus were rejected. eEPSPs were recorded by stimulating the appropriate nerve at 0.2 Hz, 10 times (6 V, 300 μs pulse). An average eEPSP amplitude was calculated from the 10 traces. mEPSPs were analyzed with a genotype specific template. Quantal contents were calculated by dividing the mean eEPSP by mean mEPSP for each cell.
In vivo two-photon live calcium imaging and analysis. Two-photon imaging of odor-evoked calcium responses was conducted in 3-5-day-old mixed-sex flies expressing LexAop-GCaMP6f in VT1211-LexA. For imaging, flies were briefly anesthetized on ice and mounted in a custom made chamber by immobilizing wings, head and proboscis with wax. The head capsule was opened in sugar-free HL3-like extracellular saline 87 . Odor stimulation consisted of a 1.5 s OCT pulse followed by a 30 s break and then a 1.5 s MCH pulse followed again by another 30 s break. This alternating odor pulse protocol was consecutively repeated five times (odor dilution in mineral oil 1/1000). Odors were delivered on a clean air carrier stream and image acquisition and odor stimulation was synchronized temporally using a custom-designed system. Fluorescence was centered on 910 nm generated by a Ti-Sapphire laser (Chameleon Ultra II, Coherent, CA, USA). Images with a pixel size of 0.3 × 0.3 μm were acquired at 70 Hz using two-photon microscopy (Femto2D-Resonant by Femtonics Ltd., Hungary) with a 20×, 1.0 NA waterimmersion objective, controlled by MESc v3.5 software (Femtonics Ltd., Hungary). For each animal, a single hemisphere was analyzed. All OCT and MCH responses of a fly were averaged respectively and resulting traces were averaged between flies.
Mean intensity values of two-photon fluorescence were calculated, while F0 was defined as the mean F from 0 to 1.5 s at the beginning of a recording (MESc v 3.5 software). Image processing for single frames were manually performed using ImageJ (images did not require registration).
Odor avoidance conditioning. All flies were 3 to 5 days old, raised in 12 h:12 h, light:dark cycle and at 65% relative humidity. One day before the experiment, the flies were transferred to fresh food vials. One hour prior to the experiment, flies were pre-conditioned to experimental conditions (dim red light, 25°C, humidity of 80%). The aversive odors 3-Octanol (OCT) and Methylcyclohexanol (MCH) were diluted 1:100 in paraffin oil and presented in 14 mm cups. A current of 120 V AC was used as a behavioral reinforcer. The associative training was performed as previously described 88 . In a single-cycle training, nearly 100 flies were presented with one odor (CS + ) paired with electrical shock (US; 12 times for 1 min). After one minute of pure air-flow, the second odor was presented without the shock (CS − ) for another minute. The flies were then immediately tested for short-term memory performance by presenting them the two odors together. A performance index (PI) was calculated as the number of flies choosing the odor without shock (CS − ), minus the number of flies choosing the odor paired with shock (CS + ), divided by the total number of flies, multiplied by 100. The values of PI ranges from 0 to 100 where 0 means no learning (50:50 distribution of flies) and a value of 100 means complete learning (all flies avoided the conditioned odor). The final learning index was calculated as the average of both reciprocal indices for the two odors. Odor Avoidance experiments were used to test innate behavior where each odor was presented to the flies without conditioning. The PIs were calculated as stated above.
Quantification and statistical analysis. Data were analyzed using Prism (Graph-Pad Software, CA, USA). Per default Student's T test was performed to compare the means of two groups unless the data were either non-normally distributed (as assessed by D'Agostino-Pearson omnibus normality test) or if variances were unequal (assessed by F test) in which case they were compared by a Mann-Whitney U Test. However, in the Source Data file both tests are provided for all relevant cases. For comparison of more than two groups, one-way analysis of variance (ANOVA) tests were used, followed by a Tukey's multiple comparison test. P values and N values are given in the Source Data file. Means are annotated ± s.e.m. Asterisks are used to denote significance: *P < 0.05; **P < 0.01; ***P < 0.001; n.s. (not significant), P > 0.05.
Reporting summary. Further information on experimental design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The data that support the findings of this study, additional information and requests for resources and reagents as well as MATLAB and ImageJ codes used in this study are available from Alexander M. Walter (awalter@fmp-berlin.de) upon request.