Olfactory coding in the antennal lobe of the bumble bee Bombus terrestris

Sociality is classified as one of the major transitions in evolution, with the largest number of eusocial species found in the insect order Hymenoptera, including the Apini (honey bees) and the Bombini (bumble bees). Bumble bees and honey bees not only differ in their social organization and foraging strategies, but comparative analyses of their genomes demonstrated that bumble bees have a slightly less diverse family of olfactory receptors than honey bees, suggesting that their olfactory abilities have adapted to different social and/or ecological conditions. However, unfortunately, no precise comparison of olfactory coding has been performed so far between honey bees and bumble bees, and little is known about the rules underlying olfactory coding in the bumble bee brain. In this study, we used in vivo calcium imaging to study olfactory coding of a panel of floral odorants in the antennal lobe of the bumble bee Bombus terrestris. Our results show that odorants induce reproducible neuronal activity in the bumble bee antennal lobe. Each odorant evokes a different glomerular activity pattern revealing this molecule’s chemical structure, i.e. its carbon chain length and functional group. In addition, pairwise similarity among odor representations are conserved in bumble bees and honey bees. This study thus suggests that bumble bees, like honey bees, are equipped to respond to odorants according to their chemical features.

www.nature.com/scientificreports/ functional studies of odor coding in bumble bees are scarce. In the 1980's, two studies described bumble bees' peripheral equipment in cuticular sensilla on the antennae and performed electroantennogram (EAG) recordings of their antenna, showing that it responds to a wide range of volatiles, including both floral and pheromonal odorants 39,40 . These approaches were used again later to show that sensillar equipment and olfactory sensitivity increase with worker size in bumble bees 41 as well as to study left-right asymetries 42 . With regards to neural odor coding, one study demonstrated the existence of glomerulus-size odor-induced oscillations in the bumble bee AL 43 . More recently, extracellular recordings of AL neurons showed reproducible responses to odorants, and observed a specific response pattern for a pheromonal compound compared to other odorants 32 . Apart from these findings, little is known about the rules underlying olfactory coding in the bumble bee brain.
In the present work, we used in vivo calcium imaging to study olfactory coding by l-ALT PNs in the AL of the bumble bee Bombus terrestris. To compare odor-coding rules in bumble bees and honey bees, we presented a panel of floral odorants previously used in studies on olfactory processing and perception in honey bees 13,38,44 . Our results show that odorants induce reproducible neuronal activity in the bumble bee antennal lobe. Each odorant evokes a different glomerular activity pattern depending on the molecules' chemical structure, i.e. carbon chain length and functional group. Odor-similarity relationships in the bumble bee AL are highly correlated to those found in the honey bee AL.

