Quantitative trait loci that influence the expression of guarding and stinging behaviors of individual honey bees

Seasonal variation in defense behavior in European and scutellata-hybrid honey bees (Apis mellifera) in Southern California

Nest defense in the honey bee (Apis mellifera) is a complex collective behavior modulated by various interacting social, environmental, and genetic factors. Scutellata-hybrid ("Africanized") honey bees are usually considered to be far more defensive than European honey bees which are therefore preferred for commercial and hobbyist beekeeping. In the most recent zone of scutellata hybridization, the southern USA, the degree to which this defensiveness differs among current strains, and the extent to which defensiveness varies across a season has not been measured. We quantified the levels of A. m. scutellata ancestry in colonies and conducted a seasonal assessment (May through November) of colony nest defensiveness in feral scutellata-hybrid and a popular lineage of European honey bee commonly used in managed environments (sold as A. mellifera ligustica) hives at two apiaries in Southern California. Standard measures of defensiveness were low in both scutellata-hybrid and European colonies during May. Defensiveness increased during the later months of the study in scutellata-hybrid colonies. Most measures of defensiveness did not increase in managed colonies. Defensiveness in the scutellata-hybrids appears lower than what has been previously documented in Brazil and Mexico, possibly due to their lower proportion of A. m. scutellata ancestry.

Identification of quantitative trait loci associated with calmness and gentleness in honey bees using whole-genome sequences

The identification of quantitative trait loci (QTL) through genome-wide association studies (GWAS) is a powerful method for unravelling the genetic background of selected traits and improving early-stage predictions. In honey bees (Apis mellifera), past genetic analyses have particularly focused on individual queens and workers. In this study, we used pooled whole-genome sequences to ascertain the genetic variation of the entire colony. In total, we sampled 216 Apis mellifera mellifera and 28 Apis mellifera carnica colonies. Different experts subjectively assessed the gentleness and calmness of the colonies using a standardised protocol. Conducting a GWAS for calmness on 211 purebred A. m. mellifera colonies, we identified three QTL, on chromosomes 8, 6, and 12. The two first QTL correspond to LOC409692 gene, coding for a disintegrin and metalloproteinase domain-containing protein 10, and to Abscam gene, coding for a Dscam family member Abscam protein, respectively. The last gene has been reported to be involved in the domestication of A. mellifera. The third QTL is located 13 kb upstream of LOC102655631, coding for a trehalose transporter. For gentleness, two QTL were identified on chromosomes 4 and 3. They are located within gene LOC413669, coding for a lap4 protein, and gene LOC413416, coding for a bicaudal C homolog 1-B protein, respectively. The identified positional candidate genes of both traits mainly affect the olfaction and nervous system of honey bees. Further research is needed to confirm the results and to better understand the genetic and phenotypic basis of calmness and gentleness.

Selection and hybridization shaped the rapid spread of African honey bee ancestry in the Americas

Recent biological invasions offer 'natural' laboratories to understand the genetics and ecology of adaptation, hybridization, and range limits. One of the most impressive and well-documented biological invasions of the 20th century began in 1957 when Apis mellifera scutellata honey bees swarmed out of managed experimental colonies in Brazil. This newly-imported subspecies, native to southern and eastern Africa, both hybridized with and out-competed previously-introduced European honey bee subspecies. Populations of scutellata-European hybrid honey bees rapidly expanded and spread across much of the Americas in less than 50 years. We use broad geographic sampling and whole genome sequencing of over 300 bees to map the distribution of scutellata ancestry where the northern and southern invasions have presently stalled, forming replicated hybrid zones with European bee populations in California and Argentina. California is much farther from Brazil, yet these hybrid zones occur at very similar latitudes, consistent with the invasion having reached a climate barrier. At these range limits, we observe genome-wide clines for scutellata ancestry, and parallel clines for wing length that span hundreds of kilometers, supporting a smooth transition from climates favoring scutellata-European hybrid bees to climates where they cannot survive winter. We find no large effect loci maintaining exceptionally steep ancestry transitions. Instead, we find most individual loci have concordant ancestry clines across South America, with a build-up of somewhat steeper clines in regions of the genome with low recombination rates, consistent with many loci of small effect contributing to climate-associated fitness trade-offs. Additionally, we find no substantial reductions in genetic diversity associated with rapid expansions nor complete dropout of scutellata ancestry at any individual loci on either continent, which suggests that the competitive fitness advantage of scutellata ancestry at lower latitudes has a polygenic basis and that scutellata-European hybrid bees maintained large population sizes during their invasion. To test for parallel selection across continents, we develop a null model that accounts for drift in ancestry frequencies during the rapid expansion. We identify several peaks within a larger genomic region where selection has pushed scutellata ancestry to high frequency hundreds of kilometers past the present cline centers in both North and South America and that may underlie high-fitness traits driving the invasion.

