Olfactory responses of Drosophila are encoded in the organization of projection neurons

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Appendix 1
Appendix 2
Data availability
References
Article and author information
Metrics

Abstract

The projection neurons (PNs), reconstructed from electron microscope (EM) images of the Drosophila olfactory system, offer a detailed view of neuronal anatomy, providing glimpses into information flow in the brain. About 150 uPNs constituting 58 glomeruli in the antennal lobe (AL) are bundled together in the axonal extension, routing the olfactory signal received at AL to mushroom body (MB) calyx and lateral horn (LH). Here we quantify the neuronal organization in terms of the inter-PN distances and examine its relationship with the odor types sensed by Drosophila. The homotypic uPNs that constitute glomeruli are tightly bundled and stereotyped in position throughout the neuropils, even though the glomerular PN organization in AL is no longer sustained in the higher brain center. Instead, odor-type dependent clusters consisting of multiple homotypes innervate the MB calyx and LH. Pheromone-encoding and hygro/thermo-sensing homotypes are spatially segregated in MB calyx, whereas two distinct clusters of food-related homotypes are found in LH in addition to the segregation of pheromone-encoding and hygro/thermo-sensing homotypes. We find that there are statistically significant associations between the spatial organization among a group of homotypic uPNs and certain stereotyped olfactory responses. Additionally, the signals from some of the tightly bundled homotypes converge to a specific group of lateral horn neurons (LHNs), which indicates that homotype (or odor type) specific integration of signals occurs at the synaptic interface between PNs and LHNs. Our findings suggest that before neural computation in the inner brain, some of the olfactory information are already encoded in the spatial organization of uPNs, illuminating that a certain degree of labeled-line strategy is at work in the Drosophila olfactory system.

Editor's evaluation

Choi et al. explore how olfactory information flows across the three major neuropils in the Drosophila brain – the antennal lobe (AL), mushroom Body (MB), and the lateral horn (LH). They use the two connectomes of adult Drosophila and 'inter-PN distances' to do this. Using this neuroanatomy based approach, they find support for a labeled-line strategy, which they subsequently test for with synaptic connectivity data for a subset of PNs. They find that while some labelled lines may exist, PNs generally participate in multi-channel integration at the MB and LH. This manuscript will be of interest to neuroscientists interested in olfactory processing and to those working on connectomic-level circuit analysis.

https://doi.org/10.7554/eLife.77748.sa0

Introduction

Anatomical details of neurons obtained based on a full connectome of the Drosophila hemisphere reconstructed from electron microscope (EM) image datasets (Bates et al., 2020; Scheffer et al., 2020) offer the wiring diagram of the brain, shedding light on the origin of brain function. Out of the immense amount of data, we study the second-order neurons, known as the projection neurons (PNs) of the olfactory system. It is the PNs that bridge the olfactory receptor neurons (ORNs) in the antenna and maxillary palp to higher olfactory centers where neural computation occurs for Drosophila to sense and perceive the environment (Hallem and Carlson, 2004a). The three neuropils, namely the antennal lobe (AL), mushroom body (MB) calyx, and lateral horn (LH), are the regions that abound with an ensemble of axonal branches of PNs and synapses (Figure 1). PNs can be classified as uniglomerular and multiglomerular PNs based on their structure and connectivity to other PNs. The uniglomerular PNs (uPNs) in AL constitute glomeruli that collect olfactory signals from ORNs of the same receptor type (Gao et al., 2000; Couto et al., 2005). uPNs innervating MB calyx and LH relay the signals further inside the brain through synaptic junctions with the Kenyon cells (KCs) and lateral horn neurons (LHNs), respectively. Multiglomerular PNs (mPNs), on the other hand, innervate multiple glomeruli, often contributing to the inhibitory regulation of signals relayed from ORNs to third-order olfactory neurons (Berck et al., 2016). PNs can functionally be categorized into either excitatory (cholinergic) or inhibitory (GABAergic), where a many GABAergic PNs tend to bypass MB calyx while innervating multiple glomeruli in AL (and hence are mPNs) (Schultzhaus et al., 2017; Shimizu and Stopfer, 2017).

Figure 1

Download asset Open asset

A schematic of the *Drosophila* olfactory system.

uPNs comprising each glomerulus in AL collect input signals from ORNs of the same receptor type and relay the signals to MB calyx and LH. uPNs in MB calyx synapse onto KCs; and uPNs in LH synapse onto LHNs.

Since the seminal work by Cajal, 1911, who recognized neurons as the basic functional units of the nervous system, there have been a series of attempts at classifying neurons using different representations of neuronal morphologies and at associating the classified anatomies with their electrophysiological responses and functions (Uylings and van Pelt, 2002; Scorcioni et al., 2008; Jefferis et al., 2007; Seki et al., 2010; Gillette and Ascoli, 2015; Lu et al., 2015; Li et al., 2017; Kanari et al., 2018; Mihaljević et al., 2018; Gouwens et al., 2019; Laturnus et al., 2020). Systematic and principled analyses of neuronal anatomy would be a prerequisite for unveiling a notable link between the PN organization and olfactory representations. Several different metrics involving spatial projection patterns (Jefferis et al., 2007), electrophysiological properties (Seki et al., 2010; Gouwens et al., 2019), topological characteristics (e.g. morphometrics) (Uylings and van Pelt, 2002; Scorcioni et al., 2008; Lu et al., 2015; Mihaljević et al., 2018; Gouwens et al., 2019), intersection profiles (Gouwens et al., 2019), and NBLAST scores (Jeanne et al., 2018; Zheng et al., 2018; Bates et al., 2020; Scheffer et al., 2020) have been utilized in the past. More recently, machine learning approaches have been popularized as a tool for classification tasks (Vasques et al., 2016; Buccino et al., 2018; Mihaljević et al., 2018; Zhang et al., 2021).

Among a multitude of information that can be extracted from the neural anatomy associated with uPNs, the inter-PN organization draws our attention. To compare spatial characteristics of uPNs across each neuropil and classify them based on the odor coding information, we confine ourselves to uPNs innervating all three neuropils, most of which are cholinergic and follow the medial antennal lobe tract (mALT). Within this scope, we first calculate inter-PN distance matrices in each neuropil and study them in reference to the glomerular types (homotypes) to discuss how the inter-PN organization changes as the PNs extend from AL to MB calyx and from AL to LH.

In this study, we utilize two representative EM-based reconstruction datasets for the analysis (the latest FAFB Bates et al., 2020 and the hemibrain datasets Scheffer et al., 2020). The FAFB dataset specifically encompasses the Drosophila olfactory system, while the hemibrain dataset aims for a reconstruction of the entire right hemisphere of the Drosophila brain. The results based on the two datasets are largely consistent and interchangeable, which generalizes our findings.

We have conducted statistical analyses to unravel potential associations between the uPN organization and the behavioral responses of Drosophila to external stimuli encoded by glomerular homotypes, finding that certain odor types and behavioral responses are linked to a characteristic inter-neuronal organization. The map of synaptic connectivity between uPNs and the third-order neurons (KCs and LHNs in MB calyx and LH, respectively) complements the functional implication of the association between the inter-PN organization and olfactory processing. A ‘labeled-line design’ in olfaction is generally considered to exhibit a chain of neurons dedicated to encoding a single olfactory feature with no direct integration with other features as the signal is passed onto higher-order neurons. While we do not demonstrate the full architecture of labeled-line design in the Drosophila olfactory system as the signals from odor-sensing by ORNs are passed down to the inner brain for perception, our analysis shows that homotypic uPNs encoding particular odor types not only maintain their spatially localized and bundled structure throughout all three neuropils but also display synaptic connections that converge to a narrow subset of third-order neurons. The Drosophila olfactory system leverages the efficiency of the labeled-line design in sensory information processing (Min et al., 2013; Howard and Gottfried, 2014; Andersson et al., 2015; Galizia, 2014).

Results

Spatial organization of neurons inside neuropils

The inter-PN distance $d_{α β}$

First, we define a metric with which to quantify the spatial proximity between neurons. Specifically, the inter-PN distance $d_{α β}$ represents the average taken over the minimum Euclidean distances between two uPNs $α$ and $β$ , such that $d_{α β}$ is small when two uPNs are tightly bundled together (see Equation 1 and Figure 2—figure supplement 1A). Although metrics such as the NBLAST score (Costa et al., 2016) and others (Kohl et al., 2013) can be used to study the PN organization, these metrics take both the morphological similarity and the spatial proximity into account. The distance $d_{α β}$ only measures the pairwise distance but not the dot product term (which measures the similarity of two neuronal morphologies), whereas the NBLAST score considers both. Therefore, while the distance $d_{α β}$ is computationally comparable to the NBLAST score, it only measures the spatial proximity between two neurons. We notice that the features of the uPN organization captured by the NBLAST distance are not necessarily aligned with $d_{α β}$ (see Figure 2—figure supplement 1B). The two distances are correlated but with significant dispersion, indicating that these two metrics are not the same. Since we are solely interested in the spatial proximity (or co-location) between two uPN innervations but not the morphological similarity between them (which the NBLAST score accounts for, a point also noted by Zheng et al., 2018), we deliberately chose the metric $d_{α β}$ instead of the NBLAST score for our analyses.

The distances $d_{α β}$ (Equation 1) between all the possible pairs ( $α$ and $β$ ) of 135 uPNs from the FAFB dataset are visualized in the form of a matrix (Figure 2). We perform hierarchical clustering on the distance matrix for uPNs in each neuropil (see the outcomes of $d_{α β}$ -based clustering analysis in Figure 2—figure supplement 2 and Materials and methods for the details). Individual clusters from the hierarchical clustering of uPNs in MB calyx and LH from the FAFB dataset are visualized in Figure 3 and Figure 4 with the colors denoting the odor types encoded by the individual uPNs, which will be discussed in detail later.

Figure 2 with 3 supplements see all

Download asset Open asset

The three matrices representing the pairwise distances $d_{α β}$ in units of $μ m$ between individual uPN in AL, MB calyx, and LH.

The matrices are calculated based on uPNs available in the FAFB dataset. The diagonal blocks reflect the homotypic uPNs comprising the 57 glomerular homotypes defined in the FAFB dataset (Bates et al., 2020), labeled at the edges.

Figure 3

Download asset Open asset

The $d_{α β}$ -based clustering on uPNs based on the FAFB dataset in MB calyx resulting in 10 clusters.

The individual uPNs are color-coded based on the encoded odor types (Dark green: decaying fruit, lime: yeasty, green: fruity, gray: unknown/mixed, cyan: alcoholic fermentation, red: general bad/unclassified aversive, beige: plant matter, brown: animal matter, purple: pheromones, pink: hygro/thermo) (Mansourian and Stensmyr, 2015; Bates et al., 2020). The first and second columns illustrate the dorsal and the anterior view, respectively (D: dorsal, M: medial, P: posterior). The black line denotes the approximate boundary of MB calyx.

Figure 4 with 2 supplements see all

Download asset Open asset

The $d_{α β}$ -based clustering on uPNs based on the FAFB dataset in LH resulting in 11 clusters.

(inset) A cartoon illustrating the relative position between clusters $C_{I}^{LH} = C_{3}^{LH}$ and $C_{O}^{LH} = C_{4, 5, 6, 9}^{LH}$ . The individual uPNs are color-coded based on the encoded odor types (Dark green: decaying fruit, lime: yeasty, green: fruity, gray: unknown/mixed, cyan: alcoholic fermentation, red: general bad/unclassified aversive, beige: plant matter, brown: animal matter, purple: pheromones, pink: hygro/thermo). The first and second columns illustrate the dorsal and the anterior view, respectively (D: dorsal, M: medial, P: posterior). The black line denotes the approximate boundary of LH.

Spatial proximity-based clustering results

In MB calyx, the hierarchical clustering divides the uPNs from the FAFB dataset into 10 clusters (Figure 3). Clusters $C_{2}^{MB}$ and $C_{10}^{MB}$ largely encompass the dorsal region and clusters $C_{6}^{MB}$ and $C_{7}^{MB}$ encompass the ventral region of the neuropil. The cluster $C_{7}^{MB}$ shows a characteristic biforked pattern projecting to the lateral and medial regions. The cluster $C_{3}^{MB}$ also exhibits the same structural pattern but is composed of a tight bundle of uPNs that are part of DL2d and DL2v (both of which encodes food-related odors). The cluster $C_{8}^{MB}$ is located between the biforked innervation pattern of clusters $C_{6}^{MB}$ and $C_{7}^{MB}$ , and predominantly innervates the posterior region. Lastly, clusters $C_{1}^{MB}$ , $C_{4}^{MB}$ , and $C_{5}^{MB}$ , innervate the anterior region of MB calyx, spatially separated from other uPNs.

In LH, 11 clusters are identified in the FAFB dataset (Figure 4). The cluster $C_{3}^{LH}$ is the largest, which mainly innervates the dorsal posterior region of LH. Clusters $C_{4}^{LH}$ , $C_{5}^{LH}$ , $C_{6}^{LH}$ , and $C_{9}^{LH}$ display variable biforked projection patterns along the coronal plane, enveloping the boundary of the cluster $C_{3}^{LH}$ . This creates a spatial pattern where a large blob of uPNs ( $C_{I}^{LH}$ ) are surrounded by a claw-like structure ( $C_{O}^{LH}$ ) (Figure 4, inset). Clusters $C_{1}^{LH}$ , $C_{2}^{LH}$ , and $C_{7}^{LH}$ innervate the anterior-ventral region and display clear segregation from the other uPNs. Another group composed of clusters $C_{10}^{LH}$ and $C_{11}^{LH}$ innervates the posterior-ventral-medial region.

We use Pearson’s $χ^{2}$ -test (see Materials and methods for the details) to assess the likelihood of dependence between the $d_{α β}$ -based clustering outputs for MB calyx, LH, and the glomerular labels (homotypes) statistically significant correlations are found in terms of both the p-value and the Cramér’s $V$ (see Appendix 1—table 1 and Methods for a detailed explanation of the meaning behind the p-value and the Cramér’s $V$ ), the latter of which is analogous to the correlation coefficient for the $χ^{2}$ -test. The mutual information between the same set of nominal variables, which is calculated to verify our $χ^{2}$ -tests (see Materials and methods), offers a similar conclusion (see Appendix 1 and Appendix 1—table 3).

We also categorize the spatial organization of uPNs in reference to the glomerular labels. The homotypic uPNs constituting a tightly bundled glomerulus in AL manifest themselves as the block diagonal squares in the $d_{α β}$ -matrix (Figure 2). This is apparent in the dendrogram constructed from the distance matrix for the uPNs at AL (Figure 2—figure supplement 3), where uPNs sharing the same glomerular label are grouped under a common branch, thereby demonstrating the spatial proximity between uPNs forming the same glomerulus. The $d_{α β}$ -matrix indicates that such organizations are also preserved in MB calyx and LH. However, clear differences are found in the off-diagonal part of $d_{α β}$ matrices (Figure 2).

The same hierarchical clustering analysis performed on the hemibrain dataset results in 14 clusters for uPN innervation in MB calyx and 13 clusters in LH. Despite the differences in the number of clusters, we find that spatial and structural characteristics of individual clusters observed from the FAFB dataset are well translated and comparable to those from the hemibrain dataset (see the clustering result in Figure 8—figure supplement 1). Furthermore, various statistical tests used in this study (e.g. Pearson’s $χ^{2}$ -test) on the hemibrain dataset lead to the same conclusion (see Appendix 1, Appendix 1—table 2, and Appendix 1—table 4).

The degrees of bundling, packing, and overlapping

