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

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.

Beyond conflict: Kinship theory of intragenomic conflict predicts individual variation in altruistic behaviour

Behavioural variation is essential for animals to adapt to different social and environmental conditions. The Kinship Theory of Intragenomic Conflict (KTIC) predicts that parent-specific alleles can support different behavioural strategies to maximize allele fitness. Previous studies, including in honey bees (Apis mellifera), supported predictions of the KTIC for parent-specific alleles to promote selfish behaviour. Here, we test the KTIC prediction that for altruism-promoting genes (i.e. those that promote behaviours that support the reproductive fitness of kin), the allele with the higher altruism optimum should be selected to be expressed while the other is silenced. In honey bee colonies, workers act altruistically when tending to the queen by performing a 'retinue' behaviour, distributing the queen's mandibular pheromone (QMP) throughout the hive. Workers exposed to QMP do not activate their ovaries, ensuring they care for the queen's brood instead of competing to lay unfertilized eggs. Due to the haplodiploid genetics of honey bees, the KTIC predicts that response to QMP is favoured by the maternal genome. We report evidence for parent-of-origin effects on the retinue response behaviour, ovarian development and gene expression in brains of worker honey bees exposed to QMP, consistent with the KTIC. Additionally, we show enrichment for genes with parent-of-origin expression bias within gene regulatory networks associated with variation in bees' response to QMP. Our study demonstrates that intragenomic conflict can shape diverse social behaviours and influence expression patterns of single genes as well as gene networks.

Histone acetyltransferases and external demands influence task switching in Temnothorax ants

In social hymenopterans, workers specialize in different tasks. Whether a worker nurses the brood or forages is influenced by the responsiveness for task-related cues which in turn is determined by gene expression. Task choice is dynamic and changes throughout a worker's life, e.g. with age or in response to increased demands for certain tasks. Behavioural switches require the ability to adjust gene expression but the mechanisms regulating such transcriptional adaptations remain elusive. We investigated the role of histone acetylation in task specialization and behavioural flexibility in Temnothorax longispinosus ants. By inhibiting p300/CBP histone acetyltransferases (HAT) and manipulating colony composition, we found that HAT inhibition impairs the ability of older workers to switch to brood care. Yet, HAT inhibition increased the ability of young workers to accelerate their behavioural development and switch to foraging. Our data suggest that HAT in combination with social signals indicating task demands play an important role in modulating behaviour. Elevated HAT activity may contribute to keeping young brood carers from leaving the nest, where they would be exposed to high mortality. These findings shed light on the epigenetic processes underlying behavioural flexibility in animals and provide insight into the mechanisms of task specialization in social insects.

A comparative transcriptome analysis of the head of 1 and 9 days old worker honeybees (Apis mellifera)

The role of bees in the environment, economic, biodiversity and pharmaceutical industries is due to its social behavior, which is oriented from the brain and hypopharyngeal gland that is the center of royal jelly (RJ) production. Limited studies have been performed on the head gene expression profile at the RJ production stage. The aim of this study was to compare the gene expressions in 9 and 1-day-old (DO) honeybee workers in order to achieve better understanding about head gene expression pattern. After sequencing of RNAs, transcriptome and their networks were compared. The head expression profile undergoes various changes. 1662 gene transcripts had differential expressions which 1125 and 537 were up and down regulated, respectively, in 9_DO compared with 1_DO honey bees. The day 1th had more significant role in the expression of genes related to RJ production as major RJ protein 1, 2, 3, 5, 6 and 9 encoding genes, but their maximum secretion occurred at day 9th. All process related to hypopharyngeal glands activities as CYP450 gene, fatty acid synthase gene, vitamin B6 metabolism and some of genes involved in fatty acid elongation and degradation process had an upward trend from 1_DO and were age-dependent. By increasing the age, the activity of pathways related to immune system increased for keeping the health of bees against the chemical compound. The expression of aromatic amino acid genes involved in Phenylalanine, tyrosine and tryptophan biosynthesis pathway are essential for early stage of life. In 9_DO honeybees, the energy supplying, reducing stress, protein production and export pathways have a crucial role for support the body development and the social duties. It can be stated that the activity of honeybee head is focused on energy supply instead of storage, while actively trying to improve the level of cell dynamics for increasing the immunity and reducing stress. Results of current study identified key genes of certain behaviors of honeybee workers. Deeper considering of some pathways will be evaluated in future studies.

Evolutionary Invariant of the Structure of DNA Double Helix in RNAP II Core Promoters

