Search

Little is known about the combined effects of stressors on social immunity of honey bees (Apis mellifera) and related gene expression. The interaction between sublethal doses of a neurotoxin, clothianidin, and the ectoparasite, Varroa destructor, was examined by measuring differentially expressed genes (DEGs) in brains, deformed wing virus (DWV) and the proportion and intensity of self-grooming. Evidence for an interaction was observed between the stressors in a reduction in the proportion of intense groomers. Only the lowest dose of clothianidin alone reduced the proportion of self-groomers and increased DWV levels. V. destructor shared a higher proportion of DEGs with the combined stressors compared to clothianidin, indicating that the effects of V. destructor were more pervasive than those of clothianidin when they were combined. The number of up-regulated DEGs were reduced with the combined stressors compared to clothianidin alone, suggesting an interference with the impacts of clothianidin. Clothianidin and V. destructor affected DEGs from different biological pathways but shared impacts on pathways related to neurodegenerative disorders, like Alzheimer’s, which could be related to neurological dysfunction and may explain their negative impacts on grooming. This study shows that the combination of clothianidin and V. destructor resulted in a complex and non-additive interaction.

Sexual reproduction brings genes from two parents (matrigenes and patrigenes) together into one individual. These genes, despite being unrelated, should show nearly perfect cooperation because each gains equally through the production of offspring. However, an individual's matrigenes and patrigenes can have different probabilities of being present in other relatives, so kin selection could act on them differently. Such intragenomic conflict could be implemented by partial or complete silencing (imprinting) of an allele by one of the parents. Evidence supporting this theory is seen in offspring-mother interactions, with patrigenes favoring acquisition of more of the mother's resources if some of the costs fall on half-siblings who do not share the patrigene. The kinship theory of intragenomic conflict is little tested in other contexts, but it predicts that matrigene-patrigene conflict may be rife in social insects. We tested the hypothesis that honey bee worker reproduction is promoted more by patrigenes than matrigenes by comparing across nine reciprocal crosses of two distinct genetic stocks. As predicted, hybrid workers show reproductive trait characteristics of their paternal stock, (indicating enhanced activity of the patrigenes on these traits), greater patrigenic than matrigenic expression, and significantly increased patrigenic-biased expression in reproductive workers. These results support both the general prediction that matrigene-patrigene conflict occurs in social insects and the specific prediction that honey bee worker reproduction is driven more by patrigenes. The success of these predictions suggests that intragenomic conflict may occur in many contexts where matrigenes and patrigenes have different relatednesses to affected kin.

Parent-specific gene expression (PSGE) is little known outside of mammals and plants. PSGE occurs when the expression level of a gene depends on whether an allele was inherited from the mother or the father. Kin selection theory predicts that there should be extensive PSGE in social insects because social insect parents can gain inclusive fitness benefits by silencing parental alleles in female offspring. We searched for evidence of PSGE in honey bees using transcriptomes from reciprocal crosses between European and Africanized strains. We found 46 transcripts with significant parent-of-origin effects on gene expression, many of which overexpressed the maternal allele. Interestingly, we also found a large proportion of genes showing a bias toward maternal alleles in only one of the reciprocal crosses. These results indicate that PSGE may occur in social insects. The nonreciprocal effects could be largely driven by hybrid incompatibility between these strains. Future work will help to determine if these are indeed parent-of-origin effects that can modulate inclusive fitness benefits.

In order to identify genes that are influencing defensive behaviors, we have taken a new approach by dissecting colony-level defensive behavior into individual behavioral measurements using two families containing backcross workers from matings involving European and Africanized bees. We removed the social context from stinging behavior by using a laboratory assay to measure the stinging response of individual bees. A mild shock was given to bees using a constant-current stimulator. The time it took bees to sting in response to this stimulus was recorded. In addition, bees that were seen performing guard behaviors at the hive entrance were collected. We performed QTL mapping in two backcross families with SNP probes within genes and identified two new QTL regions for stinging behavior and another QTL region for guarding behavior. We also identified several candidate genes involved in neural signaling, neural development and muscle development that may be influencing stinging and guarding behaviors. The lack of overlap between these regions and previous defensive behavior QTL underscores the complexity of this behavior and increases our understanding of its genetic architecture.

Invertebrate models have greatly furthered our understanding of ethanol sensitivity and alcohol addiction. The honey bee (Apis mellifera), a widely used behavioral model, is valuable for comparative studies. A quantitative trait locus (QTL) mapping experiment was designed to identify QTL and genes influencing ethanol vapor sensitivity. A backcross mating between ethanol-sensitive and resistant lines resulted in worker offspring that were tested for sensitivity to the sedative effects of alcohol. A linkage map was constructed with over 500 amplified fragment length polymorphism (AFLP) and sequence-tagged site (STS) markers. Four QTL were identified from three linkage groups with log of odds ratio (LOD) scores of 2.28, 2.26, 2.23, and 2.02. DNA from markers within and near QTL were cloned and sequenced, and this data was utilized to integrate our map with the physical honey bee genome. Many candidate genes were identified that influence synaptic transmission, neuronal growth, and detoxification. Others affect lipid synthesis, apoptosis, alcohol metabolism, cAMP signaling, and electron transport. These results are relevant because they present the first search for QTL that affect resistance to acute ethanol exposure in an invertebrate, could be useful for comparative genomic purposes, and lend credence to the use of honey bees as biomedical models of alcohol metabolism and sensitivity.

