Characterization of quantitative trait loci for the age of first foraging in honey bee workers

Alternative double strand break repair pathways shape the evolution of high recombination in the honey bee, Apis mellifera

Social insects, particularly honey bees, have exceptionally high genomic frequencies of genetic recombination. This phenomenon and underlying mechanisms are poorly understood. To characterise the patterns of crossovers and gene conversion in the honey bee genome, a recombination map of 187 honey bee brothers was generated by whole-genome resequencing. Recombination events were heterogeneously distributed without many true hotspots. The tract lengths between phase shifts were bimodally distributed, indicating distinct crossover and gene conversion events. While crossovers predominantly occurred in G/C-rich regions and seemed to cause G/C enrichment, the gene conversions were found predominantly in A/T-rich regions. The nucleotide composition of sequences involved in gene conversions that were associated with or distant from crossovers corresponded to the differences between crossovers and gene conversions. These combined results suggest two types of DNA double-strand break repair during honey bee meiosis: non-canonical homologous recombination, leading to gene conversion and A/T enrichment of the genome, and the canonical homologous recombination based on completed double Holliday Junctions, which can result in gene conversion or crossover and is associated with G/C bias. This G/C bias may be selected for to balance the A/T-rich base composition of eusocial hymenopteran genomes. The lack of evidence for a preference of the canonical homologous recombination for double-strand break repair suggests that the high genomic recombination rate of honey bees is mainly the consequence of a high rate of double-strand breaks, which could in turn result from the life history of honey bees and their A/T-rich genome.

The Genomic Basis of Adaptation to High Elevations in Africanized Honey Bees

A range of different genetic architectures underpin local adaptation in nature. Honey bees (Apis mellifera) in the Eastern African Mountains harbor high frequencies of two chromosomal inversions that likely govern adaptation to this high-elevation habitat. In the Americas, honey bees are hybrids of European and African ancestries and adaptation to latitudinal variation in climate correlates with the proportion of these ancestries across the genome. It is unknown which, if either, of these forms of genetic variation governs adaptation in honey bees living at high elevations in the Americas. Here, we performed whole-genome sequencing of 29 honey bees from both high- and low-elevation populations in Colombia. Analysis of genetic ancestry indicated that both populations were predominantly of African ancestry, but the East African inversions were not detected. However, individuals in the higher elevation population had significantly higher proportions of European ancestry, likely reflecting local adaptation. Several genomic regions exhibited particularly high differentiation between highland and lowland bees, containing candidate loci for local adaptation. Genes that were highly differentiated between highland and lowland populations were enriched for functions related to reproduction and sperm competition. Furthermore, variation in levels of European ancestry across the genome was correlated between populations of honey bees in the highland population and populations at higher latitudes in South America. The results are consistent with the hypothesis that adaptation to both latitude and elevation in these hybrid honey bees are mediated by variation in ancestry at many loci across the genome.

Major royal jelly proteins influence the neurobiological regulation of the division of labor among honey bee workers

Age-based division of labor among workers is a fundamental life-history trait of many social insects, including the Western honey bee, Apis mellifera L. Extensive studies of the causation of the most pronounced transition from performing tasks in the nest to outside foraging indicate hormonal regulation of complex physiological changes. However, the proximate neurobiological mechanisms that cause the behavioral repertoire to change are still not understood and require novel approaches to be fully characterized. Thus, we established the first comprehensive monoclonal antibody microarray in honey bees with 16,320 antibodies to directly identify proteins in the brain that regulate the transition to foraging. Major royal jelly protein (MRJP) 1 and MRJP3 were identified as potential protein effectors and further investigated. A series of experimental manipulations of the workers' behavioral transition led to changes in MRJP1 and MRJP3 quantities in accordance with their presumed functional role. Injection of MRJPs into the brain resulted in increased task-reversal from foraging to nursing and decreased task-progression from nursing to foraging, while the latter was increased by injection with MRJP antibodies. Finally, down-regulation of MRJP1 and MRJP3 expression via RNAi injection into the brain increased the transition from in-hive nursing to outside foraging, confirming a causal role of these two proteins in the proximate regulation of behavior and life-history of honey bee workers. Interaction partners of MRJP1 and MRJP3 in the honey bee brain included other regulators of honey bee behavior and life history. Thus, our transformative methodological advancement of proteome analysis in honey bees reveals novel regulators of honey bee behavior, extends our understanding of the functional pleiotropy of MRJPs, and supports a general nutrition-based model of the regulation of the age-based division of labor in honey bees.

