Effects of microRNA-305 knockdown on brain gene expression associated with division of labor in honey bee colonies (Apis mellifera)

Division of labor in honey bee colonies is based on the behavioral maturation of adult workers that involves a transition from working in the hive to foraging. This behavioral maturation is associated with distinct task-related transcriptomic profiles in the brain and abdominal fat body that are related to multiple regulatory factors including juvenile hormone (JH) and queen mandibular pheromone (QMP). A prominent physiological feature associated with behavioral maturation is a loss of abdominal lipid mass as bees transition to foraging. We used transcriptomic and physiological analyses to study whether microRNAs (miRNAs) are involved in the regulation of division of labor. We first identified two miRNAs that showed patterns of expression associated with behavioral maturation, ame-miR-305-5p and ame-miR-375-3p. We then downregulated the expression of these two miRNAs with sequence-specific antagomirs. Neither ame-miR-305-5p nor ame-miR-375-3p knockdown in the abdomen affected abdominal lipid mass on their own. Similarly, knockdown of ame-miR-305-5p in combination with JH or QMP also did not affect lipid mass. By contrast, ame-miR-305-5p knockdown in the abdomen caused substantial changes in gene expression in the brain. Brain gene expression changes included genes encoding transcription factors previously implicated in behavioral maturation. The results of these functional genomic experiments extend previous correlative associations of microRNAs with honey bee division of labor and point to specific roles for ame-miR-305-5p.

The genomic case against genetic determinism

Animal studies reveal that the molecular wiring of the brain can be altered by heredity, the environment, and their interaction. A deeper molecular understanding of these interactions could be a potent antidote to societal concerns of genetic determinism for human behavior, but this requires a paradigm that extends beyond traditional genome-wide association study (GWAS).

Gut microbiota contribute to variations in honey bee foraging intensity

Gut microbiomes are increasingly recognized for mediating diverse biological aspects of their hosts, including complex behavioral phenotypes. Although many studies have reported that experimental disruptions to the gut microbial community result in atypical host behavior, studies that address how gut microbes contribute to adaptive behavioral trait variation are rare. Eusocial insects represent a powerful model to test this, because of their simple gut microbiota and complex division of labor characterized by colony-level variation in behavioral phenotypes. Although previous studies report correlational differences in gut microbial community associated with division of labor, here, we provide evidence that gut microbes play a causal role in defining differences in foraging behavior between European honey bees (Apis mellifera). We found that gut microbial community structure differed between hive-based nurse bees and bees that leave the hive to forage for floral resources. These differences were associated with variation in the abundance of individual microbes, including Bifidobacterium asteroides, Bombilactobacillus mellis, and Lactobacillus melliventris. Manipulations of colony demography and individual foraging experience suggested that differences in gut microbial community composition were associated with task experience. Moreover, single-microbe inoculations with B. asteroides, B. mellis, and L. melliventris caused effects on foraging intensity. These results demonstrate that gut microbes contribute to division of labor in a social insect, and support a role of gut microbes in modulating host behavioral trait variation.

Single-cell dissection of aggression in honeybee colonies

Understanding how genotypic variation results in phenotypic variation is especially difficult for collective behaviour because group phenotypes arise from complex interactions among group members. A genome-wide association study identified hundreds of genes associated with colony-level variation in honeybee aggression, many of which also showed strong signals of positive selection, but the influence of these 'colony aggression genes' on brain function was unknown. Here we use single-cell (sc) transcriptomics and gene regulatory network (GRN) analyses to test the hypothesis that genetic variation for colony aggression influences individual differences in brain gene expression and/or gene regulation. We compared soldiers, which respond to territorial intrusion with stinging attacks, and foragers, which do not. Colony environment showed stronger influences on soldier-forager differences in brain gene regulation compared with brain gene expression. GRN plasticity was strongly associated with colony aggression, with larger differences in GRN dynamics detected between soldiers and foragers from more aggressive relative to less aggressive colonies. The regulatory dynamics of subnetworks composed of genes associated with colony aggression genes were more strongly correlated with each other across different cell types and brain regions relative to other genes, especially in brain regions involved with olfaction and vision and multimodal sensory integration, which are known to mediate bee aggression. These results show how group genetics can shape a collective phenotype by modulating individual brain gene regulatory network architecture.

The Influences of Illumination Regime on Egg-laying Rhythms of Honey Bee Queens

Honey bee queens show extreme fecundity, commonly laying more than a thousand eggs in a single day. It has proven challenging to study the temporal organization of egg-laying behavior because queens are typically active around the clock in the dark cavity of a densely populated nest. To contend with this challenge, we developed two novel methods allowing detailed monitoring of queen activity and egg laying. We first adapted a high-resolution, continuous, tracking system allowing to track the position of barcode-tagged queens in observation hives with colonies foraging outside. We found that the queen is active ~96% of the day with typically no diurnal rhythm. Next, we developed a new laboratory procedure to monitor egg laying at single egg resolution under different light regimes. We found that under constant darkness (DD) and temperature conditions, queens laid eggs with no circadian rhythms. Queen fecundity was severely reduced under constant light (LL). Under a 12:12 illumination regime, queen fecundity was comparable to under constant darkness, with a higher number of eggs during the light phase. These daily rhythms in egg laying continued when these queens were released to DD conditions, suggesting that egg-laying rhythms are influenced by endogenous circadian clocks. These results suggest that honey bee queens are active and lay eggs around the clock with no diurnal rhythms. Light has complex influences on these behaviors, but more studies are needed to determine whether these effects reflect the influence of light directly on the queen or indirectly by affecting workers that interact with the queen.

Behavioral genetics and genomics: Mendel's peas, mice, and bees

