Nitric oxide acts as a cotransmitter in a subset of dopaminergic neurons to diversify memory dynamics
Abstract
Animals employ diverse learning rules and synaptic plasticity dynamics to record temporal and statistical information about the world. However, the molecular mechanisms underlying this diversity are poorly understood. The anatomically defined compartments of the insect mushroom body function as parallel units of associative learning, with different learning rates, memory decay dynamics and flexibility (Aso and Rubin, 2016). Here, we show that nitric oxide (NO) acts as a neurotransmitter in a subset of dopaminergic neurons in Drosophila. NO’s effects develop more slowly than those of dopamine and depend on soluble guanylate cyclase in postsynaptic Kenyon cells. NO acts antagonistically to dopamine; it shortens memory retention and facilitates the rapid updating of memories. The interplay of NO and dopamine enables memories stored in local domains along Kenyon cell axons to be specialized for predicting the value of odors based only on recent events. Our results provide key mechanistic insights into how diverse memory dynamics are established in parallel memory systems.
Introduction
An animal’s survival in a dynamically changing world depends on storing distinct sensory information about their environment as well as the temporal and probabilistic relationship between those cues and punishment or reward. Thus, it is not surprising that multiple distributed neuronal circuits in the mammalian brain have been shown to process and store distinct facets of information acquired during learning (White and McDonald, 2002). Even a simple form of associative learning such as fear conditioning induces enduring changes, referred to as memory engrams, in circuits distributed across different brain areas (Herry and Johansen, 2014). Do these multiple engrams serve different mnemonic functions, what molecular and circuit mechanisms underlie these differences, and how are they integrated to control behavior? Localizing these distributed engrams, understanding what information is stored in each individual memory unit and how units interact to function as one network are important but highly challenging problems.
The Drosophila mushroom body (MB) provides a well-characterized and experimentally tractable system to study parallel memory circuits. Olfactory memory formation and retrieval in insects requires the MB (de Belle and Heisenberg, 1994; Dubnau et al., 2001; Erber et al., 1980; Heisenberg, 2003; McGuire et al., 2001). In associative olfactory learning, exposure to an odor paired with a reward or punishment results in formation of a positive- or negative-valence memory, respectively (Quinn et al., 1974; Tempel et al., 1983; Tully and Quinn, 1985). In the MB, sensory stimuli are represented by the sparse activity of ~2000 Kenyon cells (KCs). Each of 20 types of dopaminergic neurons (DANs) innervates compartmental regions along the parallel axonal fibers of the KCs. Similarly, 22 types of mushroom body output neurons (MBONs) arborize their dendrites in specific axonal segments of the KCs; together, the arbors of the DANs and MBONs define the compartmental units of the MB (Aso et al., 2014a; Mao and Davis, 2009; Tanaka et al., 2008). Activation of individual MBONs can cause behavioral attraction or repulsion, depending on the compartment in which their dendrites arborize, and MBONs appear to use a population code to govern behavior (Aso et al., 2014b; Owald et al., 2015).
A large body of evidence indicates that these anatomically defined compartments of the MB are also the units of associative learning (Aso et al., 2012; Aso et al., 2014b; Aso et al., 2010; Berry et al., 2018; Blum et al., 2009; Bouzaiane et al., 2015; Burke et al., 2012; Claridge-Chang et al., 2009; Huetteroth et al., 2015; Ichinose et al., 2015; Isabel et al., 2004; Krashes et al., 2009; Lin et al., 2014; Liu et al., 2012; Owald et al., 2015; Pai et al., 2013; Plaçais et al., 2013; Schwaerzel et al., 2003; Séjourné et al., 2011; Trannoy et al., 2011; Yamagata et al., 2015; Zars et al., 2000). Despite the long history of behavioral genetics in fly learning and memory, many aspects of the signaling pathways governing plasticity—especially whether they differ between compartments—remain poorly understood. Nevertheless, dopaminergic neurons and signaling play a key role in all MB compartments, and flies can be trained to form associative memories by pairing the presentation of an odor with stimulation of a single dopaminergic neuron (Aso et al., 2010). Punishment or reward activates distinct sets of DANs that innervate specific compartments of the MB (Das et al., 2014; Galili et al., 2014; Kirkhart and Scott, 2015; Liu et al., 2012; Mao and Davis, 2009; Riemensperger et al., 2005; Tomchik, 2013). Activation of the DAN innervating an MB compartment induces enduring depression of KC-MBONs synapses in those specific KCs that were active in that compartment at the time of dopamine release (Berry et al., 2018; Bouzaiane et al., 2015; Cohn et al., 2015; Hige et al., 2015; Owald et al., 2015; Séjourné et al., 2011). Thus, which compartment receives dopamine during training appears to determine the valence of the memory, while which KCs were active during training determines the sensory specificity of the memory (Figure 1A).
Compartments operate with distinct learning rules. Selective activation of DANs innervating specific compartments has revealed that they can differ extensively in their rates of memory formation, decay dynamics, storage capacity, and flexibility to learn new associations (Aso et al., 2012; Aso and Rubin, 2016; Huetteroth et al., 2015; Yamagata et al., 2015). For instance, the dopaminergic neuron PAM-α1 can induce a 24 hr memory with a single 1 min training session, whereas PPL1-α3 requires 10 repetitions of the same training to induce a 24 hr memory. PPL1-γ1pedc (aka MB-MP1) can induce a robust short-lasting memory with a single 10 s training, but cannot induce long-term memories even after 10 repetitions of a 1 min training. PAM-α1 can write a new memory without compromising an existing memory, whereas PPL1-γ1pedc extinguishes the existing memory when writing a new memory (Aso and Rubin, 2016). What molecular and cellular differences are responsible for the functional diversity of these compartments? Some differences might arise from differences among KC cell types (reviewed in Keene and Waddell, 2007; McGuire et al., 2005), but memory dynamics are different even between compartments that lie along the axon bundles of the same Kenyon cells (for example, α1 and α3). In this paper, we show that differences in memory dynamics between MB compartments can arise from the deployment of distinct cotransmitters by the DAN cell types that innervate them.
Results
Dopaminergic neurons can induce memories without dopamine, but with inverted valence
DANs release diverse cotransmitters in the mammalian brain (Maher and Westbrook, 2008; Stuber et al., 2010; Sulzer et al., 1998; Tecuapetla et al., 2010; Tritsch et al., 2012). In Drosophila, the terminals of the MB DANs contain both clear and dense-core vesicles (Takemura et al., 2017), prompting us to ask if the DAN cell types innervating different MB compartments might use distinct cotransmitters that could play a role in generating compartment-specific learning rules. We individually activated several DAN cell types in a tyrosine hydroxylase (TH) mutant background that eliminates dopamine synthesis in the nervous system (Cichewicz et al., 2017; Riemensperger et al., 2011) and assayed their ability to induce associative learning when paired with an odor stimulus.
We first examined the PPL1 cluster of DANs, which innervate several MB compartments involved in aversive learning, driven by stimuli such as electric shock, noxious temperature, and bitter taste (Galili et al., 2014; Kirkhart and Scott, 2015; Mao and Davis, 2009; Riemensperger et al., 2005; Tomchik, 2013). Using an optogenetic olfactory arena (Aso and Rubin, 2016), we trained flies by pairing odor exposure with optogenetic activation of these DANs using CsChrimson-mVenus and then immediately tested memory (Figure 1B,C). In flies with a wild-type TH allele, this training protocol induced robust negative-valence memory of the paired odor (Figure 1D), as observed previously (Claridge-Chang et al., 2009; Schroll et al., 2006). In the dopamine-deficient background, activation of the same DANs still induced a robust odor memory, but its valence was now positive (Figure 1E). This result is consistent with previous findings that TH mutant flies show a weak positive-valence memory after odor-shock conditioning (Riemensperger et al., 2011), although the positive-valence memory we observed is much stronger.
Arguing against the possibility that this valence-inversion phenotype resulted from a developmental defect caused by the constitutive absence of dopamine (Niens et al., 2017), feeding L-DOPA and carbidopa to adult-stage flies fully restored normal valence memory (Figure 1F and Figure 1—figure supplement 1). Nor did the valence-inversion phenotype result from lack of dopamine signaling outside the MB, as restoring TH expression specifically in the PPL1 DANs was sufficient to restore formation of negative-valence memory (Figure 1G; Figure 1—figure supplement 2). Moreover, valence-inversion in the absence of dopamine was not limited to punishment-representing DANs. Activation of reward-representing PAM cluster DANs (Burke et al., 2012; Liu et al., 2012) in the TH mutant background also induced odor memory of opposite valence, in this case negative rather than positive (Figure 1D,E); as we found for PPL1-induced memories, either L-DOPA plus carbidopa feeding or TH expression in reward-representing DANs restored the ability to form a memory of the valence that is observed in wild-type flies (Figure 1F,G). These observations suggested the possible presence of a cotransmitter in these DANs that exerts an opposite effect from dopamine on synaptic plasticity and memory.
Putative cotransmitter effects differ among DAN cell types
If DAN cell types use different cotransmitters, we might expect the effects of activating DANs in the TH mutant background to vary with cell type. We tested this idea by comparing the associative memories formed in wild-type and TH mutant backgrounds when an odor was paired with optogenetic activation of different subsets of DAN cell types using seven driver lines (Figure 2 and Figure 2—figure supplement 1). We also demonstrated that valence inversion is not limited to training using direct DAN stimulation with CsChrimson; activating bitter taste sensory neurons using Gr66a-GAL4, which activate PPL1-DANs (Das et al., 2014; Kirkhart and Scott, 2015), likewise induced memories of inverted valence (Figure 2). We identified two DAN cell types that exhibited the valence-inversion phenotype: PPL1-γ1pedc and PAM-γ5/PAM-β'2a. With PPL1-γ1pedc stimulation, memory valence switched from negative to positive in TH mutant animals. Conversely, with PAM-γ5 and PAM-β'2a stimulation, memory valence flipped in the opposite direction, from positive to negative in TH mutants.
The valence-inversion phenotype was not, however, observed in all compartments. Activation of PPL1-γ2α'1 resulted in negative-valence memory in both TH mutant and wild type, suggesting a cotransmitter with the same sign of action as dopamine. Activation of PPL1-α3 or PAM-β'1 in the TH mutant background did not induce significant memory, indicating that these cells do not express a cotransmitter capable of inducing memory without dopamine.
Identification of nitric oxide synthase in a subset of DAN cell types by transcript profiling
To identify potential cotransmitters, we profiled (using RNA-Seq methods) the transcriptomes of the DAN cell types in these seven split-GAL4 lines, and looked for candidates whose expression correlated with the valence-inversion effect. Isolation of pure populations of specific DAN cell types is challenging because of their low abundance. For example, PPL1-γ1pedc is a single cell in each brain hemisphere and thus requires ~50,000 fold enrichment. We used a collection of split-GAL4 driver lines (Aso et al., 2014a) to fluorescently mark the soma of specific DAN cell types and confirmed that the nuclear-targeted reporters we used for sorting visualized the same restricted set of cells as the membrane-targeted reporters in the original study. In this way, we selected a combination of reporter and driver lines that provided the most specific labeling of the targeted cell type. The number of detected genes and the correlation across biological replicates of RNA profiling experiments has been observed to be highly dependent on cDNA yield during library construction (Davis et al., 2018). Due to different soma sizes, the amount of mRNA per cell is expected to differ across cell types. To estimate the number of cells necessary for our experimental condition, we started with MBON-γ1pedc>α/β cells, a cell type that occurs once per hemisphere. Three replicates with ~350 cells yielded on average 5.3 µg of cDNA, and we observed a high correlation across biological replicates on these dataset (Pearson R = 0.90). Thus, we aimed for similar cDNA yields by sorting more cells for cell types with a smaller soma (i.e. KCs and PAM cluster DANs). We collected data from 10 driver lines for DAN cell types, with two to four biological replicates per line. We also examined three classes of KCs, DPM, APL and five additional MBON cell types. On average, we collected 1,000–2,500 cells for KCs, 56+ /- 18 for DPM and APL, 546 + /- 60 cells for PAM DANs, and 296 + /- 15 cells for PPL1 DAN and MBON cell types per replicate, by hand or fluorescence activated cell sorting (FACS), yielding 3.39 + /- 0.21 µg of cDNA, 17.7 + /- 1.2 million mapped reads per replicate, and 0.80 + /- 0.02 Pearson R across biological replicates (Figure 3—figure supplement 1; Supplementary file 1).
We analyzed these data for different splicing isoforms (see Materials and methods). Using DE-seq2 (Love et al., 2014), we searched for transcripts that were differentially expressed among DAN cell types and, in particular, for those commonly expressed in DAN cell types that showed the valence-inversion phenotype (PPL1-γ1pedc and PAM-γ5), but not in other cell types (Figure 3A,B). Through this analysis, we identified nitric oxide synthase (NOS) as a strong candidate for an enzyme synthesizing a cotransmitter. Similar enrichment in PPL1-γ1pedc and PAM-γ5 was found in only five other transcripts, none of which are likely to encode a neurotransmitter: (i) epac-RG, a cAMP-activated guanine nucleotide exchange factor. (ii) br-RO, br-RI, both transcripts of the zinc finger transcription factor broad, (iii) CG32547-RD, a G-protein-coupled receptor, and (iv) CG12717-RA, a SUMO-specific isopeptidase (data reviewed in FlyBase) (Thurmond et al., 2019). In addition to these transcripts that matched our criteria, we found other potential candidates whose expression was not a precise match. We detected a high level of the DH44 neuropeptide in PPL1-γ1pedc, but not in PAM-γ5. A receptor for DH44, DH44-R1, was expressed in PAM-γ4 and/or γ4<γ1γ2 and to lower extent in α/β Kenyon cells, but neither of the two known receptors for DH44 was detected in γ Kenyon cells (Figure 3—figure supplement 9). Transcripts of the neuropeptide gene Nplp1 were detected in PPL1-γ1pedc and PAM-γ5, but other DANs and MBONs also expressed this gene (Figure 3—figure supplement 8) (Croset et al., 2018). Expression of Gyc76C, the receptor for Nplp1, was barely detectable in KCs, DANs, and MBONs. Complete transcript data were deposited to NCBI Gene Expression Omnibus (accession number GSE139889) and presented in Supplementary file 1, and the expression of genes encoding neurotransmitters, neuropeptides and their receptors, as well as components of gap junctions, is summarized in Figure 3—figure supplement 2–10. Although we cannot formally rule out a contribution of other genes and pathways, we chose to pursue NOS, as it was the most promising candidate gene for producing a cotransmitter that might be responsible for the valence-inversion effect.
Drosophila has only one gene encoding nitric oxide synthase (Nos), but this gene has multiple splicing isoforms (Figure 3C) (Rabinovich et al., 2016; Regulski and Tully, 1995; Stasiv et al., 2001). Only NOS1, the full-length isoform, is functional, while the truncated isoforms can function as a dominant-negative. Thus, identifying the expressed splicing isoform of Nos was crucial for understanding NOS functions in DANs. NOS1 was the most abundantly expressed isoform as judged by RNA profiling. We confirmed NOS1 expression by combining fluorescent in situ hybridization (FISH) and antibody staining. For whole-brain FISH (Long et al., 2017), we used 40 probes against c-terminus exons that are present in NOS1 and NOS4, but not NOS-RK, transcripts (Figure 3C and Figure 3—figure supplement 1D). PPL1-γ1pedc was labeled with these FISH probes (Figure 3—figure supplement 1D), confirming expression of NOS1 or NOS4 in these cells. For immunohistochemistry, we used an antibody raised against exon 16 of NOS (Kuntz et al., 2017; Yakubovich et al., 2010) that is present in NOS1 and NOS-RK, but not in NOS4 (Figure 3C). We validated its specificity by demonstrating a loss of the staining that accompanied RNAi-mediated knockdown of NOS (Ni et al., 2011) in PPL1-γ1pedc (Figure 3D and Figure 3—figure supplement 1B). In the MB lobes, the γ1pedc and γ5 compartments showed enriched anti-NOS immunoreactivity, as expected from the RNA-Seq data (Figure 3E). Furthermore, we confirmed colocalization of NOS with γ1pedc’s terminals by expansion microscopy and lattice light sheet microscopy (ExM LLSM; Figure 3F and G, Video 1) (Chen et al., 2015; Gao et al., 2019; Tillberg et al., 2016). Although the majority of anti-NOS immunoreactivity in γ1 compartment was diminished by expressing NOS-shRNA in PPL1-γ1pedc, we noticed sparse but large NOS-positive puncta remained (Figure 3D). Similar large structures were observed outside the DAN’s terminals with ExM LSSM and were found to colocalize with terminals of octopaminergic neurons (Figure 3F right). In addition, γ3 and γ4 also showed significant anti-NOS immunoreactivity (Figure 3E). However, the low NOS transcript levels observed via RNA-Seq in PAM-γ3 and PAM-γ4 are most consistent with anti-NOS immunoreactivity in γ3 and γ4 arising from non-DAN cell types or developmental expression. The cell-type-specific expression and localization of NOS1 in a subset of DANs associated with compartments that display the valence-inversion phenotype prompted us to test the hypothesis that NO plays a direct role in generating the diversity of memory dynamics observed in different compartments.
NOS in dopaminergic neurons contributes to memory formation
We next evaluated the role of NOS in memory formation in the absence of dopamine biosynthesis. If NO is indeed the cotransmitter that supports the valence-inverted memory in TH mutant flies, we would predict that inhibiting NOS should block this effect, and that flies would show no memory. To assess the requirement for NO synthesis, we fed flies the competitive NOS inhibitor L-NNA for one day before training and then measured the memory induced by optogenetic training using PPL1-γ1pedc in a TH mutant background. We found that this treatment reduced valence-inverted memory formation in a dose-dependent manner (Figure 4A). Whereas NOS is broadly expressed in the brain (Kuntz et al., 2017), two lines of evidence suggest that L-NNA fed flies are capable of olfactory learning. First, when we bypassed the TH mutant by feeding L-DOPA and carbidopa to restore dopamine levels, the L-NNA fed flies showed a normal level of negative-valence odor memory formation in response to pairing an odor with PPL1- γ1pedc activation (Figure 4A). Second, the effect of L-NNA feeding was cell-type specific; memory formed by activation of either PPL1-γ1pedc or PAM-γ5/β'2a was affected but that formed by activation of PPL1-γ2α'1, was not (Figure 4B). We obtained consistent results in knockdown experiments where we expressed NOS-RNAi in all PPL1 DANs (Figure 4C). In this case, the negative-valence memory observed with NOS-RNAi in the TH mutant background is consistent with the presence of a different cotransmitter in PPL1-γ2α'1 (Figure 2).
We also examined whether we could transfer this valence-inversion property to another compartment by ectopically expressing NOS. We expressed NOS in PPL1-α3, a compartment where we observed no intrinsic NOS expression, and where extended optogenetic training induces a negative-valence memory (Aso and Rubin, 2016). TH mutant flies formed no odor association with PPL1-α3 activation, but when NOS was ectopically expressed in the α3 compartment, the same training protocol induced a positive-valence memory (Figure 4D and E). Similarly, ectopic expression of NOS in PPL1-γ2α'1 also resulted in valence inversion phenotype (Figure 4—figure supplement 1). In other words, NO was able to form an association of opposite valence to that formed by dopamine. These results demonstrate the functional significance of NOS in DANs, but leave open its mechanism of action.
Soluble guanylate cyclase in Kenyon cells is required to form NO-dependent memories
In the MB, dopamine induces synaptic plasticity by binding to dopamine receptors on the axons of KCs and activating the rutabaga-encoded adenylyl cyclase in these cells (Gervasi et al., 2010; Tomchik and Davis, 2009). Does NO released from DANs also act on receptors in the KCs? RNA-seq data revealed expression of the subunits of soluble guanylate cyclase (sGC), Gycα99B and Gycβ100B, in KCs, DANs, and MBONs (Figure 5A–B). Similar to its mammalian homologs, Drosophila sGC is activated upon binding NO (Morton et al., 2005). Transcripts for other guanylyl cyclases were found at lower levels, if at all (Figure 5A–B; Supplementary file 1), suggesting that sGC formed by Gycα99B and Gycβ100B is the primary source of cGMP in KCs. Consistent with these transcriptomic data, we observed that a protein trap Gycβ100B-EGFP fusion protein (MI08892; Venken et al., 2011) was broadly expressed in the MB lobes (Figure 5—figure supplement 1).
We tested the role of Gycβ100B with acute RNAi knockdown in KCs using the MB-switch system (McGuire et al., 2003). This RNAi knockdown effectively reduced expression of endogenous Gycβ100B tagged with GFP in the MB (Figure 5—figure supplement 1C). We found that a reduction in Gycβ100B produced in this way abolished the valence-inverted memory induced by activation of PPL1-γ1pedc in a TH mutant background (Figure 5C and Figure 5—figure supplement 1D). Although MB-switch showed significant leaky expression without RU-486 feeding under our experimental conditions (Figure 5—figure supplement 1E and F), we could restore normal aversive memory by feeding TH mutant flies L-DOPA (Figure 5C), indicating that flies carrying MB-switch driven Gycβ100B-RNAi do not have a general defect in learning. Taken together, our data argue strongly that NO functions as cotransmitter that is released by DANs and acts on sGC in postsynaptic KCs to regulate cGMP levels, although our results do not exclude the possibility that NO has other targets. Indeed, both MBONs and DANs express sGC (Figure 5B) and we found that optogenetic stimulation of PPL1-γ1pedc in a TH-mutant background can induce a prolonged increase in calcium levels in MBON-γ1pedc>α/β, which was suppressed by L-NNA (Figure 5—figure supplement 2A). However, DANs expressing sGC in other compartments did not respond to stimulation of PPL1-γ1pedc in a TH-mutant background (Figure 5—figure supplement 2B), raising the possibility that NO may not be able to diffuse over the 30 µm distances required to reach those compartments.
NO-dependent and dopamine-dependent memories have different kinetics
The above results establish the functional significance of NOS in a TH mutant background, but leave open the question of what role NO signaling plays in memory formation in wild-type flies. PPL1-γ1pedc induces a valence inverted memory when its activity precedes an odor such that the odor predicts relief from punishment (Aso and Rubin, 2016; Tanimoto et al., 2004). NOS is dispensable for this timing dependent inversion of memory valence (Figure 5—figure supplement 3), an observation consistent with a recent report that timing dependent inversion of long-term-depression of KC-to-MBON synapses to long-term-potentiation is mediated by different dopamine receptors (Handler et al., 2019). Next, we examined if NO-dependent memory requires, scribble (scrib), a gene encoding a scaffold protein which is important for forgetting; loss of scribble in KCs prolongs memory retention of odor-electric shock associative memories (Cervantes-Sandoval et al., 2016). RNAi knockdown of scrib impaired NO-dependent memory in TH-mutant flies (Figure 5D), suggesting a role of NO in regulating memory stability.
To further explore the interplay of dopamine and NO signaling on memory dynamics we examined the consequences of their combined action. First, we examined memory acquisition rates when flies were trained by activation of PPL1-γ1pedc (Figure 6A). In wild-type flies, PPL1-γ1pedc activation as brief as 10 s can induce significant negative-valence memory. Blocking NOS activity with L-NNA did not affect the memory scores observed shortly after a wide range of training protocols (Figure 6B). In contrast, longer and repetitive training was required to induce robust NO-dependent positive-valence memory in the absence of dopamine (Figure 6A–B).
Riemensperger et al. (2011) reported that dopamine-deficient flies developed weak positive-valence memory after odor-shock conditioning, but this memory was not detectable until 2 hr after the training. Motivated by this observation, we examined the kinetics of NO-dependent memory formation and the role of NOS in memory retention. When we used a single cycle of training, we found that NO-dependent memory develops slowly over time. Memory scores were not significantly different from zero at 1, 3, 5 min after training, and only became significant after 10 min. Once formed, however, these NO-dependent memories were long lasting and were still more than half maximal after 6 hr (Figure 6C). This result contrasts with the time course of memory formation by PPL1-γ1pedc activation in a wild-type background, where memory is detectable within 1 min after training but has a half-life of only 2–3 hr (Aso et al., 2012). These observations raised the possibility that NO signaling, with its opposite valence and slower dynamics, might serve to limit memory retention in a wild-type background. Indeed, we found that expression of NOS-RNAi or L-NNA feeding prolonged the retention of memories induced by either optogenetic training with PPL1-γ1pedc or odor-shock conditioning (Figure 6D, Figure 6—figure supplement 1A–B).
Memory persistence is often enhanced by repetitive training. However, PPL1-γ1pedc fails to induce long-lasting memory even after 10x repetitive training at spaced intervals (Aso and Rubin, 2016). In contrast, other DANs from the PPL1 cluster that do not exhibit significant NOS expression, PPL1-α3 or a combination of PPL1-γ2α'1 and PPL1-α'2α2, are able to induce stable memory lasting for 4 days after 10x spaced training (Aso and Rubin, 2016). Our results strongly suggest that NO signaling is responsible for this difference in memory retention. Spaced training with PPL1-γ1pedc induced memory lasting 1 day when NOS signaling was compromised, either by knockdown with RNAi (Figure 6E) or inhibition by L-NNA (Figure 6—figure supplement 1C). The valence-inverted memories formed in the γ1pedc compartment following repetitive training in TH mutants also lasted 1 day after training (Figure 6E). Thus, the effects of NOS signaling accumulate slowly, but can be long lasting. These effects are antagonistic to memories formed by dopamine signaling, and serve to sculpt the time course of memory retention. As discussed below, NO-signaling also contributes to other features of memory dynamics.
Nitric oxide enhances fast updating of memory
We designed behavioral experiments to examine memory dynamics when flies that had been previously trained encounter a new experience. We tested three different types of new experience: (1) switching which odor is paired with DAN activation during odor conditioning (reversal conditioning); (2) exposing flies to odors without DAN activation; and (3) activating DANs without odor exposure (Figure 7A). In wild-type flies, odor preference can be altered by a single trial of reversal conditioning with PPL1-γ1pedc activation (Figure 7B, left), whereas this process became slower and required more training when NOS was inhibited (Figure 7B; center), such that switching odor preference required a repetition of reversal conditioning. In TH mutants, NO-dependent memory was also altered by reversal conditioning, but with an even slower time course (Figure 7B; right). The second type of new experience, exposure to odors alone, did not change the existing memory in this assay (Figure 7C). The third type of new experience, DAN activation alone, quickly reduced conditioned odor response in wild-type flies. Inhibiting NOS slowed this process (Figure 7D; center). NO-dependent memory was also reduced by unpaired activation of DANs, but it took five trials to detect significant reduction (Figure 7D; right). These results suggest that both the slow formation and the persistence of NO-dependent memory facilitate the fast updating of memories stored in NOS-positive MB compartments in response to changing conditions.
Modeling the function of NO and dopamine in memory dynamics
To understand the interplay between dopamine (DA) and NO-dependent plasticity, we fit a minimal model to our data that accounts for the observed effects of these two pathways on the formation of odor memories. We then used this model both to infer the parameters of a synaptic plasticity rule consistent with the data as well as to test hypotheses about the mechanisms of DA and NO-mediated synaptic modifications that would be able to generate the memory dynamics we observed. Imaging and physiology experiments have demonstrated that DANs induce intracellular signaling cascades in KCs and bidirectionally modulate the weights of KC-MBON synapses, with the direction of the plasticity determined by whether each KC is active or inactive (Cohn et al., 2015; Gervasi et al., 2010; Hattori et al., 2017; Hige et al., 2015; Berry et al., 2018; Boto et al., 2014; Tomchik and Davis, 2009; Owald et al., 2015; Bouzaiane et al., 2015; Cassenaer and Laurent, 2012; Handler et al., 2019; Okada et al., 2007). This dependence on presynaptic KC firing ensures the odor-specificity of memories formed following DAN activation. NO-dependent memories, like DA-dependent memories, are odor-specific, suggesting that plasticity induced by NO also depends on KC activity. Based on this observation and the results of previous studies, we constructed our model by assuming that: (1) both DA and NO bidirectionally modulate KC-MBON synapses depending on KC activity, (2) memory decay is due to background levels of DAN activity following conditioning (Plaçais et al., 2012; Berry et al., 2018; Sitaraman et al., 2015), and (3) the effects of DA and NO occur via independent pathways and can coexist (Figure 6).
Specifically, we denoted the effects of DA and NO-dependent synaptic plasticity at time t by two quantities, D(t) and N(t), that lie between 0 and 1. We assumed that coincident KC-DAN activation increases D and N with a timescale of 30 s and 10 min, respectively, to account for the slower induction of NO-mediated effects (Figure 8A,B). Based on previous observations in TH wild-type flies that pairing of activation of the PPL1-γ1pedc DAN with odor induces synaptic depression between odor activated KCs and MBON-γ1pedc>α/β (Hige et al., 2015), we assumed that the effect of an increase in D is a reduction in the weight of the corresponding KC-MBON synapse (Figure 8B, left). As the activity of this MBON promotes approach behavior (Aso et al., 2014b; Owald et al., 2015), its reduced response to the conditioned odor after DA-dependent synaptic depression results in avoidance during subsequent odor presentations. In the TH-null background, in contrast, we have shown that PPL1-γ1pedc activation leads to a positive-valence memory (Figure 2). This is most readily explained in our model by assuming that, in such flies, NO induces potentiation of synapses between odor-activated KCs and MBONs (Figure 8B, right). Thus, in the model, the effect of an increase in N is an increase in synaptic weight, opposite to the effect of D. Finally, based on observations that activation of DANs alone can reverse synaptic depression induced by KC-DAN pairing (Cohn et al., 2015; Hattori et al., 2017; Berry et al., 2018), we assumed that DAN activation in the absence of KC activity causes a reduction in D and N, recovering the synapse to its baseline weight.
We fit the model by assuming that the performance index (PI) is determined by the odor-evoked activation of the MBON and then determining the parameters that best match the behavioral data reported in Figure 6B–D. We used data that isolates the effects of DA and NO-dependent plasticity mechanisms to fit the parameters for the two pathways separately (see Materials and methods). In the resulting model, NO-dependent plasticity develops more gradually and requires more KC-DAN pairings to produce a memory of equal magnitude, compared to DA-dependent plasticity (Figure 8C).
We next asked how these plasticity mechanisms interact to determine effective KC-to-MBON synaptic weights when both DA and NO pathways are active. The synaptic weight is a function of both DA and NO-mediated effects, w = f(D,N). One possibility for the function f is a difference between terms corresponding to DA-dependent depression and NO-dependent facilitation; that is, w ∝ N – D. When we fit a model with this functional form to our data, we found that it incorrectly predicts a reduction in memory strength after repeated pairings, because of the slower accumulation of NO-dependent facilitation after DA-dependent potentiation has saturated (Figure 8D, gray curve; Figure 8C). Another possibility for f is a multiplicative form, for example w ∝ (1 - D) × (1 + N). While we cannot unambiguously determine the identity of the biophysical parameters underlying DA and NO-mediated effects, such a form would arise naturally if the two pathways modulated parameters such as quantal size and the probability of synaptic vesicle release from KCs. We found that the multiplicative model provides a good match to our data (Figure 8C, blue curve). In this model, strong DA-dependent depression (i.e. D close to 1) leads to a small synaptic weight even in the presence of NO-dependent facilitation.
We further tested our model by investigating its behavior in other paradigms. Assuming that spontaneous activity in DANs leads to memory decay (Plaçais et al., 2012; Berry et al., 2018; Sitaraman et al., 2015) accounted for the NO-dependent reduction in memory lifetime that we observed (Figure 8E). Fitting the magnitudes of DA-dependent facilitation and NO-dependent depression in the absence of KC activity using the data of Figure 7D also predicted the dynamics of reversal learning and its facilitation by NO with no additional free parameters (Figure 8F; Figure 8—figure supplement 1). On the other hand, the additive model failed to accurately predict these dynamics (Figure 8—figure supplement 2). In total, modeling a multiplicative interaction between DA- and NO-dependent plasticity accounts for the immediate effects of these pathways on odor memories that we observed experimentally. A notable exception is that this model cannot account for the enhanced persistence of memories after 10x training for DA or NO-null conditions (Figure 6E), suggesting a recruitment of additional consolidation mechanisms after spaced conditioning, as previously proposed (Pagani et al., 2009; Tully et al., 1994; Pai et al., 2013; Scheunemann et al., 2018; Cervantes-Sandoval et al., 2013; Miyashita et al., 2018; Huang et al., 2012; Akalal et al., 2011). Also, we found that spaced training with PPL-α3 can induce long-lasting memory but it requires at least 3 hr after training to develop (Figure 9A). This observation is consistent with the report that protein synthesis in MBON-α3 in the first 3 hr time window after odor-shock spaced training is required for LTM formation (Pai et al., 2013; Wu et al., 2017), but the current model does not account for the postsynaptic parameters.
Discussion
Evidence from a wide range of organisms establishes that dopaminergic neurons often release a second neurotransmitter, but the role of such cotransmitters in diversifying neuronal signaling is much less clear. In rodents, subsets of dopaminergic neurons co-release glutamate or GABA (Maher and Westbrook, 2008; Stuber et al., 2010; Sulzer et al., 1998; Tecuapetla et al., 2010; Tritsch et al., 2012). In mice and Drosophila, single-cell expression profiling reveals expression of diverse neuropeptides in dopaminergic neurons (Croset et al., 2018; Poulin et al., 2014). EM connectome studies of the mushroom body in adult and larval Drosophila reveal the co-existence of small-clear-core and large-dense-core synaptic vesicles in individual terminals of dopaminergic neurons (Eichler et al., 2017; Takemura et al., 2017); moreover, the size of the observed large-dense-core vesicles differs between DAN cell types (Takemura et al., 2017).
We found that NOS, the enzyme that synthesizes NO, was located in the terminals of a subset of DAN cell types. NOS catalyzes the production of nitric oxide (NO) from L-arginine. Drosophila NOS is regulated by Ca2+/calmodulin (Regulski and Tully, 1995), as is the neuronal isoform of NOS in the mammalian brain (Abu-Soud and Stuehr, 1993), raising the possibility that NO synthesis might be activity dependent. Furthermore, the localization of the NOS1 protein in the axonal terminals of DANs (Figure 3D) is consistent with NO serving as a cotransmitter. Our conclusion that NO acts as a neurotransmitter is supported by the observation that NO signaling requires the presence of a putative receptor, soluble guanylate cyclase, in the postsynaptic Kenyon cells. This role contrasts with the proposed cell-autonomous action of NOS in the ellipsoid body, in which NO appears to target proteins within the NOS-expressing ring neurons themselves, rather than conveying a signal to neighboring cells (Kuntz et al., 2017).
The valence-inversion phenotype we observed when PPL1-γ1pedc was optogenetically activated in a dopamine-deficient background can be most easily explained if NO induces synaptic potentiation between odor-activated KCs and their target MBONs. Our modeling work is consistent with this idea, but testing this idea and other possible mechanisms for NO action will require physiological experiments.
Antagonistic functions of dopamine and nitric oxide
During olfactory learning, the concentration of Ca2+ in KC axons represents olfactory information. The coincidence of a Ca2+ rise in spiking KCs and activation of the G-protein-coupled Dop1R1 dopamine receptor increases adenylyl cyclase activity (Abrams et al., 1998; Boto et al., 2014; Byrne et al., 1991; Kim et al., 2007; McGuire et al., 2003; Tomchik and Davis, 2009). The resultant cAMP in turn activates protein kinase A (Gervasi et al., 2010; Skoulakis et al., 1993), a signaling cascade that is important for synaptic plasticity and memory formation throughout the animal phyla (e.g. Davis et al., 1995). In contrast, when DANs are activated without KC activity, and thus during low intracellular Ca2+ in the KCs, molecular pathways involving the Dop1R2 receptor, Rac1 and Scribble facilitate decay of memory (Berry et al., 2012; Cervantes-Sandoval et al., 2016; Kim et al., 2007; Shuai et al., 2010).
We found that NOS in PPL1-γ1pedc shortens memory retention, while facilitating fast updating of memories in response to new experiences. These observations could be interpreted as indicating that NO regulates forgetting. Indeed, NO-dependent effect requires scribble in KCs, a gene previously reported as a component of active forgetting (Figure 5D). However, it is an open question whether the signaling pathways for forgetting, which presumably induce recovery from synaptic depression (Berry et al., 2018; Cohn et al., 2015), are related to signaling cascades downstream of NO and guanylate cyclase, which appear to be able to induce memory without prior induction of synaptic depression by dopamine. Lack of detectable 1-day memory formation after spaced training with PPL1-γ1pedc can be viewed as a balance between two distinct, parallel biochemical signals, one induced by dopamine and the other by NO (Figure 6E), rather than the loss of information (that is, forgetting). Confirming this interpretation will require better understanding of the signaling pathways downstream of dopamine and NO. The search for such pathways will be informed by the prediction from our modeling that dopamine and NO may alter two independent parameters that define synaptic weights with a multiplicative interaction.
In the vertebrate cerebellum, which has many architectural similarities to the MB (Farris, 2011; Marr, 1969; Medina et al., 2002), long-term-depression at parallel fiber-Purkinje cell synapses (equivalent to KC-MBON synapses) induced by climbing fibers (equivalent to DANs) can coexist with long-term-potentiation by NO (Bredt et al., 1990; Lev-Ram et al., 2002; Shibuki and Okada, 1991). In this case, the unaltered net synaptic weight results from a balance between coexisting LTD and LTP rather than recovery from LTD. This balance was suggested to play an important role in preventing memory saturation in the cerebellum and allowing reversal of motor learning. In the Drosophila MB, we observed a similar facilitation of reversal learning by NO (Figure 7B). The antagonistic roles of NO and synaptic depression may be a yet another common feature of the mushroom body and the cerebellum.
Distinct dynamics of dopamine and nitric oxide
Opposing cotransmitters have been observed widely in both invertebrate and vertebrate neurons (Nusbaum et al., 2017). A common feature in these cases is that the transmitters have distinct time courses of action. For instance, hypothalamic hypocretin-dynorphin neurons that are critical for sleep and arousal synthesize excitatory hypocretin and inhibitory dynorphin. When they are released together repeatedly, the distinct kinetics of their receptors result in an initial outward current, then little current, and then an inward current in the postsynaptic cells (Li and van den Pol, 2006). In line with these observations, we found that dopamine and NO show distinct temporal dynamics: NO-dependent memory requires repetitive training and takes longer to develop than dopamine-dependent memory. What molecular mechanisms underlie these differences? Activation of NOS may require stronger or more prolonged DAN activation than does dopamine release. Alternatively, efficient induction of the signaling cascade in the postsynaptic KCs might require repetitive waves of NO input. Direct measurements of release of dopamine and NO, and downstream signaling events by novel sensors will be needed to address these open questions (Chen et al., 2014b; Eroglu et al., 2016; Patriarchi et al., 2018; Sun et al., 2018; Tang and Yasuda, 2017).
Toward subcellular functional mapping of memory genes
Decades of behavioral genetic studies have identified more than one hundred genes underlying olfactory conditioning in Drosophila (Keene and Waddell, 2007; McGuire et al., 2005; Thurmond et al., 2019; Walkinshaw et al., 2015). Mutant and targeted rescue studies have been used to map the function of many memory-related genes encoding synaptic or intracellular signaling proteins (for example, rutabaga, DopR1/dumb, DopR2/DAMB, PKA-C1/DC0, Synapsin, Bruchpilot, Orb2 and Rac1) to specific subsets of Kenyon cells (Akalal et al., 2006; Berry et al., 2012; Blum et al., 2009; Gervasi et al., 2010; Han et al., 1996; Kim et al., 2007; Knapek et al., 2011; Krüttner et al., 2015; McGuire et al., 2003; Niewalda et al., 2015; Pavot et al., 2015; Qin et al., 2012; Shuai et al., 2010; Skoulakis et al., 1993; Trannoy et al., 2011; Vanover et al., 2019; Zars et al., 2000). However, it is largely unknown if these proteins physically colocalize at the same KC synapses to form intracellular signaling cascades. Some of these proteins might preferentially localize to specific MB compartments. Alternatively, they may distribute uniformly along the axon of Kenyon cells, but be activated in only specific compartments. Our identification of cell-type-specific cotransmitters in DANs enabled us to begin to explore this question.
We used optogenetic activation of specific DANs to induce memory in specific MB compartments, while manipulating genes in specific types of KCs. This approach allowed us to map and characterize the function of memory-related genes at a subcellular level. For example, the Gycbeta100B gene, which encodes a subunit of sGC, has been identified as ‘memory suppressor gene’ that enhances memory retention when pan-neuronally knocked down (Walkinshaw et al., 2015), but the site of its action was unknown. Gycbeta100B appears to be broadly dispersed throughout KC axons, based on the observed distribution of a Gycbeta100B-EGFP fusion protein (Figure 5—figure supplement 1). Our experiments ectopically expressing NOS in PPL1-α3 DANs that do not normally signal with NO is most easily explained if sGC is available for activation in all MB compartments (Figure 4D–E, Figure 4—figure supplement 1).
What are the molecular pathways downstream to cGMP? How do dopamine and NO signaling pathways interact in regulation of KC synapses? Previous studies and RNA-Seq data suggest several points of possible crosstalk. In cultured KCs from cricket brains, cGMP-dependent protein kinase (PKG) mediates NO-induced augmentation of a Ca2+ channel current (Kosakai et al., 2015). However, we failed to detect expression of either of the genes encoding Drosophila PKGs (foraging and Pkg21D) in KCs in our RNA-Seq studies (Supplementary file 1). On the other hand, cyclic nucleotide-gated channels and the cGMP-specific phosphodiesterase Pde9 are expressed in KCs. Biochemical studies have shown that the activity of sGC is calcium dependent and that PKA can enhance the NO-induced activity of sGC by phosphorylating sGC; sGC isolated from flies mutant for adenylate cyclase, rutabaga, show lower activity than sGC from wild-type brains (Morton et al., 2005; Shah and Hyde, 1995), suggesting crosstalk between the cAMP and cGMP pathways.
The benefits of parallel memory units with heterogeneous dynamics
All memory systems must contend with a tension between the strength and longevity of the memories they form. The formation of a strong immediate memory interferes with and shortens the lifetimes of previously formed memories, and reducing this interference requires a reduction in initial memory strength that can only be overcome through repeated exposure (Amit and Fusi, 1994). Theoretical studies have argued that this tension can be resolved by memory systems that exhibit a heterogeneity of timescales, balancing the need for both fast, labile memory and slow consolidation of reliable memories (Fusi et al., 2005; Lahiri and Ganguli, 2013; Benna and Fusi, 2016). The mechanisms responsible for this heterogeneity, and whether they arise from complex signaling within synapses themselves (Benna and Fusi, 2016), heterogeneity across brain areas (Roxin and Fusi, 2013), or both, have not been identified.
We found that NO acts antagonistically to dopamine and reduces memory retention (Figure 6) while facilitating the rapid updating of memory following a new experience (Figure 7). Viewed in isolation, the NO-dependent reduction in memory retention within a single compartment may seem disadvantageous, but in the presence of parallel learning pathways, this shortened retention may represent a key mechanism for the generation of multiple memory timescales that are crucial for effective learning. During shock conditioning, for example, multiple DANs respond to the aversive stimulus, including PPL1-γ1pedc, PPL1-γ2α'1, PPL1-α3 (Mao and Davis, 2009; Riemensperger et al., 2005). We have shown that optogenetic activation of these DAN cell types individually induces negative-valence olfactory memory with distinct learning rates (Aso and Rubin, 2016). The NOS-expressing PPL1-γ1pedc induces memory with the fastest learning rate in a wild-type background, and we show here that it induces an NO-dependent memory trace when dopamine synthesis is blocked, with a much slower learning rate and opposite valence. Robust and stable NO-dependent effects were only observed when training was repeated 10 times (Figure 6E). Under such repeated training, compartments with slower learning rates, such as α3, form memory traces in parallel to those formed in γ1pedc (Pai et al., 2013; Séjourné et al., 2011; Wu et al., 2017). Thus, flies may benefit from the fast and labile memory formed in γ1pedc without suffering the potential disadvantages of shortened memory retention, as long-term memories are formed in parallel in other compartments (Figure 9B). The Drosophila MB provides a tractable experimental system to study the mechanisms and benefits of diversifying learning rate, retention, and flexibility in parallel memory units, as well as exploring how the outputs from such units are integrated to drive behavior.
Materials and methods
Flies
Drosophila strains used in this study are listed in the Supplementary file 2 KEY RESOURCES TABLE. Crosses listed in Supplementary file 3 were kept on standard cornmeal food supplemented with retinal (0.2 mM all-trans-retinal prior to eclosion and then 0.4 mM) at 21°C at 60% relative humidity in the dark. Female flies were sorted on cold plates at least 1 d prior to the experiments, and 4–10 day old flies were used for experiments. Additional drugs were administered by feeding with retinal containing fly food mixed with drugs. The L‐DOPA (D9628, Sigma) or L-NNA were mixed directly into melted fly food at final concentrations of 1 mg/ml or 1–100 mM, respectively. S‐(−)‐Carbidopa (C1335, Sigma) was dissolved in 1 ml of water at 10x the final concentration and mixed with 9 ml of melted fly food.
Olfactory learning assay
Request a detailed protocolGroups of approximately 20 female flies, 4–10 day post-eclosion, were trained and tested at 25°C at 50% relative humidity in the fully automated olfactory arena for optogenetics experiments (Aso and Rubin, 2016; Pettersson, 1970; Vet et al., 1983). The odors were diluted in paraffin oil (Sigma–Aldrich): 3-octanol (OCT; 1:1000; Merck) and 4-methylcyclohexanol (MCH; 1:750; Sigma–Aldrich). Videography was performed at 30 frames per second and analyzed using Fiji (Schindelin et al., 2012). Statistical comparisons were performed (Prism; Graphpad Inc, La Jolla, CA 92037) using the Kruskal Wallis test followed by Dunn's post-test for multiple comparison, except those in Figures 4C, E, 6B, C and E and E, 9A, Figure 4—figure supplement 1 and Figure 5—figure supplement 3C, which used Wilcoxon signed-rank test with Bonferroni correction to compare from zero. Appropriate sample size for olfactory learning experiment was estimated based on the standard deviation of performance index in previous study using the same assay (Aso and Rubin, 2016). We set the effect size, power and significance as 0.15, 0.8 and 0.05, respectively.
Immunohistochemistry
Request a detailed protocolDissection and immunohistochemistry of fly brains were carried out as previously described with minor modifications (Jenett et al., 2012) using the antibodies listed in KEY RESOURCES TABLE. Brains and VNCs of 3- to 10-day-old female flies were dissected in Schneider's insect medium and fixed in 2% paraformaldehyde in Schneider's medium for 55 min at room temperature (RT). After washing in PBT (0.5% Triton X-100 in PBS), tissues were blocked in 5% normal goat serum (or normal donkey serum, depending on the secondary antibody) for 90 min. Subsequently, tissues were incubated in primary antibodies diluted in 5% serum in PBT for 2–4 days on a nutator at 4°C, washed four times in PBT for 15 min or longer, then incubated in secondary antibodies diluted in 5% serum in PBT for 2–4 days on a Nutator at 4°C. Tissues were washed thoroughly in PBT four times for 15 min or longer and mounted on glass slides with DPX.
For immunolabeling of NOS, the serum against NOS exon 16 was obtained from N. Yakuobovich and P. H. O’Farrell (Yakubovich et al., 2010), and then affinity purified as described below. In order to minimize non-specific signals, we absorbed 200 µL of anti-NOS antibody (1:1000) for 1 day with 30 fly brains in which NOS was knocked down panneuronally using elav-GAL4 and NOS-RNAi strain#1 (TRiP.HMC03076), and the supernatant was used for subsequent immunohistochemistry.
Purification of dNOS proteins and antibody
Request a detailed protocolThe pRSET-dNOS exon 16 construct, containing an N-terminal His tag and T7 gene 10 leader RBS site, was assembled as follows. The NEBuilder Assembly tool was used to design primers for the NEBuilder HiFi DNA assembly (New England Biolabs # E2621S) of dNOS exon 16 as contained in pET28a-dNOS exon 16 (gift of Nikita Yakubovich, O-Farrell lab, UCSF) into the backbone vector pRSET (gift of Ariana Tkachuk, Janelia) which was digested with BamHI/EcoRI. The assembled product, pRSET-dNOS exon 16, was first transformed into NEB 5-alpha competent cells (New England Biolabs #E2621S) and plated on LB plus ampicillin (60 µl/ml).
For protein purification, pRSET-dNOS exon 16 DNA was then transformed into T7 Express lysY/lq E. coli protein expression cells (New England Biolabs #C3013) and plated on LB plus ampicillin (60 µl/ml). For large-scale growth, 500 ml of Miller’s LB plus ampicillin (60 µl/ml) was inoculated from 5 ml of a starter culture and grown for ~3 hr (~O.D. 0.5–0.7) at 37°C and then induced by adding 0.5 mM IPTG. The culture was allowed to grow at 18–20°C overnight before spinning down and freezing the recovered pellets which were divided in two 250 ml bottles.
To resuspend the thawed pellets (frozen overnight), 10–12 ml of the nonionic detergent- based lysis reagent B-PER (Thermo Scientific #78266) in phosphate buffer containing 1 mg/ml of lysozyme (Thermo Scientific #89833), nuclease (0.1 µl/ml, Thermo Scientific # 88701), and 1X HALT protease inhibitors (Thermo Scientific #1861279). The suspension, divided into two 50 ml conical tubes, was gently shaken for 15–20 min at 30°C before spinning down at 8000 x g for 15 mins at 4°C. We found that that the majority of the dNOS exon 16 protein was in inclusion bodies and therefore we carried out purification starting with the pellet.
Each pellet was resuspended in ~10 ml B-PER containing 200 µg/ml lysozyme. The suspension was then mixed with 100 ml of a wash buffer containing 1:10 B-PER in ice cold 1xPBS (diluted from 10X PBS stock, Fisher Scientific #BP3994) by pipetting the mixture up and down and gentle agitation. The mixture was spun down at 15,000 rpm for 15 min at 4°C. The pellet was similarly washed four more times. The washed pellet was either stored at −20°C overnight or resuspended in 7–12 ml of inclusion body solubilization reagent (Thermo Scientific #78115). The protein suspension was shaken for 30–40 min at 20°C and then ultracentrifuged at 35,000 x g for 20–30 min at 4°C.
For affinity purification of anti-dNOS antibody, the supernatant fraction was concentrated using 50 ml conical tubes Vivaspin 20, 10,000 MWCO concentrators (Sartorius # VS2002) and then dialyzed in 3–12 ml dialysis cassettes 10,000 MWCO (Thermo Scientific #66810) against 1 L of 4 M guanidine HCL (diluted from 6 M stock, Sigma #SRE0066) in 1X PBS pH 8.0 for ~6 hr at 4°C. The medium was further diluted to 2 M guanidine HCL and the protein continued to dialyze overnight. Affinity purification of the protein-antibody complex was performed using the AminoLink Immobilization Kit (Thermo Scientific #44890). Approximately 6 mg of soluble dNOS exon 16 protein in 2 M guanidine in 1X PBS was bound to the agarose beads in the column as antigen, and 1.9 ml of crude rabbit anti sera to dNOS exon 16 (gift of Nikita Yakubovich and Patrick O-Farrell) was run through the column. The purified dNOS exon 16 rabbit anti-antibody fractions were eluted with IgG Elution buffer (Thermo Scientific # 21004) and then concentrated with Vivaspin 20 tube concentrators before dialyzed in 1X PBS at 4°C for 2.5 days with one change of fresh 1X PBS.
Tissue expansion
Request a detailed protocolTissue expansion was performed as described in Tillberg et al. (2016). All solutions were prepared in milliQ-grade water unless otherwise specified. AcX stock: acryloyl-X, SE (ThermoFisher, A20770) at 10 mg/mL in anhydrous DMSO. PLL solution: poly-l-lysine solution (Ted Pella 18026) with Photo-Flo detergent (EMS 74257) added 1:500. Acrylate stock: 4 M, prepared by neutralizing 5.5 mL acrylic acid (99% purity; Sigma, 147230) with 10 N NaOH using a water bath and fume hood, in a total volume of 20 mL. Acrylamide stock: 50% (w/v) (Sigma, A9099). Bisacrylamide stock: 1% (w/v) (Sigma, M7279). Monomer stock: 11.5 mL sodium acrylate stock, 2.5 mL acrylamide stock, 7.5 mL bisacrylamide stock, 18 mL 5 M NaCl (Sigma, S5150-1L), 5 mL 10xPBS (ThermoFisher, 70011044), and 2.5 mL water for a total volume of 47 mL. 4-HT stock: 4-hydroxy TEMPO at 0.5% (w/v) (Sigma, 176141). TEMED stock: N,N,N′,N′-Tetramethylethylenediamine at 10% (v/v) (Sigma, T7024). APS stock: ammonium persulfate at 10% (w/v) (Sigma, A3678). ProK digestion buffer: 0.5% Triton X-100, 500 mM NaCl, 1 mM EDTA, 50 mM Tris pH8. Appropriate caution was exercised when handling acrylamide, a known toxin.
Dissected, fixed, and antibody-stained Drosophila brains were treated with AcX stock solution diluted 1:100 in 1xPBS, with shaking, overnight. Brains were then washed with 1xPBS. A gelation chamber was created by applying a Press-to-Seal silicone gasket (ThermoFisher, P24740) to a glass slide, which was then coated with PLL solution. AcX-treated brains were immobilized on the PLL surface, up to nine per gasket. Gelation solution was created by adding 10 µL each of 4HT, TEMED, and APS stock solutions to 470 µL of monomer stock solution on ice. Brains were washed with gelation solution and then the gelation chamber was filled with ~200 µL of gelation solution and incubated on ice for 25 min. The gelation chamber was then sealed with a cover slip and placed in a 37°C incubator to gel and cure for 2 hr.
Gelation chambers were disassembled and individual gels trimmed close to each brain. Gels were trimmed to a right trapezoid shape to facilitate keeping track of specimen orientation. Gels were incubated with proteinase K digestion enzyme (NEB, P8107S) diluted 1:100 in proK digestion buffer with shaking, overnight. Digested gels were washed with water 4 × 30 min followed by equilibration overnight.
Image acquisition and analysis
Request a detailed protocolWe used a LSM710 confocal microscope (Zeiss; 20x/0.8 M27 or 63x/1.40 oil immersion objective) for imaging brains and a custom-made lattice light sheet microscope for imaging expanded brains as previously described (Chen et al., 2014a; Gao et al., 2019). Images were analyzed using Fiji (http://fiji.sc/), and visualized with VVD_Viewer (https://github.com/takashi310/VVD_Viewer/blob/master/README.md), a modified version of Fluorender (http://www.sci.utah.edu/software/13-software/127-fluorender.html; Wan et al., 2012).
Fluorescence in situ hybridization (FISH)
Request a detailed protocolFISH probe libraries were designed based on transcript sequences and were purchased from Biosearch Technologies. The FISH protocol and dye labeling procedures were described previously (Long et al., 2017). FISH probes for detecting tyrosine hydroxylase transcripts were described (Meissner et al., 2019). FISH probes for NOS are listed in KEY RESOURCES TABLE. Each probe contains a 3’-end amine-modified nucleotide that allows direct coupling to an NHS-ester Cy3 dye (GE Healthcare, PA23001) according to the manufacturer’s instructions. The brains of 3–5 day old adult flies were dissected in 1xPBS and fixed in 2% paraformaldehyde diluted PBS at room temperature for 55 min. Brain tissues were washed in 0.5% PBT, dehydrated, and stored in 100% ethanol at 4°C. After exposure to 5% acetic acid at 4°C for 5 min, the tissues were fixed in 2% paraformaldehyde in 1xPBS for 55 min at 25°C. The tissues were then washed in 1 × PBS with 1% of NaBH4 at 4°C for 30 min. Following a 2 hr incubation in prehybridization buffer (15% formamide, 2 × SSC, 0.1% Triton X-100) at 50°C, the brains were introduced to hybridization buffer (10% formamide, 2x SSC, 5x Denhardt's solution, 1 mg/ml yeast tRNA, 100 μg/ml, salmon sperm DNA, 0.1% SDS) containing FISH probes at 50°C for 10 hr and then at 37°C for an additional 10 hr. After a series of wash steps, the brains were dehydrated, cleared in xylene, and mounted in DPX. Image Z-stacks were collected using an LSM880 confocal microscope fitted with an LD LCI Plan-Apochromat 25x/0.8 oil or Plan-Apochromat 63x/1.4 oil objective after the tissue cured for 24 hr.
Calcium imaging
Request a detailed protocolFlies were reared at 25°C on cornmeal medium supplemented with retinal (0.2 mM) that was shielded from light. Flies were transferred to cornmeal media mixed with 100 mM L-NG-Nitroargenine (L-NNA) 0.4 mM retinal while siblings used for comparison were transferred to the same media without L-NNA 24 hr before testing. All experiments were performed on female flies, 3–5 days post-eclosion. Brains were dissected in a saline bath (103 mM NaCl, 3 mM KCl, 2 mM CaCl2, 4 mM MgCl2, 26 mM NaHCO3, 1 mM NaH2PO4, 8 mM trehalose, 10 mM glucose, 5 mM TES, bubbled with 95% O2/5% CO2). After dissection, the brain was positioned anterior side up on a coverslip in a Sylgard dish submerged in 2 ml saline at 20°C.
The sample was imaged with a resonant scanning 2-photon microscope with near-infrared excitation (920 nm, Spectra-Physics, INSIGHT DS DUAL) and a 25X objective (Nikon MRD77225 25XW). The microscope was controlled by using ScanImage 2015.v3 (Vidrio Technologies)15. Volumes were acquired with 141 μm × 141 μm field of view at 512 × 512 pixel resolution at 2 μm steps over 42 slices, at approximately 1 Hz. The excitation power for Ca2+ imaging measurement was 12 mW.
For photostimulation, the light-gated ion channel CsChrimson was activated with a 660 nm LED (M660L3 Thorlabs) coupled to a digital micromirror device (Texas Instruments DLPC300 Light Crafter) and combined with the imaging path with a FF757-DiO1 dichroic (Semrock). On the emission side, the primary dichroic was Di02-R635 (Semrock), the detection arm dichroic was 565DCXR (Chroma), and the emission filters were FF03-525/50 and FF01-625/90 (Semrock). Photostimulation from PPL1-γ1pedc to MBON-γ1pedc occurred 5 times at 30 s intervals at maximum intensity of 1.2 mW/mm2 as show Figure 5—figure supplement 1. When testing the PPL1-γ1pedc to the PAM-DAN connection there were 9 stimulations at 30 s intervals with an intensity of 0.89 mW/mm2, as measured using Thorlabs S170C power sensor. The stimulation duration was 1 s in all trials.
After testing, laser power was increased to take two color high-resolution images containing fluorescence from both red and green channels. Using custom python scripts and features in the images, regions of interest (ROIs) were draw corresponding to the mushroom body compartments. Fluorescence in a background ROI, that contained no endogenous fluorescence, was subtracted from the mushroom body ROIs. In the ΔF/F calculations, baseline fluorescence is the median fluorescence over a 5-s time period before stimulation started. The ΔF is the fluorescence minus the baseline. Then the ΔF is divided by baseline to normalize the signal (ΔF/F). The final signal was run through a gaussian filter (sigma = 1). When more than one trial was performed, the mean ΔF/F was calculated over the trials and then the mean between animals.
RNA-Seq
Expression checks
Request a detailed protocolNeurons of interest were isolated by expressing a fluorescent protein, either mCD8-GFP or tdTomato, using split-Gal4 drivers specific for particular cell types and then manually picking the fluorescent neurons from dissociated brain tissue. As a preliminary to the sorting process, each driver/reporter combination was ‘expression checked’ to determine if the marked cells were sufficiently bright to be sorted effectively and if there was any off-target expression in neurons other than those of interest. Drivers that met both these requirements were used in sorting experiments as described below.
Sorting of fluorescent-labeled neurons
Request a detailed protocolDrosophila adults were collected daily as they eclosed, and aged 3–5 days prior to dissection. For each sample, 60–100 brains were dissected in freshly prepared, ice cold Adult Hemolymph Solution (AHS; 108 mM NaCl, 5 mM KCl, 2 mM CaCl2, 8.2 mM MgCl2, 4 mM NaHCO3, 1 mM NaH2PO4, 5 mM HEPES, 6 mM Trehalose, 10 mM Sucrose), and the major tracheal branches removed. The brains were transferred to an 1.5 ml Eppendorf tube containing 500 microliters 1 mg/ml Liberase DH (Roche, prepared according to the manufacturer’s recommendation) in AHS, and digested for 1 hr at room temperature. The Liberase solution was removed and the brains washed three times with ice cold AHS. The final wash was removed completely and 400 μl of AHS+2% Fetal Bovine Serum (FBS, Sigma) were added. The brain samples were gently triturated with a series of fire-polished, FBS-coated Pasteur pipettes of descending pore sizes until the tissue was homogenized, after which the tube was allowed to stand for 2–3 min so that the larger debris could settle.
For hand sorting, the cell suspension was transferred to a Sylgard-lined Glass Bottom Dish (Willco Wells), leaving the debris at the bottom of the Eppendorf tube, and distributed evenly in a rectangular area in the center of the plate with the pipet tip. The cells were allowed to settle for 10–30 min prior to picking. Fluorescent cells were picked with a mouth aspirator consisting of a 0.8 mm Nalgene Syringe Filter (Thermo), a short stretch of tubing, a plastic needle holder, and a pulled Kwik-Fil Borosilicate Glass capillary (Fisher). Cells picked off the primary plate were transferred to a Sylgard-lined 35 mm Mat Tek Glass Bottom Microwell Dishes (Mat Tek) filled with 170 microliters AHS+2%FBS, allowed to settle, and then re-picked. Three washes were performed in this way before the purified cells were picked and transferred into 50 μl buffer XB from the PicoPure RNA Isolation Kit (Life Technologies), lysed for 5 m at 42°C, and stored at −80°C.
For FACS sorting, the cell suspension was passed through a Falcon 5 ml round-bottom tube fitted with a 35 micrometer cell strainer cap (Fisher), and sorted on a Becton Dickson FACSAria II cell sorter, gated for single cells with a fluorescence intensity exceeding that of a non-fluorescent control. Positive events were sorted directly into 50 microliters PicoPure XB buffer, the sample lysed for 5 m at 42°C, and stored at −80°C.
Library preparation and sequencing
Request a detailed protocolTotal RNA was extracted from 100 to 500 pooled cells using the PicoPure kit (Life Technologies) according to the manufacturer’s recommendation, including the on-column DNAse step. The extracted RNA was converted to cDNA and amplified with the Ovation RNA-Seq System V2 (NuGEN), and the yield quantified by NanoDrop (Thermo). The cDNA was fragmented and the sequencing adaptors ligated onto the fragments using the Ovation Rapid Library System (NuGEN). Library quality and concentration was determined with the Kapa Illumina Library Quantification Kit (KK4854, Kapa Biosystems), and the libraries were pooled and sequenced on an Illumina NextSeq 550 with 75 base pair reads. Sequencing adapters were trimmed from the reads with Cutadapt (Martin, 2011) prior to alignment with the STARsolo aligner (https://github.com/alexdobin/STAR/blob/master/docs/STARsolo.md) (Dobin et al., 2013; Dobin, 2020) to the Drosophila r6.34 genome assembly on Flybase (Thurmond et al., 2019). The resulting transcript alignments were passed to RSEM (Li and Dewey, 2011) to generate gene expression counts. The data set was deposited to NCBI Gene Expression Omnibus (accession number GSE139889).
Modeling
Inferring parameters of DA and NO plasticity
Request a detailed protocolWe assume that the immediate effects of DA-mediated and NO-mediated changes at a KC-to-MBON synapse are described by two variables d(t) and n(t), respectively. When the KC and corresponding DAN are coactive, these variables are modified according to:
When the DAN is active but the KC is inactive,
A and B determine how quickly the variables approach their maximum value of 1 or minimum value of 0.
To model the time it takes for the effects of synaptic plasticity to occur, we also define two additional variables
D(t) and N(t) follow the values of d(t) and n(t), but with slower timescales τDand τNrespectively. Based on the data, we assume τD = 30 s, and τN = 10 min. This accounts for the slower induction of NO-mediated effects.
We start by inferring AD and AN from data by relating the values of D(t) and N(t) in our model to the activation of the MBON, and finally to the PI. We assume that the normalized KC-to-MBON synaptic weight is given by w(t) ∝ 1 − D(t) (DA-mediated depression) if NO is absent, or w(t) ∝ 1 + N(t) (NO-dependent facilitation) if DA is absent. In our model, odors A and B activate a random 10% of KCs (the results do not depend on the total number of KCs NKC), and the activation of the MBON is given by , where si = 1 if the ith KC is active and 0 otherwise. At the beginning of each trial, D(t) = N(t) = 0, so wi = 1. If rA and rB are the MBON activations for odors A and B, then we assume the probability of the fly choosing odor A is equal to a softmax function of this activation:
We also infer g, which determines how strongly the MBON activation influences the decision, with g = 0 corresponding to random choices.
We infer AD, AN, and g separately for DA-null and NO-null conditions. To do so, we determine the values of the parameters that minimize the mean squared distance between model prediction and experimentally measured PIs for the 1 × 10 s, 1 × 1 min, 3 × 1 min, and 10 × 1 min protocols (Figure 6B). This is accomplished by simulating the model defined by Equations Equation 1 to Equation 3 and calculating the resulting preference index using Equation Equation 4. Optimal values for AD, AN, and g are found using a grid search. The optimization leads to AD = 4.3 min−1 and AN = 0.96 min−1 (the inferred values of g are similar for the two cases; 14.8 and 12.6, respectively). These values indicate that the DA effect saturates more quickly than the NO effect, nearly reaching its maximum effect after 1 × 1 min pairing.
Next, we model synaptic weights when both DA and NO-dependent changes occur. We consider two forms of interaction, an additive one with w(t) ∝ N(t) − D(t) and a multiplicative one with w(t) ∝ (1 − D(t))(1 + N(t)). We use the values of AD and AN inferred previously but allow g to be readjusted to best match the data (for the multiplicative model, g = 14.3, similar to above, while for the additive model g = 24.6). Only the multiplicative model qualitatively matches the experimental data (Figure 8D).
Modeling memory decay after 1 × 1 min pairing
Request a detailed protocolNext, we ask how memory decays after a pairing protocol. We assume that, after the pairing protocol is complete, there is a background level of activity in the DANs, which leads to depression according to Eq. 2. We infer the values of and (where the superscript denotes decay due to background DAN activity) using a grid search to minimize the mean squared distance between the predicted and actual PI from Figure 6D. Other parameters are set to the previously determined values using the multiplicative model above. This leads to = 2.7 x10−3 min−1) and = 1.6 x10−3 min−1).
Predicting memory dynamics in reversal learning
Request a detailed protocolFinally, we also model the effects of different pairing protocols. We start by considering DAN activation in the absence of odor. We infer values for BD and BN (different from the background levels above, since DANs are now activated rather than at their background levels of activation) by minimizing the mean squared difference between predicted and actual PI for NO-null and DA-null conditions (Figure 7F), and with the remaining parameters determined previously for the multiplicative model. This yields BD = 0.26 (min−1) and BN = 0.16 (min−1). These parameters are used to predict the behavior during reversal learning (Figure 7F).
Data availability
Complete transcript data were deposited to NCBI Gene Expression Omnibus (accession number GSE139889).
-
NCBI Gene Expression OmnibusID GSE139889. Bulk RNA-seq data from the dopaminergic neurons, MB output neurons, Kenyon cells, APL and DPM neurrons in adult Drosophila.
References
-
Roles for Drosophila mushroom body neurons in olfactory learning and memoryLearning & Memory 13:659–668.https://doi.org/10.1101/lm.221206
-
Learning in neural networks with material synapsesNeural Computation 6:957–982.https://doi.org/10.1162/neco.1994.6.5.957
-
Specific dopaminergic neurons for the formation of labile aversive memoryCurrent Biology 20:1445–1451.https://doi.org/10.1016/j.cub.2010.06.048
-
Computational principles of synaptic memory consolidationNature Neuroscience 19:1697–1706.https://doi.org/10.1038/nn.4401
-
Neural and molecular bases of nonassociative and associative learning in AplysiaAnnals of the New York Academy of Sciences 627:124–149.https://doi.org/10.1111/j.1749-6632.1991.tb25918.x
-
System-like consolidation of olfactory memories in DrosophilaJournal of Neuroscience 33:9846–9854.https://doi.org/10.1523/JNEUROSCI.0451-13.2013
-
A new brain dopamine-deficient Drosophila and its pharmacological and genetic rescueGenes, Brain and Behavior 16:394–403.https://doi.org/10.1111/gbb.12353
-
Drosophila learn opposing components of a compound food stimulusCurrent Biology 24:1723–1730.https://doi.org/10.1016/j.cub.2014.05.078
-
The cyclic AMP system and Drosophila learningMolecular and Cellular Biochemistry 150:271–278.https://doi.org/10.1007/978-1-4615-2015-3_31
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Localization of short-term memory in the brain of the bee, Apis melliferaPhysiological Entomology 5:343–358.https://doi.org/10.1111/j.1365-3032.1980.tb00244.x
-
Are mushroom bodies cerebellum-like structures?Arthropod Structure & Development 40:368–379.https://doi.org/10.1016/j.asd.2011.02.004
-
Mushroom body memoir: from maps to modelsNature Reviews Neuroscience 4:266–275.https://doi.org/10.1038/nrn1074
-
Encoding of fear learning and memory in distributed neuronal circuitsNature Neuroscience 17:1644–1654.https://doi.org/10.1038/nn.3869
-
Drosophila olfactory memory: single genes to complex neural circuitsNature Reviews Neuroscience 8:341–354.https://doi.org/10.1038/nrn2098
-
Gustatory learning and processing in the Drosophila mushroom bodiesThe Journal of Neuroscience 35:5950–5958.https://doi.org/10.1523/JNEUROSCI.3930-14.2015
-
Bruchpilot, a synaptic active zone protein for anesthesia-resistant memoryJournal of Neuroscience 31:3453–3458.https://doi.org/10.1523/JNEUROSCI.2585-10.2011
-
ConferenceA memory frontier for complex synapsesAdvances in Neural Information Processing Systems 26 (NIPS 2013). pp. 1034–1042.
-
Neural correlates of water reward in thirsty DrosophilaNature Neuroscience 17:1536–1542.https://doi.org/10.1038/nn.3827
-
Quantitative mRNA imaging throughout the entire Drosophila brainNature Methods 14:703–706.https://doi.org/10.1038/nmeth.4309
-
Co-transmission of dopamine and GABA in periglomerular cellsJournal of Neurophysiology 99:1559–1564.https://doi.org/10.1152/jn.00636.2007
-
A theory of cerebellar cortexThe Journal of Physiology 202:437–470.https://doi.org/10.1113/jphysiol.1969.sp008820
-
Thirty years of olfactory learning and memory research in Drosophila MelanogasterProgress in Neurobiology 76:328–347.https://doi.org/10.1016/j.pneurobio.2005.09.003
-
Parallels between cerebellum- and amygdala-dependent conditioningNature Reviews Neuroscience 3:122–131.https://doi.org/10.1038/nrn728
-
Dopamine modulates serotonin innervation in the Drosophila brainFrontiers in Systems Neuroscience 11:76.https://doi.org/10.3389/fnsys.2017.00076
-
Synapsin determines memory strength after punishment- and relief-learningJournal of Neuroscience 35:7487–7502.https://doi.org/10.1523/JNEUROSCI.4454-14.2015
-
Functional consequences of neuropeptide and small-molecule co-transmissionNature Reviews Neuroscience 18:389–403.https://doi.org/10.1038/nrn.2017.56
-
Learning-related plasticity in PE1 and other mushroom body-extrinsic neurons in the honeybee brainJournal of Neuroscience 27:11736–11747.https://doi.org/10.1523/JNEUROSCI.2216-07.2007
-
An aphid sex attractantInsect Systematics & Evolution 1:63–73.https://doi.org/10.1163/187631270X00357
-
Punishment prediction by dopaminergic neurons in DrosophilaCurrent Biology 15:1953–1960.https://doi.org/10.1016/j.cub.2005.09.042
-
Efficient partitioning of memory systems and its importance for memory consolidationPLOS Computational Biology 9:e1003146.https://doi.org/10.1371/journal.pcbi.1003146
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Dopamine and octopamine differentiate between aversive and appetitive olfactory memories in DrosophilaThe Journal of Neuroscience 23:10495–10502.https://doi.org/10.1523/JNEUROSCI.23-33-10495.2003
-
Two Drosophila genes that encode the alph and beta subunits of the brain soluble guanylyl cyclaseThe Journal of Biological Chemistry 270:15368–15376.https://doi.org/10.1074/jbc.270.25.15368
-
The Drosophila Nitric-oxide Synthase Gene ( dNOS ) Encodes a Family of Proteins That Can Modulate NOS Activity by Acting as Dominant Negative RegulatorsJournal of Biological Chemistry 276:42241–42251.https://doi.org/10.1074/jbc.M105066200
-
Dopaminergic terminals in the nucleus accumbens but not the dorsal striatum corelease glutamateJournal of Neuroscience 30:8229–8233.https://doi.org/10.1523/JNEUROSCI.1754-10.2010
-
Dopamine neurons make glutamatergic synapses in vitroThe Journal of Neuroscience 18:4588–4602.https://doi.org/10.1523/JNEUROSCI.18-12-04588.1998
-
Neuronal assemblies of the Drosophila mushroom bodyThe Journal of Comparative Neurology 508:711–755.https://doi.org/10.1002/cne.21692
-
Glutamatergic signaling by mesolimbic dopamine neurons in the nucleus accumbensJournal of Neuroscience 30:7105–7110.https://doi.org/10.1523/JNEUROSCI.0265-10.2010
-
FlyBase 2.0: the next generationNucleic Acids Research 47:D759–D765.https://doi.org/10.1093/nar/gky1003
-
Dopaminergic neurons encode a distributed, asymmetric representation of temperature in DrosophilaJournal of Neuroscience 33:2166–2176.https://doi.org/10.1523/JNEUROSCI.3933-12.2013
-
Classical conditioning and retention in normal and mutant Drosophila MelanogasterJournal of Comparative Physiology A 157:263–277.https://doi.org/10.1007/BF01350033
-
FluoRender: an application of 2D image space methods for 3D and 4D confocal microscopy data visualization in neurobiology researchIEEE Pacific Visualization Symposium : [Proceedings]. IEEE Pacific Visualisation Symposium pp. 201–208.https://doi.org/10.1109/pacificvis.2012.6183592
-
Multiple parallel memory systems in the brain of the ratNeurobiology of Learning and Memory 77:125–184.https://doi.org/10.1006/nlme.2001.4008
-
Nitric oxide synthase is not essential for Drosophila developmentCurrent Biology 20:R141–R142.https://doi.org/10.1016/j.cub.2009.12.011
Article and author information
Author details
Funding
Howard Hughes Medical Institute
- Yoshinori Aso
- Robert P Ray
- Xi Long
- Daniel Bushey
- Teri-TB Ngo
- Brandi Sharp
- Christina Christoforou
- Amy Hu
- Andrew L Lemire
- Paul Tillberg
- Gerald M Rubin
National Institutes of Health (R01 GM84128)
- Karol Cichewicz
- Jay Hirsh
Simons Foundation (Global Brain)
- Yoshinori Aso
- Ashok Litwin-Kumar
- Gerald M Rubin
Burroughs Wellcome Fund
- Ashok Litwin-Kumar
Gatsby Charitable Foundation
- Ashok Litwin-Kumar
National Science Foundation (DBI-1707398)
- Ashok Litwin-Kumar
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Jui-Chun Kao, and the Janelia Fly Facility for help with fly husbandry. We thank Eizaburo Doi, Joshua Dudman, Yichun Shuai, Mehrab Modi, Karen Hibbard Gowan Tervo, Adam Hantman, Daisuke Hattori, Larry Abbott, Vanessa Ruta, Stott Waddell, Hiromu Tanimoto, Toshihide Hige, Brett Mensh and Glenn Turner for stimulating discussions and for comments on earlier drafts of the manuscript. Andrew Lin for fly stocks. Nikita Yakubovich and Patrick O'Farrell shared anti-NOS 16 serum pET28a-dNOS exon 16 plasmid. Phil Borden, Gaby Paez, Silvia Sanchez Martinez, Ariana Tkachuk, Tim Brown, Yan Zhang, Lauren Porter, Jonathan Marvin, and Loren Looger provided guidance on affinity purification of the anti-NOS antibody. Burkhard Poeck and Roland Strauß shared confocal images of fly brains stained with anti-NOS. Daniel E Milkie, Jennifer Colonell, Ruixuan Gao, Srigokul Upadhyayula, Eric Betzig, Michael DeSantis, Damien Alcor and Janelia Advanced Imaging Center for supporting the lattice light sheet imaging. Igor Pisarev, Stephan Saalfeld, Takashi Kawase and Hideo Otsuna provided help with image analysis. Rebecca Vorimo, Allison Sowell, Kari Close, Project Technical Resources, and the FlyLight Project Team at Janelia Research Campus provided brain dissection and histological preparation. The laboratory of JH was supported by NIH R01 GM GM84128.
Copyright
© 2019, Aso 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
-
- 7,783
- views
-
- 1,073
- downloads
-
- 108
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Sour taste, which is elicited by low pH, may serve to help animals distinguish appetitive from potentially harmful food sources. In all species studied to date, the attractiveness of oral acids is contingent on concentration. Many carboxylic acids are attractive at ecologically relevant concentrations but become aversive beyond some maximal concentration. Recent work found that Drosophila ionotropic receptors IR25a and IR76b expressed by sweet-responsive gustatory receptor neurons (GRNs) in the labellum, a peripheral gustatory organ, mediate appetitive feeding behaviors toward dilute carboxylic acids. Here, we disclose the existence of pharyngeal sensors in Drosophila melanogaster that detect ingested carboxylic acids and are also involved in the appetitive responses to carboxylic acids. These pharyngeal sensors rely on IR51b, IR94a, and IR94h, together with IR25a and IR76b, to drive responses to carboxylic acids. We then demonstrate that optogenetic activation of either Ir94a+ or Ir94h+ GRNs promotes an appetitive feeding response, confirming their contributions to appetitive feeding behavior. Our discovery of internal pharyngeal sour taste receptors opens up new avenues for investigating the internal sensation of tastants in insects.