Societies to genes: can we get there from here?

Understanding the organization and evolution of social complexity is a major task because it requires building an understanding of mechanisms operating at different levels of biological organization from genes to social interactions. I discuss here, a unique forward genetic approach spanning more than 30 years beginning with human-assisted colony-level selection for a single social trait, the amount of pollen honey bees (Apis mellifera L.) store. The goal was to understand a complex social trait from the social phenotype to genes responsible for observed trait variation. The approach combined the results of colony-level selection with detailed studies of individual behavior and physiology resulting in a mapped, integrated phenotypic architecture composed of correlative relationships between traits spanning anatomy, physiology, sensory response systems, and individual behavior that affect individual foraging decisions. Colony-level selection reverse engineered the architecture of an integrated phenotype of individuals resulting in changes in the social trait. Quantitative trait locus (QTL) studies combined with an exceptionally high recombination rate (60 kb/cM), and a phenotypic map, provided a genotype-phenotype map of high complexity demonstrating broad QTL pleiotropy, epistasis, and epistatic pleiotropy suggesting that gene pleiotropy or tight linkage of genes within QTL integrated the phenotype. Gene expression and knockdown of identified positional candidates revealed genes affecting foraging behavior and confirmed one pleiotropic gene, a tyramine receptor, as a target for colony-level selection that was under selection in two different tissues in two different life stages. The approach presented here has resulted in a comprehensive understanding of the structure and evolution of honey bee social organization.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Tyramine and its receptor TYR1 linked behavior QTL to reproductive physiology in honey bee workers (Apis mellifera)

Honey bees (Apis mellifera) provide an excellent model for studying how complex social behavior evolves and is regulated. Social behavioral traits such as the division of labor have been mapped to specific genomic regions in quantitative trait locus (QTL) studies. However, relating genomic mapping to gene function and regulatory mechanism remains a big challenge for geneticists. In honey bee workers, division of labor is known to be regulated by reproductive physiology, but the genetic basis of this regulation remains unknown. In this case, QTL studies have identified tyramine receptor 1 (TYR1) as a candidate gene in region pln2, which is associated with multiple worker social traits and reproductive anatomy. Tyramine (TA), a neurotransmitter, regulates physiology and behavior in diverse insect species including honey bees. Here, we examine directly the effects of TYR1 and TA on worker reproductive physiology, including ovariole number, ovary function and the production of vitellogenin (VG, an egg yolk precursor). First, we used a pharmacology approach to demonstrate that TA affects ovariole number during worker larval development and increases ovary maturation during the adult stage. Second, we used a gene knockdown approach to show that TYR1 regulates vg transcription in adult workers. Finally, we estimated correlations in gene expression and propose that TYR1 may regulate vg transcription by coordinating hormonal and nutritional signals. Taken together, our results suggest TYR1 and TA play important roles in regulating worker reproductive physiology, which in turn regulates social behavior. Our study exemplifies a successful forward-genetic strategy going from QTL mapping to gene function.

Insulin Receptor Substrate Gene Knockdown Accelerates Behavioural Maturation and Shortens Lifespan in Honeybee Workers

