In many organisms, interactions among genes lead to multiple functional states, and changes to interactions can lead to transitions into new states. These transitions can be related to bifurcations (or critical points) in dynamical systems theory. Characterizing these collective transitions is a major challenge for systems biology. Here, we develop a statistical method for identifying bistability near a continuous transition directly from high-dimensional gene expression data. We apply the method to data from honey bees, where a known developmental transition occurs between bees performing tasks in the nest and leaving the nest to forage. Our method, which makes use of the expected shape of the distribution of gene expression levels near a transition, successfully identifies the emergence of bistability and links it to genes that are known to be involved in the behavioral transition. This proof of concept demonstrates that going beyond correlative analysis to infer the shape of gene expression distributions might be used more generally to identify collective transitions from gene expression data.

Variation in endocrine signaling is proposed to underlie the evolution and regulation of social life histories, but the genetic architecture of endocrine signaling is still poorly understood. An excellent example of a hormonally influenced set of social traits is found in the honey bee (Apis mellifera): a dynamic and mutually suppressive relationship between juvenile hormone (JH) and the yolk precursor protein vitellogenin (Vg) regulates behavioral maturation and foraging of workers. Several other traits cosegregate with these behavioral phenotypes, comprising the pollen hoarding syndrome (PHS) one of the best-described animal behavioral syndromes. Genotype differences in responsiveness of JH to Vg are a potential mechanistic basis for the PHS. Here, we reduced Vg expression via RNA interference in progeny from a backcross between 2 selected lines of honey bees that differ in JH responsiveness to Vg reduction and measured JH response and ovary size, which represents another key aspect of the PHS. Genetic mapping based on restriction site-associated DNA tag sequencing identified suggestive quantitative trait loci (QTL) for ovary size and JH responsiveness. We confirmed genetic effects on both traits near many QTL that had been identified previously for their effect on various PHS traits. Thus, our results support a role for endocrine control of complex traits at a genetic level. Furthermore, this first example of a genetic map of a hormonal response to gene knockdown in a social insect helps to refine the genetic understanding of complex behaviors and the physiology that may underlie behavioral control in general.

Foraging exposes organisms to rewarding and aversive events, providing a selective advantage for maximizing the former while minimizing the latter. Honey bees (Apis mellifera) associate environmental stimuli with appetitive or aversive experiences, forming preferences for scents, locations, and visual cues. Preference formation is influenced by inter-individual variation in sensitivity to rewarding and aversive stimuli, which can be modulated by pharmacological manipulation of biogenic amines. We propose that foraging experiences act on biogenic amine pathways to induce enduring changes to stimulus responsiveness. To simulate varied foraging conditions, freely-moving bees were housed in cages where feeders offered combinations of sucrose solution, floral scents, and aversive electric shock. Transient effects were excluded by providing bees with neutral conditions for three days prior to all subsequent assays. Sucrose responsiveness was reduced in bees that had foraged for scented rather than unscented sucrose under benign conditions. This was not the case under aversive foraging conditions, suggesting an adaptive tuning process which maximizes preference for high quality, non-aversive floral sites. Foraging conditions also influenced antennal lobe octopamine and serotonin, neuromodulators involved in stimulus responsiveness and foraging site evaluation. Our results suggest that individuals' foraging experiences durably modify neurochemistry and shape future foraging behaviour.

Iron granules containing superparamagnetic magnetite act as magnetoreceptor for magnetoreception in honey bees. Biomineralization of iron granules occurs in the iron deposition vesicles of trophocytes and requires the participation of actin, myosin, ferritin2, and ATP synthase. The mechanism of magnetoreception in honey bees can be explored by suppressing the formation of iron granules. Toward this goal, we injected double-stranded RNA of ferritin2 and ferritin1 into newly emerged worker honey bees to knock down these genes via RNA interference. We confirmed that mRNA and protein production of the ferritins was inhibited, leading to immature iron granules. Downregulating ferritin2 and ferritin1, moreover, leads to different deposition morphology of 7.5-nm diameter iron particles, indicating that the two genes play different roles in the formation of iron granules in worker honey bees.