Several candidate genes identified from quantitative trait loci (QTL) for defensive behavior in honey bees (Apis mellifera L.) are homologous to genes known to influence ethanol sensitivity in other organisms. To investigate this possible link between aggression/defense and ethanol sensitivity, assays were developed to evaluate ethanol vapor responses in worker bees from a low-defensive (gentle) colony and a high-defensive colony. Defensive workers exhibited characteristic signs of ethanol-induced sedation significantly faster than gentle workers upon exposure to ethanol vapor. Backcross workers displayed ethanol sensitivity intermediate to the parental defensive and gentle lines, suggesting a genetic basis for the trait. Workers were screened with sequence-tagged site markers linked to three defensive-behavior QTL and their genotypes were tested for associations with ethanol sensitivity. There were no significant associations, indicating that the defensive QTL were not having a pleiotropic effect on ethanol sensitivity. It is possible that gentle-source alleles at these QTL are dominant with respect to sensitivity, one or more of these QTL were not segregating in the backcross family, or unidentified QTL are influencing alcohol sensitivity.

Background Varroa mites are widely considered the biggest honey bee health problem worldwide. Until recently, Varroa jacobsoni has been found to live and reproduce only in Asian honey bee (Apis cerana) colonies, while V. destructor successfully reproduces in both A. cerana and A. mellifera colonies. However, we have identified an island population of V. jacobsoni that is highly destructive to A. mellifera, the primary species used for pollination and honey production. The ability of these populations of mites to cross the host species boundary potentially represents an enormous threat to apiculture, and is presumably due to genetic variation that exists among populations of V. jacobsoni that influences gene expression and reproductive status. In this work, we investigate differences in gene expression between populations of V. jacobsoni reproducing on A. cerana and those either reproducing or not capable of reproducing on A. mellifera, in order to gain insight into differences that allow V. jacobsoni to overcome its normal species tropism. Results We sequenced and assembled a de novo transcriptome of V. jacobsoni. We also performed a differential gene expression analysis contrasting biological replicates of V. jacobsoni populations that differ in their ability to reproduce on A. mellifera. Using the edgeR, EBSeq and DESeq R packages for differential gene expression analysis, we found 287 differentially expressed genes (FDR ≤ 0.05), of which 91% were up regulated in mites reproducing on A. mellifera. In addition, mites found reproducing on A. mellifera showed substantially more variation in expression among replicates. We searched for orthologous genes in public databases and were able to associate 100 of these 287 differentially expressed genes with a functional description. Conclusions There is differential gene expression between the two mite groups, with more variation in gene expression among mites that were able to reproduce on A. mellifera. A small set of genes showed reduced expression in mites on the A. mellifera host, including putative transcription factors and digestive tract developmental genes. The vast majority of differentially expressed genes were up-regulated in this host. This gene set showed enrichment for genes associated with mitochondrial respiratory function and apoptosis, suggesting that mites on this host may be experiencing higher stress, and may be less optimally adapted to parasitize it. Some genes involved in reproduction and oogenesis were also overexpressed, which should be further studied in regards to this host shift. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3130-3) contains supplementary material, which is available to authorized users.

This study was conducted to analyze the stinging response thresholds of individual European and Africanized worker honeybees (Apis mellifera L.) to electrical stimulation. Newly emerged workers were identified, and either were placed into an incubator, into their natal colonies, or cross-fostered in common colonies of European or Africanized ancestry. Nest and guard bees of each type were collected and exposed to an electric stimulus of 0.5 mA, and the time they took to sting a leather substrate was recorded. Africanized bees consistently had significant lower thresholds of defensive response than European bees across all of the environments tested. Guards were faster to sting than nest bees only for the Africanized genotype, suggesting that alleles of African origin have pleiotropic effects on guarding and stinging. This is the first study that shows that single individuals specialized in guarding also may have a lower response threshold for stinging. Environmental effects were also evident. In all cases, bees responded faster to the electrical stimulation after being kept in environments other than their natal nest. Moreover, significant genotype by environment and genotype by task specialization interactions were found. Our results fit a model of division of labor based on differences in response thresholds to stimuli among workers of different genotypes and task groups that result in non-additive effects on colony behavior.