In animals, dietary restriction or suppression of genes involved in nutrient sensing tends to increase lifespan. In contrast, food restriction in honeybees (Apis mellifera) shortens lifespan by accelerating a behavioural maturation program that culminates in leaving the nest as a forager. Foraging is metabolically demanding and risky, and foragers experience increased rates of aging and mortality. Food-deprived worker bees forage at younger ages and are expected to live shorter lives. We tested whether suppression of a molecular nutrient sensing pathway is sufficient to accelerate the behavioural transition to foraging and shorten worker life. To achieve this, we reduced expression of the insulin receptor substrate (irs) gene via RNA interference in two selected lines of honeybees used to control for behavioural and genetic variation. irs encodes a membrane-associated protein in the insulin/insulin-like signalling (IIS) pathway that is central to nutrient sensing in animals. We measured foraging onset and lifespan and found that suppression of irs reduced worker bee lifespan in both genotypes, and that this effect was largely driven by an earlier onset of foraging behaviour in a genotype-conditional manner. Our results provide the first direct evidence that an IIS pathway gene influences behavioural maturation and lifespan in honeybees and highlight the importance of considering social environments and behaviours when investigating the regulation of aging and lifespan in social animals.

Adaptive maintenance of European alleles in the Brazilian Africanized honeybee

The Anthropocene is an epoch hallmarked by intensified human intrusion across ecosystems. One such intrusion is the movement and re-introduction of long-separated populations. By facilitating introgression - intraspecific genetic admixture - secondary contact can facilitate range expansion and the establishment of invasive species. The proximate mechanisms through which introgression facilitates expansion are rarely known (Bock et al., ; Rius & Darling, ), but managed species provide a useful avenue for exploration. Bee-keepers have been interbreeding highly diverged honeybee clades for centuries, often to introduce "useful" phenotypic variation to their stocks. Across the Western honeybee's (Apis mellifera) European range, this practice has not resulted in range expansion (Moritz, Härtel, & Neumann, ). In the Americas, however, introgression of European with African subspecies resulted in a widely publicized invasive population: The Africanized honeybee (AHB). In this issue of Molecular Ecology, Nelson, Wallberg, Simões, Lawson, and Webster () have made the first step towards understanding how this invasive species successfully spread across the Americas.

Recent advances in population and quantitative genomics of honey bees

The increase in the availability of individual Apis mellifera genomes has resulted in significant progress toward understanding the evolution and adaptation of the honey bee. These efforts have identified new subspecies, evolutionary lineages, and a significant number of genes involved with adaptations and colony-level quantitative traits. Many studies have also developed genetic assays that are being used to monitor the movement and admixture of honey bee populations. These resources are valuable for conservation and breeding programs that seek to improve the economic value of colonies or preserve locally adapted populations and subspecies. This review provides a brief discussion on how population and quantitative genomic studies has improved our understanding of the honey bee.

Expression Characterization and Localization of the foraging Gene in the Chinese Bee, Apis cerana cerana (Hymenoptera: Apidae)

In social insects, the foraging gene (for) regulates insect age- and task-based foraging behaviors. We studied the expression and localization of the for gene (Acfor) in Apis cerana cerana workers to explore whether the differential regulation of this gene is associated with the behaviors of nurses and foragers. The expression profiles of Acfor in different tissues and at different ages were examined using real-time quantitative reverse transcription polymerase chain reaction. Cellular localization in the brain was detected using in situ hybridization. Acfor transcripts in different ages workers showed that Acfor expression was detected in all the heads of 1- to 30-d-old worker bees. Acfor expression reached a peak at 25 d of age, and then declined with increasing age. The results showed that Acfor gene expression in five tissues was respectively significantly higher in foragers than in nurses. In nurses, the relative expression of Acfor was the highest in the antennae. There was a highly significant difference in expression between antennae, legs, and the other three tissues. In foragers, Acfor expression was the highest in the thorax, which was significantly different from all other tissues. In situ hybridization showed that Acfor was highly expressed in the lamina of the optic lobes, in a central column of Kenyon cells in the mushroom bodies of the brain of workers, and in the antennal lobes. This suggested that Acfor expression affects age-related foraging behavior in Apis cerana cerana, and that it may be related to flight activity.

Genomewide analysis of admixture and adaptation in the Africanized honeybee

