Butterfly eyespots evolved via cooption of an ancestral gene-regulatory network that also patterns antennae, legs, and wings.

Butterfly eyespots are beautiful novel traits with an unknown developmental origin. Here we show that eyespots likely originated via cooption of parts of an ancestral appendage gene-regulatory network (GRN) to novel locations on the wing. Using comparative transcriptome analysis, we show that eyespots cluster most closely with antennae, relative to multiple other tissues. Furthermore, three genes essential for eyespot development, Distal-less (Dll), spalt (sal), and Antennapedia (Antp), share similar regulatory connections as those observed in the antennal GRN. CRISPR knockout of cis-regulatory elements (CREs) for Dll and sal led to the loss of eyespots, antennae, legs, and also wings, demonstrating that these CREs are highly pleiotropic. We conclude that eyespots likely reused an ancient GRN for their development, a network also previously implicated in the development of antennae, legs, and wings.

Identification and preliminary characterization of chemosensory-related proteins in the gall fly, Procecidochares utilis by transcriptomic analysis.

Chemoreception is critical for insect behaviors such as foraging, host searching and oviposition. The process of chemoreception is mediated by a series of proteins, including odorant-binding proteins (OBPs), gustatory receptors (GRs), odorant receptors (ORs), ionotropic receptors (IRs), chemosensory proteins (CSPs) and sensory neuron membrane proteins (SNMPs). The tephritid stem gall fly, Procecidochares utilis Stone, is a type of egg parasitic insect, which is an effective biological control agent for the invasive weed Ageratina adenophora in many countries. However, the study of molecular components related to the olfactory system of P. utilis has not been investigated. Here, we conducted the developmental transcriptome (egg, first-third instar larva, pupa, female and male adult) of P. utilis using next-generation sequencing technology and identified a total of 133 chemosensory genes, including 40 OBPs, 29 GRs, 24 ORs, 28 IRs, 6 CSPs, and 6 SNMPs. The sequences of these candidate chemosensory genes were confirmed by BLAST, and phylogenetic analysis was performed. Quantitative real-time PCR (qRT-PCR) confirmed that the expression levels of the candidate OBPs varied at the different developmental stages of P. utilis with most OBPs expressed mainly in the pupae, female and male adults but scarcely in eggs and larvae, which was consistent with the differentially expressed genes (DEGs) analysis using the fragments per kilobase per million fragments (FPKM) value. Our results provide a significant contribution towards the knowledge of the set of chemosensory proteins and help advance the use of P. utilis as biological control agents.

Interplay between axonal Wnt5-Vang and dendritic Wnt5-Drl/Ryk signaling controls glomerular patterning in the Drosophila antennal lobe.

Despite the importance of dendritic targeting in neural circuit assembly, the mechanisms by which it is controlled still remain incompletely understood. We previously showed that in the developing Drosophila antennal lobe, the Wnt5 protein forms a gradient that directs the ~45˚ rotation of a cluster of projection neuron (PN) dendrites, including the adjacent DA1 and VA1d dendrites. We report here that the Van Gogh (Vang) transmembrane planar cell polarity (PCP) protein is required for the rotation of the DA1/VA1d dendritic pair. Cell type-specific rescue and mosaic analyses showed that Vang functions in the olfactory receptor neurons (ORNs), suggesting a codependence of ORN axonal and PN dendritic targeting. Loss of Vang suppressed the repulsion of the VA1d dendrites by Wnt5, indicating that Wnt5 signals through Vang to direct the rotation of the DA1 and VA1d glomeruli. We observed that the Derailed (Drl)/Ryk atypical receptor tyrosine kinase is also required for the rotation of the DA1/VA1d dendritic pair. Antibody staining showed that Drl/Ryk is much more highly expressed by the DA1 dendrites than the adjacent VA1d dendrites. Mosaic and epistatic analyses showed that Drl/Ryk specifically functions in the DA1 dendrites in which it antagonizes the Wnt5-Vang repulsion and mediates the migration of the DA1 glomerulus towards Wnt5. Thus, the nascent DA1 and VA1d glomeruli appear to exhibit Drl/Ryk-dependent biphasic responses to Wnt5. Our work shows that the final patterning of the fly olfactory map is the result of an interplay between ORN axons and PN dendrites, wherein converging pre- and postsynaptic processes contribute key Wnt5 signaling components, allowing Wnt5 to orient the rotation of nascent synapses through a PCP mechanism.

Identification and phylogenetics of Spodoptera frugiperda chemosensory proteins based on antennal transcriptome data.