The honeybee has been the most important insect species for study of social behavior. The recently released draft genomic sequence for the bee will accelerate honeybee behavioral genetics. Although we lack sufficient tools to manipulate this genome easily, quantitative trait loci (QTLs) that influence natural variation in behavior have been identified and tested for their effects on correlated behavioral traits. We review what is known about the genetics and physiology of two behavioral traits in honeybees, foraging specialization (pollen versus nectar), and defensive behavior, and present evidence that map-based cloning of genes is more feasible in the bee than in other metazoans. We also present bioinformatic analyses of candidate genes within QTL confidence intervals (CIs). The high recombination rate of the bee made it possible to narrow the search to regions containing only 17-61 predicted peptides for each QTL, although CIs covered large genetic distances. Knowledge of correlated behavioral traits, comparative bioinformatics, and expression assays facilitated evaluation of candidate genes. An overrepresentation of genes involved in ovarian development and insulin-like signaling components within pollen foraging QTL regions suggests that an ancestral reproductive gene network was co-opted during the evolution of foraging specialization. The major QTL influencing defensive/aggressive behavior contains orthologs of genes involved in central nervous system activity and neurogenesis. Candidates at the other two defensive-behavior QTLs include modulators of sensory signaling (Am5HT(7) serotonin receptor, AmArr4 arrestin, and GABA-B-R1 receptor). These studies are the first step in linking natural variation in honeybee social behavior to the identification of underlying genes.

The first draft of the honey bee genome sequence and improved genetic maps are utilized to analyze a genome displaying 10 times higher levels of recombination (19 cM/Mb) than previously analyzed genomes of higher eukaryotes. The exceptionally high recombination rate is distributed genome-wide, but varies by two orders of magnitude. Analysis of chromosome, sequence, and gene parameters with respect to recombination showed that local recombination rate is associated with distance to the telomere, GC content, and the number of simple repeats as described for low-recombining genomes. Recombination rate does not decrease with chromosome size. On average 5.7 recombination events per chromosome pair per meiosis are found in the honey bee genome. This contrasts with a wide range of taxa that have a uniform recombination frequency of about 1.6 per chromosome pair. The excess of recombination activity does not support a mechanistic role of recombination in stabilizing pairs of homologous chromosome during chromosome pairing. Recombination rate is associated with gene size, suggesting that introns are larger in regions of low recombination and may improve the efficacy of selection in these regions. Very few transposons and no retrotransposons are present in the high-recombining genome. We propose evolutionary explanations for the exceptionally high genome-wide recombination rate.

Background Meiotic recombination has traditionally been explained based on the structural requirement to stabilize homologous chromosome pairs to ensure their proper meiotic segregation. Competing hypotheses seek to explain the emerging findings of significant heterogeneity in recombination rates within and between genomes, but intraspecific comparisons of genome-wide recombination patterns are rare. The honey bee (Apis mellifera) exhibits the highest rate of genomic recombination among multicellular animals with about five cross-over events per chromatid. Results Here, we present a comparative analysis of recombination rates across eight genetic linkage maps of the honey bee genome to investigate which genomic sequence features are correlated with recombination rate and with its variation across the eight data sets, ranging in average marker spacing ranging from 1 Mbp to 120 kbp. Overall, we found that GC content explained best the variation in local recombination rate along chromosomes at the analyzed 100 kbp scale. In contrast, variation among the different maps was correlated to the abundance of microsatellites and several specific tri- and tetra-nucleotides. Conclusions The combined evidence from eight medium-scale recombination maps of the honey bee genome suggests that recombination rate variation in this highly recombining genome might be due to the DNA configuration instead of distinct sequence motifs. However, more fine-scale analyses are needed. The empirical basis of eight differing genetic maps allowed for robust conclusions about the correlates of the local recombination rates and enabled the study of the relation between DNA features and variability in local recombination rates, which is particularly relevant in the honey bee genome with its exceptionally high recombination rate. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1281-2) contains supplementary material, which is available to authorized users.

The complete mitochondrial genome from an Africanized honey bee population (AHB, derived from Apis mellifera scutellata) was assembled and analyzed. The mitogenome is 16,411 bp long and contains the same gene repertoire and gene order as the European honey bee (13 protein coding genes, 22 tRNA genes and 2 rRNA genes). ND4 appears to use an alternate start codon and the long rRNA gene is 48 bp shorter in AHB due to a deletion in a terminal AT dinucleotide repeat. The dihydrouracil arm is missing from tRNA-Ser (AGN) and tRNA-Glu is missing the TV loop. The A + T content is comparable to the European honey bee (84.7%), which increases to 95% for the 3rd position in the protein coding genes.

We demonstrate the effects of a new quantitative trait locus (QTL), designated pln3, that was mapped in a backcross population derived from strains of bees selected for the amount of pollen they store in combs. We independently confirmed pln3 by demonstrating its effects on individual foraging behavior, as we did previously for QTLs pln1 and pln2 (Hunt et al. 1995). QTL pln2 is very robust in its effects on foraging behavior. In this study, pln2 was again shown to affect individual foraging behavior of workers derived from a hybrid backcross of the selected strains. In addition, pln2 was shown to affect the amount of pollen stored in combs of colonies derived from a wide cross of European and Africanized honeybees. This is noteworthy because it demonstrates that we can map QTLs for behavior in interstrain crosses derived from selective breeding and study their effects in unselected, natural populations. The results we present also demonstrate the repeatability of finding QTLs with measurable effects, even after outcrossing selected strains, suggesting that there is a relatively small subset of QTLs with major effects segregating in the population from which we selected our founding breeding populations. The different QTLs, pln1, pln2, and pln3, appear to have different effects, revealing the complex genetic architecture of honeybee foraging behavior.