Genetic exchange by hybridization or admixture can make an important contribution to evolution, and introgression of favourable alleles can facilitate adaptation to new environments. A small number of honeybees (Apis mellifera) with African ancestry were introduced to Brazil ~60 years ago, which dispersed and hybridized with existing managed populations of European origin, quickly spreading across much of the Americas in an example of a massive biological invasion. Here, we analyse whole-genome sequences of 32 Africanized honeybees sampled from throughout Brazil to study the effect of this process on genome diversity. By comparison with ancestral populations from Europe and Africa, we infer that these samples have 84% African ancestry, with the remainder from western European populations. However, this proportion varies across the genome and we identify signals of positive selection in regions with high European ancestry proportions. These observations are largely driven by one large gene-rich 1.4-Mbp segment on chromosome 11 where European haplotypes are present at a significantly elevated frequency and likely confer an adaptive advantage in the Africanized honeybee population. This region has previously been implicated in reproductive traits and foraging behaviour in worker bees. Finally, by analysing the distribution of ancestry tract lengths in the context of the known time of the admixture event, we are able to infer an average generation time of 2.0 years. Our analysis highlights the processes by which populations of mixed genetic ancestry form and adapt to new environments.

Genetic architecture of a hormonal response to gene knockdown in honey bees

Variation in endocrine signaling is proposed to underlie the evolution and regulation of social life histories, but the genetic architecture of endocrine signaling is still poorly understood. An excellent example of a hormonally influenced set of social traits is found in the honey bee (Apis mellifera): a dynamic and mutually suppressive relationship between juvenile hormone (JH) and the yolk precursor protein vitellogenin (Vg) regulates behavioral maturation and foraging of workers. Several other traits cosegregate with these behavioral phenotypes, comprising the pollen hoarding syndrome (PHS) one of the best-described animal behavioral syndromes. Genotype differences in responsiveness of JH to Vg are a potential mechanistic basis for the PHS. Here, we reduced Vg expression via RNA interference in progeny from a backcross between 2 selected lines of honey bees that differ in JH responsiveness to Vg reduction and measured JH response and ovary size, which represents another key aspect of the PHS. Genetic mapping based on restriction site-associated DNA tag sequencing identified suggestive quantitative trait loci (QTL) for ovary size and JH responsiveness. We confirmed genetic effects on both traits near many QTL that had been identified previously for their effect on various PHS traits. Thus, our results support a role for endocrine control of complex traits at a genetic level. Furthermore, this first example of a genetic map of a hormonal response to gene knockdown in a social insect helps to refine the genetic understanding of complex behaviors and the physiology that may underlie behavioral control in general.

The Architecture of the Pollen Hoarding Syndrome in Honey Bees: Implications for Understanding Social Evolution, Behavioral Syndromes, and Selective Breeding

Social evolution has influenced every aspect of contemporary honey bee biology, but the details are difficult to reconstruct. The reproductive ground plan hypothesis of social evolution proposes that central regulators of the gonotropic cycle of solitary insects have been coopted to coordinate social complexity in honey bees, such as the division of labor among workers. The predicted trait associations between reproductive physiology and social behavior have been identified in the context of the pollen hoarding syndrome, a larger suite of interrelated traits. The genetic architecture of this syndrome is characterized by a partially overlapping genetic architecture with several consistent, pleiotropic QTL. Despite these central QTL and an integrated hormonal regulation, separate aspects of the pollen hoarding syndrome may evolve independently due to peripheral QTL and additionally segregating genetic variance. The characterization of the pollen hoarding syndrome has also demonstrated that this syndrome involves many non-behavioral traits, which may be the case for numerous "behavioral" syndromes. Furthermore, the genetic architecture of the pollen hoarding syndrome has implications for breeding programs for improving honey health and other desirable traits: If these traits are comparable to the pollen hoarding syndrome, consistent pleiotropic QTL will enable marker assisted selection, while sufficient additional genetic variation may permit the dissociation of trade-offs for efficient multiple trait selection.

Multiple mating but not recombination causes quantitative increase in offspring genetic diversity for varying genetic architectures

