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

  1. Kiri Choi  Is a corresponding author
  2. Won Kyu Kim  Is a corresponding author
  3. Changbong Hyeon  Is a corresponding author
  1. School of Computational Sciences, Korea Institute for Advanced Study, Republic of Korea
14 figures, 8 tables and 2 additional files

Figures

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.

Figure 2 with 3 supplements
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 2—figure supplement 1
An illustration of inter-PN distance calculation and comparison with the NBLAST distance.

(A) A schematic showing how the ‘distance’ between two neurons α and β (dαβ) is calculated. (B) The comparisons between dαβ and the normalized NBLAST distance (which is equal to 1-NBLAST score) for uPN innervation to AL, MB calyx, and LH. While the two distances are correlated, significant dispersion is also present. Unlike NBLAST score, dαβ measures the spatial proximity between two neurons α and β only, but not their morphological similarity.

Figure 2—figure supplement 2
Two 135×135 matrices representing the inter-neuronal distances in (A) MB calyx and (B) LH that are reorganized based on the clustering results.

The analysis was done on the uPNs in the FAFB dataset. The diagonal block represents each cluster (see Figure 3, Figure 4, and Figure 2—figure supplement 3 for detailed information on the clustering labels). The indices of uPNs are reorganized based on the results from the clustering analysis.

Figure 2—figure supplement 3
The dendrograms of dαβ-based clustering on uPNs innervation in (A) AL, (B) MB calyx, and (C) LH.

The analysis was done on the uPNs in the FAFB dataset. In (B) and (C), the different colored leaves correspond to a cluster generated from the tree-cutting method. A leaf represents an individual uPN and the label depicts the corresponding homotype for each uPN.

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
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 CILH=C3LH and COLH=C4,5,6,9LH. 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.

Figure 4—figure supplement 1
15 clusters from the dαβ-based clustering on the entire PN innervation in LH including those that do not innervate all three neuropils such as GABAergic mPNs.

The 15 clusters from the dαβ-based clustering on the entire uPN innervation in LH using the FAFB dataset. 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.

Figure 4—figure supplement 2
Comparison of innervation pattern of PNs in LH between the uPNs innervating all three neuropils (gray, most of which follow mALT) and those that bypass MB calyx (black, most of which follow mlALT) for the FAFB dataset.

Shown are 14 homotypes consisting of uPNs whose innervation is localized to LH.

Figure 5 with 2 supplements
Organization of homotypic uPNs in the three neuropils.

Plots of d¯intra (blue, degree of bundling), 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).

Figure 5—figure supplement 1
Comparison of the intra- (d¯intra,X) and inter-PN (d¯inter,X) distances of X-th homotype.

Comparison of the intra- (d¯intra,X, degree of bundling) and inter-PN (d¯inter,X, degree of packing) distances of X-th homotype in AL, MB calyx, and LH (from lighter to darker colors) from (A) the FAFB dataset and (B) the hemibrain dataset. The homotype label is color-coded based on the odor types obtained from the literature (Dark green: decaying fruit, lime: yeasty, green: fruity, gray: unknown, cyan: alcoholic fermentation, red: general bad/unclassified aversive, beige: plant matter, brown: animal matter, purple: pheromones, pink: hygro/thermo). Homotypes are ordered based on both the odor type and the values of λX in 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 and X denote attractive and aversive odors, respectively.

Figure 5—figure supplement 2
A plot depicting d¯intra (blue, degree of bundling), d¯inter (orange, degree of packing), and the ratio between the two distances λ (red, degree of overlapping) of 15 homotypes without (left) and with (right) 27 additional uPNs added to the FAFB dataset, which are mostly GABAergic and follow mlALT.
Figure 6 with 1 supplement
Degree of overlapping between inter-homotypic uPNs, λX (X=VM4,VC5,,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 6—figure supplement 1
Selected morphologies of uPN innervation at each neuropil for single uPN homotypes and multiple uPN homotypes.
Figure 7 with 1 supplement
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.

Figure 7—figure supplement 1
The same graph as Figure 7 based on the hemibrain dataset.

Scatter plots depicting the relationships between λX s at two different neuropils calculated from the uPNs based on the hemibrain 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.

Figure 8 with 1 supplement
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).