To conduct a quantitative and concise analysis of $d_{α β}$ matrices, we define the mean intra- and inter-homotypic uPN distances, ${\bar{d}}_{intra, X}$ and ${\bar{d}}_{inter, X}$ (see Methods for detailed formulation). The ${\bar{d}}_{intra, X}$ is the average distance between uPNs in the same homotype and measures the degree of uPNs in the homotype $X$ being bundled. Therefore, a smaller ${\bar{d}}_{intra, X}$ signifies a tightly bundled structure of $X$ -th homotypic uPNs (see Figure 5—figure supplement 1 for raw ${\bar{d}}_{intra, X}$ values). Similarly, ${\bar{d}}_{inter, X}$ , which measures the degree of packing (or segregation), is defined as the average distance between the neurons comprising the $X$ -th homotype and neurons comprising other homotypes. Thus, a small value of ${\bar{d}}_{inter, X}$ signifies tight packing of heterotypic uPNs around $X$ -th homotype, while a large value indicates that the homotypic uPNs comprising the homotype $X$ are well segregated from other homotypes (see Figure 5—figure supplement 1 for raw ${\bar{d}}_{inter, X}$ values).

The degrees of bundling averaged over all homotypes ( ${\bar{d}}_{i n t r a} = N_{X}^{- 1} \sum_{X}^{N_{X}} {\bar{d}}_{i n t r a, X} \approx 4 μ m$ ) are comparable over all three neuropils (blue dots in Figure 5A and B). On the other hand, from ${\bar{d}}_{inter}$ , which is defined as the mean inter-homotype distance averaged over all $X$ s, we find that homotypic uPNs are well segregated from others in AL as expected, whereas spatial segregation among homotypes is only weakly present in MB calyx (orange dots in Figure 5A and B and the cartoon of Figure 5C). We also observe that the ${\bar{d}}_{intra}$ and ${\bar{d}}_{inter}$ are comparable for the two different datasets. A minor difference is observed in ${\bar{d}}_{intra}$ , indicating a slightly tighter bundling structure for the hemibrain dataset.

Figure 5 with 2 supplements see all

Download asset Open asset

Organization of homotypic uPNs in the three neuropils.

Plots of ${\bar{d}}_{intra}$ (blue, degree of bundling), ${\bar{d}}_{inter}$ (orange, degree of packing), and the ratio between the two distances $λ$ (red, degree of overlapping) calculated based on (A) the FAFB dataset and (B) the hemibrain dataset. Error bars depict the standard deviation. (C) Diagram illustrating the overall organization of uPNs at each neuropil. Homotypic uPNs are tightly bundled and segregated in AL. Several groups of homotypic uPNs form distinct heterotypic spatial clusters at higher olfactory centers, extensively overlapping in MB calyx (see Figure 3).

Next, we take the ratio of mean intra- to inter-PN distances of $X$ -th homotype as $λ_{X}$ to quantify the degree of overlapping around $X$ -th homotype (see Materials and methods). The term ‘overlapping’ is specifically chosen to describe the situation where different homotypes are occupying the same space. A large value of $λ_{X}$ (particularly $λ_{X} > 0.4$ ) suggests that the space occupied by the uPNs of the $X$ -th homotype is shared with the uPNs belonging to other homotypes. The value $λ (= N_{X}^{- 1} \sum_{X}^{N_{X}} λ_{X})$ averaged over all the homotypes (red in Figure 5A and B) suggests that the extent of overlapping between uPNs is maximal in MB calyx and minimal in AL (Figure 5C).

Figure 6A and B, Figure 7, and Figure 7—figure supplement 1 show individual values of $λ_{X}$ for all homotypes in the three neuropils. We identify the following features: (i) In AL, $λ_{X} \leq 0.4$ for all homotypes except DL5 (a homotype encoding aversive odors), indicating that homotypic uPNs are tightly bundled and segregated from uPNs in other glomeruli. The same trend is observable in the hemibrain dataset (Figure 6B), but with $λ_{D L 5} \leq 0.4$ .; (ii) In MB calyx, a large portion ( $\approx 65 %$ ) of $λ_{X}$ ’s exceed 0.4 and even the cases with $λ_{X} > 1$ are found (VC5, DL5), implying that there is a substantial amount of overlap between different homotypes. In the hemibrain dataset, $\approx 76 %$ of $λ_{X}$ ’s exceed 0.4.; (iii) Although not as significant as those in AL, many of uPNs projecting to LH are again bundled and segregated in comparison to those in MB calyx (see Figure 7B). (iv) The scatter plot of $λ_{X}$ between MB calyx and LH (Figure 7C) indicates that there exists a moderate positive correlation ( $r = 0.642, p < 0.0001$ ) between $λ_{X}$ at MB calyx and LH. This implies that a higher degree of overlapping in MB calyx carries over to the uPN organization in LH. The result from the hemibrain dataset is similar ( $r = 0.677$ , $p < 0.0001$ , see Figure 7—figure supplement 1).

Figure 6 with 1 supplement see all

Download asset Open asset

Degree of overlapping between inter-homotypic uPNs, $λ_{X}$ ( $X = VM4, VC5, \dots, VP5$ ).

The degree of overlapping ( $λ_{X}$ ) for $X$ -th homotype in AL, MB calyx, and LH (from lighter to darker colors) calculated from the uPNs based on (A) the FAFB dataset and (B) the hemibrain dataset. The homotype label is color-coded based on the odor types associated with the glomerulus obtained from the literature and is sorted based on the value of $λ_{X}$ for each odor type at LH. Asterisks (*) mark homotypes composed of a single uPN while plus (+) mark homotypes composed of a single uPN under our selection criterion but are actually a multi-uPN homotype, whose intra-homotype uPN distance is not available. O (attractive) and X (aversive) indicate the putative valence information collected from the literature. The blue horizontal line denotes $λ_{X} = 0.4$ . (C) Two homotypes taken from the FAFB dataset, DL3 (purple) and DL5 (red), which are indicated by yellow triangles in (A), are highlighted along with other uPNs (gray).

Figure 7 with 1 supplement see all

Download asset Open asset

Scatter plots depicting the relationships between $λ_{X}$ s at two different neuropils calculated from the uPNs based on the FAFB dataset; (A) AL versus MB calyx, (B) AL versus LH, and (C) LH versus MB calyx.

The color code is the same as in Figure 6. The blue lines in (A) and (B) denote $λ_{X} = 0.4$ .

The entire neuron morphologies of uPNs from two homotypes with a small ( $X =$ DL3, which largely responds to pheromones) and a large ( $X =$ DL5) $λ_{X}$ s in LH are visualized along with the other uPNs (gray) (Figure 6C). The homotype DL3, which seldom overlaps with others in AL ( $λ_{DL3} \approx 0.07$ ) and LH ( $λ_{DL3} \approx 0.17$ ), displays an increased overlapping in MB calyx ( $λ_{DL3} \approx 0.31$ ). Therefore, DL3 is tightly packed in AL and LH, whereas it is relatively dispersed in MB calyx. Meanwhile, the homotype DL5 displays a significant dispersion in all three neuropils, although the dispersion is the smallest in AL ( $λ_{DL5} \approx 0.74$ ) compared to that in MB calyx ( $λ_{DL5} \approx 1.1$ ) and LH ( $λ_{DL5} \approx 1.5$ ).

There are minor variations between the FAFB and the hemibrain datasets in terms of ${\bar{d}}_{intra, X}$ , ${\bar{d}}_{inter, X}$ , and $λ_{X}$ , and they likely arise from the factors such as a minor mismatch in the glomerulus label annotations that sometimes affects the number of uPNs constituting a given homotype, and the difference in the number of uPNs between two datasets as a result of our selection criterion. Regardless, still present in both datasets are the spatial and organizational trends described above. Taken together, the organization of olfactory uPNs varies greatly in the three neuropils. The clear homotype-to-homotype segregation in AL no longer holds in MB calyx. Instead, the $d_{α β}$ -based clustering suggests the presence of clusters made of multiple different homotypic uPNs (Figure 5C). For some homotypes, the well-segregated organizations in AL are recovered when they reach LH (compare Figure 7A and Figure 7B).

Relationship between neuronal organization and olfactory features

Now we explore how the structural features identified from our clustering outputs are associated with odor types and valences (behavioral responses). As briefly mentioned earlier, the color codes in Figure 3, Figure 4, Figure 6, and Figure 7 depict odor types encoded by corresponding homotypic uPNs, which follow the same categorical convention used by Mansourian and Stensmyr, 2015 and Bates et al., 2020. The O and X represent the putative valence, which indicates whether Drosophila is attracted to or repelled from the activation of specific homotypic uPNs. For example, DA2 responds to geosmin, a chemical generated from harmful bacteria and mold, which evokes a strong repulsion in Drosophila (Stensmyr et al., 2012). Similarly, VM3 is suggested to encode repulsive odors, while VM2 and VM7d encode attractive odors (Mansourian and Stensmyr, 2015; Bates et al., 2020). Overall, the following information is acquired from the literature (Hallem et al., 2004b; Galizia and Sachse, 2010; Mansourian and Stensmyr, 2015; Badel et al., 2016; Bates et al., 2020) and labeled accordingly:

DA1, DA3, DL3, DM1, DM4, VA1v, VA2, VA3, VC1, VC2, VM1, VM2, VM4, VM5d, VM5v, VM7d, and VM7v (17 homotypes) encode attractive (O) odor.
D, DA2, DA4l, DA4m, DC1, DC2, DC3, DC4, DL1, DL4, DL5, DM2, DM3, DM5, DM6, DP1m, V, VA5, VA6, VA7l, VA7m, VC3, VL2a, VL2p, and VM3 (25 homotypes) encode aversive (X) odor.
The remaining homotypes are characterized as either unknown, non-preferential, or conflicting valence information.

Collecting the glomerular types of tightly bundled homotypic uPNs with $λ_{X} < 0.4$ in LH (Figure 6, Figure 7, and Figure 7—figure supplement 1), we explore the presence of any organizational trend.

In LH, out of 37 homotypes composed of multiple uPNs ( $2 \leq n \leq 8$ ) based on our selection criterion, 29 homotypes (DL2v, DL2d, VM1, VL1, DM6, VM7d, VA3, VM5v, DA3, VM2, DL1, VA7m, VC3, VM7v, VC4, V, DM2, VM3, DA2, D, DC2, VA5, VA1v, DA1, DC3, DL3, VA1d, VP1d, and VP1l) satisfy the condition of $λ_{X} < 0.4$ . In the hemibrain dataset, a couple of homotypes (VM1, DL5, DL3, and VP1l) are suggested to be single uPN homotypes based on our selection criterion.
Homotypes VA1v, DA1, DC3, DL3, and VA1d (colored purple in Figure 3, Figure 4, Figure 6, and Figure 7) encode pheromones involved with reproduction (Grabe et al., 2016; Bates et al., 2020; Dweck et al., 2015), and VM4, VM1, VM7d, DM1, DM4, VC2, VM5d, VA3, VM5v, DA3, and VM2 encode odors presumed to be associated with identifying attractive food sources (Couto et al., 2005; Semmelhack and Wang, 2009; Mohamed et al., 2019; Bates et al., 2020) (see Figure 6A). A previous work (Grosjean et al., 2011) has identified a group of glomeruli that co-process food stimuli and pheromones via olfactory receptor gene knock-in coupled with behavioral studies. The list of homotypes mentioned above is largely consistent with those glomeruli reported by Grosjean et al., 2011.
Homotypes DM6, DM2, VM3, VL2p, DA2, and D are likely associated with aversive food odors. DA2 responds to bacterial growth/spoilage; VL2p, DM2, and VM3 to the alcoholic fermentation process; DM6 and D to flowers (Galizia and Sachse, 2010; Bates et al., 2020).
Many homotypes responding to odors which can be described as kairomones, a type of odors emitted by other organisms (Kohl et al., 2015), are part of the 29 homotypes with $λ_{X} < 0.4$ . This includes the pheromone encoding groups (VA1v, DA1, DC3, DL3, and VA1d) and others such as DA2, VC3, and VA5, which respond to geosmin, 1-hexanol, and 2-methyl phenol, respectively (Hallem et al., 2004b; Galizia and Sachse, 2010).

Figure 8 recapitulates the cluster information from $d_{α β}$ -based analysis along with homotypes, odor types (color-codes), and putative valence (attractive (O) and aversive (X) odors). Some points are worth making:

Even though uPNs innervating MB calyx exhibit large $λ_{X}$ s, the hierarchical clustering grouped homotypic uPNs together. This suggests the homotypic uPNs are still proximal in MB calyx, indicating the reduction in $d_{inter}$ is what is driving the increase in overlapping. This is already shown through ${\bar{d}}_{intra}$ in Figure 5A, B and is supported by our statistical tests (see Appendix 1—table 1 and Appendix 1—table 3). The same is true for LH. The grouping of homotypic uPNs is also observable from the hemibrain dataset (Figure 8—figure supplement 1).
In the FAFB dataset, 13 out of 57 homotypes are made of a single uPN ( $n = 1$ , the asterisked glomeruli in Figure 6A and Figure 8), which tend to be characterized by comparatively dense branched structures (see Figure 6—figure supplement 1), suggestive of homotypic uPN number dependence for the neuron morphology. Among the 13 homotypes, 7 encode aversive stimuli (X), 4 encode attractive stimuli (O), and 2 have no known valence information (see Appendix 1—table 5). In the hemibrain dataset, 7 encode aversive stimuli (X), 5 encode attractive stimuli (O), and 1 has no known valence information (see Appendix 1—table 6). The relative prevalence of single-uPN homotypes encoding aversive stimuli is noteworthy.
In LH, the cluster $C_{1}^{LH}$ , located in the anterior-ventral region of the neuropil, is composed only of pheromone-encoding homotypic uPNs, DA1 and DC3. The cluster $C_{2}^{LH}$ is also mostly composed of pheromone-encoding homotypic uPNs, DL3 and VA1d (Figure 4 and Figure 8), which is consistent with the results by Jefferis et al., 2007. In MB calyx, the majority of the uPNs encoding pheromones, except DL3, are grouped into the cluster $C_{8}^{MB}$ (see Figure 3 and Figure 8). A similar trend is observed in the hemibrain dataset, although the arbitrary cluster labels differ (see clusters $C_{4}^{LH}$ , $C_{8}^{LH}$ , and $C_{10}^{MB}$ in Figure 8—figure supplement 1).
Hygro/thermo-sensing homotypes such as VP2 and VP4 are spatially segregated from other odor-encoding uPNs, which is observable through clusters composed predominantly of hygro/thermo-sensing homotypes (see Figure 8 and Figure 8—figure supplement 1). In MB calyx, these neurons rarely project to anterior region and are distributed along the base of the neuropil. This is in line with previous literature (Li et al., 2020). In LH, they are clustered in the posterior-ventral-medial region, hardly innervating the neuropil but covering the medial side of the neuropil (Figure 3 and Figure 4).
Along with the clusters of uPNs visualized in Figure 3 and Figure 4, of particular note are the clusters formed by a combination of several homotypic uPNs. A large portion of uPNs innervating LH that encodes potentially aversive responses are grouped into clusters $C_{4}^{LH}$ , $C_{5}^{LH}$ , $C_{6}^{LH}$ , and $C_{9}^{LH}$ , which envelop the cluster $C_{3}^{LH}$ where mostly food-related homotypes converge (Figure 4). In the hemibrain dataset, these correspond to $C_{10}^{LH}$ and $C_{11}^{LH}$ for the aversive responses and $C_{6}^{LH}$ and $C_{13}^{LH}$ for the food-related homotypes (Figure 8—figure supplement 1).

Figure 8 with 1 supplement see all

Download asset Open asset

A diagram summarizing how the clusters of uPNs in MB calyx (10 clusters) and LH (11 clusters) from the FAFB dataset are associated with the odor types (Dark green: decaying fruit, olive: yeasty, green: fruity, cyan: alcoholic fermentation, red: general bad/unclassified aversive, beige: plant matter, brown: animal matter, purple: pheromones, gray: unknown, pink: hygro/thermo).