Results
Anatomy of the bumble bee antennal lobe. Olfactory sensory neuron innervation. Using fluorescent tracers, we performed mass staining of olfactory sensory neurons (OSNs) in the bumble bee Bombus terrestris (Fig. 1A). The tracers migrated along the antennal nerve until the OSNs' axonal projections in the cortex (outer layer) of the glomeruli in the AL. As previously reported 33 , we found a similar arrangement of sensory tracts in the bumble bee antennal lobe as in the honey bee. The most prominent tract, T1, is easily identifiable, crossing the center of the antennal lobe from the antennal nerve caudally (Fig. 1B) to the most ventral and rostral part of the antennal lobe where it innervates many glomeruli. The T3 tract is also prominent, leaving the antennal nerve on the caudal side of the antennal lobe, propagating medially on its outskirts and innervating many glomeruli on the dorso-caudal region. T3 divides itself into at least 3 sub-branches: two running medially (T3a and T3b) and innervating many medial glomeruli, and one running laterally (T3c) innervating caudo-lateral glomeruli. Tract T2 is a much smaller tract that goes from the nerve entrance through the medial part of the lobe neuropil at approximately half depth and innervates only a few medial glomeruli. Tract T4 is another smaller tract, which runs laterally along the outer side of the glomerular region, and innervates a set of tear-shaped glomeruli on the most dorsal part of the antennal lobe, close to the dorsal lobe. Contrary to other glomeruli with a clearly stained cortex, these T4 glomeruli are characterized by a homogeneous staining of sensory neurons. A conspicuous tract of neurons bypasses the antennal lobe completely on its dorso-lateral side and forms the two tracts innervating the dorsal lobe (T5) and the subesophageal zone (T6) that transmit mechanosensory and gustatory information respectively. Single glomeruli from the confocal images were reconstructed (Fig. 1C) and we found 158 ± 4 glomeruli in the antennal lobe of bumble bees (n = 4 bumble bees), a slightly lower number compared to honey bees (~ 160-166 glomeruli) 15,45,46 , which roughly corresponds to the number of OR genes found in bumble bees 1 . We observed that in bumble bees, as in honey bees, the outer surface of the antennal lobe consists of a single layer of glomeruli. This arrangement is particularly well adapted to optical measurements of glomerular activity (see below).
Projection neuron innervation. Further similarities between honey bee and bumble bee olfactory systems were observed at the level of projection neurons innervation 32 . Using the classical technique used in honey bees 38,47 , we stained the lateral antennal-lobe tract (l-ALT) of projection neurons (Fig. 1A). By introducing tracers into the protocerebrum at a location lateral to the α-lobe of the mushroom bodies and rostral to the lateral horn, we obtained clear staining of l-ALT PNs ( Fig. 2A) in rostro-ventral glomeruli of the antennal lobe (Fig. 1C, red glomeruli innervated by T1 mainly). In contrast to anterograde staining of OSNs, PN staining was found to be homogeneous in the whole volume of the glomeruli, with PN somata visible on the edge of the glomerular area of the AL ( Fig. 2A).
All these observations confirm that OSN and l-ALT PN glomerular innervations are highly similar in bumble bees compared to the ones observed in the honey bee olfactory system. We then wondered whether similar olfactory coding rules are found in the two olfactory systems.
In vivo calcium imaging. We performed in vivo calcium imaging measurements in 14 bumble bee AL using the calcium indicator Fura-2 dextran, and recorded calcium responses from the dendrites of l-ALT PNs in rostro-ventral glomeruli (T1 region). We studied the coding of floral odorants. Inspired by previous work on honey bee olfactory perception and coding 13,38,44,47 , we presented to the bumble bees a set of 16 odorants differing systematically in their functional group and chain length.
Intensity of odor-induced responses. All odorants induced remarkable activity in a combination of AL glomeruli, while air control stimulation had no effect (Figs. 2B, 3A; n = 14; RM-ANOVA, F 16, 208 = 17.1, p < 0.0001, comparisons to the control: Dunnett test, p < 0.01). As the odorants systematically varied in terms of chemical group and carbon chain length, we evaluated the effect of these properties on the intensity of calcium responses. Odorants with different functional groups induced different activity levels (Fig. 3B, RM-ANOVA F 3,39 = 22.1, p < 0.0001). Among functional groups, the weakest responses were evoked by primary alcohols, which induced significantly lower responses than the other chemical groups (Tukey HSD test: p < 0.01 compared to secondary alcohols and p < 0.001 compared to ketones and aldehydes), which did not differ from each other. Odorants Generally, global response intensity decreased with increasing chain length, i.e. odorant molecules with 6 and 7 carbons induced stronger neural activity than odorants with 8 and 9 carbons (Tukey HSD test: a vs. b: p < 0.01). This pattern of results recapitulates the observations made in honey bees 38,44,47 and can be explained by the volatility of the odorants, as measured by their individual vapor pressure. Indeed, AL response was highly correlated to vapor pressure (Fig. 3D, R 2 = 0.88, F 1,14 = 106.7, p < 0.0001) confirming that the more volatile the odorant (i.e. the larger its vapor pressure), the more molecules were present in headspace in the sample and the larger was the recorded AL response to this odorant. In the presented odorant panel, alcohols and molecules with longer carbon chains possess lower volatility and thus induced lower responses.
Similarity among odor response maps. We then evaluated how chemical characteristics of odorants affected similarity relationships among AL response maps. We thus calculated a measure of (dis-)similarity between response maps (pixelwise Euclidian distance) for all possible pairs of the 16 tested odorants, and produced a distance matrix, which provides an overview of similarity relationships among these odorants (Fig. 4). The more similar odor responses were between two odorants, the smaller are the Euclidian distances and the more intense is the color in the matrix. The matrix reveals a strong effect of the odorant's carbon chain length on similarity relationships, as shown by the red diagonal lines in the matrix (e.g. for primary alcohols vs. secondary alcohols). Generally, distances between any two odorants of the same carbon chain length were smaller than distances between odorants with different carbon chain lengths. A remarkable exception was observed with secondary alcohols which appeared to be closer to the primary alcohol with a shorter chain by one C atom (see 44,47 for a similar effect in honey bees). The matrix also suggests that odor pairs with longer carbon chains (C8 vs. C9) evoke more similar activation patterns (i.e. smaller Euclidian distances) than odor pairs with shorter carbon chain length (C6 vs. C7). This more pronounced similarity is also visible in single recordings, as for example shown in Fig. 2B, where a distinct change in the glomerular activity map can be seen between C7 and C8 odorants, but not between C6 and C7 or between C8 and C9 molecules. Odorants' functional group also plays a role in similarity relationships, although this effect is less easily visible in the matrix. Some pairs of functional groups show higher similarity than others, for instance most primary and secondary alcohols show a high similarity (low distance). We confirmed these observations by performing multidimensional analyses using these Euclidian distance measures (Fig. 5). A hierarchical cluster analysis using Ward's classification method (Fig. 5A) showed that the odorants formed three main clusters. Odorants primarily segregated along two branches. The upper branch predominantly grouped odorants with short carbon chain lengths (C6 and C7). Within this branch, odorants were www.nature.com/scientificreports/ grouped according to their functional groups, with primary and secondary alcohols in one subgroup (C-OH functional group) and aldehydes and ketones in the other (C=O functional group). The lower branch exclusively contained odorants with longer carbon chain lengths (C8 and C9). Within this branch, odorants also tended to be distributed according to their functional group, apart from 1-nonanol. This analysis shows, as can be seen in the matrix (Fig. 4) as well as in individual recordings (Fig. 2B), that long-chain molecules evoke highly similar activity patterns, which are less dependent on the functional group than shorter molecules. www.nature.com/scientificreports/ We next performed a multidimensional scaling analysis 48 , using the Euclidian distance matrix, to understand the most meaningful dimensions underlying similarity relationships among odorants, focusing on the 3 main dimensions (Fig. 5B). Dimension 1 mostly provided information about odorants' carbon chain length, as odorants are represented along this axis by increasing carbon chain length for all functional groups (Fig. 5B left). Dimension 2 contained both functional group and chain length information. Primary and secondary alcohols were not separated from each other, but these C-OH functional groups were clearly separated from both ketones and aldehydes with a C=O functional group (Fig. 5B left). Dimension 2 also contains carbon chain length information for primary and secondary alcohols, as the odorants are represented along this axis by increasing carbon chain length. Lastly, dimension 3 clearly separates aldehydes (lower values) from ketones, primary and secondary alcohols (higher values, Fig. 5B right). To summarize, the three main dimensions of the multidimensional scaling analysis represented odorants' chain length and functional group information, distinguishing alcohols, ketones and aldehydes from each other.
The observations made on the distance matrix (Fig. 4) and the multidimensional analysis (Fig. 5) are supported by statistical analyses (Fig. 6). First, odor-specific coding is demonstrated by the fact that odor response maps for presentations of the same odorant were more similar (smaller Euclidian distances) than odor response maps for presentations of two different odorants (Fig. 6A, Paired t-test, t = 6.94, p < 0.0001, 13 df). Second, odorants with the same functional group induced significantly more similar odor response maps compared to odorants with different functional groups (Fig. 6B, Paired t-test, t = 4.69, p < 0.001). Lastly, odorants with the same carbon chain length induced more similar response maps than odorants with different carbon chain lengths (Fig. 6C, Paired t-test, t = 4.99, p < 0.001). This effect increased with the difference in the number of carbon atoms between the odorant molecules (Fig. 6D). The difference between odor maps was thus stronger when the molecules differed by at least 2 carbons, i.e. C6 vs. C8 or C6 vs. C9 (ANOVA F 3,39 = 25.86, p < 0.0001; Tukey HSD test: a vs. b p < 0.01; a vs. c p < 0.001). These analyses thus demonstrate that odor coding in the bumble bee AL relies on both odorants' chain length and odorant's functional group.
Euclidian distances measure (dis-)similarity between odor-response maps by taking into account differences both in the pattern of activated regions (glomeruli) and in intensity between responses. To concentrate more on the activity patterns, these analyses were also performed by using the correlation coefficient between pixel responses as a measure of similarity between odor-response maps. These analyses provided exactly the same results. First, correlation coefficients were higher (i.e. the maps were more similar) for presentations of the same www.nature.com/scientificreports/ odorant than for presentations of different odorants (Fig. S1A, Paired t-test, t = 11.0, p < 0.0001, 13 df). Second, correlation coefficients were higher for odorants with the same functional group than for odorants with different functional groups (Fig. S1B, Paired t-test, t = 5.88, p < 0.0001). Lastly, correlation coefficients were also higher for odorants with the same carbon chain length compared to odorants with different carbon chain lengths (Fig. S1C, Paired t-test, t = 4.46, p < 0.001). As for Euclidian distances, this difference was significant when the odorants differed by at least 2 carbons (Fig. S1D, ANOVA F 3,39 = 47.05, p < 0.0001; Tukey HSD test: a vs. b p < 0.01; a vs. c p < 0.001).
Comparison of honey bee and bumble bee data. The results we have described so far for bumble bees are generally very similar to the data obtained when imaging the homologous region of the honey bee AL 38,47 . We thus assessed the similarity of odor coding in bumble bees and honey bees by comparing odor-evoked intensity and similarity relationships between the two species. We performed linear regression analyses of response intensity (Fig. 7A) and similarity measures between bumble bees and honey bees (Fig. 7B). Response intensities measured for the 16 odorants were highly correlated (R 2 = 0.57; F 1,14 = 18.64, p < 0.001) showing that odorants inducing strong responses in bumble bees also induce strong activity in honey bees. This expected observation is a direct effect of response intensity directly depending on odorant vapor pressure in both species. . Short linkage distance between branches indicates odorants with similar response maps. Functional groups are shown in different colors: primary alcohols in blue, secondary alcohols in green, aldehydes in black, and ketones in red. The analysis shows a first separation (node 1) between odorants with short and long carbon chain lengths. Odorants with a short carbon chain are then subdivided (node 2) into alcohols (primary and secondary, C-OH functional group) and ketones/aldehydes (C = O functional group). (B) Multidimensional scaling analysis based on the Euclidian distance matrix for the 16 odorants. The first dimension (left panel) explains 34% of overall variance and orders molecules according to their chain length from short (on the right, C6 and C7) to long (on the left, C8 and C9). The second dimension explains 22.5% of variance and distinctly separates alcohols (blue, green) from ketones (red) and aldehydes (black). Functional group separation is clearer for short-chain than for long-chain molecules. The third dimension (right panel) explains 9.7% of variance and separates aldehydes (black) from other molecules. Altogether, odorants' chain length and functional group represent main coding dimensions for odorants in l-ALT projection neurons. www.nature.com/scientificreports/  www.nature.com/scientificreports/ With regards to similarity relationships, Euclidian distances between odor response maps for the 120 odor pairs were strongly correlated between honey bees and bumble bees (R 2 = 0.55; Mantel test p < 0.0001). This result is illustrated by the similar arrangements of odorant representations observed in the multidimensional scaling analysis performed on bumble bee data (Fig. 5B) and in that previously performed on honey bee data (see Fig. 4A, l-ALT PNS, in 38 ). Similarly, similarity measures based on correlation coefficients also yielded a significant correlation between honey bee and bumble bee data (data not shown; R 2 = 0.21; Mantel test p < 0.005). These results indicate that odorants inducing similar activity patterns in the honey bee AL also induce similar activity patterns in the bumble bee AL.