Understanding the interaction between the insect olfactory system and the environment is crucial for fully explaining the molecular mechanisms underlying insect behavior, and providing new strategies for integrated pest management. Although there is good evidence that olfactory proteins play a vital role in mediating insect behaviors, the olfactory mechanism of insects remains poorly understood. We identified a total of 71 chemosensory genes; 25 odorant-binding proteins (OBPs), 27 odorant receptors (ORs), 8 ionotropic receptors (IRs), 8 chemosensory proteins (CSPs) and 3 sensory neuron membrane proteins (SNMPs), in the antennae of male and female fall armyworms, Spodoptera frugiperda, an invasive global pest that causes significant economic damage worldwide. We used differential gene expression (DGE) and fragments per kilobase per million fragments (FPKM) values to compare the transcript levels of candidate chemosensory genes, and qRT-PCR to compare the expression levels of the OR gene, in male and female antennae. The expression of candidate OR genes in male and female antennae was consistent with the DGE data, and the expression of the SfruCL4419.Contig1-All and SfruUnigene1070-All genes was sex-biased. These results not only provide new information on the olfactory mechanism of S. frugiperda, and insects in general, but also suggest new gene targets for pest control.

Analysis of the antennal transcriptome and odorant-binding protein expression profiles of the parasitoid wasp Encarsia formosa.

The parasitoid of whiteflies Encarsia formosa has been widely applied to reduce whitefly-mediated damage on vegetables and ornamental plants grown in greenhouses. Although its chemosensory behavior has been described, the mechanism by which E. formosa recognizes chemical volatiles at the molecular level remains unknown. In this study, we obtained 66,632 unigenes from antennae transcriptomic architecture of E. formosa, of which 19,473 (29.2%) were functionally annotated. All that matters is that we manually identified 39 odorant-binding proteins (OBPs) from above dataset, and further investigated the tissue and stage-specific expression profiles of all identified OBP genes by real-time quantitative PCR. Among these OBP genes, 32 were enriched in antennae, and 2 in body. In addition, 4 OBPs were highly expressed in pupae, and 32 in 6-hour-age adults after eclosion. In addition to identifying OBP genes from E. formosa, this study provides a molecular basis for further functional studies of OBPs and the interactions of hosts and parasitic wasps.

The function of spineless in antenna and wing development of the brown planthopper, Nilaparvata lugens.

A wide array of sensilla are distributed on insect antennae, and they play a variety of important roles. Rice planthoppers, destructive pests on rice, have a unique antenna sensilla structure called the 'sensory plaque organ'. The spineless (ss) gene encodes a bHLH-PAS transcription factor and plays a key role in antenna development. In the current study, a 3029 bp full-length cDNA of the Nilaparvata lugens ss gene (Nlss) was cloned, and it encodes 654 amino acid residues. The highest level of Nlss expression was detected in the thorax of fourth-instar nymphs. Knockdown of Nlss in nymphs led to a decrease in the number and size of plaque organs. Moreover, the flagella of the treated insects were poorly developed, wilted, and even dropped off from the pedicel. Nlss-knockdown also resulted in twisted wings in both long-winged and short-winged brown planthoppers. Y-type olfactometer analyses indicated that antenna defects originating from Nlss depletion resulted in less sensitivity to host volatiles. This study represents the first report of the characteristics and functions of Nlss in N. lugens antenna and wing development and illuminates the function of the plaque organ of N. lugens in host volatile perception.

Embryonic origin and genetic basis of cave associated phenotypes in the isopod crustacean Asellus aquaticus.

Characteristics common to animals living in subterranean environments include the reduction or absence of eyes, lessened pigmentation and enhanced sensory systems. How these characteristics have evolved is poorly understood for the majority of cave dwelling species. In order to understand the evolution of these changes, this study uses an invertebrate model system, the freshwater isopod crustacean, Asellus aquaticus, to examine whether adult differences between cave and surface dwelling individuals first appear during embryonic development. We hypothesized that antennal elaboration, as well as eye reduction and pigment loss, would be apparent during embryonic development. We found that differences in pigmentation, eye formation, and number of segments of antenna II were all present by the end of embryonic development. In addition, we found that cave and surface hatchlings do not significantly differ in the relative size of antenna II and the duration of embryonic development. To investigate whether the regions responsible for eye and pigment differences could be genetically linked to differences in article number, we genotyped F2 hybrids for the four previously mapped genomic regions associated with eye and pigment differences and phenotyped these F2 hybrids for antenna II article number. We found that the region previously known to be responsible for both presence versus absence of pigment and eye size also was significantly associated with article number. Future experiments will address whether pleiotropy and/or genetic linkage play a role in the evolution of cave characteristics in Asellus aquaticus.