Honey bees (Apis mellifera) provide a system for studying social and food-related behavior. A caste of workers performs age-related tasks: young bees (nurses) usually feed the brood and other adult bees inside the nest, while older bees (foragers) forage outside for pollen, a protein/lipid source, or nectar, a carbohydrate source. The workers' transition from nursing to foraging and their foraging preferences correlate with differences in gustatory perception, metabolic gene expression, and endocrine physiology including the endocrine factors vitellogenin (Vg) and juvenile hormone (JH). However, the understanding of connections among social behavior, energy metabolism, and endocrine factors is incomplete. We used RNA interference (RNAi) to perturb the gene network of Vg and JH to learn more about these connections through effects on gustation, gene transcripts, and physiology. The RNAi perturbation was achieved by single and double knockdown of the genes ultraspiracle (usp) and vg, which encode a putative JH receptor and Vg, respectively. The double knockdown enhanced gustatory perception and elevated hemolymph glucose, trehalose, and JH. We also observed transcriptional responses in insulin like peptide 1 (ilp1), the adipokinetic hormone receptor (AKHR), and cGMP-dependent protein kinase (PKG, or "foraging gene" Amfor). Our study demonstrates that the Vg-JH regulatory module controls changes in carbohydrate metabolism, but not lipid metabolism, when worker bees shift from nursing to foraging. The module is also placed upstream of ilp1, AKHR, and PKG for the first time. As insulin, adipokinetic hormone (AKH), and PKG pathways influence metabolism and gustation in many animals, we propose that honey bees have conserved pathways in carbohydrate metabolism and conserved connections between energy metabolism and gustatory perception. Thus, perhaps the bee can make general contributions to the understanding of food-related behavior and metabolic disorders.

Various bioactive food compounds may confer health and longevity benefits, possibly through altering or preserving the epigenome. While bioactive food compounds are widely being marketed for human consumption as 'improving health and longevity' by counteracting harmful effects of poor nutrition and lifestyle, claimed effects are often not adequately documented. Using the honey bee (Apis mellifera) as a model species, we here employed a multi-step screening approach to investigate seven compounds for effects on lifespan and DNA methylation using ELISA and whole genome bisulfite sequencing (WGBS). A positive longevity effect was detected for valproic acid, isovaleric acid, and cyanocobalamin. For curcumin, we found that lifespan shortening caused by ethanol intake, was restored when curcumin and ethanol were co-administered. Furthermore, we identified region specific DNA methylation changes as a result of ethanol intake. Ethanol specific changes in DNA methylation were fully or partially blocked in honey bees receiving ethanol and curcumin together. Ethanol-affected and curcumin-blocked differentially methylated regions covered genes involved in fertility, temperature regulation and tubulin transport. Our results demonstrate fundamental negative effects of low dose ethanol consumption on lifespan and associated DNA methylation changes and present a proof-of-principle on how longevity and DNA methylation changes can be negated by the bioactive food component curcumin. Our findings provide a fundament for further studies of curcumin in invertebrates.

Loss of brain function is one of the most negative and feared aspects of aging. Studies of invertebrates have taught us much about the physiology of aging and how this progression may be slowed. Yet, how aging affects complex brain functions, e.g., the ability to acquire new memory when previous experience is no longer valid, is an almost exclusive question of studies in humans and mammalian models. In these systems, age related cognitive disorders are assessed through composite paradigms that test different performance tasks in the same individual. Such studies could demonstrate that afflicted individuals show the loss of several and often-diverse memory faculties, and that performance usually varies more between aged individuals, as compared to conspecifics from younger groups. No comparable composite surveying approaches are established yet for invertebrate models in aging research. Here we test whether an insect can share patterns of decline similar to those that are commonly observed during mammalian brain aging. Using honey bees, we combine restrained learning with free-flight assays. We demonstrate that reduced olfactory learning performance correlates with a reduced ability to extinguish the spatial memory of an abandoned nest location (spatial memory extinction). Adding to this, we show that learning performance is more variable in old honey bees. Taken together, our findings point to generic features of brain aging and provide the prerequisites to model individual aspects of learning dysfunction with insect models.

The reproductive ground plan hypothesis of social evolution suggests that reproductive controls of a solitary ancestor have been co-opted during social evolution, facilitating the division of labor among social insect workers. Despite substantial empirical support, the generality of this hypothesis is not universally accepted. Thus, we investigated the prediction of particular genes with pleiotropic effects on ovarian traits and social behavior in worker honey bees as a stringent test of the reproductive ground plan hypothesis. We complemented these tests with a comprehensive genome scan for additional quantitative trait loci (QTL) to gain a better understanding of the genetic architecture of the ovary size of honey bee workers, a morphological trait that is significant for understanding social insect caste evolution and general insect biology.