The question of the heritability of behavior has been of long fascination to scientists and the broader public. It is now widely accepted that most behavioral variation has a genetic component, although the degree of genetic influence differs widely across behaviors. Starting with Mendel's remarkable discovery of "inheritance factors," it has become increasingly clear that specific genetic variants that influence behavior can be identified. This goal is not without its challenges: Unlike pea morphology, most natural behavioral variation has a complex genetic architecture. However, we can now apply powerful genome-wide approaches to connect variation in DNA to variation in behavior as well as analyses of behaviorally related variation in brain gene expression, which together have provided insights into both the genetic mechanisms underlying behavior and the dynamic relationship between genes and behavior, respectively, in a wide range of species and for a diversity of behaviors. Here, we focus on two systems to illustrate both of these approaches: the genetic basis of burrowing in deer mice and transcriptomic analyses of division of labor in honey bees. Finally, we discuss the troubled relationship between the field of behavioral genetics and eugenics, which reminds us that we must be cautious about how we discuss and contextualize the connections between genes and behavior, especially in humans.

Context-dependent influence of threat on honey bee social network dynamics and brain gene expression

Adverse social experience affects social structure by modifying the behavior of individuals, but the relationship between an individual's behavioral state and its response to adversity is poorly understood. We leveraged naturally occurring division of labor in honey bees and studied the biological embedding of environmental threat using laboratory assays and automated behavioral tracking of whole colonies. Guard bees showed low intrinsic levels of sociability compared with foragers and nurse bees, but large increases in sociability following exposure to a threat. Threat experience also modified the expression of caregiving-related genes in a brain region called the mushroom bodies. These results demonstrate that the biological embedding of environmental experience depends on an individual's societal role and, in turn, affects its future sociability.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Neural and Molecular Mechanisms of Biological Embedding of Social Interactions

Animals operate in complex environments, and salient social information is encoded in the nervous system and then processed to initiate adaptive behavior. This encoding involves biological embedding, the process by which social experience affects the brain to influence future behavior. Biological embedding is an important conceptual framework for understanding social decision-making in the brain, as it encompasses multiple levels of organization that regulate how information is encoded and used to modify behavior. The framework we emphasize here is that social stimuli provoke short-term changes in neural activity that lead to changes in gene expression on longer timescales. This process, simplified-neurons are for today and genes are for tomorrow-enables the assessment of the valence of a social interaction, an appropriate and rapid response, and subsequent modification of neural circuitry to change future behavioral inclinations in anticipation of environmental changes. We review recent research on the neural and molecular basis of biological embedding in the context of social interactions, with a special focus on the honeybee.

The genomic basis of evolutionary differentiation among honey bees

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

Neurodevelopmental and transcriptomic effects of CRISPR/Cas9-induced somatic orco mutation in honey bees

In insects, odorant receptors facilitate olfactory communication and require the functionality of the highly conserved co-receptor gene orco. Genome editing studies in a few species of ants and moths have revealed that orco can also have a neurodevelopmental function, in addition to its canonical role in adult olfaction, discovered first in Drosophila melanogaster. To extend this analysis, we determined whether orco mutations also affect the development of the adult brain of the honey bee Apis mellifera, an important model system for social behavior and chemical communication. We used CRISPR/Cas9 to knock out orco and examined anatomical and molecular consequences. To increase efficiency, we coupled embryo microinjection with a laboratory egg collection and in vitro rearing system. This new workflow advances genomic engineering technologies in honey bees by overcoming restrictions associated with field studies. We used Sanger sequencing to quickly select individuals with complete orco knockout for neuroanatomical analyses and later validated and described the mutations with amplicon sequencing. Mutant bees had significantly fewer glomeruli, smaller total volume of all the glomeruli, and higher mean individual glomerulus volume in the antennal lobe compared to wild-type controls. RNA-Sequencing revealed that orco knockout also caused differential expression of hundreds of genes in the antenna, including genes related to neural development and genes encoding odorant receptors. The expression of other types of chemoreceptor genes was generally unaffected, reflecting specificity of CRISPR activity in this study. These results suggest that neurodevelopmental effects of orco are related to specific insect life histories.

Assessing Agrochemical Risk to Mated Honey Bee Queens

Current risk assessment strategies for honey bees rely heavily upon laboratory tests performed on adult or immature worker bees, but these methods may not accurately capture the effects of agrochemical exposure on honey bee queens. As the sole producer of fertilized eggs inside a honeybee colony, the queen is arguably the most important single member of a functioning colony unit. Therefore, understanding how agrochemicals affect queen health and productivity should be considered a critical aspect of pesticide risk assessment. Here, an adapted method is presented to expose honey bee queens and worker queen attendants to agrochemical stressors administered through a worker diet, followed by tracking egg production in the laboratory and assessing first instar eclosion using a specialized cage, referred to as a Queen Monitoring Cage. To illustrate the method's intended use, results of an experiment in which worker queen attendants were fed diet containing sublethal doses of imidacloprid and effects on queens were monitored are described.

Individual differences in honey bee behavior enabled by plasticity in brain gene regulatory networks

Understanding the regulatory architecture of phenotypic variation is a fundamental goal in biology, but connections between gene regulatory network (GRN) activity and individual differences in behavior are poorly understood. We characterized the molecular basis of behavioral plasticity in queenless honey bee (Apis mellifera) colonies, where individuals engage in both reproductive and non-reproductive behaviors. Using high-throughput behavioral tracking, we discovered these colonies contain a continuum of phenotypes, with some individuals specialized for either egg-laying or foraging and 'generalists' that perform both. Brain gene expression and chromatin accessibility profiles were correlated with behavioral variation, with generalists intermediate in behavior and molecular profiles. Models of brain GRNs constructed for individuals revealed that transcription factor (TF) activity was highly predictive of behavior, and behavior-associated regulatory regions had more TF motifs. These results provide new insights into the important role played by brain GRN plasticity in the regulation of behavior, with implications for social evolution.

Juvenile hormone regulates brain-reproduction tradeoff in bumble bees but not in honey bees