Eukaryotic and archaeal RNA polymerase II (POL II) machinery is highly conserved, regardless of the extreme changes in promoter sequences in different organisms. The goal of our work is to find the cause of this conservatism. The representative sets of aligned promoter sequences of fifteen organisms belonging to different evolutional stages were studied. Their textual profiles, as well as profiles of the indexes that characterize the secondary structure and the mechanical and physicochemical properties, were analyzed. The evolutionarily stable, extremely heterogeneous special secondary structure of POL II core promoters was revealed, which includes two singular regions—hexanucleotide “INR” around TSS and octanucleotide “TATA element” of about −28 bp upstream. Such structures may have developed at some stage of evolution. It turned out to be so well matched for the pre-initiation complex formation and the subsequent initiation of transcription for POL II machinery that in the course of evolution there were selected only those nucleotide sequences that were able to reproduce these structural properties. The individual features of specific sequences representing the singular region of the promoter of each gene can affect the kinetics of DNA-protein complex formation and facilitate strand separation in double-stranded DNA at the TSS position.

Gene expression and chromatin conformation differs between worker bees performing different tasks

BACKGROUND

Revealing the effect of transcriptomic regulation on behavioral differences is a fundamental goal in biology, but the relationship between gene regulatory networks and individual behavior differences remains largely unknown. Honey bees are considered as good models for studying the mechanisms underlying gene expression changes and behavioral differences since they exhibit strong and obvious differences in tasks between individuals. The cis-regulatory regions usually contain the binding sites of diverse transcription factor (TFs) influencing bee behavior. Thus, the identification of cis-regulatory elements in the brains across different behavioral states is important for understanding how genomic and transcriptomic variations affect different tasks in honeybees.

METHODS

In this study, we employed transcriptome and genome-wide chromatin accessibility assays to analyze brain tissues of honey bees in different behavioral states for examining the relationship between individual behavior differences and brain gene expression changes. We also used the obtained open chromatin regions to identify cis-motifs associating differentially expressed TFs and genes in order to reveal the transcriptional regulatory mechanism related to the different tasks.

RESULTS

We identified genetic regulatory modules regulating different tasks that contained key TFs (CTCF, Trl and schlank) associated with open chromatin regions enriched for DNA sequence motifs belonging to the family of the corresponding TFs. The most prominent transcriptomic changes, which correlated with chromatin accessibility modifications within their proximal promoter regions, occurred in nervous system development, and were associated with behavior switch.

CONCLUSIONS

Our results revealed the regulatory landscape among three behavioral stages in honeybees and identified interactive molecular networks regulating different tasks. These results provide a comprehensive insight into behavioral differences of honeybees, which offers reference for future study.

Social divergence: molecular pathways underlying castes and longevity in a facultatively eusocial small carpenter bee

Unravelling the evolutionary origins of eusocial life is a longstanding endeavour in the field of evolutionary-developmental biology. Descended from solitary ancestors, eusocial insects such as honeybees have evolved ontogenetic division of labour in which short-lived workers perform age-associated tasks, while a long-lived queen produces brood. It is hypothesized that (i) eusocial caste systems evolved through the co-option of deeply conserved genes and (ii) longevity may be tied to oxidative damage mitigation capacity. To date, however, these hypotheses have been examined primarily among only obligately eusocial corbiculate bees. We present brain transcriptomic data from a Japanese small carpenter bee, Ceratina japonica (Apidae: Xylocopinae), which demonstrates both solitary and eusocial nesting in sympatry and lives 2 or more years in the wild. Our dataset captures gene expression patterns underlying first- and second-year solitary females, queens and workers, providing an unprecedented opportunity to explore the molecular mechanisms underlying caste-antecedent phenotypes in a long-lived and facultatively eusocial bee. We find that C. japonica's queens and workers are underpinned by divergent gene regulatory pathways, involving many differentially expressed genes well-conserved among other primitively eusocial bee lineages. We also find support for oxidative damage reduction as a proximate mechanism of longevity in C. japonica.

Tissue-specific transcriptional patterns underlie seasonal phenotypes in honey bees (Apis mellifera)