The effect of larval nutrition on female fertility in honey bees is a focus for both scientific studies and for practical applications in beekeeping. In general, morphological traits are standards for classifying queens and workers and for evaluating their quality. In recent years, in vitro rearing techniques have been improved and used in many studies; they can produce queen-like and worker-like bees. Here, we questioned whether queens and workers reared in vitro are the same as queens and workers reared in a natural hive environment. We reared workers and queens both in vitro and naturally in beehives to test how these different environments affect metabolic physiology and candidate genes in newly emerged queens and workers. We found that sugar (glucose and trehalose) levels differed between queens and workers in both in vitro and in-hive-reared bees. The in vitro-reared bees had significantly higher levels of lipids in the abdomen. Moreover, hive reared queens had almost 20 times higher levels of vitellogenin than in vitro-reared queens, despite similar morphologies. In addition, hive-reared bees had significantly higher levels of expression of mrjp1 In conclusion, in vitro rearing produces queens and workers that differ from those reared in the hive environment at physiological and gene expression levels.This article has an associated First Person interview with the first author of the paper.

Eusocial insects divide their labour so that individuals working inside the nest are affected by external conditions through a cascade of social interactions. Honey bees (Apis mellifera) transfer food and information via mouth-to-mouth social feeding, ie trophallaxis, a process known to be modulated by the rate of food flow at feeders and familiarity of food's scent. Little is understood about how aversive foraging conditions such as predation and con-specific competition affect trophallaxis. We hypothesized that aversive conditions have an impact on food transfer inside the colony. Here we explore the effect of foragers' aversive experience on downstream trophallaxis in a cage paradigm. Each cage contained one group of bees that was separated from feeders by mesh and allowed to feed only through trophallaxis, and another group that had access to feeders and self-specialized to either forage or distribute food. Our results show that aversive foraging conditions increase non-foragers' trophallaxis with bees restricted from feeder access when food is scented, and have the opposite effect when food is unscented. We discuss potential behavioural mechanisms and implications for the impact of aversive conditions such as malaise inducing toxins, predation, and con-specific competition.

Food choice and eating behavior affect health and longevity. Large-scale research efforts aim to understand the molecular and social/behavioral mechanisms of energy homeostasis, body weight, and food intake. Honey bees (Apis mellifera) could provide a model for these studies since individuals vary in food-related behavior and social factors can be controlled. Here, we examine a potential role of peripheral insulin receptor substrate (IRS) expression in honey bee foraging behavior. IRS is central to cellular nutrient sensing through transduction of insulin/insulin-like signals (IIS). By reducing peripheral IRS gene expression and IRS protein amount with the use of RNA interference (RNAi), we demonstrate that IRS influences foraging choice in two standard strains selected for different food-hoarding behavior. Compared with controls, IRS knockdowns bias their foraging effort toward protein (pollen) rather than toward carbohydrate (nectar) sources. Through control experiments, we establish that IRS does not influence the bees' sucrose sensory response, a modality that is generally associated with food-related behavior and specifically correlated with the foraging preference of honey bees. These results reveal a new affector pathway of honey bee social foraging, and suggest that IRS expressed in peripheral tissue can modulate an insect's foraging choice between protein and carbohydrate sources.

Our interest in healthy aging and in evolutionarily conserved mechanisms of lifespan extension prompted us to investigate whether features of age-related decline in the honey bee could be attenuated with resveratrol. Resveratrol is regarded as a caloric restriction mimetic known to extend lifespan in some but not all model species. The current, prevailing view is that resveratrol works largely by activating signaling pathways. It has also been suggested that resveratrol may act as an antioxidant and confer protection against nervous system impairment and oxidative stress. To test whether honey bee lifespan, learning performance, and food perception could be altered by resveratrol, we supplemented the diets of honey bees and measured lifespan, olfactory learning, and gustatory responsiveness to sucrose. Furthermore, to test the effects of resveratrol under metabolic challenge, we used hyperoxic environments to generate oxidative stress. Under normal oxygen conditions, two resveratrol treatments-30 and 130 μM-lengthened average lifespan in wild-type honey bees by 38% and 33%, respectively. Both resveratrol treatments also lengthened maximum and median lifespan. In contrast, hyperoxic stress abolished the resveratrol life-extension response. Furthermore, resveratrol did not affect learning performance, but did alter gustation. Honey bees that were not fed resveratrol exhibited greater responsiveness to sugar, while those supplemented with resveratrol were less responsive to sugar. We also discovered that individuals fed a high dose of resveratrol-compared to controls-ingested fewer quantities of food under ad libitum feeding conditions.

