A versatile recording device for the analysis of continuous daily external activity in colonies of highly eusocial bees.

As pollinators, bees are key to maintaining the biodiversity of angiosperm plants, and for agriculture they provide a billion-dollar ecosystem service. But they also compete for resources (primarily nectar and pollen), especially the highly social bees that live in perennial colonies. So, how do they organize their daily temporal activities? Here, we present a versatile, low-cost device for the continuous, automatic recording and data analysis of the locomotor activity in the colony-entrance tube of highly eusocial bees. Consisting of an in-house built block containing an infrared detector, the passage of bees in the colony entrance tunnel is registered and automatically recorded in an Arduino environment, together with concomitant recordings of temperature and relative humidity. With a focus on the highly diverse Neotropical stingless bees (Meliponini), we obtained 10-day consecutive recordings for two colonies each of the species Melipona quadrifasciata and Frieseomelitta varia, and also for the honey bee. The Lomb-Scargle periodogram analysis identified a predominant circadian rhythmicity for all three species, but also indications of ultradian rhythms. For M. quadrifasciata, which is comparable in size to the honey bee, we found evidence for a possibly anticipatory activity already before sunrise. As all three species also presented activity at night in the colony entrance tube, this also raises questions about sleep organization in social insects. The cost and versatility of the device and the open-source options for data analysis make this an attractive system for conducting studies on circadian rhythms in social bees under natural conditions, complementing studies on flower visits by these important pollinators.

Sex- and caste-specific transcriptomes of larval honey bee (Apis mellifera L.) gonads: DMRT A2 and Hsp83 are differentially expressed and regulated by juvenile hormone.

The gonads of honey bee, Apis mellifera, queens and drones are each composed of hundreds of serial units, the ovarioles and testioles, while the ovaries of the adult subfertile workers consist of only few ovarioles. We performed a comparative RNA-seq analysis on early fifth-instar (L5F1) larval gonads, which is a critical stage in gonad development of honey bee larvae. A total of 1834 genes were identified as differentially expressed (Padj < 0.01) among the three sex and caste phenotypes. The Gene Ontology analysis showed significant enrichment for metabolism, protein or ion binding, and oxidoreductase activity, and a KEGG analysis revealed metabolic pathways as enriched. In a principal component analysis for the total transcriptomes and hierarchical clustering of the DEGs, we found higher similarity between the queen and worker ovary transcriptomes compared to the drone testis, despite the onset of programmed cell death in the worker ovaries. Four DEGs were selected for RT-qPCR analyses, including their response to juvenile hormone (JH), which is a critical factor in the caste-specific development of the ovaries. Among these, DMRT A2 and Hsp83 were found upregulated by JH and, thus, emerged as potential molecular markers for sex- and caste-specific gonad development in honey bees.

DNA methylation is not a driver of gene expression reprogramming in young honey bee workers.

The presence of DNA methylation marks within genic intervals, also called gene body methylation, is an evolutionarily-conserved epigenetic hallmark of animal and plant methylomes. In social insects, gene body methylation is thought to contribute to behavioural plasticity, for example between foragers and nurse workers, by modulating gene expression. However, recent studies have suggested that the majority of DNA methylation is sequence-specific, and therefore cannot act as a flexible mediator between environmental cues and gene expression. To address this paradox, we examined whole-genome methylation patterns in the brains and ovaries of young honey bee workers that had been subjected to divergent social contexts: the presence or absence of the queen. Although these social contexts are known to bring about extreme changes in behavioral and reproductive traits through differential gene expression, we found no significant differences between the methylomes of workers from queenright and queenless colonies. In contrast, thousands of regions were differentially methylated between colonies, and these differences were not associated with differential gene expression in the subset of genes examined. Methylation patterns were highly similar between brain and ovary tissues and only differed in nine regions. These results strongly indicate that DNA methylation is not a driver of differential gene expression between tissues or behavioral morphs. Finally, despite the lack of difference in methylation patterns, queen presence affected the expression of all four DNA methyltransferase genes, suggesting that these enzymes have roles beyond DNA methylation. Therefore, the functional role of DNA methylation in social insect genomes remains an open question.