Faced with adverse conditions, such as winter in temperate regions or hot and dry conditions in tropical regions, many insect species enter a state of diapause, a period of dormancy associated with a reduction or arrest of physical activity, development and reproduction. Changes in common physiological pathways underlie diapause phenotypes in different insect species. However, most transcriptomic studies of diapause have not simultaneously evaluated and compared expression patterns in different tissues. Honey bees (Apis mellifera) represent a unique model system to study the mechanisms underpinning diapause-related phenotypes. In winter, honey bees exhibit a classic diapause phenotype, with reduced metabolic activity, increased physiological nutritional resources and altered hormonal profiles. However, winter bees actively heat their colony by vibrating their wing muscles; thus, this tissue is not quiescent. Here, we evaluated the transcriptional profiles of flight muscle tissue and fat body tissue (involved in nutrient storage, metabolism and immune function) of winter bees. We also evaluated two behavioural phenotypes of summer bees: nurses, which exhibit high nutritional stores and low flight activity, and foragers, which exhibit low nutritional stores and high flight activity. We found winter bees and nurses have similar fat body transcriptional profiles, whereas winter bees and foragers have similar flight muscle transcriptional profiles. Additionally, differentially expressed genes were enriched in diapause-related gene ontology terms. Thus, honey bees exhibit tissue-specific transcriptional profiles associated with seasonal phenotypes, laying the groundwork for future studies evaluating the mechanisms, evolution and consequences of this tissue-specific regulation.

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.

Missing Nurse Bees-Early Transcriptomic Switch From Nurse Bee to Forager Induced by Sublethal Imidacloprid

The environmental residue/sublethal doses of neonicotinoid insecticides are believed to generate a negative impact on pollinators, including honey bees. Here we report our recent investigation on how imidacloprid, one of the major neonicotinoids, affects worker bees by profiling the transcriptomes of various ages of bees exposed to different doses of imidacloprid during the larval stage. The results show that imidacloprid treatments during the larval stage severely altered the gene expression profiles and may induce precocious foraging. Differential expression of foraging regulators was found in 14-day-old treated adults. A high transcriptome similarity between larvae-treated 14-day-old adults and 20-day-old controls was also observed, and the similarity was positively correlated with the dose of imidacloprid. One parts per billion (ppb) of imidacloprid was sufficient to generate a long-term impact on the bee's gene expression as severe as with 50 ppb imidacloprid. The disappearance of nurse bees may be driven not only by the hive member constitution but also by the neonicotinoid-induced precocious foraging behavior.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

BACKGROUND

Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.

RESULTS

The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.

CONCLUSIONS

The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

Functional genomics of parental care of insects

Parental care was likely the first step most lineages made towards sociality. However, the molecular mechanisms that generate parental care are not broadly characterized. Insects are important as an evolutionary independent group from classic models of parental care, such as, house mice. They provide an opportunity to test the generality of our understanding. With this review, I survey the functional genomics of parental care of insects, summarize several recent advances in the broader framework for studying and understanding parental care, and finish with suggested priorities for further research. Although there are too few studies to draw definitive conclusions, I argue that natural selection appears to be rewiring existing gene networks to produce parental care, that the epigenetic mechanisms influencing parental care are not well understood, and, as an interesting early consensus, that genes strongly associated with carer/offspring interactions appear biased towards proteins that are secreted. I summarize the studies that have functionally validate candidate genes and highlight the increasing need to perform this work. I finish with arguments for both conceptual and practical changes moving forward. I argue that future work can increase the use of predictive frameworks, broaden its definition of conservation of mechanism to gene networks rather than single genes, and increase the use of more established comparative methods. I further highlight the practical considerations of standardizing analyses and reporting, increasing the sampling of both carers and offspring, better characterizing gene regulatory networks, better characterizing taxonomically restricted genes and any consistent role they have underpinning parental care, and using factorial designs to disentangle the influence of multiple variables on the expression of parental care.

Immediate early gene kakusei potentially plays a role in the daily foraging of honey bees

kakusei is a non-coding RNA that is overexpressed in foraging bee brain. This study describes a possible role of the IEG kakusei during the daily foraging of honey bees. kakusei was found to be transiently upregulated within two hours during rewarded foraging. Interestingly, during unrewarded foraging the gene was also found to be up-regulated, but immediately lowered when food was not rewarded. Moreover, the kakusei overexpression was diminished within a very short time when the time schedule of feeding was changed. This indicates the potential role of kakusei on the motivation of learned reward foraging. These results provide evidence for a dynamic role of kakusei during for aging of bees, and eventually its possible involvement in learning and memory. Thus the kakusei gene could be used as search tool in finding distinct molecular pathways that mediate diverse behavioral components of foraging.

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.

Social modularity: conserved genes and regulatory elements underlie caste-antecedent behavioural states in an incipiently social bee