Discussion
This study shows that bumble bees are a suitable model organism for studying central olfactory processing and coding. Using neuroanatomical and neurophysiological approaches, we described AL architecture and measured glomerular activity patterns in response to odorants. We found that bumble bee l-ALT PNs provide clear information about odorants' chemical features, here their functional group and chain length. Although odor-evoked activity patterns were not topologically the same in bumble bees and honey bees, the general rules of olfactory coding in terms of intensity and similarity relationships were similar.
Neuroanatomical staining and 3D reconstructions indicated that the structure of the bumble bee AL greatly resembles that of the honey bee 45,49 . Both consist of a single layer of glomeruli around an inner coarse neuropil characterized by the presence of numerous local interneurons and projection neurons (see Fig. 1). The restricted innervation of the glomerular cortex by olfactory sensory neurons (OSN) seen in bumble bees is also reminiscent of the honey bee AL 45 . To note, in both species, a set of the most dorsal glomeruli innervated by the T4 tract present a homogeneous innervation compared with the exclusively peripheral innervation of the other glomeruli. The existence of the 4 OSN tracts in the bumble bee antennal lobe (Fig. 1B, plus 2 bypassing tracts towards dorsal lobe and subesophageal zone), and their similar arrangement to that observed in the honey bee 45,49 suggest a strong homology between the olfactory systems of both insects. The most prominent T1 tract innervates a large proportion of glomeruli on the rostro-ventral surface of the AL, which are directly accessible when opening the brain capsule. As in honey bees, these glomeruli could be stained retrogradely by placing dye crystals on the l-ALT tract of projection neurons (Fig. 2) 32,50 . We are thus confident that the group of glomeruli that we imaged in bumble bees is structurally homologous to the glomeruli usually imaged in honey bees with the same preparation 20,47,50 .
Using in vivo optical recordings, we demonstrated that a panel of 16 aliphatic odorants evokes reproducible neuronal activity in the glomeruli of the bumble bee antennal lobe. In particular, we found that an odorant's functional group and chain length influence the intensity of odor-evoked signals, with primary alcohols inducing significantly lower activity than other functional groups, and short chain molecules activating glomeruli more strongly than molecules with longer chain lengths (Fig. 3). These effects, also found in honey bees 38,44,47 are explained by the strong correlation found between response intensity and odorant's vapor pressure. Thus, the bumble bee antennal lobe, as its honey bee counterpart, does not display any specific sensitivity for any of the odorants in our panel, and the intensity of AL activity mainly reflects odor concentration in vapor phase.
We then showed that each odorant evokes a specific glomerular activity pattern (Fig. 6A), which is different from that evoked by other odorants. Olfactory coding was influenced by both tested chemical features, carbon chain length and functional group, as previously observed both in invertebrates 38,44,47,51,52 and in mammals 53,54 . Odor-evoked activity within PNs is the product of OSN activity entering the AL and of local inhibitory networks carrying out local computations 38,50,55 . Despite a highly similar organization of their olfactory pathways, the main difference between honey bee and bumble bee systems lies in the repertoire of ORs expressed at the periphery. A recent study analyzed the genomes of two Bombus species, Bombus terrestris and Bombus impatiens 1 and aiming to identify key genes in the evolution of sociality, they compared the genomes of these species with that of the honey bee Apis mellifera. Concerning chemoreception, they found that Bombus genomes contain a slightly less numerous OR family than Apis mellifera, with 159 intact OR genes (excluding 5 pseudogenes). The number of glomeruli that we found in the AL of Bombus terrestris in our reconstructions (n = 158 ± 4) corresponds well to the number of OR proteins found in the genome of this species, fitting with the general hypothesis in insects that each OSN expresses one type of odor-specific receptor, while all OSNs carrying the same receptor project to the same glomerulus in the AL 56 . This hypothesis is especially appealing as in honey bees the number of olfactory receptor genes largely coincides with the number of glomeruli in the AL (~ 165 glomeruli and ~ 163 intact OR genes excluding pseudogenes 57 ). Note however, that in Drosophila melanogaster, several AL glomeruli are not innervated by OSNs expressing OR family genes, but rather by neurons expressing ionotropic receptors 58,59 . As bumble bee and honey bee genomes each contain about ~ 20 IR genes, with several orthologous genes, a proportion of their ALs may be innervated by IR expressing sensory neurons.
The comparison of the honey bee and the bumble bee OR family genes also showed duplications of genes in one or both species, several large species-specific gene lineage expansions, and at least 22 gene losses, reflecting the typical birth-and-death evolution of these receptors 7 . Recent evidence in different species of the genus Drosophila suggests that the number of olfactory receptor genes has remained quite similar for the entire period of Drosophila evolution (63 million years 60 ), but that frequent gains and losses of genes occurred in each evolutionary lineage 61 . This may have changed the sequence of olfactory receptor neurons leading to different glomerular wiring patterns. The most recent common ancestors of honey bees and bumble bees are estimated to have lived between 70 and 90 million years ago [62][63][64] . This long time of separate evolution suggests that profound changes could also have taken place in the sequences of olfactory receptor genes in both species, modifying each receptor's sensitivity spectrum to odorant molecules, as well as the localization of its corresponding glomerulus in the AL. In accordance with these observations, direct comparison of odor-evoked activity patterns in bumble bees www.nature.com/scientificreports/ and honey bees was difficult and we did not manage to identify possibly homologous glomeruli with a sufficient level of confidence (see Fig. S2 for an example). Our finding of highly similar olfactory coding rules, supported by a clear coding of chain length and functional group information in bumble bees as in honey bees does not imply that ORs/glomeruli should be conserved between the two species. In our view, this finding is simply a by-product of the joined sensitivities and selectivities of the numerous OSN/ORs imaged simultaneously within each species. In a previous study, we showed that similarity relationships among inter-odorant maps measured in the ALs of an ant (Camponotus fellah) and the honey bee using calcium imaging were similar to those measured in the rat olfactory bulb using an utterly different recording technique (2-deoxyglucose autoradiography) 52 . The general rule was simple: odorants with a similar molecular structure (chain length and/or functional group) induced similar activity patterns in each insect's antennal lobe as well as in this mammal's olfactory bulb. More generally, it was observed that similarity relationships in a range of different species, including invertebrates and rodents, could be predicted to some extent based on purely molecular descriptors of odorants 65 . Thus, an emerging property of the multiple coding channels of each system (the ORs) which each detect a different but overlapping range of odorant molecule features is that neural representations mirror the chemical characteristics of the molecules. Importantly, this emerging property does not depend on the type of receptors expressed at the periphery since it is well established that olfactory receptor (OR) proteins in insects and vertebrates are unrelated 66,67 . Coming back to our study, honeybee and bumble bee ORs may well have evolved independently for a long time, our study shows that the two neural ensembles that were recorded in their ALs perform a reliable depiction of odorants' chemical features, granting these insects with a clear representation of odorants' structure. In honey bees, we previously showed that inter-odorant similarity relationships in the AL could predict bees' behavioral responses in a generalization protocol, so that similar odorants in the AL were treated as similar by the bees in their behavior 13,38,44 . The high correlation we found between odor-evoked response maps in the bumble bee and honeybee ALs suggests that bumble bees are indeed equipped, like honey bees, to respond to odorants according to chemical dimensions. We thus predict that future behavioral experiments in bumble bees shall reveal a similar organization of their olfactory perceptual space based on odorants chemical dimensions, as found in honey bees 13 .
In conclusion, our study unravels a high similarity in the general organization of the primary olfactory processing center of bumble bees and honey bees. In addition, it shows similar olfactory coding rules conveying each system with a reliable depiction of odorants' chemical structure. While we concentrated here on the coding of general odorant features, we expect that future studies devoted to the coding of species-specific odorants, like social pheromones, may reveal more remarkable differences between both systems.