Asterisks (*) mark homotypes composed of a single uPN while plus (+) mark homotypes composed of a single uPN under our selection criterion but are actually a multi-uPN homotype, whose intra-homotype uPN distance is not available. O and X represent the putative valence information collected from the literature (O: attractive, X: aversive).

Given that the synaptic communications with KCs and LHNs are critical for neural computation in the inner brain, the specific type of uPN organization in each neuropil should be of great relevance. Indeed, it has been suggested that the spatial convergence, segregation, and overlapping of different homotypic uPNs within neuropil influence the information processing in higher olfactory centers (Grosjean et al., 2011).

According to previous studies (Jefferis et al., 2007; Liang et al., 2013; Kohl et al., 2013; Fişek and Wilson, 2014), uPN innervation in LH and LHNs are highly stereotyped in terms of connectivity and response. Homotypic uPNs are spatially organized in AL, and to a certain degree, in LH, based on the odor type and valence information (Min et al., 2013; Huoviala et al., 2020). The presence of tightly bundled anatomy of homotypic uPNs ( $λ_{X} < 0.4$ ) in both AL and LH (Figure 7B and Figure 7—figure supplement 1B) may imply that the Drosophila olfactory system dedicates a part of the second-order neural circuit on behalf of the ‘labeled-line’ design, which enables the organism to sense urgent chemical stimuli at the early stage of information processing without going through more sophisticated neural computation in the inner brain (Howard and Gottfried, 2014; Andersson et al., 2015; Min et al., 2013).

Labeled-line design of the higher order olfactory neurons

The concept of labeled-line design is widely considered at work at the ORN-PN interface (AL) as the signal generated from specific olfactory receptors converges to a single glomerulus (Vosshall et al., 2000; Couto et al., 2005; Fishilevich and Vosshall, 2005). A potential labeled-line strategy or separated olfactory processing of aversive odors encoded by DA2 has been extensively discussed (Stensmyr et al., 2012; Seki et al., 2017; Huoviala et al., 2020). It has been shown that pheromone-encoding homotypes in LH (Jefferis et al., 2007; Ruta et al., 2010; Kohl et al., 2013; Frechter et al., 2019; Bates et al., 2020; Das Chakraborty and Sachse, 2021) are at work in specific third-order olfactory neurons. So far, we have shown that the labeled-line design is present in the architecture of higher olfactory centers of second-order neurons, that is, MB calyx and LH, where homotypic uPNs are tightly bundled together despite the lack of glomerular structure. In this section, we will conduct a comprehensive analysis of the synaptic connectivity between PNs and third-order olfactory neurons (KCs and LHNs) using three demonstrations. We ask (i) whether the labeled-line strategy implied in the uPN organization is translated over to the third-order olfactory neurons, (ii) to what extent the signals encoded by different homotypic uPNs are integrated at synaptic interfaces with the third-order neurons, and (iii) whether the spatial properties of pre-synaptic neurons (PNs) play any role in signal integration by the third-order neurons.

Homotype-specific connections

For the analysis of the interface between homotypic uPNs and third-order neurons, we study the connectivity matrices $C^{P N - K C}$ and $C^{P N - L H N}$ (see Figure 9, Figure 9—figure supplement 1), which are extracted from the hemibrain dataset (Scheffer et al., 2020). The $C^{ξ}$ ( $ξ =$ PN-KC or PN-LHN) is a binary matrix ( $C_{X, i}^{ξ} = 0$ or 1 dictating the connectivity) of synaptic connectivity between $X$ -th homotypic uPNs and $i$ -th third-order neuron (KC or LHN). It is observed that most of the KCs and LHNs integrate information from multiple homotypes, but that there are also a small number of KCs and LHNs that synapse only with a single homotype (Figure 10).

Figure 9 with 1 supplement see all

Download asset Open asset

A schematic illustrating the connectivity between homotypes ( $X = A, B, \dots, E$ ) and third-order neurons ( $i = 1, 2, \dots, 7$ ).

(A) The connectivity matrix $C$ , where $C_{X, i} = 1$ when any uPNs in the $X$ -th homotype and $i$ -th third-order neuron synapses and $C_{X, i} = 0$ otherwise. (B) The number of $X$ -th homotype-specific connections ( $N_{X, sp}$ ) and the total number of third-order neurons synapsed to any uPNs in the $X$ -th homotype. (C) The common synapse matrix ( $S$ ) whose element specifies the number of third-order neurons commonly connected between two homotypes. The homotype A is connected to three third-order neurons 1, 2, and 3 ( $N_{A, tot} = 3$ ). Neuron 1 is not synapsing with any other homotype but A, and hence $N_{A, sp} = 1$ ; similarly, $N_{D, sp} = 2$ (the blue lines depict specific connections). The signals from the two homotypes B and C are shared by the third-order neurons 2, 3, and 5; therefore, $S_{B C} = 3$ in the common synapse matrix $S$ .

Figure 10 with 1 supplement see all

Download asset Open asset

Bar graphs depicting the number of KCs/LHNs that synapse with a specific homotype $X$ ( $N_{X, sp}$ , blue) and the percentage of KCs/LHNs that synapse with a specific homotype $X$ ( $f_{X} = N_{X, sp} / N_{X, tot}$ , red) at (A) PN-KC and (B) PN-LHN interfaces, with the synaptic weight of $N = 3$ used as the threshold.

The ‘homotype-specific’ connections, defined as the number of third-order neurons that only synapses with a specific homotype but not with the others (see Figure 10 and Methods for more information) can be quantified in terms of the total number of third-order neurons in contact with $X$ -th homotypic uPNs, and it can be obtained by counting the non-zero elements of the matrix $C$ with fixed $X$ . For the case of the PN-KC interface, this number can be obtained from $N_{X, t o t}^{P N - K C} = \sum_{i = 1}^{1754} C_{X, i}^{P N - K C}$ . Specifically, Figure 10A shows $N_{X, sp}^{PN - KC}$ and those normalized by $N_{X, t o t}^{P N - K C}$ ( $f_{X} = N_{X, s p}^{P N - K C} / N_{X, t o t}^{P N - K C}$ , see Materials and methods for the detailed algorithms behind the calculation), for all homotypes ( $X = V M 4, V C 5, \dots, V P 4$ ). Compared to those in KCs, the ‘homotype-specific’ connections are much more prevalent in LHNs (Figure 10). Certain homotypic uPNs, in particular, the hygro/thermo-sensing homotypes are connected to the LHNs which are dedicated to process the signals from hygro/thermo-sensing homotypes ( $\geq 10 %$ of PN-LHN connections made by homotypes).

To address the concern with potential false positives in the detected synapses, we reexamine our results based on the synaptic connectivity with a higher threshold ( $N = 8$ ). Figure 10—figure supplement 1 demonstrates that the homotype-specific connections tend to increase under a more stringent synapse selection criterion, especially in LH. This is most notable in homotypes DM1, DM4, DP1l, and VM6. The existence of these ‘homotype-specific’ third-order neurons suggests that a subset of olfactory processing may rely on the labeled-line strategy that extends beyond the layer of second-order neurons to the higher brain center.

Third-order neurons mediate signal integration

Figure 11A, B show the ‘common synapse matrices’ representing the number of commonly connected third-order neurons between two homotypes $X$ and $Y$ ( $S_{X Y}^{η}$ with $η =$ KC or LHN), which provide glimpses into the extent of signal integration mediated by KCs and LHNs (see Figure 9C and the caption for how these matrices are constructed from the connectivity matrix).

Figure 11 with 2 supplements see all

Download asset Open asset

Common synapse matrices (A) $S^{K C}$ and (B) $S^{L H N}$ , each of which represents the extent of signal integration from homotypic uPNs to KCs and LHNs.

The color code represents $S_{X Y}$ , which is the number of the third-order neurons (KCs or LHNs) synapsing with both homotypes $X$ and $Y$ . The black color is used when there is no third-order neuron-mediated signal integration ( $S_{X Y} = 0$ ) happening between two homotypes $X$ and $Y$ . See Figure 9C and its caption for how the common synapse matrices are calculated from the connectivity matrices provided in Figure 9—figure supplement 1.

Overall, the number of synaptic connections between uPNs and KCs is greater than that between uPNs and LHNs ( $S_{X Y}^{K C} > S_{X Y}^{L H N}$ , see Figure 11).
In MB calyx, the signals from food-related odors-encoding homotypes (e.g. Yeasty, Fruity, or Alcoholic Fermentation odor types) are shared by a large number of KCs, which constitute a few large clusters in $S^{K C}$ matrix, depicted in red ( $S_{X Y}^{K C} ≳ 35$ ) and indicated by the blue arrows on the top in Figure 11A. Some KCs process signals almost exclusively from the hygro/thermo-sensing homotypes without sharing any signal from other homotypes ( $S_{X Y}^{K C} = 0$ for the cases of $X$ and $Y$ homotype pairs without any signal integration, which are depicted in black in Figure 11). There are also homotypes with significantly less number of overall synaptic connections to KCs, dictated by the diagonal element of the matrix $S^{K C}$ (see Figure 11—figure supplement 2A). In comparison with $S^{L H N}$ , the $S^{K C}$ suggests a stronger but less organized signal integration between heterotypic uPNs by KCs and lends support to the previous literature pointing to the random synapsing of KCs with uPNs at MB calyx (Caron et al., 2013; Stevens, 2015; Eichler et al., 2017; Zheng et al., 2020).
$S^{L H N}$ , on the other hand, demonstrates LHN-mediated signal integration localized to subsets of homotypes. When we collect LHNs connected to a particular homotype and check which other homotypes these LHNs are also synapsing (thereby analyzing the scope of signal integration happening at LH), we find a strong tendency of signals from pheromone and hygro/thermo-sensing uPNs to be integrated within the given odor/signal type (Figure 11). The fact that the pheromone-encoding and hygro/thermo-sensing homotypes share the synaptic connections to LHNs among themselves are demonstrated as the homotype-specific block patterns along the diagonal of the $S^{L H N}$ matrix (see purple and pink arrows on the side in Figure 11B). The $S^{L H N}$ matrix also shows that signals from various food-related odor encoding homotypes, such as DP1l, DP1m, VA2, and VL2p or DM1, DM4, and VA4 are also integrated (see blue arrows in Figure 11B). Many of these homotypes encode signals originating from esters, which is intriguing given the ester-encoding LHN cluster shown by Frechter et al., 2019. The results suggest that certain odor types are processed through common channels of LHNs that are largely dedicated to encoding a particular odor type.

A more stringent selection criterion for synaptic connectivity does not affect our conclusion on the signal integration by the third-order olfactory neurons (Figure 11—figure supplement 2). The only notable change is the general increase in the cases with no integration ( $S_{X Y} = 0$ ) in $S^{L H N}$ , especially for hygro/thermo-sensing homotypes. Thus, the extent of signal integration from homotypic uPNs to KCs and LHNs summarized in $S^{K C}$ and $S^{L H N}$ is robust.

Spatial proximity-based versus connectivity-based clustering

Next, we study the relationship between spatial proximity-based clustering and connectivity-based clustering results. Upon visual inspection, the connectivity-based clustering at MB calyx (Figure 12A on the right) appears less structured than the spatial proximity-based clustering (Figure 12A on the left). Specifically, many homotypic uPNs are grouped under a common branch in the tree structure obtained from the spatial proximity-based clustering, whereas such a feature is largely absent in the output of the connectivity-based clustering. Therefore, the spatially well-clustered uPNs at MB calyx (or stereotyped structure) do not necessarily translate into structured connectivity patterns (or stereotyped connectivity), which is consistent with the notion of randomized PN-KC connections (Caron et al., 2013; Stevens, 2015; Eichler et al., 2017; Zheng et al., 2020). In stark contrast to the outcomes for MB calyx, most homotypic uPNs are grouped in the connectivity-based clustering for LH (Figure 12B). This suggests that the spatially proximal uPNs synapse with a similar group of LHNs. The stereotyped organization and stereotyped connectivity of uPNs in LH have been suggested before (Jefferis et al., 2007; Liang et al., 2013; Kohl et al., 2013; Fişek and Wilson, 2014), and we demonstrate such stereotypies are, in reality, expressed throughout LH over all uPNs. In LH, spatial and organizational characteristics of uPNs are well translated to connectivity to LHNs.

Figure 12

Download asset Open asset

Tanglegrams comparing the tree structures generated from the inter-PN distances-based (left) and the connectivity-based clustering (right) (A) between uPNs and KCs, and (B) between uPNs and LHNs.

The same uPNs in the two tree structures are connected with lines, which visualize where the uPNs clustered by one method end up in the clustering results of another. The labels for uPNs are representative of the homotype and are color-coded based on the encoded odor types (Dark green: decaying fruit, lime: yeasty, green: fruity, gray: unknown/mixed, cyan: alcoholic fermentation, red: general bad/unclassified aversive, beige: plant matter, brown: animal matter, purple: pheromones, pink: hygro/thermo).

A quantitative comparison of two trees based on statistical tests lends support to the notion that the spatial organization of uPNs can be indicative of connective properties, most evident in LH (see Appendix 2 for Baker’s Gamma index, entanglement, and cophenetic distance correlation).

Discussion

The inter-PN organization revealed in this study and its association with odor type/valence are reminiscent of the generally accepted notion that the form determines the function in biology. Previously observed stereotypes of neurons in the Drosophila olfactory system were largely based on the differentiation between pheromones and non-pheromones (Ruta et al., 2010; Kohl et al., 2013; Frechter et al., 2019; Das Chakraborty and Sachse, 2021), the whole-cell patch-clamp recording (Seki et al., 2017), and imaging studies suggestive of stimulus-dependent arrangement of neurons in LH (Marin et al., 2002; Wong et al., 2002; Jefferis et al., 2007). Our results are generally consistent with the previous studies, which suggest that a level of stereotypy in uPN organization in MB calyx and LH is universal throughout Drosophila, which can be captured through different metrics and methodologies. In line with Lin et al., 2007, our study finds that homotypes DL2v and DL2d constitute a bilateral cluster in MB calyx ( $C_{3}^{MB}$ ), and that the dual organization of uPNs is present in MB calyx and LH, such that homotypes DC2, DL1, and VA5 are sorted into the same cluster in LH while sharing similar innervation pattern in MB calyx. Our clustering results in LH share similarities with the NBLAST score-based LH clusters (Bates et al., 2020). The uPNs that ended up in the same cluster or nearby clusters, such as homotypes DM1, DM3, DM4, VA4, and VM3 in the cluster $C_{3}^{LH}$ , are also grouped in the NBLAST score-based clustering analysis (Bates et al., 2020). We find a significant correlation of $d_{α β}$ with NBLAST score (see Figure 2—figure supplement 1) despite the fact that two metrics prioritize different aspects of neuronal anatomy.

Our inter-PN distances and clustering results suggest the spatial organization of uPNs differs greatly in each neuropil (Figure 5). Some of the tightly bundled organization of uPN homotypes are well preserved throughout the neuropils despite the lack of glomerulus in MB calyx and LH. The spatial segregation between different homotypes is, however, practically not present in MB calyx, leading to a high degree of overlapping. Therefore, the heterogeneity of homotypes at the PN-KC synaptic interface may physically assist the randomized sampling known to exist between uPNs and KCs (Caron et al., 2013; Stevens, 2015; Eichler et al., 2017; Zheng et al., 2020).