The evolutionary origins of advanced eusociality, one of the most complex forms of phenotypic plasticity in nature, have long been a focus within the field of sociobiology. Although eusocial insects are known to have evolved from solitary ancestors, sociogenomic research among incipiently social taxa has only recently provided empirical evidence supporting theories that modular regulation and deeply conserved genes may play important roles in both the evolutionary emergence and elaboration of insect sociality. There remains, however, a paucity of data to further test the biological reality of these and other evolutionary theories among taxa in the earliest stages of social evolution. Here, we present brain transcriptomic data from the incipiently social small carpenter bee, Ceratina calcarata, which captures patterns of cis-regulation and gene expression associated with female maturation, and underlying two well-defined behavioural states, foraging and guarding, concurrently demonstrated by mothers and daughters during early autumn. We find that an incipiently social nest environment may dramatically affect gene expression. We further reveal foraging and guarding behaviours to be putatively caste-antecedent states in C. calcarata, and offer strong empirical support for the operation of modular regulation, involving deeply conserved and differentially expressed genes in the expression of early social forms.

Honey bee (Apis mellifera) larval pheromones may regulate gene expression related to foraging task specialization

Background

Foraging behavior in honey bees (Apis mellifera) is a complex phenotype that is regulated by physiological state and social signals. How these factors are integrated at the molecular level to modulate foraging behavior has not been well characterized. The transition of worker bees from nursing to foraging behaviors is mediated by large-scale changes in brain gene expression, which are influenced by pheromones produced by the queen and larvae. Larval pheromones can also stimulate foragers to leave the colony to collect pollen. However, the mechanisms underpinning this rapid behavioral plasticity in foragers that specialize in collecting pollen over nectar, and how larval pheromones impact these different behavioral states, remains to be determined. Here, we investigated the patterns of gene expression related to rapid behavioral plasticity and task allocation among honey bee foragers exposed to two larval pheromones, brood pheromone (BP) and (E)-beta-ocimene (EBO). We hypothesized that both pheromones would alter expression of genes in the brain related to foraging and would differentially impact brain gene expression depending on foraging specialization.

Results

Combining data reduction, clustering, and network analysis methods, we found that foraging preference (nectar vs. pollen) and pheromone exposure are each associated with specific brain gene expression profiles. Furthermore, pheromone exposure has a strong transcriptional effect on genes that are preferentially expressed in nectar foragers. Representation factor analysis between our study and previous landmark honey bee transcriptome studies revealed significant overlaps for both pheromone communication and foraging task specialization.

Conclusions

Our results suggest that, as social signals, pheromones alter expression patterns of foraging-related genes in the bee’s brain to increase pollen foraging at both long and short time scales. These results provide new insights into how social signals and task specialization are potentially integrated at the molecular level, and highlights the possible role that brain gene expression may play in honey bee behavioral plasticity across time scales.

Electronic supplementary material

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

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.

lncRNA profile of Apis mellifera and its possible role in behavioural transition from nurses to foragers

BACKGROUND

The behavioural transition from nurses to foragers in honey bees is known to be affected by intrinsic and extrinsic factors, including colony demography, hormone levels, brain chemistry and structure, and gene expression in the brain. However, the molecular mechanism underlying this behavioural transition of honey bees is still obscure.

RESULTS

Through RNA sequencing, we performed a comprehensive analysis of lncRNAs and mRNAs in honey bee nurses and foragers. Nurses and foragers from both typical colonies and single-cohort colonies were used to prepare six libraries to generate 49 to 100 million clear reads per sample. We obtained 6863 novel lncRNAs, 1480 differentially expressed lncRNAs between nurses and foragers, and 9308 mRNAs. Consistent with previous studies, lncRNAs showed features distinct from mRNAs, such as shorter lengths, lower exon numbers, and lower expression levels compared to mRNAs. Bioinformatic analysis showed that differentially expressed genes were mostly involved in the regulation of sensory-related events, such as olfactory receptor activity and odorant binding, and enriched Wnt and FoxO signaling pathways. Moreover, we found that lncRNAs TCONS_00356023, TCONS_00357367, TCONS_00159909 and mRNAs dop1, Kr-h1 and HR38 may play important roles in behavioural transition in honey bees.

CONCLUSION

This study characterized the expression profile of lncRNAs in nurses and foragers and provided a framework for further study of the role of lncRNAs in honey bee behavioural transition.

Mass Spectrometric Quantification of Arousal Associated Neurochemical Changes in Single Honey Bee Brains and Brain Regions

