# Effects of Insemination Quantity on Honey Bee Queen Physiology

Freddie-Jeanne Richard, David R Tarpy, Christina M Grozinger

► To cite this version:

Freddie-Jeanne Richard, David R Tarpy, Christina M Grozinger. Effects of Insemination Quantity on Honey Bee Queen Physiology. PLoS ONE, 2007, 2 (10), pp.e980. 10.1371/journal.pone.0000980. hal-05285604

HAL Id: hal-05285604

https://hal.science/hal-05285604v1

Submitted on 26 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Effects of Insemination Quantity on Honey Bee Queen Physiology

Freddie-Jeanne Richard¹,²,³, David R. Tarpy¹,³, Christina M. Grozinger¹,²,³*

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

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.

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

## 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

Academic Editor: Robert Brooks, The University of New South Wales, Australia

Received July 17, 2007; Accepted September 12, 2007; Published October 3, 2007

**Copyright:** © 2007 Richard et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by a USDA-NRI grant to CMG and DRT (2006-35607-16625) and the Eastern Apiculture Society.

**Competing Interests:** The authors have declared that no competing interests exist.

* To whom correspondence should be addressed. E-mail: cmgrozin@ncsu.edu

<figure><img src="image_2.png" /></figure>
(9-HDA), methyl *p*-hydroxybenzoate (HOB), and 4-hydroxy-3-methoxyphenylethanol (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-old 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).

**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.

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 chroma-
tography. A total of 27 compounds were apparent in the GC
analysis and used in the subsequent comparisons. There were

<figure><img src="image_3.png" /><figcaption>**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).</figcaption></figure>

doi:10.1371/journal.pone.0000980.g001

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>**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</figcaption></figure>

doi:10.1371/journal.pone.0000980.g002

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⁻⁴; 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.

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.

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.

*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).

<caption>**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⁻⁴; All group distances <0.005). Ellipses have been drawn to emphasize the categories, but have no specific statistical meaning</caption>

doi:10.1371/journal.pone.0000980.g003

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 inter-