Our analysis suggests that LH is compartmentalized into four regions: (1) Posterior-dorsal region primarily occupied by food-related uPNs; (2) Anterior-ventral region occupied by pheromone-encoding uPNs; (3) Biforked bundle surrounding posterior-dorsal region largely occupied by food-related uPNs with an aversive response; (4) Posterior-ventral-medial region occupied by hygro/thermo-sensing uPNs. Previous attempts at identifying regions of odorant space in LH revealed compatible results. The three domains (LH-PM, LH-AM, and LH-AL) suggested by Strutz et al., 2014 seem to be a different combination of our clustering result (LH-PM and LH-AM correspond to the posterior-dorsal region and LH-AL corresponds to a combination of anterior-ventral region and the biforked bundle). Although not perfect, the study of the axo-axonic communities in LH yields results with comparable characteristics (Bates et al., 2020), understandably due to the necessity of inter-neuronal proximity to form synapses. For example, the community 12 by Bates et al., 2020 is predominantly composed of homotypes VP1l and DL5, which resembles our cluster $C_{10}^{LH}$ . The community 6 contains a mixture of homotypes VA5, VC1, D, DA4l, DC2, DA3, and VA7m, which is reminiscent of our cluster $C_{6}^{LH}$ .

Many homotypic uPNs that are spatially localized in LH can be associated with key survival features and a strong innate response (Seki et al., 2017). In this sense, the stereotyped localization of pheromone-encoding uPNs in $C_{8}^{MB}$ , $C_{1}^{LH}$ , and $C_{2}^{LH}$ is of great interest. Our study not only lends support to the existing studies pointing to the labeled-line strategy in the Drosophila olfactory system but also suggests that there may exist an even more sophisticated level of spatial organization, which supersedes the pheromone versus non-pheromone segregation. Interestingly, while the spatial organization of uPNs in LH has a basis on the functionality of the odor type encoded, it does not seem to be directly translated to segregated chemical features seen in LHNs (Frechter et al., 2019). The apparent divergence observed at the PN-LHN interface, coupled with strongly stereotyped connectivity may contribute to a higher resolution of odor categorization.

The uPN organizations from FAFB and hemibrain datasets are consistent

Our analyses of both the FAFB and the hemibrain dataset (Scheffer et al., 2020) find that the results from both datasets are generally consistent. For example, ${\bar{d}}_{intra}$ , ${\bar{d}}_{inter}$ , and $λ$ analyzed based on two different datasets are almost identical (see Figure 5A and B). ${\bar{d}}_{intra,X}$ , ${\bar{d}}_{inter,X}$ , and $λ_{X}$ show slight differences due to a mismatch between the FAFB and the hemibrain dataset (on glomerulus labels and the number of uPNs based on our selection criterion) leading to a different number of uPNs per homotype (Figure 6A and B and Figure 5—figure supplement 1), but the correlation between $λ_{X}$ s at MB calyx and LH are still observed (Figure 7C and Figure 7—figure supplement 1C). Most importantly, the clustering results are similar, where many spatial clusters in both datasets share the same set of homotypes. Additionally, odor type-dependent spatial properties are retained (Figure 8 and Figure 8—figure supplement 1), with all statistical tests supporting our hypothesis. In conclusion, the outcomes from our analyses of the two EM datasets lend support to the previous claims of stereotypy in the Drosophila brain and neuronal structures (Jenett et al., 2012; Jeanne et al., 2018; Schlegel et al., 2021).

Odor signal processing and labeled-lines

Our study suggests that while the primary connectivity motif of third-order olfactory neurons indeed integrates signals, there still exist several labeled lines. The synaptic connections at the PN-KC interface in MB calyx are largely integrative and randomized - with an exception of hygro/thermo-sensing homotypes that display stereotypy even in terms of the connectivity to the KCs. A similar observation has been made by Li et al., 2020, who employed NBLAST score to identify a structural segregation between odor-encoding and hygro/thermo-sensing homotypes. They found that specific KC types are preferentially targeted by hygro/thermo-sensing homotypes. Marin et al., 2020, who carried out connectome analysis specific to hygro/thermo-sensing homotypes, also identified that lateral accessory calyx (lACA), the anterior-dorsal part of the calyx, are primarily targeted by hygro/thermo-sensing homotypes (analogous to our clusters $C_{1}^{MB}$ and $C_{4}^{MB}$ in Figure 3), and found that a number of KCs are dedicated to encoding signals from these homotypes. The uPNs in LH are spatially segregated, which translates to connectivity in three different levels. First, certain LHNs are dedicated to encoding signals from a specific homotype. The number of these ‘homotype-specific’ LHNs varies across the homotype and can make up a significant portion of PN-LHN connections depending on the homotype (Figure 10). Second, synaptic connectivity maps between uPNs and LHNs indicate odor type-dependent integration occurs at LH (Figure 11B). Channels of LHNs predominantly encoding specific odor types are observed; one primarily integrates responses from certain food-related homotypes, one integrates pheromone-encoding homotypes, and another integrates hygro/thermo-sensing homotypes. Third, homotypic uPNs share similar connectivity to LHNs, unlike those in MB calyx. The signals relayed from the spatially well-organized (or tightly bundled) homotypes are localized into a specific group of LHNs, thereby forming a ‘homotype-specific’ connectivity motif (Figure 10, Figure 11, and Figure 12).

In our study of the labeled-line strategy, we made several interesting observations, which are worth comparing with the concept of ‘fovea’ introduced by Zheng et al., 2020. A ‘fovea’ delineates deviations between experimentally observed connectivity matrices and connectivity under the assumption of random synapses in MB calyx, specifically for certain food-related uPNs (Zheng et al., 2020). A group of common KCs predominantly sampling ’food-related’ uPNs manifest themselves in the common synapse matrix $S^{K C}$ (see the group of homotypes comprising the clusters, highlighted by the blue arrows in Figure 11A). A subset of homotypic uPNs under the food-related ‘fovea’ reported by Zheng et al., 2020 are also spatially clustered (e.g. DM1, DM4, DP1m, DP1l, VA2, and VA4). While most of these homotypes are spatially proximal (the vast majority of the uPNs are located in clusters $C_{6}^{MB}$ and $C_{7}^{MB}$ ), some homotypes under the food-related ‘fovea’ such as VA2 are sampled from spatially disparate clusters. Thus, it is likely that factors other than the spatial organization of uPNs in neuropils contribute to creating the ‘fovea’. Interestingly, the spatial proximity of pheromone-encoding homotypes in MB calyx may suggest the existence of pheromone-encoding ‘fovea,’ but most uPNs in these homotypes do not converge in connectivity-based clustering with an exception of VA1d. In fact, we suspect the spatial organization of pheromone-encoding homotypes in MB calyx, which is placed at the center of the neuropil, to facilitate the observed randomization of connections by increasing the accessibility of KCs to these homotypes. There is, however, a potential hygro/thermo ‘fovea,’ where homotypes such as VP1d and VP2 are spatially clustered together and the signals from these homotypes are relayed by the same set of KCs. Curiously, VL1 is part of this hygro/thermo ‘fovea’ (Figure 12A).

Multiglomerular PNs are spatially distinctive

Apart from uPNs primarily explored in this study, a host of local neurons (LNs) and multiglomerular PNs (mPNs) also constitute sophisticated neural circuits to regulate the signals received from ORNs (Sudhakaran et al., 2012; Bates et al., 2020), playing a significant role in the olfactory signal processing (Olsen et al., 2010; Jeanne and Wilson, 2015; Seki et al., 2017). A large portion of these mPNs is GABAergic and inhibitory (Berck et al., 2016; Tobin et al., 2017), whereas the role of interneurons can be both inhibitory and excitatory (Wilson et al., 2004; Turner et al., 2008). Electrophysiological measurements indicate that mPNs are narrowly tuned to a specific set of odor stimuli (Berck et al., 2016), which is significant given that PNs are generally thought to be more broadly tuned than presynaptic ORNs (Wilson et al., 2004). Several PNs do not follow the typical mALT, but mediolateral antennal lobe (mlALT) or lateral antennal lobe tracts (lALT) instead, thereby bypassing innervation through one of the higher olfactory centers (Schultzhaus et al., 2017; Zheng et al., 2018; Bates et al., 2020). As stated previously, we confined ourselves to uPNs innervating all three neuropils to compare the spatial organization across neuropils for each uPN. As a result, 28 uPNs present in the FAFB dataset are not explored in our study. In MB calyx, only two uPNs constituting VP3 were dropped as a result of our selection criterion, which ended up in an almost identical clustering output once hierarchical clustering was performed on the entire 137 uPNs that innervate MB calyx. Two missing uPNs were grouped into clusters $C_{4}^{MB}$ and $C_{6}^{MB}$ , along with other hygro/thermo-sensing homotypes. On the other hand, the addition of 27 uPNs constituting 15 homotypes innervating LH but not MB calyx created four new clusters when hierarchical clustering was performed (Figure 4—figure supplement 1). The additional uPNs changed the content of the individual clusters; that is, the tree-cutting algorithm broke down a few clusters that became larger due to the additional uPNs. Furthermore, when we calculate the ${\bar{d}}_{intra}$ , ${\bar{d}}_{inter}$ , and $λ$ in LH for the 15 homotypes that included the 27 uPNs, we find that the ${\bar{d}}_{intra}$ values increased when the 27 uPNs were included (see Figure 5—figure supplement 2). This suggests that the previously removed uPNs, most of which follow mlALT, are significantly different in terms of spatial and organizational characteristics and thus should be analyzed separately. Out of 27 additional uPNs in LH, 21 were in mlALT, 5 were in trans-lALT, and 1 was in mALT. Figure 4—figure supplement 2 illustrates how these 27 uPNs innervate LH which demonstrates the reason behind increased ${\bar{d}}_{intra}$ values. These 27 uPNs are mostly GABAergic (21 are labeled as GABAergic, 1 as cholinergic, and 4 as unknown neurotransmitter type), covering 84% of all GABAergic uPNs available in the FAFB dataset. These uPNs innervate LH differently from other uPNs in the same homotype that follow mALT (see homotypes such as DA1, DC4, DL2d, DL2v, DP1l, VA1d, VA1v, VL2a, VL2p, and VP5 in Figure 4—figure supplement 2). Morphologically, inhibitory GABAergic neurons are often considered ‘smooth’ and aspiny (Douglas et al., 1989; Bopp et al., 2014; Gouwens et al., 2019), which are discernible from Figure 4—figure supplement 2.

The single-uPN homotypes may have different morphological properties

It is of great interest that many of the single-uPN homotypes, characterized by densely branched morphology, encode signals with aversive responses. Direct transmission of the associated signals across the three neuropils via a single PN might simplify the overall processing of the olfactory signals as well as reduce the energetic cost. Similarly, the morphological characteristics of uPN innervation at each neuropil are intriguing. Even though a structural difference exists between the single-uPN and multi-uPN homotypes, all uPN innervations within neuropil share a similar morphology regardless of the homotype (see Figure 6—figure supplement 1; Choi et al., 2022). A localized morphological diversity within a neuron may be a characteristic aspect of pseudo-unipolar neurons like uPN and suggests a fundamentally multi-scale characteristic of neuron morphology.

The Drosophila brain EM reconstruction project has evolved to its near completion since the EM image dataset was first released (Dorkenwald et al., 2022). The reconstruction of the majority of the Drosophila central brain as well as the corresponding connectome with detailed information of the individual synapses has become publicly available (Scheffer et al., 2020). Our analysis of the second-order neurons inside the Drosophila olfactory system may be translated to other parts of the nervous system in Drosophila as well as different organisms including the central nervous system (CNS) of humans. For the mammalian olfactory system, the details of analyses must be adapted, however, since the wiring scheme is much more complex than that of an insect (Maresh et al., 2008). For example, multiple glomeruli encoding the same olfactory signal exist in humans (Mombaerts et al., 1996). When analyzing the spatial properties, this can be accounted for by prioritizing the individual glomerulus over the homotypes. Then, homotypic PNs forming different glomeruli may be compared or averaged if one were to consider the homotype-dependent characteristics. According to the neurotransmitter map from a recent study (Dolan et al., 2019), sophisticated processes beyond neuronal anatomy are apparently at work in the olfactory signal processing. Thus, functional studies incorporating odor response profiles in PNs (Badel et al., 2016) and ORNs (Münch and Galizia, 2016; Bak et al., 2018) would supplement our findings. The extension of our study to the other regulatory interneurons and mPNs, morphological studies of second-order neurons, and spatial analysis of third-order neurons will be of great interest for a better understanding of the olfactory signal processing beyond the implication of the neural anatomy and connectivity studied here.

Materials and methods

Data preparation

Request a detailed protocol

We used the neuron morphology reconstruction of 346 Drosophila olfactory neurons from the FAFB dataset (Bates et al., 2020) traced from EM images. The neurons were extracted from the right hemisphere of the female Drosophila. Out of 346 neurons in the FAFB dataset, 164 neurons were uPNs. One uPN in the dataset (neuron ID = 1356477 forming VP3) did not have an associated reconstruction (.swc file) available and was therefore ignored. For this study, uPNs that innervate all three neuropils were chosen because our aim is (1) to compare spatial characteristics of the uPN innervation across each neuropil and (2) to classify each uPN based on the odor encoding information. Thus, out of 164 uPNs, a total of 135 uPNs constituting 57 homotypes were collected under this criteria, resulting in mostly cholinergic uPNs that follow mALT. Rest of the uPNs that did not innervate all three neuropils are collected for the supplementary analysis. In MB calyx, a total of 137 PNs are identified with two PNs constituting VP3 that do not innervate all three neuropils. On the other hand, in LH, a total of 162 PNs are identified, indicating that 27 PNs constituting 15 homotypes do not innervate all three neuropils. The morphological information of each neuron is stored as a set of 3D coordinates with the connectivity specified with the parent nodes. Complete reconstruction of neuron morphology was made by connecting data points based on their parent-child relationship.

The hemibrain dataset (Scheffer et al., 2020) was taken from the neuPrint database (Clements et al., 2020), where we collected from the right hemisphere of the female Drosophila a total of 120 uPNs forming 58 glomeruli based on the same criterion we used for the FAFB dataset (uPNs that innervate all three neuropils). Unlike the FAFB dataset, the neurons in the hemibrain dataset are labeled with regions of interest (ROI), which are used to query uPNs conforming to our selection criterion. The discrepancy in the number of uPNs between the two datasets most likely resulted from the difference between the neuropil boundary we used and the region defined by the hemibrain dataset. In fact, we find that the total number of uPNs in both datasets is comparable, with 164 uPNs in the FAFB dataset and 162 uPNs in the hemibrain dataset. The two datasets also had a minor mismatch in the glomerulus label annotations, sometimes affecting the number of uPNs constituting a given homotype. Among the 120 uPNs from the hemibrain dataset, five uPNs had ambiguity in terms of their glomerulus labels, which is presumably due to poorly formed glomerular structures. For these uPNs, we adopted the glomerulus labels of the FAFB dataset with the matching hemibrain neuron IDs.

Additionally, a recent community-led effort identified three glomeruli in both databases with conflicting glomerulus labels, which have been a source of confusion (Schlegel et al., 2021). After an extensive study, the community agreed to rename the glomeruli in both datasets labeled as VC3l, VC3m, and VC5 as VC3, VC5, and VM6, respectively (Schlegel et al., 2021). Thus, we have manually incorporated these labels into our analyses for both the FAFB and the hemibrain dataset.

Next, we systematically demarcated the regions of AL, LH, and MB calyx. The density of data points projected to each axis was used for the identification since the neuropils are featured with a much higher density of data points than the rigid backbone connecting them. The boundaries defining each neuropil were systematically chosen from local minima that separate neuropils from rigid backbones. Due to the unique structure of uPNs, sometimes the projection along a given axis cannot fully differentiate two neuropils. To resolve this issue, projections along each axis were sampled while rotating the data points along the reference axes at $5^{\circ}$ increments to obtain multiple snapshots. The densities were analyzed to choose the optimal degrees of rotation along the reference axes that could best segment the neuropils. We used the smallest average and deviation value of density at the local minima as the criteria to choose the optimal rotation. The process has been repeated for each neuropil to identify a set of boundaries along multiple transformed axes with various degrees of rotations that optimally confine each neuropil. This information has been combined to create a set of conditions per neuropil for segmentation. The resulting neuropils were confirmed through visual inspection. We compared our neuropil segmentation boundaries with neuropil volume surface coordinates provided by Ito et al., 2014 via CATMAID (Saalfeld et al., 2009) and found the boundaries are comparable (data not shown). An overview of the segmentation process is available in Figure 13.