Honey bee foragers show a strong diurnal rhythm of foraging activity, and such behavioral changes are likely under the control of specific neuromodulators. To identify and quantify neuromodulators involved in regulating rest and arousal in honey bees, we established a mass spectrometric method for quantifying 14 different neurochemicals and precursor molecules. We measured forager type and brain region specific differences in amine levels from individual honey bee brains and brain regions. The observed differences in amine levels between resting and aroused foragers resemble findings in other species indicating a conserved molecular mechanism by glutamate and GABA in regulating arousal. Subesophageal ganglion specific changes in the histaminergic system and global increases in aspartate during arousal suggest a possible role of histamine and aspartate in feeding and arousal, respectively. More aminergic systems were significantly affected due to arousal in nectar foragers than in pollen foragers, implying that forager phenotypes differ not only in their food preference but also in their neuromodulatory signaling systems (brain states). Finally, we found that neurotransmitter precursors were better at distinguishing brain states in the central brain, while their end products correlated with arousal associated changes in sensory regions like the optic and antennal lobes.

Gene expression is more strongly associated with behavioural specialization than with age or fertility in ant workers

The ecological success of social insects is based on division of labour, not only between queens and workers, but also among workers. Whether a worker tends the brood or forages is influenced by age, fertility and nutritional status, with brood carers being younger, more fecund and more corpulent. Here, we experimentally disentangle behavioural specialization from age and fertility in Temnothorax longispinosus ant workers and analyse how these parameters are linked to whole-body gene expression. A total of 3,644 genes were associated with behavioural specialization which is ten times more than associated with age and 50 times more than associated with fertility. Brood carers were characterized by an upregulation of three Vitellogenin (Vg) genes, one of which, Vg-like A, was the most differentially expressed gene that was recently shown experimentally to control the switch from brood to worker care. The expression of Conventional Vg was unlinked to behavioural specialization, age or fertility, which contrasts to studies on bees and some ants. Diversity in Vg/Vg-like copy number and expression bias across ants supports subfunctionalization of Vg genes and indicates that some regulatory mechanisms of division of labour diverged in different ant lineages. Simulations revealed that our experimental dissociation of co-varying factors reduced transcriptomic noise, suggesting that confounding factors could potentially explain inconsistencies across transcriptomic studies of behavioural specialization in ants. Thus, our study reveals that worker gene expression is mainly linked to the worker's function for the colony and provides novel insights into the evolution of sociality in ants.

Immediate early genes in social insects: a tool to identify brain regions involved in complex behaviors and molecular processes underlying neuroplasticity

Social insects show complex behaviors and master cognitive tasks. The underlying neuronal mechanisms, however, are in most cases only poorly understood due to challenges in monitoring brain activity in freely moving animals. Immediate early genes (IEGs) that get rapidly and transiently expressed following neuronal stimulation provide a powerful tool for detecting behavior-related neuronal activity in vertebrates. In social insects, like honey bees, and in insects in general, this approach is not yet routinely established, even though these genes are highly conserved. First studies revealed a vast potential of using IEGs as neuronal activity markers to analyze the localization, function, and plasticity of neuronal circuits underlying complex social behaviors. We summarize the current knowledge on IEGs in social insects and provide ideas for future research directions.

Recent Advances in Behavioral (Epi)Genetics in Eusocial Insects

Eusocial insects live in societies in which distinct family members serve specific roles in maintaining the colony and advancing the reproductive ability of a few select individuals. Given the genetic similarity of all colony members, the diversity of morphologies and behaviors is surprising. Social communication relies on pheromones and olfaction, as shown by mutants of orco, the universal odorant receptor coreceptor, and through electrophysiological analysis of neuronal responses to pheromones. Additionally, neurohormonal factors and epigenetic regulators play a key role in caste-specific behavior, such as foraging and caste switching. These studies start to allow an understanding of the molecular mechanisms underlying social behavior and provide a technological foundation for future studies of eusocial insects. In this review, we highlight recent findings in eusocial insects that advance our understanding of genetic and epigenetic regulations of social behavior and provide perspectives on future studies using cutting-edge technologies.

Egr-1: A Candidate Transcription Factor Involved in Molecular Processes Underlying Time-Memory

