Thrice out of Asia and the adaptive radiation of the western honey bee.

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world's most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic "hotspots." We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Defense Response in Brazilian Honey Bees (Apis mellifera scutellata × spp.) Is Underpinned by Complex Patterns of Admixture.

In 1957, an invasive and highly defensive honey bee began to spread across Brazil. In the previous year, Brazilian researchers hoped to produce a subtropical-adapted honey bee by crossing local commercial honey bees (of European origin) with a South African honey bee subspecies (Apis mellifera scutellata; an A-lineage honey bee subspecies). The resulting cross-African hybrid honey bees (AHBs)-escaped from their enclosure and spread through the Americas. Today, AHB is the most common honey bee from Northern Argentina to the Southern United States. AHBs are much more likely to sting nest intruders than managed European-derived honey bee colonies. Previous studies have explored how genetic variation contributes to differences in defense response between European-derived honey bee and AHB. Although this work demonstrated very strong genetic effects on defense response, they have yet to pinpoint which genes influence variation in defense response within AHBs, specifically. We quantified defense response for 116 colonies in Brazil and performed pooled sequencing on the most phenotypically divergent samples. We identified 65 loci containing 322 genes that were significantly associated with defense response. Loci were strongly associated with metabolic function, consistent with previous functional genomic analyses of this phenotype. Additionally, defense-associated loci had nonrandom and unexpected patterns of admixture. Defense response was not simply the product of more A-lineage honey bee ancestry as previously assumed, but rather an interaction between A-lineage and European alleles. Our results suggest that a combination of A-lineage and European alleles play roles in defensive behavior in AHBs.

Colonization history and population differentiation of the Honey Bees (Apis mellifera L.) in Puerto Rico.

Honey bees (Apis mellifera L.) are the primary commercial pollinators across the world. The subspecies A. m. scutellata originated in Africa and was introduced to the Americas in 1956. For the last 60 years, it hybridized successfully with European subspecies, previous residents in the area. The result of this hybridization was called Africanized honey bee (AHB). AHB has spread since then, arriving to Puerto Rico (PR) in 1994. The honey bee population on the island acquired a mosaic of features from AHB or the European honey bee (EHB). AHB in Puerto Rico shows a major distinctive characteristic, docile behavior, and is called gentle Africanized honey bees (gAHB). We used 917 SNPs to examine the population structure, genetic differentiation, origin, and history of range expansion and colonization of gAHB in PR. We compared gAHB to populations that span the current distribution of A. mellifera worldwide. The gAHB population is shown to be a single population that differs genetically from the examined populations of AHB. Texas and PR groups are the closest genetically. Our results support the hypothesis that the Texas AHB population is the source of gAHB in Puerto Rico.

Finding the missing honey bee genes: lessons learned from a genome upgrade.

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.

Genes associated with honey bee behavioral maturation affect clock-dependent and -independent aspects of daily rhythmic activity in fruit flies.

BACKGROUND: In the honey bee, the age-related and socially regulated transition of workers from in-hive task performance (e.g., caring for young) to foraging (provisioning the hive) is associated with changes in many behaviors including the 24-hour pattern of rhythmic activity. We have previously shown that the hive-bee to forager transition is associated with extensive changes in brain gene expression. In this study, we test the possible function of a subset of these genes in daily rhythmic activity pattern using neural-targeted RNA interference (RNAi) of an orthologous gene set in Drosophila melanogaster. PRINCIPAL FINDINGS: Of 10 genes tested, knockdown of six affected some aspect of locomotor activity under a 12 h:h light:dark regime (LD). Inos affected anticipatory activity preceding lights-off, suggesting a possible clock-dependent function. BM-40-SPARC, U2af50 and fax affected peak activity at dawn without affecting anticipation or overall inactivity (proportion of 15-min intervals without activity), suggesting that these effects may depend on the day-night light cycle. CAH1 affected overall inactivity. The remaining gene, abl, affected peak activity levels but was not clearly time-of-day-specific. No gene tested affected length of period or strength of rhythmicity in constant dark (DD), suggesting that these genes do not act in the core clock. SIGNIFICANCE: Taking advantage of Drosophila molecular genetic tools, our study provides an important step in understanding the large set of gene expression changes that occur in the honey bee transition from hive bee to forager. We show that orthologs of many of these genes influence locomotor activity in Drosophila, possibly through both clock-dependent and -independent pathways. Our results support the importance of both circadian clock and direct environmental stimuli (apart from entrainment) in shaping the bee's 24-hour pattern of activity. Our study also outlines a new approach to dissecting complex behavior in a social animal.