Figure 13

Download asset Open asset

A diagram depicting the neuropil segmentation process.

The data points from skeletal reconstruction are projected to each axis to generate distributions from which local minima are obtained. The process is repeated while rotating the uPNs along each axis. A collection of histograms and corresponding local minima are surveyed to generate a set of optimal rotations and boundaries for individual neuropil. The resulting parameters are combined to form a collection of conditions to segment each neuropil.

The odor type and odor valence information were extracted from various literature (Hallem et al., 2004b; Galizia and Sachse, 2010; Stensmyr et al., 2012; Mansourian and Stensmyr, 2015; Badel et al., 2016; Bates et al., 2020) and we closely followed the categorical convention established by Mansourian and Stensmyr, 2015 and Bates et al., 2020. However, we note that the categorization of a uPN under a specific odor category may overshadow the complete spectrum of odorants a uPN might encode, especially if the uPN encodes ORs that are broadly tuned. Therefore, we focused on the well-separated pheromone/non-pheromone encoding types and valence information.

To test our labeled-line hypothesis, the connectivity information between uPNs and higher olfactory neurons such as KCs and LHNs was necessary. Since only the hemibrain dataset contains detailed connectivity information, all of our connectivity studies are done using uPNs, KCs, and LHNs queried from the hemibrain dataset. We chose KCs and LHNs that made at least three synaptic connections with any of the 120 uPNs from the hemibrain dataset. This resulted in 1754 KCs and 1295 LHNs, creating bipartite connectivity matrices at each neuropil.

Inter-PN distance

Request a detailed protocol

The ‘distance’ $d_{α β}$ between two neurons, $α$ and $β$ , with different lengths ( $N_{α} \leq N_{β}$ ) is quantified by calculating.

d_{α β}^{2} = \frac{1}{N_{α}} \sum_{i = 1}^{N_{α}} min [(r_{i}^{α} - r_{j}^{β})^{2}],

where $r_{i}^{α}$ is an i-th coordinate forming the neuron $α$ . Equation 1 is evaluated over all pairs of $r_{i}^{α}$ and $r_{j}^{β}$ with $j = 1, \dots, N_{β}$ that gives rise to the minimum value. This means that when $N_{α} \leq N_{β}$ , for every i-th coordinate in the neuron $α$ ( $r_{i}^{α}$ ), we find j-th coordinate in the neuron $β$ ( $r_{j}^{β}$ ) that is the closest to $r_{i}^{α}$ . Then, the spatial proximity of a given pair of neurons is assessed by the $d_{α β}$ that denotes the average of all the minimum Euclidean distances between the pair of coordinates.

The degree of bundling, packing, and overlapping

Request a detailed protocol

We define the mean intra- and inter-homotype distances as.

{\bar{d}}_{i n t r a, X} \equiv \frac{1}{N} \sum_{α, β \in X}^{N} d_{α β}

and

{\bar{d}}_{i n t e r, X} \equiv \frac{1}{N} \sum_{α \in X, β \notin X}^{N} d_{α β},

where $X$ denotes a homotype and $N$ is the total number of uPN pairs to be averaged. The ${\bar{d}}_{intra, X}$ is calculated over all the pairs of uPNs in the $X$ -th homotype, quantifying the tightness of bundling of uPNs constituting the $X$ -th homotype. On the other hand, ${\bar{d}}_{inter, X}$ is calculated over the pairs of uPNs between $α$ -th uPN belonging to the $X$ -th homotype and $β$ -th uPN in the $Y$ -th homotype ( $Y \neq X$ ), such that it measures the extent of packing of uPNs around the $X$ -th homotype. The degree of overlapping for the $X$ -th homotype, $λ_{X}$ , is defined as the ratio of average intra- and inter-homotype distances,

λ_{X} = \frac{{\bar{d}}_{i n t r a, X}}{{\bar{d}}_{i n t e r, X}},

which represents how clearly the $X$ -th homotype is segregated from other homotypes in a given space. A large value of $λ_{X}$ ( $λ_{X} ≫ 1$ ) implies that the space spanned by the $X$ -th homotype is not clearly discerned from other homotypes.

Spatial clustering of projection neurons

Request a detailed protocol

Hierarchical/agglomerative clustering was used to cluster the uPN innervation at each neuropil using the pairwise $d_{α β}$ matrices. First, the linkage was decided based on the pairwise distance matrix built with the Farthest Point Algorithm (or ‘complete’ method), where uses the maximum distance between neurons to define the distance between two clusters. This criterion is used to build hierarchical relations (or nested clusters) in a bottom-up approach where each neuron is treated as a cluster at the beginning. The result is a fixed tree structure of individual neurons from which the finalized clusters are formed using an optimal tree-cutting algorithm. In the dendrogram from AL (Figure 2—figure supplement 3), homotypic uPNs are grouped together with high accuracy, suggesting our distance metric $d_{α β}$ is adequate. We tested various tree-cutting criteria such as elbow method, gap statistics, maximum average silhouette coefficient, and dynamic hybrid cut tree method (Langfelder et al., 2008) to determine the optimal number of clusters. Among them, we selected the dynamic hybrid cut tree method, since it performed the best in giving the cluster number closest to the number of different odor types (which is 10) (Table 1). We deployed the dynamic hybrid cut tree method with the minimum cluster size of 4 neurons for the tree-cutting, following the neuron clustering procedure used by Gouwens et al., 2019.

Table 1

The optimal number of clusters of uPNs in the FAFB dataset determined by employing the dynamic hybrid cut tree method, elbow method, gap statistics, and maximum average silhouette coefficient.

	Dynamic hybrid	Elbow	Gap	Silhouette
AL	19	14	8	54
MB calyx	10	11	7	2
LH	11	9	7	2

Pearson’s $x^{2}$ -test of independence

Request a detailed protocol

The association between two categorical variables is assessed using Pearson’s $χ^{2}$ -test. For the test, a contingency table, which lists the categorical frequency of two variables, is created. For example, $O_{i j}$ of the $i$ - and $j$ -th element of the contingency table shown below is the frequency counting the putative valence $i = 1$ (attractive), 2 (aversive), 3 (unknown), and the number of uPNs in one of the 10 clusters in MB calyx with $j = 1$ ( $C_{1}^{MB}$ ), 2 ( $C_{2}^{MB}$ ), ... , 10 ( $C_{10}^{MB}$ ).

	$C_{1}^{M B}$	$C_{2}^{M B}$	$C_{3}^{M B}$	$C_{4}^{M B}$	$C_{5}^{M B}$	$C_{6}^{M B}$	$C_{7}^{M B}$	$C_{8}^{M B}$	$C_{9}^{M B}$	$C_{10}^{M B}$	Total
Attractive	0	4	0	1	0	5	4	11	11	8	44
Aversive	1	2	0	0	4	12	9	8	8	3	47
Unknown	4	7	8	5	6	5	1	2	3	3	44
Total	5	13	8	6	10	22	14	21	22	14	135

Then the $χ^{2}$ value is evaluated based on the table using.

χ^{2} = \sum_{i = 1}^{R} \sum_{j = 1}^{C} \frac{(O_{i j} - E_{i j})^{2}}{E_{i j}},

where $R$ and $C$ are the numbers of rows and columns, and $O_{i j}$ and $E_{i j}$ are the observed and expected frequencies of the event in the $i$ -th row and $j$ -th column, respectively. $E_{i j}$ is calculated from $O_{i j}$ as.

E_{i j} = N p_{i \cdot} p_{\cdot j},

where $p_{i, \cdot} = \sum_{j}^{C} O_{i j} / N$ and $p_{\cdot j} = \sum_{i}^{R} O_{i j} / N$ with $N$ being the total count. Thus, $E_{i j}$ is the frequency expected by assuming that the two categorical data are statistically independent. Pearson’s $χ^{2}$ test aims to check whether there is a significant difference between $O_{i j}$ and $E_{i j}$ .

In the $χ^{2}$ -test, the p-values are estimated using $f_{k} (x)$ , the $χ^{2}$ -distribution with the degree of freedom $k = (R - 1) (C - 1)$ . If the test returns a $χ^{2}$ value that gives rise to a p-value smaller than the defined significance level ( $α = 0.01$ ), the null hypothesis of independence between the two data sets should be rejected. As a result, the distribution of the categorical data is deemed significantly different from a randomly generated distribution, which concludes that the association between two sets of data is statistically significant.

For the above contingency table with $k = 18$ , which leads to $χ^{2} \approx 66.1$ (Equation 5), we get a p-value much smaller than the significance level ( $α = 0.01$ ), $p = 1 - \int_{0}^{χ^{2}} f_{k = 18} (x) d x \approx 2.016 \times 10^{- 7} ≪ α = 0.01$ .

When Pearson’s $χ^{2}$ statistics are available, one can calculate Cramér’s $V$ with bias correction, a measure of association between two categorical variables, as follows.