In honey bees, continuous foraging is accompanied by a sustained up-regulation of the immediate early gene Egr-1 (early growth response protein-1) and candidate downstream genes involved in learning and memory. Here, we present a series of feeder training experiments indicating that Egr-1 expression is highly correlated with the time and duration of training even in the absence of the food reward. Foragers that were trained to visit a feeder over the whole day and then collected on a day without food presentation showed Egr-1 up-regulation over the whole day with a peak expression around 14:00. When exposed to a time-restricted feeder presentation, either 2 h in the morning or 2 h in the evening, Egr-1 expression in the brain was up-regulated only during the hours of training. Foragers that visited a feeder in the morning as well as in the evening showed two peaks of Egr-1 expression. Finally, when we prevented time-trained foragers from leaving the colony using artificial rain, Egr-1 expression in the brains was still slightly but significantly up-regulated around the time of feeder training. In situ hybridization studies showed that active foraging and time-training induced Egr-1 up-regulation occurred in the same brain areas, preferentially the small Kenyon cells of the mushroom bodies and the antennal and optic lobes. Based on these findings we propose that foraging induced Egr-1 expression can get regulated by the circadian clock after time-training over several days and Egr-1 is a candidate transcription factor involved in molecular processes underlying time-memory.

Vitellogenin-like A-associated shifts in social cue responsiveness regulate behavioral task specialization in an ant

Division of labor and task specialization explain the success of human and insect societies. Social insect colonies are characterized by division of labor, with workers specializing in brood care early and foraging later in life. Theory posits that this task switching requires shifts in responsiveness to task-related cues, yet experimental evidence is weak. Here, we show that a Vitellogenin (Vg) ortholog identified in an RNAseq study on the ant T. longispinosus is involved in this process: using phylogenetic analyses of Vg and Vg-like genes, we firstly show that this candidate gene does not cluster with the intensively studied honey bee Vg but falls into a separate Vg-like A cluster. Secondly, an experimental knockdown of Vg-like A in the fat body caused a reduction in brood care and an increase in nestmate care in young ant workers. Nestmate care is normally exhibited by older workers. We demonstrate experimentally that this task switch is at least partly based on Vg-like A-associated shifts in responsiveness from brood to worker cues. We thus reveal a novel mechanism leading to early behavioral maturation via changes in social cue responsiveness mediated by Vg-like A and associated pathways, which proximately play a role in regulating division of labor.

NanoCAGE: A Method for the Analysis of Coding and Noncoding 5'-Capped Transcriptomes

Transcripts in all eukaryotes are characterized by the 5'-end specific cap structure in mRNAs. Cap Analysis Gene Expression or CAGE makes use of these caps to specifically obtain cDNA fragments from the 5'-end of RNA and sequences those at high throughput for transcript identification and genome-wide mapping of transcription start sites for coding and noncoding genes. Here, we provide an improved version of our nanoCAGE protocol that has been developed for preparing CAGE libraries from as little as 50 ng of total RNA within three standard working days. Key steps in library preparation have been improved over our previously published protocol to obtain libraries having a good 5'-end selection and a more equal size distribution for higher sequencing efficiency on Illumina MiSeq and HiSeq sequencers. We recommend nanoCAGE as the method of choice for transcriptome profiling projects even from limited amounts of RNA, and as the best approach for genome-wide mapping of transcription start sites within promoter regions.

Honey bee foraging induces upregulation of early growth response protein 1, hormone receptor 38 and candidate downstream genes of the ecdysteroid signalling pathway

In honey bees, continuous foraging at an artificial feeder induced a sustained upregulation of the immediate early genes early growth response protein 1 (Egr-1) and hormone receptor 38 (Hr38). This gene expression response was accompanied by an upregulation of several Egr-1 candidate downstream genes: ecdysone receptor (EcR), dopamine/ecdysteroid receptor (DopEcR), dopamine decarboxylase and dopamine receptor 2. Hr38, EcR and DopEcR are components of the ecdysteroid signalling pathway, which is highly probably involved in learning and memory processes in honey bees and other insects. Time-trained foragers still showed an upregulation of Egr-1 when the feeder was presented at an earlier time of the day, suggesting that the genomic response is more dependent on the food reward than training time. However, presentation of the feeder at the training time without food was still capable of inducing a transient increase in Egr-1 expression. Thus, learnt feeder cues, or even training time, probably affect Egr-1 expression. In contrast, whole brain Egr-1 expression changes did not differ between dancing and nondancing foragers. On the basis of our results we propose that food reward induced continuous foraging ultimately elicits a genomic response involving Egr-1 and Hr38 and their downstream genes. Furthermore this genomic response is highly probably involved in foraging-related learning and memory responses.

