No retraction, correction or concern found when last checked on Oct 21, 2026.
A Bidirectional Circuit Switch Reroutes Pheromone Signals in Male and Female Brains
Johannes Kohl, Aaron D. Ostrovsky, Shahar Frechter and 1 more
Cell · 2013 · CC BY
doi:10.1016/j.cell.2013.11.025
Abstract
The Drosophila sex pheromone cVA elicits different behaviors in males and females. First- and second-order olfactory neurons show identical pheromone responses, suggesting that sex genes differentially wire circuits deeper in the brain. Using in vivo whole-cell electrophysiology, we now show that two clusters of third-order olfactory neurons have dimorphic pheromone responses. One cluster responds in females; the other responds in males. These clusters are present in both sexes and share a common input pathway, but sex-specific wiring reroutes pheromone information. Regulating dendritic position, the fruitless transcription factor both connects the male-responsive cluster and disconnects the female-responsive cluster from pheromone input. Selective masculinization of third-order neurons transforms their morphology and pheromone responses, demonstrating that circuits can be functionally rewired by the cell-autonomous action of a switch gene. This bidirectional switch, analogous to an electrical changeover switch, provides a simple circuit logic to activate different behaviors in males and females.
Full text
From PubMed Central, U.S. National Library of Medicine, as current on Oct 21, 2026. Paragraph numbers are ours, so an answer can point at one.
Introduction
Many species exhibit sexually dimorphic behaviors, typically as part of their reproductive repertoire. These behaviors, which often have a substantial unlearned component, provide highly tractable paradigms to explore the genetic and neural circuit basis of behavior. As potent releasers of specific dimorphic behavior, sex pheromones are particularly experimentally advantageous. Nevertheless, even here the neural mechanisms underlying differential processing within the brain remain largely unknown.
Several models have been proposed for how pheromones can elicit different behavior in males and females. One model is exemplified by classic work on the attraction of male silkmoths to bombykol. Here, one sex expresses a pheromone receptor, while the other is blind to this cue. However, this peripheral change cannot account for situations in which a common stimulus produces behavior in both sexes. These are likely due to circuit differences within the brain. For instance, in mice, only males show courtship behavior toward females, but after ablation of the vomeronasal organ females show female-directed courtship. This leads to a second model in which both sexes can express male behaviors, but these are normally repressed in females by sex-specific circuits downstream of pheromone detection. However, these circuit differences remain unknown, because the relevant receptors and downstream pathways have yet to be identified.
A simpler paradigm is offered by analogous results in flies and mice, in which a monomolecular pheromone can activate identified sensory neurons in both sexes. In the mouse, the male pheromone ESP1 activates V2Rp5 sensory neurons in both sexes but produces distinct patterns of immediate-early gene expression in deeper brain nuclei. ESP1 triggers lordosis in females, but no effect on male behavior has been reported.
In Drosophila, the male pheromone 11-cis-vaccenyl acetate (cVA) stimulates courtship in females but decreases courtship and increases aggression in males. Because both first- and second-order olfactory neurons show similar cVA responses in males and females, it is likely that some circuit difference deeper in the brain results in sex-specific behavioral output. Two further studies have characterized downstream elements of this pathway. used an elegant tracing approach based on sequential photoactivation of green fluorescent protein to identify candidate third- and fourth-order neurons, some of which were shown to be cVA responsive in males. However, they were unable to characterize these neurons anatomically or functionally in females, so the presence or nature of any circuit dimorphism remained unclear. In a parallel study, used a genetic mosaic technique to carry out an exhaustive analysis of sexually dimorphic neurons in male and female brains. In the olfactory system they found two groups of third-order neurons, present in both sexes, that appeared to be differentially connected, suggesting a precise circuit hypothesis for differential pheromone processing in male and female brains (Figure 1A).
We now combine targeted in vivo whole-cell electrophysiology, high-resolution neuroanatomy, and genetic analysis to analyze cVA processing in male, female, and sex mosaic flies. We first establish a simple but efficient circuit motif: a bidirectional (or changeover) switch in which a common input is routed to different active outputs in each sex. We then demonstrate that the fruitless gene sets the state of this switch, specifying both the dendritic placement and pheromone responses of third-order olfactory neurons in a cell-autonomous manner.
Sex-Specific Pheromone Responses in fru+ LHNs
cVA processing in the first three layers of the fly olfactory system provides an ideal model to investigate the logic of neural circuit switches. Or67d olfactory receptor neurons (ORNs) are narrowly tuned to cVA and send axons to the DA1 glomerulus in the brain, where they synapse with DA1 projection neurons (PNs). First-order ORNs and second-order PNs both express the terminal sex determination gene fruitless (henceforth fru+ neurons) but appear functionally isomorphic. Recent anatomical work on fru+ neurons has identified five clusters of candidate third-order neurons of the lateral horn that may receive cVA pheromone information (summarized in Table 1). Each cluster descends from a different neuroblast (neural stem cell). highlighted two neuronal clusters that were present in both sexes but had dendrites in sex-specific locations: aSP-f neurons had dendritic overlap with DA1 PN axon terminals in males, but not females, whereas aSP-g dendrites overlapped in females, but not in males (Figure 1B). Although suggestive, these purely anatomical results provided no functional evidence for a wiring difference that altered pheromone processing.
also characterized aSP-f (DC1) neurons in males. However, negative anatomical observations led to a conclusion that these neurons were absent in females. Critically, then demonstrated that male aSP-f/DC1 neurons receive input from the DA1 glomerulus and respond to cVA in males. However, in the absence of positive anatomical data or physiological recordings in females, it remained unclear whether these responses were sex-specific. Furthermore, this study did not identify neurons that might selectively receive pheromone information in females.
Our anatomical data prompted us to make in vivo recordings from fru+ lateral horn neurons (LHNs) in males and females. We obtained stable whole-cell patch clamp recordings (most >1 hr), giving access to subthreshold responses and morphology of every recorded neuron. Cells were filled, classified, reconstructed (see Figure 1C), and registered to a template brain, allowing us to compare the overlap of LHN dendrites with incoming PN axons. One important technical point quickly became clear: fruGal4 is too weakly expressed in females to target some cells for recording (e.g., aSP-f neurons). However, in the course of a large enhancer trap screen (S.F., J.K., and G.S.X.E.J., unpublished data; Experimental Procedures), we obtained two new driver lines, JK1029 and JK56, that label subsets of fru+ neurons, including the aSP-f cluster, in both sexes (Figure S1 available online).
Because aSP-f dendrites only overlap with DA1 PNs in males (Figure 1D), we expected them to show male-specific responses. Indeed, about half of the male—but almost no female—aSP-f neurons showed significant spiking cVA responses (Figure 1M). cVA-responsive aSP-f neurons in males were narrowly tuned to cVA (Figures 1J and 2G), matching the narrow tuning of Or67d ORNs and DA1 PNs.
Because aSP-g dendrites overlap with DA1 PNs only in females (Figure 1E), we expected them to show female-specific responses. Indeed, the majority of female—but almost no male—aSP-g neurons responded to cVA (Figure 1N). Female aSP-g neurons showed weaker cVA responses and broader odor tuning than did male aSP-f neurons (Table S1C; Figures 1K and 2G); this is likely due to the partial overlap of aSP-g dendrites with incoming DA1 PN axons and suggests that their dendrites also receive information from other PN classes.
We also recorded from a third cluster of fru+ LHNs. aSP-h neurons were examined anatomically in males by (who referred to them as DC2 neurons and proposed that they were absent from females), whereas examined both sexes and observed a difference in the density of dendritic arbors in the ventral lateral horn (see Table 1). Examining single aSP-h neurons filled during recording, we found more dendritic overlap with DA1 PN axons in males than in females (Figures 1F and 1I). Functionally, cVA spiking responses in these broadly tuned neurons were stronger and more frequent in males but were still occasionally present in females (Figures 1L, 1O, and 2C). Given the quantitative nature of this difference and the broader tuning of these neurons, our subsequent analysis focused on aSP-f and aSP-g neurons.
Morphological and Functional Correlations
The difference in dendritic location for male and female aSP-f neurons (Figure 1D) provides a simple circuit hypothesis for the origin of functional differences between these neurons. We confirmed this relationship by examining the three-dimensional (3D) morphology of 37 male aSP-f neurons and 36 female aSP-f neurons that were filled during recording. Almost all male neurons had dendrites in the ventral lateral horn, whereas female neurons never did. In addition to clear differences in dendritic arborization, there was a consistent difference in axonal morphology. Male aSP-f axons terminate in the male-enlarged arch and lateral junction neuropil regions, whereas female aSP-f axons project to the arch and the superior protocerebrum (Figure 2D). aSP-g and aSP-h neurons show a similar axonal dimorphism with axons in females targeting the same superior protocerebral region (Figures 2E and 2F). This region is the female-enlarged neuropil described by that appears to match a focus for female receptivity described by early gynandromorph studies.
The lack of spiking responses in half of the male aSP-f neurons (Figure 1M) was initially surprising, because all but one of these nonresponders had dendrites in the ventral lateral horn (see Figure S2A) with the potential to form synapses with DA1 PNs. However, 9/17 of these neurons showed significant subthreshold cVA responses (see Extended Experimental Procedures), indicating that they do receive input, but that it is unable to drive a spiking response. Morphological analysis of aSP-f neurons revealed two major classes in males, unilateral neurons and bilateral neurons, whose axons project through the arch to the contralateral protocerebrum (Figures 2A and 2D). Intriguingly, cross-comparison of morphology and physiology revealed that all unilateral neurons in our study showed strong spiking responses, whereas responses from bilateral neurons were infrequent and weaker when present (Figures 2A and 2H; see Table S1C). aSP-f neurons therefore have distinct functionally and morphologically related subtypes. Furthermore, these subtypes are genetically heterogeneous, because the JK56 driver line exclusively labels bilateral neurons. This difference may be functionally significant because bilateral male aSP-f neurons have additional dendritic arborizations ventral to the lateral horn (Figure 2H) and may therefore integrate both olfactory and nonolfactory stimuli; coincident inputs would likely result in stronger responses.
Analysis of individual aSP-g (Figures 2B and 2E) and aSP-h (Figures 2C and 2F) neurons clearly revealed the correlated morphological and functional differences between the sexes. However, although aSP-g neurons showed clear morphological subtypes within each sex (Figures 2B and 2E), no strong structure/function correlations were obvious for these subtypes. aSP-h neurons appeared morphologically homogeneous (Figures 2C and 2F).
cVA Responses in fru+ LHNs Depend on a Common Input
Our bidirectional switch model (Figure 1A) predicts that pheromone responses depend on a common sensory pathway. cVA detection has been linked to two classes of sensory neurons that express either olfactory receptor: Or67d or Or65a. However, the available anatomical data suggest that aSP-f neurons in males and aSP-g neurons in females are postsynaptic to DA1 PNs, which receive input from Or67d sensory neurons. We therefore recorded from fru+ LHNs in flies lacking Or67d (Figure 3D). In these Or67d−/− flies, cVA-evoked spiking and subthreshold responses were abolished in both male aSP-f (Figures 3A and 3C) and female aSP-g neurons (Figures 3B and 3C; Table S1C). The absence of even subthreshold cVA responses indicates that Or65a ORNs provide minimal, if any, input to these neurons. Responses of female aSP-g neurons to other odorants were preserved in Or67d−/− flies (Figure 3B). This suggests that female aSP-g neurons integrate cVA information from the Or67d/DA1-labeled line along with general odor information encoded by other ORN/PN classes. In conclusion, the same Or67d sensory pathway is necessary for cVA responses in both aSP-f and aSP-g LHNs, consistent with the bidirectional switch hypothesis.
DA1 PNs Form Sex-Specific Connections with fru+ LHNs
The Or67d receptor is necessary for pheromone responses in fru+ LHNs. Is stimulating this pathway also sufficient to excite these neurons? Or67d sensory neurons project to the DA1 glomerulus, synapsing with DA1 PN dendrites. We used local acetylcholine iontophoresis to stimulate the dendrites of DA1 PNs, while simultaneously recording intracellularly from fru+ LHNs (Figure 3E). DA1 stimulation produced both spiking responses and large depolarizations in all male aSP-f and almost all female aSP-g neurons (Figures 3F–3H; Figure S3A); male aSP-g neurons were unresponsive. Control stimulation in neighboring glomeruli produced minimal responses (Figures 3F and 3H), confirming the specificity of stimulation; this also suggests that glomeruli in more distant parts of the antennal lobe are the origin of non-cVA responses in female aSP-g neurons. We previously noted that only half of aSP-f neurons showed cVA spiking responses, whereas all male aSP-f neurons responded to glomerular stimulation. This suggests that all aSP-f neurons receive input from DA1 PNs, but the strength of this input varies across different morphological classes. It also appears that stimulation can reveal functional connections that are too weak to generate spiking responses to odor stimuli.
The Or67d/DA1 pathway is sufficient for sex-specific excitation of aSP-f and aSP-g LHNs, but do these LHNs receive direct input from DA1 PNs? We measured the latency between presynaptic stimulation and postsynaptic response. Latencies to the first spike were variable though sometimes as low as 4.5 ms. However, whole-cell recordings allowed us to measure the latency to the start of the evoked postsynaptic response (see Figure S3B; Experimental Procedures) in male aSP-f and female aSP-g neurons. We found values of 1.8 ± 0.4 ms and 1.8 ± 0.3 ms (n = 7 each), respectively. This is less than half the reported latency for ORN to PN connections but is consistent with measurements from paired recordings of connected nonpheromonal PNs and LHNs (1.5 ms; ).
Our recordings therefore provide conclusive evidence for sex-specific input from DA1 PNs to aSP-f and aSP-g LHNs, exactly as predicted from their anatomy and odor responses. Furthermore, the very short latency to subthreshold response is strong evidence for a monosynaptic connection.
FruM Is Necessary for the Male Form of the Switch
Our results so far identify a bidirectional switch in pheromone processing (Figure 1A), where a common sensory pathway is wired to different target neurons in male and female animals. What is the genetic basis of this circuit switch? In Drosophila, an alternative splicing cascade converts sex chromosome status into sex-specific action of two terminal transcription factors, fruitless and doublesex. The action of fruitless is confined to the male nervous system, where the protein products of male-specific fruitless transcripts (collectively termed FruM) present in about 2,000 neurons are critical for male behavior. Numerous studies have shown that fruitless loss-of-function mutations can change the morphology of both central and sensory neurons and the survival of central neurons. However, in only one case has fruitless been shown to be necessary for a sexually dimorphic neuronal connection: fruitless is required for survival of the Mind motor neuron in males, which in turn induces the formation of its target, the muscle of Lawrence. We now show a direct effect of fruitless on brain wiring: a functionally validated change in connectivity between identified neurons.
First-, second-, and third-order olfactory neurons associated with pheromone signaling all express FruM protein in males but do not express doublesex. Does FruM therefore specify the male form of the bidirectional circuit switch? We used a heteroallelic loss-of-function combination fruF/ fru4-40 (henceforth fru−/−; Experimental Procedures) to remove FruM from all fruitless-expressing neurons (Figure S4A) and examined the morphology of the 6–7 aSP-f neurons labeled by the sparse JK56 driver line (Figure S1A) in wild-type males, females, and fru−/− males.
In fru−/− males the ventral lateral horn lacked male-specific aSP-f dendrites (Figure S4B), closely resembling the female pattern. We performed whole-cell recordings to examine whether single aSP-f neurons were morphologically and functionally feminized in these fru−/− males. Wild-type female JK56 aSP-f neurons showed no dendritic overlap with DA1 PNs and had unilateral axonal projections, whereas wild-type male aSP-f neurons contacted DA1 PNs and had bilateral axonal projections (Figure 4A). fru−/− neurons had minimal overlap with DA1 PNs and no contralateral projections (Figure 4A) and therefore resembled wild-type female neurons. Morphological clustering confirmed this impression (Figure S4E). Most (10/14) reconstructed fru−/− neurons coclustered with female neurons, with the remaining neurons displaying unusual projections in the dorsal lateral horn, which is never innervated by wild-type neurons of either sex. cVA elicited very weak spiking responses (range 2–6 Hz) in 3/14 fru−/− aSP-f neurons, two of which had aberrant morphology. In contrast 8/20 neurons in control males showed responses (range 4–36 Hz) (Figure 4C). This difference in response magnitude was statistically significant (Table S1C). The few odor-responsive fru−/− aSP-f neurons were also more broadly tuned than were their wild-type male counterparts (Figure 4D; Table S1C). Thus, although this fru−/− heteroallelic combination does not result in feminization of all neurons, FruM is clearly necessary to establish the male form of the circuit switch (Figure 4E). This parallels the observation that the same fruitless loss-of-function combination disrupts normal male courtship but does not lead to full behavioral feminization.
FruM Specifies the Male Form of the Circuit Switch
Many studies have demonstrated that fruitless mutations selectively disrupt male behavior, leading to the influential hypothesis that FruM builds the potential for male sexual behavior into the fly nervous system. This hypothesis was dramatically validated by studies in which misexpression of FruM in females was sufficient to recapitulate many steps of male courtship behavior or to produce masculinized courtship song. If FruM can specify male behavior, can it also specify the male form of the circuit switch in higher olfactory neurons?
We used mosaic analysis with a repressible cell marker (MARCM) to label fru+ neuronal clusters in females expressing FruM in all fru+ neurons (fruM females). Examining MARCM clones labeling 16 sexually dimorphic clusters, we found that nine clusters were indistinguishable from wild-type males, both in morphology and cell number; four clusters were not transformed. Completely masculinized clones included aSP-f (Figure 5A) and aSP-g (Figure 5B), components of the circuit switch. We recorded and filled LHNs in fruM females and found that single reconstructed neurons were morphologically indistinguishable from male neurons (Figures 5C, 5D, and 5G; Figure S5A). cVA elicited spiking responses in about half of fruM female aSP-f neurons (Figures 5E, 5G, and 5H), whereas no fruM aSP-g neurons showed cVA spiking responses (Figures 5F–5H; Table S1). FruM therefore specifies the male form of the circuit switch, coordinating both the connection of aSP-f neurons to, and the disconnection of aSP-g neurons from, pheromone input (Figure 5I).
Selectively Masculinizing fru+ LHNs Can Flip the Switch
The experiments in Figures 4 and 5 demonstrate that FruM is both necessary and sufficient for the male form of the bidirectional circuit switch. FruM is expressed in <5% of the neurons in the fly brain, but these include pheromone-responsive first-, second-, and third-order olfactory neurons. We therefore asked whether selectively masculinizing fru+ LHNs in an otherwise female brain is sufficient to transform these neurons. We used null mutants in the transformer (tra) gene: tra1 mutant females are morphologically and behaviorally completely masculinized and loss of tra can masculinize individual somatic cells in a cell-autonomous manner. Because fru+ LHNs do not express doublesex, any transformation should depend on FruM.
We generated MARCM clones homozygous mutant for tra1, masculinizing these neurons in female brains. aSP-f and aSP-g (Figures 6A and 6B) tra1 clones were morphologically indistinguishable from wild-type male clones, even when aSP-f or aSP-g were the only tra-deficient fru+ clones in the brain. This indicates a cell-autonomous effect. tra1 aSP-h neurons were also masculinized (Figures S5D and S5E).
In order to determine whether this morphological transformation was reflected functionally, we performed whole-cell recordings in females containing labeled (and therefore masculinized) LHN clusters (Figure S6; Experimental Procedures). Because of the stochastic nature of MARCM, each animal was examined on the electrophysiology rig to determine if GFP-labeled mutant clones were present in the lateral horn (n = 297 flies). In total we observed 15 aSP-f, 9 aSP-g, and 17 aSP-h clones, from which we recorded 14, 9, and 8 single neurons, respectively. Morphologically, individual tra1 aSP-f and aSP-g (Figures 6C and 6D) neurons were completely masculinized, coclustering with their wild-type male counterparts (Figure 6G).
Transformed aSP-g neurons had dendrites outside the ventral lateral horn, so we strongly predicted that they would lose their cVA responses. Indeed, no tra1 aSP-g neurons responded to cVA (Figures 6F–6H). Conversely, transformed aSP-f neurons have dendrites close to the axon terminals of female DA1 PNs. Is this sufficient for them to form functional connections? Strikingly, we observed cVA spiking responses in 5/14 tra1 mutant aSP-f neurons (Figures 6E, 6G, and 6H). tra1 mutant aSP-h neurons also gained male-type responses (Figures S5F and S5G; see Table S1C for full statistical tests). We conclude that the cell-autonomous transformation of fru+ LHNs is sufficient to recapitulate male form and function (Figure 6I).
Discussion
Our study reveals principles of neural circuit organization and development that are of general significance. First, we show that two populations of neurons, present in both sexes, show reciprocal, sex-specific responses to the same stimulus. Second, we demonstrate that these responses result from differential wiring of a common input to different outputs. Together, these results define an elegant principle of neural circuit organization: a developmental circuit switch directly analogous to an electrical changeover (or single pole, double throw, SPDT) switch that efficiently reroutes a common input signal to one of two possible outputs. This model appears directly applicable to sex-specific processing of mouse pheromones, including ESP1 and Darcin, but not to Caenorhabditis elegans ascarosides, where recent data suggest wiring differences may not be required. The electrical changeover switch is the prototype for a wide-range of electrical switches in which concerted changes involving three or more contacts reroute signals; it is very likely that neural circuits, including those involved in pheromone processing, contain more complex switches or assemblies of multiple switches that elaborate on the basic mechanism that we have described here. Indeed, we previously identified over 700 sites of dimorphic neuronal overlap that may form such switches in other sensory pathways, multimodal interneurons, or motor circuits across the fly brain.
Third, we identify sex-specific placement of target neuron dendrites as the primary cellular basis of the switch that we have described. This contrasts with earlier studies of this circuit that proposed that axonal dimorphism or neurons present only in one sex were the key dimorphic element. Regarding axonal dimorphism, hypothesized that a male-specific extension of DA1 PN axon terminals is the basis of differential wiring in this system, and subsequently proposed that this extension synapses with the dendrites of aSP-f LHNs in males. The large shifts in dendritic position that we observe in aSP-f and aSP-g neurons mean the male-specific extension of DA1 PNs cannot be sufficient for rewiring. Is it necessary? In our mosaic masculinization experiments, aSP-f and aSP-h neurons adopt male morphology and pheromone responses in a brain in which other neurons (including DA1 PNs) are female. Therefore, the male-specific ventral extension is either not necessary for differential wiring or is a secondary consequence of changes in the dendrites of post-synaptic LHNs. Of course, this extension may increase contact between DA1 PNs and aSP-f and aSP-h LHNs, strengthening responses of those LHNs in males. All three mechanisms (dendritic and axonal dimorphisms, dimorphic cell numbers) are likely relevant to different degrees in different circuits.
Fourth, having defined this bidirectional switch, we demonstrate that its male form is specified by the fruitless gene. We show that this transcription factor has a dual function, coordinating the disconnection of one group of target neurons and the connection of the other. Fifth, we show that masculinization of third-order neurons alone is sufficient for functional rewiring. Although previous studies have demonstrated a cell-autonomous effect of fruitless on neuronal morphology, we now demonstrate a difference in functional connectivity. This is surprising because many would predict that connectivity changes would depend on coordinate regulation of genes in synaptic partner neurons. Such simplicity has evolutionary implications: it may allow variation in circuit structure and ultimately in behavior, through evolution of cis-regulatory elements, as previously shown for somatic characters, such as wing spots.
Sixth, studies of pheromone processing in general and cVA processing in particular have emphasized a labeled line processing model. However, our data indicate that both narrowly (aSP-f) and broadly tuned (aSP-h) cVA-responsive neurons coexist in males. Likewise in females, aSP-g neurons respond to cVA and general odors, such as vinegar, but only cVA responses depend on the Or67d receptor. It will be very interesting to determine the circuit origin and behavioral significance of this integration of odor channels. For example, it seems reasonable to speculate that coincidence of cVA and food odors could interact in a supralinear way to promote female courtship or egg laying. This parallels the convergence in the lateral horn of a labeled line responsive to non-cVA fly odors (Or47b/VA1lm neurons) and one responsive to a specific food odorant, phenylacetic acid, that acts as a male aphrodisiac.
Our study naturally raises additional questions. The action of fruitless within fewer than 5% of the neurons in the fly brain can specify behavior, and we now show that it can reroute pheromone signals within those neurons. But what is the behavioral relevance of this particular bidirectional switch? Testing this will require the development of sensitive behavioral assays of cVA processing and a reliable genetic approach to control this switch without affecting the many other dimorphic elements in sensory and motor circuits. Indeed, it remains to be seen whether flipping a single switch in sensory processing is sufficient to engage motor behavior typical of the opposite sex without masculinizing downstream circuitry. We note that could force the production of courtship song by activating fruitless-positive neurons in headless females but were almost never successful in intact females.
Another open question concerns the functional significance of female aSP-f and male aSP-g neurons, which do not respond to cVA or other tested odors. Do they receive input at all? One possibility, based on our in silico analysis of the brain-wide 3D maps in, is that they receive gustatory input, perhaps from contact pheromones, although further work is necessary to test this hypothesis. Finally, which genes does fruitless regulate in order to differentially wire the switch? Our clonal transformation experiments strongly support our earlier proposal that male and female aSP-f/g/h clusters are generated by neuroblasts common to both sexes but that those neurons develop in a sex-specific manner. Therefore, cell-surface molecules required for dendritic guidance are plausible targets. It will be intriguing to see if the same fru-dependent factor(s) direct(s) male aSP-f and female aSP-g dendrites to the ventral lateral horn and, more generally, whether fruitless acts on conserved downstream targets across all the dimorphic neurons in the fly brain.
Fly Stocks
The fruitlessGal4 (fruGal4), fruF, fruM, tra1, and Or67dGal4 stocks were as described previously (see Extended Experimental Procedures). JK56 and JK1029 are Split Gal4 enhancer trap P-element insertions of the Herpes Simplex VP16 activation domain (VP16-AD), which were identified in a screen of 2,000 new insertions generated by our group (Extended Experimental Procedures). MARCM labeling of tra1 mutant clones used y w hs-FLP UAS-mCD8-GFP / +; UAS-mCD8-GFP FRTG13 / +; tra1 FRT2A fruGal4 / tubP-Gal80 FRT2A flies. For anatomical experiments, MARCM clones were generated by heat shock of first-instar larvae for 17 min (males) or 23 min (females) at 37°C 0 hr–3 hr after larval hatching. For whole-cell recordings of tra1 mutant clones, heatshock time was extended to 1.5 hr (males and females), increasing clone frequency.
Immunochemistry
Immunochemistry was as described previously, except that blocking was overnight at 4°C. For FruM staining, fixation was in 2% PFA for 30 min on ice. Primary antibodies included mouse anti-nc82 (DSHB, University of Iowa) 1:20–1:40, chicken anti-GFP (Abcam, ab13970), and rabbit anti-FruM (rabbit polyclonal against male-specific 101 amino acids of FruM, gift of S. Goodwin) 1:400. Secondary antibodies (all from Life Technologies) included Alexa-568 anti-mouse (A-11004) 1:1,200, Alexa-633 anti-mouse (A-21052) 1:1,200, Alexa-488 anti-chicken (A-11039) 1:1,200, and Alexa-568 anti-rabbit (A-11011) 1:1,200. Filled neurons were visualized with Streptavidin Alexa-568 (S-11226) 1:1,300.
Image Acquisition and Analysis
Confocal stacks were acquired on a Zeiss 710 with a 40× NA1.3 oil objective, voxel resolution 0.46 × 0.46 × 1 μm. Images were registered to the IS2 template brain with the Computational Morphometry Toolkit (CMTK, http://www.nitrc.org/projects/cmtk). Neuron tracing used the skeletonize module in Amira (VSG). Tracings were transformed to the left brain hemisphere using the AnalysisSuite package (https://github.com/jefferis/AnalysisSuite) written in R (http://www.r-project.org). Amira was used for 3D visualization. See http://jefferislab.org/si/frulhns for details and data download. Morphological analysis of traced neurons in R used an algorithm that scores the similarity of the local geometry of two neurons by calculating the distance between matching points and the dot products of the tangent vectors (see Extended Experimental Procedures for details and links to R code).
Electrophysiology
Recordings were made from 2- to 3-day-old flies essentially as described previously, with the changes indicated in the Extended Experimental Procedures. A different protocol was developed for recording mutant clones in tra1 MARCM females (Extended Experimental Procedures). Single glomerulus stimulation was performed largely as described by, with modifications indicated in the Extended Experimental Procedures. Field recordings were performed to ensure that animals were odor-responsive (Extended Experimental Procedures). Data acquisition and initial analysis were carried out in Igor Pro with the NeuroMatic analysis software package (Jason Rothman, University College London; see http://neuromatic.thinkrandom.com); subsequent analysis was in R (Extended Experimental Procedures).
We quantified odor responses by finding the mean spike number in 500 ms window starting 150 ms after valve opening, subtracting the mean spike number for control stimulus. We assessed significance by an exact one-sided Poisson test of the number of spikes to odor and control stimuli using data from four trials per cell. We adjusted raw p values to control the false discovery rate using R’s p.adjust function; cells were declared significant for FDR adjusted p < 0.01.
Odor Stimulation
Odorant delivery used a custom odor delivery device (ODD; Extended Experimental Procedures and http://jefferislab.org/si/odd). All odorants were of the highest purity available and were prepared 1:100 v/v in mineral oil (Sigma, M8410), except propionic, butyric, and acetic acid, which were dissolved 1:100 v/v in water, and phenylacetic acid, which was diluted 1:200 w/v in water. cVA was undiluted.
Odorant abbreviations include ctr, mineral oil control; cVA, 11-cis-vaccenyl acetate; 4ol, butanol; PAA, phenylacetic aldehyde; IAA, isoamyl acetate; pro, propionic acid; far, farnesol; vin, apple cider vinegar; pac, phenylacetic acid; aac, acetic acid; ger, geranyl acetate; lin, linalool; bty, butyric acid; hxe, E2-hexenal; ben, benzaldehyde; met, methyl salicylate; pra, propyl acetate; hxa, 1-hexanol; ehb, ethyl 3-hydroxybutyrate; eta, ethyl acetate; cit, b-citronellol. (See http://jefferislab.org/si/frulhns for detailed odorant descriptions.)
Additional data and computer code are available at http://jefferislab.org/si/frulhns. All 3D image data are available on a hard drive on request to GSXEJ.The fruitlessGal4 (fruGal4) allele is a targeted insertion of the yeast transcription factor Gal4 into the P1 promoter of the fruitless gene. In fruM mutants, the entire 1,601 bp female-specific part of the S exon is deleted, enforcing male-specific splicing and thus expression of male-specific FruM isoforms. In fruF mutants, point mutations introduced at the male splice donor site of the S exon of fruitless abolish splicing at this site without altering the coding potential of the unspliced transcripts. fru4-40 is a lethal deletion in the fruitless locus, extending distally from fru4 for > 70 kb. Heteroallelic fruF / fru4-40 flies are referred to as fru−/− in this study. In tra1 loss-of-function mutants, the entire coding region of transformer is deleted. The JK56 insertion was mapped by inverse PCR to cytological location 91B5, genomic coordinate 3R: 14410045, ∼40 kb upstream of the fruitless locus. The JK1029 insertion is located on 3R, close to the fruitless locus, but a genomic location is not yet available. Both JK56 and JK1029 were recombined with Cha-Gal4-DBD, an insertion that expresses a Gal4 DNA binding domain-leucine zipper fusion protein under control of the Cha promoter. Or67dGal4 is a mutant knock-in allele in which the open reading frame of Or67d is replaced with that of the yeast transcriptional activator Gal4. Flies homozygous for Or67dGal4 are referred to as Or67d−/− in this study.Wild-type neuron clusters were targeted for whole-cell recordings in flies of the following genotypes:UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / FRTG13 UAS-mCD8-GFP; fruGal4 / +UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK56-VP16AD Cha-Gal4-DBD / +UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD / +Single glomerulus stimulation experiments were performed in flies of the following genotypes:UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / Mz19-Gal4 UAS-mCD8-GFP; JK1029-VP16AD Cha-Gal4-DBD / +UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD Or67dGal4 / +Or67d mutant recordings were performed in flies of the genotype:UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD Or67dGal4 / Or67dGal4fruM mutant clusters were targeted in female flies of the genotypes:UAS-mCD8-GFP / +; UAS-mCD8-GFP / +; fruGal4 / fruMUAS-mCD8-GFP / +; UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD / fruMfruF mutant clusters were targeted in male flies (fru−/− males) of the genotype:UAS-mCD8-GFP / Y; UAS-mCD8-GFP / +; fru4-40 Cha-Gal4-DBD / fruF JK56-VP16ADtra1 mutant clones were targeted in female flies of the genotype:y w hs-FLP UAS-mCD8-GFP / +; UAS-mCD8-GFP FRTG13 / +; tra1 FRT2A fruGal4 / tubP-Gal80 FRT2AMARCM labeling of wild-type LHN clones used flies of the genotype:y w hs-FLP UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP FRTG13 / FRTG13 tubP-GAL80; fruGal4 / +.MARCM labeling of fruM mutant clones used flies of the genotype:y w hs-FLP UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP FRTG13 / FRTG13 tubP-GAL80; fruGal4 / fruMProlonged incubation (2–3 days rotating at 4°C) with primary and secondary antibodies was required for homogeneous staining. Specimens were whole mounted in Vectashield (Vector Laboratories) on charged slides to avoid movement.Brains were imaged at 768 × 768 pixel resolution and 0.6 zoom factor. Images of dye-filled neurons were acquired with 2x (frame) averaging. Detail images were taken with a Plan-Apochromat 63x/1.4 Oil objective at 2–3x zoom and contained about 30 768 × 768 pixel slices with a voxel size of 0.06 × 0.06 × 0.15 μm. All images were taken using 16bit color depths.For morphological analysis, the fully connected reconstruction was transformed into a “DotProperties” representation that retains only the 3D position of points along with the local heading (tangent vector) of the neurite. Neurons in this representation could then be compared in pairwise fashion. For each point on a neuron designated as the query neuron, the closest point on the target neuron was identified using a nearest neighbor algorithm (R package RANN, http://cran.r-project.org/web/packages/RANN). Each point pair was then given a score which was a function of the distance between them multiplied by the absolute dot product of the two tangent vectors. The similarity score for these two neurons was then simply the average of all the scores for all the individual point pairs. The similarity score S(Q,T) (and the corresponding distance score, D(Q,T)) for a query neuron, Q, and a target neuron, T, are precisely defined as follows:where n is the number of points in the query neuron, dij is the distance between point i in the query neuron and its nearest neighbor (point j) in the target neuron and qi and tj are the tangent vectors at these points. There is one free parameter, σ, which determines how close in space points must be to be considered similar; we set this value to 3 μm based on previous estimates of registration accuracy/biological variability in the fly brain. Similarity scores will range between 1 (identical) and 0 (completely different). This approach is implemented in R in the AnalysisSuite function WeightedNNBasedLinesetMatching (https://github.com/jefferis/AnalysisSuite/blob/master/R/Code/NeuroBlast/NeuriteMatchingFunctions.R#L164). For clustering, pairwise similarity scores were converted to distances by subtraction from 1 and clustering was performed using R function hclust using Ward’s method.Flies were reared on standard fly food at a 12-12 light-dark cycle at 25°C and 60% humidity. Fly preparations and recordings were essentially performed as described, with the following changes. UV-curable glue (Kemxert, KOA-300) was used to fix the fly in custom holders which were photofabricated from a 0.025 mm thick sheet of 302 grade stainless steel (Photofabrication Limited). After the forelegs were removed, the proboscis was fixed to a piece of hair running perpendicular to the longitudinal axis of the body of the fly with wax. To prevent motion in the brain, muscle 16 (paired muscle dorsal to the esophagus) was removed with fine forceps. Gentle mechanical desheathing of the exposed brain was performed using fine forceps. Enzymatic desheathing was necessary in most cases to access cell bodies for whole-cell recordings. For this purpose, wide (∼15 μm) desheathing pipettes were pulled from thin-filament (1.5 mm O.D. x 1.17 mm I.D.) borosilicate glass capillaries (Harvard Apparatus) and filled with a 0.22 μm syringe-filtered solution of 0.5 mg ml−1 collagenase type IV (Worthington or Abnova) in saline (see below). Gentle positive (∼30 mm Hg) pressure was applied to rupture the glial sheath and then collagenase was allowed to diffuse into the tissue under reduced pressure (∼10 mm Hg) for 1–2 min. Perfusion rate during experiments was 0.25 ml min-1, controlled by a peristaltic pump. Saline was continuously bubbled with 95% O2/ 5% CO2 and frequently exchanged during the dissection and recording procedure. Only one neuron per brain was recorded and only neurons for which both physiological and morphological data (i.e., biocytin fill) could be obtained were included in this study (see Experimental Procedures; Tables S1A and S1B).Because of the low probability of finding a transformed (i.e., GFP-positive) LHN clone in tra1 MARCM females, typically 10–20 animals had to be screened before the recording procedure could begin. Therefore, as soon as a fly was inserted into the holder, the next vial was put on ice in order to avoid unnecessary delays. Only 2–4 females were kept per vial to avoid repeated anesthesia over short periods of time. The proboscis was immobilized (see above) and tracheae were removed, but muscle 16 was left intact and no mechanical desheathing was initially performed. Flies were then transferred onto the microscope platform and only if labeled LHN clones were present in the right hemisphere the remaining steps were performed and recordings initiated. Recordings were only performed from neuroblast clones since reliable identification of single-cell clones was challenging.Patch-clamp recordings were obtained using an Olympus BX51W1 upright microscope with IR-DIC illumination and a LUMPlanFL/IR 60x water immersion objective. A Grasshopper 143S camera (Point Grey Research) controlled by Micro-Manager software (http://www.micro-manager.org) was used to guide the recording electrode. Patch-clamp electrodes (7–9 MΩ) were pulled from thin-filament (1.5 mm O.D. x 1.17 mm I.D.) borosilicate glass capillaries (Harvard Apparatus) using a Zeitz DMZ Universal Puller (Zeitz - Instruments Vertriebs GmbH) and pressure-polished with a CPM-2 Coating and Polishing Microforge (ALA Scientific Instruments) mounted onto an Olympus CKX41 inverted microscope with a 50x LWD lens. Patch pipettes were filled with a solution of (in mM): Potassium Aspartate 125, CaCl2 0.1, HEPES 10, MgATP 4, NaGTP 0.5, EGTA 1.1, biocytin hydrazide (Life Technologies, B-1603) 13 (pH = 7.2, adjusted to 265 mOsm), that was filtered through a 0.22 μm syringe filter. Saline (external solution) was composed as follows (in mM): NaCl 103, KCl 3, NaH2PO4 1, MgCl2.6H20 4, CaCl2.2H20 1.5, NaHCO3 26, TES 5, glucose 10, trehalose 10 (pH = 7.25, adjusted to 275 mOsm). All chemicals were purchased from Sigma at the highest purity grade available, unless indicated otherwise. Signals were acquired on an Axoclamp 200B amplifier (Molecular Devices) via a CV-7B headstage and low-pass filtered at 10 kHz, digitized at 11.1kHz. Note that although the digitization frequency was below the Nyquist frequency for the amplifier’s 10 kHz filter, we confirmed that any aliasing artifact in the 100-1000Hz range was not detectable (i.e., of much lower amplitude than other sources of noise) and had no impact on spike finding. Recording protocols were delivered via Igor Pro software (Wavemetrics) in conjunction with Nclamp (Jason Rothman, University College London, UK, see http://neuromatic.thinkrandom.com).Field recordings were performed as follows: a patch pipette was filled with 0.22 μm syringe-filtered saline and lowered into the cell body layer dorsal of the mushroom body calyx under positive pressure (∼10 mm Hg). Odor-evoked field responses were then recorded in current clamp mode. Field odor responses were typically visible as 0.1–0.3 mV downward potential deflections (Figure S2C).Initial analysis was carried out in Igor Pro using the NeuroMatic analysis software package (Jason Rothman, University College London, UK, see http://neuromatic.thinkrandom.com). Cells were accepted for analysis if spikes could be elicited by current injection and extracted using the following protocol. Current clamp traces were smoothed with a binomial algorithm at 180–550 Hz and action potentials (‘spikes’) were extracted by applying a negative threshold to the second derivative of the smoothed trace. The resultant spike times were then imported into the R statistical environment for further processing, statistical analysis and production of final graphics (our open source R code is available for download at http://jefferislab.org/si/frulhns). All quantities reported in the paper are mean ± s.d. unless otherwise indicated.The odor tuning of responsive neurons was quantified using a measure of tuning curve sharpness (S, lifetime sparseness, ). Lifetime sparseness ranges from 0 (for a neuron that responds equally to all odorants) to 1 (for a neuron that is excited by one odorant while showing no response to all others). In order to compute peak firing rates for aSP-f and aSP-g neurons, we calculated response rates in five 200 ms windows starting at 150 ms after valve opening, reporting the peak value for each cell.When comparing odor responses between groups of cells we used a permutation approach implemented in the R library coin, specifically the oneway_test function. This is a nonparametric two-sample test analogous to a t test or Wilcoxon text to detect a shift between two populations, with the advantage over a Wilcoxon test that the magnitude (not just the rank) of responses is considered. The permutation approach means that there are no distributional assumptions (including equality of sample variances).In order to assess the significance of subthreshold responses to cVA we first smoothed the membrane potential with a 45ms boxcar filter. We then calculated the peak response in a window from +200 to +800 ms after valve opening, subtracting a baseline value of the mean response between −1000 to 0 ms. We repeated the same procedure for a response window preceding odor presentation, generating a null distribution of pseudo responses from all of the traces available for each cell. The cVA responses were then compared with the pseudo responses using the non parametric oneway_test of the R coin package. Raw p values were adjusted to control the False Discovery Rate using the p.adjust function in R; responses were declared significant for FDR adjusted p < 0.01.Odorants were kept in glass vials with glass wool (Acros Organics) to increase the effective surface area, thus facilitating odor headspace formation in the vial. A constant airstream of 2 l min-1 was directed at the fly throughout the recording. This airstream was composed of a carrier stream (1.75 l min-1) and an odor stream (0.25 l min-1) that was directed from a mineral oil control vial to an odor vial by a trigger-controlled solenoid valve. The rate between carrier and odor stream was set via two mass flow controllers (MFCs). The two streams were joined 70 cm from the end of the delivery tube, which measured 4 mm in diameter and was positioned 8 mm from the fly. All odor concentrations, as they leave the delivery tube, are therefore ∼800x dilutions. Typical odor response latencies with this setup were 150–200 ms, as measured by spiking onset in our recordings and with a photoionization detector (mini-PID, Aurora Scientific) that was placed near the end of the delivery tube. In addition, we sometimes observed a quasi-instantaneous spiking onset in neurons (latency < 5 ms) that we ascribe to a change in pressure.Several precautions were taken to account for false negative results, i.e., recordings from flies that were unable to respond to odors per se: (a) successful odor delivery was monitored with a photoionization detector in series with the animal (see above), (b) great care was taken during the preparation to fully expose the olfactory appendages of the fly without damaging them or the associated olfactory nerves, (c) in order to determine whether neurons were capable of generating action potentials (“spikes”), positive current was injected, and (d) neurons were inspected to establish if they exhibited at least one supra- or subthreshold response to odor. Recordings were only included if criteria (a–c) were fulfilled; if (d) was not fulfilled, a field recording (Figure S2C) was performed to ensure that the animal was odor-responsive (see above).For single glomerulus stimulations, theta glass stimulation electrodes (Warner Instruments) were pulled with a 3–5 μm tip and filled with 2 mM acetylcholine chloride (Sigma) in saline. With the exception of the experiment in Figure S3D, both sides of the theta electrode were filled with identical solution (i.e., effectively a one-barreled electrode), using a separate bath electrode as ground. While there is therefore no liquid junction potential between the two barrels of the theta electrode, there is a liquid junction potential between the bath electrode (in saline) and the theta electrode (saline + acetylcholine). However, since the concentration of acetylcholine used in these experiments was very low (2 mM), this liquid junction potential is less than 0.2 mV.Stimulation electrodes were inserted into the GFP-positive DA1 glomerulus. These recordings were performed in flies with an additional copy of Mz19-Gal4 or Or67dGal4 in order to reliably identify the DA1 glomerulus. In order to reduce spontaneous activity in DA1 PNs, both antennae were removed with forceps. Square voltage pulses (500 ms), ranging from 0.2–10 V were delivered via the DAC output of the ITC-18 data acquisition interface (HEKA Elektronik) using a common bath electrode. This configuration results in both transient and sustained stimulation artifacts; the sustained artifacts (typically 1 mV/V stimulation) were fitted for each cell and subtracted before quantification. Stimulation depolarizes PNs both chemically and electrically (see Figure S3D). For control stimulations the stimulation electrode was moved to neighboring glomeruli in the dorsal antennal lobe. Note that although this controls for stimulation spreading ’across’ glomeruli, it does not rule out input from more ventral glomeruli to (broadly odor-tuned) female aSP-g neurons. After control stimulations of non-DA1 glomeruli, the stimulation electrode was repositioned in DA1 to ensure that there was no drop in responsivity throughout the recording procedure. A small number of experiments investigated the minimum stimulation voltage required to elicit a postsynaptic response, using stimulation intensities as low as 10 mV.To measure stimulation latency we first applied a digital boxcar filter at 2.2 kHz to the signals. We normalized the voltage level by subtracting the mean value between 0 and 0.5 ms after stimulation. Responses to both DA1 and control glomerulus stimulation were plotted (Figure S3B) and binned in 0.5 ms intervals. We then tested for divergence between voltage traces for DA1 and control glomerulus stimulation by carrying out a Wilcoxon rank sum test within each 0.5 ms bin. The latency was then defined as the time after which the DA1 trace remained significantly larger (p < 0.05, one-tailed) than the control trace for the remainder of the stimulation period.
Additional data and computer code are available at http://jefferislab.org/si/frulhns. All 3D image data are available on a hard drive on request to GSXEJ.
Fly Stocks
The fruitlessGal4 (fruGal4) allele is a targeted insertion of the yeast transcription factor Gal4 into the P1 promoter of the fruitless gene. In fruM mutants, the entire 1,601 bp female-specific part of the S exon is deleted, enforcing male-specific splicing and thus expression of male-specific FruM isoforms. In fruF mutants, point mutations introduced at the male splice donor site of the S exon of fruitless abolish splicing at this site without altering the coding potential of the unspliced transcripts. fru4-40 is a lethal deletion in the fruitless locus, extending distally from fru4 for > 70 kb. Heteroallelic fruF / fru4-40 flies are referred to as fru−/− in this study. In tra1 loss-of-function mutants, the entire coding region of transformer is deleted. The JK56 insertion was mapped by inverse PCR to cytological location 91B5, genomic coordinate 3R: 14410045, ∼40 kb upstream of the fruitless locus. The JK1029 insertion is located on 3R, close to the fruitless locus, but a genomic location is not yet available. Both JK56 and JK1029 were recombined with Cha-Gal4-DBD, an insertion that expresses a Gal4 DNA binding domain-leucine zipper fusion protein under control of the Cha promoter. Or67dGal4 is a mutant knock-in allele in which the open reading frame of Or67d is replaced with that of the yeast transcriptional activator Gal4. Flies homozygous for Or67dGal4 are referred to as Or67d−/− in this study.
Wild-type neuron clusters were targeted for whole-cell recordings in flies of the following genotypes:
UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / FRTG13 UAS-mCD8-GFP; fruGal4 / +
UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK56-VP16AD Cha-Gal4-DBD / +
UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD / +
Single glomerulus stimulation experiments were performed in flies of the following genotypes:
UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / Mz19-Gal4 UAS-mCD8-GFP; JK1029-VP16AD Cha-Gal4-DBD / +
UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD Or67dGal4 / +
Or67d mutant recordings were performed in flies of the genotype:
UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD Or67dGal4 / Or67dGal4
fruM mutant clusters were targeted in female flies of the genotypes:
UAS-mCD8-GFP / +; UAS-mCD8-GFP / +; fruGal4 / fruM
UAS-mCD8-GFP / +; UAS-mCD8-GFP / +; JK1029-VP16AD Cha-Gal4-DBD / fruM
fruF mutant clusters were targeted in male flies (fru−/− males) of the genotype:
UAS-mCD8-GFP / Y; UAS-mCD8-GFP / +; fru4-40 Cha-Gal4-DBD / fruF JK56-VP16AD
tra1 mutant clones were targeted in female flies of the genotype:
y w hs-FLP UAS-mCD8-GFP / +; UAS-mCD8-GFP FRTG13 / +; tra1 FRT2A fruGal4 / tubP-Gal80 FRT2A
MARCM
MARCM labeling of wild-type LHN clones used flies of the genotype:
y w hs-FLP UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP FRTG13 / FRTG13 tubP-GAL80; fruGal4 / +.
MARCM labeling of fruM mutant clones used flies of the genotype:
y w hs-FLP UAS-mCD8-GFP / (+ or Y); UAS-mCD8-GFP FRTG13 / FRTG13 tubP-GAL80; fruGal4 / fruM
Immunohistochemistry
Prolonged incubation (2–3 days rotating at 4°C) with primary and secondary antibodies was required for homogeneous staining. Specimens were whole mounted in Vectashield (Vector Laboratories) on charged slides to avoid movement.
Image Acquisition and Analysis
Brains were imaged at 768 × 768 pixel resolution and 0.6 zoom factor. Images of dye-filled neurons were acquired with 2x (frame) averaging. Detail images were taken with a Plan-Apochromat 63x/1.4 Oil objective at 2–3x zoom and contained about 30 768 × 768 pixel slices with a voxel size of 0.06 × 0.06 × 0.15 μm. All images were taken using 16bit color depths.
For morphological analysis, the fully connected reconstruction was transformed into a “DotProperties” representation that retains only the 3D position of points along with the local heading (tangent vector) of the neurite. Neurons in this representation could then be compared in pairwise fashion. For each point on a neuron designated as the query neuron, the closest point on the target neuron was identified using a nearest neighbor algorithm (R package RANN, http://cran.r-project.org/web/packages/RANN). Each point pair was then given a score which was a function of the distance between them multiplied by the absolute dot product of the two tangent vectors. The similarity score for these two neurons was then simply the average of all the scores for all the individual point pairs. The similarity score S(Q,T) (and the corresponding distance score, D(Q,T)) for a query neuron, Q, and a target neuron, T, are precisely defined as follows:where n is the number of points in the query neuron, dij is the distance between point i in the query neuron and its nearest neighbor (point j) in the target neuron and qi and tj are the tangent vectors at these points. There is one free parameter, σ, which determines how close in space points must be to be considered similar; we set this value to 3 μm based on previous estimates of registration accuracy/biological variability in the fly brain. Similarity scores will range between 1 (identical) and 0 (completely different). This approach is implemented in R in the AnalysisSuite function WeightedNNBasedLinesetMatching (https://github.com/jefferis/AnalysisSuite/blob/master/R/Code/NeuroBlast/NeuriteMatchingFunctions.R#L164). For clustering, pairwise similarity scores were converted to distances by subtraction from 1 and clustering was performed using R function hclust using Ward’s method.
Electrophysiology
Flies were reared on standard fly food at a 12-12 light-dark cycle at 25°C and 60% humidity. Fly preparations and recordings were essentially performed as described, with the following changes. UV-curable glue (Kemxert, KOA-300) was used to fix the fly in custom holders which were photofabricated from a 0.025 mm thick sheet of 302 grade stainless steel (Photofabrication Limited). After the forelegs were removed, the proboscis was fixed to a piece of hair running perpendicular to the longitudinal axis of the body of the fly with wax. To prevent motion in the brain, muscle 16 (paired muscle dorsal to the esophagus) was removed with fine forceps. Gentle mechanical desheathing of the exposed brain was performed using fine forceps. Enzymatic desheathing was necessary in most cases to access cell bodies for whole-cell recordings. For this purpose, wide (∼15 μm) desheathing pipettes were pulled from thin-filament (1.5 mm O.D. x 1.17 mm I.D.) borosilicate glass capillaries (Harvard Apparatus) and filled with a 0.22 μm syringe-filtered solution of 0.5 mg ml−1 collagenase type IV (Worthington or Abnova) in saline (see below). Gentle positive (∼30 mm Hg) pressure was applied to rupture the glial sheath and then collagenase was allowed to diffuse into the tissue under reduced pressure (∼10 mm Hg) for 1–2 min. Perfusion rate during experiments was 0.25 ml min-1, controlled by a peristaltic pump. Saline was continuously bubbled with 95% O2/ 5% CO2 and frequently exchanged during the dissection and recording procedure. Only one neuron per brain was recorded and only neurons for which both physiological and morphological data (i.e., biocytin fill) could be obtained were included in this study (see Experimental Procedures; Tables S1A and S1B).
Because of the low probability of finding a transformed (i.e., GFP-positive) LHN clone in tra1 MARCM females, typically 10–20 animals had to be screened before the recording procedure could begin. Therefore, as soon as a fly was inserted into the holder, the next vial was put on ice in order to avoid unnecessary delays. Only 2–4 females were kept per vial to avoid repeated anesthesia over short periods of time. The proboscis was immobilized (see above) and tracheae were removed, but muscle 16 was left intact and no mechanical desheathing was initially performed. Flies were then transferred onto the microscope platform and only if labeled LHN clones were present in the right hemisphere the remaining steps were performed and recordings initiated. Recordings were only performed from neuroblast clones since reliable identification of single-cell clones was challenging.
Patch-clamp recordings were obtained using an Olympus BX51W1 upright microscope with IR-DIC illumination and a LUMPlanFL/IR 60x water immersion objective. A Grasshopper 143S camera (Point Grey Research) controlled by Micro-Manager software (http://www.micro-manager.org) was used to guide the recording electrode. Patch-clamp electrodes (7–9 MΩ) were pulled from thin-filament (1.5 mm O.D. x 1.17 mm I.D.) borosilicate glass capillaries (Harvard Apparatus) using a Zeitz DMZ Universal Puller (Zeitz - Instruments Vertriebs GmbH) and pressure-polished with a CPM-2 Coating and Polishing Microforge (ALA Scientific Instruments) mounted onto an Olympus CKX41 inverted microscope with a 50x LWD lens. Patch pipettes were filled with a solution of (in mM): Potassium Aspartate 125, CaCl2 0.1, HEPES 10, MgATP 4, NaGTP 0.5, EGTA 1.1, biocytin hydrazide (Life Technologies, B-1603) 13 (pH = 7.2, adjusted to 265 mOsm), that was filtered through a 0.22 μm syringe filter. Saline (external solution) was composed as follows (in mM): NaCl 103, KCl 3, NaH2PO4 1, MgCl2.6H20 4, CaCl2.2H20 1.5, NaHCO3 26, TES 5, glucose 10, trehalose 10 (pH = 7.25, adjusted to 275 mOsm). All chemicals were purchased from Sigma at the highest purity grade available, unless indicated otherwise. Signals were acquired on an Axoclamp 200B amplifier (Molecular Devices) via a CV-7B headstage and low-pass filtered at 10 kHz, digitized at 11.1kHz. Note that although the digitization frequency was below the Nyquist frequency for the amplifier’s 10 kHz filter, we confirmed that any aliasing artifact in the 100-1000Hz range was not detectable (i.e., of much lower amplitude than other sources of noise) and had no impact on spike finding. Recording protocols were delivered via Igor Pro software (Wavemetrics) in conjunction with Nclamp (Jason Rothman, University College London, UK, see http://neuromatic.thinkrandom.com).
Field recordings were performed as follows: a patch pipette was filled with 0.22 μm syringe-filtered saline and lowered into the cell body layer dorsal of the mushroom body calyx under positive pressure (∼10 mm Hg). Odor-evoked field responses were then recorded in current clamp mode. Field odor responses were typically visible as 0.1–0.3 mV downward potential deflections (Figure S2C).
Analysis
Initial analysis was carried out in Igor Pro using the NeuroMatic analysis software package (Jason Rothman, University College London, UK, see http://neuromatic.thinkrandom.com). Cells were accepted for analysis if spikes could be elicited by current injection and extracted using the following protocol. Current clamp traces were smoothed with a binomial algorithm at 180–550 Hz and action potentials (‘spikes’) were extracted by applying a negative threshold to the second derivative of the smoothed trace. The resultant spike times were then imported into the R statistical environment for further processing, statistical analysis and production of final graphics (our open source R code is available for download at http://jefferislab.org/si/frulhns). All quantities reported in the paper are mean ± s.d. unless otherwise indicated.
The odor tuning of responsive neurons was quantified using a measure of tuning curve sharpness (S, lifetime sparseness, ). Lifetime sparseness ranges from 0 (for a neuron that responds equally to all odorants) to 1 (for a neuron that is excited by one odorant while showing no response to all others). In order to compute peak firing rates for aSP-f and aSP-g neurons, we calculated response rates in five 200 ms windows starting at 150 ms after valve opening, reporting the peak value for each cell.
When comparing odor responses between groups of cells we used a permutation approach implemented in the R library coin, specifically the oneway_test function. This is a nonparametric two-sample test analogous to a t test or Wilcoxon text to detect a shift between two populations, with the advantage over a Wilcoxon test that the magnitude (not just the rank) of responses is considered. The permutation approach means that there are no distributional assumptions (including equality of sample variances).
In order to assess the significance of subthreshold responses to cVA we first smoothed the membrane potential with a 45ms boxcar filter. We then calculated the peak response in a window from +200 to +800 ms after valve opening, subtracting a baseline value of the mean response between −1000 to 0 ms. We repeated the same procedure for a response window preceding odor presentation, generating a null distribution of pseudo responses from all of the traces available for each cell. The cVA responses were then compared with the pseudo responses using the non parametric oneway_test of the R coin package. Raw p values were adjusted to control the False Discovery Rate using the p.adjust function in R; responses were declared significant for FDR adjusted p < 0.01.
Odor Stimulation
Odorants were kept in glass vials with glass wool (Acros Organics) to increase the effective surface area, thus facilitating odor headspace formation in the vial. A constant airstream of 2 l min-1 was directed at the fly throughout the recording. This airstream was composed of a carrier stream (1.75 l min-1) and an odor stream (0.25 l min-1) that was directed from a mineral oil control vial to an odor vial by a trigger-controlled solenoid valve. The rate between carrier and odor stream was set via two mass flow controllers (MFCs). The two streams were joined 70 cm from the end of the delivery tube, which measured 4 mm in diameter and was positioned 8 mm from the fly. All odor concentrations, as they leave the delivery tube, are therefore ∼800x dilutions. Typical odor response latencies with this setup were 150–200 ms, as measured by spiking onset in our recordings and with a photoionization detector (mini-PID, Aurora Scientific) that was placed near the end of the delivery tube. In addition, we sometimes observed a quasi-instantaneous spiking onset in neurons (latency < 5 ms) that we ascribe to a change in pressure.
Several precautions were taken to account for false negative results, i.e., recordings from flies that were unable to respond to odors per se: (a) successful odor delivery was monitored with a photoionization detector in series with the animal (see above), (b) great care was taken during the preparation to fully expose the olfactory appendages of the fly without damaging them or the associated olfactory nerves, (c) in order to determine whether neurons were capable of generating action potentials (“spikes”), positive current was injected, and (d) neurons were inspected to establish if they exhibited at least one supra- or subthreshold response to odor. Recordings were only included if criteria (a–c) were fulfilled; if (d) was not fulfilled, a field recording (Figure S2C) was performed to ensure that the animal was odor-responsive (see above).
Single Glomerulus Stimulation
For single glomerulus stimulations, theta glass stimulation electrodes (Warner Instruments) were pulled with a 3–5 μm tip and filled with 2 mM acetylcholine chloride (Sigma) in saline. With the exception of the experiment in Figure S3D, both sides of the theta electrode were filled with identical solution (i.e., effectively a one-barreled electrode), using a separate bath electrode as ground. While there is therefore no liquid junction potential between the two barrels of the theta electrode, there is a liquid junction potential between the bath electrode (in saline) and the theta electrode (saline + acetylcholine). However, since the concentration of acetylcholine used in these experiments was very low (2 mM), this liquid junction potential is less than 0.2 mV.
Stimulation electrodes were inserted into the GFP-positive DA1 glomerulus. These recordings were performed in flies with an additional copy of Mz19-Gal4 or Or67dGal4 in order to reliably identify the DA1 glomerulus. In order to reduce spontaneous activity in DA1 PNs, both antennae were removed with forceps. Square voltage pulses (500 ms), ranging from 0.2–10 V were delivered via the DAC output of the ITC-18 data acquisition interface (HEKA Elektronik) using a common bath electrode. This configuration results in both transient and sustained stimulation artifacts; the sustained artifacts (typically 1 mV/V stimulation) were fitted for each cell and subtracted before quantification. Stimulation depolarizes PNs both chemically and electrically (see Figure S3D). For control stimulations the stimulation electrode was moved to neighboring glomeruli in the dorsal antennal lobe. Note that although this controls for stimulation spreading ’across’ glomeruli, it does not rule out input from more ventral glomeruli to (broadly odor-tuned) female aSP-g neurons. After control stimulations of non-DA1 glomeruli, the stimulation electrode was repositioned in DA1 to ensure that there was no drop in responsivity throughout the recording procedure. A small number of experiments investigated the minimum stimulation voltage required to elicit a postsynaptic response, using stimulation intensities as low as 10 mV.
To measure stimulation latency we first applied a digital boxcar filter at 2.2 kHz to the signals. We normalized the voltage level by subtracting the mean value between 0 and 0.5 ms after stimulation. Responses to both DA1 and control glomerulus stimulation were plotted (Figure S3B) and binned in 0.5 ms intervals. We then tested for divergence between voltage traces for DA1 and control glomerulus stimulation by carrying out a Wilcoxon rank sum test within each 0.5 ms bin. The latency was then defined as the time after which the DA1 trace remained significantly larger (p < 0.05, one-tailed) than the control trace for the remainder of the stimulation period.
References 60
- 1.
Baker B.S., Ridge K.A. Sex and the single cell. I. On the action of major loci affecting sex determination in Drosophila melanogaster. Genetics. 1980;94:383-423.
Not in our library10.1093/genetics/94.2.383
- 2.
Baker B.S., Taylor B.J., Hall J.C. Are complex behaviors specified by dedicated regulatory genes? Reasoning from Drosophila. Cell. 2001;105:13-24.
Not in our library10.1016/s0092-8674(01)00293-8
- 3.
Benjamini Y., Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Stat. Methodol. 1995;57:289-300.
Not in our library
- 4.
Billeter J.C., Rideout E.J., Dornan A.J., Goodwin S.F. Control of male sexual behavior in Drosophila by the sex determination pathway. Curr. Biol. 2006;16:R766-R776.
Not in our library10.1016/j.cub.2006.08.025
- 5.
Billeter J.C., Villella A., Allendorfer J.B., Dornan A.J., Richardson M., Gailey D.A., Goodwin S.F. Isoform-specific control of male neuronal differentiation and behavior in Drosophila by the fruitless gene. Curr. Biol. 2006;16:1063-1076.
Not in our library10.1016/j.cub.2006.04.039
- 6.
Butler B., Pirrotta V., Irminger-Finger I., Nöthiger R. The sex-determining gene tra of Drosophila: molecular cloning and transformation studies. EMBO J. 1986;5:3607-3613.
Not in our library10.1002/j.1460-2075.1986.tb04689.x
- 7.
- 8.
Chiang A.S., Lin C.Y., Chuang C.C., Chang H.M., Hsieh C.H., Yeh C.W., Shih C.T., Wu J.J., Wang G.T., Chen Y.C. Three-dimensional reconstruction of brain-wide wiring networks in Drosophila at single-cell resolution. Curr. Biol. 2011;21:1-11.
Not in our library10.1016/j.cub.2010.11.056
- 9.
Clyne J.D., Miesenböck G. Sex-specific control and tuning of the pattern generator for courtship song in Drosophila. Cell. 2008;133:354-363.
Not in our library10.1016/j.cell.2008.01.050
- 10.
Couto A., Alenius M., Dickson B.J. Molecular, anatomical, and functional organization of the Drosophila olfactory system. Curr. Biol. 2005;15:1535-1547.
Not in our library10.1016/j.cub.2005.07.034
- 11.
Datta S.R., Vasconcelos M.L., Ruta V., Luo S., Wong A., Demir E., Flores J., Balonze K., Dickson B.J., Axel R. The Drosophila pheromone cVA activates a sexually dimorphic neural circuit. Nature. 2008;452:473-477.
Not in our library10.1038/nature06808
- 12.
Demir E., Dickson B.J. fruitless splicing specifies male courtship behavior in Drosophila. Cell. 2005;121:785-794.
Not in our library10.1016/j.cell.2005.04.027
- 13.
Dickson B.J. Wired for sex: the neurobiology of Drosophila mating decisions. Science. 2008;322:904-909.
Not in our library10.1126/science.1159276
- 14.
Dulac C., Kimchi T. Neural mechanisms underlying sex-specific behaviors in vertebrates. Curr. Opin. Neurobiol. 2007;17:675-683.
Not in our library10.1016/j.conb.2008.01.009
- 15.
Ejima A., Smith B.P.C., Lucas C., van der Goes van Naters W., Miller C.J., Carlson J.R., Levine J.D., Griffith L.C. Generalization of courtship learning in Drosophila is mediated by cis-vaccenyl acetate. Curr. Biol. 2007;17:599-605.
Not in our library10.1016/j.cub.2007.01.053
- 16.
Evers J.F., Schmitt S., Sibila M., Duch C. Progress in functional neuroanatomy: precise automatic geometric reconstruction of neuronal morphology from confocal image stacks. J. Neurophysiol. 2005;93:2331-2342.
Not in our library10.1152/jn.00761.2004
- 17.
Fisek M., Wilson R.I. Stereotyped connectivity and computations in higher order olfactory neurons. Nat. Neurosci. 2013.
Not in our library10.1038/nn.3613
- 18.
Fishilevich E., Vosshall L.B. Genetic and functional subdivision of the Drosophila antennal lobe. Curr. Biol. 2005;15:1548-1553.
Not in our library10.1016/j.cub.2005.07.066
- 19.
Grosjean Y., Rytz R., Farine J.P., Abuin L., Cortot J., Jefferis G.S.X.E., Benton R. An olfactory receptor for food-derived odours promotes male courtship in Drosophila. Nature. 2011;478:236-240.
Not in our library10.1038/nature10428
- 20.
Ha T.S., Smith D.P. A pheromone receptor mediates 11-cis-vaccenyl acetate-induced responses in Drosophila. J. Neurosci. 2006;26:8727-8733.
Not in our library10.1523/jneurosci.0876-06.2006
- 21.
Haga S., Hattori T., Sato T., Sato K., Matsuda S., Kobayakawa R., Sakano H., Yoshihara Y., Kikusui T., Touhara K. The male mouse pheromone ESP1 enhances female sexual receptive behaviour through a specific vomeronasal receptor. Nature. 2010;466:118-122.
Not in our library10.1038/nature09142
- 22.
Horowitz P., Hill W. The Art of Electronics. 1989.
Not in our library
- 23.
Hothorn T., Hornik K., Van De Wiel M.A., Zeileis A. A lego system for conditional inference. Am. Stat. 2006;60:257-263.
Not in our library
- 24.
Ito H., Fujitani K., Usui K., Shimizu-Nishikawa K., Tanaka S., Yamamoto D. Sexual orientation in Drosophila is altered by the satori mutation in the sex-determination gene fruitless that encodes a zinc finger protein with a BTB domain. Proc. Natl. Acad. Sci. USA. 1996;93:9687-9692.
Not in our library10.1073/pnas.93.18.9687
- 25.
Ito H., Sato K., Koganezawa M., Ote M., Matsumoto K., Hama C., Yamamoto D. Fruitless recruits two antagonistic chromatin factors to establish single-neuron sexual dimorphism. Cell. 2012;149:1327-1338.
Not in our library10.1016/j.cell.2012.04.025
- 26.
Jang H., Kim K., Neal S.J., Macosko E., Kim D., Butcher R.A., Zeiger D.M., Bargmann C.I., Sengupta P. Neuromodulatory state and sex specify alternative behaviors through antagonistic synaptic pathways in C. elegans. Neuron. 2012;75:585-592.
Not in our library10.1016/j.neuron.2012.06.034
- 27.
- 28.
Kazama H., Wilson R.I. Homeostatic matching and nonlinear amplification at identified central synapses. Neuron. 2008;58:401-413.
Not in our library10.1016/j.neuron.2008.02.030
- 29.
Kimchi T., Xu J., Dulac C. A functional circuit underlying male sexual behaviour in the female mouse brain. Nature. 2007;448:1009-1014.
Not in our library10.1038/nature06089
- 30.
Kimura K.I., Ote M., Tazawa T., Yamamoto D. Fruitless specifies sexually dimorphic neural circuitry in the Drosophila brain. Nature. 2005;438:229-233.
Not in our library10.1038/nature04229
- 31.
Kimura K.I., Hachiya T., Koganezawa M., Tazawa T., Yamamoto D. Fruitless and doublesex coordinate to generate male-specific neurons that can initiate courtship. Neuron. 2008;59:759-769.
Not in our library10.1016/j.neuron.2008.06.007
- 32.
Kohatsu S., Koganezawa M., Yamamoto D. Female contact activates male-specific interneurons that trigger stereotypic courtship behavior in Drosophila. Neuron. 2011;69:498-508.
Not in our library10.1016/j.neuron.2010.12.017
- 33.
Kurtovic A., Widmer A., Dickson B.J. A single class of olfactory neurons mediates behavioural responses to a Drosophila sex pheromone. Nature. 2007;446:542-546.
Not in our library10.1038/nature05672
- 34.
Lee T., Luo L. Mosaic analysis with a repressible cell marker for studies of gene function in neuronal morphogenesis. Neuron. 1999;22:451-461.
Not in our library10.1016/s0896-6273(00)80701-1
- 35.
Lee G., Foss M., Goodwin S.F., Carlo T., Taylor B.J., Hall J.C. Spatial, temporal, and sexually dimorphic expression patterns of the fruitless gene in the Drosophila central nervous system. J. Neurobiol. 2000;43:404-426.
Not in our library10.1002/1097-4695(20000615)43:4<404::aid-neu8>3.0.co;2-d
- 36.
Luan H., Peabody N.C., Vinson C.R., White B.H. Refined spatial manipulation of neuronal function by combinatorial restriction of transgene expression. Neuron. 2006;52:425-436.
Not in our library10.1016/j.neuron.2006.08.028
- 37.
Manoli D.S., Foss M., Villella A., Taylor B.J., Hall J.C., Baker B.S. Male-specific fruitless specifies the neural substrates of Drosophila courtship behaviour. Nature. 2005;436:395-400.
Not in our library10.1038/nature03859
- 38.
Masse N.Y., Cachero S., Ostrovsky A.D., Jefferis G.S.X.E. A mutual information approach to automate identification of neuronal clusters in Drosophila brain images. Front. Neuroinform. 2012;6:21.
Not in our library10.3389/fninf.2012.00021
- 39.
Mellert D.J., Knapp J.M., Manoli D.S., Meissner G.W., Baker B.S. Midline crossing by gustatory receptor neuron axons is regulated by fruitless, doublesex and the Roundabout receptors. Development. 2010;137:323-332.
Not in our library10.1242/dev.045047
- 40.
Nojima T., Kimura K.i., Koganezawa M., Yamamoto D. Neuronal synaptic outputs determine the sexual fate of postsynaptic targets. Curr. Biol. 2010;20:836-840.
Not in our library10.1016/j.cub.2010.02.064
Cited by papers we hold 10
- Information flow, cell types and stereotypy in a full olfactory connectome · 2021No notice
- Achieving functional neuronal dendrite structure through sequential stochastic growth and retraction · 2020Not checked
- Complete Connectomic Reconstruction of Olfactory Projection Neurons in the Fly Brain · 2020No notice
- The natverse, a versatile toolbox for combining and analysing neuroanatomical data · 2020No notice
- Neurogenetic dissection of the Drosophila lateral horn reveals major outputs, diverse behavioural functions, and interactions with the mushroom body · 2019No notice
- NeuroMatic: An Integrated Open-Source Software Toolkit for Acquisition, Analysis and Simulation of Electrophysiological Data · 2018No notice
- P1 interneurons promote a persistent internal state that enhances inter-male aggression in Drosophila · 2015No notice
- The Olfactory Transcriptomes of Mice · 2014No notice
- Decoding odor quality and intensity in the Drosophila brain · 2014No notice
- The neuronal architecture of the mushroom body provides a logic for associative learning · 2014No notice
Reads
Read by people on this page
None yet.
Checked or read by AI through naulon
None yet.
A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.
Are you an author?
Is this your paper?
Claim it to add your note and be paid when agents cite it. Free.