Effects of Insemination Quantity on Honey Bee Queen Physiology

Abstract

Mating has profound effects on the physiology and behavior of female insects, and in honey bee (Apis mellifera) queens, these changes are permanent. Queens mate with multiple males during a brief period in their early adult lives, and shortly thereafter they initiate egg-laying. Furthermore, the pheromone profiles of mated queens differ from those of virgins, and these pheromones regulate many different aspects of worker behavior and colony organization. While it is clear that mating causes dramatic changes in queens, it is unclear if mating number has more subtle effects on queen physiology or queen-worker interactions; indeed, the effect of multiple matings on female insect physiology has not been broadly addressed. Because it is not possible to control the natural mating behavior of queens, we used instrumental insemination and compared queens inseminated with semen from either a single drone (single-drone inseminated, or SDI) or 10 drones (multi-drone inseminated, or MDI). We used observation hives to monitor attraction of workers to SDI or MDI queens in colonies, and cage studies to monitor the attraction of workers to virgin, SDI, and MDI queen mandibular gland extracts (the main source of queen pheromone). The chemical profiles of the mandibular glands of virgin, SDI, and MDI queens were characterized using GC-MS. Finally, we measured brain expression levels in SDI and MDI queens of a gene associated with phototaxis in worker honey bees (Amfor). Here, we demonstrate for the first time that insemination quantity significantly affects mandibular gland chemical profiles, queen-worker interactions, and brain gene expression. Further research will be necessary to elucidate the mechanistic bases for these effects: insemination volume, sperm and seminal protein quantity, and genetic diversity of the sperm may all be important factors contributing to this profound change in honey bee queen physiology, queen behavior, and social interactions in the colony.

About the Authors:

Freddie-Jeanne Richard

Affiliations Department of Entomology, North Carolina State University, Raleigh, North Carolina, United States of America, Department of Genetics, North Carolina State University, Raleigh, North Carolina, United States of America, W.M. Keck Center for Behavioural Biology, North Carolina State University, Raleigh, North Carolina, United States of America

David R. Tarpy

Affiliations Department of Entomology, North Carolina State University, Raleigh, North Carolina, United States of America, W.M. Keck Center for Behavioural Biology, North Carolina State University, Raleigh, North Carolina, United States of America

Christina M. Grozinger

*To whom correspondence should be addressed. E-mail: [email protected]

Affiliations Department of Entomology, North Carolina State University, Raleigh, North Carolina, United States of America, Department of Genetics, North Carolina State University, Raleigh, North Carolina, United States of America, W.M. Keck Center for Behavioural Biology, North Carolina State University, Raleigh, North Carolina, United States of America

Introduction

While many previous studies have considered the effects of mating on the physiology and behavior of female insects (i.e., reference 1), the effects of mating number and semen quantity have not been broadly considered. Honey bees are an excellent system in which to study this process: colonies typically consist of a single reproductive queen and thousands of sterile worker bees, and the queen mates with several males (12 males on average; reviewed in reference 2) during a brief period in her early adult life. However, studies of the effects of mating number have focused almost exclusively on the consequences of worker genetic diversity on overall colony health [3]–[5] and social interactions [6]. Relatively little is known about how queen mating number may alter queen physiology and thereby alter intracolony social interactions. Understanding the physiological effects of multiple mating could offer insights into the mechanisms that govern mating behavior and how the process is regulated in honey bees and female insects in general, as well as what consequences this may have for colony organization in bees and other social insects.

Virgin honey bee queens initiate mating very early in their lives, when they are approximately 1–2 weeks old, by taking multiple mating flights and mating with numerous males (drones) on each flight. On average, queens are successfully inseminated by 12 males (based on molecular genotyping of workers [2]), but mating number is highly variable among queens (range from 1 to 28). The factors that determine mating frequency are not fully understood [7], [8]. Mating has profound and permanent effects on queen behavior, physiology, and resultant queen-worker interactions. Once they begin to oviposit, mated queens will never mate again and will remain in the colony for the rest of their lives (unless they depart during colony swarming). Mating stimulates vitellogenesis and oocyte-maturation in the ovaries [9], which prompts the initiation of egg-laying of up to 1500 eggs/day [10]. There are also profound changes in a queen's brain after mating, where levels of dopamine significantly decrease [11] and the ratio of the neuropil/cell body volume in the mushroom bodies significantly increases [12].

Finally, mated and virgin queens differ dramatically in their pheromone profiles, and these pheromones are important for regulating colony organization and worker behaviour [13], [14]. Queens produce pheromones from a variety of glands, and the complete “queen pheromone” has not been fully characterized [15]–[17]. However, five components in particular, termed “queen mandibular pheromone” or QMP, are produced by the mandibular gland and elicit many similar worker behavioral and physiological responses as a live queen. QMP consists of 9-oxo-(E)-2-decenoic acid (ODA), (R)- and (S)-9-hydroxy-(E)-2-decenoic acid (9-HDA), methyl p-hydroxybenzoate (HOB), and 4-hydroxy-3-methyoxyphenylethanol (HVA) [18]. QMP has both primer (long-term) and releaser (short-term) effects on worker behaviour and physiology. It elicits a “retinue response”, in which workers are attracted from a distance (several cm) and then antennate and groom the queen [18]–[21]. QMP also inhibits worker ovary activation [22], [23], inhibits queen cell production by workers [24], [25], stimulates pollen foraging and brood rearing in small, newly founded colonies [26], increases nectar foraging [27], [28], and delays the age-at-onset of foraging and reduces juvenile hormone secretion [27]. Some components of QMP are involved in drone attraction and mating [29], [30]. However, while the five component QMP blend is a potent modulator of worker behaviour and physiology, it is still not as potent as a live queen, suggesting that additional pheromonal components may exist. Indeed, four additional compounds were identified that synergize with QMP to increase the retinue response, but even those nine compounds do not elicit as robust a response as a live queen [17].

Although it is clear that mandibular gland pheromone profiles differ dramatically between young virgins and naturally mated, laying queens, the short-term effects of mating on pheromone profiles are largely unknown. Plettner et al. [14] compared the quantities of the QMP components between 6-day-old virgins and 1-year-old mated laying queens, and found that mated queens had significantly higher levels of 9-ODA, 9-HDA, HOB and HVA than did virgins. However, a similar study by Slessor et al [13], which compared 6- and 12-day-old virgins with mated queens that had been laying for 1 day, 5 weeks, or 2 years, found different results. In that study, levels of 9-ODA were not significantly different between any of these groups, 9-HDA levels were significantly higher in mated versus virgin queens, HOB levels were significantly higher in the 2-year-old mated queens than in both the virgins and 1-day laying queens and intermediate in the 5-week laying queens, and HVA levels were significantly higher in the 2-year mated queens compared to all other groups. Thus, pheromone profiles are strongly affected by mating, age (or perhaps time since mating), and potentially by egg-laying. Finally, Apsegaite and Skirkevicius [31] found slight differences in the quantity of 9-ODA between different strains of honey bees, suggesting that genotypic differences might also alter pheromone profiles.