Juvenile hormone affects age polyethism, ovarian status and cuticular hydrocarbon profile in workers of the wasp Polybia occidentalis.

Division of labor is one of the most striking features in the evolution of eusociality. Juvenile hormone (JH) mediates reproductive status and aggression among nestmates in primitively eusocial Hymenoptera (species without morphologically distinct castes). In highly social species it has apparently lost its gonadotropic role and primarily regulates the division of labor in the worker caste. Polybia occidentalis, a Neotropical swarm-founding wasp, is an ideal model to understand how JH levels mirror social context and reproductive opportunities because of the absence of a clear morphological caste dimorphism. In this study, we tested the hypothesis that JH influences division of labor, ovary activation and cuticular hydrocarbon profiles of workers. Our observations confirmed that a JH analog (methoprene) and an inhibitor of JH biosynthesis (precocene) affected the cuticular chemical profile associated with age polyethism. Also, methoprene and precocene-I treatment of females influenced ovarian activation differently (individuals treated with methoprene expressed more activated ovaries while precocene treatment did not have significant effect). These results suggest that different hormonal levels induce a differential expression of cuticular chemicals associated with workers' age polyethism, which may be essential for keeping the social cohesion among workers throughout their lives in the colony. Furthermore, JH is likely to play a gonadotropic role in P. occidentalis. JH has apparently undergone certain modifications in social Hymenoptera, presenting multifaceted functions in different species.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste.

BACKGROUND: Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste. RESULTS: The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication. CONCLUSIONS: The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

Queen pheromone modulates the expression of epigenetic modifier genes in the brain of honeybee workers.

Pheromones are used by many insects to mediate social interactions. In the highly eusocial honeybee (Apis mellifera), queen mandibular pheromone (QMP) is involved in the regulation of the reproductive and other behaviour of workers. The molecular mechanisms by which QMP acts are largely unknown. Here, we investigate how genes responsible for epigenetic modifications to DNA, RNA and histones respond to the presence of QMP in the environment. We show that several of these genes are upregulated in the honeybee brain when workers are exposed to artificial QMP. We propose that pheromonal communication systems, such as those used by social insects, evolved to respond to environmental signals by making use of existing epigenomic machineries.

Differential expression of antioxidant system genes in honey bee (Apis mellifera L.) caste development mitigates ROS-mediated oxidative damage in queen larvae.

The expression of morphological differences between the castes of social bees is triggered by dietary regimes that differentially activate nutrient-sensing pathways and the endocrine system, resulting in differential gene expression during larval development. In the honey bee, Apis mellifera, mitochondrial activity in the larval fat body has been postulated as a link that integrates nutrient-sensing via hypoxia signaling. To understand regulatory mechanisms in this link, we measured reactive oxygen species (ROS) levels, oxidative damage to proteins, the cellular redox environment, and the expression of genes encoding antioxidant factors in the fat body of queen and worker larvae. Despite higher mean H2O2 levels in queens, there were no differences in ROS-mediated protein carboxylation levels between the two castes. This can be explained by their higher expression of antioxidant genes (MnSOD, CuZnSOD, catalase, and Gst1) and the lower ratio between reduced and oxidized glutathione (GSH/GSSG). In worker larvae, the GSG/GSSH ratio is elevated and antioxidant gene expression is delayed. Hence, the higher ROS production resulting from the higher respiratory metabolism in queen larvae is effectively counterbalanced by the up-regulation of antioxidant genes, avoiding oxidative damage. In contrast, the delay in antioxidant gene expression in worker larvae may explain their endogenous hypoxia response.

The Behavior and Reproductive Physiology of a Solitary Progressive Provisioning Vespid Wasp: Evidence for a Solitary-Cycle Origin of Reproductive Castes.