Background: The first generation of genome sequence assemblies and annotations have had a significant impact upon our understanding of the biology of the sequenced species, the phylogenetic relationships among species, the study of populations within and across species, and have informed the biology of humans. As only a few Metazoan genomes are approaching finished quality (human, mouse, fly and worm), there is room for improvement of most genome assemblies. The honey bee (Apis mellifera) genome, published in 2006, was noted for its bimodal GC content distribution that affected the quality of the assembly in some regions and for fewer genes in the initial gene set (OGSv1.0) compared to what would be expected based on other sequenced insect genomes. Results: Here, we report an improved honey bee genome assembly (Amel_4.5) with a new gene annotation set (OGSv3.2), and show that the honey bee genome contains a number of genes similar to that of other insect genomes, contrary to what was suggested in OGSv1.0. The new genome assembly is more contiguous and complete and the new gene set includes 5000 more protein-coding genes, 50% more than previously reported. About 1/6 of the additional genes were due to improvements to the assembly, and the remaining were inferred based on new RNAseq and protein data. Conclusions: Lessons learned from this genome upgrade have important implications for future genome sequencing projects. Furthermore, the improvements significantly enhance genomic resources for the honey bee, a key model for social behavior and essential to global ecology through pollination. Funding for the project was provided by a grant to RG from the National Human Genome Research Institute, National Institutes of Health (NHGRI, NIH) U54 HG003273. Contributions from members of the CGE lab were supported by Agriculture and Food Research Initiative Competitive grant no. 2010-65205-20407 from the USDA National Institute of Food Agriculture. AKB was supported by a Clare Luce Booth Fellowship at Georgetown University

A linkage map was constructed for the honey bee based on the segregation of 365 random amplified polymorphic DNA (RAPD) markers in haploid male progeny of a single female bee. The X locus for sex determination and genes for black body color and malate dehydrogenase were mapped to separate linkage groups. RAPD markers were very efficient for mapping, with an average of about 2.8 loci mapped for each 10-nucleotide primer that was used in polymerase chain reactions. The mean interval size between markers on the map was 9.1 cM. The map covered 3110 cM of linked markers on 26 linkage groups. We estimate the total genome size to be approximately 3450 cM. The size of the map indicated a very high recombination rate for the honey bee. The relationship of physical to genetic distance was estimated at 52 kb/cM, suggesting that map-based cloning of genes will be feasible for this species.

Varroa mites (V. destructor) are a major threat to honey bees (Apis melilfera) and beekeeping worldwide and likely lead to colony decline if colonies are not treated. Most treatments involve chemical control of the mites; however, Varroa has evolved resistance to many of these miticides, leaving beekeepers with a limited number of alternatives. A non-chemical control method is highly desirable for numerous reasons including lack of chemical residues and decreased likelihood of resistance. Varroa sensitive hygiene behavior is one of two behaviors identified that are most important for controlling the growth of Varroa populations in bee hives. To identify genes influencing this trait, a study was conducted to map quantitative trait loci (QTL). Individual workers of a backcross family were observed and evaluated for their VSH behavior in a mite-infested observation hive. Bees that uncapped or removed pupae were identified. The genotypes for 1,340 informative single nucleotide polymorphisms were used to construct a high-resolution genetic map and interval mapping was used to analyze the association of the genotypes with the performance of Varroa sensitive hygiene. We identified one major QTL on chromosome 9 (LOD score = 3.21) and a suggestive QTL on chromosome 1 (LOD = 1.95). The QTL confidence interval on chromosome 9 contains the gene 'no receptor potential A' and a dopamine receptor. 'No receptor potential A' is involved in vision and olfaction in Drosophila, and dopamine signaling has been previously shown to be required for aversive olfactory learning in honey bees, which is probably necessary for identifying mites within brood cells. Further studies on these candidate genes may allow for breeding bees with this trait using marker-assisted selection.

