Differential gene expression in the mandibular glands of queen and worker honeybees, Apis mellifera L.: implications for caste-selective aldehyde and fatty acid metabolism

Hydroxy Fatty Acid Synthesis-Related mRNA as the Biomarker for Detecting Mislabeling of Honey Entomological Origin

The authentication of the entomological origin of honey is a widespread concern, necessitating the prompt establishment of an effective approach for distinguishing between Apis cerana cerana honey (ACH) and Apis mellifera ligustica honey (AMH). Hydroxy fatty acids (HFAs) found in honey are bee-derived components synthesized by the mandibular glands of worker bees. We previously discovered significant variations in the hydroxy fatty acid composition between ACH and AMH, suggesting their potential as indicators for identifying the authenticity of the entomological origin of honey. Herein, we identified differentially expressed genes associated with HFA synthesis by conducting transcriptome sequencing of the mandibular glands of AC and AM honeybees. Subsequently, we proposed a method for the relative quantitative analysis of bee-derived RNA components using real-time fluorescence quantitative polymerase chain reaction, which was supplemented by multivariate statistical analysis to further discern differences in HFA synthesis-related mRNA between ACH and AMH. The results showed that the mRNAs of FAXDC2 (fatty acid hydroxylase domain-containing protein 2) and FAS (fatty acid synthase) may serve as indicators to discern the entomological origin of honey. This study presents two novel biomarkers for detecting mislabeling of the entomological origin in ACH and AMH based on variations in bee-derived components.

The diverging epigenomic landscapes of honeybee queens and workers revealed by multiomic sequencing

The role of the epigenome in phenotypic plasticity is unclear presently. Here we used a multiomics approach to explore the nature of the epigenome in developing honey bee (Apis mellifera) workers and queens. Our data clearly showed distinct queen and worker epigenomic landscapes during the developmental process. Differences in gene expression between workers and queens become more extensive and more layered during the process of development. Genes known to be important for caste differentiation were more likely to be regulated by multiple epigenomic systems than other differentially expressed genes. We confirmed the importance of two candidate genes for caste differentiation by using RNAi to manipulate the expression of two genes that differed in expression between workers and queens were regulated by multiple epigenomic systems. For both genes the RNAi manipulation resulted in a decrease in weight and fewer ovarioles of newly emerged queens compared to controls. Our data show that the distinct epigenomic landscapes of worker and queen bees differentiate during the course of larval development.

A Comparison of RNA Interference via Injection and Feeding in Honey Bees

RNA interference (RNAi) has been used successfully to reduce target gene expression and induce specific phenotypes in several species. It has proved useful as a tool to investigate gene function and has the potential to manage pest populations and reduce disease pathogens. However, it is not known whether different administration methods are equally effective at interfering with genes in bees. Therefore, we compared the effects of feeding and injection of small interfering RNA (siRNA) on the messenger RNA (mRNA) levels of alpha-aminoadipic semialdehyde dehydrogenase (ALDH7A1), 4-coumarate-CoA ligase (4CL), and heat shock protein 70 (HSP70). Both feeding and injection of siRNA successfully knocked down the gene but feeding required more siRNA than the injection. Our results suggest that both feeding and injection of siRNA effectively interfere with brain genes in bees. The appropriateness of each method would depend on the situation.

Non-targeted lipidomics and transcriptomics analysis reveal the molecular underpinnings of mandibular gland development in Apis mellifera ligustica

The mandibular gland is an important exocrine gland of worker bees, which mainly secretes fatty acids and pheromones. Lipids have important roles in energy storage, membrane structure stabilization, and signaling. However, molecular underpinnings of mandibular gland development and lipid remodeling at the different physiological stages of worker bees is still lacking. In this study, we used scanning and transmission electron microscopy to reveal the morphological changes in secretory cells, and liquid chromatography-mass spectrometry and RNA-seq to investigate the lipidome and gene transcripts during development. The morphology of secretory cells was flat in newly emerged workers, becoming vacuolated and turgid when they were activated in nurse bees and foragers. Transport vesicles became denser from newly emerged bees to 21-day worker bees. Concentrations of 10-HDA reached a maximum within 15d workers and changes in genes expression were consistent with 10-HDA content. Non-targeted lipidomics analysis of newly emerged, 6d, and 15d worker bees revealed that PC and TAG were the main lipids in mandibular gland, and lipids dramatically altered across developmental stages. TAG 54:4 was increased most strongly at 6d and 15d worker bees, meanwhile, the abundances of TAG 54:1 and TAG 54:2 were decreased sharply. Further, transcriptomics analysis showed that differentially expressed genes were significantly enriched in key nutrient metabolic pathways, particularly lipid metabolism, in 6d and 15d bees. This multi-omic perspective provides a unique resource and deeper insight into bee mandibular gland development and baseline data for further study of the mandibular gland in worker bees.