V = \sqrt{\frac{ϕ^{' 2} / N}{min (R^{'} - 1, C^{'} - 1)}},

where $ϕ^{', 2} = \max (0, χ^{2} / N - (R - 1) (C - 1) / (N - 1))$ , $R^{'} = R - {(R - 1)}^{2} / (N - 1)$ , and $C^{'} = C - {(C - 1)}^{2} / (N - 1)$ . Similar to the Pearson correlation coefficient, the value $V$ ranges between 0 and 1 where 0 indicates no correlation and 1 indicates a complete correlation between two categorical variables.

Mutual information

Request a detailed protocol

Mutual information ( $I$ ) is used to verify the significance of association between nominal variables observed in Pearson’s $χ^{2}$ -test for independence. The $I$ measures the information transfer or the similarity between two data. The concept can be extended to clustering outputs to check how two different clustering labels from the same data are similar to each other. Traditionally, the $I$ between two jointly discrete variables $A$ and $B$ is given by.

I (A; B) = \sum_{i = 1}^{n_{A}} \sum_{j = 1}^{n_{B}} P (A_{i}, B_{j}) \log [\frac{P (A_{i}, B_{j})}{P (A_{i}) P (B_{j})}],

where $n_{A}$ (or $n_{B}$ ) is the number of clusters in $A$ (or $B$ ). Numerically, the $I$ between two clustering outputs $A$ and $B$ is calculated by evaluating $P (A_{i}) = N_{A_{i}} / N$ , $P (B_{i}) = N_{B_{i}} / N$ , and $P (A_{i}, B_{j}) = N_{A_{i} \cap B_{j}} / N$ where $N$ is the total count and $N_{A_{i} \cap B_{j}}$ is the number of elements common in both clusters $A_{i}$ and $B_{j}$ .

The significance was assessed by comparing the observed $I$ with the distribution of $I$ s from randomly sampled variables. Specifically, the cluster label was randomly sampled 1000 times to generate a distribution of $I$ under the assumption of independence. The value of observed $I$ is considered significant if the approximated p-value is below 0.01 (p< 0.01).

Analysis of synaptic interfaces

Request a detailed protocol

We conducted three different analyses on the synaptic interfaces of uPNs with the third-order neurons (KCs or LHNs) from the hemibrain dataset.

(i) The ‘homotype-specific’ connections ( $N_{X, sp}^{ξ}$ with $ξ =$ PN-KC or PN-LHN) are obtained by counting the number of third-order neurons that synapse with a homotype $X$ but do not synapse with any other homotypes, the information of which is provided by the binarized connectivity matrix $C$ . The total number of synaptic connections for a homotype $X$ is simply the sum of the row of the connectivity matrix $C$ ( $N_{X, t o t}^{ξ} = Σ_{i = 1}^{N_{ξ}} C_{X i}$ ).

(ii) To generate the $S$ matrices, we counted the number of third-order neurons synapsing with a given homotype $X$ that also synapses with other homotypes.

(iii) The tanglegram study required a hierarchical clustering of uPNs based on their connectivity to third-order neurons. The distances between uPNs in the connectivity matrix $C$ represent the similarity of the connectivity patterns to third-order neurons between two uPNs. We utilized the metric of cosine distance, which is widely used for analyzing the connectivity matrix (Bates et al., 2019; Bates et al., 2020; Li et al., 2020; Eschbach et al., 2020; Schlegel et al., 2021). The cosine distance is defined as.

d_{c o s} = 1 - \frac{u \cdot v}{| u | | v |},

where $u$ and $v$ are two vectors to be compared. After calculating the distances, we performed hierarchical clustering by Ward’s criterion, which minimizes the variance of merged clusters, to generate the tree structure. The results of hierarchical clustering using the spatial proximity ( $d_{α β}$ ) and connectivity ( $d_{c o s}$ ) are compared using a tanglegram (Figure 12) after untangling two trees using the ’step1side’ method (Galili, 2015).

Appendix 1

Monte Carlo approach to independence test

In this section, we describe an alternative method to the independence test inspired by the Monte Carlo significance test (Hope, 1968) to further support our Pearson’s $χ^{2}$ -test of independence. The procedure is as follows: (1) For a given contingency table, randomize the observation such that the marginal sum of each row remains the same as the observed contingency table. That is, for each row, randomize the vector with integers while the sum of the vector stays the same as the observed contingency table. This procedure randomly shuffles the distribution of the clusters while keeping the distribution of a particular categorical variable intact. (2) Calculate the $χ^{2}$ value from the randomized contingency table. (3) Repeat steps 1 and 2 for 1000 times to generate a distribution of the $χ^{2}$ values. (4) Obtain the mean and the standard deviation of $χ^{2}$ values. The distribution of $χ^{2}$ values is approximately normal. (5) If the $χ^{2}$ value from the observed contingency table is more than $4 σ$ different from the distribution, we consider the observed $χ^{2}$ value statistically significant and reject the null hypothesis. Whenever we ran a Pearson’s $χ^{2}$ -test, we performed the above procedure alongside. The output of this procedure supported whichever conclusion we drew from Pearson’s $χ^{2}$ -test.

Identifying the agreement between two categorical data via mutual information

We verified our Pearson’s $χ^{2}$ -test of independence of two categorical variables by calculating the mutual information $I$ (see Methods). In the FAFB dataset, the mutual information between glomerular labels and $d_{α β}$ -based clustering output in MB calyx was equal to $I (glo; C^{MB}) = 1.892$ , which is significantly (more than $4 σ$ ) different from the mean of randomly sampled $I$ distribution, $I {(glo; C^{MB})}_{rand} = 1.386 \pm 0.035$ . This result is consistent with our $χ^{2}$ -test, as the mutual information of the observed variables is significantly larger than the mutual information under the assumption of random sampling, suggestive of a statistically significant association between glomerular labels and MB calyx cluster labels. In LH, the mutual information between glomerular labels and $d_{α β}$ -based clustering output was $I (glo; C^{LH}) = 2.128$ which deviated $4 σ$ or more from the mean of the randomly sampled $I$ distribution, $I {(glo; C^{LH})}_{rand} = 1.466 \pm 0.035$ .

The same method is applied to confirm that a statistically significant association exists between odor type and the clustering outputs, with $I (odor; C^{MB}) = 0.819$ and $I (odor; C^{LH}) = 0.963$ , all of which differ by more than $4 σ$ from the means of the randomly sampled $I$ distributions, $I {(odor; C^{MB})}_{rand} = 0.337 \pm 0.044$ and $I {(odor; C^{LH})}_{rand} = 0.372 \pm 0.043$ . For odor valence, we obtain $I (val; C^{MB}) = 0.277$ and $I (val; C^{LH}) = 0.326$ , where both $I (val; C^{MB})$ and $I (val; C^{LH})$ differ significantly from the means of the randomly sampled $I$ distributions, $I {(val; C^{MB})}_{rand} = 0.073 \pm 0.026$ and $I {(val; C^{LH})}_{rand} = 0.081 \pm 0.026$ . Overall, the conclusion drawn from the association study based on mutual information is identical to Pearson’s $χ^{2}$ -test.

Appendix 1—table 1

Pearson’s $χ^{2}$ tests of independence of variables in the FAFB dataset.

$C^{Z}$ indicates cluster labels from $d_{α β}$ -based clustering in $Z$ neuropil. Cramér’s V values are displayed on each cell and the corresponding p-values are shown in parentheses.

	$C^{L H}$	Glomerular Labels	Odor Type	Odor Valence
$C^{M B}$	0.502 (1.149E-36)	0.610 (1.255E-27)	0.401 (3.303E-21)	0.425 (2.016E-07)
$C^{L H}$		0.671 (2.266E-40)	0.416 (1.980E-22)	0.455 (2.586E-08)

Appendix 1—table 2

Pearson’s $χ^{2}$ tests of independence of variables in the hemibrain dataset.

$C^{Z}$ indicates cluster labels from $d_{α β}$ -based clustering in $Z$ neuropil. Cramér’s V values are displayed on each cell and the corresponding p-values are shown in parentheses.

	$C^{L H}$	Glomerular Labels	Odor Type	Odor Valence
$C^{M B}$	0.495 (6.635E-34)	0.577 (3.461E-25)	0.425 (9.400E-18)	0.463 (6.283E-07)
$C^{L H}$		0.685 (1.523E-40)	0.502 (6.072E-29)	0.521 (2.932E-09)

Appendix 1—table 3

Mutual information (observed mutual information (top), randomly sampled mutual information (bottom) in each cell) from the association study using the FAFB dataset.

$C^{Z}$ is cluster labels from $d_{α β}$ -based clustering at $Z$ neuropil. The observed mutual information differs from the randomly sampled mutual information by more than $4 σ$ .

	$C^{L H}$	Glomerular Labels	Odor Type	Odor Valence
$C^{M B}$	1.076 0.397±0.045	1.892 1.386±0.035	0.819 0.337±0.044	0.277 0.073±0.026
$C^{L H}$		2.128 1.466±0.035	0.963 0.372±0.043	0.326 0.081±0.026

Appendix 1—table 4

Mutual information (observed mutual information (top), randomly sampled mutual information (bottom) in each cell) from the association study using the hemibrain dataset.

$C^{Z}$ is cluster labels from $d_{α β}$ -based clustering at $Z$ neuropil. The observed mutual information differs from the randomly sampled mutual information by more than $4 σ$ .

	$C^{L H}$	Glomerular Labels	Odor Type	Odor Valence
$C^{M B}$	1.371 0.710±0.048	2.244 1.783±0.033	1.036 0.527±0.047	0.336 0.124±0.035
$C^{L H}$		2.344 1.717±0.034	1.211 0.493±0.048	0.434 0.116±0.033

Appendix 1—table 5

Statistics of homotypes composed of a single uPN (or multiple uPNs) in the FAFB dataset and the corresponding putative valence.

	Aversive	Attractive	Unknown	Total
Single uPN Homotypes Count	7	4	2	13
Multiple uPN Homotypes Count	18	13	13	44
Total	25	17	15	57

Appendix 1—table 6

Statistics of homotypes composed of a single uPN (or multiple uPNs) in the hemibrain dataset and the corresponding putative valence.

	Aversive	Attractive	Unknown	Total
Single uPN Homotypes Count	7	5	1	13
Multiple uPN Homotypes Count	18	12	15	45
Total	25	17	16	58

Appendix 2

Testing the labeled-line hypothesis

We detail the analyses performed on the tanglegram and the respective outputs (Figure 12). First, we applied the dynamic hybrid cut tree method on the dendrogram generated from connectivity and conducted Pearson’s $χ^{2}$ test. The results are shown in Table 1. The p-values for the connectivity-based clustering between uPNs and LHNs for glomerular labels, odor types, and odor valence were very small. For the connectivity between uPNs and KCs, we see a moderate to no association for the given categorical variables (Appendix 2—table 1).

The similarity between two tree structures from spatial proximity-based and connectivity-based clustering at a given synaptic interface is measured in several different ways to provide a comprehensive comparison. First, we quantified the similarity using Baker’s Gamma index (Baker, 1974), which is a measure of rank correlation (or ordinal relation) calculated from concordant and discordant pairs given by.

G_{B a k e r} = \frac{N_{c o n} - N_{d i s}}{N_{c o n} + N_{d i s}},

where $N_{con}$ is the number of concordant pairs (the ordering of elements in two trees match) and $N_{dis}$ is the number of discordant pairs (the ordering of elements in two trees do not match). Baker’s Gamma index ranges from -1 to 1 where 0 represents the ordering of two trees is completely dissimilar and 1 or -1 indicate the ordering of two trees match. We find $G_{Baker}^{MB} = 0.286$ and $G_{Baker}^{LH} = 0.219$ (which we double-checked using both the in-house code and ‘dendextend’ package in R). Baker’s Gamma index for LH is very similar to the one obtained by Bates et al., 2020 ( $G_{Baker}^{LH} = 0.21$ ), who conducted a similar study using the NBLAST score and connectivity. However, the fact that $G_{B a k e r}^{M B} > G_{B a k e r}^{L H}$ when the tanglegram of MB calyx is seemingly more incoherent (Figure 12A) raises a question of whether Baker’s Gamma index alone is enough to describe the tanglegram.

Apart from the ordinal relations between two sets of leaves, we employed two additional metrics to compare the two trees: (1) entanglement, a measure spanning from 0 to 1 quantifying the number of lines crossing, and (2) cophenetic distance correlation, a measure spanning from 0 to 1 quantifying how similar the two branching structures are. The entanglement between two trees for MB calyx was 0.35 (higher entanglement), while the entanglement for LH was 0.26 (lower entanglement), which agrees with Figure 12. To calculate cophenetic distance correlation, we measured the pairwise cophenetic distances within each tree and calculated the Pearson correlation coefficient. The cophenetic distance between two leaves in the dendrogram is equal to the minimum distance (or height) to the branching point that contains both leaves. The Pearson correlation coefficient between cophenetic distances of the spatial proximity-based and connectivity-based tree structures was $r = - 0.032$ ( $p > 0.001$ ) for MB calyx and $r = 0.236$ ( $p ≪ 0.001$ ) for LH, reflecting the less disrupted tree structure in LH compared to MB calyx.

Appendix 2—table 1

Pearson’s $χ^{2}$ tests of independence of variables on the connectivity-based clustering results.

Cramér’s V values are displayed on each cell and the corresponding p-values are shown in parentheses. Bold entries are used to specify statistically significant results.

	Glomerular Labels	Odor Type	Odor Valence
$C^{P N - K C}$	0.433 (2.472E-08)	0.316 (9.978E-09)	0.271 (0.012)
$C^{P N - L H N}$	0.765 (1.410E-67)	0.630 (1.519E-48)	0.604 (4.055E-12)

Data availability

All data generated during this study and Python script are available in Drosophila Olfaction-main.zip included as the supporting file. They are also available at https://github.com/kirichoi/DrosophilaOlfaction, (copy archived at swh:1:rev:91dd60f4231a58590e2571e72b660c5dfee261b6).

The following previously published data sets were used

1. Bates AS
2. Schlegel P
3. Roberts RJV
4. Drummond N
5. Tamimi IFM
6. Turnbull R
7. Zhao X
8. Marin EC
9. Popovici PD
10. Dhawan S
11. Jamasb A
12. Javier A
13. Capdevila LS
14. Li F
15. Rubin GM
16. Waddell S
17. Bock DD
18. Costa M
19. Jefferis GSXE
(2020) FAFB
ID catmaid. Complete Connectomic Reconstruction of Olfactory Projection Neurons in the Fly Brain.

https://fafb.catmaid.virtualflybrain.org/
1. Scheffer LK
2. Xu CS
3. Januszewski M
4. Lu Z
5. Takemura SY
6. Hayworth KJ
7. Huang GB
8. Shinomiya K
9. Maitlin-Shepard J
10. Berg S
11. Clements J
12. Hubbard PM
13. Katz WT
14. Umayam L
15. Zhao T
16. Ackerman D
17. Blakely T
18. Bogovic J
19. Dolafi T
20. Kainmueller D
21. Kawase T
22. Khairy KA
23. Leavitt L
24. Li PH
25. Lindsey L
26. Neubarth N
27. Olbris DJ
28. Otsuna H
29. Trautman ET
30. Ito M
31. Bates AS
32. Goldammer J
33. Wolff T
34. Svirskas R
35. Schlegel P
36. Neace E
37. Knecht CJ
38. Alvarado CX
39. Bailey DA
40. Ballinger S
41. Borycz JA
42. Canino BS
43. Cheatham N
44. Cook M
45. Dreher M
46. Duclos O
47. Eubanks B
48. Fairbanks K
49. Finley S
50. Forknall N
51. Francis A
52. Hopkins GP
53. Joyce EM
54. Kim S
55. Kirk NA
56. Kovalyak J
57. Lauchie SA
58. Lohff A
59. Maldonado C
60. Manley EA
61. McLin S
62. Mooney C
63. Ndama M
64. Ogundeyi O
65. Okeoma N
66. Ordish C
67. Padilla N
68. Patrick CM
69. Paterson T
70. Phillips EE
71. Phillips EM
72. Rampally N
73. Ribeiro C
74. Robertson MK
75. Rymer JT
76. Ryan SM
77. Sammons M
78. Scott AK
79. Scott AL
80. Shinomiya A
81. Smith C
82. Smith K
83. Smith NL
84. Sobeski MA
85. Suleiman A
86. Swift J
87. Takemura S
88. Talebi I
89. Tarnogorska D
90. Tenshaw E
91. Tokhi T
92. Walsh JJ
93. Yang T
94. Horne JA
95. Li F
96. Parekh R
97. Rivlin PK
98. Jayaraman V
99. Costa M
100. Jefferis GS
101. Ito K
102. Saalfeld S
103. George R
104. Meinertzhagen IA
105. Rubin GM
106. Hess HF
107. Jain V
108. Plaza SM
(2020) figshare
A connectome and analysis of the adult Drosophila central brain.

https://doi.org/10.25378/janelia.12818645.v1

References

(2015) Insect olfaction and the evolution of receptor tuning
Frontiers in Ecology and Evolution 3:53.

https://doi.org/10.3389/fevo.2015.00053
- Google Scholar
1. Badel L
2. Ohta K
3. Tsuchimoto Y
4. Kazama H
(2016) Decoding of context-dependent olfactory behavior in Drosophila
Neuron 91:155–167.

https://doi.org/10.1016/j.neuron.2016.05.022
- PubMed
- Google Scholar
1. Bak JH
2. Jang SJ
3. Hyeon C
(2018) Implications for human odor sensing revealed from the statistics of odorant-receptor interactions
PLOS Computational Biology 14:e1006175.

https://doi.org/10.1371/journal.pcbi.1006175
- PubMed
- Google Scholar
1. Baker FB
(1974) Stability of two hierarchical grouping techniques case i: sensitivity to data errors
Journal of the American Statistical Association 69:440–445.

https://doi.org/10.1080/01621459.1974.10482971
- Google Scholar
(2019) Neuronal cell types in the fly: single-cell anatomy meets single-cell genomics
Current Opinion in Neurobiology 56:125–134.

https://doi.org/10.1016/j.conb.2018.12.012
- PubMed
- Google Scholar
1. Bates AS
2. Schlegel P
3. Roberts RJV
4. Drummond N
5. Tamimi IFM
6. Turnbull R
7. Zhao X
8. Marin EC
9. Popovici PD
10. Dhawan S
11. Jamasb A
12. Javier A
13. Serratosa Capdevila L
14. Li F
15. Rubin GM
16. Waddell S
17. Bock DD
18. Costa M
19. Jefferis G
(2020) Complete connectomic reconstruction of olfactory projection neurons in the fly brain
Current Biology 30:3183–3199.

https://doi.org/10.1016/j.cub.2020.06.042
- PubMed
- Google Scholar
1. Berck ME
2. Khandelwal A
3. Claus L
4. Hernandez-Nunez L
5. Si G
6. Tabone CJ
7. Li F
8. Truman JW
9. Fetter RD
10. Louis M
11. Samuel AD
12. Cardona A
(2016) The wiring diagram of a glomerular olfactory system
eLife 5:e14859.

https://doi.org/10.7554/eLife.14859
- PubMed
- Google Scholar
(2014) Pyramidal cells make specific connections onto smooth (gabaergic) neurons in mouse visual cortex
PLOS Biology 12:e1001932.

https://doi.org/10.1371/journal.pbio.1001932
- PubMed
- Google Scholar
Conference
(2018) A deep learning approach for the classification of neuronal cell types
Annual International Conference of the IEEE Engineering in Medicine and Biology Society. IEEE Engineering in Medicine and Biology Society. Annual International Conference. pp. 999–1002.

https://doi.org/10.1109/EMBC.2018.8512498
- Google Scholar
Book
1. Cajal SR
(1911) Histologie du système nerveux de l’homme & des vertébrés: cervelet, cerveau moyen, rétine, couche optique, corps strié, écorce cérébrale générale & régionale, grand sympathique
In: Cajal SR, editors. BHL. Springer. pp. 63–90.

https://doi.org/10.5962/bhl.title.48637
- Google Scholar
1. Caron SJC
2. Ruta V
3. Abbott LF
4. Axel R
(2013) Random convergence of olfactory inputs in the Drosophila mushroom body
Nature 497:113–117.

https://doi.org/10.1038/nature12063
- PubMed
- Google Scholar
1. Choi K
2. Kim WK
3. Hyeon C
(2022) Polymer physics-based classification of neurons
Neuroinformatics 1:e3.

https://doi.org/10.1007/s12021-022-09605-3
- Google Scholar
Preprint
1. Clements J
2. Dolafi T
3. Umayam L
4. Neubarth NL
5. Berg S
6. Scheffer LK
7. Plaza SM
(2020) Neu Print: Analysis Tools for EM Connectomics
bioRxiv.

https://doi.org/10.1101/2020.01.16.909465
- Google Scholar
(2016) NBLAST: rapid, sensitive comparison of neuronal structure and construction of neuron family databases
Neuron 91:293–311.

https://doi.org/10.1016/j.neuron.2016.06.012
- PubMed
- Google Scholar
(2005) Molecular, anatomical, and functional organization of the Drosophila olfactory system
Current Biology 15:1535–1547.

https://doi.org/10.1016/j.cub.2005.07.034
- PubMed
- Google Scholar
1. Das Chakraborty S
2. Sachse S
(2021) Olfactory processing in the lateral horn of Drosophila
Cell and Tissue Research 383:113–123.

https://doi.org/10.1007/s00441-020-03392-6
- PubMed
- Google Scholar
1. Dolan M-J
2. Frechter S
3. Bates AS
4. Dan C
5. Huoviala P
6. Roberts RJ
7. Schlegel P
8. Dhawan S
9. Tabano R
10. Dionne H
11. Christoforou C
12. Close K
13. Sutcliffe B
14. Giuliani B
15. Li F
16. Costa M
17. Ihrke G
18. Meissner GW
19. Bock DD
20. Aso Y
21. Rubin GM
22. Jefferis GS
(2019) Neurogenetic dissection of the Drosophila lateral horn reveals major outputs, diverse behavioural functions, and interactions with the mushroom body
eLife 8:e43079.

https://doi.org/10.7554/eLife.43079
- PubMed
- Google Scholar
1. Dorkenwald S
2. McKellar CE
3. Macrina T
4. Kemnitz N
5. Lee K
6. Lu R
7. Wu J
8. Popovych S
9. Mitchell E
10. Nehoran B
11. Jia Z
12. Bae JA
13. Mu S
14. Ih D
15. Castro M
16. Ogedengbe O
17. Halageri A
18. Kuehner K
19. Sterling AR
20. Ashwood Z
21. Zung J
22. Brittain D
23. Collman F
24. Schneider-Mizell C
25. Jordan C
26. Silversmith W
27. Baker C
28. Deutsch D
29. Encarnacion-Rivera L
30. Kumar S
31. Burke A
32. Bland D
33. Gager J
34. Hebditch J
35. Koolman S
36. Moore M
37. Morejohn S
38. Silverman B
39. Willie K
40. Willie R
41. Yu S-C
42. Murthy M
43. Seung HS
(2022) FlyWire: online community for whole-brain connectomics
Nature Methods 19:119–128.

https://doi.org/10.1038/s41592-021-01330-0
- PubMed
- Google Scholar
(1989) A canonical microcircuit for neocortex
Neural Computation 1:480–488.

https://doi.org/10.1162/neco.1989.1.4.480
- Google Scholar
1. Dweck HKM
2. Ebrahim SAM
3. Thoma M
4. Mohamed AAM
5. Keesey IW
6. Trona F
7. Lavista-Llanos S
8. Svatoš A
9. Sachse S
10. Knaden M
11. Hansson BS
(2015) Pheromones mediating copulation and attraction in Drosophila
PNAS 112:E2829–E2835.

https://doi.org/10.1073/pnas.1504527112
- PubMed
- Google Scholar
1. Eichler K
2. Li F
3. Litwin-Kumar A
4. Park Y
5. Andrade I
6. Schneider-Mizell CM
7. Saumweber T
8. Huser A
9. Eschbach C
10. Gerber B
11. Fetter RD
12. Truman JW
13. Priebe CE
14. Abbott LF
15. Thum AS
16. Zlatic M
17. Cardona A
(2017) The complete connectome of a learning and memory centre in an insect brain
Nature 548:175–182.

https://doi.org/10.1038/nature23455
- PubMed
- Google Scholar
1. Eschbach C
2. Fushiki A
3. Winding M
4. Schneider-Mizell CM
5. Shao M
6. Arruda R
7. Eichler K
8. Valdes-Aleman J
9. Ohyama T
10. Thum AS
11. Gerber B
12. Fetter RD
13. Truman JW
14. Litwin-Kumar A
15. Cardona A
16. Zlatic M
(2020) Recurrent architecture for adaptive regulation of learning in the insect brain
Nature Neuroscience 23:544–555.

https://doi.org/10.1038/s41593-020-0607-9
- PubMed
- Google Scholar
1. Fişek M
2. Wilson RI
(2014) Stereotyped connectivity and computations in higher-order olfactory neurons
Nature Neuroscience 17:280–288.

https://doi.org/10.1038/nn.3613
- PubMed
- Google Scholar
1. Fishilevich E
2. Vosshall LB
(2005) Genetic and functional subdivision of the Drosophila antennal lobe
Current Biology 15:1548–1553.

https://doi.org/10.1016/j.cub.2005.07.066
- PubMed
- Google Scholar
1. Frechter S
2. Bates AS
3. Tootoonian S
4. Dolan MJ
5. Manton J
6. Jamasb AR
7. Kohl J
8. Bock D
9. Jefferis G
(2019) Functional and anatomical specificity in a higher olfactory centre
eLife 8:e44590.

https://doi.org/10.7554/eLife.44590
- PubMed
- Google Scholar
1. Galili T
(2015) Dendextend: an R package for visualizing, adjusting and comparing trees of hierarchical clustering
Bioinformatics 31:3718–3720.

https://doi.org/10.1093/bioinformatics/btv428
- PubMed
- Google Scholar
1. Galizia CG
2. Sachse S
(2010) Odor coding in insects
The Neurobiology of Olfaction 35:70.

https://doi.org/10.1201/9781420071993
- Google Scholar
1. Galizia CG
(2014) Olfactory coding in the insect brain: data and conjectures
The European Journal of Neuroscience 39:1784–1795.

https://doi.org/10.1111/ejn.12558
- PubMed
- Google Scholar
1. Gao Q
2. Yuan B
3. Chess A
(2000) Convergent projections of Drosophila olfactory neurons to specific glomeruli in the antennal lobe
Nature Neuroscience 3:780–785.

https://doi.org/10.1038/77680
- PubMed
- Google Scholar
1. Gillette TA
2. Ascoli GA
(2015) Topological characterization of neuronal arbor morphology via sequence representation: I--motif analysis
BMC Bioinformatics 16:1–15.

https://doi.org/10.1186/s12859-015-0604-2
- PubMed
- Google Scholar
1. Gouwens NW
2. Sorensen SA
3. Berg J
4. Lee C
5. Jarsky T
6. Ting J
7. Sunkin SM
8. Feng D
9. Anastassiou CA
10. Barkan E
11. Bickley K
12. Blesie N
13. Braun T
14. Brouner K
15. Budzillo A
16. Caldejon S
17. Casper T
18. Castelli D
19. Chong P
20. Crichton K
21. Cuhaciyan C
22. Daigle TL
23. Dalley R
24. Dee N
25. Desta T
26. Ding S-L
27. Dingman S
28. Doperalski A
29. Dotson N
30. Egdorf T
31. Fisher M
32. de Frates RA
33. Garren E
34. Garwood M
35. Gary A
36. Gaudreault N
37. Godfrey K
38. Gorham M
39. Gu H
40. Habel C
41. Hadley K
42. Harrington J
43. Harris JA
44. Henry A
45. Hill D
46. Josephsen S
47. Kebede S
48. Kim L
49. Kroll M
50. Lee B
51. Lemon T
52. Link KE
53. Liu X
54. Long B
55. Mann R
56. McGraw M
57. Mihalas S
58. Mukora A
59. Murphy GJ
60. Ng L
61. Ngo K
62. Nguyen TN
63. Nicovich PR
64. Oldre A
65. Park D
66. Parry S
67. Perkins J
68. Potekhina L
69. Reid D
70. Robertson M
71. Sandman D
72. Schroedter M
73. Slaughterbeck C
74. Soler-Llavina G
75. Sulc J
76. Szafer A
77. Tasic B
78. Taskin N
79. Teeter C
80. Thatra N
81. Tung H
82. Wakeman W
83. Williams G
84. Young R
85. Zhou Z
86. Farrell C
87. Peng H
88. Hawrylycz MJ
89. Lein E
90. Ng L
91. Arkhipov A
92. Bernard A
93. Phillips JW
94. Zeng H
95. Koch C
(2019) Classification of electrophysiological and morphological neuron types in the mouse visual cortex
Nature Neuroscience 22:1182–1195.

https://doi.org/10.1038/s41593-019-0417-0
- PubMed
- Google Scholar
(2016) Elucidating the neuronal architecture of olfactory glomeruli in the Drosophila antennal lobe
Cell Reports 16:3401–3413.

https://doi.org/10.1016/j.celrep.2016.08.063
- PubMed
- Google Scholar
1. Grosjean Y
2. Rytz R
3. Farine J-P
4. Abuin L
5. Cortot J
6. Jefferis GSXE
7. Benton R
(2011) An olfactory receptor for food-derived odours promotes male courtship in Drosophila
Nature 478:236–240.

https://doi.org/10.1038/nature10428
- PubMed
- Google Scholar
1. Hallem EA
2. Carlson JR
(2004a) The odor coding system of Drosophila
TRENDS in Genetics 20:453–459.

https://doi.org/10.1016/j.tig.2004.06.015
- PubMed
- Google Scholar
(2004b) The molecular basis of odor coding in the Drosophila antenna
Cell 117:965–979.

https://doi.org/10.1016/j.cell.2004.05.012
- PubMed
- Google Scholar
1. Hope ACA
(1968) A simplified monte carlo significance test procedure
Journal of the Royal Statistical Society 30:582–598.

https://doi.org/10.1111/j.2517-6161.1968.tb00759.x
- Google Scholar
1. Howard JD
2. Gottfried JA
(2014) Configural and elemental coding of natural odor mixture components in the human brain
Neuron 84:857–869.

https://doi.org/10.1016/j.neuron.2014.10.012
- PubMed
- Google Scholar
Preprint
1. Huoviala P
2. Dolan MJ
3. Love F
4. Myers P
5. Frechter S
6. Namiki S
7. Pettersson L
8. Roberts RJ
9. Turnbull R
10. Mitrevica Z
(2020) Neural circuit basis of aversive odour processing in Drosophila from sensory input to descending output
bioRxiv.

https://doi.org/10.1101/394403
- Google Scholar
(2014) A systematic Nomenclature for the insect brain
Neuron 81:755–765.

https://doi.org/10.1016/j.neuron.2013.12.017
- PubMed
- Google Scholar
1. Jeanne JM
2. Wilson RI
(2015) Convergence, divergence, and reconvergence in a feedforward network improves neural speed and accuracy
Neuron 88:1014–1026.

https://doi.org/10.1016/j.neuron.2015.10.018
- PubMed
- Google Scholar
(2018) The organization of projections from olfactory glomeruli onto higher-order neurons
Neuron 98:1198–1213.

https://doi.org/10.1016/j.neuron.2018.05.011
- PubMed
- Google Scholar
1. Jefferis GSXE
2. Potter CJ
3. Chan AM
4. Marin EC
5. Rohlfing T
6. Maurer CR Jr
7. Luo L
(2007) Comprehensive maps of Drosophila higher olfactory centers: spatially segregated fruit and pheromone representation
Cell 128:1187–1203.

https://doi.org/10.1016/j.cell.2007.01.040
- PubMed
- Google Scholar
1. Jenett A
2. Rubin GM
3. Ngo T-TB
4. Shepherd D
5. Murphy C
6. Dionne H
7. Pfeiffer BD
8. Cavallaro A
9. Hall D
10. Jeter J
11. Iyer N
12. Fetter D
13. Hausenfluck JH
14. Peng H
15. Trautman ET
16. Svirskas RR
17. Myers EW
18. Iwinski ZR
19. Aso Y
20. DePasquale GM
21. Enos A
22. Hulamm P
23. Lam SCB
24. Li H-H
25. Laverty TR
26. Long F
27. Qu L
28. Murphy SD
29. Rokicki K
30. Safford T
31. Shaw K
32. Simpson JH
33. Sowell A
34. Tae S
35. Yu Y
36. Zugates CT
(2012) A gal4-driver line resource for Drosophila neurobiology
Cell Reports 2:991–1001.

https://doi.org/10.1016/j.celrep.2012.09.011
- PubMed
- Google Scholar
1. Kanari L
2. Dłotko P
3. Scolamiero M
4. Levi R
5. Shillcock J
6. Hess K
7. Markram H
(2018) A topological representation of branching neuronal morphologies
Neuroinformatics 16:3–13.

https://doi.org/10.1007/s12021-017-9341-1
- PubMed
- Google Scholar
(2013) A bidirectional circuit switch reroutes pheromone signals in male and female brains
Cell 155:1610–1623.

https://doi.org/10.1016/j.cell.2013.11.025
- PubMed
- Google Scholar
(2015) Pheromone processing in Drosophila
Current Opinion in Neurobiology 34:149–157.

https://doi.org/10.1016/j.conb.2015.06.009
- PubMed
- Google Scholar
(2008) Defining clusters from a hierarchical cluster tree: the dynamic tree cut package for R
Bioinformatics 24:719–720.

https://doi.org/10.1093/bioinformatics/btm563
- PubMed
- Google Scholar
(2020) A systematic evaluation of interneuron morphology representations for cell type discrimination
Neuroinformatics 18:591–609.

https://doi.org/10.1007/s12021-020-09461-z
- PubMed
- Google Scholar
1. Li Y
2. Wang D
3. Ascoli GA
4. Mitra P
5. Wang Y
(2017) Metrics for comparing neuronal tree shapes based on persistent homology
PLOS ONE 12:e0182184.

https://doi.org/10.1371/journal.pone.0182184
- PubMed
- Google Scholar
1. Li F
2. Lindsey JW
3. Marin EC
4. Otto N
5. Dreher M
6. Dempsey G
7. Stark I
8. Bates AS
9. Pleijzier MW
10. Schlegel P
11. Nern A
12. Takemura S-Y
13. Eckstein N
14. Yang T
15. Francis A
16. Braun A
17. Parekh R
18. Costa M
19. Scheffer LK
20. Aso Y
21. Jefferis GS
22. Abbott LF
23. Litwin-Kumar A
24. Waddell S
25. Rubin GM
(2020) The connectome of the adult Drosophila mushroom body provides insights into function
eLife 9:e62576.

https://doi.org/10.7554/eLife.62576
- PubMed
- Google Scholar
1. Liang L
2. Li Y
3. Potter CJ
4. Yizhar O
5. Deisseroth K
6. Tsien RW
7. Luo L
(2013) Gabaergic projection neurons route selective olfactory inputs to specific higher-order neurons
Neuron 79:917–931.

https://doi.org/10.1016/j.neuron.2013.06.014
- PubMed
- Google Scholar
1. Lin HH
2. Lai JSY
3. Chin AL
4. Chen YC
5. Chiang AS
(2007) A map of olfactory representation in the Drosophila mushroom body
Cell 128:1205–1217.

https://doi.org/10.1016/j.cell.2007.03.006
- PubMed
- Google Scholar
1. Lu Y
2. Carin L
3. Coifman R
4. Shain W
5. Roysam B
(2015) Quantitative arbor analytics: unsupervised harmonic co-clustering of populations of brain cell arbors based on l-measure
Neuroinformatics 13:47–63.

https://doi.org/10.1007/s12021-014-9237-2
- PubMed
- Google Scholar
1. Mansourian S
2. Stensmyr MC
(2015) The chemical ecology of the fly
Current Opinion in Neurobiology 34:95–102.

https://doi.org/10.1016/j.conb.2015.02.006
- PubMed
- Google Scholar
(2008) Principles of glomerular organization in the human olfactory bulb -- implications for odor processing
PLOS ONE 3:e2640.

https://doi.org/10.1371/journal.pone.0002640
- PubMed
- Google Scholar
1. Marin EC
2. Jefferis GSXE
3. Komiyama T
4. Zhu H
5. Luo L
(2002) Representation of the glomerular olfactory map in the Drosophila brain
Cell 109:243–255.

https://doi.org/10.1016/s0092-8674(02)00700-6
- PubMed
- Google Scholar
1. Marin EC
2. Büld L
3. Theiss M
4. Sarkissian T
5. Roberts RJV
6. Turnbull R
7. Tamimi IFM
8. Pleijzier MW
9. Laursen WJ
10. Drummond N
11. Schlegel P
12. Bates AS
13. Li F
14. Landgraf M
15. Costa M
16. Bock DD
17. Garrity PA
18. Jefferis GSXE
(2020) Connectomics analysis reveals first-, second-, and third-order thermosensory and hygrosensory neurons in the adult Drosophila brain
Current Biology 30:3167–3182.

https://doi.org/10.1016/j.cub.2020.06.028
- PubMed
- Google Scholar
(2018) Towards a supervised classification of neocortical interneuron morphologies
BMC Bioinformatics 19:1–22.

https://doi.org/10.1186/s12859-018-2470-1
- PubMed
- Google Scholar
1. Min S
2. Ai M
3. Shin SA
4. Suh GSB
(2013) Dedicated olfactory neurons mediating attraction behavior to ammonia and amines in Drosophila
PNAS 110:E1321–E1329.

https://doi.org/10.1073/pnas.1215680110
- PubMed
- Google Scholar
(2019) Odor mixtures of opposing valence unveil inter-glomerular crosstalk in the Drosophila antennal lobe
Nature Communications 10:1–17.

https://doi.org/10.1038/s41467-019-09069-1
- PubMed
- Google Scholar
1. Mombaerts P
2. Wang F
3. Dulac C
4. Chao SK
5. Nemes A
6. Mendelsohn M
7. Edmondson J
8. Axel R
(1996) Visualizing an olfactory sensory map
Cell 87:675–686.

https://doi.org/10.1016/s0092-8674(00)81387-2
- PubMed
- Google Scholar
1. Münch D
2. Galizia CG
(2016) Door 2.0 -- comprehensive mapping of Drosophila melanogaster odorant responses
Scientific Reports 6:1–14.

https://doi.org/10.1038/srep21841
- PubMed
- Google Scholar
(2010) Divisive normalization in olfactory population codes
Neuron 66:287–299.

https://doi.org/10.1016/j.neuron.2010.04.009
- PubMed
- Google Scholar
1. Ruta V
2. Datta SR
3. Vasconcelos ML
4. Freeland J
5. Looger LL
6. Axel R
(2010) A dimorphic pheromone circuit in Drosophila from sensory input to descending output
Nature 468:686–690.

https://doi.org/10.1038/nature09554
- PubMed
- Google Scholar
(2009) CATMAID: collaborative annotation toolkit for massive amounts of image data
Bioinformatics 25:1984–1986.

https://doi.org/10.1093/bioinformatics/btp266
- PubMed
- Google Scholar
1. Scheffer LK
2. Xu CS
3. Januszewski M
4. Lu Z
5. Takemura S-Y
6. Hayworth KJ
7. Huang GB
8. Shinomiya K
9. Maitlin-Shepard J
10. Berg S
11. Clements J
12. Hubbard PM
13. Katz WT
14. Umayam L
15. Zhao T
16. Ackerman D
17. Blakely T
18. Bogovic J
19. Dolafi T
20. Kainmueller D
21. Kawase T
22. Khairy KA
23. Leavitt L
24. Li PH
25. Lindsey L
26. Neubarth N
27. Olbris DJ
28. Otsuna H
29. Trautman ET
30. Ito M
31. Bates AS
32. Goldammer J
33. Wolff T
34. Svirskas R
35. Schlegel P
36. Neace E
37. Knecht CJ
38. Alvarado CX
39. Bailey DA
40. Ballinger S
41. Borycz JA
42. Canino BS
43. Cheatham N
44. Cook M
45. Dreher M
46. Duclos O
47. Eubanks B
48. Fairbanks K
49. Finley S
50. Forknall N
51. Francis A
52. Hopkins GP
53. Joyce EM
54. Kim S
55. Kirk NA
56. Kovalyak J
57. Lauchie SA
58. Lohff A
59. Maldonado C
60. Manley EA
61. McLin S
62. Mooney C
63. Ndama M
64. Ogundeyi O
65. Okeoma N
66. Ordish C
67. Padilla N
68. Patrick CM
69. Paterson T
70. Phillips EE
71. Phillips EM
72. Rampally N
73. Ribeiro C
74. Robertson MK
75. Rymer JT
76. Ryan SM
77. Sammons M
78. Scott AK
79. Scott AL
80. Shinomiya A
81. Smith C
82. Smith K
83. Smith NL
84. Sobeski MA
85. Suleiman A
86. Swift J
87. Takemura S
88. Talebi I
89. Tarnogorska D
90. Tenshaw E
91. Tokhi T
92. Walsh JJ
93. Yang T
94. Horne JA
95. Li F
96. Parekh R
97. Rivlin PK
98. Jayaraman V
99. Costa M
100. Jefferis GS
101. Ito K
102. Saalfeld S
103. George R
104. Meinertzhagen IA
105. Rubin GM
106. Hess HF
107. Jain V
108. Plaza SM
(2020) A connectome and analysis of the adult Drosophila central brain
eLife 9:e57443.

https://doi.org/10.7554/eLife.57443
- PubMed
- Google Scholar
1. Schlegel P
2. Bates AS
3. Stürner T
4. Jagannathan SR
5. Drummond N
6. Hsu J
7. Serratosa Capdevila L
8. Javier A
9. Marin EC
10. Barth-Maron A
11. Tamimi IF
12. Li F
13. Rubin GM
14. Plaza SM
15. Costa M
16. Jefferis GSXE
(2021) Information flow, cell types and stereotypy in a full olfactory connectome
eLife 10:e66018.

https://doi.org/10.7554/eLife.66018
- PubMed
- Google Scholar
(2017) The role of the Drosophila lateral horn in olfactory information processing and behavioral response
Journal of Insect Physiology 98:29–37.

https://doi.org/10.1016/j.jinsphys.2016.11.007
- PubMed
- Google Scholar
(2008) L-measure: a web-accessible tool for the analysis, comparison and search of digital reconstructions of neuronal morphologies
Nature Protocols 3:866–876.

https://doi.org/10.1038/nprot.2008.51
- PubMed
- Google Scholar
1. Seki Y.
2. Rybak J
3. Wicher D
4. Sachse S
5. Hansson BS
(2010) Physiological and morphological characterization of local interneurons in the Drosophila antennal lobe
Journal of Neurophysiology 104:1007–1019.

https://doi.org/10.1152/jn.00249.2010
- PubMed
- Google Scholar
1. Seki Y
2. Dweck HKM
3. Rybak J
4. Wicher D
5. Sachse S
6. Hansson BS
(2017) Olfactory coding from the periphery to higher brain centers in the Drosophila brain
BMC Biology 15:1–20.

https://doi.org/10.1186/s12915-017-0389-z
- PubMed
- Google Scholar
1. Semmelhack JL
2. Wang JW
(2009) Select Drosophila glomeruli mediate innate olfactory attraction and aversion
Nature 459:218–223.

https://doi.org/10.1038/nature07983
- PubMed
- Google Scholar
1. Shimizu K
2. Stopfer M
(2017) A population of projection neurons that inhibits the lateral horn but excites the antennal lobe through chemical synapses in Drosophila
Frontiers in Neural Circuits 11:30.

https://doi.org/10.3389/fncir.2017.00030
- PubMed
- Google Scholar
1. Stensmyr MC
2. Dweck HKM
3. Farhan A
4. Ibba I
5. Strutz A
6. Mukunda L
7. Linz J
8. Grabe V
9. Steck K
10. Lavista-Llanos S
11. Wicher D
12. Sachse S
13. Knaden M
14. Becher PG
15. Seki Y
16. Hansson BS
(2012) A conserved dedicated olfactory circuit for detecting harmful microbes in Drosophila
Cell 151:1345–1357.

https://doi.org/10.1016/j.cell.2012.09.046
- PubMed
- Google Scholar
1. Stevens CF
(2015) What the fly’s nose tells the fly’s brain
PNAS 112:9460–9465.

https://doi.org/10.1073/pnas.1510103112
- PubMed
- Google Scholar
1. Strutz A
2. Soelter J
3. Baschwitz A
4. Farhan A
5. Grabe V
6. Rybak J
7. Knaden M
8. Schmuker M
9. Hansson BS
10. Sachse S
(2014) Decoding odor quality and intensity in the Drosophila brain
eLife 3:e04147.

https://doi.org/10.7554/eLife.04147
- PubMed
- Google Scholar
(2012) Plasticity of recurrent inhibition in the Drosophila antennal lobe
The Journal of Neuroscience 32:7225–7231.

https://doi.org/10.1523/JNEUROSCI.1099-12.2012
- PubMed
- Google Scholar
(2017) Wiring variations that enable and constrain neural computation in a sensory microcircuit
eLife 6:e24838.

https://doi.org/10.7554/eLife.24838
- PubMed
- Google Scholar
(2008) Olfactory representations by Drosophila mushroom body neurons
Journal of Neurophysiology 99:734–746.

https://doi.org/10.1152/jn.01283.2007
- PubMed
- Google Scholar
1. Uylings HBM
2. van Pelt J
(2002) Measures for quantifying dendritic arborizations
Network 13:397–414.

https://doi.org/10.1088/0954-898X_13_3_309
- PubMed
- Google Scholar
1. Vasques X
2. Vanel L
3. Villette G
4. Cif L
(2016) Morphological neuron classification using machine learning
Frontiers in Neuroanatomy 10:102.

https://doi.org/10.3389/fnana.2016.00102
- PubMed
- Google Scholar
(2000) An olfactory sensory map in the fly brain
Cell 102:147–159.

https://doi.org/10.1016/s0092-8674(00)00021-0
- PubMed
- Google Scholar
(2004) Transformation of olfactory representations in the Drosophila antennal lobe
Science 303:366–370.

https://doi.org/10.1126/science.1090782
- PubMed
- Google Scholar
1. Wong AM
2. Wang JW
3. Axel R
(2002) Spatial representation of the glomerular map in the Drosophila protocerebrum
Cell 109:229–241.

https://doi.org/10.1016/s0092-8674(02)00707-9
- PubMed
- Google Scholar
1. Zhang T
2. Zeng Y
3. Zhang Y
4. Zhang X
5. Shi M
6. Tang L
7. Zhang D
8. Xu B
(2021) Neuron type classification in rat brain based on integrative convolutional and tree-based recurrent neural networks
Scientific Reports 11:1–14.

https://doi.org/10.1038/s41598-021-86780-4
- PubMed
- Google Scholar
1. Zheng Z
2. Lauritzen JS
3. Perlman E
4. Robinson CG
5. Nichols M
6. Milkie D
7. Torrens O
8. Price J
9. Fisher CB
10. Sharifi N
11. Calle-Schuler SA
12. Kmecova L
13. Ali IJ
14. Karsh B
15. Trautman ET
16. Bogovic JA
17. Hanslovsky P
18. Jefferis GSXE
19. Kazhdan M
20. Khairy K
21. Saalfeld S
22. Fetter RD
23. Bock DD
(2018) A complete electron microscopy volume of the brain of adult Drosophila melanogaster
Cell 174:730–743.

https://doi.org/10.1016/j.cell.2018.06.019
- PubMed
- Google Scholar
Preprint
1. Zheng Z
2. Li F
3. Fisher C
4. Ali IJ
5. Sharifi N
6. Calle-Schuler S
7. Hsu J
8. Masoodpanah N
9. Kmecova L
10. Kazimiers T
11. Perlman E
12. Nichols M
13. Li PH
14. Jain V
15. Bock DD
(2020) Structured Sampling of Olfactory Input by the Fly Mushroom Body
bioRxiv.

https://doi.org/10.1101/2020.04.17.047167
- Google Scholar

Article and author information

Author details

Kiri Choi

School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea

Contribution
Conceptualization, Data curation, Software, Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing – review and editing

For correspondence
ckiri0315@kias.re.kr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-0156-8410
Won Kyu Kim

School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea

Contribution
Conceptualization, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing – review and editing

For correspondence
wonkyukim@kias.re.kr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-6286-0925
Changbong Hyeon

School of Computational Sciences, Korea Institute for Advanced Study, Seoul, Republic of Korea

Contribution
Conceptualization, Resources, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing – review and editing

For correspondence
hyeoncb@kias.re.kr

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4844-7237

Funding

KIAS individual grant (CG077001)

Kiri Choi

KIAS individual grant (CG076002)

Won Kyu Kim

KIAS individual grant (CG035003)

Changbong Hyeon

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dr. Ji Hyun Bak for helpful discussions. This study was supported by KIAS Individual Grants CG077001 (KC), CG076002 (WKK), and CG035003 (CH). We thank the Center for Advanced Computation in KIAS for providing the computing resources.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,532

views
362

downloads
8

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Kiri Choi
Won Kyu Kim
Changbong Hyeon

(2022)

Olfactory responses of Drosophila are encoded in the organization of projection neurons

eLife 11:e77748.

https://doi.org/10.7554/eLife.77748

Categories and tags

Research organism

D. melanogaster

Share this article

Cite this article

A schematic of the Drosophila olfactory system.

The three matrices representing the pairwise distances dα⁢β in units of μ⁢m between individual uPN in AL, MB calyx, and LH.

The dα⁢β-based clustering on uPNs based on the FAFB dataset in MB calyx resulting in 10 clusters.

The dα⁢β-based clustering on uPNs based on the FAFB dataset in LH resulting in 11 clusters.

Organization of homotypic uPNs in the three neuropils.

Degree of overlapping between inter-homotypic uPNs, λX (X=VM4,VC5,…,VP5).

Scatter plots depicting the relationships between λX s at two different neuropils calculated from the uPNs based on the FAFB dataset; (A) AL versus MB calyx, (B) AL versus LH, and (C) LH versus MB calyx.

A schematic illustrating the connectivity between homotypes (X=A,B,…,E) and third-order neurons (i=1,2,…,7).

Bar graphs depicting the number of KCs/LHNs that synapse with a specific homotype X (NX,sp, blue) and the percentage of KCs/LHNs that synapse with a specific homotype X (fX=NX,sp/NX,tot, red) at (A) PN-KC and (B) PN-LHN interfaces, with the synaptic weight of N=3 used as the threshold.

Common synapse matrices (A) SKC and (B) SLHN, each of which represents the extent of signal integration from homotypic uPNs to KCs and LHNs.

Tanglegrams comparing the tree structures generated from the inter-PN distances-based (left) and the connectivity-based clustering (right) (A) between uPNs and KCs, and (B) between uPNs and LHNs.

A diagram depicting the neuropil segmentation process.

The optimal number of clusters of uPNs in the FAFB dataset determined by employing the dynamic hybrid cut tree method, elbow method, gap statistics, and maximum average silhouette coefficient.

Pearson’s χ2 tests of independence of variables in the FAFB dataset.

Pearson’s χ2 tests of independence of variables in the hemibrain dataset.

Mutual information (observed mutual information (top), randomly sampled mutual information (bottom) in each cell) from the association study using the FAFB dataset.

Mutual information (observed mutual information (top), randomly sampled mutual information (bottom) in each cell) from the association study using the hemibrain dataset.

Statistics of homotypes composed of a single uPN (or multiple uPNs) in the FAFB dataset and the corresponding putative valence.

Statistics of homotypes composed of a single uPN (or multiple uPNs) in the hemibrain dataset and the corresponding putative valence.

Pearson’s χ2 tests of independence of variables on the connectivity-based clustering results.

Author details

Kiri Choi

Contribution

For correspondence

Competing interests

Won Kyu Kim

Contribution

For correspondence

Competing interests

Changbong Hyeon

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organism

Further reading

The three matrices representing the pairwise distances $d_{α β}$ in units of $μ m$ between individual uPN in AL, MB calyx, and LH.

The $d_{α β}$ -based clustering on uPNs based on the FAFB dataset in MB calyx resulting in 10 clusters.

The $d_{α β}$ -based clustering on uPNs based on the FAFB dataset in LH resulting in 11 clusters.

Degree of overlapping between inter-homotypic uPNs, $λ_{X}$ ( $X = VM4, VC5, \dots, VP5$ ).

Scatter plots depicting the relationships between $λ_{X}$ s at two different neuropils calculated from the uPNs based on the FAFB dataset; (A) AL versus MB calyx, (B) AL versus LH, and (C) LH versus MB calyx.

A schematic illustrating the connectivity between homotypes ( $X = A, B, \dots, E$ ) and third-order neurons ( $i = 1, 2, \dots, 7$ ).

Common synapse matrices (A) $S^{K C}$ and (B) $S^{L H N}$ , each of which represents the extent of signal integration from homotypic uPNs to KCs and LHNs.

Pearson’s $χ^{2}$ tests of independence of variables in the FAFB dataset.

Pearson’s $χ^{2}$ tests of independence of variables in the hemibrain dataset.

Pearson’s $χ^{2}$ tests of independence of variables on the connectivity-based clustering results.