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.

Developmental Transcriptome for a Facultatively Eusocial Bee, Megalopta genalis

Transcriptomes provide excellent foundational resources for mechanistic and evolutionary analyses of complex traits. We present a developmental transcriptome for the facultatively eusocial bee Megalopta genalis, which represents a potential transition point in the evolution of eusociality. A de novo transcriptome assembly of Megalopta genalis was generated using paired-end Illumina sequencing and the Trinity assembler. Males and females of all life stages were aligned to this transcriptome for analysis of gene expression profiles throughout development. Gene Ontology analysis indicates that stage-specific genes are involved in ion transport, cell-cell signaling, and metabolism. A number of distinct biological processes are upregulated in each life stage, and transitions between life stages involve shifts in dominant functional processes, including shifts from transcriptional regulation in embryos to metabolism in larvae, and increased lipid metabolism in adults. We expect that this transcriptome will provide a useful resource for future analyses to better understand the molecular basis of the evolution of eusociality and, more generally, phenotypic plasticity.

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.

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.

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.