GABA immunostaining in the central complex of dicondylian insects.

The central complex is a group of midline-crossing neuropils in the insect brain involved in head direction coding, sky compass navigation, and spatial visual memory. To compare the neuroarchitecture and neurochemistry of the central complex in insects that differ in locomotion, ways of orientation, time of activity (diurnal, nocturnal), and evolutionary history, we studied the distribution of γ-aminobutyric acid (GABA) immunostaining in the central complex of 29 species, ranging from Zygentoma to Diptera. In all species, the lower division of the central body was densely innervated by GABA-immunoreactive tangential neurons. These neurons had additional arborizations in the bulb, a distinct region of synaptic complexes in the lateral complex, and somata in a cell cluster mediodorsally to the antennal lobe. Differences in the appearance of GABA immunostaining in the lower division of the central body corresponded to differences in neuropil architecture, such as transformation of the lower division into a toroid in certain Diptera and Heteroptera. In nearly all species two additional systems of tangential neuron of the upper division of the central body were GABA-immunoreactive. One of these systems diffusely invaded a superior layer, while the second system showed fan-like projections in an inferior layer. Sparse immunostaining in the protocerebral bridge was detected in cockroaches, a cricket, and two hemipteran species. The data show that three systems of GABA-immunoreactive tangential neurons of the central body are highly conserved and suggest that the layered organization of the upper division of the central body is, likewise, largely maintained from basal to advanced species.

Diverse populations of local interneurons integrate into the Drosophila adult olfactory circuit.

Drosophila olfactory local interneurons (LNs) in the antennal lobe are highly diverse and variable. How and when distinct types of LNs emerge, differentiate, and integrate into the olfactory circuit is unknown. Through systematic developmental analyses, we found that LNs are recruited to the adult olfactory circuit in three groups. Group 1 LNs are residual larval LNs. Group 2 are adult-specific LNs that emerge before cognate sensory and projection neurons establish synaptic specificity, and Group 3 LNs emerge after synaptic specificity is established. Group 1 larval LNs are selectively reintegrated into the adult circuit through pruning and re-extension of processes to distinct regions of the antennal lobe, while others die during metamorphosis. Precise temporal control of this pruning and cell death shapes the global organization of the adult antennal lobe. Our findings provide a road map to understand how LNs develop and contribute to constructing the olfactory circuit.

Cell fate changes induced by a Distal-less enhancer-trap transgene in the Drosophila antennal imaginal disc.

The imaginal discs of the genetically tractable model organism Drosophila melanogaster have been used to study cell-fate specification and plasticity, including homeotic changes and regeneration-induced transdetermination. The identity of the reprogramming mechanisms that induce plasticity has been of great interest in the field. Here we identify a change from antennal fate to eye fate induced by a Distal-less-GAL4 (DllGAL4) P-element insertion that is a mutant allele of Dll and expresses GAL4 in the antennal imaginal disc. While this fate change is not induced by tissue damage, it appears to be a hybrid of transdetermination and homeosis as the GAL4 expression causes upregulation of Wingless, and the Dll mutation is required for the fate change. Neither GAL4 expression nor a Dll mutation on its own is able to induce antenna-to-eye fate changes. This plasticity appears to be unique to the DllGAL4 line, possibly due to cellular stress induced by the high GAL4 expression combined with the severity of the Dll mutation. Thus, we propose that even in the absence of tissue damage, other forms of cellular stress caused by high GAL4 expression can induce determined cell fates to change, and selector gene mutations can sensitize the tissue to these transformations.

Transcriptome analysis in different developmental stages of Batocera horsfieldi (Coleoptera: Cerambycidae) and comparison of candidate olfactory genes.