Gonadotropic hormones coordinate processes in diverse tissues regulating animal reproductive physiology and behavior. Juvenile hormone (JH) is the ancient and most common gonadotropin in insects, but not in advanced eusocial honey bees and some ants. To start probing the evolutionary basis of this change, we combined endocrine manipulations, transcriptomics, and behavioral analyses to study JH regulated processes in a bumble bee showing a relatively simple level of eusociality. We found that in worker fat body, more JH-regulated genes were up- rather than down-regulated, and enriched for metabolic and biosynthetic pathways. This transcriptomic pattern is consistent with earlier evidence that JH is the major gonadotropin in bumble bees. In the brain, more JH-regulated genes were down- rather than up-regulated and enriched for protein turnover pathways. Brain ribosomal protein gene expression shows a similar trend of downregulation in dominant workers, which naturally have high JH titers. In other species, similar downregulation of protein turnover is found in aging brains or under stress, associated with compromised long-term memory and health. These findings suggest a previously unknown gonadotropin-mediated tradeoff. Analysis of published data reveals no such downregulation of protein turnover pathways in the brain of honey bee workers, which exhibit more complex eusociality and in which JH is not a gonadotropin but rather regulates division of labor. These results suggest that the evolution of complex eusociality in honey bees was associated with modifications in hormonal signalling supporting extended and extremely high fertility while reducing the ancient costs of high gonadotropin titers to the brain.

Behavior-related gene regulatory networks: A new level of organization in the brain

Neuronal networks are the standard heuristic model today for describing brain activity associated with animal behavior. Recent studies have revealed an extensive role for a completely distinct layer of networked activities in the brain-the gene regulatory network (GRN)-that orchestrates expression levels of hundreds to thousands of genes in a behavior-related manner. We examine emerging insights into the relationships between these two types of networks and discuss their interplay in spatial as well as temporal dimensions, across multiple scales of organization. We discuss properties expected of behavior-related GRNs by drawing inspiration from the rich literature on GRNs related to animal development, comparing and contrasting these two broad classes of GRNs as they relate to their respective phenotypic manifestations. Developmental GRNs also represent a third layer of network biology, playing out over a third timescale, which is believed to play a crucial mediatory role between neuronal networks and behavioral GRNs. We end with a special emphasis on social behavior, discuss whether unique GRN organization and cis-regulatory architecture underlies this special class of behavior, and review literature that suggests an affirmative answer.

Involving Urban Single Low-Income African American Mothers in Genomic Research: Giving Voice to How Place Matters in Health Disparities and Prevention Strategies

This article describes the process of using principles from community-based participatory action research to involve low-income, single, African American mothers on the south side of Chicago in genomic research, including as citizen scientists. The South Chicago Black Mothers' Resiliency Project used a mixed methods design to investigate how the stress of living in neighborhoods with high levels of violence affects mothers' mental and physical health. This article seeks to serve as a model for physicians and scholars interested in successfully involving low-income African American mothers in genomic research, and other health-related activities in ways that are culturally sensitive and transformative. The lives of Black mothers who struggle under interlocking systems of oppression that are often hidden from view of most Americans are at the center of this article. Therefore, we provide extensive information about the procedures used to collect the various types of data, the rationale for our procedures, the setting, the responses of mothers in our sample and methodological challenges. This study also has implications for the current COVID-19 pandemic and the need to train a corps of citizen scientists in health and wellness to avoid future extreme loss of life such as the 106,195 lives lost in the United States as of June 1, 2020.

Meta-analysis of honey bee neurogenomic response links Deformed wing virus type A to precocious behavioral maturation

Crop pollination by the western honey bee Apis mellifera is vital to agriculture but threatened by alarmingly high levels of colony mortality, especially in Europe and North America. Colony loss is due, in part, to the high viral loads of Deformed wing virus (DWV), transmitted by the ectoparasitic mite Varroa destructor, especially throughout the overwintering period of a honey bee colony. Covert DWV infection is commonplace and has been causally linked to precocious foraging, which itself has been linked to colony loss. Taking advantage of four brain transcriptome studies that unexpectedly revealed evidence of covert DWV-A infection, we set out to explore whether this effect is due to DWV-A mimicking naturally occurring changes in brain gene expression that are associated with behavioral maturation. Consistent with this hypothesis, we found that brain gene expression profiles of DWV-A infected bees resembled those of foragers, even in individuals that were much younger than typical foragers. In addition, brain transcriptional regulatory network analysis revealed a positive association between DWV-A infection and transcription factors previously associated with honey bee foraging behavior. Surprisingly, single-cell RNA-Sequencing implicated glia, not neurons, in this effect; there are relatively few glial cells in the insect brain and they are rarely associated with behavioral plasticity. Covert DWV-A infection also has been linked to impaired learning, which together with precocious foraging can lead to increased occurrence of infected bees from one colony mistakenly entering another colony, especially under crowded modern apiary conditions. These findings provide new insights into the mechanisms by which DWV-A affects honey bee health and colony survival.

Valence of social information is encoded in different subpopulations of mushroom body Kenyon cells in the honeybee brain

Over 600 Myr of evolutionary divergence between vertebrates and invertebrates is associated with considerable neuroanatomical variation both across and within these lineages. By contrast, valence encoding is an important behavioural trait that is evolutionarily conserved across vertebrates and invertebrates, and enables individuals to distinguish between positive (potentially beneficial) and negative (potentially harmful) situations. We tested the hypothesis that social interactions of positive and negative valence are modularly encoded in the honeybee brain (i.e. encoded in different cellular subpopulations) as in vertebrate brains. In vertebrates, neural activation patterns are distributed across distinct parts of the brain, suggesting that discrete circuits encode positive or negative stimuli. Evidence for this hypothesis would suggest a deep homology of neural organization between insects and vertebrates for valence encoding, despite vastly different brain sizes. Alternatively, overlapping localization of valenced social information in the brain would imply a 're-use' of circuitry in response to positive and negative social contexts, potentially to overcome the energetic constraints of a tiny brain. We used immediate early gene expression to map positively and negatively valenced social interactions in the brain of the western honeybee Apis mellifera. We found that the valence of a social signal is represented by distinct anatomical subregions of the mushroom bodies, an invertebrate sensory neuropil associated with social behaviour, multimodal sensory integration, learning and memory. Our results suggest that the modularization of valenced social information in the brain is a fundamental property of neuroanatomical organization.

