|Home | About | Journals | Submit | Contact Us | Français|
Local interneurons play essential roles in information processing by neural circuits. Here we present a comprehensive genetic, anatomical, and electrophysiological analysis of local interneurons (LNs) in the Drosophila antennal lobe, the first olfactory processing center in the brain. We find that LNs are diverse in their neurotransmitter profiles, connectivity, and physiological properties. Analysis of >1500 individual LNs reveals major morphological classes characterized by coarsely stereotyped glomerular innervation patterns. Some of these morphological classes exhibit distinct physiological properties. However, the finer-scale connectivity of an individual LN varies considerably across brains and there is notable physiological variability within each morphological or genetic class. Finally, we show that LN innervation requires interaction with olfactory receptor neurons during development, and some individual variability also likely reflects LN-LN interactions. Our results reveal an unexpected degree of complexity and individual variation in an invertebrate neural circuit, a result that creates challenges for solving the Drosophila connectome.
Neurons can be divided into two general categories: projection neurons whose axons connect discrete regions of neural tissue, and local interneurons whose processes are restricted to a single region. Local interneurons play an important role in neural computation. For example, horizontal and amacrine cells in the vertebrate retina are instrumental in transforming incoming signals from photoreceptors, such that information regarding contrast, motion and color can be encoded by retinal ganglion cells and transmitted into the brain. Every vertebrate brain region contains a multitude of local interneurons, as defined by morphology, histological markers, and physiology, and likely distinct functions1. However, it is less clear whether the same kind of diversity exists in simpler organisms such as Drosophila. Systematic identification and characterization of local interneurons in the fly is a prerequisite to understanding information processing in neural circuits2.
The fly antennal lobe is an attractive model neural circuit because of its elegant anatomical and functional organization and its genetic tractability. All olfactory receptor neurons (ORNs) that express the same odorant receptor project their axons to the same antennal lobe glomerulus, and most antennal lobe projection neurons (PNs) send dendrites into a single glomerulus to receive direct input from one ORN class3. Genetic access is available to many classes of ORNs and PNs3. In contrast to these input and output neurons, local interneurons (LNs) in the fly antennal lobe are much less well characterized.
In the olfactory bulb, the vertebrate analog of the antennal lobe, there are ~100 local interneurons for each projection neuron4, one of the highest local interneuron/projection neuron ratios in any mammalian brain region. Olfactory bulb interneurons are morphologically diverse and are hypothesized to play important roles in gain control, spike synchronization, and decorrelating the representations of similar odorants4-6. Studies of antennal lobe LNs in larger insects have described a variety of LN morphologies7-11. In comparison, the limited datasets available suggested that Drosophila antennal lobe LNs have relatively simple morphologies, with each LN innervating all glomeruli, although LNs which are not pan-glomerular have been described in more recent studies12-18. Understanding how ORN input is transformed into PN output by the antennal lobe requires a deeper understanding of antennal lobe LNs.
Here we take a systematic approach to identify and characterize LNs in the Drosophila antennal lobe. Using specific Gal4 lines in conjunction with the MARCM (Mosaic Analysis with Repressible Cell Marker) method19 and dye fills, we analyzed the morphology of >1,500 LNs at a single cell level. This analysis reveals a striking diversity of glomerular innervation patterns, which can be used to define distinct morphological classes of LNs. Major morphological LN classes have distinctive electrophysiological properties that help define their potential functions within the circuit. At a fine scale, our analysis also demonstrates that the innervation pattern of an individual LN varies across brains. Similarly, even within a morphologically- or genetically-defined class of LNs, we find considerable variability in physiological properties. This individual variation may in part reflect the role of cell-cell interactions during development. Our findings imply that LNs play unexpectedly diverse roles in olfactory processing in the fly antennal lobe. The wiring variability of LNs suggests that the wiring diagram differs considerably between individual fly brains, and that individual variations in behavior may in part reflect these kinds of fine-scale variations at the cellular level.
We used 10 Gal4 lines to genetically access antennal lobe LNs. The expression patterns of these lines span a range from hundreds of neurons near the antennal lobe to just a handful of LNs (Fig. 1). For simplicity, we refer to these Gal4 lines as Lines 1-10 (see Methods). In addition to LNs—which we define as neurons whose cell bodies are near the antennal lobe and whose processes are restricted to the antennal lobes—several lines also label ORNs, PNs and/or neurons near the antennal lobe that project to other areas of the brain (Supplementary Fig. 1 online).
LN cell bodies are located in two clusters: a large but continuous cluster lateral and dorsolateral to the antennal lobe (most LNs labeled by Lines 1, 3-9), and a separate cluster ventral to the antennal lobe (most LNs labeled by Line10) (Fig. 1, Supplementary Fig. 1 online). To determine the number and the potential overlap of LNs labeled by these Gal4 lines, we counted the number of nuclei labeled by UAS-nuclear marker driven from individual Gal4s, as well as from their combinations in the same fly (Supplementary Table 1 online). These data suggest that some lines label largely non-overlapping LN populations, whereas other lines overlap partly or completely.
Line5 appears to include all cells that are labeled by Lines 6-9. All ~56 cells labeled by Line5 are LNs because 1) Line5 driven mCD8-GFP does not label any neurons that send processes out of the antennal lobe (Fig. 1), and 2) all 578 single cells from systematic MARCM analysis and biocytin fill (see below) from Line5 have their entire arborization within the antennal lobe. Based on the total number of LNs labeled by Lines 3, 4 and 5 in the same fly (Supplementary Table 1 online) and cell body positions from neuroblast clone analysis (Supplementary Fig. 1 online), we estimate that there are ~100 LNs labeled by Lines 1 and 3-9, which are mostly ipsilaterally projecting LNs in the lateral cluster (see below). Line10 labels mostly bilaterally projecting LNs (see below) in the ventral cluster. Therefore, the lower bound for Drosophila LNs (labeled by our Gal4 lines) is ~100 ipsilaterally projecting and ~100 bilaterally projecting LNs for each antennal lobe.
Although LNs are traditionally considered to be inhibitory neurons that release GABA as their neurotransmitter, recent studies have suggested that some LNs can be excitatory14, 20 and cholinergic14, 17. We examined the neurotransmitter type of different Gal4 lines by co-staining with antibodies against GABA and Choline acetyltransferase. Antibody specificity was validated by loss of staining in MARCM clones homozygous mutant for Glutamic acid decarboxylase 1 and Choline acetyltransferase (Supplementary Fig. 2a-h online). We found that the vast majority of LNs are GABAergic; however, there are a few cholinergic cells, and a larger minority that are neither GABAergic nor cholinergic (Supplementary Table 1 online; Supplementary Fig. 2i-l online). Staining for additional neurotransmitter types identifies a few dopaminergic and a large number of (mostly ventral) glutamatergic LNs (Supplementary Fig. 3 online). These results indicate that although the large majority of LNs are GABAergic, some LNs use different neurotransmitters or their combinations, emphasizing the potential diversity in the functional roles of different LNs.
The MARCM method can be used to label single neurons born at specific developmental times21, 22. We performed systematic MARCM labeling of LNs from these 10 Gal4 lines throughout development. Figure 2 shows representative images from 1,439 single cell LN clones. As is evident, LNs exhibit remarkable morphological diversity.
We identified 5 broad categories of ipsilaterally projecting LNs based on the extent of their arborization in the antennal lobe: the entire antennal lobe (pan-glomerular; Fig. 2a), all but a few glomeruli (Fig. 2b), a continuous region of the antennal lobe (Fig. 2c), discontinuous or patchy regions of the antennal lobe (Fig. 2d), 1-3 glomeruli (few-glomerular; Fig. 2e). Within each category, however, there is considerable diversity. For example, although many LNs arborize nearly the entire antennal lobe, they differ in the density of arborization and the thickness of their processes (Fig. 2a-b; Supplementary Fig. 4a-c online). For LNs that skip a few glomeruli or arborize in a continuous region, different individuals skip different glomeruli or arborize in different regions of the antennal lobe (Fig. 2b-c; Supplementary Fig. 5 online). Finally, bilaterally projecting LNs exhibit heterogeneous ipsilateral and contralateral arborization patterns (Fig. 2f).
Since all neurotransmitter release of LNs occurs within the antennal lobe, does each LN releases neurotransmitter throughout its arbor, or only in a restricted zone? To examine how presynaptic terminals are distributed within individual LNs, we performed a large subset of MARCM single cell labeling with an additional synaptotagmin-HA marker that labels presynaptic terminals23. In general, synaptotagmin-HA labels puncta throughout an LN's processes (red staining in Fig. 2, Supplementary Fig. 4 online). These data suggest that as a general rule, LNs broadcast transmitter release across their arborization without notable glomerular selectivity. However, in exception to this rule, some few-glomerular and bilateral LNs have a more selective synaptotagmin-HA distribution (Fig. 2e-f, Supplementary Fig. 4f online). These LNs therefore appear to deliver information to only a subset of glomeruli they innervate.
Given the diversity of morphologies, we next analyzed these data systematically with the aim of identifying organizational principles of the LN network. An important feature of an LN is its glomerular innervation pattern, as this determines which olfactory channels the LN may receive information from and send information to. We therefore started by analyzing the glomerular innervation patterns of individual LNs.
We scored 1,532 individual LNs (1,439 from MARCM single cell clones; 93 from dye fills during whole cell patch clamp recording) for their innervation of 54 glomeruli in the antennal lobe (Fig. 3b). We employed hierarchical clustering to organize the 1,532 LNs according to their binary ipsilateral glomerular innervation patterns. This analysis revealed several distinct morphological classes (Fig. 3a). For example, neurons within subcluster C tend to avoid certain glomeruli innervated by trichoid ORN classes that sense pheromones24 (Fig. 3c). Neurons within subcluster D innervate the central antennal lobe and skip the dorsal and ventral glomeruli; we call these dumbbell cells after their shape (Fig. 3d). Neurons within subcluster E exhibit patchy innervation (Fig. 3e).
Clustering all cells labeled by the same Gal4 (Supplementary Fig. 6 online) reveals that some Gal4 lines label a relatively uniform morphological population of cells (e.g. Lines 6-8), whereas others label a more heterogeneous group (e.g. Lines 1, 2, 3, and 5). This suggests that some Gal4 lines predominantly label LNs belonging to just a few functional classes, whereas others label a functionally heterogeneous population. We will examine this idea in more detail using electrophysiological recordings (see below).
Clustering all GABA+ vs. GABA- cells from our collection revealed that both subtypes include a diversity of morphological patterns (Supplementary Fig. 7 online). Because most LNs are GABAergic, we considered the hypothesis that LN innervation might be denser in glomeruli that receive more ORN input. We estimated the average lifetime activity of each ORN type based on a published dataset25 by averaging the responses of each ORN type across 110 odor stimuli. We found a significant positive correlation between LN innervation probability and mean odor-evoked ORN firing rate (Fig. 3f). This would be consistent with the idea that the density of inhibitory innervation is adapted to the average level of afferent input to each glomerulus26, assuming that these stimuli are representative of natural odors.
In addition to hierarchical clustering, we used principal component analysis (PCA) to identify emergent relationships in the glomerular innervation dataset without any arbitrary threshold of significance. The first and second principal components (PC1 and PC2; Fig. 3g and Supplementary Fig. 8 online) for the entire set of innervation patterns accounted for 40% and 10% of the variance in the data, respectively, but 13 additional principal components were needed to account for an additional 25%. The high number of dimensions necessary to account for variations indicates that innervation patterns are not well described as the linear combination of any small number of basic glomerular relationships. PC1 is essentially a proxy for the number of glomeruli innervated (correlation coefficient=0.9993). Thus, a defining feature of an LN's innervation pattern is the number of glomeruli it connects to. LNs were bimodally distributed on PC2, with a distinctive subcluster containing mainly dumbbell cells (Fig. 3d), confirming that the dumbbell cells constitute a distinctive morphological class.
Next, we considered the LNs that arborize bilaterally. We compared ipsilateral innervation patterns (Supplementary Fig. 9a online) with corresponding contralateral innervation patterns (Supplementary Fig. 9b online). We found that many bilateral LNs have rather symmetric glomerular innervation patterns and overall, ipsi- and contralateral patterns were significantly correlated (Pearson's correlation coefficient, r=0.46±0.05 (SEM), p<0.01, n=38) (Supplementary Fig. 9c online). In general, the exceptions were LNs with nearly pan-glomerular ipsilateral innervation.
Finally, based on clustering we identified a minimal number of LN classes with identifiable innervation patterns and further investigated birth timing and Gal4 expression for those cell types (Supplementary Fig. 5 online, Supplementary Table 2 online). We found that some (although not all) morphological types of LNs are labeled by only a few Gal4 lines and are born in limited time windows. Together, these analyses indicate that LNs are composed of genetically distinct groups of cells defined by a combination of morphology, birth timing and Gal4 expression.
Do LN genetic categories (Gal4 lines) and morphological categories (glomerular innervation patterns) have any correlation with physiological properties? To address this, we selected five Gal4 lines that label relatively small subsets of LNs (Lines 5-9) and made whole-cell patch-clamp recordings from GFP-labeled somata in each line. We recorded the responses of each cell to a panel of diverse odors (see Methods), and filled each cell with biocytin via the patch pipette for visualization after recording. We successfully filled and scored the innervation patterns of 93 LNs.
All these LNs fired spontaneous action potentials, and their spiking was always modulated by odors. LN odor responses were remarkably diverse and typically varied more across cells than across odors within a cell (Fig. 4a,b). Some LNs recorded in the same Gal4 line responded to odors in a similar manner (Fig. 4a). For example, the LNs we recorded in Line7 were relatively stereotyped in their functional properties, as were Line8 LNs. In other cases, a Gal4 line could label LNs having very different response profiles (Fig. 4b). Line5 and Line9 were particularly diverse in this respect.
In all these LNs, we measured spontaneous firing rate, mean odor response, and maximum odor response. To quantify the time course of the odor response, we also measured the percentage of spikes fired during the first 100 msec of the odor response. All four of these properties were significantly dependent on Gal4 lines (Fig. 4c). This implies that differences in Gal4 expression reflect physiological differences between LNs, in spite of the variability within each Gal4 line and the partial overlap in the cells that are labeled by these lines.
Next, we asked whether LN physiological properties are correlated with their morphological class. Because the principal axis for morphological variation is the number of glomeruli each LN innervates (Fig. 3g), we began by asking whether any response properties were correlated with this axis. Unexpectedly, LNs that innervate many glomeruli have lower odor-evoked firing rates than LNs that innervate fewer glomeruli. The number of glomeruli innervated by an LN is significantly (although weakly) anti-correlated with the average strength of its odor responses (Pearson's r=−0.2, p<0.05).
This finding motivated us to examine the physiology of pan-glomerular LNs in particular, because these cells innervate the largest number of glomeruli. This class comprised 28% of the cells we recorded from. On average, pan-glomerular LNs had significantly higher spontaneous firing rates than other LNs (Fig. 5b). In the presence of an odor, spontaneous spiking in many pan-glomerular cells shut down completely, sometimes after a brief burst at odor onset, while others modestly increased their firing rate in the presence of odors (Fig. 5a). Overall, odor-evoked changes in firing rate were significantly weaker in pan-glomerular cells that in other LNs (Fig. 5b).
A second major class of LNs that we recorded from were LNs that selectively avoid glomerulus VA1d and frequently also avoid DL3 (Fig. 3c). These glomeruli are innervated by trichoid ORNs thought to be selective for pheromones24. This innervation pattern (with slight variations) accounts for ~15% of all LNs in our data set from Lines 1-7, notably those labeled by Line6 (Fig. 3a; see below), and 10% of cells we filled in electrophysiological recordings. These cells differed from other LNs in having especially transient bursts of excitation at odor onset (Fig. 5c). Overall, this morphological class fired a higher percentage of their spikes in the first 100 msec of the response, as compared to other LNs (Fig. 5d). Thus, these LNs could create a transient pulse of GABAergic inhibition at odor onset. If so, then some pheromone glomeruli would evidently be selectively excluded from this transient pulse of inhibition.
Does the diversity we observe in LN morphology and physiology arise from variability of the same LNs across different individual animals, or from many diverse LN types each with stereotypic patterns in all individuals? The first alternative is implied by our finding of 847 distinct innervation patterns in 1,489 ipsilateral-projecting LNs (Fig. 3a), far exceeding the number of estimated total ipsilateral-projecting LNs (~100) within an individual antennal lobe. LN arborization patterns therefore cannot be completely stereotyped across animals.
We sought to address this question by dense sampling of a Gal4 line (Line6) that labels a small population (7 LNs per antennal lobe, Supplementary Table 1 online). Of 131 labeled single cells for Line6, we observed 76 distinct glomerular innervation patterns (Fig. 6a). Since there are many more innervation patterns than the number of cells per animal, individual Line6 labeled LNs cannot be identical across animals.
Electrophysiological recordings showed that Line6 LNs also have diverse functional properties (Fig. 6b-e and and5c).5c). Some cells (Fig. 6b,c) showed a transient burst at the onset of almost every odor, followed by inhibition. Other cells (Fig. 6d,e) exhibited sustained excitation, off-excitation, and more odor-specific tuning. Thus, neither the morphology nor the physiology of Line6 LNs is highly stereotyped.
Despite the clear evidence that individual LNs labeled by Line6 are not rigidly stereotyped, the properties of these LNs are far from randomly distributed. For example, the odor response properties of Line6 cells are significantly less varied across cells than the odor response properties of all LNs. Mean and maximum odor-evoked firing rates were less variable within the Line6 population (n=16) than across the LN population as a whole (n=92, p<5×10-5, F-tests). In addition, spontaneous firing rates were also less variable within the Line6 population than across all LNs (p<5×10-6, F-test).
The glomerular innervation patterns of Line6 LNs are also far from random. Each Line6 LN innervates 51.1 glomeruli (n=131, SEM=0.22) on average. If the innervation pattern were completely random, then every glomerulus would have a 94.4% probability of being innervated. Assuming a binomial probability distribution, the number of samples innervating a glomerulus should vary by the standard error, ±2%. The actual distribution clearly deviates from this prediction (Fig. 6f). Most glomeruli are almost always innervated, whereas a subset of glomeruli is often missed. Among this latter group of glomeruli, VA1d is most often omitted. Indeed, most of the biocytin-filled LNs that selectively avoided VA1d (Fig. 5d) were Line6 LNs (7 of 9).
Since our binary scores for innervation patterns are only a coarse measure of a LN's morphology, we additionally measured the innervation density of Line6 LNs for selected glomeruli that are either always innervated or often missed. We found that DM1 (always innervated) has a significantly higher innervation density compared to VA1d or DA1 (occasionally innervated), whereas randomly selected pan-glomerular LNs exhibit similar innervation density of DM1, VA1d and DA1. Innervation density of Line6 neurons in DA1 and VA1d glomeruli is significantly reduced compared with that of pan-glomerular LNs (Fig. 6g). These data strengthen the distinction of these two groups of glomeruli with regard to Line6 LN innervation. One group is always innervated, while the other group has highly reduced innervation if not completely avoided.
In summary, analysis of Line6 indicated a principle likely applicable to most other morphologically and genetically-defined classes of LNs. LNs exhibit marked variability within a class, and some of this variability must represent variations in the same cells across different brains. Yet the properties of these LNs are clearly not drawn randomly from the entire distribution of LN properties. This picture is consistent with the idea that the coarse properties of these cells are genetically pre-programmed, whereas their finer-scale properties may also reflect factors such as developmental plasticity and sensory experience.
The variability of glomerular innervation patterns is exemplified in a small subset of LNs from Lines 1, 3, and 5 that exhibit a distinctive pattern of patchy innervation (Fig. 2d, ,3e).3e). Strikingly, comparison of the glomerular innervation patterns of 161 patchy cells (Fig. 7a) indicates that none of these single cells innervate identical sets of glomeruli. Given that these patchy cells represent only a small fraction of cells labeled by LN Gal4s and a small fraction of the ~100 ipsilaterally projecting LNs, it is clear that the innervation patterns of patchy LNs are highly variable across different animals.
The patchy nature of these cells suggested that these innervation patterns might be established through cell-cell interactions among LNs. To examine the relationships between innervation patterns of different LNs in the same animal, we designed a genetic method to visualize a pair of sister patchy cells with two different colors in the same animal (Fig. 7b). Although the efficiency of this method was extremely low, yielding only 5 pairs of patchy cells after thousands of dissected brains, it is evident that processes of sister patchy cells avoid each other to carve out their innervation domains, sometimes by splitting a single glomerulus (outlined in Fig. 7c-e, Supplementary Fig. 10 online). By contrast, differentially labeled sister cells of other LN types inter-mix at much finer scales (Fig. 7f-h). We confirmed the qualitative difference by quantifying the degree of overlap between the signals of the two cells after the signals had been systematically dilated to varying degrees (see Methods). Non-patchy cells exhibited a great deal of overlap with just a small amount of dilation, while patchy cells exhibited similar degrees of overlap only after they had been dilated much more (Fig. 7i). These observations suggest that sister patchy LNs may tile the antennal lobe the way ganglion cells tile the vertebrate retina27 and sensory neurons tile the Drosophila body wall28. We also estimated that each patchy LN occupies 13.1±1.6% (mean±SD, n=8) of the antennal lobe volume, implying that eight patchy LNs could tile an entire antennal lobe. These findings suggest that LN-LN interactions play a role in establishing patchy LN morphologies.
Finally, we tested the contributions of ORNs to the glomerular arborization pattern of LNs. During the construction of the antennal lobe circuit, PN dendrites initiate pattern formation by targeting dendrites to positions similar to where they are found in the adult antennal lobe before the arrival of pioneering axons. ORN axon invasion into the antennal lobe and subsequent ORN axon-PN dendrite recognition results in the formation of distinguishable glomeruli29, 30. The development of LNs has not previously been described.
We first tested whether LN arborization during pupal development requires ORN axons. We made use of the fact that Hedgehog signaling is required for antennal disc proliferation, and clonal loss of Smoothened, a component essential for Hedgehog signaling, results in occasional loss of maxillary palps while leaving the antennae intact30. Under this condition, ORN axons from the antennae would innervate glomeruli that are normal targets of antennal ORN axons, but glomeruli that are normal targets of maxillary palp ORNs are devoid of ORN input30. By introducing Line5 Gal4 driven mCD8-GFP into this genetic background and screening for flies with bilateral loss of maxillary palps (~1 in 250; Fig. 8a), we found that LN processes are not present in the glomeruli that are normal targets of maxillary palp ORNs. LN processes still innervate the nearby glomeruli that are targets of antennal ORNs (Fig. 8b, bottom panels, compared to Fig. 8b, top panels; Fig. 8d). By contrast, PN dendritic processes are mostly still present in glomeruli devoid of ORN axons (Fig. 8c, d). This experiment suggests that glomerular innervation of LNs requires the presence of ORN axons.
To test whether the ORN axons are required for the maintenance of LN processes, we bilaterally removed the maxillary palps in adults, well after antennal lobe wiring is completed. We examined the processes of Line6 LNs (which are a subset of Line5 LNs) at least 10 days after removal of the maxillary palps to allow complete degeneration of ORN axons. We found that LN processes still innervate maxillary palp glomeruli (Fig. 8e). The volume of one of the two glomeruli quantified is significantly reduced (Fig. 8f, top panel), likely because ORN axons contribute to the glomerular volume31. However, the total length of LN processes and the total number of presynaptic terminals as marked by synaptotagmin-HA puncta are unaffected by the adult removal of ORNs (Fig. 8f, middle and bottom panels). We also compared variability of glomerular innervation patterns of Line6 LNs between wild-type flies and those with ORNs removed. We found no significant differences in the number of glomeruli innervated or unique innervation patterns (Supplementary Fig. 11 online), suggesting that variability of Line6 LNs is not dependent on the presence of ORN axons in adult. Taken together, these experiments indicate that ORN axons are essential for LN innervation during development, but are not required for the maintenance of their glomerular innervation in adulthood.
Knowledge of local interneuron properties is essential to understand how a neural circuit processes information. Here we present a systematic single-cell analysis of >1500 antennal lobe LNs in Drosophila. The scale of our analysis is an order of magnitude greater than the most comprehensive analysis of insect LNs published to date8. Our genetic strategy also avoided sampling biases involved with dye filling, which favors large cells with easy access. Despite the caveat that our sampling of LNs is not entirely even (see Methods), our study revealed a striking diversity of antennal lobe LNs beyond that described previously in any insect. LNs are diverse in their neurotransmitter types, glomerular innervation patterns, fine dendritic structures, densities and distribution of presynaptic terminals, and odor response properties. Furthermore, our analysis of glomerular innervation patterns and odor responses revealed coarse stereotypy but also a surprising degree of fine-scale variability in genetically-defined subclasses of LNs. Below we discuss the implications of these findings for the development and function of the antennal lobe circuit, and for the complexity of the Drosophila brain as a model for neural circuit analysis.
The final glomerular innervation patterns of LNs likely result from some combination of: 1) genetic specification, 2) cell-cell interactions during development, and 3) neuronal activity and plasticity. Our study suggests that these mechanisms are used differentially in specifying different aspects of LN diversification. We also cannot rule out additional contributions of sex, age and rearing conditions (see Methods).
We envision that the distinctions between major LN categories are specified largely by genetic mechanisms similar to those used to specify other types of central neurons, such as antennal lobe PNs. These include neuroblast lineage, birth timing, and specific gene expression. For example, most bilateral projecting LNs have cell bodies located ventral to the antennal lobe, are labeled by Line10, and derive from a distinct neuroblast (Supplementary Fig. 1 online). All patchy LNs are born during a well-defined time window in the mid-larval stage (Supplementary Fig. 5 online, Supplementary Table 2 online) and are derived exclusively from Lines 1, 3, and 5. Dumbbell cells derive exclusively from Lines 1, 2, and 9, and are also born during a defined larval stage (Supplementary Fig. 5 online, Supplementary Table 2 online). At a finer level, even though Line6 LNs exhibit a large number of glomerular innervation patterns, these patterns are biased against certain pheromonal glomeruli. Furthermore, the density of arborization and presynaptic terminals within glomeruli that are always innervated is higher than those that are only occasionally innervated by Line6. One can imagine that mechanisms that translate lineage and birth timing to wiring specificity of PNs and other glomerular targeting mechanisms32 can also be used to specify glomerular innervation specificity of certain types of LNs.
Compared to ORNs and PNs, however, LN glomerular innervation patterns are much more variable and less deterministic. The most striking example is patchy LNs: the 161 single cells we examined exhibited 161 distinct glomerular innervation patterns. Our two-cell clone analysis suggests a developmental mechanism for achieving this diversity. We found that processes of sister patchy cells tend to avoid each other even when they innervate common glomeruli. Our quantitative analysis is consistent with the idea that these complex dendrites tile the three-dimensional volume of the antennal lobe. These data suggest that homophilic repulsion might play a role in establishing final innervation patterns of some LN morphological types. From a functional perspective, it seems likely that the specific glomerular innervation patterns of individual patchy LNs are unimportant as long as a group of similar types of patchy LNs cover the entire antennal lobe.
Local interneuron classes defined by morphology and molecular markers are thought to have distinct functions1. The Drosophila antennal lobe represents an unusual opportunity to assess in vivo the functional correlates of LN classes. Overall, we found that LNs varied widely in their odor-evoked firing rates and their temporal dynamics. The selective modulation of LNs with particular dynamics could be a way to modulate the temporal pattern of PN activity.
On a coarse level, some physiological diversity correlated with differences in morphology or genetic markers. For example, pan-glomerular LNs had unusually high rates of spontaneous activity, suggesting that they may be particularly important in broadcasting tonic inhibition throughout the antennal lobe. Tonic inhibition is thought to regulate the sensitivity of odor detection in the mammalian olfactory bulb33, 34. Pan-glomerular LNs were either inhibited by odors or were only weakly excited. Because they are mainly GABAergic, inhibition of pan-glomerular LNs would tend to disinhibit the entire antennal lobe in response to odors.
To take another example, LNs that avoid certain pheromone-selective glomeruli fired unusually transient bursts at odor onset. These LNs may have a special role in inhibiting the initial 100 msec of PN odor responses, when PN firing rates peak35. Some pheromone glomeruli are evidently excluded from this inhibition. This may reflect the fact that pheromones elicit weaker peak firing rates as compared to other odors24, and may help explain why weak pheromone responses are strongly amplified in postsynaptic PNs36.
Although we found a coarse correlation between some functional properties and morphological or genetic classifications, we also found that each LN class we examined in finer detail contained substantial diversity. For example, both pan-glomerular and pheromone-avoiding LNs had varied odor responses. Similarly, Line6 LNs were also heterogeneous, especially considering that there are only ~7 such cells. Thus, we are unlikely to be able to accurately predict every aspect of an LN's physiology based solely on its morphology or a few molecular markers. Interestingly, interneurons of the mammalian cortex are also remarkably varied, prompting the remark that “inhibitory cells seem to be so diverse as perhaps to endanger the notion of cell types in the central nervous system”37. Here we have shown that interneuron diversity in the Drosophila brain rivals that of the mammalian brain, and we have begun to establish developmental principles underlying this diversity.
It is widely believed that, compared to vertebrates38, invertebrate nervous systems are much more stereotyped across different members of the same species, with individually identifiable neurons and stereotyped connection patterns. Indeed, previous analysis of projection patterns of individual antennal lobe PNs provided a striking example in support of this notion39, 40. Thus, the degree of variability of glomerular innervation patterns of individual LNs came as a surprise. Such variability should translate into variability in the functional connectivity of LNs, ORNs and PNs.
While the glomerular organization of the antennal lobe provides a convenient and unequivocal demonstration of variability in the Drosophila “connectome”, it remains to be determined whether variable connectivity is a general feature of the fly brain. Interestingly, even the axons of antennal lobe PNs, which exhibit strikingly stereotyped terminal arborization patterns in the lateral horn, display less stereotypy in the mushroom body23, 39, 40. Indeed, a recent physiology study suggests non-stereotypic connections between PN axons and their postsynaptic targets in the mushroom body41. It is possible that, upon detailed examination, such lack of stereotypy is the norm rather than an exception even in the relatively simple insect brain.
We and others42 believe that a wiring diagram is a necessary (although certainly not sufficient) prerequisite for understanding the function of neural circuits. However, our results imply that the complete reconstruction of the wiring diagram of a single Drosophila brain will not yield a general wiring diagram for all Drosophila brains. At the same time, our findings identify a fundamental similarity between insect and vertebrate brains in the variability of their connectomes. Although it may prove more challenging than previously assumed to map the connectivity of the fly brain, solving the Drosophila connectome may turn out to be more relevant to understanding our own brain.
Line1 (HB4-93) and Line6 (HB8-145) were identified from an anatomical screen of an unpublished Gal4 library from Ulrike Heberlein. Line2 (NP6277) and Line5 (NP3056) were ordered from the NP Gal4 library44 based on our various unpublished screens. Line7 (LCCH3) is an unpublished line from Julie Simpson. Other lines have been described before for their expression in LNs or in other brain neurons: Line3=H2445; Line4=GH29846; Line8=Nan-Gal447; Line9=krasavietz-Gal448; Line10=OK10721.
Our 1,532 single cells from 10 Gal4 lines were collected over the course of two years so experimental conditions vary. Most MARCM analysis used 0-14 day old female flies except Lines 2, 3 and 10 where a subset of experiments were performed with male flies, and Line 9 used exclusively males. A subset of Line6 single cell clones were collected after antennae and/or maxillary palps were removed, so these flies as well as their untreated controls were dissected between 10-20 days after eclosion. All electrophysiology experiments were performed using 2-3 day old females. Flies were reared at 25°C for MARCM experiments and at room temperature for physiology experiments.
It is important to note that the relative frequency of single cells in our collection does not necessarily reflect the fraction of these cells in a given animal because of our sampling bias with regard to Gal4 lines and heatshock windows—the number of MARCM-labeled single cells is not proportional to the number of LNs in a given Gal4 line, and not all developmental windows are equally sampled.
MARCM analysis was performed according to established protocols19. Briefly, MARCM ready stocks were crossed, eggs collected for various lengths of time (ranging from 2-24 hours), and samples heat-shocked at time points of interest during embryonic, larval and early pupal development. Heatshocks at middle and late pupal stages rarely induced LN clones in our early experiments and thus were not employed further. Heat-shock time windows for all cells are listed in Supplementary Table 2. Adult brains were dissected, fixed and stained as described before22. After staining, brains were mounted in SlowFade Gold (Molecular Probes) and imaged on a Zeiss Meta-510 confocal. Images were processed using LSM software (Zeiss) and Image J (http://rsbweb.nih.gov/ij/). Further measurements including glomerular volume, dendrite length and synaptotagmin-HA density (Fig. 6, Fig. 8) were made in Imaris (Bitplane) and statistical analysis performed in Excel and Graphpad Prism 5.0 software. To determine the volume of the antennal lobe occupied by individual patchy cells, LSM software was used to trace the outline of the antennal lobe (marked by N-cadherin) or patchy cell processes (labeled by GFP or CD2) in each confocal section collected at 1 mm intervals. The volume was determined as the sum of areas from each 1 mm section. Cell counts were performed manually or with the aid of Imaris.
Primary antibodies used include: rat anti-mCD8a (1:100, Invitrogen), chicken anti-GFP (1:1000; Aves Labs), rabbit anti-GFP (1:1000, Molecular Probes), mouse anti-nc82 (1:30, Developmental Studies Hybridoma Bank [DSHB]), mouse anti-rCD2 (1:100, OX34, Serotec), rat anti-N-cadherin (1:40, DN-EX #8, DSHB), mouse anti-HA (1:100, 12CA5, gift of Dr. K. Wehner), rabbit anti-HA (1:1000, Abcam), rabbit anti-βgal (1:1000, CAPPEL/MB Biomedicals), chicken anti-βgal (1:1000, Abcam), mouse anti-Cha (1:100, ChAT4B1, DSHB), rabbit anti-GABA (1:200, Sigma), rabbit anti-DVGLUT(1:500; gift of Dr. A. DiAntonio). Secondary antibodies from Molecular Probes (MP) or Jackson ImmunoResearch Laboratories (JIRL) and used at 1:500 include: Zenon mouse Cy5 (MP), goat anti-mouse Cy3 (JIRL), goat anti-rat Alexa488 (MP), goat anti-rat Cy5 (JIRL), goat anti-rat Cy3 (JIRL), goat anti-rabbit Alexa488 (MP), goat anti-rabbit Cy5 (JIRL), donkey anti-chicken FITC (JIRL), donkey anti-chicken Cy3 (1:300, JIRL).
Whole-cell patch-clamp recordings from LN somata were performed as previously described13. External saline contained (in mM): 103 NaCl, 3 KCl, 5 N-tris(hydroxymethyl) methyl-2-aminoethane-sulfonic acid, 8 trehalose, 10 glucose, 26 NaHCO3, 1 NaH2PO4, 1.5 CaCl2, and 4 MgCl2 (osmolarity adjusted to 270-275 mOsm). The saline was bubbled with 95% O2/5% CO2 and reached a pH of 7.3. The internal patch-pipette solution contained (in mM): 140 potassium aspartate, 10 HEPES, 4 MgATP, 0.5 Na3GTP, 1 EGTA, 1 KCl, and 13 biocytin hydrazide. The pH of the internal solution was adjusted to 7.3 and the osmolarity was adjusted to ~ 265 mOsm. LN somata were targeted by means of the GFP signal under epifluorescent illumination in an Olympus BX51WI microscope with a 40X water-immersion objective. Spontaneous firing rates were measured in cell-attached mode prior to whole-cell break-in. After break-in, a small amount of hyperpolarizing current was injected constantly to compensate for the seal conductance and thereby to return the cell to its native spontaneous firing rate49. Recordings were acquired with an A-M Systems Model 2400 amplifier in current clamp mode (10MΩ headstage), low-pass filtered at 2 kHz, and digitized at 10 kHz. One neuron was recorded per brain and the morphology of each cell was visualized post hoc with biocytin-streptavidin and nc82 histochemistry as described previously35 using a Zeiss LSM 510 META confocal microscope.
The total number of LNs recorded in each Gal4 line (and the number whose biocytin fills were successfully reconstructed) was: Line5 55 recorded (51 reconstructed), Line6 18 (17), Line7 21 (11), Line8 8 (5), Line9 11 (8). In some of these cells, we did not complete the full set of physiological measurements, which is why the number experiments cited in connection with some statistical tests is smaller than the numbers listed here. We emphasize that this distribution is not necessarily representative of the actual underlying distribution of these LNs. For example, we deliberately over-sampled Line5 in an attempt to try to hit some of the unusual morphological types labeled by this line (including patchy, dumbbell, and “few-glomerular” LNs). Despite these attempts, we were not able to record from any of these particular morphological classes of LNs, either for technical reasons (e.g., their somata may be unusually small, or they may express low levels of GFP) or potentially due to a low abundance of these morphological classes. Thus, a systematic comparison of all morphological and physiological categories was not possible. The physiological diversity we describe is therefore only a lower bound on the diversity of the LN population.
Odors (Sigma) were delivered with a custom olfactometer, as described previously35. All odors were diluted 1:100 v/v in paraffin oil, and the odor vial headspace was diluted 10-fold in air before reaching the fly. We chose a panel of 10 odors to maximize the chemical diversity of our stimuli, to maximize overlap with odors used in other studies of the same ORNs, and to activate a wide variety of ORN types to varying degrees50. Spike times were extracted and converted to PSTHs by binning them in 50-msec counting windows that overlapped by 25 msec. Each odor was presented 6 times at intervals of 40–60 sec and those trials were then averaged to generate a PSTH for that odor/cell combination. Where the odor response of a cell is expressed as a single time-averaged firing rate, these rates were measured over a 1-sec period beginning at odor onset, and are expressed as a change in firing rate relative to the spontaneous firing rate of that cell. Analyses were performed in IgorPro and in MATLAB.
Scoring data in clustering analyses were collected manually by examining individual planes of 1 μM z-stacks. 54 antennal glomeruli, as depicted in Fig. 3b, were identified in the nc82 or N-cadherin channel. A glomerulus was considered innervated if any green-colored LN process entered the nc82/N-cadherin dense region of a glomerulus.
Hierarchical clustering of innervation data was performed using freely available Cluster (http://bonsai.ims.u-tokyo.ac.jp/~mdehoon/software/cluster/software.htm) or R (http://www.r-project.org/) software. Data were clustered using complete linkage and a Euclidean distance measure in either program. Results were visualized as dendrograms and heatmaps using JavaTreeView (http://jtreeview.sourceforge.net/) or native R graphics facilities. The coherence of clusters was assayed in R using bootstrap analysis available through the pvclust package with n≥5000 at p-value>95% and at p-value>90%. Hierarchical clustering was verified by comparison to iterative k-means based clustering. k was determined for each iteration by plotting the within groups sum of squares or by using the GAP statistic. Clusters identified by hierarchical clustering that correlate with those revealed by iterative k-means based clustering were judged significant.
Innervation data were also imported into Matlab R2007a (Mathworks), where principal component analysis (PCA) was performed using customized software. Data were centered by subtracting off column means but were not rescaled.
Theoretical probabilities for LN6 innervation were calculated based on the assumption that all glomeruli have an equal chance of being innervated. Thus, for a given set of n innervation patterns, the chance a given glomerulus is innervated = (the total number of innervated glomeruli in the dataset/54)/n. Assuming that the collection of innervation probabilities is binomially distributed, the variation in observations can be calculated as the standard error = √[p*(1-p)/n], where p = probability of innervation and n = number of samples. These theoretical probabilities are presented in Fig. 6f as the average % innervation by glomerulus ± standard error as the upper and lower bounds.
We use the following intuition to quantify the overlap between patchy and non-patchy sister cells. For non-patchy sister cells colored red and green, a given volume of the antennal lobe of small size should contain both red and green signal, while for patchy sister cells, the same volume of antennal lobe could conceivably contain just red or just green signal. Only if the volume of interest were much larger should the strength of the overall red and green signals within the volume of interest become more comparable. To test this idea, cell bodies and brain regions outside the antennal lobe were first excluded from the analysis. An intensity threshold was then applied to sister cell image stacks to isolate the innervation patterns from background with all pixels below the 97.5th percentile set to 0. Such thresholded images were then dilated using cube shaped kernels of varying size. For image vectors and , an index of overlap OI was computed as follows:
using the Euclidean norm and where θ is the angle between and . Finally, we plotted the overlap index as a function of dilation kernel size.
To eliminate both maxillary palps in early development, smo3 mutant clones were induced by eyFLP in flies expressing Line5 or GH146-Gal4 driven mCD8GFP. Adults missing both MPs but retaining both antennae were dissected and their brains were analyzed for LN or PN innervation in glomeruli targeted by MP ORNs.
To remove ORN input in adults, MARCM ready Line6 flies were heat-shocked between 24-72h after egg laying to generate single cell LN clones. Antennae, maxillary palps or both were severed and removed between 0-2 days after eclosion using a #5 forceps. Flies were aged between 10-20 days before dissection and immunohistochemistry. Confocal images were loaded into Imaris where individual glomeruli were traced to obtain volumetric measurements. LN projection length within individual glomeruli as well as Syt-HA puncta were also quantified using Imaris.
We thank Ulrike Heberlein and Lisa Marin respectively for providing and screening unpublished Gal4 lines which lead to identification of Line1 and Line6; Julie Simpson for providing unpublished LCCH3 (Line7) Gal4; the Bloomington Stock Center, Kyoto Stock Center, NIG (GETDB), and Developmental Studies Hybridoma Bank for other reagents. We are grateful for the help of Jess Brooks in data collection and the Stanford Department of Statistics Consulting Service for technical help with statistical analyses. We thank Tom Clandinin, Greg Jefferis, and members of the Luo and Wilson labs for helpful comments on the manuscript. This work was supported by NIH grants to L.L. (R01-DC005982) and R.I.W. (R01-DC008174), a Pew Scholar Award, a McKnight Scholar Award, a Sloan Foundation Research Fellowship, and Beckman Young Investigator Award (to R.I.W). M.L.S. is an NRSA predoctoral fellow. E.Y. is supported by a Human Frontiers Science Program Long Term Fellowship. J.C.S.L. is supported by the Medical Scientist Training Program at Stanford University. L.L. is an investigator of the Howard Hughes Medical Institute.
AUTHOR CONTRIBUTIONS: Y. C. and M.L.S. performed the anatomical and developmental experiments. E.Y. performed the physiological experiments. L.C.S.L. helped with statistical analysis. L.L. and R.I.W. supervised the project and wrote the paper.