The white-striped longhorn beetle Batocera horsfieldi (Coleoptera: Cerambycidae) is a polyphagous wood-boring pest that causes substantial damage to the lumber industry. Moreover olfactory proteins are crucial components to function in related processes, but the B. horsfieldi genome is not readily available for olfactory proteins analysis. In the present study, developmental transcriptomes of larvae from the first instar to the prepupal stage, pupae, and adults (females and males) from emergence to mating were built by RNA sequencing to establish a genetic background that may help understand olfactory genes. Approximately 199 million clean reads were obtained and assembled into 171,664 transcripts, which were classified into 23,380, 26,511, 22,393, 30,270, and 87, 732 unigenes for larvae, pupae, females, males, and combined datasets, respectively. The unigenes were annotated against NCBI's non-redundant nucleotide and protein sequences, Swiss-Prot, Gene Ontology (GO), Pfam, Clusters of Eukaryotic Orthologous Groups (KOG), and KEGG Orthology (KO) databases. A total of 43,197 unigenes were annotated into 55 sub-categories under the three main GO categories; 25,237 unigenes were classified into 26 functional KOG categories, and 25,814 unigenes were classified into five functional KEGG Pathway categories. RSEM software identified 2,983, 3,097, 870, 2,437, 5,161, and 2,882 genes that were differentially expressed between larvae and males, larvae and pupae, larvae and females, males and females, males and pupae, and females and pupae, respectively. Among them, genes encoding seven candidate odorant binding proteins (OBPs) and three chemosensory proteins (CSPs) were identified. RT-PCR and RT-qPCR analyses showed that BhorOBP3, BhorCSP2, and BhorOBPC1/C3/C4 were highly expressed in the antenna of males, indicating these genes may may play key roles in foraging and host-orientation in B. horsfieldi. Our results provide valuable molecular information about the olfactory system in B. horsfieldi and will help guide future functional studies on olfactory genes.

Anticipatory flexibility: larval population density in moths determines male investment in antennae, wings and testes.

Developmental plasticity provides individuals with a distinct advantage when the reproductive environment changes dramatically. Variation in population density, in particular, can have profound effects on male reproductive success. Females may be easier to locate in dense populations, but there may be a greater risk of sperm competition. Thus, males should invest in traits that enhance fertilization success over traits that enhance mate location. Conversely, males in less dense populations should invest more in structures that will facilitate mate location. In Lepidoptera, this may result in the development of larger antennae to increase the likelihood of detecting female sex pheromones, and larger wings to fly more efficiently. We explored the effects of larval density on adult morphology in the gum-leaf skeletonizer moth, Uraba lugens, by manipulating both the number of larvae and the size of the rearing container. This experimental arrangement allowed us to reveal the cues used by larvae to assess whether absolute number or density influences adult responses. Male investment in testes size depended on the number of individuals, while male investment in wings and antennae depended upon larval density. By contrast, the size of female antennae and wings were influenced by an interaction of larval number and container size. This study demonstrates that male larvae are sensitive to cues that may reveal adult population density, and adjust investment in traits associated with fertilization success and mate detection accordingly.

Notch-dependent epithelial fold determines boundary formation between developmental fields in the Drosophila antenna.

Compartment boundary formation plays an important role in development by separating adjacent developmental fields. Drosophila imaginal discs have proven valuable for studying the mechanisms of boundary formation. We studied the boundary separating the proximal A1 segment and the distal segments, defined respectively by Lim1 and Dll expression in the eye-antenna disc. Sharp segregation of the Lim1 and Dll expression domains precedes activation of Notch at the Dll/Lim1 interface. By repressing bantam miRNA and elevating the actin regulator Enable, Notch signaling then induces actomyosin-dependent apical constriction and epithelial fold. Disruption of Notch signaling or the actomyosin network reduces apical constriction and epithelial fold, so that Dll and Lim1 cells become intermingled. Our results demonstrate a new mechanism of boundary formation by actomyosin-dependent tissue folding, which provides a physical barrier to prevent mixing of cells from adjacent developmental fields.

Comparative analysis of behavioral and transcriptional variation underlying CO2 sensory neuron function and development in Drosophila.

Carbon dioxide is an important environmental cue for many insects, regulating many behaviors including some that have direct human impacts. To further improve our understanding of how this system varies among closely related insect species, we examined both the behavioral response to CO2 as well as the transcriptional profile of key developmental regulators of CO2 sensory neurons in the olfactory system across the Drosophila genus. We found that CO2 generally evokes repulsive behavior across most of the Drosophilids we examined, but this behavior has been lost or reduced in several lineages. Comparisons of transcriptional profiles from the developing and adult antennae for subset these species suggest that behavioral differences in some species may be due to differences in the expression of the CO2 co-receptor Gr63a. Furthermore, these differences in Gr63a expression are correlated with changes in the expression of a few genes known to be involved in the development of the CO2 circuit, namely dac, an important regulator of sensilla fate for sensilla that house CO2 ORNs, and mip120, a member of the MMB/dREAM epigenetic regulatory complex that regulates CO2 receptor expression. In contrast, most of the other known structural, molecular, and developmental components of the peripheral Drosophila CO2 olfactory system seem to be well-conserved across all examined lineages. These findings suggest that certain components of CO2 sensory ORN development may be more evolutionarily labile, and may contribute to differences in CO2-evoked behavioral responses across species.