<figure><img src="image_6.png" /></figure>
<table><caption>Table 1. Absolute quantity of the mandibular gland compounds (mean±STD) in µg of virgin, SDI and MDI queens.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Retention Index</th><th colspan="2">Virgin (N = 12)</th><th colspan="2">Inseminated (n = 25)</th><th colspan="2">Mann Whitney (Virgin vs Inseminated)</th><th colspan="2">SDI (n = 15)</th><th colspan="2">MDI (n = 10)</th><th colspan="2">Mann Whitney (SDI vs MDI)</th></tr><tr><th>mean</th><th>STD</th><th>mean</th><th>STD</th><th>U</th><th>p</th><th>mean</th><th>STD</th><th>mean</th><th>STD</th><th>U</th><th>p</th></tr></thead><tbody><tr><td>Unidentified 1</td><td>1408</td><td>0.70</td><td>±0.09</td><td>0.70</td><td>±0.04</td><td>143</td><td>0.835</td><td>0.69</td><td>±0.04</td><td>0.71</td><td>±0.09</td><td>69</td><td>0.765</td></tr><tr><td>HOB</td><td>1499</td><td>0.20</td><td>±0.01</td><td>0.20</td><td>±0.02</td><td>126</td><td>0.451</td><td>0.22</td><td>±0.03</td><td>0.16</td><td>±0.04</td><td>49</td><td>0.160</td></tr><tr><td>4-hydroxy-3methoxybenzoic acid</td><td>1512</td><td>4.58</td><td>±0.44</td><td>2.32</td><td>±0.20</td><td>31</td><td>0.000</td><td>2.70</td><td>±0.21</td><td>1.74</td><td>±0.33</td><td>34</td><td>0.023</td></tr><tr><td>8-HOAA</td><td>1627</td><td>12.12</td><td>±1.08</td><td>6.19</td><td>±0.80</td><td>42</td><td>0.000</td><td>7.71</td><td>±0.96</td><td>3.92</td><td>±1.06</td><td>35</td><td>0.026</td></tr><tr><td>4-hydroxybenzoic acid</td><td>1635</td><td>8.09</td><td>±0.64</td><td>3.77</td><td>±0.28</td><td>14</td><td>0.000</td><td>3.92</td><td>±0.29</td><td>3.55</td><td>±0.55</td><td>63</td><td>0.531</td></tr><tr><td>3-HDAA</td><td>1658</td><td>1.07</td><td>±0.07</td><td>0.53</td><td>±0.05</td><td>18</td><td>0.000</td><td>0.59</td><td>±0.05</td><td>0.43</td><td>±0.08</td><td>42</td><td>0.071</td></tr><tr><td>4-hydroxydihydrocinnamyl alcohol</td><td>1687</td><td>0.07</td><td>±0.02</td><td>0.22</td><td>±0.04</td><td>89</td><td>0.049</td><td>0.22</td><td>±0.05</td><td>0.22</td><td>±0.09</td><td>61</td><td>0.461</td></tr><tr><td>9-ODA</td><td>1714</td><td>252.45</td><td>±20.15</td><td>125.16</td><td>±11.30</td><td>23</td><td>0.000</td><td>144.70</td><td>±13.18</td><td>95.84</td><td>±16.94</td><td>35</td><td>0.026</td></tr><tr><td>HVA</td><td>1716</td><td>0.51</td><td>±0.04</td><td>0.29</td><td>±0.02</td><td>24</td><td>0.000</td><td>0.32</td><td>±0.02</td><td>0.24</td><td>±0.03</td><td>46</td><td>0.115</td></tr><tr><td>9-oxodecanoic acid</td><td>1750</td><td>0.77</td><td>±0.09</td><td>0.44</td><td>±0.04</td><td>53</td><td>0.001</td><td>0.51</td><td>±0.05</td><td>0.34</td><td>±0.04</td><td>32</td><td>0.016</td></tr><tr><td>4-hydroxy-3methoxybenzoic acid</td><td>1785</td><td>2.66</td><td>±0.21</td><td>1.48</td><td>±0.13</td><td>33</td><td>0.000</td><td>1.73</td><td>±0.14</td><td>1.12</td><td>±0.19</td><td>33</td><td>0.019</td></tr><tr><td>9-HDA</td><td>1798</td><td>37.65</td><td>±3.87</td><td>19.56</td><td>±1.68</td><td>35</td><td>0.000</td><td>22.91</td><td>±2.04</td><td>14.54</td><td>±2.12</td><td>32</td><td>0.016</td></tr><tr><td>Unidentified 2</td><td>1817</td><td>0.91</td><td>±0.09</td><td>0.63</td><td>±0.05</td><td>66</td><td>0.006</td><td>0.66</td><td>±0.07</td><td>0.58</td><td>±0.09</td><td>63</td><td>0.531</td></tr><tr><td>10-HDAA</td><td>1820</td><td>2.84</td><td>±0.23</td><td>1.57</td><td>±0.13</td><td>28</td><td>0.000</td><td>1.75</td><td>±0.15</td><td>1.30</td><td>±0.21</td><td>47</td><td>0.129</td></tr><tr><td>10-HDA</td><td>1873</td><td>34.87</td><td>±3.95</td><td>23.90</td><td>±2.63</td><td>79</td><td>0.021</td><td>24.64</td><td>±3.57</td><td>22.81</td><td>±4.04</td><td>73</td><td>0.935</td></tr><tr><td>Unidentified 3</td><td>1885</td><td>0.49</td><td>±0.04</td><td>0.32</td><td>±0.03</td><td>45</td><td>0.000</td><td>0.33</td><td>±0.02</td><td>0.29</td><td>±0.07</td><td>56</td><td>0.311</td></tr><tr><td>C10:0 DA</td><td>1907</td><td>0.17</td><td>±0.01</td><td>0.16</td><td>±0.02</td><td>130</td><td>0.532</td><td>0.19</td><td>±0.02</td><td>0.13</td><td>±0.02</td><td>47</td><td>0.129</td></tr><tr><td>Dihydroferulic acid</td><td>1907</td><td>0.30</td><td>±0.02</td><td>0.16</td><td>±0.01</td><td>19</td><td>0.000</td><td>0.18</td><td>±0.01</td><td>0.13</td><td>±0.02</td><td>41</td><td>0.062</td></tr><tr><td>Alkane 1</td><td>1918</td><td>2.17</td><td>±0.15</td><td>1.22</td><td>±0.10</td><td>31</td><td>0.000</td><td>1.40</td><td>±0.11</td><td>0.95</td><td>±0.17</td><td>65</td><td>0.605</td></tr><tr><td>C10:1 DA</td><td>1957</td><td>3.29</td><td>±0.32</td><td>2.35</td><td>±0.20</td><td>71</td><td>0.009</td><td>2.48</td><td>±0.23</td><td>2.16</td><td>±0.36</td><td>40</td><td>0.055</td></tr><tr><td>Unidentified 4</td><td>1986</td><td>2.86</td><td>±0.31</td><td>1.36</td><td>±0.16</td><td>43</td><td>0.000</td><td>1.59</td><td>±0.21</td><td>1.03</td><td>±0.24</td><td>28</td><td>0.008</td></tr><tr><td>Unidentified 5</td><td>2023</td><td>6.88</td><td>±0.68</td><td>4.31</td><td>±0.45</td><td>62</td><td>0.003</td><td>5.23</td><td>±0.52</td><td>2.92</td><td>±0.61</td><td>64</td><td>0.567</td></tr><tr><td>Palmitic acid</td><td>2052</td><td>0.79</td><td>±0.04</td><td>0.42</td><td>±0.03</td><td>14</td><td>0.000</td><td>0.39</td><td>±0.03</td><td>0.45</td><td>±0.06</td><td>45</td><td>0.103</td></tr><tr><td>Unidentified 6</td><td>2068</td><td>4.74</td><td>±0.40</td><td>3.48</td><td>±0.33</td><td>74</td><td>0.013</td><td>3.91</td><td>±0.43</td><td>2.83</td><td>±0.45</td><td>59</td><td>0.397</td></tr><tr><td>Alkane 2</td><td>2119</td><td>0.36</td><td>±0.03</td><td>0.27</td><td>±0.02</td><td>60</td><td>0.003</td><td>0.25</td><td>±0.02</td><td>0.29</td><td>±0.05</td><td>66</td><td>0.643</td></tr><tr><td>Octadecenoic acid C18:1</td><td>2022</td><td>7.55</td><td>±0.34</td><td>2.60</td><td>±0.35</td><td>10</td><td>0.000</td><td>2.94</td><td>±0.52</td><td>2.08</td><td>±0.35</td><td>69</td><td>0.765</td></tr></tbody></table>