Mitochondrial capacity, oxidative damage and hypoxia gene expression are associated with age-related division of labor in honey bee (Apis mellifera L.) workers

During adult life, honey bee workers undergo a succession of behavioral states. Nurse bees perform tasks inside the nest, and when they are about 2-3 weeks old they initiate foraging. This switch is associated with alterations in diet, and with the levels of juvenile hormone and vitellogenin circulating in hemolymph. It is not clear whether this behavioral maturation involves major changes at the cellular level, such as mitochondrial activity and the redox environment in the head, thorax and abdomen. Using high-resolution respirometry, biochemical assays and RT-qPCR, we evaluated the association of these parameters with this behavioral change. We found that tissues from the head and abdomen of nurses have a higher oxidative phosphorylation capacity than those of foragers, while for the thorax we found the opposite situation. As higher mitochondrial activity tends to generate more H₂O₂, and H₂O₂ is known to stabilize HIF-1α, this would be expected to stimulate hypoxia signaling. The positive correlation that we observed between mitochondrial activity and hif-1α gene expression in abdomen and head tissue of nurses would be in line with this hypothesis. Higher expression of antioxidant enzyme genes was observed in foragers, which could explain their low levels of protein carbonylation. No alterations were seen in nitric oxide (NO) levels, suggesting that NO signaling is unlikely to be involved in behavioral maturation. We conclude that the behavioral change seen in honey bee workers is reflected in differential mitochondrial activities and redox parameters, and we consider that this can provide insights into the underlying aging process.

The Neuropeptide Corazonin Controls Social Behavior and Caste Identity in Ants

Social insects are emerging models to study how gene regulation affects behavior because their colonies comprise individuals with the same genomes but greatly different behavioral repertoires. To investigate the molecular mechanisms that activate distinct behaviors in different castes, we exploit a natural behavioral plasticity in Harpegnathos saltator, where adult workers can transition to a reproductive, queen-like state called gamergate. Analysis of brain transcriptomes during the transition reveals that corazonin, a neuropeptide homologous to the vertebrate gonadotropin-releasing hormone, is downregulated as workers become gamergates. Corazonin is also preferentially expressed in workers and/or foragers from other social insect species. Injection of corazonin in transitioning Harpegnathos individuals suppresses expression of vitellogenin in the brain and stimulates worker-like hunting behaviors, while inhibiting gamergate behaviors, such as dueling and egg deposition. We propose that corazonin is a central regulator of caste identity and behavior in social insects.

RNA-sequence analysis of gene expression from honeybees (Apis mellifera) infected with Nosema ceranae

Honeybees (Apis mellifera) are constantly subjected to many biotic stressors including parasites. This study examined honeybees infected with Nosema ceranae (N. ceranae). N. ceranae infection increases the bees energy requirements and may contribute to their decreased survival. RNA-seq was used to investigate gene expression at days 5, 10 and 15 Post Infection (P.I) with N. ceranae. The expression levels of genes, isoforms, alternative transcription start sites (TSS) and differential promoter usage revealed a complex pattern of transcriptional and post-transcriptional gene regulation suggesting that bees use a range of tactics to cope with the stress of N. ceranae infection. N. ceranae infection may cause reduced immune function in the bees by: (i)disturbing the host amino acids metabolism (ii) down-regulating expression of antimicrobial peptides (iii) down-regulation of cuticle coatings and (iv) down-regulation of odorant binding proteins.

Molecular Evolution of Insect Sociality: An Eco-Evo-Devo Perspective

The evolution of eusociality is a perennial issue in evolutionary biology, and genomic advances have fueled steadily growing interest in the genetic changes underlying social evolution. Along with a recent flurry of research on comparative and evolutionary genomics in different eusocial insect groups (bees, ants, wasps, and termites), several mechanistic explanations have emerged to describe the molecular evolution of eusociality from solitary behavior. These include solitary physiological ground plans, genetic toolkits of deeply conserved genes, evolutionary changes in protein-coding genes, cis regulation, and the structure of gene networks, epigenetics, and novel genes. Despite this proliferation of ideas, there has been little synthesis, even though these ideas are not mutually exclusive and may in fact be complementary. We review available data on molecular evolution of insect sociality and highlight key biotic and abiotic factors influencing social insect genomes. We then suggest both phylogenetic and ecological evolutionary developmental biology (eco-evo-devo) perspectives for a more synthetic view of molecular evolution in insect societies.

Involvement of phosphorylated Apis mellifera CREB in gating a honeybee's behavioral response to an external stimulus