The emergence of queens and workers from solitary antecedents mark a major evolutionary transition in the history of life. The solitary progressive provisioning wasp Synagris cornuta, a member of the subfamily Eumeninae (basal to eusocial vespid wasps), alternates between behavioral states characterized as queenlike and worker-like. Akin to a queen in eusocial wasps, a S. cornuta female initiates construction of a cell into which she oviposits and then, similar to a worker, cares for the brood as it develops. The ovarian groundplan (OGP) hypothesis for caste origins predicts that these behavioral states are associated with cyclical changes in ovarian status, where females performing queenlike tasks have eggs and those performing worker-like tasks possess only small oocytes. Our findings show strong support for the OGP hypothesis: the ovaries of S. cornuta females undergo differential oogenesis depending on the behavioral phase: the largest oocyte in the ovaries of females building a cell progresses faster compared to that of females attending brood. Yet contrary to the OGP hypothesis, neither juvenile hormone nor ecdysteroids is associated with the reproductive cycle. Finally, the cuticular hydrocarbon profile showed no link with ovarian status, suggesting that fertility signals evolved subsequent to the emergence of group living.

Toward an Understanding of Divergent Compound Eye Development in Drones and Workers of the Honeybee (Apis mellifera L.): A Correlative Analysis of Morphology and Gene Expression.

Eye development in insects is best understood in Drosophila melanogaster, but little is known for other holometabolous insects. Combining a morphological with a gene expression analysis, we investigated eye development in the honeybee, putting emphasis on the sex-specific differences in eye size. Optic lobe development starts from an optic lobe anlage in the larval brain, which sequentially gives rise to the lobula, medulla, and lamina. The lamina differentiates in the last larval instar, when it receives optic nerve projections from the developing retina. The expression analysis focused on seven genes important for Drosophila eye development: eyes absent, sine oculis, embryonic lethal abnormal vision, minibrain, small optic lobes, epidermal growth factor receptor, and roughest. All except small optic lobes were more highly expressed in third-instar drone larvae, but then, in the fourth and fifth instar, their expression was sex-specifically modulated, showing shifts in temporal dynamics. The clearest differences were seen for small optic lobes, which is highly expressed in the developing eye of workers, and minibrain and roughest, which showed a strong expression peak coinciding with retina differentiation. A microarray analysis for optic lobe/retina complexes revealed the differential expression of several metabolism-related genes, as well as of two micro-RNAs. While we could not see major morphological differences in the developing eye structures before the pupal stage, the expression differences observed for the seven candidate genes and in the transcriptional microarray profiles indicate that molecular signatures underlying sex-specific optic lobe and retina development become established throughout the larval stages.

Mitochondrial capacity, oxidative damage and hypoxia gene expression are associated with age-related division of labor in honey bee (Apis mellifera L.) workers.

During adult life, honey bee workers undergo a succession of behavioral states. Nurse bees perform tasks inside the nest, and when they are about 2-3 weeks old they initiate foraging. This switch is associated with alterations in diet, and with the levels of juvenile hormone and vitellogenin circulating in hemolymph. It is not clear whether this behavioral maturation involves major changes at the cellular level, such as mitochondrial activity and the redox environment in the head, thorax and abdomen. Using high-resolution respirometry, biochemical assays and RT-qPCR, we evaluated the association of these parameters with this behavioral change. We found that tissues from the head and abdomen of nurses have a higher oxidative phosphorylation capacity than those of foragers, while for the thorax we found the opposite situation. As higher mitochondrial activity tends to generate more H2O2, and H2O2 is known to stabilize HIF-1α, this would be expected to stimulate hypoxia signaling. The positive correlation that we observed between mitochondrial activity and hif-1α gene expression in abdomen and head tissue of nurses would be in line with this hypothesis. Higher expression of antioxidant enzyme genes was observed in foragers, which could explain their low levels of protein carbonylation. No alterations were seen in nitric oxide (NO) levels, suggesting that NO signaling is unlikely to be involved in behavioral maturation. We conclude that the behavioral change seen in honey bee workers is reflected in differential mitochondrial activities and redox parameters, and we consider that this can provide insights into the underlying aging process.