Explaining the evolution of sex and recombination is particularly intriguing for some species of eusocial insects because they display exceptionally high mating frequencies and genomic recombination rates. Explanations for both phenomena are based on the notion that both increase colony genetic diversity, with demonstrated benefits for colony disease resistance and division of labor. However, the relative contributions of mating number and recombination rate to colony genetic diversity have never been simultaneously assessed. Our study simulates colonies, assuming different mating numbers, recombination rates, and genetic architectures, to assess their worker genotypic diversity. The number of loci has a strong negative effect on genotypic diversity when the allelic effects are inversely scaled to locus number. In contrast, dominance, epistasis, lethal effects, or limiting the allelic diversity at each locus does not significantly affect the model outcomes. Mating number increases colony genotypic variance and lowers variation among colonies with quickly diminishing returns. Genomic recombination rate does not affect intra- and inter-colonial genotypic variance, regardless of mating frequency and genetic architecture. Recombination slightly increases the genotypic range of colonies and more strongly the number of workers with unique allele combinations across all loci. Overall, our study contradicts the argument that the exceptionally high recombination rates cause a quantitative increase in offspring genotypic diversity across one generation. Alternative explanations for the evolution of high recombination rates in social insects are therefore needed. Short-term benefits are central to most explanations of the evolution of multiple mating and high recombination rates in social insects but our results also apply to other species.

Evolution of multiple additive loci caused divergence between Drosophila yakuba and D. santomea in wing rowing during male courtship

In Drosophila, male flies perform innate, stereotyped courtship behavior. This innate behavior evolves rapidly between fly species, and is likely to have contributed to reproductive isolation and species divergence. We currently understand little about the neurobiological and genetic mechanisms that contributed to the evolution of courtship behavior. Here we describe a novel behavioral difference between the two closely related species D. yakuba and D. santomea: the frequency of wing rowing during courtship. During courtship, D. santomea males repeatedly rotate their wing blades to face forward and then back (rowing), while D. yakuba males rarely row their wings. We found little intraspecific variation in the frequency of wing rowing for both species. We exploited multiplexed shotgun genotyping (MSG) to genotype two backcross populations with a single lane of Illumina sequencing. We performed quantitative trait locus (QTL) mapping using the ancestry information estimated by MSG and found that the species difference in wing rowing mapped to four or five genetically separable regions. We found no evidence that these loci display epistasis. The identified loci all act in the same direction and can account for most of the species difference.

Genetics of reproduction and regulation of honeybee (Apis mellifera L.) social behavior

Honeybees form complex societies with a division of labor for reproduction, nutrition, nest construction and maintenance, and defense. How does it evolve? Tasks performed by worker honeybees are distributed in time and space. There is no central control over behavior and there is no central genome on which selection can act and effect adaptive change. For 22 years, we have been addressing these questions by selecting on a single social trait associated with nutrition: the amount of surplus pollen (a source of protein) that is stored in the combs of the nest. Forty-two generations of selection have revealed changes at biological levels extending from the society down to the level of the gene. We show how we constructed this vertical understanding of social evolution using behavioral and anatomical analyses, physiology, genetic mapping, and gene knockdowns. We map out the phenotypic and genetic architectures of food storage and foraging behavior and show how they are linked through broad epistasis and pleiotropy affecting a reproductive regulatory network that influences foraging behavior. This is remarkable because worker honeybees have reduced reproductive organs and are normally sterile; however, the reproductive regulatory network has been co-opted for behavioral division of labor.

Complex pleiotropy characterizes the pollen hoarding syndrome in honey bees (Apis mellifera L.)