Local Variation in Recombination Rates of the Honey Bee (Apis mellifera) Genome among Samples from Six Disparate Populations

Meiotic recombination is an essential component of eukaryotic sexual reproduction but its frequency varies within and between genomes. Although it is well-established that honey bees have a high recombination rate with about 20 cM/Mbp, the proximate and ultimate causes of this exceptional rate are poorly understood. Here, we describe six linkage maps of the Western Honey Bee Apis mellifera that were produced with consistent methodology from samples from distinct parts of the species' near global distribution. We compared the genome-wide rates and distribution of meiotic crossovers among the six maps and found considerable differences. Overall similarity of local recombination rates among our samples was unrelated to geographic or phylogenetic distance of the populations that our samples were derived from. However, the limited sampling constrains the interpretation of our results because it is unclear how representative these samples are. In contrast to previous studies, we found only in two datasets a significant relation between local recombination rate and GC content. Focusing on regions of particularly increased or decreased recombination in specific maps, we identified several enriched gene ontologies in these regions and speculate about their local adaptive relevance. These data are contributing to an increasing comparative effort to gain an understanding of the intra-specific variability of recombination rates and their evolutionary role in honey bees and other social insects.

Comparative transcriptomics of social insect queen pheromones

Queen pheromones are chemical signals that mediate reproductive division of labor in eusocial animals. Remarkably, queen pheromones are composed of identical or chemically similar compounds in some ants, wasps and bees, even though these taxa diverged >150MYA and evolved queens and workers independently. Here, we measure the transcriptomic consequences of experimental exposure to queen pheromones in workers from two ant and two bee species (genera: Lasius, Apis, Bombus), and test whether they are similar across species. Queen pheromone exposure affected transcription and splicing at many loci. Many genes responded consistently in multiple species, and the set of pheromone-sensitive genes was enriched for functions relating to lipid biosynthesis and transport, olfaction, production of cuticle, oogenesis, and histone (de)acetylation. Pheromone-sensitive genes tend to be evolutionarily ancient, positively selected, peripheral in the gene coexpression network, hypomethylated, and caste-specific in their expression. Our results reveal how queen pheromones achieve their effects, and suggest that ants and bees use similar genetic modules to achieve reproductive division of labor.

Proteomics Improves the New Understanding of Honeybee Biology

The honeybee is one of the most valuable insect pollinators, playing a key role in pollinating wild vegetation and agricultural crops, with significant contribution to the world's food production. Although honeybees have long been studied as model for social evolution, honeybee biology at the molecular level remained poorly understood until the year 2006. With the availability of the honeybee genome sequence and technological advancements in protein separation, mass spectrometry, and bioinformatics, aspects of honeybee biology such as developmental biology, physiology, behavior, neurobiology, and immunology have been explored to new depths at molecular and biochemical levels. This Review comprehensively summarizes the recent progress in honeybee biology using proteomics to study developmental physiology, task transition, and physiological changes in some of the organs, tissues, and cells based on achievements from the authors' laboratory in this field. The research advances of honeybee proteomics provide new insights for understanding of honeybee biology and future research directions.

Shotgun proteomics deciphered age/division of labor-related functional specification of three honeybee (Apis mellifera L.) exocrine glands

The honeybee (Apis mellifera L.) uses various chemical signals produced by the worker exocrine glands to maintain the functioning of its colony. The roles of worker postcerebral glands (PcGs), thoracic glands (TGs), and mandibular glands (MGs) and the functional changes they undergo according to the division of labor from nursing to foraging are not as well studied. To comprehensively characterize the molecular roles of these glands in workers and their changes according to the division of labor of workers, we analyzed the proteomes of PcGs, TGs, and MGs from nurse bees and foragers using shotgun proteomics technology. We identified approximately 2000 proteins from each of the nurse bee or forager glands and highlighted the features of these glands at the molecular level by semiquantitative enrichment analyses of frequently detected, gland-selective, and labor-selective proteins. First, we found the high potential to produce lipids in PcGs and MGs, suggesting their relation to pheromone production. Second, we also found the proton pumps abundant in TGs and propose some transporters possibly related to the saliva production. Finally, our data unveiled candidate enzymes involved in labor-dependent acid production in MGs.