Epigenetic modifications and their relation to caste and sex determination and adult division of labor in the stingless bee Melipona scutellaris.

Stingless bees of the genus Melipona, have long been considered an enigmatic case among social insects for their mode of caste determination, where in addition to larval food type and quantity, the genotype also has a saying, as proposed over 50 years ago by Warwick E. Kerr. Several attempts have since tried to test his Mendelian two-loci/two-alleles segregation hypothesis, but only recently a single gene crucial for sex determination in bees was evidenced to be sex-specifically spliced and also caste-specifically expressed in a Melipona species. Since alternative splicing is frequently associated with epigenetic marks, and the epigenetic status plays a major role in setting the caste phenotype in the honey bee, we investigated here epigenetic chromatin modification in the stingless bee Melipona scutellaris. We used an ELISA-based methodology to quantify global methylation status and western blot assays to reveal histone modifications. The results evidenced DNA methylation/demethylation events in larvae and pupae, and significant differences in histone methylation and phosphorylation between newly emerged adult queens and workers. The epigenetic dynamics seen in this stingless bee species represent a new facet in the caste determination process in Melipona bees and suggest a possible mechanism that is likely to link a genotype component to the larval diet and adult social behavior of these bees.

The scent of royalty: a p450 gene signals reproductive status in a social insect.

Cooperation requires communication; this applies to animals and humans alike. The main communication means differ between taxa and social insects (ants, termites, and some bees and wasps) lack the cognitive abilities of most social vertebrates. Central to the regulation of the reproductive harmony in insect societies is the production of a royalty scent which signals the fertility status of the reproducing queen to the nonreproducing workers. Here, we revealed a central genetic component underlying this hallmark of insect societies in the termite Cryptotermes secundus. Communication between queens and workers relied upon the expression of a gene, Neofem4, which belongs to the cytochrome P450 genes. We inhibited Neofem4 in queens by RNA interference. This resulted in the loss of the royalty scent in queens and the workers behaved as though the queen were absent. The queen's behavior was not generally affected by silencing Neofem4. This suggests that the lack of the royalty scent lead to workers not recognizing her anymore as queen. P450 genes are known to be involved in the production of chemical signals in cockroaches and their expression has been linked to a major fertility regulator, juvenile hormone. This makes P450 genes, both a suitable and available evolutionary substrate in the face of natural selection for production of a queen substance. Our data suggest that in an organism without elaborate cognitive abilities communication has been achieved by the exploitation of a central gene that links the fertility network with the chemical communication pathway. As termites and social Hymenoptera seem to share the same class of compounds in signaling fertility, this role of P450 genes might be more widespread across social insects.

Insights into the dynamics of hind leg development in honey bee (Apis mellifera L.) queen and worker larvae - A morphology/differential gene expression analysis.

Phenotypic plasticity is a hallmark of the caste systems of social insects, expressed in their life history and morphological traits. These are best studied in bees. In their co-evolution with angiosperm plants, the females of corbiculate bees have acquired a specialized structure on their hind legs for collecting pollen. In the highly eusocial bees (Apini and Meliponini), this structure is however only present in workers and absent in queens. By means of histological sections and cell proliferation analysis we followed the developmental dynamics of the hind legs of queens and workers in the fourth and fifth larval instars. In parallel, we generated subtractive cDNA libraries for hind leg discs of queen and worker larvae by means of a Representational Difference Analysis (RDA). From the total of 135 unique sequences we selected 19 for RT-qPCR analysis, where six of these were confirmed as differing significantly in their expression between the two castes in the larval spinning stage. The development of complex structures such as the bees' hind legs, requires diverse patterning mechanisms and signaling modules, as indicated by the set of differentially expressed genes related with cell adhesion and signaling pathways.

Identification of differentially expressed genes in female Drosophila antonietae and Drosophila meridionalis in response to host cactus odor.

