- Research article
- Open access
- Published:
Morphological and Transcriptomic Analysis of a Beetle Chemosensory System Reveals a Gnathal Olfactory Center
BMC Biology volume 14, Article number: 90 (2016)
Abstract
Background
The red flour beetle Tribolium castaneum is an emerging insect model organism representing the largest insect order, Coleoptera, which encompasses several serious agricultural and forest pests. Despite the ecological and economic importance of beetles, most insect olfaction studies have so far focused on dipteran, lepidopteran, or hymenopteran systems.
Results
Here, we present the first detailed morphological description of a coleopteran olfactory pathway in combination with genome-wide expression analysis of the relevant gene families involved in chemoreception. Our study revealed that besides the antennae, also the mouthparts are highly involved in olfaction and that their respective contribution is processed separately. In this beetle, olfactory sensory neurons from the mouthparts project to the lobus glomerulatus, a structure so far only characterized in hemimetabolous insects, as well as to a so far non-described unpaired glomerularly organized olfactory neuropil in the gnathal ganglion, which we term the gnathal olfactory center. The high number of functional odorant receptor genes expressed in the mouthparts also supports the importance of the maxillary and labial palps in olfaction of this beetle. Moreover, gustatory perception seems equally distributed between antenna and mouthparts, since the number of expressed gustatory receptors is similar for both organs.
Conclusions
Our analysis of the T. castaneum chemosensory system confirms that olfactory and gustatory perception are not organotopically separated to the antennae and mouthparts, respectively. The identification of additional olfactory processing centers, the lobus glomerulatus and the gnathal olfactory center, is in contrast to the current picture that in holometabolous insects all olfactory inputs allegedly converge in the antennal lobe. These findings indicate that Holometabola have evolved a wider variety of solutions to chemoreception than previously assumed.
Background
Insects use chemical cues for most tasks they encounter during their life history. Over long distances, airborne chemical stimuli guide insects to food sources, mates, and places for oviposition [1–6]. Within close range, olfaction as well as gustation are used to discriminate between different food qualities, to avoid toxins or harmful microbes, to communicate intra- or interspecifically, to identify suitable mating partners, and to find appropriate egg-laying sites [6–15]. Because of insects’ devastating impact on agriculture and stored food products, as well as their ability to serve as vectors for detrimental diseases, insect olfaction has become an important research field in biology [4].
Chemical signals are typically perceived within specialized antennal and palpal cuticular structures, the olfactory or gustatory sensilla. These chemosensory sensilla form a hollow structure filled with aqueous lymph and harbor the dendritic branches of the chemosensory neurons (CSNs), namely the olfactory (OSNs) or gustatory sensory neurons (GSNs) [16, 17]. They are divided into several sub-types according to their different morphology [16]. The volatile molecules enter the cavity through wall pores finally to reach and activate the chemoreceptors on the dendrites of the OSNs. To enhance olfactory sensitivity and specificity, odorant binding proteins (OBPs) or potentially chemosensory proteins (CSPs) facilitate the translocation of many, mostly hydrophobic, chemicals through the aqueous lymph [18]. In insects, typically three different receptor families are involved in chemoreception [4]: the ionotropic glutamate-like receptors (IRs) [19, 20], the gustatory receptors (GRs) [13, 21], and the odorant receptors (ORs) [22–24]. The IRs are evolutionarily highly conserved chemoreceptors involved in protostome olfaction [19]. They contain three transmembrane domains and form functional heteromers between an odor-specific IR and a co-receptor (IR8a and IR25a) The GRs are seven transmembrane receptors found across arthropods [25–28] whose quaternary structure [13, 29–31], as well as the signal transduction mechanism [32, 33], are still under debate. The typical ORs are seven transmembrane receptors found in pterygote insects [24] that form functional heteromers with the atypical (general) odorant receptor co-receptor (Orco) [22, 23, 34–36]. Their signal transduction mechanism is currently discussed and they may either form an ionotropic receptor complex that is regulated by second messengers or be functional metabotropic receptors [22, 23, 34, 37–40]. The described influence of G-proteins and affiliated second messengers on insect olfaction supports both mechanisms [41–47]. Moreover, sensitive pheromone detection requires the OR/Orco complex to interact with a sensory neuron membrane protein (SNMP) related to the scavenger receptor CD36 [48–50]. Besides the perireceptor events involved in effective activation, the high temporal resolution of olfactory reception probably also requires signal termination, which is supposedly mediated by secreted or membrane-bound odorant-degrading enzymes (ODEs) [51–55].
Activation of the described chemoreceptors elicits action potentials in the CSNs that are further transmitted via the antennal nerve to the antennal lobe (AL), the first integration center of the olfactory pathway in the brain, or for GSNs, to the primary gustatory center of the gnathal ganglion (GNG) [56]. The AL of insects consists typically of spherical sub-compartments, the olfactory glomeruli [57]. Usually OSNs express only one typical (specific) OR gene and all antennal OSNs expressing the same typical OR converge into the same olfactory glomerulus, creating a chemotropic map-like representation of chemical coding in the AL [58–60], known as the central dogma of olfaction [61, 62]. In Drosophila melanogaster, the OR/Orco and IR derived sensory information from the antennae and the maxillary palps is processed in the AL [63], whereas in several hemimetabolous insects, CSNs from the palps converge typically in the lobus glomerulatus (LG), next to but outside the AL [57, 64–66]. In the AL, olfactory information from the OSNs, is processed by a complex network of local interneurons [57, 67–69]. The processed odor information is further relayed by distinct antennal lobe tracts (ALTs) formed by the projection neurons (PNs) to the mushroom body (MB) and the lateral horn (LH) [57, 70]. The MBs are higher-order integration centers for multiple processed sensory information and are responsible for odor discrimination, associative learning, as well as memory storage and retrieval. The LHs receive odor input directly from the ALs or indirectly from the MBs, decode the quality and intensity of the information, and finally trigger immediate odor-driven behavior [71–77].
Despite the evolutionary success and ecological as well as economic importance of beetles [78, 79], little is known on the neuroanatomy, genetics, or biochemistry of their olfactory pathway. Within the Coleoptera, the red flour beetle Tribolium castaneum has become the prime model organism for developmental biology and pest management [80]. With its fully annotated genome [81–83] and the multiple powerful genetic tools – such as systemic RNA interference [84, 85], insertional mutagenesis [86], and transgene-based misexpression systems [87, 88] – T. castaneum represents an eligible beetle model organism for olfaction. In the current study, we present a substantial overview of the olfactory pathway in T. castaneum, covering the morphology of the sensilla and the antenna, all major neuropils including AL, MB, LH, LG, and the gnathal olfactory center (GOC), a previously undescribed glomerularly organized neuropil in the GNG. Additional support for the importance of the gnathal input into olfaction is provided by genome-wide expression analysis of gene families involved in chemoreception (e.g., ORs, GRs, IRs, SMNPs, and ODEs) and CSPs and OBPs, which have recently been published [89].
Results
The Antenna of Tribolium castaneum
To determine the distribution and number of CSNs, we used immunohistochemistry (IHC) with a cross-reactive antibody against Orco, fluorescent in situ hybridization with an Orco-specific probe, and a transgenic line, EF1-B-DsRed, that labels almost all and only CSNs in the adult antenna (see ‘Methods’ for a detailed characterization). Moreover, we generated an Orco-Gal4 line that partially covers the Orco pattern, which we refer to as the partial Orco-Gal4 line (see ‘Methods’ for a detailed characterization of reagents). These different approaches unequivocally confirm that CSNs are restricted to the distal three segments (9–11) that form the enlarged club of the antenna [90] (Fig. 1a; Additional file 1: Figure S1a; and Additional file 2: Figure S2a). To improve on previous data in respect to the characterization, location, and exact number of antennal sensilla [90], we used in addition to the confocal laser-scanning microscopy (CLSM) approaches also scanning electron microscopy (SEM) (Figs. 1b–h and 2a–g). This morphologically verified the presence of chemosensory sensilla exclusively on the three club segments [90], with the highest number and diversity on the apical part of the terminal segment 11 (Fig. 1b–b''; Additional file 3: Figure S3).
Sensilla types and distribution on Tribolium castaneum antennae I. a Chemosensory sensilla are restricted to the distal three segments (9–11) of the T. castaneum antenna, which is composed of scape (S), pedicel (P), and flagellum, and the last labial palp (LP) and maxillary palp (MP) segment. CLSM-stack voltex projection of a transgenic beetle head (ventral view, green: partial Orco-Gal4/UAS-tGFP; yellowish eye, brownish cuticle: autofluorescence). b–b'' SEM images of the club segments with close-up of segments 9 (b') and 11 (b''). Single sensilla: CLSM maximum intensity projection overlays (c–h) of antibody-enhanced EF1-B-DsRed reporter signal (magenta, c'–h') and cuticle autofluorescence (green, c''–h''). c'''–h''' SEM analysis. Mechanoreceptive sensilla: SCam are small, smooth, and dome-shaped sensilla restricted to segment 11 (Additional file 3: Figure S3a); SCha – previously described as spines [90] – are longitudinally corrugated, connected to a neuron at the socket (c'; blue), jointed (c'''; arrow), and solid (c''''; arrowhead). d–d'''' SpaB – in T. brevicornis called sensilla squamiformium [95] – resemble modified (slightly thicker tip) SCha [96] restricted to segment 11 (Additional file 3: Figure S3b). e–e'''' mSTri (structurally similar to SCha but smaller more hair-like appearance) have previously been described in other species [24, 244]. CLSM analysis showed joint-like structures at the base (c–e, c''–e'', open squares) of the mechanoreceptive sensilla and SEM revealed a small gap at their base (c'''–e''', arrow). Chemoreceptive sensilla: f–f''' cSTri are hair-like structures restricted to segment 11 (Additional file 3: Figure S3d) with a rounded tip and a smooth transition of the base; g–g'''' SBas are smooth-surfaced pegs with rounded tips and smooth transitions at the base (g'''; arrow). h–h''' SCoe are short and corrugated, and their transition into the antennal cuticle shows a typical elevation (b'', h'''). All chemoreceptive sensilla (f, g, f'–h') house dendritic branches of CSNs labeled by DsRed. The close-up in c' shows a non-CSN fiber entering only the base of a SCha labeled with phalloidin (blue). Chemoreceptive sensilla show a smooth transition into the antennal cuticle (f'''–h''', arrow). Whereas all mechanoreceptive sensilla are solid cuticular structures (fractured in c''''–e''''), chemoreceptive SBas appear hollow (fractured in g''''). CLSM confocal laser-scanning microscopy, CSN chemosensory neuron, cSTri chemosensilla trichoidea, LP labial palp, MP maxillary palp, mSTri mechanosensilla trichoidea, P pedicel, S scape, SBas sensilla basiconica, SCam sensilla campaniformes, SCha sensilla chaetica, SCoe sensilla coeloconica, SEM scanning electron microscopy, SpaB spaculate bristle
Sensilla types and distribution on Tribolium castaneum antennae II. a–e SEM images of SBas with one to five prongs. f, f' SEM image of the tenth segment of the antenna with a close-up of the lateral corner (f') containing SCoe, SBas, and mSTri. g, g' SEM image of the ninth segment with a close-up of the lateral corner (g') showing SCoe and mSTri. Chemoreceptive SCoe were previously described as “minute spicule-like sensilla trichoidea” [90], are relatively rare (Additional file 3: Figure S3e), and located besides the lateral corners of segments 9 and 10 (f',g') mostly at the apical side of segment 11 (Fig. 1b''). Chemoreceptive SBas are arranged in an axial ring at the distal margins of all three club segments (f, g, Fig. 1b–b''). For mechanoreceptive mSTri, we identified about 37 on the apical side of segment 11 (Fig. 1b'') and four in lateral corners of segments 9 and 10 (f, f', g, g', and Additional file 3: Figure S3c, h). h Voltex projection based on a CLSM image stack of the tenth segment from the EF1-B-DsRed line displaying CSNs (orange) and autofluorescence of the cuticle (green). The dendrites of the CSNs converge into the SBas (on average, six per prong), while the axons unite at the center of the segment and join the antennal nerve (AN). i–i'' Overlay of the signals of the DsRed reporter (magenta, i') and the Orco antibody (green, i'') together with DAPI staining (light blue) in the EF1-B-DsRed line, demonstrating a high level of colocalization between DsRed and Orco in segments 9 and 10, but not in 11, where some DsRed-immunoreactive CSNs are spared (compare with Additional file 1: Figure S1a). AN antennal nerve, CLSM confocal laser-scanning microscopy, CSN chemosensory neuron, mSTri mechanosensilla trichoidea, Orco odorant receptor co-receptor, SBas sensilla basiconica, SCoe sensilla coeloconica, Seg segment, SEM scanning electron microscopy
Four mechanoreceptive and three chemoreceptive sensilla types could be confirmed by the combination of these techniques (Fig. 1b–b'') and the respective number of contained CSNs was identified. The mechanoreceptive sensilla include the spatulate bristles (SpaB; Fig. 1d–d''''), the mechanosensilla trichoidea (mSTri; Fig. 1e–e''''), the sensilla campaniformes (SCam; Fig. 1b''), and the sensilla chaetica (SCha; Fig. 1c–c''''), which are the most dominant sensilla type present on the lateral sites of all 11 segments (Fig. 1a). The chemoreceptive sensilla subdivide into chemo-sensilla trichoidea (cSTri, Fig. 1f–f'''), sensilla basiconica (SBas; Fig. 1g–g''''), and sensilla coeloconica (SCoe, Fig. 1h–h'''). For the chemoreceptive sensilla, segments 9 and 10 carry mostly SBas (about 15) arranged in an axial ring at the apical edge of each segment (Fig. 2f, g) and two SCoe (Fig. 2f, g'), whereas the terminal segment 11 harbors SBas (about 25), some SCoe (about 7), and many cSTri (about 87) (Fig. 1b–b''). A detailed analysis of the number and distribution of the different sensilla types in males and females revealed no sexual dimorphism (Additional file 3: Figure S3).
The number of CSNs per antenna was estimated based on the number of CSNs per sensillum or prong and the number of the respective sensilla per antenna. cSTri contain typically one Orco-immunoreactive OSN (Additional file 1: Figure S1b; Additional file 3: Figure S3i). This type of sensilla is known for its pheromone receiving abilities in Lepidoptera [91–93] and had been described as olfactory sensilla in D. melanogaster [59] and Culex quinquefasciatus [94]. SBas of T. castaneum consist of up to five prongs (Fig. 2a–e) like other Tenebrionidae [90, 95]. Each prong harbors about six CSNs (Additional file 3: Figure S3i) – the same number as in Tenebrio molitor [96]. Of them, four or five can be considered olfactory based on Orco-immunoreactivity (Additional file 1: Figure S1c). Findings in Tribolium brevicornis [95] suggest an additional gustatory function of SBas, leading to the conclusion that the SBas are bimodal chemosensilla. Because of this constant number of CSNs per prong and the shared lymphatic space (Fig. 1g–g''), we propose that multiple pronged SBas are derived from a fusion of single sensilla. Nonetheless, we refer to and count multiple pronged SBas as a single sensillum independent of the number of prongs. SCoe contain three CSNs (Additional file 3: Figure S3i) without Orco-immunoreactivity (Additional file 1: Figure S1d). The SCoe in T. castaneum might therefore harbor IRs as shown in D. melanogaster [19, 97]. Altogether, we found on each antenna about 100 prongs of SBas with six CSNs each, 87 cSTri with one CSN, and 11 SCoe with about three CSNs (Additional file 3: Figure S3). This leads to a total number of about 720 CSNs per antenna of T. castaneum.
Anatomy of the olfactory pathway in the red flour beetle brain
Antennal projections
To get an impression of the innervation pattern of chemosensory neuropils, we performed antennal and palpal backfills. Backfills via the antennal nerve labeled the ipsilateral AL (Fig. 3a; Additional file 4: Movie S1; Additional file 5: Movie S2; Additional file 6: Figure S4), the antennal mechanosensory and motor center (AMMC) (Fig. 3b), as well as a distinct area in the GNG (Fig. 3c). While this ipsilateral restriction is common in many insects [57], it is in contrast to D. melanogaster and Ceratitis capitata, where the majority of OSNs innervate the ipsi- and contralateral sides [59, 98]. The antennal backfills labeled all AL glomeruli except one, which is the only glomerulus labeled by backfills of the maxillary palp via the GNG (Fig. 3a; Additional file 5: Movie S2). This resembles the situation in Lepidoptera, where CO2 responsive CSNs from the palp project into a single AL glomerulus devoid of antennal innervation [99]. The descending antennal projections into the GNG (Fig. 3c) are not labeled in the partial Orco-Gal4/UAS-DsRed line and might therefore be from gustatory or mechanosensory neurons, as described in Periplaneta americana and Locusta migratoria [100, 101].
The central olfactory pathway of T. castaneum. a Backfill of one antenna (magenta) stains all glomeruli in the ipsilateral antennal lobe (AL) except one. This glomerulus is exclusively labeled by a backfill of a maxillary palp (cyan). b In addition to the AL glomeruli, backfilling (magenta) of one antenna labeled the ipsilateral antennal mechanosensory and motor center (AMMC), located n-dorsally to the AL, c as well as descending fibers to the gnathal ganglion (GNG). d Maximum intensity projection of the backfills of mouthparts (cyan) shows massive innervation of the GNG including the gnathal olfactory center (GOC) (magnified in the inset) and the primary gustatory center (PGC). e Backfill of the mouthparts (cyan) revealed in the cerebral ganglion beside innervation of a single ipsilateral AL glomerulus also projections in the ipsilateral lobus glomerulatus (LG). f Reporter expression of the partial Orco-Gal4/UAS-DsRed line (magenta) revealed two paired input tracts (black and white arrowheads) from the maxillary (white arrowhead) and labial palps (black arrowhead) that converge in a medial and n-anterodorsally located glomerular area, the GOC, and ascend to a microglomerularly organized area, the LG. See also Additional file 7: Movie S3. Orientation bars in (a) also apply for (b) and (e). AL antennal lobe, AMMC antennal mechanosensory and motor center, GNG gnathal ganglion, GOC gnathal olfactory center, L lateral, LG lobus glomerulatus, NA neuroaxis-anterior, PGC primary gustatory center, TR tritocerebrum
Orco-immunoreactive sensory neurons in the maxillary palp. a Voltex projection of a CLSM-stack showing antibody enhanced reporter expression of the EF1-B-DsRed line (a', orange) and Orco-immunoreactive cells (a'', green) in a halved maxillary palp. b–b'' Single optical section of (a) showing partial colocalization of Orco immunoreactivity and the reporter expression of the EF1-B-DsRed line (magenta). Dotted lines in (b) highlight reporter-expressing cells that are not Orco-immunoreactive. CLSM confocal laser-scanning microscopy, Orco odorant receptor co-receptor
Additional file 4: Movie S1: 3D reconstructions of the antennal nerve, antennal lobe, antennal mechanosensory and motor center, and lobus glomerulatus. Z-stack video of a phalloidin stained brain with embedded 3D reconstruction of antennal lobe (dark blue), antennal nerve (light blue), antennal mechanosensory and motor center (turquoise), and the lobus glomerulatus (magenta). Later the neuropils are embedded in a voltex projection of the brain, also based on phalloidin staining (orange). (MP4 11379 kb)
Additional file 5: Movie S2: Camera path through a confocal stack of the AL with backfills of the antenna and maxillary palp. Obtained from the same confocal stack as Fig. 3a. Antennal backfill in green and maxillary palp in magenta. (MP4 5185 kb)
Antennal lobe
For the AL of freshly eclosed adults, about 70 distinguishable olfactory glomeruli have been previously described using a synapsin antibody [102]. To evaluate the glomeruli number in ALs of beetles, 7 days after adult eclosion, we improved the analysis by deconvolution as well as using an additional antiserum against tachykinin-related peptides (TKRPs), which distinctly labels also densely packed glomeruli [103]. This more advanced analysis resulted in the 3D reconstruction of about 90 glomeruli per AL with no obvious sexual dimorphism (females: mean 89.2, standard deviation or SD = 4.9, n = 5; males: mean 89.4, SD, 7.6, n = 5).
Palpal projections into accessory olfactory centers
Whole mouthparts or maxillary palp backfills (Fig. 3d, e) revealed besides the already mentioned single AL glomerulus, innervation of three distinct neuropil areas: an unpaired glomerular organized neuropil in the GNG, the primary gustatory center also in the GNG [104], as well as an area near the AL, resembling the LG of hemimetabolous insects [57, 105]. The unpaired neuropil located n-anterodorsal in the GNG consists of 30 to 40 glomeruli (Fig. 3d, inset), which are all innervated from both sides of the mouthparts. This neuropil is also labeled by the partial Orco-Gal4/UAS-DsRed line (Fig. 3f; Additional file 7: Movie S3), which indicates innervation by OSNs originating in the maxillary or labial palps (Figs. 1a, 3f, 4a'', 4b''; Additional file 2: Figure S2f) that project via two tracts into the GNG. This neuropil, therefore, represents an olfactory processing center in the GNG that has to our knowledge never been described before and we term the gnathal olfactory center (GOC). Some of the fibers labeled by the palpal backfills, as well as the partial Orco-Gal4/UAS-DsRed line pass through the GOC, ascend via the neck connectives, and terminate ipsilaterally in an area medioventral to the AL (Fig. 3d), resembling the LG, which to date had only been described in hemimetabolous insects [57, 105]. Since the position, innervation, and glomerularly organized structure of this paired neuropil in T. castaneum is similar to the LG in cockroach, locust, and silverfish [57, 64–66], we refer to it as LG. In summary, our data suggest that in T. castaneum, odor information from the antennae and the mouthparts are processed separately. It appears that OSNs from the mouthparts do not project into the AL but into the GOC and the LG.
Additional file 7: Movie S3: Voltex projection of the gnathal ganglion and part of the brain of the Orco-Gal4/UAS-DsRed line. The video was obtained from the same confocal stack as Fig. 3f. It shows two paired input tracts from the maxillary and labial palps that converge in the GOC and ascend to the LGs, as well as the partially labeled ALs. (MP4 16700 kb)
Projection neurons
Dye injections into the AL of adult T. castaneum revealed three ALTs formed by the PNs (Fig. 5), exclusively in the ipsilateral hemisphere. The most prominent tract, the medial antennal lobe tract (mALT), connects the AL with the calyx (CA) of the MB and the LH. The mediolateral antennal lobe tract (mlALT) passes the region near the MB spur without forming sub-branches (Fig. 5) and further projects to the LH. The lateral antennal lobe tract (lALT) projects directly to the most n-posterior part of the LH. We could not observe any obvious direct projections of the mlALT and the lALT to the CA. However, since they possibly overlap with trajectories of the mALT fibers from the CA to the LH [106], we cannot exclude their existence, as described for other holometabolous insects [70]. Previously only the mALT had been clearly identified in Coleoptera and the existence of a mlALT had only been presumed [70]. Our results indicate that three ALTs are a common feature among most holometabolous insects, including beetles.
Antennal lobe tracts. Maximum intensity projection of a CLSM image stack after dye injection into the AL (magenta) revealed three antennal lobe tracts – the medial (mALT), mediolateral (mlALT), and the lateral antennal lobe tract (lALT) – as well as the calyx (CA) and the lateral horn (LH). In the CA, most fibers from the mALT form microglomeruli (inset obtained from another preparation). The staining in the optical lobe is an artifact caused by diffusion of the dye during application. Phalloidin counterstaining in green. AL antennal lobe, ALT antennal lobe tracts, CA calyx, CLSM confocal laser-scanning microscopy, lALT lateral antennal lobe tract, LH lateral horn, mALT mediolateral lobe tract, mlALT mediolateral lobe tract
Mushroom body
The detailed architecture of the MB of T. castaneum is described in [103]. The CA is innervated by the mALT (Fig. 5) and microglomerularly organized as indicated by phalloidin or synapsin antibody stainings (Fig. 5, inset). This is similar to several insects including Apis mellifera [107, 108] and D. melanogaster [109, 110] and suggests a comparable wiring with the PNs. The Kenyon cells (KCs) were identified in DAPI (4',6-diamidino-2-phenylindole) stainings based on their smaller and brighter stained nuclei [103]. The number of about 2700 KCs was determined by interpolation of volumetric data as well as by counting of the stained nuclei using MorphoGraphX [111]. Both procedures resulted in comparable numbers with the interpolation of 13 CAs from seven animals estimating about 2800 KCs (2795; SD: 214) and the counting of nine CAs from five specimen indicating approximately 2600 KCs (2613; SD: 204) per MB.
Genome-wide expression analysis of genes involved in chemoreception in T. castaneum
The fully sequenced genome of T. castaneum [81–83] led to the annotation of the major gene families involved in chemoreception. Based on genome data and computational gene predictions, the OBPs [112], CSPs [113], IRs [114], GRs [81], ORs [115], and SNMPs [116, 117] were annotated, but only for the ORs was a RT-PCR-based expression analysis performed [115]. To validate or correct the predicted gene models of these gene family members and to determine their tissue-specific expression, we performed transcriptome analyses of adult male and female antennae, heads (without antennae, but including mouthparts), mouthparts (the part of the head capsule, anterior to the antennal bases), legs, and bodies (without legs and head). In addition, we identified potential ODEs, as well as orthologs from further genes described to be involved in D. melanogaster olfaction. The detailed analysis of the OBPs and CSPs has already been published [89] and revealed that the majority of the classic OBPs and antenna binding proteins II (ABP II) seem to be involved in chemoreception while only a few of the minus-C OBPs (C-OBPs) and CSPs are enriched in antenna or mouthparts. The following results are based on this same set of transcriptome data (GEO accession number: GSE63162 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63162 [118]). Like OBPs and CSPs [89], also for the genes presented here, no significant differences on the expression level between male and female antenna samples were identified (Fig. 6). Therefore, the female and male antenna samples can serve as biological replicates and indicate that reads above 0.1 reads per kilobase per million (RPKM) are reproducible (Fig. 6). However, to minimize the rate of potential false positives in our description, we considered only genes with RPKM ≥ 0.5 as expressed. We are aware that this might lead to an underestimation of the expressed gene numbers for each class of genes. Since it is impossible to determine the exact number of genes that are functionally involved in chemoreception based only on transcriptomic expression analyses, we always present two values for expressed genes, one based on RPKM ≥ 0.5 and the other defined by statistical analysis as significantly enriched over body. All raw values and the re-annotated gene models are summarized in Additional file 8.
Comparison of expression levels in male and female antenna. Comparison of expression levels of odorant receptors (ORs, magenta), gustatory receptors (GRs, green), ionotropic glutamate-like receptors (IRs, blue), sensory neuron membrane proteins (SNMPs, orange), orthologous of candidates obtained from D. melanogaster (Dmel candidates, grey) and potential odorant degrading enzymes (ODEs, yellow) in male and female antennae. Average values based on two male and three female antennal samples. Scatter plot of the RPKM values. Dmel D. melanogaster, GRs gustatory receptors, IRs ionotropic glutamate-like receptors, ODEs odorant degrading enzymes, ORs odorant receptors, RPKM reads per kilobase per million, SNMPs sensory neuron membrane proteins
Tissue-specific expression of ionotropic glutamate-like receptors
The RNAseq based revision of the 23 previously annotated IRs [114] confirmed the sequences of three open reading frames (ORFs); 17 had to be modified, two were incompletely covered by reads, and for a single one, no expression was detected (color coded in Additional file 8: Table S1, column B). In antennae, 16 of the IRs were significantly enriched compared to body (Fig. 7; Additional file 9: Figure S5a). In the mouthparts, five IRs are expressed, two are significantly enriched.
Expression of T. castaneum ionotropic glutamate-like receptors (IRs). Heat map showing the expression level of the 23 IRs as a log2[RPKM + 1] value in different tissues [adult antennae, head (missing antennae but including mouthparts), mouthparts, legs, and body]. The candidates are ordered according to their chromosomal localization (Additional file 9: Figure S5b). Horizontal brackets above indicate clustering in the genome. The arrowheads represent the orientation of the open reading frame. The expression levels are represented by a greyscale with highest shown expression levels labeled black. The asterisks mark statistically significantly differentially expressed genes compared to body (based on biological replicates of five antennal, two head, three mouthpart, two leg, and two body samples). The red asterisks represent up- and the blue down-regulation (p values adjusted are * < 0.05, ** < 0.01, and *** < 0.001). IR ionotropic glutamate-like receptor, RPKM reads per kilobase per million
Comparing expression profiles of the IRs from T. castaneum, D. melanogaster, and Anopheles gambiae confirmed the antennal specific expression, as well as the high degree of phylogenetic conservation of the antennal IRs (Fig. 8; Additional file 10: Figure S6; highlighted in yellow) as proposed [114]. In contrast, the divergent IRs are non-antennal specifically expressed and are highly radiated within species clades as previously shown or predicted [114]. T. castaneum has a lower number of IRs compared to D. melanogaster and An. gambiae, due to lesser expansions of divergent IRs, but maintains the basic repertoire of antennal IRs (Fig. 8; highlighted in yellow). The homologs of IR25a, IR93a, and IR40a, which are necessary for humidity perception in D. melanogaster [119], are significantly enriched in antennae. IR40a is exclusively expressed in antennae, which correlates with the essential role of antennae in T. castaneum hygro-perception [90]. The homolog of the highly sensitive salt receptor and possible co-receptor IR76b [20, 120] is significantly enriched in antennae, mouthparts, and legs, while the co-receptors IR8a and IR25a [20] are highly expressed in all tissues of T. castaneum (Fig. 7).
Phylogenetic tree of IRs. Based on protein sequences from T. castaneum (green branches), D. melanogaster (red branches), and An. gambiae (blue branches). The tree was rooted using the IR8/IR25 clade, according to [114]. Robustness of the tree topology was evaluated by 100 rapid bootstrap replications. Outer rings represent the expression in antennae and mouthparts (T. castaneum: palps, mandible, labrum, and labium; D. melanogaster: palp and proboscis; An. gambiae: maxillary palp) as log2-fold change compared to body corresponding to the scale in the left lower corner. The scale bars within the trees represent one amino acid substitution per site. Antennal IRs are highlighted in yellow. Basically the same figure is available with absolute values instead of fold changes to get an impression of the tissue-specific abundance of the transcripts as Additional file 10: Figure S6. IR ionotropic glutamate-like receptor
Tissue-specific expression of gustatory receptors
Of the 220 previously annotated GRs [81], only 207 genes had available gene models [82, 83]. Our transcriptome analysis verified the ORFs of 58 GRs, showed slight differences for 20 GRs, but did not or only incompletely cover 129 GRs (Additional file 8: Table S1, column B). In the antennae, 62 GRs are expressed, with 34 being significantly enriched and 10 being antennal-specific. Of the 69 mouthpart-expressed GRs, 36 are significantly enriched and 19 exclusive. Seventeen GRs are significantly enriched in both antenna and mouthparts. In legs, 18 GRs are expressed with three being significantly enriched (Fig. 9; Additional file 11: Figure S7a).
Expression of T. castaneum gustatory receptors (GRs). Heat map showing the expression level of the 207 analyzed GRs as a log2[RPKM + 1] value in different tissues [adult antennae, head (missing antennae but including mouthparts), mouthparts, legs, and body]. The candidates are ordered according to their chromosomal localization (Additional file 11: Figure S7b). Horizontal brackets above indicate clustering in the genome. The arrowheads represent the orientation of the open reading frame. The expression levels are represented by a greyscale with the highest shown expression levels labeled black. The asterisks mark statistically significantly differentially expressed genes compared to body (based on biological replicates of five antennal, two head, three mouthpart, two leg, and two body samples). The red asterisks represent up- and the blue down-regulation (p values adjusted are * < 0.05, ** < 0.01, and *** < 0.001). CO2 receptors are highlighted in orange, fructose receptor related genes in grey, and sugar receptors in yellow. GR gustatory receptor, RPKM reads per kilobase per million,
The phylogenetic comparison of the GRs in T. castaneum, D. melanogaster, and An. gambiae (Fig. 10; Additional file 12: Figure S8) confirmed that only the CO2 receptors (highlighted in orange) are highly conserved [121]. The other GRs seem to have undergone independent radiation during the transition to T. castaneum: e.g. the sugar receptor-related branch (highlighted in light yellow) contains 16 genes [122], twice the number compared to the two chosen dipterans. In addition, the single fructose receptor (highlighted in grey) found in D. melanogaster and An. gambiae is represented by eight homologs in T. castaneum. The remaining 180 GRs belong to several T. castaneum-specific expansion groups. Specific orthologs to the known bitter receptors of D. melanogaster [13] as well as to the thermo-sensitive DmelGR28bD [123] cannot be predicted based on our phylogenetic analysis.
Phylogenetic tree of gustatory receptors (GRs). Mid-point rooted tree based on protein sequences from T. castaneum (green branches), D. melanogaster (red branches), and An. gambiae (blue branches). Robustness of the tree topology was evaluated by 100 rapid bootstrap replications. Outer rings represent the expression in antennae and mouthparts (T. castaneum: palps, mandible, labrum, and labium; D. melanogaster: palp and proboscis; An. gambiae: maxillary palp) as log2-fold change compared to body corresponding to the scale in the left lower corner. The scale bars within the trees represent one amino acid substitution per site. Potential sugar receptors (highlighted in yellow), fructose receptors (highlighted in grey), and CO2 receptors (highlighted in orange) are labeled. Known bitter receptors from D. melanogaster are highlighted in green, and the thermos-sensitive GR28bD in light blue. Basically the same figure is available with absolute values instead of fold changes to get an impression of the tissue-specific abundance of the transcripts as Additional file 12: Figure S8. GR gustatory receptor
Like other insects [121, 124, 125], T. castaneum has three CO2 receptors (TcasGR1, TcasGR2, and TcasGR3), while D. melanogaster has only two that form functional heteromers [126, 127]. In T. castaneum, the expression of the CO2 receptors is not restricted to one of the chemosensory organs with TcasGR2 and TcasGR3 being significantly enriched in antennae but also being expressed together with TcasGR1 in the mouthparts (Fig. 9; highlighted in orange). This dual input is in contrast to but combines both the expression of the three An. gambiae CO2 receptors that are restricted to the maxillary palps [128, 129], as well as the two D. melanogaster CO2 receptors that are mainly expressed in the antennae [126, 127, 130].
The presence of GRs on insect antenna had previously been postulated based on physiological response to sugars [95, 131–133] and was identified by antennal expression analysis [124, 127, 128, 134–136]. Our interspecies comparison (Fig. 10) confirms the antennal enrichment of several GRs in the two analyzed dipterans. However, the high number of 34 significantly enriched GRs in the antenna of T. castaneum is unusual, but reflects the increased total number of GRs in this species. Interestingly, the GRs of T. castaneum are present in both antenna and mouthparts at similar numbers and expression levels (Fig. 9; Additional file 11: Figure S7a).
Tissue-specific expression of odorant receptors
Of the 341 previously annotated OR sequences [115], we could re-analyze 337 based on our RNAseq data. This revision confirmed 97, and 22 were re-annotated reviving eight previously indicated pseudogenes [115], namely TcasOR2, TcasOR18, TcasOR19, TcasOR22, TcasOR85, TcasOR99, TcasOR104, and TcasOR122. Moreover, 219 genes were not or only partially covered by our transcriptome data (color coded in Additional file 8: Table S1, column B). Over all samples, 170 ORs are expressed (Fig. 11; Additional file 13: Figure S9a). In antennae, 129 ORs are expressed, with 92 being significantly enriched and 99 exclusive. In the mouthparts, 49 ORs are expressed, with 28 being significantly enriched and 27 exclusive. In addition, 16 of the significantly mouthpart-enriched ORs are not enriched in the antenna (Fig. 11). The expression of typical ORs in the mouthparts is consistent with the high expression of Orco in this tissue (Figs. 1a and 4) and with observations in other insect species [128, 129, 137, 138]. In legs, ten ORs are expressed (Additional file 13: Figure S9a) but only one, namely TcasOR127, is statistically enriched (Fig. 9).
Expression of T. castaneum odorant receptors (ORs). Heat map showing the expression levels of the 337 analyzed ORs as log2[RPKM + 1] with a maximum of 8.1 (Orco has a value of 11.1 in antenna) in different tissues [adult antennae, head (missing antennae but including mouthparts), mouthparts, legs, and body]. The candidates are ordered according to their chromosomal localization (Additional file 13: Figure S9c). Horizontal brackets above indicate clustering in the genome, and the arrowheads represent the orientation of the open reading frame. ORs that are member of clades four, five, and six [115] are written in grey letters. The line labeled with Adult and Larva refers to data from [115]. The character H (respectively B) indicates that the corresponding OR was detected in head or body cDNA samples by reverse PCR of the labeled developmental stage. A black letter indicates that an amplicon was detected in the majority of replicates, a grey letter means only in a few replicates, a dash indicates no PCR product and no character means no data available. A comparison of the number of expressed genes is summarized in Additional file 13: Figure S9b. The expression levels are represented by a greyscale with highest shown expression levels (3 RPKM or higher) labeled black to make sure that also low level expression is identifiably presented. The asterisks mark statistically significantly differentially expressed genes compared to body (based on biological replicates of five antennal, two head, three mouthpart, two leg, and two body samples). The red asterisks represent up- and the blue down-regulation (p values adjusted are * < 0.05, ** < 0.01, and *** < 0.001). B body, H head, OR odorant receptor, RPKM reads per kilobase per million
The phylogenetic comparison of OR expression patterns in T. castaneum, D. melanogaster, and An. gambiae (Fig. 12; Additional file 14: Figure S10) revealed that the atypical odorant co-receptor Orco (in T. castaneum previously called TcOR1 [115]) is the highest expressed OR in all tissues of all three species. In T. castaneum, Orco is expressed highest in antenna, followed by mouthparts. Orco is the only OR of T. castaneum with clear orthologs in dipterans [115, 139]. The high expression levels, the distribution, and the evolutionary conservation of Orco are consistent with its ancestral origin [24] and its outstanding role as a chaperone and co-receptor, forming functional heteromers with all typical ORs [140, 141].
Phylogenetic tree of odorant receptors (ORs). Protein sequences (>300 amino acids) from T. castaneum (green branches), D. melanogaster (red branches), and An. gambiae (blue branches). The tree was rooted using the Orco clade, according to [24]. Robustness of the tree topology was evaluated by 100 rapid bootstrap replications. Outer rings represent the expression in antennae and mouthparts (T. castaneum: palps, mandible, labrum, and labium; D. melanogaster: palp and proboscis; An. gambiae: maxillary palp) as log2-fold change compared to body corresponding to the scale in the left upper corner. The surrounding numbers on the outer thin line indicate the expansion groups 1 to 6 [115]. TcasOR71 and TcasOR72PSE were previously assigned to expansion group 1. The scale bar within the tree represents one amino acid substitution per site. Basically the same figure is available with absolute values instead of fold changes to get an impression of the tissue-specific abundance of the transcripts as Additional file 14: Figure S10. OR odorant receptor
The exceptional high number of typical ORs (Fig. 12) in T. castaneum is the result of large gene radiations within the coleopteran and tenebrionid lineages [136], which were previously subdivided into six expansion groups (Fig. 12) [115]. Expansion groups 1, 2, and 3 are conserved in other coleopterans [136] and are mainly expressed in antennae. The ORs of the expansion groups 4, 5, and 6 are highly derived, have no described homologs in other insects, and their expression is unusually often mouthpart-enriched (Fig. 11; grey lettering). This is consistent with the elaborated role of the mouthparts in T. castaneum olfaction. Specific orthologs to deorphanized ORs of D. melanogaster [142] cannot be predicted based on our phylogenetic analysis.
Identification and expression of potential odorant degrading enzymes
The genome of T. castaneum contains 15 aldehyde dehydrogenases (ALDHs) (Fig. 13 and Additional file 15) with two of them being significantly enriched, but not exclusively expressed in antenna. We found four predicted genes encoding aldehyde oxidases (ALOXs) with one being highly enriched in antennae and mouthparts, which, in contrast to ALOX ODEs from lepidopterans [143–145], does not encode a signal peptide (Fig. 13). Five of the 54 identified carboxylesterases (CESs) are significantly enriched in antenna, with two of them also in the mouthparts. Two other CESs are significantly enriched exclusively in the mouthparts. Five of these seven candidates show a predicted signal peptide for secretion (Fig. 13). TcasCESXA shares sequence similarities with D. melanogaster Est6, and TcasCES7J with DmelJHEdup, with both D. melanogaster homologs having previously been identified as ODE candidates [54, 55]. TcasCES10C is expressed highest in antennae and related to a pheromone degrading enzyme from the Japanese beetle, Popillia japonica [146]. We identified six epoxide hydrolases (EHs), which are supposed to be membrane bound ODEs [147], with one being significantly enriched in antennae and having a predicted signal peptide (Fig. 13). The glutathione S-transferases (GSTs) of T. castaneum had already been annotated [148]. The revision confirmed most gene models, only TcasGSTd2 and TcasMGST2 had to be modified (available in Additional file 8: Table S1). Eight of the 41 GSTs are significantly enriched in antennae, with three also in the mouthparts (Fig. 13). One of these three, TcasGSTd2 represents a member of the GST delta subfamily such as GST-msolf1 from Manduca sexta, which is an olfactory-specific GST expressed specifically in the sex-pheromone-detecting sensilla [149]. Analysis of the 141 previously described cytochrome P450s (CYPs) [150] revealed that two predicted gene models (CYP347A4 and CYP351B1) were fusions of two separate genes (now termed CYP347A4A and CYP347A4B, as well as CYP351B1A and CYP351B1B, respectively). Seven other predictions had to be adjusted based on RNAseq data (sequences available in Additional file 8: Table S1). The expression analysis of these 141 genes showed that 26 are significantly enriched in the antenna, with 11 also in the mouthparts (Fig. 13). In addition, six CYPs are significantly enriched in mouthparts, but not in antennae. For the coleopteran Phyllopertha diversa, CYPs have been shown to be involved in pheromone degradation in a membrane-bound manner [51].
Expression of T. castaneum potential odorant degrading enzymes (ODEs). Heat map showing the expression level of the 263 potential ODEs as a log2[RPKM + 1] value in different tissues [adult antennae, head (missing antennae but including mouthparts), mouthparts, legs, and body]. The candidates are ordered according to their protein family and chromosomal localization. Horizontal brackets above indicate clustering in the genome (Additional file 15: Figure S11), and the arrowheads represent the orientation of the open reading frame. Underlined genes were previously found on the protein level in antennae by [89]. The expression levels are represented by a greyscale with the highest shown expression levels labeled black. The asterisks mark statistically significantly differentially expressed genes compared to body (based on biological replicates of five antennal, two head, three mouthpart, two leg, and two body samples). The red asterisks represent up- and the blue down-regulation (p values adjusted are * < 0.05, ** < 0.01, and *** < 0.001). A black dot in the lowest line indicates a predicted signal peptide according to a SignalP 4.0 [211] prediction. ODE odorant degrading enzyme, RPKM reads per kilobase per million
Expression of potential olfaction signal transduction pathway components
The orthologs of genes encoding signal transduction pathway components known to be involved in olfaction of D. melanogaster [46] were identified by BLAST and manually curated. The expression analysis revealed that four of them (rdgB, itpr, dgkd, and dgkt) are significantly enriched in the antennae (Fig. 14). However, there is no chemosensory-specific candidate exclusively expressed in antennae or mouthparts. Our data, therefore, do not indicate a chemosensory-specific metabotropic signal transduction pathway.
Expression of T. castaneum homologs of genes described to be involved in olfaction of D. melanogaster. Heat map showing the expression level of the several genes supposed to be involved in D. melanogaster olfaction, as a log2[RPKM + 1] value in different tissues [adult antennae, head (missing antennae but including mouthparts), mouthparts, legs, and body]. The expression levels are represented by a greyscale with highest shown expression levels labeled black. The asterisks mark statistically significantly differentially expressed genes compared to body (based on biological replicates of five antennal, two head, three mouthpart, two leg, and two body samples). The red asterisks represent up-regulation (p values adjusted are * < 0.05). RPKM reads per kilobase per million
Expression and distribution of sensory neuron membrane proteins
The transcriptome analysis revealed that one of the seven previously identified TcasSNMPs [116, 117], namely XP_969729 [116], was incorrectly annotated and does not encode for a CD36-related protein. Moreover, the gene model previously named SNMP1c (XM_001816389) was a fusion of two SNMPs and overlaps with SNMP1d (XM_001816391) [117]. In our re-annotation, we removed XP_969729 and separated TcasSNMP1c and TcasSNMP1d. In addition, the gene models of TcasSNMP2, TcasSNMP1a, and XP_975606 [116] had to be modified based on transcriptome and RACE-PCR data. For XP_975606, we propose the name TcasSNMP3, to reflect its unclear phylogenic relationship. Despite the more SNMP1-like expression pattern (Fig. 15) and chromosomal localization (Additional file 9: Figure S5b) of TcasSNMP3, the comparison of the amino acid composition revealed no clear affiliation to either the SNMP1 or the SNMP2 subgroup [151]. Interspecies comparison revealed no clear orthology of TcasSNMP3 to SNMPs from other species, including the so-called SNMP3 of Calliphora stygia [152], which, based on phylogeny, clearly represents an SNMP1 homolog. All six TcasSNMPs are expressed in antennae (Fig. 15), which was also confirmed by rapid amplification of cDNA-ends PCR (RACE-PCR) based on an antennae cDNA pool, but only TcasSNMP1a-d and TcasSNMP3 are significantly enriched in antennal tissue. Moreover, three of the TcasSNMP1, as well as TcasSNMP3, are also enriched in mouthparts (Fig. 15), further supporting the importance of the mouthparts for olfaction in T. castaneum. In contrast, TcasSNMP2 is expressed highest in body and significantly underrepresented in antennae and mouthparts (Fig. 15), which is similar to its ortholog in D. melanogaster [49]. Despite the observation that in most insects with a fully sequenced genome only two SNMPs were found [116, 117], the relatively high amount of six TcasSNMPs of T. castaneum is not unique, since transcriptome analysis, e.g., of other beetles, revealed four SNMPs in Dendroctonus valens [153] and Dastarcus helophoroides [154], as well as three in Ips typographus and Dendroctonus ponderosae [136]. However, T. castaneum is currently only exceeded by the hessian fly (Mayetiola destructor) with seven expressed SNMPs [155].
Expression of T. castaneum sensory neuron membrane proteins (SNMPs). Heat map showing the expression level of the six sensory neuron membrane proteins of T. castaneum, as a log2[RPKM + 1] value in different tissues [adult antennae, head (missing antennae but including mouthparts), mouthparts, legs, body, as well as larval head and body]. The candidates are ordered according to their chromosomal localization (Additional file 9: Figure S5b). Horizontal brackets above indicate clustering in the genome, and the arrowheads represent the orientation of the open reading frame. The expression levels are represented by a greyscale with highest shown expression levels labeled black. The asterisks mark statistically significantly differentially expressed genes compared to body (based on biological replicates of five antennal, two head, three mouthpart, two leg, and two body samples). The red asterisks represent up- and the blue down-regulation (p values adjusted are * < 0.05, ** < 0.01, and *** < 0.001). SNMP sensory neuron membrane protein, RPKM reads per kilobase per million
Discussion
Independent integration centers for antennal and palpal olfactory perception
In T. castaneum, odorants are mainly perceived with the last three segments of the antenna, which carries three types of chemoreceptive sensilla (SBas, cSTri, and SCoe), as well as with the maxillary and labial palps (Fig. 16). Accordingly, expression analysis revealed that ORs are mostly expressed in antennae, but also in the mouthparts (Figs. 11 and 12; Additional file 13: Figure S9a) as previously shown for several dipteran species [128, 129, 137, 156–158]. In contrast to the Diptera, where the palps are chemosensory appendages with limited odor coding complexity, the relatively high number of Orco-immunoreactive CSNs (Fig. 4) as well as the high number of expressed ORs, SNMPs, potential ODEs, and OBPs [89] in T. castaneum mouthparts (Figs. 11, 13, 15, and 16) imply a more prominent role of the palps in olfaction. The palpal ORs are possibly involved in the evaluation of the quality of food sources, like the ORs on the proboscis of Manduca sexta [159].
T. castaneum head scheme depicting the major olfactory pathway components. a Head section (dorsal view) showing the brain and the CSNs from the antenna (blue) and the mouthparts (green). b Head section (ventral view) showing the GNG. Section orientation is indicated at the upper right corner (lateral view of the head; V, ventral). Double-headed arrows indicate body (A, anterior ↔ P, posterior; black) and neuro-axis (NA, n-anterior ↔ NP, n-posterior; red). Chemical signals are sensed by about 720 CSNs located in 56 SBas, 87 cSTri, and 11 SCoe on the last three antennal segments. These CSNs express 16 IRs, 62 GRs, 129 ORs, and six SNMPs. Chemosensory information is also perceived in the palps by five IRs, 69 GRs, 49 ORs, and six SMNPs (number in brackets indicates significantly enriched members compared to body). The antennal nerve (AN) projects into the ipsilateral AL, where all except one (light green) of the about 90 GL (dark blue) are innervated. A separate antennal tract (*1) descends into the GNG (b, blue), where presumably gustatory and mechanosensory information is processed. Incoming olfactory information is processed by a complex network of local interneurons (LNs) in the AL and further relayed by projection neurons forming three ALTs. The medial ALT (mALT) projects to and arborizes in the calyx of the MB formed by about 2700 KCs (orange) to eventually innervate the LH (light blue). The mediolateral ALT (mlALT) and lateral ALT (lALT) directly innervate the LH. From the mouthparts, CSNs project via the maxillary (*3) and labial palp nerves (*4) into the GNG, where the gustatory information is processed in the PGC. The olfactory sensory input from the palps is processed in an unpaired glomerularly organized GNG structure, the GOC, as well as in the LG, which receives input from some palpal OSNs via ascending neurons (*2) passing through the GOC. Some of the palp-derived chemosensory information is processed in the single AL glomerulus, which lacks antennal innervation and is, therefore, exclusively innervated by projections from the mouthparts (light green). AL antennal lobe, ALT antennal lobe tracts, AN antennal nerve, CSN chemosensory neuron, cSTri chemosensilla trichoidea, GL antennal lobe glomeruli, GNG gnathal ganglia, GOC gnathal olfactory center, GR gustatory receptor, IR ionotropic glutamate-like receptor, KC Kenyon cells, lALT lateral antennal lobe tract, LG lobus glomerulatus, LH lateral horn, LNs local interneurons, mALT mediolateral lobe tract, MB mushroom body, mlALT mediolateral lobe tract, OR odorant receptor, OSN, olfactory sensory neuron; PGC primary gustatory center, SBas sensilla basiconica, SCoe sensilla coeloconica, SNMP sensory neuron membrane protein
Moreover, in addition to the differences on the perception level, major dissimilarities to the Diptera occur on the level of odor processing. The data from the partial Orco-Gal4 line as well as the backfills from the antenna and the mouthparts indicate that processing olfactory information at least at the level of the first central relay station occurs independently of each other (Fig. 16). This is surprising, as many of the ORs expressed on the mouthparts are also expressed on the antennae. In contrast, typical OR expression is mutually exclusive between antenna and palps in D. melanogaster and An. gambiae [128, 137], where in addition, projections from the palps innervate several AL glomeruli [63, 160, 161]. In T. castaneum, the olfactory input stemming from the antenna seems to be processed exclusively in the AL (Fig. 3a; Additional file 5: Movie S2; Additional file 6: Figure S4), whereas the palpal-derived olfactory information is essentially processed outside the AL, in the LG (Fig. 3e, f; Additional file 5: Movie S2) and the GOC, an unpaired and glomerularly organized first olfactory center in the GNG (Fig. 3d, f; Additional file 7: Movie S3). The LG had, as far as we know, previously been described only in hemimetabolous insects [57, 64–66, 105]. A glomerularly organized olfactory center in the GNG such as the GOC has, to our knowledge, not been described in any insect so far. The number of 49 ORs (with 28 being significantly enriched compared to body) that are expressed in the mouthparts is roughly consistent with the estimated 30 to 40 glomeruli in the GOC. This suggests that the wiring in the GOC may resemble the situation in the ALs with the difference being convergence into an unpaired medial structure. The only palpal projection into the AL is a mutually exclusive innervation of a single ipsilateral glomerulus (Fig. 3a; Additional file 5), which may be involved in CO2 perception, as described in several moth species [99] and proposed for some mosquitoes [161, 162].
Antennae serve also as key organs for gustatory perception
In T. castaneum, antennae and mouthparts express similar high numbers and levels of GRs, which indicates the antenna as a key gustatory organ besides the mouthparts (Fig. 9; Additional file 11: Figure S7a). This finding may reflect the beetles’ ground-dwelling life style and indicates that the scanning behavior with the antennae not only gathers tactile but also chemical stimuli. This is in contrast to the Diptera, where the labellum is the main gustatory organ [163–165].
Postulation of exceptions to the central dogma
The number of 129 ORs that we found to be expressed in T. castaneum antennae (Fig. 11; Additional file 13: Figure S9a) exceed the numbere of about 90 glomeruli in the AL. Moreover, some glomeruli are likely to get exclusive innervation by OSNs that express IRs, as described in D. melanogaster [20]. These observations do not conciliate with the central dogma postulating that OSNs express only one typical OR and all OSNs carrying this same OR converge into one and the same glomerulus, which was hypothesized to be the typical situation for insects [61, 137, 166]. However, for D. melanogaster, both co-expression of more than one typical OR per OSN as well as co-convergence due to innervation of one AL glomerulus by more than one OSN sub-type have been already described as exceptions [137, 167]. For T. castaneum, we propose that such exceptions are much more frequent.
Large repertoire of potentially functional odorant receptor genes and possible environmental regulation
The genome of T. castaneum harbors 341 OR genes [81, 115], of which 270 seem to encode for functional ORs. Of the 337 ORs with available full sequence information [115], we find in our RNAseq data 161 ORs to be expressed in adult antennae, mouthparts, and head by a threshold of 0.5 RPKM (Additional file 13: Figure S9a). In comparison to the RT-PCR-based data from [115], who found 112 ORs to be clearly expressed in adult heads, we only confirmed 82 ORs. In addition, we identified 41 ORs previously declared as not expressed and 37 ORs previously not tested [115] as expressed (Additional file 13: Figure S9b). This discrepancy may partially be due to the different type of methodology used to identify expression. However, culturing conditions and the specific genetic variations of the strain used may also be responsible for the differences.
Taking both studies together, there is clear experimental evidence for 191 ORs that are expressed in the adult head. By including adult leg and all adult body data, 223 ORs seem to be expressed in total, of which 17 actually do not encode an intact OR. However, for 64 OR functional gene models, no expression could be detected so far. This may be due to low expression in a single OSN or conditional expression under exceptional circumstances. The red flour beetle can live for up to two years [168]. During this long period in their natural environment, the beetles can encounter a variety of challenges such as food shortages, which possibly triggers flight migrations over tens of kilometers [169]. Under such exceptional circumstances, the not or low-expressed receptor genes may become active [115], as shown in studies in D. melanogaster [170] and An. gambiae [171] where up to fivefold upregulation of several ORs was triggered by temperature or feeding state.
Inter-species comparison of olfactory components
The comparison of the number of main components of the chemosensory pathway of different insect species reveals the high diversity of evolutionary strategies to enable proper chemoreception and thus, reflects the diversity of insects and the manifold adaptations to their specialized lifestyles (Table 1). In particular, T. castaneum has by far the lowest number of chemoreceptive sensilla (154) and consequently also of CSNs (720). In contrast to this low number, the number of GRs (220) and ORs (341), but not of the IRs (23) encoded in the genome is exceptionally high. The number of olfactory glomeruli in the AL is within the range of most other species (Table 1) [57, 172]. Comparing the relation of OR genes and number of glomeruli, the highest discrepancy occurs with about fourfold higher numbers of OR genes in T. castaneum. However, also in Aedes aegypti, OR gene numbers are more than double the number of glomeruli [165]. In most other analyzed insects, except ensiferan orthopterans that have hundreds of microglomeruli [173, 174], the number of OR genes is typically similar to the number of glomeruli (Table 1). Despite the relatively low number of IRs encoded in the genome of T. castaneum, the repertoire of IRs involved in olfaction is highly conserved (Fig. 8). The number of KCs is roughly the same as in D. melanogaster and seems to be independent of the OR or AL glomeruli number (Table 1) [175].
No apparent sexual dimorphism
Sexual dimorphism of the olfactory system is described in many insect species [57, 172, 176, 177]. However, in contrast to other coleopterans [178–182], our analysis revealed no apparent sexual dimorphism on antenna morphology or number and distribution of sensilla (Additional file 1: Figure S1b–d). Expression analysis of male and female antenna samples revealed only a small but not significant dimorphism in the OBP expression levels described earlier [89]. Also for IRs, GRs, ORs, and SNMPs, we could not find any significant sexual differences (Fig. 6), similar to the striped flea beetle Phyllotreta striolata [183] and in contrast to described situations in Diptera and Lepidoptera [63, 128, 184, 185]. Different numbers of glomeruli or different sized glomeruli were observed in several insect species [57, 70, 176] including the beetle Holotrichia diomphalia [172]. However, the comparison of the ALs of T. castaneum males and females disclosed no obvious dimorphism as previously described also for the small hive beetle (Aethina tumida) [186]. In summary, our study did not reveal any sexual dimorphism of the olfactory system in T. castaneum. This finding is consistent with behavioral studies that showed an attraction of both sexes to the aggregation pheromone 4,8-dimethyldecanal [187] and no sex preference in the mating choice of males [188].
Conclusions
Detailed analysis of the olfactory system in T. castaneum, a holometabolous insect of special importance for the study of coleopteran and pest biology, reveals that olfactory sensory input from the antennae is processed mostly in the antennal lobes of the brain, as observed in other insect species. However, tracing of olfactory projections from the mouthparts enabled the identification of two additional neuropils: a lobus glomerulatus described previously only in a hemimetabolous insect and an unpaired glomerularly organized olfactory neuropil in the GNG (the GOC), which has never before been described. In addition, the high number of GRs on both the antennae and mouthparts indicates that there is no organotopic separation of olfaction and gustation in this beetle. These findings are a reminder of the wide variety of solutions to chemoreception that have evolved in the holometabolous insects. This should remind us that we have much still to learn about olfactory systems in general.
Methods
Tribolium castaneum rearing and transgenic lines
Tribolium castaneum (Herbst, 1797; Insecta, Coleoptera, Tenebrionidae) wild-type strain San Bernardino, as well as the transgenic lines partial Orco-Gal4, UAS-DsRed, UAS-tGFP [87], and EF1-B-DsRed [189] were bred at about 30 °C and 40 % relative humidity on organic whole wheat flour supplemented with 5 % yeast powder [190]. The Orco-Gal4 and UAS-DsRed lines were generated by piggyBac-based insertional mutagenesis [191]. The donor plasmids used were assembled by a versatile two-step cloning procedure [192].
For the partial Orco-Gal4 line, a donor plasmid was generated by cloning a blunted and BamHI (Fermentas, Vilnius, Lithuania) digested PCR < product containing Gal4delta-SV40pA (amplified with primers Gal4deltafor and SV40rev from plasmid CH#757, see Additional file 16) into the BamHI and EcoRV (Fermentas) digested pSLfa1180 vector [193]. After propagation, a BamHI and BfuAI digested PCR product containing 2.5 kb upstream of the TcasOrco (amplified with TcOR1upfor and TcOR1uprev from San Bernardino gDNA) was cloned into the corresponding restriction sites to generate pSLfa1180[2.5kbOrcoUp_GAL4delta]. The whole cassette was shuttled with AscI and FseI (New England Biolabs, Ipswich, MA, USA) into the pBac[3XP3-Tcv] [194] donor plasmid. The tissue-specific expression of Gal4 in the Orco-Gal4 line was determined by crossing it with an UAS-tGFP [87] line and performing IHC on the antennae with α-tGFP and α-Orco antibody or by staining of the whole brain with α-tGFP and an α-synapsin counterstaining. These stainings revealed that only Orco-immunoreactive neurons are labeled in antennae (Additional file 2: Figure S2d), which indicates the specificity of the Orco-Gal4 driver line. However, only half of the Orco-immunoreactive neurons in the antenna express tGFP (Additional file 2: Figure S2d), which implies that the Orco-Gal4 line only partially covers the Orco pattern resulting in labelling of only half of the AL glomeruli (Additional file 2: Figure S2e). The same approach with an UAS-dsRed line and an α-RFP antibody was used to characterize the palps, in which the reporter is also exclusively expressed in Orco-immunoreactive neurons, but in only 10–20 % of the cells (Additional file 2: Figure S2f). We, therefore, refer to it as the partial Orco-Gal4 line.
For UAS-DsRed, the donor plasmid pBac[3XP3-eYFP_UAS-Tchsp68bP-DsRedex-SV40] was generated by cloning the DsRed express ORF (Clontech laboratories Inc., Mountain View, CA, USA; catalog no. 632412) into the pSLfa[UAS-Tc'Hsp-p-tGFP-SV40]fa shuttle vector [87] by using KpnI and NotI, which was followed by transferring the UAS-hsp-DsRed-SV40 cassette into the pBac[3XP3-eYFP] [193] using AscI and FseI. The UAS-DsRed line as well as the UAS-tGFP line were analyzed by confocal microscopy to ensure that no reporter expression was present in the relevant tissues in the absence of a Gal4 driver line (Additional file 17: Figure S12).
The EF1-B-DsRed line (elongation factor1-alpha regulatory region-DsRedExpress; kindly provided by Michalis Averof, Institut de Génomique Fonctionnelle de Lyon, France) has been described to label most neurons in the central nervous system of first instar larvae [189] and also shows high expression in the adult central nervous system. However, clearly not all neurons are labeled in the peripheral nervous system. We, therefore, re-analyzed adult antennae of this line using confocal microscopy in combination with antibody stainings. The labeled neurons in the antenna resemble the typical morphology of CSNs with the dendrites being embedded in the sensilla cavities (Fig. 2h; Fig. 1f, g, and h) and the axons converging to the antennal nerve (Fig. 2h; Additional file 1: Figure S1a). No labelling was detected at mechanosensory sensilla (Fig. 1c, d, and e) except the scolopidia cells of Johnston's organ (Additional file 1: Figure S1a). In addition to almost all Orco-immunoreactive ORNs (Fig. 2i), this line labels also non Orco-immunoreactive neurons that are affiliated with sensilla coeloconica (Fig. 1h) and sensilla basiconica (Fig. 1g). Whereas in the palps only about 30 to 50 % of the DsRed-immunoreactive cells are also Orco-immunoreactive (Fig. 4), in antennal segment 11 a higher percentage of CSNs is double labeled, and in segments 9 and 10, the vast majority of CSNs are double labeled (Fig. 2i–i''; Additional file 1: Figure S1a). This suggests that almost all and only CSNs are labeled by this line in the adult antenna.
Tissue preparation for SEM
Antennae of sex-separated adults were dissected and immediately fixed for at least 2 h in 5 % glutaraldehyde in 0.1 M phosphate buffered saline (PBS), pH 7.1, washed and post-fixed in osmium-tetroxide (1 % in 0.1 M Sörensen buffer, pH 7.2). Fixed samples were washed in water, dehydrated overnight in ethyleneglycolmonoethylether, and then transferred into acetone via at least three 10-min changes with 100 % acetone as described in [195]. The samples were critical-point-dried by using a Polaron E 3000 (Balzers Union, Quorum Technologies Ltd, Darmstadt, Germany). After being sputtered with gold (Balzers Union Sputter Coater, Balzers, Liechtenstein; Quorum Technologies Ltd, Ringmer, UK), the material was examined using a Hitachi S-530 SEM (Hitachi High-Technologies Europe GmbH, Krefeld, Germany). Micrographs (Figs. 1b–b'', c''', c'''', d''', d'''', e''–e'''', f'', g''', g'''', h'', and 2a–g) were taken by digital image acquisition (DISS 5, point electronic, Halle, Germany).
Immunohistochemistry
Whole mount brain IHC was performed as described in [102]. The animals were cold anesthetized, their brains were dissected in cold PBS (0.01 M, pH 7.4), and fixed subsequently overnight at 4 °C or for 1–2 h at room temperature in PBS containing 4 % paraformaldehyde (Roth, Karlsruhe, Germany). The tissue was rinsed four times for 10 min with PBS. and pre-incubated with 5 % normal goat serum (NGS, Jackson ImmunoResearch, Westgrove, PA, USA) in PBT (PBS containing 0.3 % Triton X-100; Sigma-Aldrich, Steinheim, Germany) for 1–3 days at 4 °C. After pre-incubation, nervous tissue was transferred to the primary antibody solution containing 2 % NGS in PBT and incubated for 2–4 days at 4 °C. To label neuropil regions selectively, a monoclonal primary antibody from mouse against synapsin was used in combination with specific additional antibodies and various dyes (for an overview of the antibodies and dyes employed, see Additional file 18: Table S2). After rinsing five times for 10 min with PBT, the brains were incubated with appropriate secondary antibodies and various dyes (Additional file 18: Table S2) diluted in PBT containing 2 % NGS for 1–3 days at 4 °C, followed by 3 to 5 washing steps for 10 min each with PBT. Brains and ganglia were dehydrated in an ascending ethanol series (50 %, 70 %, 90 %, 95 %, 100 %, and 100 % for 2.5 min each) and cleared with methyl salicylate (Merck, Gernsheim, Germany). Finally, they were mounted on coverslips using Permount mounting medium (Fisher Scientific, Pittsburgh, PA, USA) and a stack of two reinforcement rings (Zweckform, Oberlaindern, Germany) as spacers to prevent compression. Brains and ganglia of some of the backfills were not dehydrated and directly mounted in Aqua-Poly/Mount (Polysciences Europe Inc., Eppelheim, Germany).
Antennae and palps of the EF1-B-DsRed, the Orco-Gal4/UAS-tGFP, or Orco-Gal4/UAS-dsRed lines were dissected and fixed overnight at 4 °C in 4 % paraformaldehyde and 10 % methanol in PBT. Afterwards, they were transferred into silicone molds, embedded in tissue-freezing media (Leica, Wetzlar, Germany), and frozen for at least 1 hour at -80 °C, followed by cutting into 50 μm sections at -23 °C on a Cryotome (Cryotome CM 1959, Leica Microsystems, Wetzlar, Germany) resulting in longitudinally halved antennae. The half mounts were collected in a tube and rinsed four times for 20 min each at room temperature in PBT. The samples were pre-incubated with 5 % NGS in PBT overnight at 4 °C followed by incubation with primary antibodies and dyes together with 5 % NGS in PBT overnight. After washing four times for 20 min with PBT, the samples were incubated with appropriate secondary antibodies (Additional file 18: Table S2) overnight at 4 °C. Finally, the antennae were rinsed four times with PBT for 20 min and embedded on coverslips in Aqua-Poly/Mount with one layer of reinforcement rings as spacers.
The specificity of the Orco-antiserum (Moth-R2, kindly provide by Jürgen Krieger) in T. castaneum could be demonstrated by IHC on antennae of animals with RNA interference-mediated knock-down of Orco [115]. To circumvent problems during dsRNA synthesis previously observed with the full length CDS of TcasOrco, we cloned a 476 bp fragment from San Bernadino cDNA containing only a part of CDS and the majority of the 3' untranslated region amplified by Advantage2 Taq Polymerase and primers TcasOrco3UTRrev and TcasOrco3for (see Additional file 16) into PCRII vector (Invitrogen). Using PCR, a bidirectional template was generated followed by dsRNA synthesis with the MEGAscript T7 transcription kit (Ambion, Austin, USA) [196]. Orco dsRNA was injected into pupa of the strain San Bernardino. About 7 days after adult eclosion, the antennae of the treated animals were collected together with antennae of untreated beetles of the black strain, which can be easily discriminated based on the cuticle color, and thus, they served as internal staining controls. A maximal projection of a confocal stack of the Orco-antiserum (Moth-R2) treated antennae shows no detectable antibody staining in RNAi-treated animals (Additional file 2: Figure S2b) in contrast to the black beetle internal control (Additional file 2: Figure S2c).
In vivo backfills of the antenna, single maxillary palps, and whole mouthparts
Cold anesthetized animals were mounted with dental wax (S-U-wax wire, 2.0 mm, hard; Schuler Dental, Ulm, Germany) and modelling clay (Das grosse Dino-Knet-Set; Moses, Verlag GMBH, Kempen, Germany) using a low-temperature soldering iron (Solder-Unit ST 081; Star Tec Products, Bremen, Germany) or with rubber cement (Fixogum, Marabu, Tamm, Germany) with their dorsal side on a microscope slide. The last three segments of the antenna and the most distal segment of the maxillary palp were removed and 4 % neurobiotin in 1 M KCl (Vector Laboratories, Burlingame, UK) for the antenna and Texas Red coupled dextran 50 mg/ml in PBS (3000 MW; Molecular Probes, Invitrogen) for the maxillary palps were used as neuronal tracers. Glass micropipettes were drawn (Model P-97, Sutter Instrument, Novato, USA) from borosilicate glass (inner diameter, 0.75 mm; outer diameter, 1.5 mm; Hilgenberg, Malsfeld, Germany) and broken to a tip diameter matching to the antenna/maxillary palp stump. The dye-filled glass micropipette was put on the antenna/maxillary palp stump for about 4–6 hours in a moist chamber at 4 °C. For the backfills of the whole mouthparts, the maxillary and labial palps were cut and the antennae were protected from unintentional dye-filling by covering them with dental wax (S-U-wax wire, 2.0 mm, hard). A crystal of biotin-conjugated dextran (3000 MW; Molecular Probes, Invitrogen) was placed onto the prepared mouthparts, covered with a drop of distilled water, and stored for about 4 h in a moist chamber at 4 °C. Brains and ganglia were dissected, fixed, washed, and stained as described above. Neurobiotin was visualized with Cy3 conjugated streptavidin (Dianova, Hamburg, Germany) diluted 1/200 in PBT (0.3 % TrX). The staining solution contained in addition Alexa Fluor 488-coupled phalloidin (1/200), DAPI (1/20,000) and 2 % NGS. The incubation time was 2–3 days at 4 °C. Biotin-coupled dextran was visualized with Alexa Fluor 488-coupled streptavidin (Molecular Probes, Invitrogen) diluted 1/200 in PBT (0.3 % TrX and 2 % NGS) and applied together with synapsin (1/300) for 2–3 days at 4 °C.
In vivo dye injection into the antennal lobes
Cold anesthetized animals with fluorescent labeled ALs (partial Orco-Gal4/UAS-DsRed) were mounted with their ventral side pointing upside down with dental wax on a microscope slide. The pronotum and the head capsule were opened using a piece of a razor blade held by a blade breaker, with two parallel longitudinal cuts along the compound eyes. The cuticle, fat tissue, and tracheae were removed. Afterwards, the head capsule and pronotum were covered with ringer solution [197]. A tungsten needle was sharpened in 2 M KOH with 5–8 volts as described in [198], followed by coating with Texas Red conjugated dextran (3000 MW; Molecular Probes, Invitrogen) dissolved in NGS and air-dried. The injection of dye into the DsRed-labeled AL was performed manually under a fluorescence stereomicroscope (SteREO Lumar.V12, Carl Zeiss MicroImaging, Jena, Germany) by careful perforation. The treated animals were kept in a moist chamber for about 1 h at room temperature to let the dye diffuse. Afterwards, the brains were dissected, fixed, washed, and pre-incubated with NGS as described previously and afterwards incubated with Alexa Fluor 488-coupled phalloidin (1/200), DAPI (1/20,000), and 2 % NGS for 2 days at 4 °C. Subsequently, the brains were washed, dehydrated, cleared, and mounted in Permount as described above.
Microscopic image acquisition, processing, and analysis
The fluorescent-labeled microscopic samples were scanned with CLSM (TCS SP5, Leica Microsystems) at 1024 × 1024 or 2048 × 2048 pixel resolution, a scanning speed between 100 and 200 Hz, a pinhole of size 1 airy, a line average of 2–4, and a step size between 0.5 and 2.5 μm. Confocal images and image stacks were analyzed with the Amira 5.3.3 graphics software (FEI, Hillsboro, OR, USA). The final image processing and figure arrangements were processed using Corel Draw X3 (Corel, Ottawa, Ontario, Canada), Adobe Photoshop CS3 (Adobe Systems, San Jose, CA, USA), or Inkscape [199].
The number of CSNs per sensillum was determined based on high-resolution CLSM stacks taken from antennae of the EF1-B-DsRed line after antibody enhancement of the DsRed reporter signal in combination with Orco antibody staining. To determine the number of CSNs and Orco-immunoreactive OSNs, we traced the stained dendrites of the CSNs to their associated soma of several sensilla and calculated their average number (Additional file 3: Figure S3i; Fig. 2h).
AL glomeruli were separately labeled in the AMIRA Segmentation Editor and 3D reconstructed [102, 200] based on CLSM stacks of brains labeled with synapsin and TKRP antibodies of five male and five female A7 beetles (one AL from a random hemisphere for each brain). To optimize data quality, the CLSM stacks were previously deconvoluted in AMIRA using the blind method with initial estimation set to input data, with a border width of 10, 10, and 10, and an iteration of ten cycles.
KCs were identified based on their position, size, and density in DAPI stainings [103]. The total volumes of the whole CAs (13 CAs of seven A7 males), as well as the volumes of three randomly assigned clusters of 20 KCs per CA, were measured using 3D reconstruction. For the segmentation and reconstruction details, we refer to [201, 202]. Briefly, different layers of a structure were labeled in the Segmentation Editor and wrapped. Volumes of reconstructed structures were taken from Material Statistics. Based on the ratios between whole CA volume and the volumes of the three clusters of 20 KCs, the total number of KCs per CA was interpolated. In addition, we counted the KCs by an independent method using MorphoGraphX (www.MorphoGraphX.org). The CLSM stacks were processed with the arithmetic tool of AMIRA to mask the CAs and consequently to remove the remaining materials. The resulting stacks were converted to TIFF files with FIJI [203] by preserving the image properties. These files were analyzed with the “Local Maxima” tool of MorphoGraphX [111] with the following parameter X-/Y-/Z-radius = 1 μm, Start Label -1, and Min Color 1. Because of the inhomogeneous intensity distribution within some CLSM stacks, only nine CAs from five specimens were analyzed automatically.
Statistical analysis included determination of mean values, standard deviations, and independent/unpaired two-tailed Student’s t-tests, which were performed in Excel XP (Microsoft, Redmond, WA, USA). Bar charts were created in Excel XP and imported and revised in Corel Draw.
RNA isolation and sequencing
Total RNA of about 1000 antennae, 150 mouthparts (a piece of the head capsule anterior of the antennae), 600 legs, 50 heads (without antennae but including mouthparts), and 20 remaining bodies of sex-separated adults was isolated using the ZR Tissue and Insect RNA Micro Prep Kit (Zymo Research, Irvine, CA, USA), following the manufacturer´s protocol. The library preparations for RNA-Seq were performed using the TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA) and cDNA libraries were amplified and sequenced using the cBot and HiSeq2000 from Illumina (paired end; 2 × 100 bp). Biological triplicates of female and male antennae were sequenced. For the other tissues, one male and one female sample were used and one additional mouthpart sample was obtained from unsexed beetles. The number of biological replicates was chosen based on previous publications using similar approaches [128, 170, 171]. Each biological replicate came from independently prepared and processed tissue. No technical replicates were performed. For details, see Dippel et al. [89].
Re-annotation of olfactory genes
For manual inspection the obtained reads were mapped against the T. castaneum 4.0 genome using BLAT [204] and a genome browser was set up (http://bioinf.uni-greifswald.de/tcas/). In a genome independent approach, a de novo assembly was built with Trinity (release 2013_08_14) [205] as described in [89]. The previously published OR [115], GR [81], IR [114], and SNMP [116, 117] sequences were used for further analysis. To identify the potential ODEs, the official (OGS3) [81–83], the preliminary AU2 and AU3, and the NCBI [206, 207] gene sets were used and a protein functional analysis was conducted using InterProScan [208]. All genes belonging to a protein family containing known ODEs in other insect species were collected (namely, ALDH, ALOX, CES, EH, GST, and CYP [4]). The redundant genes were removed and the sequences were reviewed. The identified GSTs and CYPs were collated to already published sequences [148, 150] and the names were adapted. For all other candidates, a genome-based name was built reflecting the protein family and the chromosomal localization (e.g., CES2D is the fourth CES on the second chromosome). The genes supposed to be involved in olfactory transduction of D. melanogaster were taken from [46], the corresponding sequences were downloaded from the Flybase [209], and the T. castaneum orthologs were identified by pBLAST embedded in the genome browser (http://bioinf.uni-greifswald.de/tcas/).
The revision of the olfactory genes was performed in an iterative process based on sequence comparison with the de novo assembly and the RNA-seq based gene annotations (AU3), a conserved domain search [210], and manual inspection of the aligned reads in the genome browser. For discrepancies, the gene models were manually curated. Finally, the chromosomal localization of the olfactory genes was determined by pBLAST against the genome assembly Tcas4.0. The ODE candidates were searched for signal peptides using the SignalP4.1 server [211]. The sequences and read numbers are summarized in Additional file 8: Table S1. The complete dataset including all relevant parameters has been deposited in the National Center for Biotechnology Information (NCBI) database repository ‘Gene Expression Omnibus’ (GEO accession number: GSE63162).
Tribolium castaneum expression profiling
The olfactory genes were identified in the AU3 gene set by pBLAST and the corresponding gene models were replaced with the re-annotated candidate sequences. The resulting enhanced AU3 gene set was used to map the RNAseq data with Bowtie2 [212] using the very sensitive presetting. The mapped reads were counted with Samtools [213] and normalized as RPKM values. The RPKMs were visualized (matrix2png interface, version 1.2.1; [214]) and the figures were arranged in Inkscape [199]. One male antennae sample was excluded from the subsequent analysis due to massive differences with the other five antennae samples in the principal component as well as cluster analysis. Because no significant differences for the genes of interest were observed between male and female reads, the sequenced tissues were considered as biological replicates. Statistical analysis was performed in R [215] with the DESeq package (version 1.12.0) [216] from Bioconductor [217]. p values were calculated with a negative binomial test using raw read counts and adjusted for multiple testing with the Benjamini–Hochberg method. All tissues were compared to body as reference. Significant differentially expressed genes (false discovery rate <0.05) are marked with asterisks in the heat maps. However, since only two body samples served as controls, all conclusions from this should be treated as preliminary. Genes with an RPKM ≥ 0.5 were considered as specifically expressed in the tissue. The tissue comparison was visualized as Venn diagrams (e.g., Additional file 9: Figure S5a) (http://bioinformatics.psb.ugent.be/webtools/Venn/). For details, see Dippel et al. [89].
Phylogenetic analysis and interspecies comparison
We compared the T. castaneum IR, GR, and OR sequences independent from each other on a protein level with data from D. melanogaster [46, 209] and An. gambiae [128, 218]. The sequences were aligned using MAFFT (v7.040b [219]) (--genafpair --maxiterate 1000 --bl 62 --op 1.53 --ep 0.123) and the phylogeny was calculated using RAxML (version 7.8.6 [220]), with the LG substitution model and GAMMA correction. The robustness of the tree topology was evaluated by 100 rapid bootstrap replications. The relative expression levels were calculated as log2-fold changes of antenna/body and palp (mouthpart)/body as described in [89]. The D. melanogaster data set was downloaded from the EMBL gene expression atlas [221], originally published in [170], and the An. gambiae data were obtained from [128]. The phylogenetic tree was visualized by iTOL [222] and descriptions were added using Inkscape [199].
Abbreviations
- AL:
-
Antennal lobe
- ALDH:
-
Aldehyde dehydrogenase
- ALOX:
-
Aldehyde oxidase
- ALT:
-
Antennal lobe tract
- AMMC:
-
antennal mechanosensory and motor center
- AN:
-
antennal nerve
- CA:
-
Calyx of the MB
- CES:
-
Carboxylesterase
- CLSM:
-
Confocal laser-scanning microscopy
- CSN:
-
Chemosensory neuron
- CSP:
-
Chemosensory protein
- cSTri:
-
Chemosensilla trichoidea
- CYP:
-
Cytochrome P450
- EH:
-
Epoxide hydrolase
- GABA:
-
Gamma-aminobutyric acid
- GL:
-
Antennal lobe glomeruli
- GNG:
-
Gnathal ganglia
- GOC:
-
Gnathal olfactory center
- GR:
-
Gustatory receptor
- GSN:
-
Gustatory sensory neuron
- GST:
-
Glutathione S-transferase
- IHC:
-
Immunohistochemistry
- IR:
-
Ionotropic glutamate-like receptor
- KC:
-
Kenyon cells
- lALT:
-
Lateral antennal lobe tract
- LG:
-
Lobus glomerulatus
- LH:
-
Lateral horn (lateral protocerebrum)
- LN:
-
Local interneuron
- mALT:
-
Mediolateral lobe tract
- MB:
-
Mushroom body
- mlALT:
-
Mediolateral antennal lobe tract
- mSTri:
-
Mechanosensilla trichoidea
- NA:
-
Neuro-axis
- NGS:
-
normal goat serum
- OBP:
-
Odorant binding protein
- ODE:
-
Odorant degrading enzyme
- OR:
-
Odorant receptor
- Orco:
-
Odorant receptor co-receptor
- ORF:
-
Open reading frame
- OSN:
-
Olfactory sensory neuron
- PBS:
-
Phosphate buffered saline
- PGC:
-
Primary gustatory center
- PN:
-
Projection neuron
- RPKM:
-
Reads per kilobase per million
- SBas:
-
Sensilla basiconica
- SCam:
-
Sensilla campaniformes
- SCha:
-
Sensilla chaetica
- SCoe:
-
Sensilla coeloconica
- SD:
-
standard deviation
- SEM:
-
Scanning electron microscopy
- SNMP:
-
Sensory neuron membrane protein
- SpaB:
-
Spaculate bristle
- TKRP:
-
Tachykinin-related peptide
References
de Bruyne M, Smart R, Zammit E, Warr CG. Functional and molecular evolution of olfactory neurons and receptors for aliphatic esters across the Drosophila genus. J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2010;196:97–109.
Dahanukar A, Hallem EA, Carlson JR. Insect chemoreception. Curr Opin Neurobiol. 2005;15:423–30.
Tegoni M, Campanacci V, Cambillau C. Structural aspects of sexual attraction and chemical communication in insects. Trends Biochem Sci. 2004;29:257–64.
Leal WS. Odorant reception in insects: roles of receptors, binding proteins, and degrading enzymes. Annu Rev Entomol. 2013;58:373–91.
Visser JH. Host odor perception in phytophagous insects. Annu Rev Entomol. 1986;31:121–44.
Whiteman NK, Pierce NE. Delicious poison: genetics of Drosophila host plant preference. Trends Ecol Evol. 2008;23:473–8.
Liu M, Yu H, Li G. Oviposition deterrents from eggs of the cotton bollworm, Helicoverpa armigera (Lepidoptera: Noctuidae): chemical identification and analysis by electroantennogram. J Insect Physiol. 2008;54:656–62.
Sun Y-L, Huang L-Q, Pelosi P, Wang C-Z. Expression in antennae and reproductive organs suggests a dual role of an odorant-binding protein in two sibling Helicoverpa species. PLoS One. 2012;7, e30040.
Dicke M. Behavioural and community ecology of plants that cry for help. Plant Cell Environ. 2009;32:654–65.
Stensmyr MC, Dweck HKM, Farhan A, Ibba I, Strutz A, Mukunda L, et al. A conserved dedicated olfactory circuit for detecting harmful microbes in Drosophila. Cell. 2012;151:1345–57.
Laska M, Galizia CG, Giurfa M, Menzel R. Olfactory discrimination ability and odor structure-activity relationships in honeybees. Chem Senses. 1999;24:429–38.
Yang C-H, Belawat P, Hafen E, Jan LY, Jan Y-N. Drosophila egg-laying site selection as a system to study simple decision-making processes. Science. 2008;319:1679–83.
Weiss LA, Dahanukar A, Kwon JY, Banerjee D, Carlson JR. The molecular and cellular basis of bitter taste in Drosophila. Neuron. 2011;69:258–72.
Linz J, Baschwitz A, Strutz A, Dweck HKM, Sachse S, Hansson BS, et al. Host plant-driven sensory specialization in Drosophila erecta. Proc R Soc Lond B Biol Sci. 2013;280:20130626.
Paczkowski S, Paczkowska M, Dippel S, Flematti G, Schütz S. Volatile combustion products of wood attract Acanthocnemus nigricans (Coleoptera: Acanthocnemidae). J Insect Behav. 2014;27:228–38.
Steinbrecht RA. Structure and function of insect olfactory sensilla. Ciba Found Symp. 1996;200:158–77.
de Bruyne M, Warr CG. Molecular and cellular organization of insect chemosensory neurons. Bioessays. 2006;28:23–34.
Pelosi P, Iovinella I, Felicioli A, Dani FR. Soluble proteins of chemical communication: an overview across arthropods. Front Physiol. 2014;5:320.
Benton R, Vannice KS, Gomez-Diaz C, Vosshall LB. Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell. 2009;136:149–62.
Rytz R, Croset V, Benton R. Ionotropic receptors (IRs): chemosensory ionotropic glutamate receptors in Drosophila and beyond. Insect Biochem Mol Biol. 2013;43:888–97.
Montell C. A taste of the Drosophila gustatory receptors. Curr Opin Neurobiol. 2009;19:345–53.
Sato K, Pellegrino M, Nakagawa T, Nakagawa T, Vosshall LB, Touhara K. Insect olfactory receptors are heteromeric ligand-gated ion channels. Nature. 2008;452:1002–6.
Wicher D, Schäfer R, Bauernfeind R, Stensmyr MC, Heller R, Heinemann SH, et al. Drosophila odorant receptors are both ligand-gated and cyclic-nucleotide-activated cation channels. Nature. 2008;452:1007–11.
Missbach C, Dweck HK, Vogel H, Vilcinskas A, Stensmyr MC, Hansson BS, et al. Evolution of insect olfactory receptors. eLife. 2014;3, e02115.
Sánchez-Gracia A, Vieira FG, Almeida FC, Rozas J. Comparative Genomics of the Major Chemosensory Gene Families in Arthropods. eLS. 2011. doi:10.1002/9780470015902.a0022848.
Cao TNP. Genome annotation and evolution of chemosensory receptors in spider mites. Dissertation. Ghent University; 2014.
Chipman AD, Ferrier DEK, Brena C, Qu J, Hughes DST, Schröder R, et al. The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima. PLoS Biol. 2014;12, e1002005.
Chyb S. Drosophila gustatory receptors: from gene identification to functional expression. J Insect Physiol. 2004;50:469–77.
Freeman EG, Wisotsky Z, Dahanukar A. Detection of sweet tastants by a conserved group of insect gustatory receptors. Proc Natl Acad Sci. 2014;111:1598–603.
Lee Y, Moon SJ, Montell C. Multiple gustatory receptors required for the caffeine response in Drosophila. Proc Natl Acad Sci. 2009;106:4495–500.
Jiao Y, Moon SJ, Wang X, Ren Q, Montell C. Gr64f is required in combination with other gustatory receptors for sugar detection in Drosophila. Curr Biol. 2008;18:1797–801.
Ishimoto H, Takahashi K, Ueda R, Tanimura T. G-protein gamma subunit 1 is required for sugar reception in Drosophila. EMBO J. 2005;24:3259–65.
Sato K, Tanaka K, Touhara K. Sugar-regulated cation channel formed by an insect gustatory receptor. Proc Natl Acad Sci. 2011;108:11680–5.
Smart R, Kiely A, Beale M, Vargas E, Carraher C, Kralicek AV, et al. Drosophila odorant receptors are novel seven transmembrane domain proteins that can signal independently of heterotrimeric G proteins. Insect Biochem Mol Biol. 2008;38:770–80.
Mukunda L, Lavista-Llanos S, Hansson BS, Wicher D. Dimerisation of the Drosophila odorant coreceptor Orco. Front Cell Neurosci. 2014;8:261.
German PF, van der Poel S, Carraher C, Kralicek AV, Newcomb RD. Insights into subunit interactions within the insect olfactory receptor complex using FRET. Insect Biochem Mol Biol. 2013;43:138–45.
Getahun MN, Olsson SB, Lavista-Llanos S, Hansson BS, Wicher D. Insect odorant response sensitivity is tuned by metabotropically autoregulated olfactory receptors. PLoS One. 2013;8, e58889.
Stengl M, Funk NW. The role of the coreceptor Orco in insect olfactory transduction. J Comp Physiol A. 2013;199:897–909.
Nolte A, Funk NW, Mukunda L, Gawalek P, Werckenthin A, Hansson BS, et al. In situ tip-recordings found no evidence for an Orco-based ionotropic mechanism of pheromone-transduction in Manduca sexta. PLoS One. 2013;8, e62648.
Martin F, Alcorta E. Regulation of olfactory transduction in the Orco channel. Front Cell Neurosci. 2011;5:21.
Riesgo-Escovar J, Raha D, Carlson JR. Requirement for a phospholipase C in odor response: overlap between olfaction and vision in Drosophila. Proc Natl Acad Sci. 1995;92:2864–8.
Miura N, Atsumi S, Tabunoki H, Sato R. Expression and localization of three G protein α subunits, Go, Gq, and Gs, in adult antennae of the silkmoth (Bombyx mori). J Comp Neurol. 2005;485:143–52.
Chatterjee A, Roman G, Hardin PE. Go contributes to olfactory reception in Drosophila melanogaster. BMC Physiol. 2009;9:22.
Deng Y, Zhang W, Farhat K, Oberland S, Gisselmann G, Neuhaus EM. The stimulatory gαs protein is involved in olfactory signal transduction in Drosophila. PLoS One. 2011;6, e18605.
Sargsyan V, Getahun MN, Llanos SL, Olsson SB, Hansson BS, Wicher D. Phosphorylation via PKC regulates the function of the Drosophila odorant co-receptor. Front Cell Neurosci. 2011;5:5.
Martin F, Boto T, Gomez-Diaz C, Alcorta E. Elements of olfactory reception in adult Drosophila melanogaster. Anat Rec. 2013;296:1477–88.
Ignatious Raja JS, Katanayeva N, Katanaev VL, Galizia CG. Role of Go/i subgroup of G proteins in olfactory signaling of Drosophila melanogaster. Eur J Neurosci. 2014;39:1245–55.
Jin X, Ha TS, Smith DP. SNMP is a signaling component required for pheromone sensitivity in Drosophila. Proc Natl Acad Sci. 2008;105:10996–1001.
Benton R, Vannice KS, Vosshall LB. An essential role for a CD36-related receptor in pheromone detection in Drosophila. Nature. 2007;450:289–93.
Li Z, Ni JD, Huang J, Montell C. Requirement for Drosophila SNMP1 for rapid activation and termination of pheromone-induced activity. PLoS Genet. 2014;10, e1004600.
Maïbèche-Coisne M, Nikonov AA, Ishida Y, Jacquin-Joly E, Leal WS. Pheromone anosmia in a scarab beetle induced by in vivo inhibition of a pheromone-degrading enzyme. Proc Natl Acad Sci. 2004;101:11459–64.
Ishida Y, Leal WS. Rapid inactivation of a moth pheromone. Proc Natl Acad Sci. 2005;102:14075–9.
Durand N, Carot-Sans G, Bozzolan F, Rosell G, Siaussat D, Debernard S, et al. Degradation of pheromone and plant volatile components by a same odorant-degrading enzyme in the cotton leafworm, Spodoptera littoralis. PLoS ONE. 2011;6, e29147.
Chertemps T, François A, Durand N, Rosell G, Dekker T, Lucas P, et al. A carboxylesterase, Esterase-6, modulates sensory physiological and behavioral response dynamics to pheromone in Drosophila. BMC Biol. 2012;10:56.
Younus F, Chertemps T, Pearce SL, Pandey G, Bozzolan F, Coppin CW, et al. Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogaster. Insect Biochem Mol Biol. 2014;53:30–43.
Scott K, Brady Jr R, Cravchik A, Morozov P, Rzhetsky A, Zuker C, et al. A chemosensory gene family encoding candidate gustatory and olfactory receptors in Drosophila. Cell. 2001;104:661–73.
Schachtner J, Schmidt M, Homberg U. Organization and evolutionary trends of primary olfactory brain centers in Tetraconata (Crustacea + Hexapoda). Arthropod Struct Dev. 2005;34:257–99.
Vosshall LB. Olfaction in Drosophila. Curr Opin Neurobiol. 2000;10:498–503.
Stocker RF. Drosophila as a focus in olfactory research: mapping of olfactory sensilla by fine structure, odor specificity, odorant receptor expression, and central connectivity. Microsc Res Tech. 2001;55:284–96.
Keller A, Vosshall LB. Decoding olfaction in Drosophila. Curr Opin Neurobiol. 2003;13:103–10.
Jefferis GSXE. Insect olfaction: a map of smell in the brain. Curr Biol. 2005;15:R668–70.
Smith CUM. Olfaction. In: Biology of sensory systems. John Wiley & Sons; 2009. doi: 10.1002/9780470694374.ch14.
Couto A, Alenius M, Dickson BJ. Molecular, anatomical, and functional organization of the Drosophila olfactory system. Curr Biol. 2005;15:1535–47.
Ernst DKD, Boeckh J, Boeckh V. A neuroanatomical study on the organization of the central antennal pathways in insects. Cell Tissue Res. 1977;176:285–308.
Ignell R, Anton S, Hansson BS. The maxillary palp sensory pathway of Orthoptera. Arthropod Struct Dev. 2000;29:295–305.
Hofer S, Dircksen H, Tollbäck P, Homberg U. Novel insect orcokinins: characterization and neuronal distribution in the brains of selected dicondylian insects. J Comp Neurol. 2005;490:57–71.
Christensen TA, Waldrop BR, Harrow ID, Hildebrand JG. Local interneurons and information processing in the olfactory glomeruli of the moth Manduca sexta. J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 1993;173:385–99.
Chou Y-H, Spletter ML, Yaksi E, Leong JCS, Wilson RI, Luo L. Diversity and wiring variability of olfactory local interneurons in the Drosophila antennal lobe. Nat Neurosci. 2010;13:439–49.
Seki Y, Rybak J, Wicher D, Sachse S, Hansson BS. Physiological and morphological characterization of local interneurons in the Drosophila antennal lobe. J Neurophysiol. 2010;104:1007–19.
Galizia CG, Rössler W. Parallel olfactory systems in insects: anatomy and function. Annu Rev Entomol. 2010;55:399–420.
Wang Y, Chiang A-S, Xia S, Kitamoto T, Tully T, Zhong Y. Blockade of neurotransmission in Drosophila mushroom bodies impairs odor attraction, but not repulsion. Curr Biol. 2003;13:1900–4.
Jefferis GSXE, Potter CJ, Chan AM, Marin EC, Rohlfing T, Maurer Jr CR, et al. Comprehensive maps of Drosophila higher olfactory centers: spatially segregated fruit and pheromone representation. Cell. 2007;128:1187–203.
Yamagata N, Nishino H, Mizunami M. Neural pathways for the processing of alarm pheromone in the ant brain. J Comp Neurol. 2007;505:424–42.
Strutz A, Soelter J, Baschwitz A, Farhan A, Grabe V, Rybak J, et al. Decoding odor quality and intensity in the Drosophila brain. eLife. 2014;3, e04147.
de Belle J, Heisenberg M. Associative odor learning in Drosophila abolished by chemical ablation of mushroom bodies. Science. 1994;263:692–5.
Connolly JB, Roberts IJ, Armstrong JD, Kaiser K, Forte M, Tully T, et al. Associative learning disrupted by impaired Gs signaling in Drosophila mushroom bodies. Science. 1996;274:2104–7.
Heimbeck G, Bugnon V, Gendre N, Keller A, Stocker RF. A central neural circuit for experience-independent olfactory and courtship behavior in Drosophila melanogaster. Proc Natl Acad Sci. 2001;98:15336–41.
Morris SC, Grimaldi D, Engel MS. Evolution of the Insects. Cambridge, New York, Melbourne: Cambridge University Press; 2005. http://www.cambridge.org/catalogue/catalogue.asp?isbn=9780521821490&ss=fro.
Hunt T, Bergsten J, Levkanicova Z, Papadopoulou A, John OS, Wild R, et al. A comprehensive phylogeny of beetles reveals the evolutionary origins of a superradiation. Science. 2007;318:1913–6.
Brown SJ, Shippy TD, Miller S, Bolognesi R, Beeman RW, Lorenzen MD, et al. The red flour beetle, Tribolium castaneum (Coleoptera): a model for studies of development and pest biology. Cold Spring Harb Protoc. 2009;2009:pdb.emo126.
Richards S, Gibbs RA, Weinstock GM, Brown SJ, Denell R, Beeman RW, et al. The genome of the model beetle and pest Tribolium castaneum. Nature. 2008;452:949–55.
Wang L, Wang S, Li Y, Paradesi MSR, Brown SJ. BeetleBase: the model organism database for Tribolium castaneum. Nucleic Acids Res. 2007;35:D476–9.
Kim HS, Murphy T, Xia J, Caragea D, Park Y, Beeman RW, et al. BeetleBase in 2010: revisions to provide comprehensive genomic information for Tribolium castaneum. Nucleic Acids Res. 2010;38:D437–42.
Bucher G, Scholten J, Klingler M. Parental RNAi in Tribolium (Coleoptera). Curr Biol. 2002;12:R85–6.
Tomoyasu Y, Denell RE. Larval RNAi in Tribolium (Coleoptera) for analyzing adult development. Dev Genes Evol. 2004;214:575–8.
Lorenzen MD, Berghammer, Brown SJ, Denell RE, Klingler M, Beeman RW. piggyBac-mediated germline transformation in the beetle Tribolium castaneum. Insect Mol Biol. 2003;12:433–40.
Schinko JB, Weber M, Viktorinova I, Kiupakis A, Averof M, Klingler M, et al. Functionality of the GAL4/UAS system in Tribolium requires the use of endogenous core promoters. BMC Dev Biol. 2010;10:53.
Schinko JB, Hillebrand K, Bucher G. Heat shock-mediated misexpression of genes in the beetle Tribolium castaneum. Dev Genes Evol. 2012;222:287–98.
Dippel S, Oberhofer G, Kahnt J, Gerischer L, Opitz L, Schachtner J, et al. Tissue-specific transcriptomics, chromosomal localization, and phylogeny of chemosensory and odorant binding proteins from the red flour beetle Tribolium castaneum reveal subgroup specificities for olfaction or more general functions. BMC Genomics. 2014;15:1141.
Roth LM, Willis ER. Hygroreceptors in adults of Tribolium (Coleoptera, Tenebrionidae). J Exp Zool. 1951;116:527–70.
Vogt RG, Riddiford LM. Pheromone binding and inactivation by moth antennae. Nature. 1981;293:161–3.
Keil TA. Fine structure of the pheromone-sensitive sensilla on the antenna of the hawkmoth, Manduca sexta. Tissue Cell. 1989;21:139–51.
Almaas TJ, Mustaparta H. Pheromone reception in tobacco budworm moth, Heliothis virescens. J Chem Ecol. 1990;16:1331–47.
Hill SR, Hansson BS, Ignell R. Characterization of antennal trichoid sensilla from female southern house mosquito, Culex quinquefasciatus Say. Chem Senses. 2009;34:231–52.
Alabi T, Marion-Poll F, Danho M, Mazzucchelli GD, De Pauw E, Haubruge E, et al. Identification of taste receptors and proteomic characterization of the antenna and legs of Tribolium brevicornis, a stored food product pest. Insect Mol Biol. 2014;23:1–12.
Harbach RE, Larsen JR. Fine structure of antennal sensilla of the adult mealworm beetle, Tenebrio molitor L. (Coleoptera: Tenebrionidae). Int J Insect Morphol Embryol. 1977;6:41–60.
Shanbhag SR, Müller B, Steinbrecht RA. Atlas of olfactory organs of Drosophila melanogaster: 1. types, external organization, innervation and distribution of olfactory sensilla. Int J Insect Morphol Embryol. 1999;28:377–97.
Solari P, Corda V, Sollai G, Kreissl S, Galizia CG, Crnjar R. Morphological characterization of the antennal lobes in the Mediterranean fruit fly Ceratitis capitata. J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2016;202:131–46.
Kent KS, Harrow ID, Quartararo P, Hildebrand DJG. An accessory olfactory pathway in Lepidoptera: the labial pit organ and its central projections in Manduca sexta and certain other sphinx moths and silk moths. Cell Tissue Res. 1986;245:237–45.
Bräunig P, Pflüger H-J, Hustert R. The specificity of central nervous projections of locust mechanoreceptors. J Comp Neurol. 1983;218:197–207.
Nishino H, Nishikawa M, Yokohari F, Mizunami M. Dual, multilayered somatosensory maps formed by antennal tactile and contact chemosensory afferents in an insect brain. J Comp Neurol. 2005;493:291–308.
Dreyer D, Vitt H, Dippel S, Goetz B, El Jundi B, Kollmann M, et al. 3D standard brain of the red flour beetle Tribolium castaneum: a tool to study metamorphic development and adult plasticity. Front Syst Neurosci. 2010;4:3.
Binzer M, Heuer CM, Kollmann M, Kahnt J, Hauser F, Grimmelikhuijzen CJP, et al. Neuropeptidome of Tribolium castaneum antennal lobes and mushroom bodies. J Comp Neurol. 2014;522:337–57.
Miyazaki T, Ito K. Neural architecture of the primary gustatory center of Drosophila melanogaster visualized with GAL4 and LexA enhancer-trap systems. J Comp Neurol. 2010;518:4147–81.
Farris SM. Tritocerebral tract input to the insect mushroom bodies. Arthropod Struct Dev. 2008;37:492–503.
Tanaka NK, Endo K, Ito K. Organization of antennal lobe-associated neurons in adult Drosophila melanogaster brain. J Comp Neurol. 2012;520:4067–130.
Frambach I, Schürmann FW. Separate distribution of deutocerebral projection neurons in the mushroom bodies of the cricket brain. Acta Biol Hung. 2004;55:21–9.
Groh C, Lu Z, Meinertzhagen IA, Rössler W. Age-related plasticity in the synaptic ultrastructure of neurons in the mushroom body calyx of the adult honeybee Apis mellifera. J Comp Neurol. 2012;520:3509–27.
Leiss F, Groh C, Butcher NJ, Meinertzhagen IA, Tavosanis G. Synaptic organization in the adult Drosophila mushroom body calyx. J Comp Neurol. 2009;517:808–24.
Caron SJC, Ruta V, Abbott LF, Axel R. Random convergence of olfactory inputs in the Drosophila mushroom body. Nature. 2013;497:113–7.
de Reuille PB, Routier-Kierzkowska A-L, Kierzkowski D, Bassel GW, Schüpbach T, Tauriello G, et al. MorphoGraphX: a platform for quantifying morphogenesis in 4D. eLife. 2015;4, e05864.
Foret S, Maleszka R. Function and evolution of a gene family encoding odorant binding-like proteins in a social insect, the honey bee (Apis mellifera). Genome Res. 2006;16:1404–13.
Forêt S, Wanner KW, Maleszka R. Chemosensory proteins in the honey bee: insights from the annotated genome, comparative analyses and expressional profiling. Insect Biochem Mol Biol. 2007;37:19–28.
Croset V, Rytz R, Cummins SF, Budd A, Brawand D, Kaessmann H, et al. Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010;6, e1001064.
Engsontia P, Sanderson AP, Cobb M, Walden KKO, Robertson HM, Brown S. The red flour beetle’s large nose: an expanded odorant receptor gene family in Tribolium castaneum. Insect Biochem Mol Biol. 2008;38:387–97.
Nichols Z, Vogt RG. The SNMP/CD36 gene family in Diptera, Hymenoptera and Coleoptera: Drosophila melanogaster, D. pseudoobscura, Anopheles gambiae, Aedes aegypti, Apis mellifera, and Tribolium castaneum. Insect Biochem Mol Biol. 2008;38:398–415.
Vogt RG, Miller NE, Litvack R, Fandino RA, Sparks J, Staples J, et al. The insect SNMP gene family. Insect Biochem Mol Biol. 2009;39:448–56.
GEO accession: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63162.
Enjin A, Zaharieva EE, Frank DD, Mansourian S, Suh GSB, Gallio M, et al. Humidity sensing in Drosophila. Curr Biol. 2016;26:1352–8.
Zhang YV, Ni J, Montell C. The molecular basis for attractive salt-taste coding in Drosophila. Science. 2013;340:1334–8.
Robertson HM, Kent LB. Evolution of the gene lineage encoding the carbon dioxide receptor in insects. J Insect Sci. 2009;9:19.
Kent LB, Robertson HM. Evolution of the sugar receptors in insects. BMC Evol Biol. 2009;9:41.
Ni L, Bronk P, Chang EC, Lowell AM, Flam JO, Panzano VC, et al. A gustatory receptor paralogue controls rapid warmth avoidance in Drosophila. Nature. 2013;500:580–4.
Robertson HM, Wanner KW. The chemoreceptor superfamily in the honey bee, Apis mellifera: expansion of the odorant, but not gustatory, receptor family. Genome Res. 2006;16:1395–403.
Wanner KW, Robertson HM. The gustatory receptor family in the silkworm moth Bombyx mori is characterized by a large expansion of a single lineage of putative bitter receptors. Insect Mol Biol. 2008;17:621–9.
Jones WD, Cayirlioglu P, Kadow IG, Vosshall LB. Two chemosensory receptors together mediate carbon dioxide detection in Drosophila. Nature. 2007;445:86–90.
Kwon JY, Dahanukar A, Weiss LA, Carlson JR. The molecular basis of CO2 reception in Drosophila. Proc Natl Acad Sci. 2007;104:3574–8.
Pitts RJ, Rinker DC, Jones PL, Rokas A, Zwiebel LJ. Transcriptome profiling of chemosensory appendages in the malaria vector Anopheles gambiae reveals tissue- and sex-specific signatures of odor coding. BMC Genomics. 2011;12:271.
Lu T, Qiu YT, Wang G, Kwon JY, Rutzler M, Kwon H-W, et al. Odor coding in the maxillary palp of the malaria vector mosquito Anopheles gambiae. Curr Biol. 2007;17:1533–44.
Hartl M, Loschek LF, Stephan D, Siju KP, Knappmeyer C, Kadow ICG. A new Prospero and microRNA-279 pathway restricts CO2 receptor neuron formation. J Neurosci. 2011;31:15660–73.
Ramaswamy SB. Behavioural responses of Heliothis virescens (Lepidoptera: Noctuidae) to stimulation with sugars. J Insect Physiol. 1987;33:755–60.
de Brito Sanchez MG, Giurfa M, de Paula Mota TR, Gauthier M. Electrophysiological and behavioural characterization of gustatory responses to antennal ‘bitter’ taste in honeybees. Eur J Neurosci. 2005;22:3161–70.
Popescu A, Couton L, Almaas T-J, Rospars J-P, Wright GA, Marion-Poll F, et al. Function and central projections of gustatory receptor neurons on the antenna of the noctuid moth Spodoptera littoralis. J Comp Physiol A. 2013;199:403–16.
Dunipace L, Meister S, McNealy C, Amrein H. Spatially restricted expression of candidate taste receptors in the Drosophila gustatory system. Curr Biol. 2001;11:822–35.
Jacquin-Joly E, Legeai F, Montagné N, Monsempes C, François M-C, Poulain J, et al. Candidate chemosensory genes in female antennae of the noctuid moth Spodoptera littoralis. Int J Biol Sci. 2012;8:1036–50.
Andersson MN, Grosse-Wilde E, Keeling CI, Bengtsson JM, Yuen MM, Li M, et al. Antennal transcriptome analysis of the chemosensory gene families in the tree killing bark beetles, Ips typographus and Dendroctonus ponderosae (Coleoptera: Curculionidae: Scolytinae). BMC Genomics. 2013;14:198.
Vosshall LB, Stocker RF. Molecular architecture of smell and taste in Drosophila. Annu Rev Neurosci. 2007;30:505–33.
Sparks JT, Bohbot JD, Dickens JC. The genetics of chemoreception in the labella and tarsi of Aedes aegypti. Insect Biochem Mol Biol. 2014;48:8–16.
Krieger J. A candidate olfactory receptor subtype highly conserved across different insect orders. J Comp Physiol A Neuroethol Sens Neural Behav Physiol. 2003;189:519–26.
Vosshall LB, Hansson BS. A unified nomenclature system for the insect olfactory coreceptor. Chem Senses. 2011;36:497–8.
Jones WD, Nguyen T-AT, Kloss B, Lee KJ, Vosshall LB. Functional conservation of an insect odorant receptor gene across 250 million years of evolution. Curr Biol. 2005;15:R119–21.
Münch D, Galizia CG. DoOR 2.0 – comprehensive mapping of Drosophila melanogaster odorant responses. Sci Rep. 2016;6:21841.
Choo Y-M, Pelletier J, Atungulu E, Leal WS. Identification and characterization of an antennae-specific aldehyde oxidase from the navel orangeworm. PLoS One. 2013;8, e67794.
Rybczynski R, Vogt RG, Lerner MR. Antennal-specific pheromone-degrading aldehyde oxidases from the moths Antheraea polyphemus and Bombyx mori. J Biol Chem. 1990;265:19712–5.
Pelletier J, Bozzolan F, Solvar M, François M-C, Jacquin-Joly E, Maïbèche-Coisne M. Identification of candidate aldehyde oxidases from the silkworm Bombyx mori potentially involved in antennal pheromone degradation. Gene. 2007;404:31–40.
Ishida Y, Leal WS. Chiral discrimination of the Japanese beetle sex pheromone and a behavioral antagonist by a pheromone-degrading enzyme. Proc Natl Acad Sci. 2008;105:9076–80.
Vogt RG. Molecular basis of pheromone detection in insects. Compr Insect Physiol Biochem Pharmacol Mol Biol. 2005;3:753–804.
Shi H, Pei L, Gu S, Zhu S, Wang Y, Zhang Y, et al. Glutathione S-transferase (GST) genes in the red flour beetle, Tribolium castaneum, and comparative analysis with five additional insects. Genomics. 2012;100:327–35.
Rogers ME, Jani MK, Vogt RG. An olfactory-specific glutathione-S-transferase in the sphinx moth Manduca sexta. J Exp Biol. 1999;202:1625–37.
Zhu F, Moural TW, Shah K, Palli SR. Integrated analysis of cytochrome P450 gene superfamily in the red flour beetle, Tribolium castaneum. BMC Genomics. 2013;14:174.
Forstner M, Gohl T, Gondesen I, Raming K, Breer H, Krieger J. Differential expression of SNMP-1 and SNMP-2 proteins in pheromone-sensitive hairs of moths. Chem Senses. 2008;33:291–9.
Leitch O, Papanicolaou A, Lennard C, Kirkbride KP, Anderson A. Chemosensory genes identified in the antennal transcriptome of the blowfly Calliphora stygia. BMC Genomics. 2015;16:255.
Gu X-C, Zhang Y-N, Kang K, Dong S-L, Zhang L-W. Antennal transcriptome analysis of odorant reception genes in the red turpentine beetle (RTB), Dendroctonus valens. PloS One. 2015;10, e0125159.
Wang J, Li D-Z, Min S-F, Mi F, Zhou S-S, Wang M-Q. Analysis of chemosensory gene families in the beetle Monochamus alternatus and its parasitoid Dastarcus helophoroides. Comp Biochem Physiol Part D Genomics Proteomics. 2014;11:1–8.
Andersson MN, Videvall E, Walden KK, Harris MO, Robertson HM, Löfstedt C. Sex- and tissue-specific profiles of chemosensory gene expression in a herbivorous gall-inducing fly (Diptera: Cecidomyiidae). BMC Genomics. 2014;15:501.
Syed Z, Leal WS. Maxillary palps are broad spectrum odorant detectors in culex quinquefasciatus. Chem Senses. 2007;32:727–38.
Rinker DC, Zhou X, Pitts RJ, Jones PL, Rokas A, Zwiebel LJ. RNAseq in the mosquito maxillary palp: a little antennal RNA goes a long way. bioRxiv. 2015:16998. http://dx.doi.org/10.1101/016998.
Dweck HK, Ebrahim SA, Khallaf MA, Koenig C, Farhan A, Stieber R, et al. Olfactory channels associated with the Drosophila maxillary palp mediate short- and long-range attraction. eLife. 2016;5:e14925.
Haverkamp A, Yon F, Keesey IW, Missbach C, Koenig C, Hansson BS, et al. Hawkmoths evaluate scenting flowers with the tip of their proboscis. eLife. 2016;5:e15039.
Ghaninia M, Hansson BS, Ignell R. The antennal lobe of the African malaria mosquito, Anopheles gambiae – innervation and three-dimensional reconstruction. Arthropod Struct Dev. 2007;36:23–39.
Anton S, van Loon JJA, Meijerink J, Smid HM, Takken W, Rospars J-P. Central projections of olfactory receptor neurons from single antennal and palpal sensilla in mosquitoes. Arthropod Struct Dev. 2003;32:319–27.
Ignell R, Dekker T, Ghaninia M, Hansson BS. Neuronal architecture of the mosquito deutocerebrum. J Comp Neurol. 2005;493:207–40.
Dahanukar A, Foster K, van der Goes van Naters WM, Carlson JR. A Gr receptor is required for response to the sugar trehalose in taste neurons of Drosophila. Nat Neurosci. 2001;4:1182–6.
Vermehren-Schmaedick A, Scudder C, Timmermans W, Morton DB. Drosophila gustatory preference behaviors require the atypical soluble guanylyl cyclases. J Comp Physiol A. 2011;197:717–27.
Bohbot J, Pitts RJ, Kwon H-W, Rützler M, Robertson HM, Zwiebel LJ. Molecular characterization of the Aedes aegypti odorant receptor gene family. Insect Mol Biol. 2007;16:525–37.
Kaupp UB. Olfactory signalling in vertebrates and insects: differences and commonalities. Nat Rev Neurosci. 2010;11:188–200.
Goldman AL, Van der Goes van Naters W, Lessing D, Warr CG, Carlson JR. Coexpression of two functional odor receptors in one neuron. Neuron. 2005;45:661–6.
Good NE. The flour beetles of the Genus Tribolium. United States Department of Agriculture; 1936
Ridley AW, Hereward JP, Daglish GJ, Raghu S, Collins PJ, Walter GH. The spatiotemporal dynamics of Tribolium castaneum (Herbst): adult flight and gene flow. Mol Ecol. 2011;20:1635–46.
Farhadian SF, Suárez-Fariñas M, Cho CE, Pellegrino M, Vosshall LB. Post-fasting olfactory, transcriptional, and feeding responses in Drosophila. Physiol Behav. 2012;105:544–53.
Rinker DC, Pitts RJ, Zhou X, Suh E, Rokas A, Zwiebel LJ. Blood meal-induced changes to antennal transcriptome profiles reveal shifts in odor sensitivities in Anopheles gambiae. Proc Natl Acad Sci. 2013;110:8260–5.
Hu J-H, Wang Z-Y, Sun F. Anatomical organization of antennal-lobe glomeruli in males and females of the scarab beetle Holotrichia diomphalia (Coleoptera: Melolonthidae). Arthropod Struct Dev. 2011;40:420–8.
Flook PK, Klee S, Rowell CHF. Combined molecular phylogenetic analysis of the Orthoptera (Arthropoda, Insecta) and implications for their higher systematics. Syst Biol. 1999;48:233–53.
Ignell R, Anton S, Hansson BS. The antennal lobe of Orthoptera – anatomy and evolution. Brain Behav Evol. 2001;57:1–17.
Farris SM, Roberts NS. Coevolution of generalist feeding ecologies and gyrencephalic mushroom bodies in insects. Proc Natl Acad Sci. 2005;102:17394–9.
Kondoh Y, Kaneshiro KY, Kimura K, Yamamoto D. Evolution of sexual dimorphism in the olfactory brain of Hawaiian Drosophila. Proc Biol Sci. 2003;270:1005–13.
Kleineidam CJ, Obermayer M, Halbich W, Rössler W. A macroglomerulus in the antennal lobe of leaf-cutting ant workers and its possible functional significance. Chem Senses. 2005;30:383–92.
Ågren L. Architecture of a lamellicorn flagellum (Phyllopertha horticola, Scarabaeidae, Coleoptera, Insecta). J Morphol. 1985;186:85–94.
Allsopp PG. Sexual dimorphism in the adult antennae of Antitrogus parvulus Britton and Lepidiota negatoria Blackburn (Coleoptera: Scarabaeidae: Melolonthinae). Aust J Entomol. 1990;29:261–6.
Ruther J, Reinecke A, Thiemann K, Tolasch T, Francke W, Hilker M. Mate finding in the forest cockchafer, Melolontha hippocastani, mediated by volatiles from plants and females. Physiol Entomol. 2000;25:172–9.
Jourdan H, Barbier R, Bernard J, Ferran A. Antennal sensilla and sexual dimorphism of the adult ladybird beetle Semiadalia undecimnotata Schn. (Coleoptera: Coccinellidae). Int J Insect Morphol Embryol. 1995;24:307–22.
Okada K, Mori M, Shimazaki K, Chuman T. Morphological studies on the antennal sensilla of the cigarette beetle, Lasioderma serricorne (F.) (Coleoptera: Anobiidae). Appl Entomol Zool. 1992;27:269–76.
Wu Z, Bin S, He H, Wang Z, Li M, Lin J. Differential expression analysis of chemoreception genes in the striped flea beetle Phyllotreta striolata using a transcriptomic approach. PLoS One. 2016;11, e0153067.
Fishilevich E, Vosshall LB. Genetic and functional subdivision of the Drosophila antennal lobe. Curr Biol. 2005;15:1548–53.
Grosse-Wilde E, Stieber R, Forstner M, Krieger J, Wicher D, Hansson BS. Sex-specific odorant receptors of the tobacco hornworm Manduca sexta. Front Cell Neurosci. 2010;4:22.
Kollmann M, Rupenthal AL, Neumann P, Huetteroth W, Schachtner J. Novel antennal lobe substructures revealed in the small hive beetle Aethina tumida. Cell Tissue Res. 2015;363:679–92.
Suzuki T. 4, 8-Dimethyldecanal: the aggregation pheromone of the flour beetles, Tribolium castaneum and T. confusum (Coleoptera: Tenebrionidae). Agric Biol Chem. 1980;44:2519–20.
Serrano JM, Castro L, Toro MA, López-Fanjul C. The genetic properties of homosexual copulation behavior in Tribolium castaneum: diallel analysis. Behav Genet. 1991;21:547–58.
Posnien N, Koniszewski NDB, Hein HJ, Bucher G. Candidate gene screen in the red flour beetle Tribolium reveals Six3 as ancient regulator of anterior median head and central complex development. PLoS Genet. 2011;7:e1002416.
Berghammer, Bucher G, Maderspacher F, Klingler M. A system to efficiently maintain embryonic lethal mutations in the flour beetle Tribolium castaneum. Dev Genes Evol. 1999;209:382–9.
Berghammer, Klingler M, Wimmer EA. A universal marker for transgenic insects. Nature. 1999;402:370–1.
Horn C, Offen N, Nystedt S, Häcker U, Wimmer EA. PiggyBac-Based insertional mutagenesis and enhancer detection as a tool for functional insect genomics. Genetics. 2003;163:647–61.
Horn C, Wimmer EA. A versatile vector set for animal transgenesis. Dev Genes Evol. 2000;210:630–7.
Siebert KS, Lorenzen MD, Brown SJ, Park Y, Beeman RW. Tubulin superfamily genes in Tribolium castaneum and the use of a Tubulin promoter to drive transgene expression. Insect Biochem Mol Biol. 2008;38:749–55.
Santos JG, Pollák E, Rexer K-H, Molnár L, Wegener C. Morphology and metamorphosis of the peptidergic Va neurons and the median nerve system of the fruit fly, Drosophila melanogaster. Cell Tissue Res. 2006;326:187–99.
Schmitt-Engel C, Schultheis D, Schwirz J, Ströhlein N, Troelenberg N, Majumdar U, et al. The iBeetle large-scale RNAi screen reveals gene functions for insect development and physiology. Nat Commun. 2015;6:7822.
Galizia CG, Joerges J, Küttner A, Faber T, Menzel R. A semi-in-vivo preparation for optical recording of the insect brain. J Neurosci Methods. 1997;76:61–9.
Pellegrino M, Nakagawa T, Vosshall LB. Single sensillum recordings in the insects Drosophila melanogaster and Anopheles gambiae. J Vis Exp. 2010;36, e1725.
Inkscape [Internet]. Available from: http://www.inkscape.org.
Huetteroth W, Schachtner J. Standard three-dimensional glomeruli of the Manduca sexta antennal lobe: a tool to study both developmental and adult neuronal plasticity. Cell Tissue Res. 2005;319:513–24.
Kurylas AE, Rohlfing T, Krofczik S, Jenett A, Homberg U. Standardized atlas of the brain of the desert locust, Schistocerca gregaria. Cell Tissue Res. 2008;333:125–45.
Brandt R, Rohlfing T, Rybak J, Krofczik S, Maye A, Westerhoff M, et al. Three-dimensional averageshape atlas of the honeybee brain and its applications. J. Comp. Neurol. 2005;492:1–19.
Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9:676–82.
Kent WJ. BLAT—the BLAST-like alignment tool. Genome Res. 2002;12:656–64.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Maglott D, Ostell J, Pruitt KD, Tatusova T. Entrez gene: gene-centered information at NCBI. Nucleic Acids Res. 2007;35:D26–31.
Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–5.
Zdobnov EM, Apweiler R. InterProScan–an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001;17:847–8.
Attrill H, Falls K, Goodman JL, Millburn GH, Antonazzo G, Rey AJ, et al. FlyBase: establishing a Gene Group resource for Drosophila melanogaster. Nucleic Acids Res. 2016;44:D786–92.
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, et al. CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011;39:D225–9.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Pavlidis P, Noble WS. Matrix2png: a utility for visualizing matrix data. Bioinformatics. 2003;19:295–6.
R Core Team. R: A Language and Environment for Statistical Computing. 2013. Available from: http://www.R-project.org.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Lawson D, Arensburger P, Atkinson P, Besansky NJ, Bruggner RV, Butler R, et al. VectorBase: a data resource for invertebrate vector genomics. Nucleic Acids Res. 2009;37:D583–7.
Katoh K, Kuma K, Toh H, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33:511–8.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.
Kapushesky M, Adamusiak T, Burdett T, Culhane A, Farne A, Filippov A, et al. Gene Expression Atlas update–a value-added database of microarray and sequencing-based functional genomics experiments. Nucleic Acids Res. 2011;40:D1077–81.
Letunic I, Bork P. Interactive Tree of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2007;23:127–8.
Wang X, Fang X, Yang P, Jiang X, Jiang F, Zhao D, et al. The locust genome provides insight into swarm formation and long-distance flight. Nat Commun. 2014;5:2957.
Grabe V, Strutz A, Baschwitz A, Hansson BS, Sachse S. Digital in vivo 3D atlas of the antennal lobe of Drosophila melanogaster. J Comp Neurol. 2015;523:530–44.
Hinke DW. Das relative postembryonale Wachstum der Hirnteile von Culex pipiens, Drosophila melanogaster und Drosophila-mutanten. Z Für Morphol Ökol Tiere. 1961;50:81–118.
Qiu YT, van Loon JJA, Takken W, Meijerink J, Smid HM. Olfactory coding in antennal neurons of the malaria mosquito, Anopheles gambiae. Chem Senses. 2006;31:845–63.
Pitts RJ, Zwiebel LJ. Antennal sensilla of two female anopheline sibling species with differing host ranges. Malar J. 2006;5:26.
McIver S. Structure of sensilla trichodea of female Aedes aegypti with comments on innervation of antennal sensilla. J Insect Physiol. 1978;24:383–90.
Kent LB, Walden KKO, Robertson HM. The Gr family of candidate gustatory and olfactory receptors in the yellow-fever mosquito Aedes aegypti. Chem Senses. 2008;33:79–93.
Koontz MA, Schneider D. Sexual dimorphism in neuronal projections from the antennae of silk moths (Bombyx mori, Antheraea polyphemus) and the gypsy moth (Lymantria dispar). Cell Tissue Res. 1987;249:39–50.
Lee JK, Strausfeld NJ. Structure, distribution and number of surface sensilla and their receptor cells on the olfactory appendage of the male moth Manduca sexta. J Neurocytol. 1990;19:519–38.
Homberg U, Christensen TA, Hildebrand JG. Structure and function of the deutocerebrum in insects. Annu Rev Entomol. 1989;34:477–501.
Kanost MR, Arrese EL, Cao X, Chen Y-R, Chellapilla S, Goldsmith MR, et al. Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta. Insect Biochem Mol Biol. 2016;76:118–47.
Rospars JP, Hildebrand JG. Sexually dimorphic and isomorphic glomeruli in the antennal lobes of the Sphinx moth Manduca sexta. Chem Senses. 2000;25:119–29.
Esslen J, Kaissling K-E. Zahl und Verteilung antennaler Sensillen bei der Honigbiene (Apis mellifera L.). Zoomorphologie. 1976;83:227–51.
Frasnelli E, Anfora G, Trona F, Tessarolo F, Vallortigara G. Morpho-functional asymmetry of the olfactory receptors of the honeybee (Apis mellifera). Behav Brain Res. 2010;209:221–5.
Galizia CG, Sachse S, Rappert A, Menzel R. The glomerular code for odor representation is species specific in the honeybee Apis mellifera. Nat Neurosci. 1999;2:473–8.
Mobbs PG. The brain of the honeybee Apis mellifera. I. The connections and spatial organization of the mushroom bodies. Philos Trans R Soc Lond B Biol Sci. 1982;298:309–54.
Strausfeld NJ. Organization of the honey bee mushroom body: representation of the calyx within the vertical and gamma lobes. J Comp Neurol. 2002;450:4–33.
Boeckh J, Ernst K-D. Contribution of single unit analysis in insects to an understanding of olfactory function. J Comp Physiol A. 1987;161:549–65.
Neder R. Allometrisches Wachstum von Hirnteilen bei drei verschieden grossen Schabenarten. Zool. Jb., Abt. allg. Zool. u. Physiol. 1959;77:411–464.
Greenwood M, Chapman RF. Differences in numbers of sensilla on the antennae of solitarious and gregarious Locusta migratoria L. (Orthoptera: Acrididae). Int J Insect Morphol Embryol. 1984;13:295–301.
Laurent G, Naraghi M. Odorant-induced oscillations in the mushroom bodies of the locust. J Neurosci. 1994;14:2993–3004.
Altner H. Insektensensillen: Bau und Funktionsprinzipien. Verhandlungen Dtsch Zool Ges. 1977;70:139–153.
Acknowledgements
We thank the transcriptome analysis laboratory of the University Medical Center Göttingen, Germany, especially Gabriela Salinas-Riester for sequencing and technical support; Jürgen Krieger for providing a cross-reactive previously unpublished Orco antibody; Gregor Bucher and Michalis Averof for sharing transgenic beetle lines and plasmids; Hugh M. Robertson and Kimberly K. O. Walden for providing unpublished sequence information; Mario Stanke und Lizzy Gerischer for the iBeetle Genome Browser; Richard Smith for support on MorphoGraphiX; Martina Kern, Marlene Binzer, and Peter Christ for technical assistance; Montserrat Torres Oliva and Jan Kropf for technical advice; Kei Ito for comments on the manuscript; as well as Uwe Homberg, and the members of the Deutsche Forschungsgemeinschaft Priority Program SPP 1392 “Integrative Analysis of Olfaction” for a fruitful discussion.
Funding
Deutsche Forschungsgemeinschaft SPP 1392: SCHA 678/13-1 (JS) and WI 1797/4-1 (EAW). The funders had no role in the study design, data collection, interpretation, or the decision to submit the work for publication.
Availability of data and materials
Previously Published Dataset: Tissue-specific transcriptomics, chromosomal localization, and phylogeny of chemosensory and odorant binding proteins from the red flour beetle Tribolium castaneum reveal subgroup specificities for olfaction or more general functions: Wimmer EA, Dippel S, Oberhofer G, 2015, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63162, GSE63162.
Authors’ contributions
SD and MK conceived and designed the study; acquired, analyzed, and interpreted the data; and drafted and revised the article. GO and RK analyzed, and interpreted the data. AM acquired, analyzed, and interpreted the data; and revised the article. CK and MK acquired the data and revised the article. KR acquired the data. SF acquired, analyzed, and interpreted the data. JS and EAW conceived and designed the study; analyzed and interpreted the data; and drafted and revised the article. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Author information
Authors and Affiliations
Corresponding authors
Additional files
Additional file 1: Figure S1.
Antibody staining against DsRed and Orco of the EF1-B-DsRed line. a Maximum projection of a confocal image stack of a halved antenna of the EF1-B-DsRed line, with an antibody staining against DsRed and Orco and in addition DAPI. Showing Orco immunoreactivity in the last three segments and particularly in the SBas. The DsRed reporter line labels in addition the scolopidia cells of Johnston's organ (JO) in the pedicellus (P). S, scapus. b Optical section of a mechano- and chemosensillum trichoideum (mSTri and cSTri) labeled with an Orco antibody (green) shows immunoreactivity only within the sensillum cavity of the cSTri. Autofluorescence of the cuticle at 560 nm is in blue. c–c'' Single optical section of a sensilla basiconica (SBas) in the EF1-B-DsRed (magenta, b') line labeled with an Orco antibody (green, b'') reveals signals of both channels particularly within the cavity and at the base of the sensillum. Both channels also show autofluorescence of the cuticle. d–d'' Optical section of two sensilla coeloconica (SCoe) in the EF1-B-DsRed (magenta, d') line labeled with an Orco antibody (green, d'') reveals no specific immunoreactivity within the sensilla cavities. Both channels also show autofluorescence of the cuticle. (TIF 8703 kb)
Additional file 2: Figure S2.
Orco localization in antennae and palp. a Fluorescent in situ hybridization against Orco in the club segments in a maximum projection of a halved antenna. b, c Specificity of the Orco antibody. IHC against Orco in antennae of b a San Bernadino beetle after Orco RNAi treatment (light cuticle, inset in the left lower corner) and c an untreated control included in the IHC (black strain identified by dark cuticle). The cross-reactive Orco antiserum results in no detectable staining in the antenna after Orco knockdown, whereas in the antenna of the untreated beetles, the odorant receptor neurons (ORNs) are clearly labeled by the antiserum (magenta). This indicates the specificity of the Orco antiserum against TcasOrco. Counterstaining with phalloidin (green) and DAPI (blue). The Orco antibody staining was labeled with a goat anti rabbit Cy3 secondary antibody. d–d'' Immunohistochemical characterization of the Orco-Gal4 line in the antenna and brain. Double immune-staining against Orco (magenta) and tGFP (green) in the antennae of the partial Orco-Gal4/UAS-tGFP line revealed that only half of the Orco-immunoreactive neurons expressing tGFP (arrow indicates colocalization and arrowhead as an example for no colocalization). e Immunohistochemical characterization of the Orco-Gal4 line in the brain. Antibody staining against tGFP in the Orco-Gal4/UAS-tGFP line (orange) labels only half of the AL glomeruli represented as a 3D reconstruction (light green, based on a phalloidin staining). AL glomeruli not labeled are not shown. f Immunohistochemical characterization of the partial Orco-Gal4 line in the palps. Double immuno-staining against Orco (green) and dsRed (magenta) in the palps of the Orco-Gal4/UAS-dsRed line reveals that all genetically labeled neurons are also Orco-immunoreactive. However, in contrast to the antennae, in which about half of the Orco-immunoreactive neurons are labeled in the partial Orco-Gal4 line, only a few of the the Orco-immunoreactive odorant receptor neurons express the reporter in the palps. The Cy2 (Orco) signal is quenched by dsRed. (TIF 6221 kb)
Additional file 3: Figure S3.
Comparison of sensilla type numbers on the antenna of Tribolium castaneum and chemosensory neurons entering the sensilla types. a–e Number of different sensilla types on the 11th segment of the antenna: a sensilla campaniformis (Scam: ♂ 5.4; SD 1.8; ♀ 5.4; SD 0.8), b spatulate bristles (SpaB: ♂ 10.2; SD 0.9; ♀ 10.3; SD 0.9), c mechanoreceptive sensilla trichoideum (mSTri: ♂ 36.9; SD 3; ♀ 37.6; SD 4.3), d chemoreceptive sensilla trichoideum (cSTri: ♂ 86.3; SD 9.3; ♀ 87.1; SD 6.9), e sensilla coeloconica (SCoe: ♂ 6.8; SD 1.4; ♀ 7.6; SD 1.1). f Amount of sensilla basiconica on the club segments (11th: ♂ 24.4; SD 1.5; ♀ 25.5; SD 1.3; tenth: ♂ 16.2; SD 1.5; ♀ 16.4; SD 1.4; ninth: ♂ 13.6; SD 0.9; ♀ 13.8; SD 1.1), regardless of the number of prongs. g Number of sensilla basiconica as in (f), but considering the prong number. h Number of different sensilla in the lateral corner of the tenth and ninth segments. i Number of chemosensory neurons (CSNs) entering the chemoreceptive sensilla: SBas 5.92 CSNs per prong (SD = 1.2; n = 73 prongs of total 48 SBas), cSTri 1.07 CSNs (SD = 0.25; n = 61), and SCoe 3.16 CSNs (SD = 1.10; n = 5). Error bars represent standard deviations; n = number of antennae. (TIF 1147 kb)
Additional file 6: Figure S4.
Ipsilateral antennal projection. Maximum intensity projection of a brain labeled with an antibody against synapsin (green) and a neurotracer resulting from an antennal backfill (magenta). The antennal backfill labels exclusively structures in the ipsilateral hemisphere, mainly the AL via the antennal nerve (*), and a tract (arrowhead) descending to the gnathal ganglion. The inset depicts a projection of only a few optical sections showing fibers interconnecting the AL and the protocerebrum with some arborizations in the accessory medulla of the optical lobe (arrow) suggesting an integration of circadian information. (TIF 3816 kb)
Additional file 8: Table S1.
Summary of the RNAseq data. In column (A) the gene name of the GRs according to [81], the ORs from [115], the SNMPs modified after [116, 117], the T. castaneum orthologous of D. melanogaster genes named after candidates obtained from [46], the GSTs named after [148], the CYPs named after [150] and the remaining ODE candidates de novo named according to their chromosomal localization (Additional file 15: Figure S11). Column B shows the sequences of the ORF based on published annotations or existing gene models, but modified if necessary. Confirmed gene models are highlighted in grey, modified ones are highlighted in yellow, only partially covered but expressed ones are highlighted in orange, genes with low and scattered coverage are highlighted in red, and not highlighted sequences were not manually checked. Columns C–I show the average RPKM values of antennae, mouthparts (piece of the head capsule anterior of the antennae), legs, head (without antennae but including mouthparts), and remaining body of sex-separated adult animals. Columns J–Q show the results of the statistical analysis conducted in R [215] with the DESeq package (version 1.12.0) [216] (from Bioconductor [217], based on five antenna samples and two replicates for the other adult tissues in comparison to body. Columns P–AB show the individual RPKM values of the biological replicates. (XLSX 677 kb)
Additional file 9: Figure S5.
IR gene tissue expression and chromosomal localization of IR and SNMP genes. a Venn diagram showing the number of IRs expressed (RPKM ≥ 0.5) in the different body parts: antennae, legs, mouthparts (as a piece of the head capsule anterior of the antennae), heads (the whole head capsule including mouthparts but excluding the antennae), and bodies (excluding head and legs). b Based on Georgia GA-2 strain genome assembly 3 [81], only chromosomal linkage groups containing an IR or SNMP are depicted. Gene clusters are indicated by a number referring to the chromosome and a letter conveys the relative position on the chromosome. The number of genes within this cluster is indicated in square brackets. (PDF 66 kb)
Additional file 10: Figure S6.
Phylogenetic tree of IRs. Outer rings represent the expression in body, mouthparts (T. castaneum: palps, mandible, labrum, and labium; D. melanogaster: palp and proboscis; An. gambiae: maxillary palp) and antenna as a percentage compared to the highest expressed gene according to the scale in the left upper corner. Note that the methods used to obtain the different expression data (RNAseq and microarray) are not directly comparable. This figure can, thus, only give an impression of the tissue-specific abundance of the transcripts. The scale bars within the trees represent one amino acid substitution per site. Antennal IRs are highlighted in yellow. (PDF 625 kb)
Additional file 11: Figure S7.
GR gene tissue expression and their chromosomal localization. a Venn diagram showing the number of GRs expressed (RPKM ≥ 0.5) in the different body parts: antennae, legs, mouthparts (as a piece of the head capsule anterior of the antennae), heads (the whole head capsule including mouthparts but excluding the antennae), and bodies (excluding head and legs). b Based on Georgia GA-2 strain genome assembly 3.0 [81], only chromosomal linkage groups containing an IR or SNMP are depicted. Gene clusters are indicated by a number referring to the chromosome and a letter conveys the relative position on the chromosome. The number of genes within this cluster is indicated in the square brackets. (PDF 176 kb)
Additional file 12: Figure S8.
Phylogenetic mid-point rooted tree of the GRs based on protein sequences. Outer rings represent the expression in body, mouthparts (T. castaneum: palps, mandible, labrum, and labium; D. melanogaster: palp and proboscis; An. gambiae: maxillary palp) and antenna as a percentage compared to the highest expressed gene according to the scale in the left upper corner. Note that the methods used to obtain the different expression data (RNAseq and microarray) are not directly comparable. This figure can, thus, only give an impression of the tissue-specific abundance of the transcripts. The scale bars within the trees represent 1 amino acid substitution per site. Potential sugar and fructose receptors are labeled and highlighted in yellow and in grey, and CO2 receptors are highlighted in orange. (PDF 1733 kb)
Additional file 13: Figure S9.
OR gene tissue expression and their chromosomal localization. a Venn diagram showing the number of ORs expressed (RPKM ≥ 0.5) in the different body parts: antennae, legs, mouthparts (as piece of the head capsule anterior of the antennae), heads (the whole head capsule including mouthparts but excluding the antennae), and bodies (excluding head and legs). b Venn diagram comparing our results (yellow, green) with data from Engsontia et al. [115] (blue, red). Number of expressed ORs, defined by RPKM ≥ 0.5 (yellow), by RT-PCR (blue), not expressed RPKM < 0.5 (green), or with no RT-PCR amplicon (red). ORs of the brown group were not previously tested by Engsontia et al. c Chromosomal localization of T. castaneum ORs. Based on the Georgia GA-2 strain genome assembly 3.0 [81], only chromosomal linkage groups containing an IR or SNMP are depicted. Gene clusters are indicated by a number referring to the chromosome and a letter conveys the relative position on the chromosome. The number of genes within this cluster is indicated in the square brackets. (PDF 277 kb)
Additional file 14: Figure S10.
Phylogenetic tree of the ORs based on protein sequences. Outer rings represent the expression in body, mouthparts (T. castaneum: palps, mandible, labrum, and labium; D. melanogaster: palp and proboscis; An. gambiae: maxillary palp) and antenna as a percentage compared to the highest expressed gene according to the scale in the left upper corner. Note that the methods used to obtain the different expression data (RNAseq and microarray) are not directly comparable. This figure can, thus, only give an impression of the tissue-specific abundance of the transcripts. The scale bars within the trees represent one amino acid substitution per site. (PDF 2302 kb)
Additional file 15: Figure S11.
Chromosomal localization of potential T. castaneum ODE genes. Based on Georgia GA-2 strain genome assembly 3.0 [81], aldehyde dehydrogenase (ALDH, in grey), aldehyde oxidase (ALOX, in orange), carboxylesterase (CES, in blue), epoxide hydrolase (EH, in green), glutathione S-transferase (GST, in purple), and cytochrome P450 (CYP, in magenta). Gene clusters are indicated by a number referring to the chromosome and a letter conveys the relative position on the chromosome. The number of genes within this cluster is indicated in the square brackets. (PDF 257 kb)
Additional file 16:
Sequences of primers and template plasmid used to generate pSLfa1180[2.5kbOrcoUp_GAL4delta]. (PDF 111 kb)
Additional file 17: Figure S12.
UAS responder lines in the absence of Gal4 driver. In the four rows, the maximum projections of head capsules from different transgenic strains (Orco-Gal4/UAS-tGFP, UAS-tGFP, Orco-Gal4/UAS-dsRed, and UAS-dsRed) are depicted. The upper row represents the overlay of both channels (GFP/YFP in green and dsRed in red). For UAS responders without the Gal4 driver, high-resolution images of the palps and the antennal club are provided. In the second and third rows, the separated channels are given as greyscale images. The UAS-tGFP and UAS-dsRed lines do not show leaky reporter expression in the absence of a Gal4 driver in the antennae and palps. The presence of the genetic constructs is indicated by the eye markers: pBac[3XP3-dsRed_UAS-Tchsp68bP-tGFP-SV40] and pBac[3XP3-eYFP_UAS-Tchsp68bP-DsRedex-SV40]. The marker signal is quenched in the crossed lines by the vermillion rescue marker of the Orco-Gal4 construct pBac[3XP3-gVerm_2.5kbOrcoUp_GAL4delta]. (TIF 6507 kb)
Additional file 18: Table S2.
Primary and secondary antibodies and dyes used with additional information such as source and specificity. n/a not available. (PDF 118 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Dippel, S., Kollmann, M., Oberhofer, G. et al. Morphological and Transcriptomic Analysis of a Beetle Chemosensory System Reveals a Gnathal Olfactory Center. BMC Biol 14, 90 (2016). https://doi.org/10.1186/s12915-016-0304-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12915-016-0304-z