All significant differences with p<0.05 are marked in bold.

doi:10.1371/journal.pone.0000980.t001

actions. 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

<figure><img src="image_7.png" /></figure>
<table><caption>**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.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Retention Index</th><th colspan="2">Virgin (N = 12)</th><th colspan="2">Inseminated (n = 25)</th><th colspan="2">Mann Whitney (Virgin vs Inseminated)</th><th colspan="2">SDI (n = 15)</th><th colspan="2">MDI (n = 10)</th><th colspan="2">Mann Whitney (SDI vs MDI)</th></tr><tr><th>mean</th><th>STD</th><th>mean</th><th>STD</th><th>U</th><th>p</th><th>mean</th><th>STD</th><th>mean</th><th>STD</th><th>U</th><th>p</th></tr></thead><tbody><tr><td>Unidentified 1</td><td>1408</td><td>0.21</td><td>±0.05</td><td>0.42</td><td>±0.06</td><td>41</td><td><strong>0.000</strong></td><td>0.33</td><td>±0.03</td><td>0.57</td><td>±0.13</td><td>49</td><td>0.149</td></tr><tr><td>HOB</td><td>1499</td><td>0.06</td><td>±0.01</td><td>0.10</td><td>±0.01</td><td>50</td><td><strong>0.001</strong></td><td>0.10</td><td>±0.01</td><td>0.10</td><td>±0.02</td><td>74</td><td>0.956</td></tr><tr><td>4-hydroxy-3-methoxybenzoic acid</td><td>1512</td><td>1.19</td><td>±0.07</td><td>1.15</td><td>±0.08</td><td>124</td><td>0.399</td><td>1.24</td><td>±0.12</td><td>1.05</td><td>±0.09</td><td>62</td><td>0.471</td></tr><tr><td>8-HOAA</td><td>1627</td><td>3.08</td><td>±0.12</td><td>2.79</td><td>±0.26</td><td>113</td><td>0.230</td><td>3.39</td><td>±0.29</td><td>2.01</td><td>±0.36</td><td>31</td><td><strong>0.015</strong></td></tr><tr><td>4-hydroxybenzoic acid</td><td>1635</td><td>2.11</td><td>±0.10</td><td>1.97</td><td>±0.12</td><td>110</td><td>0.194</td><td>1.78</td><td>±0.10</td><td>2.31</td><td>±0.21</td><td>35</td><td><strong>0.027</strong></td></tr><tr><td>3-HDAA</td><td>1658</td><td>0.28</td><td>±0.02</td><td>0.28</td><td>±0.02</td><td>113</td><td>0.230</td><td>0.26</td><td>±0.03</td><td>0.29</td><td>±0.04</td><td>72</td><td>0.868</td></tr><tr><td>4-hydroxydihydrocinnamyl alcohol</td><td>1687</td><td>0.02</td><td>±0.00</td><td>0.12</td><td>±0.03</td><td>37</td><td><strong>0.000</strong></td><td>0.10</td><td>±0.02</td><td>0.16</td><td>±0.07</td><td>70</td><td>0.782</td></tr><tr><td>9-ODA</td><td>1714</td><td>64.77</td><td>±0.63</td><td>60.38</td><td>±0.96</td><td>73</td><td><strong>0.012</strong></td><td>61.19</td><td>±1.27</td><td>58.65</td><td>±1.36</td><td>49</td><td>0.149</td></tr><tr><td>HVA</td><td>1716</td><td>0.13</td><td>±0.01</td><td>0.15</td><td>±0.01</td><td>106</td><td>0.153</td><td>0.14</td><td>±0.00</td><td>0.16</td><td>±0.01</td><td>50</td><td>0.166</td></tr><tr><td>9-oxodecanoic acid</td><td>1750</td><td>0.19</td><td>±0.01</td><td>0.23</td><td>±0.01</td><td>84</td><td><strong>0.032</strong></td><td>0.22</td><td>±0.01</td><td>0.24</td><td>±0.03</td><td>71</td><td>0.824</td></tr><tr><td>4-hydroxy-3-methoxybenzoic acid</td><td>1785</td><td>0.68</td><td>±0.01</td><td>0.72</td><td>±0.02</td><td>104</td><td>0.136</td><td>0.76</td><td>±0.02</td><td>0.68</td><td>±0.03</td><td>47</td><td>0.120</td></tr><tr><td>9-HDA</td><td>1798</td><td>9.41</td><td>±0.40</td><td>9.70</td><td>±0.28</td><td>136</td><td>0.650</td><td>9.85</td><td>±0.34</td><td>9.43</td><td>±0.50</td><td>63</td><td>0.506</td></tr><tr><td>Unidentified 2</td><td>1817</td><td>0.23</td><td>±0.01</td><td>0.32</td><td>±0.02</td><td>52</td><td><strong>0.001</strong></td><td>0.29</td><td>±0.02</td><td>0.37</td><td>±0.03</td><td>31</td><td><strong>0.015</strong></td></tr><tr><td>10-HDAA</td><td>1820</td><td>0.72</td><td>±0.02</td><td>0.77</td><td>±0.03</td><td>128</td><td>0.475</td><td>0.77</td><td>±0.04</td><td>0.80</td><td>±0.06</td><td>63</td><td>0.506</td></tr><tr><td>10-HDA</td><td>1873</td><td>8.62</td><td>±0.56</td><td>12.04</td><td>±1.01</td><td>95</td><td>0.074</td><td>10.96</td><td>±1.28</td><td>14.11</td><td>±1.48</td><td>41</td><td>0.059</td></tr><tr><td>Unidentified 3</td><td>1885</td><td>0.13</td><td>±0.00</td><td>0.16</td><td>±0.02</td><td>67</td><td><strong>0.007</strong></td><td>0.15</td><td>±0.00</td><td>0.19</td><td>±0.04</td><td>60</td><td>0.405</td></tr><tr><td>C10:0 DA</td><td>1907</td><td>0.04</td><td>±0.00</td><td>0.08</td><td>±0.01</td><td>38</td><td><strong>0.000</strong></td><td>0.09</td><td>±0.01</td><td>0.08</td><td>±0.01</td><td>74</td><td>0.956</td></tr><tr><td>Dihydroferulic acid</td><td>1907</td><td>0.08</td><td>±0.00</td><td>0.08</td><td>±0.00</td><td>114</td><td>0.243</td><td>0.08</td><td>±0.00</td><td>0.08</td><td>±0.00</td><td>53</td><td>0.222</td></tr><tr><td>Alkane 1</td><td>1918</td><td>0.57</td><td>±0.02</td><td>0.60</td><td>±0.01</td><td>113</td><td>0.230</td><td>0.61</td><td>±0.02</td><td>0.57</td><td>±0.02</td><td>54</td><td>0.244</td></tr><tr><td>C10:1 DA</td><td>1957</td><td>0.84</td><td>±0.03</td><td>1.17</td><td>±0.06</td><td>32</td><td><strong>0.000</strong></td><td>1.09</td><td>±0.06</td><td>1.32</td><td>±0.10</td><td>37</td><td><strong>0.035</strong></td></tr><tr><td>Unidentified 4</td><td>1986</td><td>0.71</td><td>±0.04</td><td>0.68</td><td>±0.06</td><td>118</td><td>0.299</td><td>0.69</td><td>±0.09</td><td>0.63</td><td>±0.06</td><td>73</td><td>0.912</td></tr><tr><td>Unidentified 5</td><td>2023</td><td>1.75</td><td>±0.07</td><td>2.05</td><td>±0.11</td><td>99</td><td>0.098</td><td>2.31</td><td>±0.15</td><td>1.72</td><td>±0.11</td><td>30</td><td><strong>0.013</strong></td></tr><tr><td>Palmitic acid</td><td>2052</td><td>0.22</td><td>±0.02</td><td>0.26</td><td>±0.04</td><td>130</td><td>0.516</td><td>0.18</td><td>±0.02</td><td>0.37</td><td>±0.09</td><td>39</td><td><strong>0.046</strong></td></tr><tr><td>Unidentified 6</td><td>2068</td><td>1.21</td><td>±0.04</td><td>1.78</td><td>±0.11</td><td>59</td><td><strong>0.003</strong></td><td>1.76</td><td>±0.16</td><td>1.86</td><td>±0.16</td><td>62</td><td>0.471</td></tr><tr><td>Alkane 2</td><td>2119</td><td>0.10</td><td>±0.01</td><td>0.15</td><td>±0.03</td><td>72</td><td><strong>0.011</strong></td><td>0.12</td><td>±0.01</td><td>0.21</td><td>±0.06</td><td>35</td><td><strong>0.027</strong></td></tr><tr><td>Octadecenoic acid C18:1</td><td>2022</td><td>2.11</td><td>±0.23</td><td>1.42</td><td>±0.19</td><td>64</td><td><strong>0.005</strong></td><td>1.20</td><td>±0.26</td><td>1.50</td><td>±0.29</td><td>58</td><td>0.346</td></tr><tr><td>Stearic acid</td><td>2249</td><td>0.54</td><td>±0.05</td><td>0.45</td><td>±0.05</td><td>85</td><td><strong>0.035</strong></td><td>0.34</td><td>±0.06</td><td>0.56</td><td>±0.09</td><td>33</td><td><strong>0.020</strong></td></tr></tbody></table>