Vitellogenin is a yolk precursor protein in most oviparous females. In the advanced eusocial honeybee, Apis mellifera (Hymenoptera: Apidae), vitellogenin has recently attracted much interest as this protein, in addition to a classical function in oocyte development in the reproductive queen caste, has evolved functions in the facultatively sterile female worker caste not documented in other species. However, research on the spatial dynamics of vitellogenin in various tissues is not easily performed with available tools. Here we present an immunogold staining procedure that visualizes honeybee vitellogenin in resin embedded tissue. To establish the protocol, we used ovaries of worker bees from colonies with and without a queen. Under the first condition, vitellogenin is assumed not to be present in the workers' ovaries. Under the second condition, the ovaries of worker bees become vitellogenic, with abundant opportunities for detection of complex patterns of vitellogenin uptake and storage. By use of this experimental setup, the staining method is shown to be both sensitive and specific. To demonstrate the functional significance of the protocol, it was subsequently used to identify vitellogenin protein in the hypopharyngeal glands (brood food producing head glands) of nursing worker bees and in adjacent head fat body cells for the first time. Localization of vitellogenin in these tissues supports previously hypothesized roles of vitellogenin in social behavior. This protocol thus provides deeper insights into the functions of vitellogenin in the honeybee.

The protein vitellogenin (Vg) plays a central role in lipid transportation in most egg-laying animals. High Vg levels correlate with stress resistance and lifespan potential in honey bees (Apis mellifera). Vg is the primary circulating zinc-carrying protein in honey bees. Zinc is an essential metal ion in numerous biological processes, including the function and structure of many proteins. Measurements of Zn2+ suggest a variable number of ions per Vg molecule in different animal species, but the molecular implications of zinc-binding by this protein are not well-understood. We used inductively coupled plasma mass spectrometry to determine that, on average, each honey bee Vg molecule binds 3 Zn2+ -ions. Our full-length protein structure and sequence analysis revealed seven potential zinc-binding sites. These are located in the β-barrel and α-helical subdomains of the N-terminal domain, the lipid binding site, and the cysteine-rich C-terminal region of unknown function. Interestingly, two potential zinc-binding sites in the β-barrel can support a proposed role for this structure in DNA-binding. Overall, our findings suggest that honey bee Vg bind zinc at several functional regions, indicating that Zn2+ -ions are important for many of the activities of this protein. In addition to being potentially relevant for other egg-laying species, these insights provide a platform for studies of metal ions in bee health, which is of global interest due to recent declines in pollinator numbers.

The honey bee has been extensively studied as a model for neuronal circuit and memory function and more recently has emerged as an unconventional model in biogerontology. Yet, the detailed knowledge of neuronal processing in the honey bee brain contrasts with the very sparse information available on glial cells. In other systems glial cells are involved in nutritional homeostasis, detoxification, and aging. These glial functions have been linked to metabolic enzymes, such as glutamine synthetase and glycogen phosphorylase. As a step in identifying functional roles and potential differences among honey bee glial types, we examined the spatial distribution of these enzymes and asked if enzyme abundance is associated with aging and other processes essential for survival. Using immunohistochemistry and confocal laser microscopy we demonstrate that glutamine synthetase and glycogen phosphorylase are abundant in glia but appear to co-localize with different glial sub-types. The overall spatial distribution of both enzymes was not homogenous and differed markedly between different neuropiles and also within each neuropil. Using semi-quantitative Western blotting we found that rapid aging, typically observed in shortest-lived worker bees (foragers), was associated with declining enzyme levels. Further, we found enzyme abundance changes after severe starvation stress, and that glutamine synthetase is associated with food response. Together, our data indicate that aging and nutritional physiology in bees are linked to glial specific metabolic enzymes. Enzyme specific localization patterns suggest a functional differentiation among identified glial types.