Given the highly variable mating number and the profound effects mating has on queens, we hypothesize that mating number could modulate queen physiology, behavior, and pheromone production. However, “mating number” comprises a host of factors, including the number and duration of mating flights, the physical act of insemination, the volume of ejaculate, the quantity and viability of sperm and seminal proteins in the ejaculate, and the genetic diversity of sperm and seminal proteins. Furthermore, given that queens mate during flights several meters above the ground, it is not possible to control natural mating behavior. However, the number of males that inseminate a queen can be precisely controlled by using instrumental insemination [32]. Previous studies have demonstrated that queens inseminated with lower quantities of semen (less than 8 µL) are more likely to continue to take mating flights [33], but the effects of insemination quantity on other aspects of queen physiology have not been considered. To begin to address these questions, we focused on single-drone inseminated (SDI, total semen from one drone) and multi-drone inseminated (MDI, total semen from 10 drones) queens. We monitored attraction of the workers to SDI or MDI queens in observation hive colonies, as well as the attraction of workers to virgin, SDI, and MDI queen mandibular gland extracts in cages. The chemical profiles of the mandibular glands of virgin, SDI, and MDI queens were characterized using GC-MS. Finally, since fully inseminated queens cease to be phototactic and no longer take mating flights, we measured brain expression levels in SDI and MDI queens of a gene that is associated with phototaxis in worker honey bees (Amfor, the foraging gene; [34].

Results

Behavioral assays

Observation hives.

One of the most measurable effects of queen pheromone is the induced retinue response, in which workers are attracted to the queen from a short distance, and lick and antennate her. Following insemination, the retinue response to SDI and MDI queens was monitored in observation hives, twice a day for 5 days. MDI queens attracted significantly more worker bees in their retinue than the SDI queens (F = 6.73, p = 0.02; Figure 1). There was no significant effect of day (F = 0.46, p = 0.76) or day*treatment interaction (F = 0.68, p = 0.61).

[Figure omitted. See PDF.]

Figure 1. Workers are more attracted to MDI than SDI queens in observation hives.

Following insemination, the retinue response to SDI and MDI queens was monitored in observation hives, twice a day for 5 days. MDI queens attracted significantly more worker bees in their retinue than the SDI queens did (F = 6.73, p = 0.02). (SDI: 7 queens; MDI: 5 queens).

https://doi.org/10.1371/journal.pone.0000980.g001

Cage assays.

Worker retinue responses to the mandibular gland extracts of virgin, SDI, and MDI queens were tested in cages containing 4-day-old bees. The retinue size was equivalent whether workers were exposed to virgin, SDI, or MDI queen mandibular gland extracts (H(2,18) = 0.615 p = 0.735; data not shown). However, worker bees exposed to two different mandibular gland extracts at the same time preferred gland extracts from SDI and MDI queens to virgins (respectively: t = 4.2, df = 36, p<0.001 and t = 3.1, df = 38, p<0.01, Figure 2), and preferred MDI extracts to SDI extracts (t = 2.7, df = 40, p<0.01). Data for the preference assay represent two trials which were not significantly different and were pooled for the subsequent analysis.

[Figure omitted. See PDF.]

Figure 2. Workers are more attracted to the mandibular gland extracts of MDI than SDI queens.

Caged workers were exposed to mandibular gland extract of two different queens simultaneously, and the retinue response to each extract was monitored. Workers were significantly more attracted to the extracts of inseminated versus virgin queens, and significantly more attracted to the extracts of MDI versus SDI queens (0.05 queen equivalents: Virgin vs SDI queens (n = 19); Virgin vs MDI (n = 20). MDI vs SDI (n = 21); all comparison with t-test: t<2.7, p<0.01). ** p<0.01 and *** p<0.001

https://doi.org/10.1371/journal.pone.0000980.g002

Chemical analysis of queen mandibular gland profile

The chemical composition of the mandibular gland extract of virgin, SDI, and MDI queens were analysed using gas chromatography. A total of 27 compounds were apparent in the GC analysis and used in the subsequent comparisons. There were significant differences in the overall chemical profiles between the three groups of queens, as revealed by a discriminant analysis (F(33, 36) = 9.33, p<10−4; Figure 3). Two discriminant variables explained 100% of the variation. The chemical profiles of the virgin queen mandibular glands were most different from the two groups of mated queens (85% of the variation), and the other 15% of the variation could be attributed to insemination quantity.

[Figure omitted. See PDF.]

Figure 3. Insemination quantity significantly alters the chemical profile of mandibular glands.

Chemical composition of the mandibular gland extracts of Virgin, SDI and MDI queens were analysed using gas-chromatography. Discriminant analysis of mandibular extract of virgin, SDI, and MDI queens was based on the relative proportion of the chemical compounds (F(33, 36) = 9.33. p<10−4; All group distances <0.005). Ellipses have been drawn to emphasize the categories, but have no specific statistical meaning

https://doi.org/10.1371/journal.pone.0000980.g003

Next, when the actual quantities of the compounds are analyzed (Table 1), mated queens had significantly lower levels of 18 compounds relative to virgins, and higher levels of one compound. MDI queens had significantly lower levels of seven compounds than SDI queens, and significantly higher levels of two compounds. Of the five QMP components, the quantities of 9-ODA, 9-HDA, and HVA were all significantly lower in mated queens compared to virgins. Levels of 9-ODA and 9-HDA were also significantly lower in the MDI queens than in the SDI queens.

[Figure omitted. See PDF.]

Table 1. Absolute quantity of the mandibular gland compounds (mean±STD) in µg of virgin, SDI and MDI queens.

https://doi.org/10.1371/journal.pone.0000980.t001

Finally, we compared the relative proportions (compared to the total gland quantity) of the 27 individual compounds between virgin and mated queens, and between SDI and MDI queens (Table 2). The relative proportion of 13 compounds was significantly different between virgins and mated queens, and 10 of these compounds were significantly higher in mated queen mandibular glands. There were fewer compounds with significantly different relative proportions between SDI and MDI queens. Two compounds (8-hydroxyoctanoic acid and unidentified 5) were present in significantly higher proportions in SDI queens than in MDI queens (Table 2), while six compounds (4-hydroxybenzoic acid, unidentified 2, (E)-dec-2-enedioic acid, palmitic acid, alkane 2 and stearic acid) were significantly higher in MDI than in SDI queens mandibular glands (Table 2). Of the QMP components, relative levels of 9-ODA were significantly lower in mated queens compared to virgins, while HOB levels were significantly higher. None of the QMP components differed between SDI and MDI queens.

[Figure omitted. See PDF.]

Table 2. Chemical identity and relative proportion (mean±STD) of the queen mandibular gland of honey bees Apis mellifera belonging to different queen group: virgin, single (SDI) and multi (MDI) drone inseminated queens.

https://doi.org/10.1371/journal.pone.0000980.t002

Brain gene expression levels of Amfor

We used quantitative real-time PCR to measure Amfor expression levels in the queens' brains to determine if insemination treatment altered brain gene expression and therefore presumably altered neuronal properties in the brain. Amfor expression was significantly higher in SDI queens than in MDI queens (U = 10, p = 0.02; Figure 4).

[Figure omitted. See PDF.]

Figure 4. Brain gene expression levels are affected by insemination quantity.

Expression levels of Amfor were monitored using quantitative real-time PCR using individual brains (SDI: N = 9; MDI: N = 9). SDI queens have significantly higher levels of Amfor than MDI queens (U = 10, p = 0.022). Data represent mean values, −/+ standard error (converted to the same arbitrary scale as the mean), relative to the SDI brains.

https://doi.org/10.1371/journal.pone.0000980.g004

Discussion

In this study, we instrumentally inseminated queen honey bees to determine if insemination quantity alters their physiology and social interactions under controlled environmental conditions. Our results clearly demonstrate that insemination quantity alters queen physiology, queen pheromone profiles, and queen-worker interactions. MDI queens elicited a stronger retinue than SDI queens in natural colony conditions, and their mandibular gland extracts were more attractive in preference assays with caged worker bees. Moreover, GC-MS analysis of the mandibular gland chemical profiles revealed significant differences between SDI and MDI queens. Finally, brain-expression levels of a behaviorally relevant gene were significantly different between the SDI and MDI queens. All of these results suggest that insemination quantity can have profound effects on queen physiology and behaviour.

Comparisons of the mandibular gland profiles between MDI and SDI queens may lead to the identification of new components of queen pheromone. Previously, studies focusing on compounds that elicited a retinue response identified the five-component QMP blend [18]. Another study identified four additional compounds that did not elicit a retinue on their own, but synergized with QMP to produce a stronger retinue response [17]. However, the retinue response was still not as strong as that elicited by the live queens, suggesting that additional active compounds may still exist. In this study, we found that MDI queens elicited a larger retinue than SDI queens in observation hive colonies. Since the mandibular glands are the main source of QMP, we tested the mandibular gland extracts of virgin, SDI, and MDI queens on worker retinue responses in cages. When tested individually, the glands were equally attractive, though all glands are qualitatively comparable and contained the main QMP compounds [18]. However, when the glands were compared in a choice test, there were clear differences. Indeed, the glands of the inseminated queens were more attractive than gland of virgins, while the gland extracts of MDI queens were more attractive than gland of SDI queens. This suggests that MDI queen mandibular gland extracts have additional compounds that synergize with 9-ODA and 9-HDA to make these extracts more attractive in a preference assay and that the relative proportions of the entire gland chemical profile affect the retinue response. Indeed, when the extracts were analyzed with GC-MS, the chemical profiles of the mandibular glands of the inseminated queens were significantly different from those of the virgins, and those of the MDI queens were significantly different from those of the SDI queens.

Our results suggest that these changes in pheromone profiles occur within days after insemination, and that additional changes occur over a long period of time. Previous studies have compared the queen mandibular gland compounds between 6-day-old virgin and 1-year-old mated laying queens, and mated queens queens had higher levels of 9-ODA, 9-HDA, HOB, HVA, and 10-HDAA, and lower levels of 10-HDA [14]. However, when virgins and mated queens of more similar ages were compared, the differences were less robust [13]. In particular, 12-day-old virgins and mated queens that were laying for only 1 day did not differ significantly between levels of 9-ODA, 9-HDA, HOB, or HVA. In our study, virgins and inseminated queens were age-matched, and all queens were collected 5 days after insemination and were not laying eggs. However, we still observed significant differences in the overall chemical profiles and in several individual compounds, suggesting that insemination quantity causes changes in pheromone production within days. Thus, our results suggest that changes in mandibular gland profiles occur immediately or shortly after insemination, but additional modifications occur over time that may be associated with age, egg laying, or both. Furthermore, the observation hives assays with live queens and the cage studies with queen mandibular gland extract demonstrate that workers can detect these differences in pheromone profiles and respond differently to virgins versus newly inseminated queens, and singly versus multiply inseminated queens.

Expression levels of a gene associated with phototaxis are also significantly altered by insemination quantity, suggesting that insemination quantity exerts effects of neuronal properties in the brain and possibly behavior. We focused our analysis on brain expression levels of the foraging gene (Amfor; GenBank Accession AF469010). Previous studies found that expression of Amfor is significantly higher in forager bees than in-hive bees [35]. Moreover, Amfor encodes a cGMP-dependent protein kinase (PKG) and treatment with a cGMP analog stimulates foraging behavior [35] and increases phototaxis [36]. Our results demonstrate that there are significant differences in gene expression levels in brains of SDI and MDI queens, and these differences may be associated with changes in flight and phototaxis, though behavioural studies are clearly necessary to determine if this is the case. However, previous work demonstrates that queens inseminated with lower quantities of semen (<8 µL) have an increased likelihood to attempt mating flights [37].

Clearly, insemination quantity has profound effects on queens, but what are the mechanisms by which insemination quantity triggers these changes in queen honey bees? As noted previously, a number of factors might play a role. The number of stored sperm may modulate queen post-mating physiology and behavior [38]. It is highly unlikely that queens can directly assess the number of sperm in their spermathecae since there is no nerve enervation into the spermathecae and the membranes of their spermathecae are rigid [39], [40], so there is likely to be some other mechanism—one that is possibly correlated with stored sperm number—that serves as the basis for these physiological and behavioral changes. One possible explanation is that insemination volume may be the cue that triggers these changes. As queens are inseminated with numerous drones, the large quantity of semen is temporarily stored in their lateral and median oviducts for several hours before a proportion (5–10%) of each male's sperm migrates into the spermathecae and the excess semen is excreted. Indeed, preliminary studies (Richard, Tarpy, Grozinger, unpublished data) suggest insemination volume may be a key contributor to the observed post-insemination changes. One possible mechanism to induce these changes is that queens may have stretch receptors in their oviducts that provide a negative stimulus for additional mating behavior. Similar mechanisms exist in other insect systems [41]–[43], and abdomen distension has been tested in honey bees [44], but it remains unclear if abdominal enlargement is a physiological mechanism for reproductive senescence in queen bees. Alternatively, queens may use a substance in their mates' ejaculate, such as a hormonal precursor or seminal protein that causes physiological changes and decreases their subsequent reproductive behavior. Male seminal proteins have shown to have significant inhibitory effects on mating in female Drosophila [45]–[52] and possibly bumble bees [53], [54], and several substances in drone honey bee ejaculate have been identified, including sugars, metal ions, and proteins [55]–[57]. Further studies will be necessary to determine which factor(s) is involved in causing the differences in queen pheromone profile and brain gene expression that we have observed between SDI and MDI queens.

Previous studies have demonstrated that there may be colony-level adaptive benefits for genetically diverse workers, in terms of increased resistance to disease [5], [58], homozygosity at the sex determination locus [59], or better regulation of colony division of labor [60], [61]. Here, we demonstrate that insemination quantity, and thereby potentially mating number, could also affect different aspects of queen physiology or quality that could affect colony fitness. Additional research will be necessary to see if the physiological and pheromonal differences of differentially inseminated or mated queens are used by the workers as honest signals of their insemination quality or fertility [62]–[64].

Materials and Methods

Bee rearing

Colonies headed by single-drone inseminated queens (Apis mellifera carnica, Glenn Apiaries, CA) were maintained at the NCSU Lake Wheeler Honey Bee Research Facility (Raleigh, NC). These colonies were used as genetic sources for producing supersister queens for the experiments. Due to the haplodiploid genome of the hymenoptera, the progeny of SDI-queens have a genetic relatedness of G = 0.75. Additionally, colonies headed by naturally mated Apis mellifera ligustica or Buckfast-SMR queens (B. Weaver, TX) were used to provide workers for the observation-hive and cage experiments. Note that while worker bees were derived from different source colonies for convenience, each individual experiment used workers from a single colony when testing differential responses to SDI or MDI queens and/or pheromone extract. Thus, though there might be differences between colonies in worker responses to queen pheromones (i.e., reference 65), this should not be a factor in our results.

Queen rearing

Supersister queens were produced by grafting young larvae (<24 h) from a single source colony and reared as queens in a queenless colony [32]. Once all queen cells were capped, they were transferred to a dark incubator at 33°C and ∼40% relative humidity. One to two days prior to emergence, each capped queen cell was placed into an individual Plexiglas cage (10×10×7 cm). Frames of brood were removed from a hive and incubated at 33°C. Day-old workers were collected from the brood frames and 100 adult worker bees were placed in the cages with the capped queen cell. Bees were provided with water, food (45% honey/45% pollen/10% water), and a solution of 50% sucrose, and all cages were maintained in an incubator (at 33°C and ∼40% RH). The food was changed every two days. Five days after emergence, queens destined for insemination were treated with CO2 for 4.0 minutes [66], while virgins were handled but left untreated, since we wanted a baseline nonreproductive control group and CO2 treatment accelerates the transition to egg-laying [66]. Seven days after emergence, queens were again treated with CO2 for 4.0 minutes, during which time they were inseminated with semen from either one drone (single-drone inseminated treatment group, or SDI) or 10 brother drones (multi-drone inseminated treatment group, or MDI) by following standard insemination protocols which readily produce laying queens [32], though in our experiments queens were reared in cages without honeycomb, so there was no opportunities for them to lay eggs. The average drone produces approximately 1 µl of semen, thus SDI queens were inseminated with approximately 1 µl and MDI queens were inseminated with approximately 10 µl. Queens were then returned to their respective cages and collected 5 days after insemination onto dry ice and stored at 80°C. A total of 12 virgin, 15 SDI and 11 MDI queens were collected. Queens were reared in three separate cohorts (Cohorts 1∶ 5 SDI and 4 MDI; Cohorts 2∶ 5 SDI and 2 MDI and Cohorts 3∶ 12 virgin, 5 SDI and 5 MDI), but queens from different cohorts were combined for the subsequent chemical analysis and behavioural assays.

Prior to dissection, queen heads were removed and partially lyophilized to facilitate dissection [67], after which their mandibular glands and brains were dissected out on dry ice and stored for future processing. In order to verify insemination success, we counted the number of sperm in each spermatheca (SDI, N = 15 and MDI, N = 8). As would be expected, spermathecae of MDI queens (65*104±25*104 sperm) contained significantly more sperm than those of SDI queens (15*104±8*104 sperm; U = 20.0, p = 0.0098). Note that these quantities of sperm are lower than typical for inseminated queens, but the lack of comb, the presence of a small number of workers, and the small size of the cages may have minimize sperm storage [38].

Behavioral assays

Observation hives.

Supersister queens were produced as described above. Capped queen cells were removed from the queen-rearing colony and individually placed into 1.5 frame observation hives (34″×5″×21.5″ LxWxH). The entrance to each hive was blocked with queen excluder material so the queens could not exit the hives to take mating flights. Observation hives were stocked with ∼1000 one-day-old bees and a half-frame of honey/pollen. The observation hives were stored in a 33°C incubator for 3 days (to ensure that workers were old enough to properly thermoregulate) prior to being placed in a room with access to the outside via a window. As the queens emerged, individual plastic numbered tags (Thorne Ltd, UK) were glued to their thoraxes. Approximately five days after emergence, queens were captured and treated with carbon dioxide for 4.0 minutes. Two days later, queens were captured again and treated with carbon dioxide for 4.0 minutes during their insemination procedure. As outlined above, queens were inseminated with semen from either a single drone (single-drone inseminated, SDI) or 10 brother drones (multi-drone inseminated, MDI). Queens were then released back into their respective hives. One hour later, the number of worker bees in the queen's retinue (surrounding the queen, licking and antennating her) was counted every minute over a period of 5.0 minutes [17], [18]. One hour later, the number of workers contacting the queen was counted again. The retinue size was calculated by summing across all of the observations, as in previous studies [17], [68]. The retinue was monitored every day, twice a day, for five days. In total, 2 SDI and 2 MDI queens were assayed in 2005 and 5 SDI and 3 MDI in 2006.

Differences in the retinue size for SDI and MDI queens in observation hives were analyzed using a mixed-model ANOVA (PROC mixed in SAS, Cary, NC), with treatment and day viewed as fixed factors, while trial, trial*treatment and replicate (trial*treatment) were random effects. Thus ‘day’ was treated as a split plot factor.

Cage assays.

We next compared worker responses to mandibular gland extracts of virgin, SDI, and MDI queens. The glands were extracted from each supersister queen produced as outlined above.

Frames of brood were removed from a hive and incubated at 33°C. Day-old workers were collected from the brood frames and 35 bees were placed in Plexiglas cages (10×10×7 cm). Bees were provisioned as described above. The cages were kept in a 33°C incubator with ∼40% relative humidity, and manipulations and observations were performed under red light to negate any potential behavioural effects. Cages were maintained for five days.

Experiment 1: Retinue response

The retinue response was monitored as described previously [17], [21], [68], [69]. Extracts of mandibular glands of individual queens were produced as described below (see chemical analysis section). 0.05 queen equivalents (Qeq) of queen mandibular gland extract from individual queens (Virgin: n = 3; SDI: n = 8; MDI: n = 7) was placed onto a glass coverslip, the solvent was allowed to evaporate, and the coverslip was then placed inside the cage at the same time every day for five days. The number of bees antennating and licking the coverslip was counted in each cage (“retinue”). The retinue was recorded 5, 10, and 15 minutes after pheromone introduction each day of the 5-day time course. A single trial was performed.

Experiment 2: Preference assay

This experiment was repeated twice (Trial 1: Fall 2005 and Trial 2: Spring 2006). Bees were exposed to 0.1 Qeq of synthetic QMP (Pherotech, Canada). Every day, 10 µl of QMP was placed on a microscope slide and allowed to evaporate before being placed in the cage. This amount of QMP mimics a live queen in assays of worker behaviour and physiology [70]. Since worker maturation is altered in the absence of queen pheromone [67], [68], [71], we reared these workers with QMP to mimic natural colony conditions. On the fifth day of the experiment, workers were presented with two slides containing equal quantities of extract (0.05 Qeq) from either a virgin and SDI queen (Trial 1: n = 9 and Trial 2: n = 10), a virgin and MDI queen (Trial 1: n = 10 and Trial 2: n = 10), or a SDI and MDI queen (Trial 1: n = 11 and Trial 2: n = 10). The number of workers contacting each slide was counted every 5.0 minutes for 15 minutes after slide presentation. Note that a total of 6 virgin, 6 SDI, and 6 MDI queens were used in this analysis. Queens were derived from 2 cohorts. Trial 1 and trial 2 were not significantly different, so the data were pooled for the subsequent analysis.

Data are presented as mean±SEM. The effect of the mated number was evaluated with a non-parametric Kruskall-Wallis ANOVA on ranks for global comparison. The worker preferences between two slides were evaluated with a parametric t-test.

Chemical analysis

The mandibular glands were dissected and immersed in 50 µl of diethyl ether containing 0.4 µg/µl of undec-10-enoic acid (as an internal standard) for minimum of 24 hours. A 5 µl portion (approx. 0.1 bee equivalents) of an extract was placed in a small glass insert and the solvent gently evaporated. The residue was sylylated overnight at room temperature in the insert with 10 µl of neat bistrimethylsilyltrifluoroacetamide (BSTFA) [17]. The derivatized sample was diluted with hexane (100 µl) and a 2 µl portion was analyzed using gas chromatography on a HP 5890 equipped with capillary column (30 m×0.25 mm ID. 0.5 µm film thickness) DB-5 (5% diphenyl-95% dimethylsiloxane) column (J&W scientific, Folsom, CA) in splitless mode. Helium was used as the carrier gas at a head pressure of 18 psi (flow rate = 1.3 ml/min). The GC temperature was held at 100°C for 1 min and then increased at 5°C/min to 200°C (5 min), followed by an increase of 10°C/min to 250°C (15 min). Injector and FID temperatures were both set at 250°C. We extracted the mandibular gland of 12 virgin queens, 15 SDI queens and 10 MDI queens. All queens used for chemical analysis were raised from the same grafting source to reduce any genetic variation in their pheromone profiles. However, the queens were reared in three different cohorts.

Compound identification was achieved by splitless capillary gas-chromatography-mass spectrometry using a Hewlett-Packard 6890 GC and a model 5973A msd with an electron impact ion source and a HP-5ms capillary column (30 m×0.25 mm ID×0.25 µm film thickness).

To examine differences in profiles related to number of inseminations based on the relative proportion of the chemical compounds, a stepwise discriminant analysis was employed, using all the chemical compounds (Statistica 6.0. StatSoft® Inc.). A MANOVA test was also used to compare treatment effects (queen insemination quantity). As the MANOVA test was significant, we compared each compound using a univariate analysis on the relative proportion for each compound. Both parametric and non-parametric tests were used for individual component comparisons, and both tests revealed comparable statistic p value. Here, we present the results from the Mann-Whitney test.

The following abbreviations are used for chemicals found in the analysis: methyl p-hydroxybenzoate (HOB), 8-hydroxyoctanoitic acid (8-HOAA), 4-hydroxy-3-methoxyphenylethanol (HVA), (E)-9-oxodec-2-enoic acid (ODA), (E)-9-hydroxydec-2-enoic acid (9-HDA), 10-hydroxydecanoic acid (10-HDAA), (E)-10-hydroxydec-2-enoic acid (10-HDA), decanedioic acid (C10∶0 DA) and (E)-dec-2-enedioic acid (C10∶1 DA)

Brain gene expression levels of Amfor

Total RNA was isolated from dissected brains using a RNeasy RNA extraction kit (Qiagen, Valencia CA), yielding 0.6-1.5 µg/brain. cDNA was synthesized from 150 ng RNA using Arrayscript reverse transcriptase (Ambion, CA). Expression levels of Amfor were measured using quantitative real-time RT-PCR (qRT-PCR) with an ABI Prism 7900 sequence detector and the SYBR green detection method (Applied Biosystems, Foster City, CA). eIF3-S8, a housekeeping gene that did not vary in expression levels in previous bee brain microarray studies [67], [72], was used as a loading control. For each sample, triplicate qRT-PCR reactions were performed and averaged. A standard curve was generated for each primer using dilutions of genomic DNA to calculate the relative quantities of mRNA levels for each sample. A dissociation curve and negative control (cDNA reaction without RT-enzyme) were used to ensure primer specificity and lack of contamination.

The sequences for the primers (5′ to 3′) used are as follows:

Amfor -F: AATATAACTTCCGGTGCAACGTATT;

Amfor -R: CGTTTGGATCACGGAAGAAAG;

eIFS8-F: TGAGTGTCTGCTATGGATTGCAA;

eIFS8-R: TCGCGGCTCGTGGTAAA;

We evaluated the brain gene-expression levels of 9 SDI queens and 9 MDI queens (derived from two cohorts). For each individual brain sample, the ratio of the expression level of Amfor to that of the control gene (eIF3-S8) was calculated.

Acknowledgments

We would like to thank Coby Schal and Michael Roe for use of their GC and GC/MS instruments, Cavell Brownie for her advice on the statistical analyses, Josh Summers, Jennifer Keller, and Joe Flowers for their beekeeping assistance, Kelly Hutcherson and Cynthia Rouf for their assistance in the field. This experiment complies with the laws of the USA.

Author Contributions

Conceived and designed the experiments: CG FR DT. Performed the experiments: FR DT. Analyzed the data: CG FR. Contributed reagents/materials/analysis tools: CG DT. Wrote the paper: CG FR DT.

Citation: Richard F-J, Tarpy DR, Grozinger CM (2007) Effects of Insemination Quantity on Honey Bee Queen Physiology. PLoS ONE2(10): e980. https://doi.org/10.1371/journal.pone.0000980

References

1. Gillott C (2003) Male accessory gland secretions: modulators of female reproductive physiology and behavior. Annu Rev Entomol 48: 163–184.C. Gillott2003Male accessory gland secretions: modulators of female reproductive physiology and behavior.Annu Rev Entomol48163184

2. Tarpy DR, Nielsen DI (2002) Sampling error, effective paternity and estimating the genetic structure of honey bee colonies (Hymenoptera: Apidae). Ann Entomol Soc Am 95: 513–528.DR TarpyDI Nielsen2002Sampling error, effective paternity and estimating the genetic structure of honey bee colonies (Hymenoptera: Apidae).Ann Entomol Soc Am95513528

3. Oldroyd BP, Clifton MJ, Parker K, Wongsiri S, Rinderer TE, et al. (1998) Evolution of mating behavior in the genus Apis and an estimate of mating frequency in A. cerana (Hymenoptera: Apidae). Ann Entomol Soc Am 91: 700–709.BP OldroydMJ CliftonK. ParkerS. WongsiriTE Rinderer1998Evolution of mating behavior in the genus Apis and an estimate of mating frequency in A. cerana (Hymenoptera: Apidae).Ann Entomol Soc Am91700709

4. Oldroyd BP, Clifton MJ, Wongsiri S, Rinderer T, Sylvester HA, et al. (1997) Polyandry in the genus Apis, particularly Apis andreniformis. Behav Ecol Sociobiol 40: 17–26.BP OldroydMJ CliftonS. WongsiriT. RindererHA Sylvester1997Polyandry in the genus Apis, particularly Apis andreniformis.Behav Ecol Sociobiol401726

5. Tarpy DR, Seeley TD (2006) Lower disease infections in honeybee (Apis mellifera) colonies headed by polyandrous vs monandrous queens. 93: 195–199.DR TarpyTD Seeley2006Lower disease infections in honeybee (Apis mellifera) colonies headed by polyandrous vs monandrous queens.93195199

6. Crozier RH, Pamilo P (1996) Evolution of social insect colonies;. In: May RM, Harvey PH, editors. Oxford: Oxford University Press. RH CrozierP. Pamilo1996Evolution of social insect colonies;RM MayPH HarveyOxfordOxford University Press306

7. Tarpy DR, Page RE (2000) No behavioral control over mating frequency in queen honey bees (Apis mellifera L.): implications for the evolution of extreme polyandry. Am Nat 155: 820–827.DR TarpyRE Page2000No behavioral control over mating frequency in queen honey bees (Apis mellifera L.): implications for the evolution of extreme polyandry.Am Nat155820827

8. Schlüns H, Moritz RFA, Neumann P, Kryger P, Koeniger G (2005) Multiple nuptial flights, sperm transfer and the evolution of extreme polyandry in honeybee queens. Anim Behav 70: 125–131.H. SchlünsRFA MoritzP. NeumannP. KrygerG. Koeniger2005Multiple nuptial flights, sperm transfer and the evolution of extreme polyandry in honeybee queens.Anim Behav70125131

9. Tanaka ED, Hartfelder K (2004) The initial stages of oogenesis and their relation to differential fertility in the honey bee (Apis mellifera) castes. Arth Struct Dev 33: 431–442.ED TanakaK. Hartfelder2004The initial stages of oogenesis and their relation to differential fertility in the honey bee (Apis mellifera) castes.Arth Struct Dev33431442

10. Winston ML (1987) The Biology of the Honey Bee. Cambridge: Havard University Press. ML Winston1987The Biology of the Honey Bee.CambridgeHavard University Press

11. Harano K, Sasaki K, Nagao T (2005) Depression of brain dopamine and its metabolite after mating in European honeybee (Apis mellifera) queens. Naturwissenschaften 92: 310–313.K. HaranoK. SasakiT. Nagao2005Depression of brain dopamine and its metabolite after mating in European honeybee (Apis mellifera) queens.Naturwissenschaften92310313

12. Fahrbach SE, Giray T, RObinson GE (1995) Volume changes in the mushroom bodies of adult honey bee queens. Neur Learn Mem 63: 181–191.SE FahrbachT. GirayGE RObinson1995Volume changes in the mushroom bodies of adult honey bee queens.Neur Learn Mem63181191

13. Slessor KN, Kaminski LA, King GGS, Winston ML (1990) Semiochemicals of the honeybee queen mandibular glands. J Chem Ecol 16: 851–860.KN SlessorLA KaminskiGGS KingML Winston1990Semiochemicals of the honeybee queen mandibular glands.J Chem Ecol16851860

14. Plettner E, Otis GW, Wimalaratne PDC, Winston ML, Slessor KN, et al. (1997) Species- and caste-determined mandibular gland signals in honeybees (Apis). J Chem Ecol 23: 363–377.E. PlettnerGW OtisPDC WimalaratneML WinstonKN Slessor1997Species- and caste-determined mandibular gland signals in honeybees (Apis).J Chem Ecol23363377

15. Blum MS (1992) The hive and the honey bee. In: Graham JM, editor. Hamilton: Dadant and Sons. MS Blum1992The hive and the honey bee.JM GrahamHamiltonDadant and Sons

16. Free JB (1987) Pheromones of social bees. London: University Press, Chapman and Hall. JB Free1987Pheromones of social bees.LondonUniversity Press, Chapman and Hall

17. Keeling C, Slessor KN, Higo HA, Winston ML (2003) New components of the honey bee (Apis mellifera L.) queen retinue pheromone. Proc Natl Acad Sci USA 100: 4486–4491.C. KeelingKN SlessorHA HigoML Winston2003New components of the honey bee (Apis mellifera L.) queen retinue pheromone.Proc Natl Acad Sci USA10044864491

18. Slessor KN, Kaminski LA, King GGS, Borden JH, Winston ML (1988) The semiochemical basis of the retinue response to queen honey bees. Nature 332: 354–356.KN SlessorLA KaminskiGGS KingJH BordenML Winston1988The semiochemical basis of the retinue response to queen honey bees.Nature332354356

19. Pankiw T, Winston M, Slessor KN (1994) Variation in worker reponse to honey bee (Apis mellifera L.) queen mandibular pheromone (Hymenoptera: Apidae). J Ins Behav 7: 1–15.T. PankiwM. WinstonKN Slessor1994Variation in worker reponse to honey bee (Apis mellifera L.) queen mandibular pheromone (Hymenoptera: Apidae).J Ins Behav7115

20. Pankiw T, Winston M, Slessor KN (1995) Queen attendance behaviour of worker honey bees (Apis mellifera L.) that are high and low responding to queen mandibular pheromone. Insects Soc 42: 371–378.T. PankiwM. WinstonKN Slessor1995Queen attendance behaviour of worker honey bees (Apis mellifera L.) that are high and low responding to queen mandibular pheromone.Insects Soc42371378

21. Kaminski LA, Slessor KN, Winston ML, Hay NW, Borden JH (1990) Honeybee response to queen mandibular pheromone in laboratory bioassays. J Chem Ecol 16: 841–850.LA KaminskiKN SlessorML WinstonNW HayJH Borden1990Honeybee response to queen mandibular pheromone in laboratory bioassays.J Chem Ecol16841850

22. Butler CG, Fairey EM (1963) The role of the queen in preventing oogenesis in worker honey bees. J Apic Res 2: 14–18.CG ButlerEM Fairey1963The role of the queen in preventing oogenesis in worker honey bees.J Apic Res21418

23. Hoover SER, Keeling CI, Winston ML, Slessor KN (2003) The effect of queen pheromones on worker honey bee ovary develpment. Naturwissenschaften 90: 477–480.SER HooverCI KeelingML WinstonKN Slessor2003The effect of queen pheromones on worker honey bee ovary develpment.Naturwissenschaften90477480

24. Melathopoulos AP, Winston M, Pettis JS, Pankiw T (1996) Effect of queen mandibular pheromone on initiation and maintenance of queen cells in the honey bee (Apis mellifera L.). Can Entomol 128: 263–272.AP MelathopoulosM. WinstonJS PettisT. Pankiw1996Effect of queen mandibular pheromone on initiation and maintenance of queen cells in the honey bee (Apis mellifera L.).Can Entomol128263272

25. Pettis JS, Winston M, Collins AM (1995) Supprression of queen rearing in European and Africanized honey bees Apis melifera L. by synthetic queen mandibular gland pheromone. Insectes Soc 42: 113–121.JS PettisM. WinstonAM Collins1995Supprression of queen rearing in European and Africanized honey bees Apis melifera L. by synthetic queen mandibular gland pheromone.Insectes Soc42113121

26. Higo HA, Colley SJ, Winston ML, Slessor KN (1992) Effects of honey bee (Apis mellifera L.) queen mandibular gland pheromone on foraging and brood rearing. Can Entomol 124: 409–418.HA HigoSJ ColleyML WinstonKN Slessor1992Effects of honey bee (Apis mellifera L.) queen mandibular gland pheromone on foraging and brood rearing.Can Entomol124409418

27. Pankiw T, Huang ZY, Winston ML, Robinson GE (1998) Queen mandibular gland pheromone influences worker honey bee (Apis mellifera L.) foraging ontogeny and juvenile hormone titers. J Ins Physiol 44: 685–692.T. PankiwZY HuangML WinstonGE Robinson1998Queen mandibular gland pheromone influences worker honey bee (Apis mellifera L.) foraging ontogeny and juvenile hormone titers.J Ins Physiol44685692

28. Pettis JS, Winston M, Slessor KN (1995) Behavior of queen and worker honey bees (Hymenoptera, Apidae) in response to exogenous queen mandibular gland pheromone. Ann Entomol Soc Amer 88: 580–588.JS PettisM. WinstonKN Slessor1995Behavior of queen and worker honey bees (Hymenoptera, Apidae) in response to exogenous queen mandibular gland pheromone.Ann Entomol Soc Amer88580588

29. Gary NE (1962) Chemical mating attractants in the queen honey bee. Science 136: 773–774.NE Gary1962Chemical mating attractants in the queen honey bee.Science136773774

30. Brockmann A, Dietz D, Spaethe J, Tautz J (2006) Beyond 9-ODA: sex pheromone communication in the european honey bee Apis mellifera L. J Chem Ecol 32: 657–667.A. BrockmannD. DietzJ. SpaetheJ. Tautz2006Beyond 9-ODA: sex pheromone communication in the european honey bee Apis mellifera L.J Chem Ecol32657667

31. Apsegaite V, Skirkevicius A (1999) Content of (E)-9-oxo-2-decenoic acid in pheromones of honeybee (Apis mellifera L.) queens. Pheromones 9: 27–32.V. ApsegaiteA. Skirkevicius1999Content of (E)-9-oxo-2-decenoic acid in pheromones of honeybee (Apis mellifera L.) queens.Pheromones92732

32. Laidlaw HHj (1977) Instrumental inseminantion of honey bee queens: its origin and development. Hamilton, IL: Dandant. HHj Laidlaw1977Instrumental inseminantion of honey bee queens: its origin and development.Hamilton, ILDandant

33. Woyke J, Jasinski A, Fliszkiewicz (1995) Further investigations on natural mating of instrumentally inseminated Apis mellifera queens. J Apic Res 32: 105–106.J. WoykeA. JasinskiFliszkiewicz1995Further investigations on natural mating of instrumentally inseminated Apis mellifera queens.J Apic Res32105106

34. Ben-Shahar Y (2005) The foraging gene, behavioral plasticity, and honeybee division of labor. J Comp Physiol A 191: 987–994.Y. Ben-Shahar2005The foraging gene, behavioral plasticity, and honeybee division of labor.J Comp Physiol A191987994

35. Ben-Shahar Y, Robichon A, Sokolowski MB, Robinson GE (2002) Behavior influenced by gene actio across different time scales. Science 296: 742–744.Y. Ben-ShaharA. RobichonMB SokolowskiGE Robinson2002Behavior influenced by gene actio across different time scales.Science296742744

36. Ben-Shahar Y, Leung HT, Pak WL, Sokolowski MB, Robinson GE (2003) cGMP-dependent changes in phototaxis: a possible role for the foraging gene in honey bee division of labor. J Exp Biol 14: 2507–2515.Y. Ben-ShaharHT LeungWL PakMB SokolowskiGE Robinson2003cGMP-dependent changes in phototaxis: a possible role for the foraging gene in honey bee division of labor.J Exp Biol1425072515

37. Woyke J (1966) Which is the most suitable volume of semen in artificial inseminantion of queen honey bees. USDA Bee Cult. Lab. Warsaw 12–15.J. Woyke1966Which is the most suitable volume of semen in artificial inseminantion of queen honey bees.USDA Bee Cult. Lab. Warsaw1215

38. Woyke J (1983) Dynamics of entry of spermatozoa into the spermatheca of instrumentally inseminated queen honeybees. J Apic Res 22: 150–154.J. Woyke1983Dynamics of entry of spermatozoa into the spermatheca of instrumentally inseminated queen honeybees.J Apic Res22150154

39. Dade HA (1978) Anatomy and dissection of the honeybee. London: International Bee Research Association. HA Dade1978Anatomy and dissection of the honeybee.LondonInternational Bee Research Association158

40. Tomanek B, Jasinski A, Sulek Z, Muszynska J, Kulinowski Pea (1996) Magnetic resonance microscopy of internal struction of drone and queen honey bees. J Apic Res 35: 3–9.B. TomanekA. JasinskiZ. SulekJ. MuszynskaPea Kulinowski1996Magnetic resonance microscopy of internal struction of drone and queen honey bees.J Apic Res3539

41. Nijhout HF (1984) Abdominal stretch reception in Dipetalogaster-Maximus (Hemiptera, Reduviidae). J Ins Physiol 30: 629–633.HF Nijhout1984Abdominal stretch reception in Dipetalogaster-Maximus (Hemiptera, Reduviidae).J Ins Physiol30629633

42. Ringo J (1996) Sexual receptivity in insects. Annu Rev Entomol 41: 473–494.J. Ringo1996Sexual receptivity in insects.Annu Rev Entomol41473494

43. Roth LM (1962) Hypersexual activity induced in females of the cockroach Nauphoeta cinerea. Science 138: 1267–1269.LM Roth1962Hypersexual activity induced in females of the cockroach Nauphoeta cinerea.Science13812671269

44. Kaftanoglu O, Peng YS (1982) Effects of inseminantion on the initiation of oviposition in the queen honey bee. J Apic Res 21: 3–6.O. KaftanogluYS Peng1982Effects of inseminantion on the initiation of oviposition in the queen honey bee.J Apic Res2136

45. Chapman T, Davies SJ (2004) Functions and analysis of the seminal fluid proteins of male Drosophila melanogaster fruit flies. Peptides 25: 1477–1490.T. ChapmanSJ Davies2004Functions and analysis of the seminal fluid proteins of male Drosophila melanogaster fruit flies.Peptides2514771490

46. Heifetz Y, Lung O, Frongillo EA, Wolfner MF (2000) The drosophila seminal fluid protein Acp26Aa stimulates release of oocytes by the ovary. Current Biology 10: 99–102.Y. HeifetzO. LungEA FrongilloMF Wolfner2000The drosophila seminal fluid protein Acp26Aa stimulates release of oocytes by the ovary.Current Biology1099102

47. Mueller JL, Ripoll DR, Aquadro CF, Wolfner MF (2004) Comparative structural modeling and inference of conserved protein classes in Drosophila seminal fluid. PNAS 101: 13542–13547.JL MuellerDR RipollCF AquadroMF Wolfner2004Comparative structural modeling and inference of conserved protein classes in Drosophila seminal fluid.PNAS1011354213547

48. Riddiford LM, Ashenhurst JB (1973) The switchover from virgin to mated behavior in female Cecropia moths: the role of the Bursa copulatrix. Biol Bull 144: 162–171.LM RiddifordJB Ashenhurst1973The switchover from virgin to mated behavior in female Cecropia moths: the role of the Bursa copulatrix.Biol Bull144162171

49. Sugawara T (1986) Oviposition behavior of the cricket Teleogryllus commodus: the site of oction of an oviposition-stimulating factor and the role of the nervous system. J Ins Physiol 32: 485–492.T. Sugawara1986Oviposition behavior of the cricket Teleogryllus commodus: the site of oction of an oviposition-stimulating factor and the role of the nervous system.J Ins Physiol32485492

50. Tram U, Wolfner MF (1998) Seminal fluid regulation of female sexual attractiveness in Drosophila melanogaster. PNAS 95: 4051–4054.U. TramMF Wolfner1998Seminal fluid regulation of female sexual attractiveness in Drosophila melanogaster.PNAS9540514054

51. Wolfner MF (1997) Tokens of love: functions and regulation of Drosophila male accessory gland products. Insect Biochem Molec Biol 27: 179–192.MF Wolfner1997Tokens of love: functions and regulation of Drosophila male accessory gland products.Insect Biochem Molec Biol27179192

52. Wolfner MF (2002) THe gifts that keep on giving: physiological functions and evolutionary dynamics of male seminal proteins in Drosophila. Heredity 88: 85–93.MF Wolfner2002THe gifts that keep on giving: physiological functions and evolutionary dynamics of male seminal proteins in Drosophila.Heredity888593

53. Baer B, Maile R, Schmid-Hempel P, Morgan ED, Jones GR (2000) Chemistry of a mating plug in bumblebees. J Chem Ecol 26: 1869–1875.B. BaerR. MaileP. Schmid-HempelED MorganGR Jones2000Chemistry of a mating plug in bumblebees.J Chem Ecol2618691875

54. Baer B, Morgan ED, Schmid-Hempel P (2001) A nonspecific fatty acid within the bumblebee mating plug prevents females from remating. PNAS 98: 3926–3928.B. BaerED MorganP. Schmid-Hempel2001A nonspecific fatty acid within the bumblebee mating plug prevents females from remating.PNAS9839263928

55. Blum MS, Glowska ZS, Taber I (1962) Chemistry of the drone honey bee repoductive system. II. Carbohydrates in the reproductive organs and semen. Ann Entomol Soc Am 55: 135–139.MS BlumZS GlowskaI. Taber1962Chemistry of the drone honey bee repoductive system. II. Carbohydrates in the reproductive organs and semen.Ann Entomol Soc Am55135139

56. Collins AM, Williams V, Evans JD (2004) Sperm storage and antioxidative enzyme expression in the honey bee, Apis mellifera. Ins Molec Biol 13: 141–146.AM CollinsV. WilliamsJD Evans2004Sperm storage and antioxidative enzyme expression in the honey bee, Apis mellifera.Ins Molec Biol13141146

57. Klenk M, Koeniger G, Koeniger N, Fasold H (2004) Proteins in spermathecal gland secretion and spermathecal fluid and the properties of a 29 kDa protein in queens of Apis mellifera. Apidologie 35: 371–381.M. KlenkG. KoenigerN. KoenigerH. Fasold2004Proteins in spermathecal gland secretion and spermathecal fluid and the properties of a 29 kDa protein in queens of Apis mellifera.Apidologie35371381

58. Seeley TD, Tarpy DR (2007) Queen promiscuity lowers disease within honeybee colonies. Proc R Soc B 274: 67–72.TD SeeleyDR Tarpy2007Queen promiscuity lowers disease within honeybee colonies.Proc R Soc B2746772

59. Tarpy DR, Page RE (2002) Sex determinantion and the evolution of polyandry in honey bees (Apis mellifera). Behav Ecol Sociobiol 52: 143–150.DR TarpyRE Page2002Sex determinantion and the evolution of polyandry in honey bees (Apis mellifera).Behav Ecol Sociobiol52143150

60. Lobo JA, Kerr WE (1993) Estimation of the number of matings in Apis mellifera: extentsions of the model and comparison of different estimates. Ethol Ecol Evol 5: 337–345.JA LoboWE Kerr1993Estimation of the number of matings in Apis mellifera: extentsions of the model and comparison of different estimates.Ethol Ecol Evol5337345

61. Oldroyd BP, Rinderer TE, Harbo JR, Buco SM (1992) Effects of intracolonial genetic diversity on honey bee (Hymenoptera: Apidae) colony performance. Ann Entomol Soc Am 85: 335–343.BP OldroydTE RindererJR HarboSM Buco1992Effects of intracolonial genetic diversity on honey bee (Hymenoptera: Apidae) colony performance.Ann Entomol Soc Am85335343

62. Keller L, Nonacs P (1993) The role of queen pheromones in social insects: queen control or queen signal? Anim Behav 45: 787–794.L. KellerP. Nonacs1993The role of queen pheromones in social insects: queen control or queen signal?Anim Behav45787794

63. Seeley T (1985) Honeybee ecology: a study of adaptation in social life. Princeton, New Jersey: Princeton University Press. T. Seeley1985Honeybee ecology: a study of adaptation in social life.Princeton, New JerseyPrinceton University Press201

64. Woyciechowskii M, Lomnicki A (1987) Multiple mating of queens and the sterility of workers among eusocial hymenoptera. J Theor Biol 128: 317–328.M. WoyciechowskiiA. Lomnicki1987Multiple mating of queens and the sterility of workers among eusocial hymenoptera.J Theor Biol128317328

65. Pankiw T, Winston ML, Fondrk MK, Slessor KN (2000) Selection on worker honeybee responses to queen pheromone (Apis mellifera L.). Naturwissenschaften 87: 487–490.T. PankiwML WinstonMK FondrkKN Slessor2000Selection on worker honeybee responses to queen pheromone (Apis mellifera L.).Naturwissenschaften87487490

66. Laidlaw HHj, Page REj (1997) Queen rearing and bee breeding;. In: Cheschire CW, editor. HHj LaidlawREj Page1997Queen rearing and bee breeding;CW Cheschire

67. Grozinger CM, Sharabash NM, Whitfield CW, Robinson GE (2003) Pheromone-mediated gene expression in the honey bee brain. Proc Natl Acad Sci USA 100: 14519–14525.CM GrozingerNM SharabashCW WhitfieldGE Robinson2003Pheromone-mediated gene expression in the honey bee brain.Proc Natl Acad Sci USA1001451914525

68. Grozinger CM, Fischer P, Hampton JE (2006) Uncoupling primer and releaser responses to pheromone in honey bees. Naturwissenschaften 94: 375–379.CM GrozingerP. FischerJE Hampton2006Uncoupling primer and releaser responses to pheromone in honey bees.Naturwissenschaften94375379

69. Hoover SER, Oldroyd BP, Wossler TC, Winston ML (2005) Anarchistic queen honey bees have normal queen mandibular pheromones. Insectes Soc 52: 6–10.SER HooverBP OldroydTC WosslerML Winston2005Anarchistic queen honey bees have normal queen mandibular pheromones.Insectes Soc52610

70. Slessor KN, Winston M, Le Conte Y (2005) Pheromone communication in the honey bee (Apis mellifera L.). J Chem Ecol 31: 2731–2745.KN SlessorM. WinstonY. Le Conte2005Pheromone communication in the honey bee (Apis mellifera L.).J Chem Ecol3127312745

71. Morgan SM, Butz Huryn VM, Downes SR, Mercer AR (1998) The effects of queenlessness on the maturation of the honey bee olfactory system. Behav Brain Res 91: 115–126.SM MorganVM Butz HurynSR DownesAR Mercer1998The effects of queenlessness on the maturation of the honey bee olfactory system.Behav Brain Res91115126

72. Whitfield CW, Cziko AM, Robinson G (2003) Gene expression profiles in the brain predict behavior in individual Honey Bees. Science 302: 296–299.CW WhitfieldAM CzikoG. Robinson2003Gene expression profiles in the brain predict behavior in individual Honey Bees.Science302296299