A genome-wide signature of positive selection in ancient and recent invasive expansions of the honey bee Apis mellifera.

Apis mellifera originated in Africa and extended its range into Eurasia in two or more ancient expansions. In 1956, honey bees of African origin were introduced into South America, their descendents admixing with previously introduced European bees, giving rise to the highly invasive and economically devastating "Africanized" honey bee. Here we ask whether the honey bee's out-of-Africa expansions, both ancient and recent (invasive), were associated with a genome-wide signature of positive selection, detected by contrasting genetic differentiation estimates (F(ST)) between coding and noncoding SNPs. In native populations, SNPs in protein-coding regions had significantly higher F(ST) estimates than those in noncoding regions, indicating adaptive evolution in the genome driven by positive selection. This signal of selection was associated with the expansion of honey bees from Africa into Western and Northern Europe, perhaps reflecting adaptation to temperate environments. We estimate that positive selection acted on a minimum of 852-1,371 genes or approximately 10% of the bee's coding genome. We also detected positive selection associated with the invasion of African-derived honey bees in the New World. We found that introgression of European-derived alleles into Africanized bees was significantly greater for coding than noncoding regions. Our findings demonstrate that Africanized bees exploited the genetic diversity present from preexisting introductions in an adaptive way. Finally, we found a significant negative correlation between F(ST) estimates and the local GC content surrounding coding SNPs, suggesting that AT-rich genes play an important role in adaptive evolution in the honey bee.

Species differences in brain gene expression profiles associated with adult behavioral maturation in honey bees.

BACKGROUND: Honey bees are known for several striking social behaviors, including a complex pattern of behavioral maturation that gives rise to an age-related colony division of labor and a symbolic dance language, by which successful foragers communicate the location of attractive food sources to their nestmates. Our understanding of honey bees is mostly based on studies of the Western honey bee, Apis mellifera, even though there are 9-10 other members of genus Apis, showing interesting variations in social behavior relative to A. mellifera. To facilitate future in-depth genomic and molecular level comparisons of behavior across the genus, we performed a microarray analysis of brain gene expression for A. mellifera and three key species found in Asia, A. cerana, A. florea and A. dorsata. RESULTS: For each species we compared brain gene expression patterns between foragers and adult one-day-old bees on an A. mellifera cDNA microarray and calculated within-species gene expression ratios to facilitate cross-species analysis. The number of cDNA spots showing hybridization fluorescence intensities above the experimental threshold was reduced by an average of 16% in the Asian species compared to A. mellifera, but an average of 71% of genes on the microarray were available for analysis. Brain gene expression profiles between foragers and one-day-olds showed differences that are consistent with a previous study on A. mellifera and were comparable across species. Although 1772 genes showed significant differences in expression between foragers and one-day-olds, only 218 genes showed differences in forager/one-day-old expression between species (p < 0.001). Principal Components Analysis revealed dominant patterns of expression that clearly distinguished between the four species but did not reflect known differences in behavior and ecology. There were species differences in brain expression profiles for functionally related groups of genes. CONCLUSION: We conclude that the A. mellifera cDNA microarray can be used effectively for cross-species comparisons within the genus. Our results indicate that there is a widespread conservation of the molecular processes in the honey bee brain underlying behavioral maturation. Species differences in brain expression profiles for functionally related groups of genes provide possible clues to the basis of behavioral variation in the genus.

Thrice out of Africa: ancient and recent expansions of the honey bee, Apis mellifera.

