Introduction {#s1}
============

Like many other species of social insect, the honey bee Apis mellifera L., is characterized by extreme reproductive division of labour. The vast majority of workers are functionally sterile, and never activate their ovaries throughout their life ([@i1536-2442-5-36-1-barron1]), leaving the production of eggs to the queen. Uncovering the genes associated with worker sterility in social insects will provide an important step forward in our understanding of a major evolutionary transition ([@i1536-2442-5-36-1-maynard1]), namely, the transition from individual organisms to integrated societies, and our best example yet of kin-selected 'genes for altruism' (sensu [@i1536-2442-5-36-1-dawkins1]): genes that reduce an individual\'s direct reproductive success while enhancing the reproductive success of relatives.

All worker honeybees develop into adults with vestigial ovaries. Although some or all wild type workers of intermediate age are able to subsequently activate their ovaries in response to a period of protracted queenlessness, and lay eggs, queenright workers almost never do so, retaining their ovaries in a vestigial state. An important exception is the selectively-bred 'anarchistic' line of honeybees maintained at The University of Sydney in which about 40% of 10 day-old workers have functionally activated ovaries and lay eggs at high frequency ([@i1536-2442-5-36-1-barron1]). This mutant line demonstrates that worker sterility has a genetic basis ([@i1536-2442-5-36-1-oldroyd1]). Furthermore, the patterns of inheritance of the anarchistic phenotype show that there are a small number of genes that regulate worker ovary activation in response to social cues ([@i1536-2442-5-36-1-oldroyd2]). Given the strong genetic component to the regulation of ovary activation in honeybees, the gene-expression profiles of workers with and without activated ovaries should reveal genes associated with the regulatory control of worker sterility.

For reasons that are not well understood, the ovaries of unmated (virgin) queens can be experimentally activated, as if mated, by subjecting queens to double CO~2~ narcosis ([@i1536-2442-5-36-1-mackensen1]). This phenomenon is routinely exploited to induce oviposition in unmated but artificially inseminated queens ([@i1536-2442-5-36-1-laidlaw1]). This idiosyncratic response has facilitated the physiological study of queen egg production ([@i1536-2442-5-36-1-engels1]; [@i1536-2442-5-36-1-engels1]; [@i1536-2442-5-36-1-engels3]). However, few studies have applied this technique to the study of worker egg production ([@i1536-2442-5-36-1-harris1]; [@i1536-2442-5-36-1-harris3]), and no study has yet examined locus-specific changes in gene-expression associated with ovary activation, in queens or workers. The effects of CO~2~ narcosis in queen honeybees indicate that the treatment triggers the up-regulation of genes that in turn activate the queen\'s ovaries. This is particularly evident in the case of vitellogenin, an egg protein, whose synthesis is dramatically stimulated in queen abdomens following CO~2~ treatment ([@i1536-2442-5-36-1-engels1]). Genes directly associated with ovary activation in queens, like vitellogenin for example, are prime candidates for genes associated with ovary activation in normally sterile workers.

We speculated that the genes associated with ovary activation in queens in response to CO~2~ narcosis might also be the genes associated with ovary activation in workers. To test this hypothesis we exposed groups of caged workers to double CO~2~ narcosis and compared rates of ovary activation and levels of gene expression against groups of control (un-narcotized) workers. Differential ovary activation and gene expression in response to CO~2~ exposure occurred. These findings are discussed in the context of the regulation of ovary activation in queens and workers.

Materials and Methods {#s2}
=====================

Experimental treatment {#s2a}
----------------------

The general approach was to use groups of caged queenless workers for comparing rates of ovary activation and levels of gene expression. To this end, we incubated emerging sealed brood from wild type colonies at 35° C overnight. The following morning, newly-emerged workers were collected and groups of 30 were transferred into wire-mesh cages (approx. 14 × 10 × 7 cm) fitted with a piece of beeswax foundation comb (6 cm × 3 cm), water and food (45% honey; 10% water; 45% royal jelly, Lifetime Health Products, Australia). After two (Experiment 1) or four (Experiment 2) days of incubation (∼35° C, 40% RH), workers were narcotized at room temperature by placing whole cages into a larger plastic container and flushing the container with compressed CO~2~ until the workers were immobilized. The workers were maintained in the narcotised state for 10 minutes. The following day, workers were narcotised again, for 3 minutes (after [@i1536-2442-5-36-1-engels1]). In each case, workers were completely anaesthetized. Then, at various time intervals following the double narcosis treatment (4-hrs, 24-hrs, 48-hrs, 96-hrs), workers were collected from incubation by snap-freezing them in liquid N~2~. Thus, collections were single pairs of cages sampled at each collection interval.

Bee abdomens were later dissected according to [@i1536-2442-5-36-1-dade1] and assigned a numerical score reflecting the state of ovary activation: 0 (ovaries thin and lacking defined ova), 1 (ovaries slightly thickened but still lacking defined ova), 2 (ovaries thick with clearly defined ova), 3 (ovaries thick with a least one fully mature ova). For each collection, differences in levels of ovary activation between treated and control cages were evaluated using t-tests for equality of mean ovary score between cages, as well as *X*^2^ tests of the proportions of each ovary activation class between treatments.

RNA extraction and cDNA synthesis {#s2b}
---------------------------------

Total RNA was extracted from the abdomens of eight workers per treated and untreated cage. A combined Trizol/Qiagen protocol was used to extract RNA under standard conditions so that it was suitable for quantifying levels of mRNA. Briefly, frozen tissues were homogenized in 1.5 ml Eppendorf tubes, containing 50--70 µl of Trizol reagent (Invitrogen, [www.invitrogen.com](www.invitrogen.com)), using disposable pestles from Sigma (Sigma-Aldrich, [www.sigmaaldrich.com](www.sigmaaldrich.com)) attached to a hand-held engraving device (Arlec, Super Tool, [www.arlec.com.au/](www.arlec.com.au/)). One abdomen was used per tube. Following the homogenization step the volume was adjusted to 500 µl with Trizol and then mixed with 100 µl of chloroform. The RNA-containing buffer phase was recovered after spinning for 10 min at 10,000 × g, mixed with an equal volume of 70% ethanol and applied on a Qiagen ([www.qiagen.com](www.qiagen.com)) RNeasy column. The remaining steps were carried out according to Qiagen\'s protocols. Standardised aliquots of individual RNA extractions were pooled by treatment prior to cDNA synthesis. We used pooled rather than individual samples to examine gene expression in order to minimize inter-individual variations, as described in other studies ([@i1536-2442-5-36-1-grozinger1]; [@i1536-2442-5-36-1-tian1]). Reverse transcription was performed according to Invitrogen\'s protocol using 5 µl of the pooled total RNA. The cDNA synthesis primer was T~(20)~MN (M=A,G, or C; N=A, G, C, or T).

Selection of candidate markers for ovarian activation in the honeybee {#s2c}
---------------------------------------------------------------------

The genes selected for transcriptional evaluation in this study encode transcripts known, or that were suspected, to be involved either in ovarian development (vitellogenin, transferrin, profilin, octopamine receptor) or in signalling pathways critical for cellular growth or cellular differentiation (flotillin, nitric oxide synthase, take-out-like). Arginine kinase was also selected because it is an important component of energy transfer and its differential expression in the honey bee has already been reported ([@i1536-2442-5-36-1-kucharski1]). Selection of these candidate markers was based on: 1) their functional gene descriptions at Interactive Fly ([www.sdbonline.org](www.sdbonline.org)) or Online Mendelian Inheritance in Man ([www.ncbi.nlm.nih.gov](www.ncbi.nlm.nih.gov)) and 2) the availability of an annotated gene sequence in the honey bee (<http://www.hgsc.bcm.tmc.edu/projects/honeybee/>). [Table 1](#i1536-2442-5-36-1-t01){ref-type="table"} provides more details regarding the biochemical and functional properties of these candidate genes.

Quantitative PCR {#s2d}
----------------

Target sequences for PCR were selected from cDNAs and ESTs using the genomic scaffolds at the Baylor College of Medicine. To eliminate genomic DNA amplification we manually designed primers that either spaned intron/exon junctions or amplified exon sequences separated by long introns. Our primers were 40-60% GC-rich and contained 20-24 bases to maximize reaction efficacy. All primers were experimentally tested in an Eppendorf gradient cycler to determine the optimal annealing temperature and to ensure that only one band of correct size was produced. These gene-specific primers are listed in [Table 2](#i1536-2442-5-36-1-t02){ref-type="table"}.

The relative amount of specific cDNA template between treated and control samples was quantified using real time quantitative PCR (qRT-PCR). *Taq* polymerase (0.2 units) from Promega ([www.promega.com](www.promega.com)) was used in 20 µl reactions (in triplicate) containing 2 mM MgCl~2~, 200 µM of each dNTP, 0.25 µM each of forward and reverse primer, and cDNA template equivalent to approximately 25 ng of total RNA. For every reaction, an RNA sample without reverse transcriptase was included to control for genomic DNA contamination. Product formation was monitored by the inclusion of SYBR Green I (Fisher Biotech, [www.fishersci.com](www.fishersci.com)) at a final dilution of 1:40,000. Thermocycling was conducted in a RotorGene 3000 Thermal Cycler (Corbett Research, [www.corbettresearch.com](www.corbettresearch.com)) for 35 cycles consisting of denaturation for 30 s at 94° C, annealing at 60° C for 30 s, extension at 72° C for 30 s, and fluorescence acquisition at 84° C for 15 s. Cycling was preceded by a 15 min 95° C activation step. Specificity of amplification was confirmed through a melt curve analysis of final PCR products by ramping the rotor temperature from 55° C to 99° C at 0.2° C ^s-1^ with fluorescence acquired after every 1° C increase.

Estimation of Changes in Transcript Abundance {#s2e}
---------------------------------------------

The amplification efficiencies of primer pairs used in qRT-PCR assays varied by less than 10% compared with that of the normalizer (ribosomal protein S8; [Table 1](#i1536-2442-5-36-1-t01){ref-type="table"}), which was ∼1.85. These values were determined by the slope of the curve generated by amplification of serially diluted cDNA over at least three orders of magnitude (data not shown). Accordingly, all primer pairs were nominated as having approximately equal efficiencies and, under this simplifying assumption, fold changes in relative transcript abundance were calculated using the method of [@i1536-2442-5-36-1-livak1].

Results {#s3}
=======

CO~2~ narcosis and ovary activation - Experiment 1 {#s3a}
--------------------------------------------------

Very young workers, bees narcotised on days 2 and 3 and examined just 4 hrs after treatment, had relatively low mean ovary scores, with no individual scored as having well-defined ova (i.e. all scores ≤ 1). Neither the mean score [(Figure 1A)](#i1536-2442-5-36-1-f01){ref-type="fig"} nor the categorical distribution of ovary activation scores (*X*^2^ ~1~= 0.069, *P* = 0.8) differed between treated and control cages. At 24 hrs post-treatment a general increase to mean score was observed in both treated and control groups. However, this increase did not differ between treatments *X*^2^ ~1~= 0.271, *P* = 0.6, [(Fig.1)](#i1536-2442-5-36-1-f01){ref-type="fig"}.

By 48 hrs post-treatment, however, an effect was apparent. Control workers showed significantly greater levels of ovary activation than did treated workers, *X*^2^ ~2~= 7.735, *P* = 0.021 [(Fig. 1A)](#i1536-2442-5-36-1-f01){ref-type="fig"}. Finally, bees examined 96 hrs post-treatment showed the greatest difference ovary activation between treated and control cages *X*^2^ ~3~= 10.764, *P* = 0.013 [(Fig. 1A)](#i1536-2442-5-36-1-f01){ref-type="fig"}, again with the control group showing significantly higher levels of ovary activation.

CO~2~ narcosis and ovary activation - Experiment 2 {#s3b}
--------------------------------------------------

For bees narcotised 4 and 5 days post emergence, the pattern was similar. At 4 hrs after treatment bees had low and similar ovary scores, *X*^2^ ~2~= 1.077, *P* = 0.6, while at 48 hrs following treatment controls showed significantly greater ovary activation than did treated workers *X*^2^ ~3~= 16.831, *P* = 0.001 [(Fig. 1B)](#i1536-2442-5-36-1-f01){ref-type="fig"}.

Quantitative real-time PCR {#s3c}
--------------------------

[Figure 2](#i1536-2442-5-36-1-f02){ref-type="fig"} shows the responsiveness of eight candidate genes to CO~2~ treatment. Bees narcotised on days 4 and 5 (Experiment 2) and collected 48 hrs later were used as the basis for examining gene expression differences. For all eight genes, the level of expression was calculated for treated bees relative to non-treated controls. Thus, genes whose relative expression was estimated at '1' were deemed to show no expression difference between treated and control groups, whereas genes significantly above or below this value, as evidenced by their 95% confidence intervals, were considered to be over- or under-expressed in treated workers, respectively. Genes encoding profilin (aka chickadee; PRF), flotillin (FLT), nitric oxide synthase (NOS), arginine kinase (ArgK), octopamine receptor (OCTR), and take-out-like carrier protein (JHBP-1) showed no consistent difference in expression [(Figure 2)](#i1536-2442-5-36-1-f02){ref-type="fig"}, even though differences in ovary activation were now apparent among the workers. Therefore these genes do not appear to be associated with experimentally induced differences in ovary activation, at least not by this assay. In contrast, the genes encoding vitellogenin (VIT) and transferrin (TRF) were down-regulated (∼ 4-5 fold; [Fig. 2](#i1536-2442-5-36-1-f02){ref-type="fig"}) in treated relative to control groups. The expression of these two genes is therefore associated with differences in functional ovary activation among workers.

Discussion {#s4}
==========

CO~2~-induced changes in ovary activation {#s4a}
-----------------------------------------

Our study confirms that double CO~2~ narcosis inhibits ovary activation in workers ([@i1536-2442-5-36-1-biedermann1]; [@i1536-2442-5-36-1-harris1]; [@i1536-2442-5-36-1-harris3]; [@i1536-2442-5-36-1-kropacova1]). This was true regardless of whether workers were narcotised when young (on days 2 and 3; [Fig. 1A](#i1536-2442-5-36-1-f01){ref-type="fig"}) or slightly older (on days 4 and 5; [Fig. 1B](#i1536-2442-5-36-1-f01){ref-type="fig"}), and all narcotised groups had mean ovary scores of less than '1'. In contrast, [@i1536-2442-5-36-1-mackensen1] showed that double CO~2~ narcosis of virgin queens actually accelerated ovary activation and egg laying. It is well known that in the absence of mating or narcosis, queens will not activate their ovaries or initiate oviposition for many weeks, if ever. Thus the contrasting effects of narcosis on ovary activation in queens versus workers provides a potentially powerful model system for studying the molecular mechanisms that regulate worker sterility in honeybees, and potentially other social insects.

Molecular correlates to ovary activation in worker honeybees {#s4b}
------------------------------------------------------------

This study shows that abdominal expression of vitellogenin and transferrin is reduced in narcotised workers compared to controls. This suggests that these two genes may be part of the network that regulates functional sterility in worker honeybees.

The process of ovary activation probably involves a hierarchy of events and hundreds of genes ([@i1536-2442-5-36-1-bownes1]), some of which may not be directly associated with the synthesis of ovarian tissue *per se*. For example, oocyte development depends on the synthetic activities of other, non-oocytic cells. For most insects, yolk is synthesised outside of the ovary and secondarily imported into the oocyte ([@i1536-2442-5-36-1-raikhel1]). Oocytes then take up the yolk proteins, such as vitellogenin, from the extracellular fluid by receptor-mediated endocytosis ([@i1536-2442-5-36-1-raikhel1]). Transferrin too can be incorporated via this mechanism into developing eggs, and likely functions to mediate iron-uptake ([@i1536-2442-5-36-1-kucharski3]; [@i1536-2442-5-36-1-kurama1]). By contrast, genes encoding proteins that are essential for the regulation of the actin cytoskeleton (i.e. flotillin) or linking the cytoskeleton with major signalling pathways (i.e. profilin) are expressed directly within the ovarian tissue and are tightly regulated during development ([@i1536-2442-5-36-1-cooley1]). By sampling whole abdomens in our assay genes involved in different functional categories related to egg production are potentially detected.

Of the eight candidate genes examined, six did not change their expression with changes in reproductive state. These genes represent different, but overlapping functional categories: receptor/sensing molecules (octopamine receptor, take-out-like carrier protein), signalling cascades (nitric oxide synthase, profilin, take-out-like), structural dynamics (profilin, flotillin), and energy transfer (arginine kinase). Broadly, these genes can be regarded as developmental regulators [(Table 1)](#i1536-2442-5-36-1-t01){ref-type="table"}, but either they are not differentially expressed between CO~2~-treated and control worker groups, or their differential expression was not detected by the assay.

Two genes, vitellogenin and transferrin, which are both involved in oocyte packaging, did show differences in expression (\> 4-fold; [Fig. 2](#i1536-2442-5-36-1-f02){ref-type="fig"}) following the double CO~2~ treatment. For vitellogenin, this response was temporally consistent with previous observations on like-treated virgin queens ([@i1536-2442-5-36-1-engels1]). Unlike queens, however, the change in worker vitellogenin expression was not conditional on bees being re-introduced into host colonies ([@i1536-2442-5-36-1-engels1]). Instead, workers responded to CO~2~ as they matured within cages, albeit as a group member rather than as lone individuals. More notably, the direction of the transcriptional response is opposite to that known for queens: workers decrease, rather than increase, vitellogenin synthesis in the abdomen following CO~2~ narcosis. This contrasting pattern of transcription in queens versus workers in response to CO~2~ treatment parallels the physiological pattern concerning ovary activation itself, and again accentuates the kin-selected difference in reproductive potential between these two castes.

Vitellogenin is a prerequisite for ovary activation in honeybees ([@i1536-2442-5-36-1-engels4]). Besides its proximate importance to egg production, however, vitellogenin is also used to synthesize proteinaceous royal jelly in the hypopharyngeal glands in the heads of workers. Nurse-age workers feed this jelly to dependent but related larvae, a behaviour that constitutes a form of kin-selected alloparental care. [@i1536-2442-5-36-1-amdam1] speculate that a key adaptation to eusociality by honeybees was the diversion of vitellogenin from its primary role as a yolk protein to a secondary role in the production of brood food. As a consequence, vitellogenin is expected to be functionally associated with colony-level traits such as alloparental care and reproductive division of labour ([@i1536-2442-5-36-1-amdam2]). Our finding that vitellogenin synthesis is down-regulated in abdomens of ovary-deactivated workers implies that this molecule is actively linked to the reproductive status of individual workers. In these experiments caged bees were not exposed to larvae and thus it is not certain whether they would show a concomitant up-regulation of vitellogenin in their heads associated with nursing and the production of royal jelly. Nontheless, these results, together with available knowledge from queens, suggests that there is a general link between CO~2~ exposure, vitellogenin synthesis, and caste-specific ovary activation in honeybees. Vitellogenin is a promising new candidate component in the regulatory pathway that controls functional sterility in workers, a pathway of great theoretical significance that has never been empirically deconstructed.

The relationship between transferrin and ovary activation is less clear because transferrin is truly multifunctional in the honey bee ([@i1536-2442-5-36-1-kucharski3]). However, its co-regulation with vitellogenin in the current study, and it\'s selective incorporation into eggs during oogenesis of Sarcophaga ([@i1536-2442-5-36-1-kurama1]) and Riptortus ([@i1536-2442-5-36-1-hirai1]), suggests that transferrin may have an important role in the activation of worker ovaries. Transferrin\'s likely function is to provide essential iron ions to developing oocytes and embryos ([@i1536-2442-5-36-1-hirai1]), and may also play a defensive role by sequestering iron away from pathogens that have entered the egg ([@i1536-2442-5-36-1-weinberg1]). Our finding that transferrin is down-regulated in ovary deactivated workers implies that, like vitellogenin, it is intimately associated with ovary activation in honey bee workers. It seems less likely that transferrin would have played an active role in the evolution of honey bee alloparental care or the reproductive division of labour, but nonetheless its known function and observed kinetics suggest it too could be a component in the regulatory pathway that controls functional sterility in workers.

Speculation on regulation of ovary activation {#s4c}
---------------------------------------------

The mechanism by which CO~2~ affects ovary activation via molecular intermediates in workers and queens is unknown, but there is some evidence to suggest that honeybees are sensitive to CO~2~ and that the gas is an important exogenous factor that modulates several aspects of social life. Firstly, honeybees are equipped with sensitive CO~2~ receptors on their antennae ([@i1536-2442-5-36-1-strange1]). Workers use these antennal receptors to tightly regulate CO~2~ concentrations within their colonies ([@i1536-2442-5-36-1-seeley1]). Second, CO~2~ narcosis is known to affect honey bee foraging ([@i1536-2442-5-36-1-ebadi1]), hoarding ([@i1536-2442-5-36-1-mardan1]), fanning behaviour ([@i1536-2442-5-36-1-seeley1]), sound production ([@i1536-2442-5-36-1-schneider1]), and some age-related polyethisms ([@i1536-2442-5-36-1-heran1]; [@i1536-2442-5-36-1-ribbands1]). Combined, these studies suggest that CO~2~ may be an important factor in the modulation of honeybee task specialization. Third, note that honeybee queens are kin-selected for extremely high fecundity and react to narcosis by accelerating reproductive development ([@i1536-2442-5-36-1-engels1]; [@i1536-2442-5-36-1-engels1]; [@i1536-2442-5-36-1-mackensen1]), whereas workers are selected for low fecundity or sterility and react to narcosis by retarding reproductive development ([@i1536-2442-5-36-1-biedermann1]: [@i1536-2442-5-36-1-harris1]; [@i1536-2442-5-36-1-harris3]; [@i1536-2442-5-36-1-kropacova1]; the present study). This contrasting effect of CO~2~ on queen versus worker reproduction suggests a caste-specific response to CO~2~ and, more generally, it suggests that CO~2~ differentially affects the regulatory mechanism underpinning honeybee reproductive division of labor.

One possibility is that exposure to CO~2~ affects the level of neurosecretions in worker brains ([@i1536-2442-5-36-1-harris3]). Levels of dopamine and serotonin, for example, co-vary with changes in worker ovarian development ([@i1536-2442-5-36-1-harris2]) and even with CO~2~-induced changes in ovarian development ([@i1536-2442-5-36-1-harris3]). It is possible therefore that CO~2~ triggers the regulatory mechanism that controls worker sterility by affecting, for example, the level of dopamine or serotonin, which are putative up-stream components in the regulatory pathway that controls ovary activation in workers. A second or dual possibility is that CO~2~ affects the titre of juvenile hormone (JH) in the hemolymph ([@i1536-2442-5-36-1-buhler1]), which in turn affects age-related behavior in workers, possibly including reproduction ([@i1536-2442-5-36-1-robinson1]). Note that, unlike for insects generally, JH has an inverse relationship with vitellogenin in honey bee workers: high JH titre turns off vitellogenin synthesis ([@i1536-2442-5-36-1-pinto1]). Thus, JH and vitellogenin are causally linked, and CO~2~ may stimulate the regulatory mechanism that normally controls worker sterility by causing an increase in JH, which causes a decrease in vitellogenin, which presumably retards ovary activation, as has observed been here.

Whatever the initial changes are that narcosis causes within the nervous system, and that ultimately gives rise to observed differences in ovary activation, we can begin to infer down-stream regulatory components via their differential expression on ovary activation or deactivation, in particular as a consequence of CO~2~ treatment. This approach will give a first indication of molecules important to the maintenance of worker sterility.

Concluding remarks {#s4d}
------------------

Gene expression involves a cascade of events, from transcription factor transactivation to RNA processing and maturation. The candidate genes approach ([@i1536-2442-5-36-1-fitzpatrick1]) used in this study is based on molecules that are probably involved in relatively down-stream events and therefore may not represent molecules that initially trigger the process of ovary activation or deactivation, in particular as a consequence of CO~2~ treatment. Nonetheless, by identifying molecular components of regulatory pathways via their differential expression, we can begin to describe the molecular circuitry of reproductive regulation. Subsequent comparative studies of honey bee castes and of other insects for which regulatory mechanisms are much better understood ([@i1536-2442-5-36-1-bownes1]) will then be possible. An association between rate of vitellogenin and transferrin transcription and worker ovary activation has not previously been demonstrated, and it is suggested that these two proteins are part of the network involved in the regulation of worker ovary functional activation, and thus involved in the regulation of functional sterility of workers.

The authors would like to thank Julianne Lim for technical assistance, and Andrew B Barron and Shelley ER Hoover for useful discussion. Support for this work was provided by a grant from the Australian Research Council to BPO and RM and by ARC Special Centre for the Molecular Genetics and Development (RM). PK was supported by grants from the Thailand Research Fund and the Royal Golden Jubilee Ph.D. Program for RGJ No. 4.C.CU/44/B.1

Figures and Tables
==================

![In *Experiment 1 Apis mellifera* workers were narcotized on days 2 and 3 of adult life, and their ovaries examined after 4 hrs, 24 hrs, 48 hrs, and 96 hrs. The mean ovary score of narcotized bees was significantly lower than controls at 48 hrs (t-test, *P* = 0.005) and 96 hrs (t-test, *P* = 0.003). **B)** In *Experiment 2* workers were narcotized on days 4 and 5 of adult life, and their ovaries examined after 4 hrs and 48 hrs. The mean ovary score of narcotized bees was significantly lower than controls at 48 hrs (t-test, *P* \< 0.001).](i1536-2442-5-36-1-f01){#i1536-2442-5-36-1-f01}

![The expression levels of eight candidate genes in Apis mellifera as estimated using real-time quantitative PCR. Relative quantities are normalized to an endogenous reference gene (S8, see text). Note that vitellogenin (VIT) and transferrin (TRF) are down-regulated following CO~2~-induced changes in reproductive state. Relative quantities represent an average from three assays of a single group of 8 individuals.](i1536-2442-5-36-1-f02){#i1536-2442-5-36-1-f02}

###### 

Description of genes used for transcriptional evaluation

![](i1536-2442-5-36-1-t01)

###### 

Sequences of gene-specific primers used for qRT-PCR assays

![](i1536-2442-5-36-1-t02)