Figure 8—figure supplement 1
The same graph as Figure 8 based on the hemibrain dataset.

A diagram summarizing how the clusters of uPNs in MB calyx (14 clusters) and LH (13 clusters) from the hemibrain 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).

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

(A) The connectivity matrix C , where CX,i=1 when any uPNs in the X-th homotype and i-th third-order neuron synapses and CX,i=0 otherwise. (B) The number of X-th homotype-specific connections (NX,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 (NA,tot=3). Neuron 1 is not synapsing with any other homotype but A, and hence NA,sp=1; similarly, ND,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, SBC=3 in the common synapse matrix S.

Figure 9—figure supplement 1
The raw connectivity matrices.

The 58×1754 and 58×1285 synaptic connectivity matrices for the synaptic interface between PNs and KCs (CPNKC) (left) and between PNs and LHNs (CPNLHN) (right). Homotypic uPNs are merged together.

Figure 10 with 1 supplement
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.
Figure 10—figure supplement 1
Homotype specific connections with the synaptic connectivity threshold of N=8.

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 connectivity threshold of N=8.

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

The color code represents SXY, 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 (SXY=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.

Figure 11—figure supplement 1
The total number of third-order neurons connected to each homotype.

The total number of third-order neurons connected to each homotype, which corresponds to the diagonal element of common synapse matrix (Sη) in MB calyx (top) and LH (bottom). Figure 11—figure supplement 2. [The same as Figure 11 but with a higher threshold (N=8) used to identify synaptic connectivity]. Common synapse matrices (A) SKC and (B) SLHN calculated with the synaptic connectivity threshold of N=8, each of which represents the extent of signal integration from homotypic uPNs to KCs and LHNs. The black color is used when there is no third-order neuron-mediated signal integration (SXY=0) between two homotypes X and Y.

Figure 11—figure supplement 2
Common Synapse matrices (A) SKC and (B) SLHN calculated with the synaptic connectivity threshold of N=8, each of which represents the extent of signal integration from homotypic uPNs to KCs and LHNs.

The black colour is when there is no third-order neuron-mediated signal integration (SXY=0) between two homotypes X and Y.

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 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.

Author response image 1
Comparison of volumes defined in CATMAID (blue) and boundaries defined by our segmentation method (black) for MB calyx (top) and LH (bottom).

Tables

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 hybridElbowGapSilhouette
AL1914854
MB calyx101172
LH11972
Appendix 1—table 1
Pearson’s χ2 tests of independence of variables in the FAFB dataset.

CZ 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.

CLHGlomerular LabelsOdor TypeOdor Valence
CMB0.502 (1.149E-36)0.610 (1.255E-27)0.401 (3.303E-21)0.425 (2.016E-07)
CLH0.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.

CZ 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.

CLHGlomerular LabelsOdor TypeOdor Valence
CMB0.495 (6.635E-34)0.577 (3.461E-25)0.425 (9.400E-18)0.463 (6.283E-07)
CLH0.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.

CZ 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σ.

CLHGlomerular LabelsOdor TypeOdor Valence
CMB1.076 

0.397±0.045
1.892 

1.386±0.035
0.819 

0.337±0.044
0.277 

0.073±0.026
CLH2.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.

CZ 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σ.

CLHGlomerular LabelsOdor TypeOdor Valence
CMB1.371

0.710±0.048
2.244 

1.783±0.033
1.036 

0.527±0.047
0.336

0.124±0.035
CLH2.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.
AversiveAttractiveUnknownTotal
Single uPN Homotypes Count74213
Multiple uPN Homotypes Count18131344
Total25171557
Appendix 1—table 6
Statistics of homotypes composed of a single uPN (or multiple uPNs) in the hemibrain dataset and the corresponding putative valence.
AversiveAttractiveUnknownTotal
Single uPN Homotypes Count75113
Multiple uPN Homotypes Count18121545
Total25171658
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 LabelsOdor TypeOdor Valence
CPNKC0.433 (2.472E-08)0.316 (9.978E-09)0.271 (0.012)
CPNLHN0.765 (1.410E-67)0.630 (1.519E-48)0.604 (4.055E-12)

Additional files

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)

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

  1. Kiri Choi
  2. Won Kyu Kim
  3. 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