Division of labor in honey bees is associated with transcriptional regulatory plasticity in the brain

Studies in evolutionary and developmental biology show that relationships between transcription factors (TFs) and their target genes can be altered to result in novel regulatory relationships that generate phenotypic plasticity. We hypothesized that context-dependent shifts in the nervous system associated with behavior may also be linked to changes in TF-target relationships over physiological time scales. We tested this hypothesis using honey bee (Apis mellifera) division of labor as a model system by performing bioinformatic analyses of previously published brain transcriptomic profiles together with new RNAi and behavioral experiments. The bioinformatic analyses identified five TFs that exhibited strong signatures of regulatory plasticity as a function of division of labor. RNAi targeting of one of these TFs (broad complex) and a related TF that did not exhibit plasticity (fushi tarazu transcription factor 1) was administered in conjunction with automated analyses of foraging behavior in the field, laboratory assays of aggression and brood care behavior, and endocrine treatments. The results showed that changes in the regulatory relationships of these TFs were associated with behavioral state, social context and endocrine state. These findings provide the first empirical evidence that TF-target relationships in the brain are altered in conjunction with behavior and social context. They also suggest that one mechanism for this plasticity involves pleiotropic TFs high up in regulatory hierarchies producing behavior-specific transcriptional responses by activating different downstream TFs to induce discrete context-dependent transcriptional cascades. These findings provide new insights into the dynamic nature of the transcriptional regulatory architecture underlying behavior in the brain.

A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds

Background

The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map.

Results

Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features.

Conclusions

The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5642-0) contains supplementary material, which is available to authorized users.

Genetic accommodation and the role of ancestral plasticity in the evolution of insect eusociality

For over a century, biologists have proposed a role for phenotypic plasticity in evolution, providing an avenue for adaptation in addition to 'mutation-first' models of evolutionary change. According to the various versions of this idea, the ability of organisms to respond adaptively to their environment through phenotypic plasticity may lead to novel phenotypes that can be screened by natural selection. If these initially environmentally induced phenotypes increase fitness, then genetic accommodation can lead to allele frequency change, influencing the expression of those phenotypes. Despite the long history of 'plasticity-first' models, the importance of genetic accommodation in shaping evolutionary change has remained controversial - it is neither fully embraced nor completely discarded by most evolutionary biologists. We suggest that the lack of acceptance of genetic accommodation in some cases is related to a lack of information on its molecular mechanisms. However, recent reports of epigenetic transgenerational inheritance now provide a plausible mechanism through which genetic accommodation may act, and we review this research here. We also discuss current evidence supporting a role for genetic accommodation in the evolution of eusociality in social insects, which have long been models for studying the influence of the environment on phenotypic variation, and may be particularly good models for testing hypotheses related to genetic accommodation. Finally, we introduce 'eusocial engineering', a method by which novel social phenotypes are first induced by environmental modification and then studied mechanistically to understand how environmentally induced plasticity may lead to heritable changes in social behavior. We believe the time is right to incorporate genetic accommodation into models of the evolution of complex traits, armed with new molecular tools and a better understanding of non-genetic heritable elements.

Quantifying the effects of pollen nutrition on honey bee queen egg laying with a new laboratory system

Honey bee populations have been declining precipitously over the past decade, and multiple causative factors have been identified. Recent research indicates that these frequently co-occurring stressors interact, often in unpredictable ways, therefore it has become important to develop robust methods to assess their effects both in isolation and in combination. Most such efforts focus on honey bee workers, but the state of a colony also depends on the health and productivity of its queen. However, it is much more difficult to quantify the performance of queens relative to workers in the field, and there are no laboratory assays for queen performance. Here, we present a new system to monitor honey bee queen egg laying under laboratory conditions and report the results of experiments showing the effects of pollen nutrition on egg laying. These findings suggest that queen egg laying and worker physiology can be manipulated in this system through pollen nutrition, which is consistent with findings from field colonies. The results generated using this controlled, laboratory-based system suggest that worker physiology controls queen egg laying behavior. Additionally, the quantitative data generated in these experiments highlight the utility of the system for further use as a risk assessment tool.

Honey bee neurogenomic responses to affiliative and agonistic social interactions

Social interactions can be divided into two categories, affiliative and agonistic. How neurogenomic responses reflect these opposing valences is a central question in the biological embedding of experience. To address this question, we exposed honey bees to a queen larva, which evokes nursing, an affiliative alloparenting interaction, and measured the transcriptomic response of the mushroom body brain region at different times after exposure. Hundreds of genes were differentially expressed at distinct time points, revealing a dynamic temporal patterning of the response. Comparing these results to our previously published research on agonistic aggressive interactions, we found both shared and unique transcriptomic responses to each interaction. The commonly responding gene set was enriched for nuclear receptor signaling, the set specific to nursing was enriched for olfaction and neuron differentiation, and the set enriched for aggression was enriched for cytoskeleton, metabolism, and chromosome organization. Whole brain histone profiling after the affiliative interaction revealed few changes in chromatin accessibility, suggesting that the transcriptomic changes derive from already accessible areas of the genome. Although only one stimulus of each type was studied, we suggest that elements of the observed transcriptomic responses reflect molecular encoding of stimulus valence, thus priming individuals for future encounters. This hypothesis is supported by behavioral analyses showing that bees responding to either the affiliative or agonistic stimulus exhibited a higher probability of repeating the same behavior but a lower probability of performing the opposite behavior. These findings add to our understanding of the biological embedding at the molecular level.