BACKGROUND: Studies of insect-plant interactions have provided critical insights into the ecology and evolution of adaptive processes within and among species. Cactophilic Drosophila species have received much attention because larval development occurs in the necrotic tissues of cacti, and both larvae and adults feed on these tissues. Such Drosophila-cactus interactions include effects of the host plant on the physiology and behavior of the flies, especially so their nutritional status, mating condition and reproduction. The aim of this work was to compare the transcriptional responses of two species, Drosophila antonietae and Drosophila meridionalis, and identify genes potentially related to responses to odors released by their host cactus, Cereus hildmannianus. The two fly species are sympatric in most of their populations and use this same host cactus in nature. RESULTS: We obtained 47 unique sequences (USs) for D. antonietae in a suppression subtractive hybridization screen, 30 of these USs had matches with genes predicted for other Drosophila species. For D. meridionalis we obtained 81 USs, 46 of which were orthologous with genes from other Drosophila species. Functional information (Gene Ontology) revealed that these differentially expressed genes are related to metabolic processes, detoxification mechanisms, signaling, response to stimuli, and reproduction. The expression of 13 genes from D. meridionalis and 12 from D. antonietae were further analyzed by quantitative real time-PCR, showing that four genes were significantly overexpressed in D. antonietae and six in D. meridionalis. CONCLUSIONS: Our results revealed the differential expression of genes related to responses to odor stimuli by a cactus, in two associated fly species. Although the majority of activated genes were similar between the two species, we also observed that certain metabolic pathways were specifically activated, especially those related to signaling pathways and detoxification mechanisms. The activation of these genes may reflect different metabolic pathways used by these flies in their interaction with this host cactus. Our findings provide insight into how the use of C. hildmannianus may have arisen independently in the two fly species, through genetic differentiation in metabolic pathways to effectively explore this cactus as a host.

The role of juvenile hormone in dominance behavior, reproduction and cuticular pheromone signaling in the caste-flexible epiponine wasp, Synoeca surinama.

BACKGROUND: The popular view on insect sociality is that of a harmonious division of labor among two morphologically distinct and functionally non-overlapping castes. But this is a highly derived state and not a prerequisite for a functional society. Rather, caste-flexibility is a central feature in many eusocial wasps, where adult females have the potential to become queens or workers, depending on the social environment. In non-swarming paper wasps (e.g., Polistes), prospective queens fight one another to assert their dominance, with losers becoming workers if they remain on the nest. This aggression is fueled by juvenile hormone (JH) and ecdysteroids, major factors involved in caste differentiation in most eusocial insects. We tested whether these hormones have conserved aggression-promoting functions in Synoeca surinama, a caste-flexible swarm-founding wasp (Epiponini) where reproductive competition is high and aggressive displays are common. RESULTS: We observed the behavioral interactions of S. surinama females in field nests before and after we had removed the egg-laying queen(s). We measured the ovarian reproductive status, hemolymph JH and ecdysteroid titers, ovarian ecdysteroid content, and analyzed the cuticular hydrocarbon (CHC) composition of females engaged in competitive interactions in both queenright and queenless contexts. These data, in combination with hormone manipulation experiments, revealed that neither JH nor ecdysteroids are necessary for the expression of dominance behaviors in S. surinama. Instead, we show that JH likely functions as a gonadotropin and directly modifies the cuticular hydrocarbon blend of young workers to match that of a reproductive. Hemolymph ecdysteroids, in contrast, are not different between queens and workers despite great differences in ovarian ecdysteroid content. CONCLUSIONS: The endocrine profile of S. surinama shows surprising differences from those of other caste-flexible wasps, although a rise in JH titers in replacement queens is a common theme. Extensive remodeling of hormone functions is also evident in the highly eusocial bees, which has been attributed to the evolution of morphologically defined castes. Our results show that hormones which regulate caste-plasticity can lose these roles even while caste-plasticity is preserved.

Reproductive status, endocrine physiology and chemical signaling in the Neotropical, swarm-founding eusocial wasp Polybia micans.