Methods
Bumble bee preparation. Medium-sized bumble bee Bombus terrestris workers were caught from an indoor colony (Koppert, Berkel en Rodenrijs, The Netherlands) and chilled on ice for 5 min until they stopped moving. Then, bumble bees were prepared following the standard preparation used to image the AL in honey bees 20,38,50 . In summary, the bumble bee's head was inserted and fixed in a plastic chamber with its antennae oriented to the front of the chamber. Using beeswax, the proboscis was flued at the front end of the holder to avoid movement of the brain during the experiment. Hairs on the top of the bumble bee head were removed and a pool was built with beeswax and pieces of plastic around the rostral part of the head capsule (behind the antennae). The pool was made waterproof with two-component epoxy glue (red Araldite, Bostik Findley, S.A.). A small window was then cut in the head cuticle from the bases of the antennae up to the ocelli, and glands as well as parts of the tracheal sheath were removed to expose the antennal lobes and parts of the protocerebrum. Finally, the pool was filled with some ringer solution (in mM: NaCl, 130; KCl, 6; MgCl 2 , 4; CaCl 2 , 5; sucrose, 160; glucose, 25; Hepes, 10; pH 6.7, 500 mOsmol; all chemicals from Sigma-Aldrich, Lyon, France), to avoid desiccation of the brain surface. Three hours prior to the experiment, a dye mixture was inserted into the brain with a broken borosilicate micropipette, aiming for the tract of l-ALT projection neurons, between the α lobe and the border of the optic lobe, rostrally from the lateral horn. The dye mixture consisted of the calcium-indicator Fura-2 dextran (10,000 kDa, Life technologies, France) and of tetramethylrhodamine dextran (10,000 kDa, Life technologies, France) for later anatomical observation, both in bovine serum albumin (2%).