Cross-species systems analysis of evolutionary toolkits of neurogenomic response to social challenge

Social challenges like territorial intrusions evoke behavioral responses in widely diverging species. Recent work has showed that evolutionary "toolkits"-genes and modules with lineage-specific variations but deep conservation of function-participate in the behavioral response to social challenge. Here, we develop a multispecies computational-experimental approach to characterize such a toolkit at a systems level. Brain transcriptomic responses to social challenge was probed via RNA-seq profiling in three diverged species-honey bees, mice and three-spined stickleback fish-following a common methodology, allowing fair comparisons across species. Data were collected from multiple brain regions and multiple time points after social challenge exposure, achieving anatomical and temporal resolution substantially greater than previous work. We developed statistically rigorous analyses equipped to find homologous functional groups among these species at the levels of individual genes, functional and coexpressed gene modules, and transcription factor subnetworks. We identified six orthogroups involved in response to social challenge, including groups represented by mouse genes Npas4 and Nr4a1, as well as common modulation of systems such as transcriptional regulators, ion channels, G-protein-coupled receptors and synaptic proteins. We also identified conserved coexpression modules enriched for mitochondrial fatty acid metabolism and heat shock that constitute the shared neurogenomic response. Our analysis suggests a toolkit wherein nuclear receptors, interacting with chaperones, induce transcriptional changes in mitochondrial activity, neural cytoarchitecture and synaptic transmission after social challenge. It shows systems-level mechanisms that have been repeatedly co-opted during evolution of analogous behaviors, thus advancing the genetic toolkit concept beyond individual genes.

Defense against territorial intrusion is associated with DNA methylation changes in the honey bee brain

BACKGROUND

Aggression is influenced by individual variation in temperament as well as behavioral plasticity in response to adversity. DNA methylation is stably maintained over time, but also reversible in response to specific environmental conditions, and may thus be a neuromolecular regulator of both of these processes. A previous study reported DNA methylation differences between aggressive Africanized and gentle European honey bees. We investigated whether threat-induced aggression altered DNA methylation profiles in the honey bee brain in response to a behavioral stimulus (aggression-provoking intruder bee or inert control). We sampled five minutes and two hours after stimulus exposure to examine the effect of time on epigenetic profiles of aggression.

RESULTS

There were DNA methylation differences between aggressive and control bees for individual cytosine-guanine dinucleotides (CpGs) across the genome. Eighteen individual CpG sites showed significant difference between aggressive and control bees 120 min post stimulus. For clusters of CpGs, we report four genomic regions differentially methylated between aggressive and control bees at the 5-min time point, and 50 regions differentially methylated at the120-minute time point following intruder exposure. Differential methylation occurred at genes involved in neural plasticity, chromatin remodeling and hormone signaling. Additionally, there was a significant overlap of differential methylation with previously published epigenetic differences that distinguish aggressive Africanized and gentle European honey bees, suggesting an evolutionarily conserved use of brain DNA methylation in the regulation of aggression. Lastly, we identified individually statistically suggestive CpGs that as a group were significantly associated with differentially expressed genes underlying aggressive behavior and also co-localize with binding sites of transcription factors involved in neuroplasticity or neurodevelopment.

CONCLUSIONS

There were DNA methylation differences in the brain associated with response to an intruder. These differences increased in number a few hours after the initial exposure and overlap with previously reported aggression-associated genes and neurobiologically relevant transcription factor binding sites. Many DNA methylation differences that occurred in association with the expression of aggression in real time also exist between Africanized bees and European bees, suggesting an evolutionarily conserved role for epigenetic regulation in aggressive behavior.

Automated monitoring of behavior reveals bursty interaction patterns and rapid spreading dynamics in honeybee social networks

Social networks mediate the spread of information and disease. The dynamics of spreading depends, among other factors, on the distribution of times between successive contacts in the network. Heavy-tailed (bursty) time distributions are characteristic of human communication networks, including face-to-face contacts and electronic communication via mobile phone calls, email, and internet communities. Burstiness has been cited as a possible cause for slow spreading in these networks relative to a randomized reference network. However, it is not known whether burstiness is an epiphenomenon of human-specific patterns of communication. Moreover, theory predicts that fast, bursty communication networks should also exist. Here, we present a high-throughput technology for automated monitoring of social interactions of individual honeybees and the analysis of a rich and detailed dataset consisting of more than 1.2 million interactions in five honeybee colonies. We find that bees, like humans, also interact in bursts but that spreading is significantly faster than in a randomized reference network and remains so even after an experimental demographic perturbation. Thus, while burstiness may be an intrinsic property of social interactions, it does not always inhibit spreading in real-world communication networks. We anticipate that these results will inform future models of large-scale social organization and information and disease transmission, and may impact health management of threatened honeybee populations.

A soft selective sweep during rapid evolution of gentle behaviour in an Africanized honeybee

Highly aggressive Africanized honeybees (AHB) invaded Puerto Rico (PR) in 1994, displacing gentle European honeybees (EHB) in many locations. Gentle AHB (gAHB), unknown anywhere else in the world, subsequently evolved on the island within a few generations. Here we sequence whole genomes from gAHB and EHB populations, as well as a North American AHB population, a likely source of the founder AHB on PR. We show that gAHB retains high levels of genetic diversity after evolution of gentle behaviour, despite selection on standing variation. We observe multiple genomic loci with significant signatures of selection. Rapid evolution during colonization of novel habitats can generate major changes to characteristics such as morphological or colouration traits, usually controlled by one or more major genetic loci. Here we describe a soft selective sweep, acting at multiple loci across the genome, that occurred during, and may have mediated, the rapid evolution of a behavioural trait.

Africanized honey bees (AHB) are notoriously aggressive, but in Puerto Rico they have a ‘gentle’ phenotype. Here, Avalos et al. show that there has been a soft selective sweep at several loci in the Puerto Rican AHB population and suggest a role in the rapid evolution of gentle behaviour.