The pollen hoarding syndrome consists of a large suite of correlated traits in honey bees that may have played an important role in colony organization and consequently the social evolution of honey bees. The syndrome was first discovered in two strains that have been artificially selected for high and low pollen hoarding. These selected strains are used here to further investigate the phenotypic and genetic links between two central aspects of the pollen hoarding syndrome, sucrose responsiveness and pollen hoarding. Sons of hybrid queen offspring of these two strains were tested for sucrose responsiveness and used to produce colonies with either a highly responsive or an unresponsive father. These two colony groups differed significantly in the amount of pollen stored on brood combs and with regards to their relationship between brood and pollen amounts. Additionally, four quantitative trait loci (QTL) for pollen hoarding behavior were assessed for their effect on sucrose responsiveness. Drone offspring of two hybrid queens were phenotyped for responsiveness and genotyped at marker loci for these QTL, identifying some pleiotropic effects of the QTL with significant QTL interactions. Both experiments thus provided corroborating evidence that the distinct traits of the pollen hoarding syndrome are mechanistically and genetically linked, and that these links are complex and dependent on background genotype. The study demonstrates genetic worker-drone correlations within the context of the pollen hoarding syndrome and establishes that an indirect selection response connects pollen hoarding and sucrose responsiveness, regardless of which trait is directly selected.

Pleiotropy of segregating genetic variants that affect honey bee worker life expectancy

In contrast to many other complex traits, the natural genetic architecture of life expectancy has not been intensely studied, particularly in non-model organisms, such as the honey bee (Apis mellifera L.). Multiple factors that determine honey bee worker lifespan have been identified and genetic analyses have been performed on some of those traits. Several of the traits are included in a suite of correlated traits that form the pollen hoarding syndrome, which was named after the behavior to store surplus pollen in the nest and is tied to social evolution. Here, seven quantitative trait loci that had previously been identified for their effects on different aspects of the pollen hoarding syndrome were studied for their genetic influence on the survival of adult honey bee workers. To gain a more comprehensive understanding of the genetic architecture of worker longevity, a panel of 280 additional SNP markers distributed across the genome was also tested. Allelic distributions were compared between young and old bees in two backcross populations of the bi-directionally selected high- and low-pollen hoarding strain. Our results suggest a pleiotropic effect of at least one of the behavioral quantitative trait loci on worker longevity and one significant and several other putative genetic effects in other genomic regions. At least one locus showed evidence for strong antagonistic pleiotropy and several others suggested genetic factors that influence pre-emergence survival of worker honey bees. Thus, the predicted association between worker lifespan and the pollen hoarding syndrome was supported at the genetic level and the magnitude of the identified effects also strengthened the view that naturally segregating genetic variation can have major effects on age-specific survival probability in the wild.

Support for the reproductive ground plan hypothesis of social evolution and major QTL for ovary traits of Africanized worker honey bees (Apis mellifera L.)

BACKGROUND

The reproductive ground plan hypothesis of social evolution suggests that reproductive controls of a solitary ancestor have been co-opted during social evolution, facilitating the division of labor among social insect workers. Despite substantial empirical support, the generality of this hypothesis is not universally accepted. Thus, we investigated the prediction of particular genes with pleiotropic effects on ovarian traits and social behavior in worker honey bees as a stringent test of the reproductive ground plan hypothesis. We complemented these tests with a comprehensive genome scan for additional quantitative trait loci (QTL) to gain a better understanding of the genetic architecture of the ovary size of honey bee workers, a morphological trait that is significant for understanding social insect caste evolution and general insect biology.

RESULTS

Back-crossing hybrid European x Africanized honey bee queens to the Africanized parent colony generated two study populations with extraordinarily large worker ovaries. Despite the transgressive ovary phenotypes, several previously mapped QTL for social foraging behavior demonstrated ovary size effects, confirming the prediction of pleiotropic genetic effects on reproductive traits and social behavior. One major QTL for ovary size was detected in each backcross, along with several smaller effects and two QTL for ovary asymmetry. One of the main ovary size QTL coincided with a major QTL for ovary activation, explaining 3/4 of the phenotypic variance, although no simple positive correlation between ovary size and activation was observed.

CONCLUSIONS

Our results provide strong support for the reproductive ground plan hypothesis of evolution in study populations that are independent of the genetic stocks that originally led to the formulation of this hypothesis. As predicted, worker ovary size is genetically linked to multiple correlated traits of the complex division of labor in worker honey bees, known as the pollen hoarding syndrome. The genetic architecture of worker ovary size presumably consists of a combination of trait-specific loci and general regulators that affect the whole behavioral syndrome and may even play a role in caste determination. Several promising candidate genes in the QTL intervals await further study to clarify their potential role in social insect evolution and the regulation of insect fertility in general.