Comparative transcriptome analysis on the synthesis pathway of honey bee (Apis mellifera) mandibular gland secretions

Secretions from mandibular glands (MGs) have important caste-specific functions that are associated with the social evolution of honey bees. To gain insights into the molecular architecture underlying these caste differences, we compared the gene expression patterns of MGs from queens, queenright workers (WQRs) and queenless workers (WQLs) using high-throughput RNA-sequencing technology. In total, we identified 46 candidate genes associated with caste-specific biosynthesis of fatty acid pheromones in the MG, including members of cytochrome P450 (CYP450) family and genes involved in fatty acid β-oxidation and ω-oxidation. For further identification of the CYP450s genes involved in the biosynthesis of MG secretions, we analyzed by means of qPCR, the expression levels of six of the CYP450 genes most abundantly expressed in the transcriptome analysis across different castes, ages, tasks and tissues. Our analysis revealed that CYP6AS8 and CYP6AS11, the most abundantly expressed CYP450 genes in worker and queen MGs, respectively, are selectively expressed in the MGs of workers and queens compared to other tissues. These results suggest that these genes might be responsible for the critical bifurcated hydroxylation process in the biosynthesis pathway. Our study contributes to the description of the molecular basis for the biosynthesis of fatty acid-derived pheromones in the MGs.

Contribution of lipids in honeybee (Apis mellifera) royal jelly to health

Honeybee (Apis mellifera) royal jelly (RJ) has a long history in human medicine because of its health-protecting properties. To develop a fundamental and comprehensive understanding of lipids in RJ, this article reviews the available literature on lipid compounds identified from RJ extracts and in vitro pharmacological effects of 10-hydroxy-2-decenoic acid in RJ and other closely related compounds, some of which are also identified as lipid compounds in RJ. Overall, the lipids in RJ are composed of mostly (aliphatic) fatty acids, almost all of which are present as free fatty acids and scarcely any as esters. Most fatty acids in RJ are medium-chain fatty acids, whether hydroxylated in terminal and/or internal positions, terminated with mono- or dicarboxylic acid groups, and saturated or monounsaturated at the 2-position. Besides these fatty acids, lipids in RJ contain sterols in minor amounts. Lipids in RJ are useful as preventive and supportive medicines with functionalities that include potential inhibitors of cancer growth, immune system modulators, alternative therapies for menopause, skin-aging protectors, neurogenesis inducers, and more. Taken together, the evidence suggests that health-protecting properties of RJ can be, in part, ascribed to actions of lipids in RJ.

Transcriptional changes associated with lack of lipid synthesis in parasitoids

Phenotypic regression of morphological, behavioral, or physiological traits can evolve when reduced trait expression has neutral or beneficial effects on overall performance. Studies on the evolution of phenotypic degradation in animals have concentrated mostly on the evaluation of resulting phenotypes, whereas much less research has been dedicated to uncovering the molecular mechanisms that underlie phenotypic regression. The majority of parasitoids (i.e., insects that develop on or inside other arthropods), do not accumulate lipid reserves during their free-living adult life-stage and represent an excellent system to study phenotypic regression in animals. Here, we study transcriptional patterns associated with lack of lipogenesis in the parasitic wasp Nasonia vitripennis. We first confirmed that N. vitripennis does not synthesize lipids by showing a reduction in lipid reserves despite ingestion of dietary sugar, and a lack of incorporation of isotopic labels into lipid reserves when fed deuterated sugar solution. Second, we investigated transcriptional responses of 28 genes involved in lipid and sugar metabolism in short- and long-term sugar-fed females relative to starved females of N. vitripennis. Sugar feeding did not induce transcription of fatty acid synthase (fas) or other key genes involved in the lipid biosynthesis pathway. Furthermore, several genes involved in carbohydrate metabolism had a lower transcription in fed than in starved females. Our results reveal that N. vitripennis gene transcription in response to dietary sugar deviates markedly from patterns typically observed in other organisms. This study is the first to identify differential gene transcription associated with lack of lipogenesis in parasitoids and provides new insights into the molecular mechanism that underlies phenotypic regression of this trait.