Nutritional Regulation of Phenotypic Plasticity in a Solitary Bee (Hymenoptera: Megachilidae)

Phenotypic plasticity involves adaptive responses to predictable environmental fluctuations and may promote evolutionary change. We studied the regulation of phenotypic plasticity in an important agricultural pollinator, the solitary alfalfa leafcutting bee (Megachile rotundata F.). Specifically, we investigated how larval nutrition affects M. rotundata diapause plasticity and how diapause plasticity affects adult female reproductive behavior. Field surveys and laboratory manipulations of aspects of larval diet demonstrated nutritional regulation of M. rotundata diapause plasticity. Manipulation of larval diet quality through the addition of royal jelly, the caste-determining substance of the honey bee Apis mellifera L., increased the probability of diapause in M. rotundata. We also found that larval nutrition and diapause status affected M. rotundata adult female reproductive behavior. Nutritional effects on larval diapause that also impact adult fitness have intriguing implications for the evolution of developmental plasticity in bees. In particular, as the solitary lifestyle of M. rotundata is considered to be the ancestral condition in bees, nutritionally regulated plasticity may have been an ancestral condition in all bees that facilitated the evolution of other forms of phenotypic plasticity, such as the castes of social bees.

Deep evolutionary conservation of autism-related genes

E. O. Wilson proposed in Sociobiology that similarities between human and animal societies reflect common mechanistic and evolutionary roots. When introduced in 1975, this controversial hypothesis was beyond science's ability to test. We used genomic analyses to determine whether superficial behavioral similarities in humans and the highly social honey bee reflect common molecular mechanisms. Here, we report that gene expression signatures for individual bees unresponsive to various salient social stimuli are significantly enriched for autism spectrum disorder-related genes. These signatures occur in the mushroom bodies, a high-level integration center of the insect brain. Furthermore, our finding of enrichment was unique to autism spectrum disorders; brain gene expression signatures from other honey bee behaviors do not show this enrichment, nor do datasets from other human behavioral and health conditions. These results demonstrate deep conservation for genes associated with a human social pathology and individual differences in insect social behavior, thus providing an example of how comparative genomics can be used to test sociobiological theory.

Transcriptomic analysis of instinctive and learned reward-related behaviors in honey bees

We used transcriptomics to compare instinctive and learned, reward-based honey bee behaviors with similar spatio-temporal components: mating flights by males (drones) and time-trained foraging flights by females (workers), respectively. Genome-wide gene expression profiling via RNA sequencing was performed on the mushroom bodies, a region of the brain known for multi-modal sensory integration and responsive to various types of reward. Differentially expressed genes (DEGs) associated with the onset of mating (623 genes) were enriched for the gene ontology (GO) categories of Transcription, Unfolded Protein Binding, Post-embryonic Development, and Neuron Differentiation. DEGs associated with the onset of foraging (473) were enriched for Lipid Transport, Regulation of Programmed Cell Death, and Actin Cytoskeleton Organization. These results demonstrate that there are fundamental molecular differences between similar instinctive and learned behaviors. In addition, there were 166 genes with strong similarities in expression across the two behaviors - a statistically significant overlap in gene expression, also seen in Weighted Gene Co-Expression Network Analysis. This finding indicates that similar instinctive and learned behaviors also share common molecular architecture. This common set of DEGs was enriched for Regulation of RNA Metabolic Process, Transcription Factor Activity, and Response to Ecdysone. These findings provide a starting point for better understanding the relationship between instincts and learned behaviors. In addition, because bees collect food for their colony rather than for themselves, these results also support the idea that altruistic behavior relies, in part, on elements of brain reward systems associated with selfish behavior.

Caste-biased gene expression in a facultatively eusocial bee suggests a role for genetic accommodation in the evolution of eusociality

Developmental plasticity may accelerate the evolution of phenotypic novelty through genetic accommodation, but studies of genetic accommodation often lack knowledge of the ancestral state to place selected traits in an evolutionary context. A promising approach for assessing genetic accommodation involves using a comparative framework to ask whether ancestral plasticity is related to the evolution of a particular trait. Bees are an excellent group for such comparisons because caste-based societies (eusociality) have evolved multiple times independently and extant species exhibit different modes of eusociality. We measured brain and abdominal gene expression in a facultatively eusocial bee, Megalopta genalis, and assessed whether plasticity in this species is functionally linked to eusocial traits in other bee lineages. Caste-biased abdominal genes in M. genalis overlapped significantly with caste-biased genes in obligately eusocial bees. Moreover, caste-biased genes in M. genalis overlapped significantly with genes shown to be rapidly evolving in multiple studies of 10 bee species, particularly for genes in the glycolysis pathway and other genes involved in metabolism. These results provide support for the idea that eusociality can evolve via genetic accommodation, with plasticity in facultatively eusocial species like M. genalis providing a substrate for selection during the evolution of caste in obligately eusocial lineages.

Conservation in Mammals of Genes Associated with Aggression-Related Behavioral Phenotypes in Honey Bees

The emerging field of sociogenomics explores the relations between social behavior and genome structure and function. An important question is the extent to which associations between social behavior and gene expression are conserved among the Metazoa. Prior experimental work in an invertebrate model of social behavior, the honey bee, revealed distinct brain gene expression patterns in African and European honey bees, and within European honey bees with different behavioral phenotypes. The present work is a computational study of these previous findings in which we analyze, by orthology determination, the extent to which genes that are socially regulated in honey bees are conserved across the Metazoa. We found that the differentially expressed gene sets associated with alarm pheromone response, the difference between old and young bees, and the colony influence on soldier bees, are enriched in widely conserved genes, indicating that these differences have genomic bases shared with many other metazoans. By contrast, the sets of differentially expressed genes associated with the differences between African and European forager and guard bees are depleted in widely conserved genes, indicating that the genomic basis for this social behavior is relatively specific to honey bees. For the alarm pheromone response gene set, we found a particularly high degree of conservation with mammals, even though the alarm pheromone itself is bee-specific. Gene Ontology identification of human orthologs to the strongly conserved honey bee genes associated with the alarm pheromone response shows overrepresentation of protein metabolism, regulation of protein complex formation, and protein folding, perhaps associated with remodeling of critical neural circuits in response to alarm pheromone. We hypothesize that such remodeling may be an adaptation of social animals to process and respond appropriately to the complex patterns of conspecific communication essential for social organization.