Genetic architecture of ovary size and asymmetry in European honeybee workers

The molecular basis of complex traits is increasingly understood but a remaining challenge is to identify their co-regulation and inter-dependence. Pollen hoarding (pln) in honeybees is a complex trait associated with a well-characterized suite of linked behavioral and physiological traits. In European honeybee stocks bidirectionally selected for pln, worker (sterile helper) ovary size is pleiotropically affected by quantitative trait loci that were initially identified for their effect on foraging behavior. To gain a better understanding of the genetic architecture of worker ovary size in this model system, we analyzed a series of crosses between the selected strains. The crossing results were heterogeneous and suggested non-additive effects. Three significant and three suggestive quantitative trait loci of relatively large effect sizes were found in two reciprocal backcrosses. These loci are not located in genome regions of known effects on foraging behavior but contain several interesting candidate genes that may specifically affect worker-ovary size. Thus, the genetic architecture of this life history syndrome may be comprised of pleiotropic, central regulators that influence several linked traits and other genetic factors that may be downstream and trait specific.

Genotype effect on regulation of behaviour by vitellogenin supports reproductive origin of honeybee foraging bias

In honeybee colonies, food collection is performed by a group of mostly sterile females called workers. After an initial nest phase, workers begin foraging for nectar and pollen, but tend to bias their collection towards one or the other. The foraging choice of honeybees is influenced by vitellogenin (vg), an egg-yolk precursor protein that is expressed although workers typically do not lay eggs. The forager reproductive ground plan hypothesis (RGPH) proposes an evolutionary path in which the behavioural bias toward collecting nectar or pollen on foraging trips is influenced by variation in reproductive physiology, such as hormone levels and vg gene expression. Recently, the connections between vg and foraging behaviour were challenged by Oldroyd and Beekman (2008), who concluded from their study that the ovary, and especially vg, played no role in foraging behaviour of bees. We address their challenge directly by manipulating vg expression by RNA interference- (RNAi) mediated gene knockdown in two honeybee genotypes with different foraging behaviour and reproductive physiology. We show that the effect of vg on the food-loading decisions of the workers occurs only in the genotype where timing of foraging onset (by age) is also sensitive to vg levels. In the second genotype, changing vg levels do not affect foraging onset or bias. The effect of vg on workers' age at foraging onset is explained by the well-supported double repressor hypothesis (DHR), which describes a mutually inhibitory relationship between vg and juvenile hormone (JH) - an endocrine factor that influences development, reproduction, and behaviour in many insects. These results support the RGPH and demonstrate how it intersects with an established mechanism of honeybee behavioural control.

The genetic basis of transgressive ovary size in honeybee workers

Ovarioles are the functional unit of the female insect reproductive organs and the number of ovarioles per ovary strongly influences egg-laying rate and fecundity. Social evolution in the honeybee (Apis mellifera) has resulted in queens with 200-360 total ovarioles and workers with usually 20 or less. In addition, variation in ovariole number among workers relates to worker sensory tuning, foraging behavior, and the ability to lay unfertilized male-destined eggs. To study the genetic architecture of worker ovariole number, we performed a series of crosses between Africanized and European bees that differ in worker ovariole number. Unexpectedly, these crosses produced transgressive worker phenotypes with extreme ovariole numbers that were sensitive to the social environment. We used a new selective pooled DNA interval mapping approach with two Africanized backcrosses to identify quantitative trait loci (QTL) underlying the transgressive ovary phenotype. We identified one QTL on chromosome 11 and found some evidence for another QTL on chromosome 2. Both QTL regions contain plausible functional candidate genes. The ovariole number of foragers was correlated with the sugar concentration of collected nectar, supporting previous studies showing a link between worker physiology and foraging behavior. We discuss how the phenotype of extreme worker ovariole numbers and the underlying genetic factors we identified could be linked to the development of queen traits.