Association between molecular markers and behavioral phenotypes in the immatures of a butterfly

Newly hatched caterpillars of the butterfly Heliconius erato phyllis routinely cannibalize eggs. In a manifestation of kin recognition they cannibalize sibling eggs less frequently than unrelated eggs. Previous work has estimated the heritability of kin recognition in H. erato phyllis to lie between 14 and 48%. It has furthermore been shown that the inheritance of kin recognition is compatible with a quantitative model with a threshold. Here we present the results of a preliminary study, in which we tested for associations between behavioral kin recognition phenotypes and AFLP and SSR markers. We implemented two experimental approaches: (1) a cannibalism test using sibling eggs only, which allowed for only two behavioral outcomes (cannibal and non-cannibal), and (2) a cannibalism test using two sibling eggs and one unrelated egg, which allowed four outcomes [cannibal who does not recognize siblings, cannibal who recognizes siblings, "super-cannibal" (cannibal of both eggs), and "super non-cannibal" (does not cannibalize eggs at all)]. Single-marker analyses were performed using χ2 tests and logistic regression with null markers as covariates. Results of the χ2 tests identified 72 associations for experimental design 1 and 73 associations for design 2. Logistic regression analysis of the markers found to be significant in the χ2 test resulted in 20 associations for design 1 and 11 associations for design 2. Experiment 2 identified markers that were more frequently present or absent in cannibals who recognize siblings and super non-cannibals; i.e. in both phenotypes capable of kin recognition.

Genomic and transcriptomic analysis of the Asian honeybee Apis cerana provides novel insights into honeybee biology

The Asian honeybee Apis cerana is one of two bee species that have been commercially kept with immense economic value. Here we present the analysis of genomic sequence and transcriptomic exploration for A. cerana as well as the comparative genomic analysis of the Asian honeybee and the European honeybee A. mellifera. The genome and RNA-seq data yield new insights into the behavioral and physiological resistance to the parasitic mite Varroa the evolution of antimicrobial peptides, and the genetic basis for labor division in A. cerana. Comparison of genes between the two sister species revealed genes specific to A. cerana, 54.5% of which have no homology to any known proteins. The observation that A. cerana displayed significantly more vigilant grooming behaviors to the presence of Varroa than A. mellifera in conjunction with gene expression analysis suggests that parasite-defensive grooming in A. cerana is likely triggered not only by exogenous stimuli through visual and olfactory detection of the parasite, but also by genetically endogenous processes that periodically activates a bout of grooming to remove the ectoparasite. This information provides a valuable platform to facilitate the traits unique to A. cerana as well as those shared with other social bees for health improvement.

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M-lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C-lineage). Effective conservation actions require reliable molecular tools to identify pure-bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q-values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.

Biased Allele Expression and Aggression in Hybrid Honeybees may be Influenced by Inappropriate Nuclear-Cytoplasmic Signaling

Hybrid effects are often exhibited asymmetrically between reciprocal families. One way this could happen is if silencing of one parent’s allele occurs in one lineage but not the other, which could affect the phenotypes of the hybrids asymmetrically by silencing that allele in only one of the hybrid families. We have previously tested for allele-specific expression biases in hybrids of European and Africanized honeybees and we found that there was an asymmetric overabundance of genes showing a maternal bias in the family with a European mother. Here, we further analyze allelic bias in these hybrids to ascertain whether they may underlie previously described asymmetries in metabolism and aggression in similar hybrid families and we speculate on what mechanisms may produce this biased allele usage. We find that there are over 500 genes that have some form of biased allele usage and over 200 of these are biased toward the maternal allele but only in the family with European maternity, mirroring the pattern observed for aggression and metabolic rate. This asymmetrically biased set is enriched for genes in loci associated with aggressive behavior and also for mitochondrial-localizing proteins. It contains many genes that play important roles in metabolic regulation. Moreover we find genes relating to the piwi-interacting RNA (piRNA) pathway, which is involved in chromatin modifications and epigenetic regulation and may help explain the mechanism underlying this asymmetric allele use. Based on these findings and previous work investigating aggression and metabolism in bees, we propose a novel hypothesis; that the asymmetric pattern of biased allele usage in these hybrids is a result of inappropriate use of piRNA-mediated nuclear-cytoplasmic signaling that is normally used to modulate aggression in honeybees. This is the first report of widespread asymmetric effects on allelic expression in hybrids and may represent a novel mechanism for gene regulation.