Calcium imaging.
In vivo optical recordings were performed as described elsewhere 22,38,68 , with a T.I.L.L.
Photonics imaging system (Martinsried, Germany), under an epifluorescence microscope (Olympus BX51WI) with a 10 × water-immersion objective (Olympus, UMPlanFL; NA 0.3), which was dipped into the ringer solution covering the brain. Only one AL was recorded in each bumble bee. Images were taken with a 640 × 480 pixels 12-bit monochrome CCD camera (T.I.L.L. Imago) cooled to -12 °C. Odor presentation. A constant clean airstream, into which odor stimuli could be presented, was directed from a distance of 2 cm to the bumble bee's antennae. Odor stimuli (see below) were given at the 15th frame for 1 s (5 frames). Odor sources consisted in exchangeable Pasteur pipettes containing a piece of filter paper (1 cm 2 ) soaked with 5 µl of pure odorant (Sigma Aldrich, France).
In a first experiment, we tested 16 different aliphatic odorants that are part of floral blends bumble bees encounter while foraging 69 . The odorants differed systematically in terms of their carbon chain lengths (between 6 and 9 carbon atoms) and their functional groups (primary alcohol, secondary alcohol, aldehyde and ketone). www.nature.com/scientificreports/ As control stimulus, we used a pipette containing a clean piece of filter paper without odor solution. This stimulus set was also used in a recent calcium imaging study of PN responses in the honey bee AL 38 allowing the comparison of odor coding in honey bees and bumble bees. The olfactory stimuli were presented three times in a pseudo-randomized order, avoiding consecutive stimuli to contain the same functional group or the same carbon chain length.
Data processing and analyses. A total of 50 bumble bees were imaged, out of which 14 presented highquality signals and were kept for further analysis. Data were analyzed using custom-made software written in IDL 6.0 (Research Systems, Boulder, CO) 37 . Each odor presentation produced a four-dimensional array consisting of the excitation wavelength (340 or 380 nm), two spatial dimensions (x-and y-coordinates) along time (50 frames). First, the fluorescence ratio between excitation wavelengths at each pixel and time point was calculated: R = F 340nm /F 380nm . The relative fluorescence changes were then computed between the recorded odor responses R at each time point compared to the background fluorescence (before any odor presentation) R 0 , defined as the average of the three images before odor stimulus onset (frames [12][13][14]. Relative fluorescence changes were thus calculated as: ∆R = (R − R 0 )/R 0 . The two spatial dimensions were then filtered with a gaussian filter of window size 7 × 7 pixels to reduce photon noise. Lastly, possible irregularities of lamp illumination were corrected by subtracting the median pixel value of each frame from each single pixel of the corresponding frame. The amplitude of the odor-induced response was calculated by subtracting the average of three consecutive frames during the odor presentation (frames [17][18][19] from the average of 3 frames before stimulus onset (frames [12][13][14]. The response intensities presented as a function of the functional group or of the chain length of the odorants (Fig. 3B,C) were calculated by averaging within each animal the intensities recorded to the 4 odorants belonging to each group. Activity maps (Fig. 2) represent the average amplitude observed over the three presentations of each odorant, in a false-color code, from dark blue (no signal) to red (maximum signal). As unambiguous identification of identical glomeruli across individual bumble bees was not feasible, odor coding was analyzed over the entire surface of the AL using a pixelwise analysis that avoids any bias due to glomerular misidentification. It was previously shown in honey bees that results based on the pixelwise method lead exactly to the same conclusions as glomerular identification 37 . For each bee, a mask was precisely drawn along the edges of the AL to limit the measure of odor-evoked responses to the glomerular area. Global glomerular activity upon odor stimulation was measured by averaging the intensity values of all pixels within the unmasked area. Evaluation of (dis-)similarity relationships between odorant representations was performed by calculating pixelwise Euclidian distances for all pairs of the 16 odorant stimuli used (120 odor pairs). For all analyses, average values for the three presentations of each odorant were used except for the comparison of Euclidian distances for the same or different odorants ( Fig. 6A and Fig. S1A), where each single odorant presentation was used.
Anatomical staining. For antennal staining of the whole antennal lobe, the scapes of the antennae were carefully opened using a microscalpel, and the antennal nerve was cut with a borosilicate micropipette coated with tetramethylrhodamine dextran (10,000 kDa, Life technologies, France). Afterwards, animals were kept in a cool place until the next day to allow the dye to migrate to the AL and to stain OSN processes within the glomeruli. The brains were removed and fixed in 4% paraformaldehyde solution for at least 24 h. They were then dehydrated in ascending concentrations of ethanol, cleared and stored in methyl salicylate (Sigma-Aldrich, Lyon, France). Images of the tetramethylrhodamine-stained glomeruli were taken using a confocal laser scanning microscope (Zeiss, LSM 700) with a W Plan-Apochromat 20x/1.0 objective and a 555 nm excitation wavelength at 2 µm optical section thickness and pixel size of 0.31 µm × 0.31 µm. Recorded stacks of images were adjusted in brightness and contrast using imageJ (Rasband; National Institutes of Health, Bethesda, MD). Segmentation and anatomical reconstruction of the antennal lobe was performed using Amira (version 4.5.1 Mercury Computer Systems, Merignac, France). The neuraxis was used for all anatomical descriptions 70,71 .
After successful calcium imaging, the brains were removed and the same techniques as described above were used.

Statistical analysis.
Normality of the data was tested and confirmed for almost all data points using Shapiro-Wilk normality test. We thus applied parametric statistics over the whole study. When normality was not achieved for all data points in an analysis, the corresponding non-parametric test was performed. In all cases, both types of tests gave the same result, and therefore the text only describes parametric results. The intensities of responses to the different odorants were compared using ANOVA for repeated measurements. When significant, Dunnett's test was applied to compare the intensity of each response to a common reference, the air control. Odor-evoked response intensities between functional groups and chain lengths were compared using ANOVA for repeated measurements, followed by Tukey post-hoc tests for further analysis of statistically significant main effects. Paired t-tests were applied to compare Euclidian distances obtained for different presentations of the same odor versus presentations of different odors, as well as for odors with the same or with a different functional group or chain length. A Pearson correlation analysis was performed between response intensity and the logarithm of odorants' vapor pressure. In some analyses (Fig. 7), data recorded in bumble bees were compared to data recorded in honey bees 38 , using exactly the same experimental and analytical procedures. A Pearson correlation analysis thus evaluated a possible correlation of odor-response intensities between the two species. A Mantel test was used to evaluate a possible correlation between bumble bee and honey bee Euclidian distance matrices. All tests were performed with GraphPad Prism (version 7, GraphPad software) or R (www.r-proje ct. org). All values are displayed as means ± SEM.