De novo analysis of the oriental armyworm Mythimna separata antennal transcriptome and expression patterns of odorant-binding proteins.

To better understand the olfactory mechanisms in the oriental armyworm Mythimna separate, one of the most serious pests of cereals, an antennal transcriptome was constructed in this study. A total of 130 olfactory related transcripts were identified. These transcripts were predicted to encode 32 odorant-binding proteins (OBPs), 16 chemosensory proteins (CSPs), 71 olfactory receptors (ORs), 8 ionotropic receptors (IRs), 1 gustatory receptor (GR) and 2 sensory neuron membrane proteins (SNMPs). Q-PCR analysis of the temporal expression profiles of seven OBPs in different tissues indicated that, except for MsepOBP19 which was highly expressed in the wings of 0-day-old adult and MsepOBP20 which was low expressed in all tissues, other tested MsepOBPs were significantly more highly expressed in the antenna than in the head (antenna excluded), thorax, abdomen, legs and wings. The expression levels of MsepOBPs were diverse in different life stages (differed on eclosion days). MsepOBP5 exhibited female-biased expression in 0- and 5-day-old adult, while no gender bias in 1- and 3-day-old adult was detected and similar expression profiles were found for MsepOBP7, 20, 24 and 26. In addition, we found that although the expression of MsepOBP22 was female biased in 0- and 5-day-old adult, in the 3-day-old adult it was male-biased. Our findings established a foundation for future studies of the functions of olfactory proteins in M. separata.

Semaphorin-1a prevents Drosophila olfactory projection neuron dendrites from mis-targeting into select antennal lobe regions.

Elucidating how appropriate neurite patterns are generated in neurons of the olfactory system is crucial for comprehending the construction of the olfactory map. In the Drosophila olfactory system, projection neurons (PNs), primarily derived from four neural stem cells (called neuroblasts), populate their cell bodies surrounding to and distribute their dendrites in distinct but overlapping patterns within the primary olfactory center of the brain, the antennal lobe (AL). However, it remains unclear whether the same molecular mechanisms are employed to generate the appropriate dendritic patterns in discrete AL glomeruli among PNs produced from different neuroblasts. Here, by examining a previously explored transmembrane protein Semaphorin-1a (Sema-1a) which was proposed to globally control initial PN dendritic targeting along the dorsolateral-to-ventromedial axis of the AL, we discover a new role for Sema-1a in preventing dendrites of both uni-glomerular and poly-glomerular PNs from aberrant invasion into select AL regions and, intriguingly, this Sema-1a-deficient dendritic mis-targeting phenotype seems to associate with the origins of PNs from which they are derived. Further, ectopic expression of Sema-1a resulted in PN dendritic mis-projection from a select AL region into adjacent glomeruli, strengthening the idea that Sema-1a plays an essential role in preventing abnormal dendritic accumulation in select AL regions. Taken together, these results demonstrate that Sema-1a repulsion keeps dendrites of different types of PNs away from each other, enabling the same types of PN dendrites to be sorted into destined AL glomeruli and permitting for functional assembly of olfactory circuitry.

Candidate chemosensory genes identified in the endoparasitoid Meteorus pulchricornis (Hymenoptera: Braconidae) by antennal transcriptome analysis.

Meteorus pulchricornis is an endoparasitoid wasp which attacks the larvae of various lepidopteran pests. We present the first antennal transcriptome dataset for M. pulchricornis. A total of 48,845,072 clean reads were obtained and 34,967 unigenes were assembled. Of these, 15,458 unigenes showed a significant similarity (E-value <10-5) to known proteins in the NCBI non-redundant protein database. Gene ontology (GO) and cluster of orthologous groups (COG) analyses were used to classify the functions of M. pulchricornis antennae genes. We identified 16 putative odorant-binding protein (OBP) genes, eight chemosensory protein (CSP) genes, 99 olfactory receptor (OR) genes, 19 ionotropic receptor (IR) genes and one sensory neuron membrane protein (SNMP) gene. BLASTx best hit results and phylogenetic analysis both indicated that these chemosensory genes were most closely related to those found in other hymenopteran species. Real-time quantitative PCR assays showed that 14 MpulOBP genes were antennae-specific. Of these, MpulOBP6, MpulOBP9, MpulOBP10, MpulOBP12, MpulOBP15 and MpulOBP16 were found to have greater expression in the antennae than in other body parts, while MpulOBP2 and MpulOBP3 were expressed predominately in the legs and abdomens, respectively. These results might provide a foundation for future studies of olfactory genes and chemoreception in M. pulchricornis.