Characterization of Genomic Variants Associated with Scout and Recruit Behavioral Castes in Honey Bees Using Whole-Genome Sequencing

Among forager honey bees, scouts seek new resources and return to the colony, enlisting recruits to collect these resources. Differentially expressed genes between these behaviors and genetic variability in scouting phenotypes have been reported. Whole-genome sequencing of 44 Apis mellifera scouts and recruits was undertaken to detect variants and further understand the genetic architecture underlying the behavioral differences between scouts and recruits. The median coverage depth in recruits and scouts was 10.01 and 10.7 X, respectively. Representation of bacterial species among the unmapped reads reflected a more diverse microbiome in scouts than recruits. Overall, 1,412,705 polymorphic positions were analyzed for associations with scouting behavior, and 212 significant (p-value < 0.0001) associations with scouting corresponding to 137 positions were detected. Most frequent putative transcription factor binding sites proximal to significant variants included Broad-complex 4, Broad-complex 1, Hunchback, and CF2-II. Three variants associated with scouting were located within coding regions of ncRNAs including one codon change (LOC102653644) and 2 frameshift indels (LOC102654879 and LOC102655256). Significant variants were also identified on the 5’UTR of membrin, and 3’UTRs of laccase 2 and diacylglycerol kinase theta. The 60 significant variants located within introns corresponded to 39 genes and most of these positions were > 1000 bp apart from each other. A number of these variants were mapped to ncRNA LOC100578102, solute carrier family 12 member 6-like gene, and LOC100576965 (meprin and TRAF-C homology domain containing gene). Functional categories represented among the genes corresponding to significant variants included: neuronal function, exoskeleton, immune response, salivary gland development, and enzymatic food processing. These categories offer a glimpse into the molecular support to the behaviors of scouts and recruits. The level of association between genomic variants and scouting behavior observed in this study may be linked to the honey bee’s genomic plasticity and fluidity of transition between castes.

Brain regions and molecular pathways responding to food reward type and value in honey bees

The ability of honey bees to evaluate differences in food type and value is crucial for colony success, but these assessments are made by individuals who bring food to the hive, eating little, if any, of it themselves. We tested the hypothesis that responses to food type (pollen or nectar) and value involve different subsets of brain regions, and genes responsive to food. mRNA in situ hybridization of c-jun revealed that brain regions responsive to differences in food type were mostly different from regions responsive to differences in food value, except those dorsal and lateral to the mushroom body calyces, which responded to all three. Transcriptomic profiles of the mushroom bodies generated by RNA sequencing gave the following results: (1) responses to differences in food type or value included a subset of molecular pathways involved in the response to food reward; (2) genes responsive to food reward, food type and food value were enriched for (the Gene Ontology categories) mitochondrial and endoplasmic reticulum activity; (3) genes responsive to only food and food type were enriched for regulation of transcription and translation; and (4) genes responsive to only food and food value were enriched for regulation of neuronal signaling. These results reveal how activities necessary for colony survival are channeled through the reward system of individual honey bees.

Laboratory Assay of Brood Care for Quantitative Analyses of Individual Differences in Honey Bee (Apis mellifera) Affiliative Behavior

Care of offspring is a form of affiliative behavior that is fundamental to studies of animal social behavior. Insects do not figure prominently in this topic because Drosophila melanogaster and other traditional models show little if any paternal or maternal care. However, the eusocial honey bee exhibits cooperative brood care with larvae receiving intense and continuous care from their adult sisters, but this behavior has not been well studied because a robust quantitative assay does not exist. We present a new laboratory assay that enables quantification of group or individual honey bee brood “nursing behavior” toward a queen larva. In addition to validating the assay, we used it to examine the influence of the age of the larva and the genetic background of the adult bees on nursing performance. This new assay also can be used in the future for mechanistic analyses of eusociality and comparative analyses of affilative behavior with other animals.

Diet and endocrine effects on behavioral maturation-related gene expression in the pars intercerebralis of the honey bee brain

Nervous and neuroendocrine systems mediate environmental conditions to control a variety of life history traits. Our goal was to provide mechanistic insights as to how neurosecretory signals mediate division of labor in the honey bee (Apis mellifera). Worker division of labor is based on a process of behavioral maturation by individual bees, which involves performing in-hive tasks early in adulthood, then transitioning to foraging for food outside the hive. Social and nutritional cues converge on endocrine factors to regulate behavioral maturation, but whether neurosecretory systems are central to this process is not known. To explore this, we performed transcriptomic profiling of a neurosecretory region of the brain, the pars intercerebralis (PI). We first compared PI transcriptional profiles for bees performing in-hive tasks and bees engaged in foraging. Using these results as a baseline, we then performed manipulative experiments to test whether the PI is responsive to dietary changes and/or changes in juvenile hormone (JH) levels. Results reveal a robust molecular signature of behavioral maturation in the PI, with a subset of gene expression changes consistent with changes elicited by JH treatment. In contrast, dietary changes did not induce transcriptomic changes in the PI consistent with behavioral maturation or JH treatment. Based on these results, we propose a new verbal model of the regulation of division of labor in honey bees in which the relationship between diet and nutritional physiology is attenuated, and in its place is a relationship between social signals and nutritional physiology that is mediated by JH.

Early-life experience affects honey bee aggression and resilience to immune challenge