We characterized Apis mellifera in both native and introduced ranges using 1136 single-nucleotide polymorphisms genotyped in 341 individuals. Our results indicate that A. mellifera originated in Africa and expanded into Eurasia at least twice, resulting in populations in eastern and western Europe that are geographically close but genetically distant. A third expansion in the New World has involved the near-replacement of previously introduced "European" honey bees by descendants of more recently introduced A. m. scutellata ("African" or "killer" bees). Our analyses of spatial transects and temporal series in the New World revealed differential replacement of alleles derived from eastern versus western Europe, with admixture evident in all individuals.

Genome scan for cis-regulatory DNA motifs associated with social behavior in honey bees.

Honey bees (Apis mellifera) undergo an age-related, socially regulated transition from working in the hive to foraging, which is associated with changes in the expression of thousands of genes in the brain. To begin to study the cis-regulatory code underlying this massive social regulation of gene expression, we used the newly sequenced honey bee genome to scan the promoter regions of eight sets of behaviorally related genes differentially expressed in the brain in the context of division of labor among worker bees, for 41 cis-regulatory motifs previously characterized in Drosophila melanogaster. Binding sites for the transcription factors Hairy, GAGA, Adf1, Cf1, Snail, and Dri, known to function in nervous system development, olfactory learning, or hormone binding in Drosophila, were significantly associated with one or more gene sets. The presence of some binding sites also predicted expression patterns for as many as 71% of the genes in some gene sets. These results suggest that there is a robust relationship between cis and social regulation of brain gene expression, especially considering that we studied <15% of all known transcription factors. These results also suggest that transcriptional networks involved in the regulation of development in Drosophila are used to regulate behavioral development in adult honey bees. However, differences in gene regulation between these two processes are suggested by the finding that the promoter regions for the behaviorally related bee genes differed in both motif occurrence and G/C content relative to their Drosophila orthologs.

Genomic dissection of behavioral maturation in the honey bee.

Honey bees undergo an age-related, socially regulated transition from working in the hive to foraging that has been previously associated with changes in the expression of thousands of genes in the brain. To understand the meaning of these changes, we conducted microarray analyses to examine the following: (i) the ontogeny of gene expression preceding the onset of foraging, (ii) the effects of physiological and genetic factors that influence this behavioral transition, and (iii) the effects of foraging experience. Although >85% of approximately 5,500 genes showed brain differences, principal component analysis revealed discrete influences of age, behavior, genotype, environment, and experience. Young bees not yet competent to forage showed extensive, age-related expression changes, essentially complete by 8 days of age, coinciding with previously described structural brain changes. Subsequent changes were not age-related but were largely related to effects of juvenile hormone (JH), suggesting that the increase in JH that influences the hive bee-forager transition may cause many of these changes. Other treatments that also influence the onset age of foraging induced many changes but with little overlap, suggesting that multiple pathways affect behavioral maturation. Subspecies differences in onset age of foraging were correlated with differences in JH and JH-target gene expression, suggesting that this endocrine system mediates the genetic differences. We also used this multifactorial approach to identify candidate genes for behavioral maturation. This successful dissection of gene expression indicates that, for social behavior, gene expression in the brain can provide a robust indicator of the interaction between hereditary and environmental information.

Semiparametric approach to characterize unique gene expression trajectories across time.

BACKGROUND: A semiparametric approach was used to identify groups of cDNAs and genes with distinct expression profiles across time and overcome the limitations of clustering to identify groups. The semiparametric approach allows the generalization of mixtures of distributions while making no specific parametric assumptions about the distribution of the hidden heterogeneity of the cDNAs. The semiparametric approach was applied to study gene expression in the brains of Apis mellifera ligustica honey bees raised in two colonies (A. m. mellifera and ligustica) with consistent patterns across five maturation ages. RESULTS: The semiparametric approach provided unambiguous criteria to detect groups of genes, trajectories and probability of gene membership to groups. The semiparametric results were cross-validated in both colony data sets. Gene Ontology analysis enhanced by genome annotation helped to confirm the semiparametric results and revealed that most genes with similar or related neurobiological function were assigned to the same group or groups with similar trajectories. Ten groups of genes were identified and nine groups had highly similar trajectories in both data sets. Differences in the trajectory of the reminder group were consistent with reports of accelerated maturation in ligustica colonies compared to mellifera colonies. CONCLUSION: The combination of microarray technology, genomic information and semiparametric analysis provided insights into the genomic plasticity and gene networks linked to behavioral maturation in the honey bee.