In the evolution of caste-based societies in Hymenoptera, the classical insect hormones juvenile hormone (JH) and ecdysteroids were co-opted into new functions. Social wasps, which show all levels of sociality and lifestyles, are an ideal group in which to study such functional changes. Virtually all studies on the physiological mechanisms underlying reproductive division of labor and caste functions in wasps have been done on independent-founding paper wasps, and the majority of these studies have focused on species specially adapted for overwintering. The relatively little-studied tropical swarm-founding wasps of the Epiponini (Vespidae) are a diverse group of permanently social wasps, with some species maintaining caste flexibility well into the adult phase. We investigated the behavior, reproductive status, JH and ecdysteroid titers in hemolymph, ecdysteroid content of the ovary and cuticular hydrocarbon (CHC) profiles in the caste-monomorphic, epiponine wasp Polybia micans Ducke. We found that the JH titer was not elevated in competing queens from established multiple-queen nests, but increased in lone queens that lack direct competition. In queenless colonies, JH titer rose transiently in young potential reproductives upon challenge by nestmates, suggesting that JH may prime the ovaries for further development. Ovarian ecdysteroids were very low in workers but higher and correlated with the number of vitellogenic oocytes in the queens. Hemolymph ecdysteroid levels were low and variable in both workers and queens. Profiles of P. micans CHCs reflected caste, age and reproductive status, but were not tightly linked to either hormone. These findings show a significant divergence in hormone function in swarm-founding wasps compared with independently founding ones.

Insulin-like peptides (AmILP1 and AmILP2) differentially affect female caste development in the honey bee (Apis mellifera L.).

The food a honey bee female larva receives determines whether she develops into a large long-lived fertile queen or a short-lived sterile worker. Through well-established nutrient-sensing and growth-promoting functions in metazoans, the insulin/insulin-like growth factor 1 signaling (IIS) pathway has become a focal topic in investigations on how differences in food environment can be translated into internal signals responsible for queen-worker determination. However, low expression levels of two insulin receptors (AmInRs) in honey bee larvae and the failure of one AmInR to influence caste differentiation are in potential conflict with such a classical growth-promoting role of IIS in queen-worker development. In view of such an apparent contradiction, and the fact that binding partners and affinities of these two AmInRs have not been worked out, we performed a functional study on insulin-like peptide genes (AmILP1 and AmILP2) in honey bee larvae by using a double-stranded RNA (dsRNA)-mediated gene knockdown approach. We found that juvenile hormone (JH) levels were diminished by AmILP1 dsRNA treatment, while the AmILP2 knockdown caused a reduction in ovary size. Blood sugar titers were not significantly affected by the treatments. From these results we conclude that AmILP2 transcript levels may influence specific organ development, such as the ovary and body mass, while more general traits of caste differentiation, such as mandibles, may require additional regulators. In addition, JH production may be regulated by AmILP1 expressed locally in the brain, similar to the function of certain ILPs in Drosophila.

Positive selection in cytochrome P450 genes is associated with gonad phenotype and mating strategy in social bees.

The honey bee, Apis mellifera differs from all other social bees in its gonad phenotype and mating strategy. Honey bee queens and drones have tremendously enlarged gonads, and virgin queens mate with several males. In contrast, in all the other bees, the male and female gonads are small, and the females mate with only one or very few males, thus, suggesting an evolutionary and developmental link between gonad phenotype and mating strategy. RNA-seq comparisons of A. mellifera larval gonads revealed 870 genes as differentially expressed in queens versus workers and drones. Based on Gene Ontology enrichment we selected 45 genes for comparing the expression levels of their orthologs in the larval gonads of the bumble bee Bombus terrestris and the stingless bee, Melipona quadrifasciata, which revealed 24 genes as differentially represented. An evolutionary analysis of their orthologs in 13 solitary and social bee genomes revealed four genes with evidence of positive selection. Two of these encode cytochrome P450 proteins, and their gene trees indicated a lineage-specific evolution in the genus Apis, indicating that cytochrome P450 genes may be involved in the evolutionary association of polyandry and the exaggerated gonad phenotype in social bees.