Social molecular pathways and the evolution of bee societies

Bees provide an excellent model with which to study the neuronal and molecular modifications associated with the evolution of sociality because relatively closely related species differ profoundly in social behaviour, from solitary to highly social. The recent development of powerful genomic tools and resources has set the stage for studying the social behaviour of bees in molecular terms. We review 'ground plan' and 'genetic toolkit' models which hypothesize that discrete pathways or sets of genes that regulate fundamental behavioural and physiological processes in solitary species have been co-opted to regulate complex social behaviours in social species. We further develop these models and propose that these conserved pathways and genes may be incorporated into 'social pathways', which consist of relatively independent modules involved in social signal detection, integration and processing within the nervous and endocrine systems, and subsequent behavioural outputs. Modifications within modules or in their connections result in the evolution of novel behavioural patterns. We describe how the evolution of pheromonal regulation of social pathways may lead to the expression of behaviour under new social contexts, and review plasticity in circadian rhythms as an example for a social pathway with a modular structure.

Proteomic analysis of honeybee (Apis mellifera L.) pupae head development

The honeybee pupae development influences its future adult condition as well as honey and royal jelly productions. However, the molecular mechanism that regulates honeybee pupae head metamorphosis is still poorly understood. To further our understand of the associated molecular mechanism, we investigated the protein change of the honeybee pupae head at 5 time-points using 2-D electrophoresis, mass spectrometry, bioinformatics, quantitative real-time polymerase chain reaction and Western blot analysis. Accordingly, 58 protein spots altered their expression across the 5 time points (13–20 days), of which 36 proteins involved in the head organogenesis were upregulated during early stages (13–17 days). However, 22 proteins involved in regulating the pupae head neuron and gland development were upregulated at later developmental stages (19–20 days). Also, the functional enrichment analysis further suggests that proteins related to carbohydrate metabolism and energy production, development, cytoskeleton and protein folding were highly involved in the generation of organs and development of honeybee pupal head. Furthermore, the constructed protein interaction network predicted 33 proteins acting as key nodes of honeybee pupae head growth of which 9 and 4 proteins were validated at gene and protein levels, respectively. In this study, we uncovered potential protein species involved in the formation of honeybee pupae head development along with their specific temporal requirements. This first proteomic result allows deeper understanding of the proteome profile changes during honeybee pupae head development and provides important potential candidate proteins for future reverse genetic research on honeybee pupae head development to improve the performance of related organs.

Differential protein expression in honeybee (Apis mellifera L.) larvae: underlying caste differentiation

Honeybee (Apis mellifera) exhibits divisions in both morphology and reproduction. The queen is larger in size and fully developed sexually, while the worker bees are smaller in size and nearly infertile. To better understand the specific time and underlying molecular mechanisms of caste differentiation, the proteomic profiles of larvae intended to grow into queen and worker castes were compared at 72 and 120 hours using two dimensional electrophoresis (2-DE), network, enrichment and quantitative PCR analysis. There were significant differences in protein expression between the two larvae castes at 72 and 120 hours, suggesting the queen and the worker larvae have already decided their fate before 72 hours. Specifically, at 72 hours, queen intended larvae over-expressed transketolase, aldehyde reductase, and enolase proteins which are involved in carbohydrate metabolism and energy production, imaginal disc growth factor 4 which is a developmental related protein, long-chain-fatty-acid CoA ligase and proteasome subunit alpha type 5 which metabolize fatty and amino acids, while worker intended larvae over-expressed ATP synthase beta subunit, aldehyde dehydrogenase, thioredoxin peroxidase 1 and peroxiredoxin 2540, lethal (2) 37 and 14-3-3 protein epsilon, fatty acid binding protein, and translational controlled tumor protein. This differential protein expression between the two caste intended larvae was more pronounced at 120 hours, with particular significant differences in proteins associated with carbohydrate metabolism and energy production. Functional enrichment analysis suggests that carbohydrate metabolism and energy production and anti-oxidation proteins play major roles in the formation of caste divergence. The constructed network and validated gene expression identified target proteins for further functional study. This new finding is in contrast to the existing notion that 72 hour old larvae has bipotential and can develop into either queen or worker based on epigenetics and can help us to gain new insight into the time of departure as well as caste trajectory influencing elements at the molecular level.