Structure and postembryonic development of the intersegmental nodules in the non-muscular joints of the antennae in Rhodnius prolixus.

The antennae of Insecta consist of two basal segments and the distal annulated flagellum lacking intrinsic muscles. Non-muscular joints are important to preserve the flexibility and structure of the long heteropteran antennae which bear an intersegmental nodule on each non-muscular joint. Little is known about their properties or function. Here we characterize the structure and postembryonic development of the non-muscular joints of Rhodnius prolixus antennae. Using Scanning Electron Microscopy, we tracked the changes in shape and size of both intersegmental nodules during the course of the hemimetabolous insect life cycle. Using Atomic Force Microscopy, we established a qualitative correlation between the topography of the surface and the rigidity of the joint between pedicel and flagellum. Our results confirmed the presence of two sub-articulations on each non-muscular joint. Also, the two intersegmental nodules have different origins: the one between the two flagellar segments (intraflagelloid) is a sclerite already present from the early nymph, while the nodule between pedicel and flagellum (prebasiflagellite) originates by gradual separation of the proximal end of the basiflagellum during postembryonic development. Various changes occur in the non-muscular joints and segments of the antenna during the life cycle of R. prolixus.

Ontogeny of pioneer neurons in the antennal nervous system of the grasshopper Schistocerca gregaria.

The nervous system of the antenna of the grasshopper Schistocerca gregaria consists of two nerve tracts in which sensory cells project their axons to the brain. Each tract is pioneered early in embryogenesis by a pair of identified cells located apically in the antennal lumen. The pioneers are thought to originate in the epithelium of the antenna and then delaminate into the lumen where they commence axogenesis. However, unambiguous molecular identification of these cells in the epithelium, of an identifiable precursor, and of their mode of generation has been lacking. In this study, we have used immunolabeling against neuron-specific horseradish peroxidase and against Lachesin, a marker for differentiating epithelial cells, in combination with the nuclear stain DAPI, to identify the pioneers within the epithelium of the early embryonic antenna. We then track their delamination into the lumen as differentiated neurons. The pioneers are not labeled by the mesodermal/mesectodermal marker Mes3, consistent with an epithelial (ectodermal) origin. Intracellular dye injection, as well as labeling against the mitosis marker phospho-histone 3, identifies precursor cells in the epithelium, each associated with a column of cells. Culturing with the S-phase label 5-ethynyl-2'-deoxyuridine (EdU) shows that both a precursor and its column have incorporated the label, confirming a lineage relationship. Each set of pioneers can be shown to belong to a separate lineage of such epithelial cells, and the precursors remain and are proliferative after generating the pioneers. Analyses of mitotic spindle orientation then enable us to propose a model in which a precursor generates its pioneers asymmetrically via self-renewal.

Identification of Complete Repertoire of Apis florea Odorant Receptors Reveals Complex Orthologous Relationships with Apis mellifera.

We developed a computational pipeline for homology based identification of the complete repertoire of olfactory receptor (OR) genes in the Asian honey bee species, Apis florea Apis florea is phylogenetically the most basal honey bee species and also the most distant sister species to the Western honey bee Apis mellifera, for which all OR genes had been identified before. Using our pipeline, we identified 180 OR genes in A. florea, which is very similar to the number of ORs identified in A. mellifera (177 ORs). Many characteristics of the ORs including gene structure, synteny of tandemly repeated ORs and basic phylogenetic clustering are highly conserved. The composite phylogenetic tree of A. florea and A. mellifera ORs could be divided into 21 clades which are in harmony with the existing Hymenopteran tree. However, we found a few nonorthologous OR relationships between both species as well as independent pseudogenization of ORs suggesting separate evolutionary changes. Particularly, a subgroup of the OR gene clade XI, which had been hypothesized to code cuticular hydrocarbon receptors showed a high number of species-specific ORs RNAseq analysis detected a total number of 145 OR transcripts in male and 162 in female antennae. Most of the OR genes were highly expressed on the female antennae. However, we detected five distinct male-biased OR genes, out of which three genes (AfOr11, AfOr18, AfOr170P) were shown to be male-biased in A. mellifera, too, thus corroborating a behavioral function in sex-pheromone communication.