Populations of honey bees in North America have been experiencing high annual colony mortality for 15–20 years. Many apicultural researchers believe that introduced parasites called Varroa mites (V. destructor) are the most important factor in colony deaths. One important resistance mechanism that limits mite population growth in colonies is the ability of some lines of honey bees to groom mites from their bodies. To search for genes influencing this trait, we used an Illumina Bead Station genotyping array to determine the genotypes of several hundred worker bees at over a thousand single-nucleotide polymorphisms in a family that was apparently segregating for alleles influencing this behavior. Linkage analyses provided a genetic map with 1,313 markers anchored to genome sequence. Genotypes were analyzed for association with grooming behavior, measured as the time that individual bees took to initiate grooming after mites were placed on their thoraces. Quantitative-trait-locus interval mapping identified a single chromosomal region that was significant at the chromosome-wide level (p

Populations of honey bees and other pollinators have declined worldwide in recent years. A variety of stressors have been implicated as potential causes, including agricultural pesticides. Neonicotinoid insecticides, which are widely used and highly toxic to honey bees, have been found in previous analyses of honey bee pollen and comb material. However, the routes of exposure have remained largely undefined. We used LC/MS-MS to analyze samples of honey bees, pollen stored in the hive and several potential exposure routes associated with plantings of neonicotinoid treated maize. Our results demonstrate that bees are exposed to these compounds and several other agricultural pesticides in several ways throughout the foraging period. During spring, extremely high levels of clothianidin and thiamethoxam were found in planter exhaust material produced during the planting of treated maize seed. We also found neonicotinoids in the soil of each field we sampled, including unplanted fields. Plants visited by foraging bees (dandelions) growing near these fields were found to contain neonicotinoids as well. This indicates deposition of neonicotinoids on the flowers, uptake by the root system, or both. Dead bees collected near hive entrances during the spring sampling period were found to contain clothianidin as well, although whether exposure was oral (consuming pollen) or by contact (soil/planter dust) is unclear. We also detected the insecticide clothianidin in pollen collected by bees and stored in the hive. When maize plants in our field reached anthesis, maize pollen from treated seed was found to contain clothianidin and other pesticides; and honey bees in our study readily collected maize pollen. These findings clarify some of the mechanisms by which honey bees may be exposed to agricultural pesticides throughout the growing season. These results have implications for a wide range of large-scale annual cropping systems that utilize neonicotinoid seed treatments.

A fundamental problem in meta-analysis is how to systematically combine information from multiple statistical tests to rigorously evaluate a single overarching hypothesis. This problem occurs in systems biology when attempting to map genomic attributes to complex phenotypes such as behavior. Behavior and other complex phenotypes are influenced by intrinsic and environmental determinants that act on the transcriptome, but little is known about how these determinants interact at the molecular level. We developed an informatic technique that identifies statistically significant meta-associations between gene expression patterns and transcription factor combinations. Deploying this technique for brain transcriptome profiles from ca. 400 individual bees, we show that diverse determinants of behavior rely on shared combinations of transcription factors. These relationships were revealed only when we considered complex and variable regulatory rules, suggesting that these shared transcription factors are used in distinct ways by different determinants. This regulatory code would have been missed by traditional gene coexpression or cis -regulatory analytic methods. We expect that our meta-analysis tools will be useful for a broad array of problems in systems biology and other fields.

Background The ectoparasitic mite Varroa destructor has emerged as the primary pest of domestic honey bees (Apis mellifera). Here we present an initial survey of the V. destructor genome carried out to advance our understanding of Varroa biology and to identify new avenues for mite control. This sequence survey provides immediate resources for molecular and population-genetic analyses of Varroa-Apis interactions and defines the challenges ahead for a comprehensive Varroa genome project. Results The genome size was estimated by flow cytometry to be 565 Mbp, larger than most sequenced insects but modest relative to some other Acari. Genomic DNA pooled from ~1,000 mites was sequenced to 4.3× coverage with 454 pyrosequencing. The 2.4 Gbp of sequencing reads were assembled into 184,094 contigs with an N50 of 2,262 bp, totaling 294 Mbp of sequence after filtering. Genic sequences with homology to other eukaryotic genomes were identified on 13,031 of these contigs, totaling 31.3 Mbp. Alignment of protein sequence blocks conserved among V. destructor and four other arthropod genomes indicated a higher level of sequence divergence within this mite lineage relative to the tick Ixodes scapularis. A number of microbes potentially associated with V. destructor were identified in the sequence survey, including ~300 Kbp of sequence deriving from one or more bacterial species of the Actinomycetales. The presence of this bacterium was confirmed in individual mites by PCR assay, but varied significantly by age and sex of mites. Fragments of a novel virus related to the Baculoviridae were also identified in the survey. The rate of single nucleotide polymorphisms (SNPs) in the pooled mites was estimated to be 6.2 × 10-5per bp, a low rate consistent with the historical demography and life history of the species. Conclusions This survey has provided general tools for the research community and novel directions for investigating the biology and control of Varroa mites. Ongoing development of Varroa genomic resources will be a boon for comparative genomics of under-represented arthropods, and will further enhance the honey bee and its associated pathogens as a model system for studying host-pathogen interactions.

A prominent theory states that animal phenotypes arise by evolutionary changes in gene regulation, but the extent to which this theory holds true for behavioral evolution is not known. Because “nature and nurture” are now understood to involve hereditary and environmental influences on gene expression, we studied whether environmental influences on a behavioral phenotype, i.e., aggression, could have evolved into inherited differences via changes in gene expression. Here, with microarray analysis of honey bees, we show that aggression-related genes with inherited patterns of brain expression are also environmentally regulated. There were expression differences in the brain for hundreds of genes between the highly aggressive Africanized honey bee compared with European honey bee (EHB) subspecies. Similar results were obtained for EHB in response to exposure to alarm pheromone (which provokes aggression) and when comparing old and young bees (aggressive tendencies increase with age). There was significant overlap of the gene lists generated from these three microarray experiments. Moreover, there was statistical enrichment of several of the same cis regulatory motifs in promoters of genes on all three gene lists. Aggression shows a remarkably robust brain molecular signature regardless of whether it occurs because of inherited, age-related, or environmental (social) factors. It appears that one element in the evolution of different degrees of aggressive behavior in honey bees involved changes in regulation of genes that mediate the response to alarm pheromone.

Background Honey bees (Apis mellifera) provide a principal example of diphenic development. Excess feeding of female larvae results in queens (large reproductives). Moderate diet yields workers (small helpers). The signaling pathway that links provisioning to female developmental fate is not understood, yet we reasoned that it could include TOR (target of rapamycin), a nutrient- and energy-sensing kinase that controls organismal growth. Methodology/Principal Findings Here, the role of Apis mellifera TOR (amTOR) in caste determination is examined by rapamycin/FK506 pharmacology and RNA interference (RNAi) gene knockdown. We show that in queen-destined larvae, the TOR inhibitor rapamycin induces the development of worker characters that are blocked by the antagonist FK506. Further, queen fate is associated with elevated activity of the Apis mellifera TOR encoding gene, amTOR, and amTOR gene knockdown blocks queen fate and results in individuals with worker morphology. Conclusions/Significance A much-studied insect dimorphism, thereby, can be governed by the TOR pathway. Our results present the first evidence for a role of TOR in diphenic development, and suggest that adoption of this ancestral nutrient-sensing cascade is one evolutionary pathway for morphological caste differentiation in social insects.

Honey bee nest defense involves guard bees that specialize in olfaction-based nestmate recognition and alarm-pheromone-mediated recruitment of nestmates to sting. Stinging is influenced by visual, tactile and olfactory stimuli. Both quantitative trait locus (QTL) mapping and behavioral studies point to guarding behavior as a key factor in colony stinging response. Results of reciprocal F1 crosses show that paternally inherited genes have a greater influence on colony stinging response than maternally inherited genes. The most active alarm pheromone component, isoamyl acetate (IAA) causes increased respiration and may induce stress analgesia in bees. IAA primes worker bees for ‘fight or flight’, possibly through actions of neuropeptides and/or biogenic amines. Studies of aggression in other species lead to an expectation that octopamine or 5-HT might play a role in honey bee defensive response. Genome sequence and QTL mapping identified 128 candidate genes for three regions known to influence defensive behavior. Comparative bioinformatics suggest possible roles of genes involved in neurogenesis and central nervous system (CNS) activity, and genes involved in sensory tuning through G-protein coupled receptors (GPCRs), such as an arrestin (AmArr4) and the metabotropic GABAB receptor (GABA-B-R1).

Ce travail résulte de la collaboration de très nombreux chercheurs. Seuls les auteurs de la rubrique Physical and Genetic Mapping sont cités explicitement.; Here we report the genome sequence of the honeybee Apis mellifera, a key model for social behaviour and essential to global ecology through pollination. Compared with other sequenced insect genomes, the A. mellifera genome has high A1T and CpG contents, lacks major transposon families, evolves more slowly, and is more similar to vertebrates for circadian rhythm, RNA interference and DNA methylation genes, among others. Furthermore, A.mellifera has fewer genes for innate immunity, detoxification enzymes, cuticle-forming proteins and gustatory receptors, more genes for odorant receptors, and novel genes for nectar and pollen utilization, consistent with its ecology and social organization. Compared to Drosophila, genes in early developmental pathways differ in Apis, whereas similarities exist for functions that differ markedly, such as sex determination, brain function and behaviour. Population genetics suggests a novel African origin for the species A.mellifera and insights into whether Africanized bees spread throughout the New World via hybridization or displacement

One of the best examples of a natural behavioral syndrome is the pollen-hoarding syndrome in honeybees that ties together multiple behavioral phenotypes, ranging from foraging behavior to behavioral ontogeny and learning performance. A central behavioral factor is the bees' responsiveness to sucrose, measured as their proboscis extension reflex. This study examines the genetics of this trait in diploid worker and haploid male honeybees (drones) to learn more about the genetic architecture of the overall behavioral syndrome, using original strains selected for pollen-hoarding behavior. We show that a significant proportion of the phenotypic variability is determined by genotype in males and workers. Second, our data present overwhelming evidence for pleiotropic effects of previously identified quantitative trait loci for foraging behavior (pln-QTL) and epistatic interactions among them. Furthermore, we report on three genomic QTL scans (two reciprocal worker backcrosses and one drone hybrid population) derived from our selection strains. We present at least one significant and two putative new QTL directly affecting the sucrose response of honeybees. Thus, this study demonstrates the modular genetic architecture of behavioral syndromes in general, and elucidates the genetic architecture of the pollen-hoarding behavioral syndrome in particular. Understanding this behavioral syndrome is important for understanding the division of labor in social insects and social evolution itself.

The defensive behavior of 52 hybrid honeybee (Apis mellifera L.) colonies from four sets of crosses was studied and compared with that of European and Africanized bee colonies. Colonies containing F(1) hybrid workers were obtained through reciprocal crosses between European and Africanized bees. The total number of stings deposited by workers in a moving leather patch in 1 min was recorded. In each of the four sets of crosses, bees from hybrid colonies of Africanized paternity left more stings in leather patches than bees from hybrid colonies of European paternity. Results strongly suggest paternal effects of African origin increasing the defensive behavior of hybrid colonies. Although some degree of dominance was observed for high-defensive behavior in one of the four sets of crosses involving European paternity, most of the dominance effects reported in the literature appear to be the result of paternal effects. Several hypotheses to explain this phenomenon, as well as the implications of these effects on the fitness and breeding of honeybees are discussed.

▪ Abstract One key advantage of eusociality is shared defense of the nest, brood, and stored food; nest defense plays an important role in the biology of eusocial bees. Recent studies on honey bees, Apis mellifera, have focused on the placement of defensive activity in the overall scheme of division of labor, showing that guard bees play a unique and important role in colony defense. Alarm pheromones function in integrating defensive responses; honey bee alarm pheromone is an excellent example of a multicomponent pheromonal blend. The genetic regulation of defensive behavior is now better understood from the mapping of quantitative trait loci (QTLs) associated with variation in defensiveness. Colony defense in other eusocial bees is less well understood, but enough information is available to provide interesting comparisons between A. mellifera and other species of Apis, as well as with allodapine, halictine, bombine, and meliponine bees. These comparative studies illustrate the wide variety of evolutionary solutions to problems in colony defense in the Apoidea.

We have sequenced an 81-kb genomic region from the honey bee, Apis mellifera, associated with a quantitative trait locus (QTL) sting-2 for aggressive behavior. This sequence represents the first extensive study of the honey-bee genome structure encompassing putative genes in a QTL for a behavioral trait. Expression of 13 putative genes, as well as two transcripts that were present in a honey-bee EST database, was confirmed through reverse transcription analysis of mRNA from the honey-bee head. Whereas most transcripts exhibited little or no variation between European and Africanized honey-bee alleles, one transcript demonstrated significant nonsynonymous substitutions, deletions, and insertions. All 13 putative genes lacked similarity to known invertebrate or vertebrate proteins or transcripts. This observation may be reflective of the processes that determine the genomic evolution of an insect with social behavior and/or haplo-diploidy and are an indication of the unique nature of the honey-bee genome. These results make this sequence an invaluable research tool for the ongoing honey-bee whole-genome sequencing effort.

This study was conducted to test for the effect of three stinging behaviors QTLs (sting-1, sting-2 and sting-3) on the expression of guarding and stinging behavior of individual honey bees, and to determine if results of defensive behavior QTLs found in studies with Africanized honey bees could be extended to other populations of bees. Samples of guards, stingers, foragers and nurse bees were taken from two backcross colonies derived from a defensive colony and a gentle colony. The genotype of each bee for both types of colonies was determined for two sequence tagged site (STS) markers linked to sting-1 and for another two STSs, one linked to sting-2 and one linked to sting-3. Results showed that sting-1 had an effect on the expression of both stinging and guarding behaviors, sting-2 and sting-3 influenced the expression of guarding behavior. These results indicate that division of labor is influenced by specific QTLs. Results also show that QTLs mapped in a population of Africanized honey bees using colony level phenotypes also influenced the expression of guarding and stinging behavior of individual bees of other populations.

We analyzed the alarm pheromone components from five colonies of Africanized honeybees and three colonies of European honeybees collected in Mexico. Analyses revealed a novel alarm pheromone component that was only present in appreciable quantities in the Africanized bee samples. Analysis of the mass spectrum and subsequent synthesis confirmed that this compound is 3-methyl-2-buten-1-yl acetate (3M2BA), an unsaturated derivative of IPA. In Africanized honeybees, sampling from stings of guards showed that 3M2BA was present at levels of 0-38% the amount of isoamyl acetate (IPA). Behavioral assays from three colonies each of Africanized and European bees showed that 3M2BA recruited worker bees from hives of both Africanized bees and European bees at least as efficiently as isopentyl acetate IPA, a compound widely reported to have the highest activity for releasing alarm and stinging behavior in honeybees. However, a mixture of of 3M2BA and IPA (1:2) recruited bees more efficiently than either of the compounds alone. None of the compounds differed in their efficacy for inducing bees to pursue the observers.

The stinging and guarding components of the defensive behavior of European, Africanized, hybrid, and backcross honeybees (Apis mellifera L.) were compared and analyzed at both colony and individual levels. Hybrid and Africanized backcross colonies stung as many times as Africanized ones. European backcross colonies stung more than European bees but not as many times as Africanized or Africanized backcross colonies. The degree of dominance for the number of times that worker bees stung a leather patch was estimated to be 84.3%, 200.8%, and 145.8% for hybrid, backcross European, and backcross Africanized colonies, respectively. Additionally, both guards at the colony entrance and fast-stinging workers of one European backcross colony had a significantly higher frequency of an Africanized DNA marker allele, located near "sting1," a QTL previously implicated in stinging behavior at the colony level. However, guards and fast-stinging bees from a backcross to the Africanized parental colony did not differ from control bees in their frequency for the Africanized and European markers, as would be expected if large genetic dominance effects for sting1 exist. These results support the hypothesis that genetic dominance influences the defensive behavior of honeybees and confirm the effect of sting1 on the defensiveness of individual worker bees.

Sex determination in Hymenoptera is controlled by haplo-diploidy in which unfertilized eggs develop into fertile haploid males. A single sex determination locus with several complementary alleles was proposed for Hymenoptera [so-called complementary sex determination (CSD)]. Heterozygotes at the sex determination locus are normal, fertile females, whereas diploid zygotes that are homozygous develop into sterile males. This results in a strong heterozygote advantage, and the sex locus exhibits extreme polymorphism maintained by overdominant selection. We characterized the sex-determining region by genetic linkage and physical mapping analyses. Detailed linkage and physical mapping studies showed that the recombination rate is

A study was conducted to identify quantitative trait loci (QTLs) that affect colony-level stinging behavior and individual body size of honey bees. An F1 queen was produced from a cross between a queen of European origin and a drone descended from an African subspecies. Haploid drones from the hybrid queen were individually backcrossed to sister European queens to produce 172 colonies with backcross workers that were evaluated for tendency to sting. Random amplified polymorphic DNA markers were scored from the haploid drone fathers of these colonies. Wings of workers and drones were used as a measure of body size because Africanized bees in the Americas are smaller than European bees. Standard interval mapping and multiple QTL models were used to analyze data. One possible QTL was identified with a significant effect on tendency to sting (LOD 3.57). Four other suggestive QTLs were also observed (about LOD 1.5). Possible QTLs also were identified that affect body size and were unlinked to defensive-behavior QTLs. Two of these were significant (LOD 3.54 and 5.15).

We identified two genomic regions that affect the amount of pollen stored in honey bee colonies and influence whether foragers will collect pollen or nectar. We selected for the amount of pollen stored in combs of honey bee colonies, a colony-level trait, and then used random amplified polymorphic DNA (RAPD) markers and interval mapping procedures with data from backcross colonies to identify two quantitative trait loci (pln1 and pln2, LOD 3.1 and 2.3, respectively). Quantitative trait loci effects were confirmed in a separate cross by demonstrating the cosegregation of marker alleles with the foraging behavior of individual workers. Both pln1 and pln2 had an effect on the amount of pollen carried by foragers returning to the colony, as inferred by the association between linked RAPD marker alleles, D8-.3f and 301-.55, and the individual pollen load weights of returning foragers. The alleles of the two marker loci were nonrandomly distributed with respect to foraging task. The two loci appeared to have different effects on foraging behavior. Individuals with alternative alleles for the marker linked to pln2 (but not pln1) differed with respect to the nectar sugar concentration of their nectar loads.

A colony-level phenotype was used to map the major sex determination locus (designated X) in the honey bee (Apis mellifera). Individual queen bees (reproductive females) were mated to single drones (fertile males) by instrumental insemination. Haploid drone progeny of an F1 queen were each backcrossed to daughter queens from one of the parental lines. Ninety-eight of the resulting colonies containing back-cross progeny were evaluated for the trait 'low brood-viability' resulting from the production of diploid drones that were homozygous at X. DNA samples from the haploid drone fathers of these colonies were used individually in polymerase chain reactions (PCR) with 10-base primers. These reactions generated random amplified polymorphic DNA (RAPD) markers that were analyzed for cosegregation with the colony-level phenotype. One RAPD marker allele was shared by 22 of 25 drones that fathered low brood-viability colonies. The RAPD marker fragment was cloned and partially sequenced. Two primers were designed that define a sequence-tagged site (STS) for this locus. The primers amplified DNA marker fragments that cosegregated with the original RAPD marker. In order to more precisely estimate the linkage between X and the STS locus, another group of bees consisting of progeny from one of the low-brood viability colonies was used in segregation analysis. Four diploid drones and 181 of their diploid sisters (workers, nonfertile females) were tested for segregation of the RAPD and STS markers. The cosegregating RAPD and STS markers were codominant due to the occurrence of fragment-length alleles. The four diploid drones were homozygous for these markers but only three of the 181 workers were homozygotes (recombinants). Therefore the distance between X and the STS locus was estimated at 1.6 cM. An additional linked marker was found that was 6.6 cM from the STS locus.

Honey bee (Apis mellifera) grooming behavior is an important mechanism of resistance against the parasitic mite Varroa destructor. This research was conducted to study associations between grooming behavior and the expression of selected immune, neural, detoxification, developmental and health-related genes. Individual bees tested in a laboratory assay for various levels of grooming behavior in response to V. destructor were also analyzed for gene expression. Intense groomers (IG) were most efficient in that they needed significantly less time to start grooming and fewer grooming attempts to successfully remove mites from their bodies than did light groomers (LG). In addition, the relative abundance of the neurexin-1 mRNA, was significantly higher in IG than in LG, no groomers (NG) or control (bees without mite). The abundance of poly U binding factor kd 68 and cytochrome p450 mRNAs were significantly higher in IG than in control bees. The abundance of hymenoptaecin mRNA was significantly higher in IG than in NG, but it was not different from that of control bees. The abundance of vitellogenin mRNA was not changed by grooming activity. However, the abundance of blue cheese mRNA was significantly reduced in IG compared to LG or NG, but not to control bees. Efficient removal of mites by IG correlated with different gene expression patterns in bees. These results suggest that the level of grooming behavior may be related to the expression pattern of vital honey bee genes. Neurexin-1, in particular, might be useful as a bio-marker for behavioral traits in bees.