All significant differences with p<0.05 are marked in bold. (E)-4-hydroxycinnamic acid (1948), (E)-coniferyl alcohol (1957), traces doi:10.1371/journal.pone.0000980.t002

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].

<figure><img src="image_8.png" /></figure>
<figure><img src="image_9.png" /><figcaption>**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. doi:10.1371/journal.pone.0000980.g004</figcaption></figure>

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 in-
seminated 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
camica*, 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 in-
semination were treated with CO₂ for 4.0 minutes [66], while virgins
were handled but left untreated, since we wanted a baseline
nonreproductive control group and CO₂ treatment accelerates the
transition to egg-laying [66]. Seven days after emergence, queens
were again treated with CO₂ 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.

<figure><img src="image_10.png" /></figure>
Prior to dissection, queen heads were removed and partially
lyophilized to facilitate dissection [67], after which their mandib-
ular 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*10⁴±25*10⁴ sperm) contained significantly more
sperm than those of SDI queens (15*10⁴±8*10⁴ 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 sylated 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

<figure><img src="image_11.png" /></figure>
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-hydroxyoctanoic 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.

## REFERENCES

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

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.

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.

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.

5. Tarpy DR, Seeley TD (2006) Lower disease infections in honeybee (*Apis mellifera*) colonies headed by polyandrous vs monandrous queens. *JQ 93* 195-199.

6. Crozier RH, Pamilo P (1996) Evolution of social insect colonies; May RM, Harvey PH, eds. Oxford: Oxford University Press. 306 p.

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.

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.

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.

10. Winston ML (1987) The Biology of the Honey Bee. Cambridge: Havard 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.

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.

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

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.

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

16. Free JB (1987) Pheromones of social bees. London: University 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.

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.

19. Pankiw T, Winston M, Slessor KN (1994) Variation in worker response to honey bee (*Apis mellifera* L.) queen mandibular pheromone (Hymenoptera: Apidae). J Ins Behav 7: 1-15.

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.

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.

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

23. Hoover SER, Keeling CI, Winston ML, Slessor KN (2003) The effect of queen pheromones on worker honey bee ovary development. Naturwissenschaften 90: 477-480.

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.

25. Pettis JS, Winston M, Collins AM (1995) Suppression of queen rearing in European and Africanized honey bees *Apis mellifera* L. by synthetic queen mandibular gland pheromone. Insectes Soc 42: 113-121.

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.

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.

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.

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

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.

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

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

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

<figure><img src="image_12.png" /></figure>
34. Ben-Shahar Y (2005) The *foraging* gene, behavioral plasticity, and honeybee division of labor. *J Comp Physiol A* 191: 987-994.

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

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.

37. Woyke J (1966) Which is the most suitable volume of semen in artificial insemination of queen honey bees. USDA Bee Cult. Lab. Warsaw. pp 12-15.

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

39. Dade HA (1978) Anatomy and dissection of the honeybee. London: International Bee Research Association. 158 p.

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.

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

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

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

44. Kaftanoglu O, Peng YS (1982) Effects of insemination on the initiation of oviposition in the queen honey bee. *J Apic Res* 21: 3-6.

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

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.

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.

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.

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

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

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

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.

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.

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.

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

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.

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.

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

59. Tarpy DR, Page RE (2002) Sex determination and the evolution of polyandry in honey bees (*Apis mellifera*). *Behav Ecol Sociobiol* 52: 143-150.

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

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.

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

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

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

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.

66. Laidlaw HHj, Page REj (1997) Queen rearing and bee breeding. *Cheshire CW*, ed.

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.

68. Grozinger CM, Fischer P, Hampton JE (2006) Uncoupling primer and releaser responses to pheromone in honey bees. *Naturwissenschaften* 94: 375-379.

69. Hoover SER, Oldroyd BP, Wossler TC, Winston ML (2005) Anarchistic queen honey bees have normal queen mandibular pheromones. *Insectes Soc* 52: 6-10.

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

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.

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

<figure><img src="image_13.png" /></figure>