Sociogenomics: social life in molecular terms.

Spectacular progress in molecular biology, genome-sequencing projects and genomics makes this an appropriate time to attempt a comprehensive understanding of the molecular basis of social life. Promising results have already been obtained in identifying genes that influence animal social behaviour and genes that are implicated in social evolution. These findings - derived from an eclectic mix of species that show varying levels of sociality - provide the foundation for the integration of molecular biology, genomics, neuroscience, behavioural biology and evolutionary biology that is necessary for this endeavour.

Pheromone-mediated gene expression in the honey bee brain.

We tested the hypothesis that queen mandibular pheromone (QMP) causes changes in gene expression in the brain of the adult worker honey bee, and that these changes can be correlated to the downstream behavioral responses induced by QMP. In support of the first hypothesis, cage experiments revealed that QMP transiently regulated expression of several hundred genes and chronically regulated the expression of 19 genes. Several of these genes were also affected by QMP in experiments with bee colonies in the field, demonstrating robust gene regulation by pheromone. To evaluate the second hypothesis, we focused on one function of QMP: delaying the transition from working in the hive (e.g., brood care, or "nursing") to foraging. We compared the list of QMP-regulated genes with the lists of genes differentially regulated in nurse and forager brains generated in a separate study. QMP consistently activated "nursing genes" and repressed "foraging genes," suggesting that QMP may delay behavioral maturation by regulating genes in the brain that produce these behavioral states. We also report here on an ortholog of the Drosophila transcription factor kruppel homolog 1 that was strongly regulated by QMP, especially in the mushroom bodies of the bee brain. These results demonstrate chronic gene regulation by a primer pheromone and illustrate the potential of genomics to trace the actions of a pheromone from perception to action, and thereby provide insights into how pheromones regulate social life.

Gene expression profiles in the brain predict behavior in individual honey bees.

We show that the age-related transition by adult honey bees from hive work to foraging is associated with changes in messenger RNA abundance in the brain for 39% of approximately 5500 genes tested. This result, discovered using a highly replicated experimental design involving 72 microarrays, demonstrates more extensive genomic plasticity in the adult brain than has yet been shown. Experimental manipulations that uncouple behavior and age revealed that messenger RNA changes were primarily associated with behavior. Individual brain messenger RNA profiles correctly predicted the behavior of 57 out of 60 bees, indicating a robust association between brain gene expression in the individual and naturally occurring behavioral plasticity.

Annotated expressed sequence tags and cDNA microarrays for studies of brain and behavior in the honey bee.

To accelerate the molecular analysis of behavior in the honey bee (Apis mellifera), we created expressed sequence tag (EST) and cDNA microarray resources for the bee brain. Over 20,000 cDNA clones were partially sequenced from a normalized (and subsequently subtracted) library generated from adult A. mellifera brains. These sequences were processed to identify 15,311 high-quality ESTs representing 8912 putative transcripts. Putative transcripts were functionally annotated (using the Gene Ontology classification system) based on matching gene sequences in Drosophila melanogaster. The brain ESTs represent a broad range of molecular functions and biological processes, with neurobiological classifications particularly well represented. Roughly half of Drosophila genes currently implicated in synaptic transmission and/or behavior are represented in the Apis EST set. Of Apis sequences with open reading frames of at least 450 bp, 24% are highly diverged with no matches to known protein sequences. Additionally, over 100 Apis transcript sequences conserved with other organisms appear to have been lost from the Drosophila genome. DNA microarrays were fabricated with over 7000 EST cDNA clones putatively representing different transcripts. Using probe derived from single bee brain mRNA, microarrays detected gene expression for 90% of Apis cDNAs two standard deviations greater than exogenous control cDNAs. [The sequence data described in this paper have been submitted to Genbank data library under accession nos. BI502708-BI517278. The sequences are also available at http://titan.biotec.uiuc.edu/bee/honeybee_project.htm.]