Genes underlying altruism

William D. Hamilton postulated the existence of 'genes underlying altruism', under the rubric of inclusive fitness theory, a half-century ago. Such genes are now poised for discovery. In this article, we develop a set of intuitive criteria for the recognition and analysis of genes for altruism and describe the first candidate genes affecting altruism from social insects and humans. We also provide evidence from a human population for genetically based trade-offs, underlain by oxytocin-system polymorphisms, between alleles for altruism and alleles for non-social cognition. Such trade-offs between self-oriented and altruistic behaviour may influence the evolution of phenotypic diversity across all social animals.

Understanding the relationship between brain gene expression and social behavior: lessons from the honey bee

Behavior is a complex phenotype that is plastic and evolutionarily labile. The advent of genomics has revolutionized the field of behavioral genetics by providing tools to quantify the dynamic nature of brain gene expression in relation to behavioral output. The honey bee Apis mellifera provides an excellent platform for investigating the relationship between brain gene expression and behavior given both the remarkable behavioral repertoire expressed by members of its intricate society and the degree to which behavior is influenced by heredity and the social environment. Here, we review a linked series of studies that assayed changes in honey bee brain transcriptomes associated with natural and experimentally induced changes in behavioral state. These experiments demonstrate that brain gene expression is closely linked with behavior, that changes in brain gene expression mediate changes in behavior, and that the association between specific genes and behavior exists over multiple timescales, from physiological to evolutionary.

A genetic analysis of the stinging and guarding behaviors of the honey bee

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.

Genetics and neurobiology of aggression in Drosophila

Aggressive behavior is widely present throughout the animal kingdom and is crucial to ensure survival and reproduction. Aggressive actions serve to acquire territory, food, or mates and in defense against predators or rivals; while in some species these behaviors are involved in establishing a social hierarchy. Aggression is a complex behavior, influenced by a broad range of genetic and environmental factors. Recent studies in Drosophila provide insight into the genetic basis and control of aggression. The state of the art on aggression in Drosophila and the many opportunities provided by this model organism to unravel the genetic and neurobiological basis of aggression are reviewed.

Honey bee aggression supports a link between gene regulation and behavioral evolution

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.

Identification of Quantitative Trait Loci and candidate genes influencing ethanol sensitivity in honey bees

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.

Characterization of honey bee sensitivity to ethanol vapor and its correlation with aggression

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.

Flight and fight: a comparative view of the neurophysiology and genetics of honey bee defensive behavior

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 GABA(B) receptor (GABA-B-R1).

Behavioral genomics of honeybee foraging and nest defense

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

Molecular marker systems in insects: current trends and future avenues

Insects comprise the largest species composition in the entire animal kingdom and possess a vast undiscovered genetic diversity and gene pool that can be better explored using molecular marker techniques. Current trends of application of DNA marker techniques in diverse domains of insect ecological studies show that mitochondrial DNA (mtDNA), microsatellites, random amplified polymorphic DNA (RAPD), expressed sequence tags (EST) and amplified fragment length polymorphism (AFLP) markers have contributed significantly for progresses towards understanding genetic basis of insect diversity and for mapping medically and agriculturally important genes and quantitative trait loci in insect pests. Apart from these popular marker systems, other novel approaches including transposon display, sequence-specific amplification polymorphism (S-SAP), repeat-associated polymerase chain reaction (PCR) markers have been identified as alternate marker systems in insect studies. Besides, whole genome microarray and single nucleotide polymorphism (SNP) assays are becoming more popular to screen genome-wide polymorphisms in fast and cost effective manner. However, use of such methodologies has not gained widespread popularity in entomological studies. The current study highlights the recent trends of applications of molecular markers in insect studies and explores the technological advancements in molecular marker tools and modern high throughput genotyping methodologies that may be applied in entomological researches for better understanding of insect ecology at molecular level.

Paternal effects on the defensive behavior of honeybees

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.

Genomic analysis in the sting-2 quantitative trait locus for defensive behavior in the honey bee, Apis mellifera

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.

Defensive behavior of honey bees: organization, genetics, and comparisons with other bees

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.