Early-life social experiences cause lasting changes in behavior and health for a variety of animals including humans, but it is not well understood how social information ''gets under the skin'' resulting in these effects. Adult honey bees (Apis mellifera) exhibit socially coordinated collective nest defense, providing a model for social modulation of aggressive behavior. Here we report for the first time that a honey bee's early-life social environment has lasting effects on individual aggression: bees that experienced high-aggression environments during pre-adult stages showed increased aggression when they reached adulthood relative to siblings that experienced low-aggression environments, even though all bees were kept in a common environment during adulthood. Unlike other animals including humans however, high-aggression honey bees were more, rather than less, resilient to immune challenge, assessed as neonicotinoid pesticide susceptibility. Moreover, aggression was negatively correlated with ectoparasitic mite presence. In honey bees, early-life social experience has broad effects, but increased aggression is decoupled from negative health outcomes. Because honey bees and humans share aspects of their physiological response to aggressive social encounters, our findings represent a step towards identifying ways to improve individual resiliency. Pre-adult social experience may be crucial to the health of the ecologically threatened honey bee.

The power and promise of applying genomics to honey bee health

New genomic tools and resources are now being used to both understand honey bee health and develop tools to better manage it. Here, we describe the use of genomic approaches to identify and characterize bee parasites and pathogens, examine interactions among these parasites and pathogens, between them and their bee hosts, and to identify genetic markers for improved breeding of more resilient bee stocks. We also discuss several new genomic techniques that can be used to more efficiently study, monitor and improve bee health. In the case of using RNAi-based technologies to mitigate diseases in bee populations, we highlight advantages, disadvantages and strategies to reduce risk. The increased use of genomic analytical tools and manipulative technologies has already led to significant advances, and holds great promise for improvements in the health of honey bees and other critical pollinator species.

Comparative brain transcriptomic analyses of scouting across distinct behavioural and ecological contexts in honeybees

Individual differences in behaviour are often consistent across time and contexts, but it is not clear whether such consistency is reflected at the molecular level. We explored this issue by studying scouting in honeybees in two different behavioural and ecological contexts: finding new sources of floral food resources and finding a new nest site. Brain gene expression profiles in food-source and nest-site scouts showed a significant overlap, despite large expression differences associated with the two different contexts. Class prediction and 'leave-one-out' cross-validation analyses revealed that a bee's role as a scout in either context could be predicted with 92.5% success using 89 genes at minimum. We also found that genes related to four neurotransmitter systems were part of a shared brain molecular signature in both types of scouts, and the two types of scouts were more similar for genes related to glutamate and GABA than catecholamine or acetylcholine signalling. These results indicate that consistent behavioural tendencies across different ecological contexts involve a mixture of similarities and differences in brain gene expression.

Big Data: Astronomical or Genomical?

Genomics is a Big Data science and is going to get much bigger, very soon, but it is not known whether the needs of genomics will exceed other Big Data domains. Projecting to the year 2025, we compared genomics with three other major generators of Big Data: astronomy, YouTube, and Twitter. Our estimates show that genomics is a “four-headed beast”—it is either on par with or the most demanding of the domains analyzed here in terms of data acquisition, storage, distribution, and analysis. We discuss aspects of new technologies that will need to be developed to rise up and meet the computational challenges that genomics poses for the near future. Now is the time for concerted, community-wide planning for the “genomical” challenges of the next decade.

This perspective considers the growth of genomics over the next ten years and assesses the computational needs that we will face relative to other "Big Data" activities such as astronomy, YouTube, and Twitter.

Insights into the Transcriptional Architecture of Behavioral Plasticity in the Honey Bee Apis mellifera

Honey bee colonies exhibit an age-related division of labor, with worker bees performing discrete sets of behaviors throughout their lifespan. These behavioral states are associated with distinct brain transcriptomic states, yet little is known about the regulatory mechanisms governing them. We used CAGEscan (a variant of the Cap Analysis of Gene Expression technique) for the first time to characterize the promoter regions of differentially expressed brain genes during two behavioral states (brood care (aka “nursing”) and foraging) and identified transcription factors (TFs) that may govern their expression. More than half of the differentially expressed TFs were associated with motifs enriched in the promoter regions of differentially expressed genes (DEGs), suggesting they are regulators of behavioral state. Strikingly, five TFs (nf-kb, egr, pax6, hairy, and clockwork orange) were predicted to co-regulate nearly half of the genes that were upregulated in foragers. Finally, differences in alternative TSS usage between nurses and foragers were detected upstream of 646 genes, whose functional analysis revealed enrichment for Gene Ontology terms associated with neural function and plasticity. This demonstrates for the first time that alternative TSSs are associated with stable differences in behavior, suggesting they may play a role in organizing behavioral state.

Voxel-based analysis of the immediate early gene, c-jun, in the honey bee brain after a sucrose stimulus

Immediate early genes (IEGs) have served as useful markers of brain neuronal activity in mammals, and more recently in insects. The mammalian canonical IEG, c-jun, is part of regulatory pathways conserved in insects and has been shown to be responsive to alarm pheromone in honey bees. We tested whether c-jun was responsive in honey bees to another behaviourally relevant stimulus, sucrose, in order to further identify the brain regions involved in sucrose processing. To identify responsive regions, we developed a new method of voxel-based analysis of c-jun mRNA expression. We found that c-jun is expressed in somata throughout the brain. It was rapidly induced in response to sucrose stimuli, and it responded in somata near the antennal and mechanosensory motor centre, mushroom body calices and lateral protocerebrum, which are known to be involved in sucrose processing. c-jun also responded to sucrose in somata near the lateral suboesophageal ganglion, dorsal optic lobe, ventral optic lobe and dorsal posterior protocerebrum, which had not been previously identified by other methods. These results demonstrate the utility of voxel-based analysis of mRNA expression in the insect brain.

Social evolution. Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.