The transcription factor cAMP-response element-binding protein (CREB) is involved in neuronal plasticity. Phosphorylation activates CREB and an increased level of phosphorylated CREB is regarded as an indicator of CREB-dependent transcriptional activation. In honeybees(Apis mellifera)we recently demonstrated a particular high abundance of the phosphorylated honeybee CREB homolog (pAmCREB) in the central brain and in a subpopulation of mushroom body neurons. We hypothesize that these high pAmCREB levels are related to learning and memory formation. Here, we tested this hypothesis by analyzing brain pAmCREB levels in classically conditioned bees and bees experiencing unpaired presentations of conditioned stimulus (CS) and unconditioned stimulus (US). We demonstrate that both behavioral protocols display differences in memory formation but do not alter the level of pAmCREB in bee brains directly after training. Nevertheless, we report that bees responding to the CS during unpaired stimulus presentations exhibit higher levels of pAmCREB than nonresponding bees. In addition, Trichostatin A, a histone deacetylase inhibitor that is thought to enhance histone acetylation by CREB-binding protein, increases the bees' CS responsiveness. We conclude that pAmCREB is involved in gating a bee's behavioral response driven by an external stimulus.

DNA base modifications in honey bee and fruit fly genomes suggest an active demethylation machinery with species- and tissue-specific turnover rates

Well-known epigenetic DNA modifications in mammals include the addition of a methyl group and a hydroxyl group to cytosine, resulting in 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) respectively. In contrast, the abundance and the functional implications of these modifications in invertebrate model organisms such as the honey bee (Apis mellifera) and the fruit fly (Drosophila melanogaster) are not well understood. Here we show that both adult honey bees and fruit flies contain 5mC and also 5hmC. Using a highly sensitive liquid chromatography/tandem mass spectrometry (LC/MS/MS) technique, we quantified 5mC and 5hmC in different tissues of adult honey bee worker castes and in adult fruit flies. A comparison of our data with reports from human and mouse shed light on notable differences in 5mC and 5hmC levels between tissues and species.

•

Reporting cytosine modifications in uncharacterized tissues, phenotypes and species.

•

Quantification of 5mC and 5hmC suggests species-specific roles and turnover.​

•

Low levels of 5hmC relative to 5mC and cytosine in honey bees compared to mammals.

•

Honey bee abdominal tissues are richer in5hmC than the brain.

•

We found a higher 5hmC to 5mC ratio in fruit flies as compared to the honey bee.

Genomic sources of phenotypic novelty in the evolution of eusociality in insects

Genomic resources are now available for closely related species that vary in social behavior, providing insight on the genomics of social evolution. Changes in the architecture of gene regulatory networks likely influence the evolutionary trajectory of social traits. Evolutionarily novel genes are likely important in the evolution of social diversity among insects, but it is unclear whether new genes played a driving role in the advent or elaboration of eusociality or if they were instead a result of other genomic features of eusociality. The worker phenotype appears to be the center of genetic novelty, but the mechanisms for this remain unresolved. Future studies are needed to understand how genetic novelty arises, becomes incorporated into existing gene regulatory networks, and the effects this has on the evolution of social traits in closely related social and solitary species.

HOCOMOCO: expansion and enhancement of the collection of transcription factor binding sites models

Models of transcription factor (TF) binding sites provide a basis for a wide spectrum of studies in regulatory genomics, from reconstruction of regulatory networks to functional annotation of transcripts and sequence variants. While TFs may recognize different sequence patterns in different conditions, it is pragmatic to have a single generic model for each particular TF as a baseline for practical applications. Here we present the expanded and enhanced version of HOCOMOCO (http://hocomoco.autosome.ru and http://www.cbrc.kaust.edu.sa/hocomoco10), the collection of models of DNA patterns, recognized by transcription factors. HOCOMOCO now provides position weight matrix (PWM) models for binding sites of 601 human TFs and, in addition, PWMs for 396 mouse TFs. Furthermore, we introduce the largest up to date collection of dinucleotide PWM models for 86 (52) human (mouse) TFs. The update is based on the analysis of massive ChIP-Seq and HT-SELEX datasets, with the validation of the resulting models on in vivo data. To facilitate a practical application, all HOCOMOCO models are linked to gene and protein databases (Entrez Gene, HGNC, UniProt) and accompanied by precomputed score thresholds. Finally, we provide command-line tools for PWM and diPWM threshold estimation and motif finding in nucleotide sequences.