Vitellogenin (Vg) has been implicated as a central protein in the immunity of egg-laying animals. Studies on a diverse set of species suggest that Vg supports health and longevity through binding to pathogens. Specific studies of honey bees (Apis mellifera) further indicate that the vitellogenin (vg) gene undergoes selection driven by local pathogen pressures. Determining the complete 3D structure of full-length Vg (flVg) protein will provide insights regarding the structure-function relationships underlying allelic variation. Honey bee Vg has been described in terms of function, and two subdomains have been structurally described, while information about the other domains is lacking. Here, we present a structure prediction, restrained by experimental data, of flVg from honey bees. To achieve this, we performed homology modeling and used AlphaFold before using a negative-stain electron microscopy map to restrict, orient, and validate our 3D model. Our approach identified a highly conserved Ca2+ -ion-binding site in a von Willebrand factor domain that might be central to Vg function. Thereafter, we used rigid-body fitting to predict the relative position of high-resolution domains in a flVg model. This mapping represents the first experimentally validated full-length protein model of a Vg protein and is thus relevant for understanding Vg in numerous species. Our results are also specifically relevant to honey bee health, which is a topic of global concern due to rapidly declining pollinator numbers.

Honey bees (Apis mellifera) provide an excellent model for studying how complex social behavior evolves and is regulated. Social behavioral traits such as the division of labor have been mapped to specific genomic regions in quantitative trait locus (QTL) studies. However, relating genomic mapping to gene function and regulatory mechanism remains a big challenge for geneticists. In honey bee workers, division of labor is known to be regulated by reproductive physiology, but the genetic basis of this regulation remains unknown. In this case, QTL studies have identified tyramine receptor 1 (TYR1) as a candidate gene in region pln2, which is associated with multiple worker social traits and reproductive anatomy. Tyramine (TA), a neurotransmitter, regulates physiology and behavior in diverse insect species including honey bees. Here, we examine directly the effects of TYR1 and TA on worker reproductive physiology, including ovariole number, ovary function and the production of vitellogenin (VG, an egg yolk precursor). First, we used a pharmacology approach to demonstrate that TA affects ovariole number during worker larval development and increases ovary maturation during the adult stage. Second, we used a gene knockdown approach to show that TYR1 regulates vg transcription in adult workers. Finally, we estimated correlations in gene expression and propose that TYR1 may regulate vg transcription by coordinating hormonal and nutritional signals. Taken together, our results suggest TYR1 and TA play important roles in regulating worker reproductive physiology, which in turn regulates social behavior. Our study exemplifies a successful forward-genetic strategy going from QTL mapping to gene function.

Many animals are characterized by declining brain function at advanced ages, including honeybees (Apis mellifera). Variation in honeybee social development, moreover, results in individual differences in the progression of aging that may be accelerated, delayed, and sometimes reversed by changes in behavior. Here, we combine manipulations of social development with a measurement of sensory sensitivity, Pavlovian (associative) learning, and a proteomic technique to study the brain of aged honeybees. First, we confirm that sensory sensitivity can remain intact during aging, and that age-associated learning deficits are specific to bees that forage, a behavior typically expressed after a period of nursing activity. These initial data go beyond previous findings by showing how foragers age in social groups of different age compositions and sizes. Thereafter, we establish that learning ability can recover in aged foragers that revert to nursing tasks. Finally, we use liquid chromatography coupled to tandem mass spectrometry (LC-MS(2)) to describe proteomic differences between central brains, from reverted former foragers that varied in recovery of learning performance, and from nurse bees that varied in learning ability but never foraged. We find that recovery is positively associated with levels of stress response/cellular maintenance proteins in the central brain, while variation in learning before aging is negatively associated with the amounts of metabolic enzymes in the brain tissue. Our work provides the strongest evidence, thus far, for reversibility of learning deficits in aged honeybees, and indicates that recovery-related brain plasticity is connected to cellular stress resilience, maintenance and repair processes.

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.

Temporal division of labor and foraging specialization are key characteristics of honeybee social organization. Worker honeybees (Apis mellifera) initiate foraging for food around their third week of life and often specialize in collecting pollen or nectar before they die. Variation in these fundamental social traits correlates with variation in worker reproductive physiology. However, the genetic and hormonal mechanisms that mediate the control of social organization are not understood and remain a central question in social insect biology. Here we demonstrate that a yolk precursor gene, vitellogenin, affects a complex suite of social traits. Vitellogenin is a major reproductive protein in insects in general and a proposed endocrine factor in honeybees. We show by use of RNA interference (RNAi) that vitellogenin gene activity paces onset of foraging behavior, primes bees for specialized foraging tasks, and influences worker longevity. These findings support the view that the worker specializations that characterize hymenopteran sociality evolved through co-option of reproductive regulatory pathways. Further, they demonstrate for the first time how coordinated control of multiple social life-history traits can originate via the pleiotropic effects of a single gene that affects multiple physiological processes.