Proteome Comparisons between Hemolymph of Two Honeybee Strains (Apis mellifera ligustica) Reveal Divergent Molecular Basis in Driving Hemolymph Function and High Royal Jelly Secretion

Sublethal acetamiprid exposure induces immunity, suppresses pathways linked to juvenile hormone synthesis in queens and affects cycle-related signaling in emerging bees

Acetamiprid is the only neonicotinoid registered in the European Union because the risks of neonicotinoids to honey bees and other pollinators are strictly regulated. Herein, we orally exposed honey bee colonies to sublethal concentrations of acetamiprid (20 μg/L) under isolated conditions. After one month of continuous exposure, the emerging bees and queens were collected and analyzed via high-throughput label-free quantitative proteomics using a data-independent acquisition strategy. Six and 34 significantly differentially expressed proteins (DEPs) were identified in the emerging bees and queens, respectively. Mrjp3 was the only DEP found in both sample types/castes, and its opposite regulation illustrated a differential response. The DEPs in the emerging bees (H/ACA RNP, Rap1GAP, Mrjp3, and JHE) suggested that sublethal exposure to acetamiprid affected cell cycle-related signaling, which may affect the life history of workers in the colony. The DEPs with increased levels in queens, such as Mrjps 1-4 and 6-7, hymenoptaecin, and apidaecin 22, indicated an activated immune response. Additionally, the level of farnesyl pyrophosphate synthase (FPPS), which is essential for the mevalonate pathway and juvenile hormone biosynthesis, was significantly decreased in queens. The impaired utilization of juvenile hormone in queens supported the identification of additional DEPs. Furthermore, the proteome changes suggested the existence of increased neonicotinoid detoxification by UDP-glucuronosyltransferase and increased amino acid metabolism. The results suggest that the continuous exposure of bee colonies to acetamiprid at low doses (nanograms per gram in feed) may pose a threat to the colonies. The different exposure routes and durations for the emerging bees and queens in our experiment must be considered, i.e., the emerging bees were exposed as larvae via feeding royal jelly and beebread provided by workers (nurse bees), whereas the queens were fed royal jelly throughout the experiment. The biological consequences of the proteomic changes resulting from sublethal/chronic exposure require future determination.

Profiling the neuroproteomics of honeybee brain: A clue for understanding the role of neuropeptides in the modulation of aggressivity

The aggressivity is modulated in honeybee brain through a series of actions in cascade mode, with the participation of the neuropeptides AmAST A (59-76) and AmTRP (254-262). The aggressivity of honeybees was stimulated by injecting both neuropeptides in the hemocoel of the worker honeybees, which were submitted to behavioral assays of aggression. The brain of stinger individuals were removed by dissection and submitted to proteomic analysis; shotgun proteomic approach of honeybee brain revealed that both neuropeptides activate a series of biochemical processes responsible by production of energy, neuronal plasticity and cell protection. In addition to this, AmTRP (254-262) elicited the expression of proteins related to the processing of the potential of action and lipid metabolism; meanwhile AmAST A (59-76) elicited the metabolism of steroids and Juvenile hormone-related metabolism, amongst others. Apparently, the most complex biochemical process seems to be the regulation of ATP production, which occurs at two levels: i) by a subgroup of proteins common to the three experimental groups, which are over-/under-regulated through glycolysis, pyruvate pathway, Krebbs cycle and oxidative phosphorylation; ii) by a subgroup of proteins unique to the each experimental group, which seems to be regulated through Protein-Protein Interactions, where the protein network regulated by AmTRP (254-262) seems to be more complex than the other two experimental groups. SIGNIFICANCE: Recently we reported the effect of the neuropeptides AmAST A (59-76) and AmTRP (254-262) in the modulation of the aggressive behavior of the worker honeybees. Up to now it is known that the simple presence of the allatostatin and tachykinin-related-peptide in bee brain, is enough for inducing the aggressive behavior. However, nothing was known about how these neuropeptides perform their action, inducing the aggressive behavior. The results of the present study elucidated some of the metabolic pathways that were activated or inhibited to support the complex defensive behavior, which includes the aggressivity. These results certainly will impact the behavioral research of honeybees, since we are paving the way for understanding the molecular base of regulation, of individual /nest defense of honeybees.

SDS-PAGE-Based Quantitative Assay of Hemolymph Proteins in Honeybees: Progress and Prospects for Field Application

In human and veterinary medicine, serum proteins are considered to be useful biomarkers for assessing the health and nutritional status of the organism. Honeybee hemolymph has a unique proteome that could represent a source of valuable biomarkers. Therefore, the aims of this study were to separate and identify the most abundant proteins in the hemolymph of worker honeybees to suggest a panel of these proteins that could represent useful biomarkers for assessing the nutritional and health status of the colonies and, finally, to analyze them in different periods of the year. Four apiaries were selected in the province of Bologna, and the bees were analyzed in April, May, July, and November. Thirty specimens from three hives of each apiary were sampled and their hemolymph was collected. The most represented bands obtained after 1D sodium-dodecyl-sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) were cut from the gel, and the proteins were identified using an LC-ESI-Q-MS/MS System. A total of twelve proteins were unmistakably identified; the two most abundant proteins were apolipophorin and vitellogenin, which are known biomarkers of bee trophic and health status. The two other proteins identified were transferrin and hexamerin 70a, the first being involved in iron homeostasis and the second being a storage protein. Most of these proteins showed an increase from April to November, mirroring the physiological changes of honeybees during the productive season. The current study suggests a panel of biomarkers from honeybee hemolymph worth testing under different physiological and pathological field conditions.

Transcriptomic Responses Underlying the High Virulence of Black Queen Cell Virus and Sacbrood Virus following a Change in Their Mode of Transmission in Honey Bees (Apis mellifera)

BACKGROUND

Over the last two decades, honey bees (Apis mellifera) have suffered high rates of colony losses that have been attributed to a variety of factors, chief among which are viral pathogens, such as deformed wing virus (DWV), whose virulence has increased because of vector-based transmission by the invasive, ectoparasitic varroa mite (Varroa destructor). A shift in the experimental mode of transmission of the black queen cell virus (BQCV) and sacbrood virus (SBV) from fecal/food-oral (direct horizontal) to vector-mediated (indirect horizontal) transmission also results in high virulence and viral titers in pupal and adult honey bees. Agricultural pesticides represent another factor that acts independently or in interaction with pathogens, and they are also thought to cause colony loss. Understanding the molecular mechanisms underlying the higher virulence following a vector-based mode of transmission provides deeper insight into honey bee colony losses, as does determining whether or not host-pathogen interactions are modulated by exposure to pesticides.

METHODS

Through an experimental design with controlled laboratory, we investigated the effects of the modes of transmission of BQCV and SBV (feeding vs. vector-mediated via injection) alone or in combination with chronic exposure to sublethal and field-realistic concentrations of flupyradifurone (FPF), a novel agricultural insecticide, on honey bee survival and transcription responses by using high-throughput RNA sequencing (RNA-seq) analysis.

RESULTS

Co-exposure to viruses via feeding (VF) or injection (VI) and FPF insecticide had no statistically significant interactive effect on their survival compared to, respectively, VF or VI treatments alone. Transcriptomic analysis revealed a distinct difference in the gene expression profiles of bees inoculated with viruses via injection (VI) and exposed to FPF insecticide (VI+FPF). The number of differentially expressed genes (DEGs) at log2 (fold-change) > 2.0 in VI bees (136 genes) or/and VI+FPF insecticide (282 genes) was very high compared to that of VF bees (8 genes) or the VF+FPF insecticide treatment (15 genes). Of these DEGs, the expression in VI and VI+FPF bees of some immune-related genes, such as those for antimicrobial peptides, Ago2, and Dicer, was induced. In short, several genes encoding odorant binding proteins, chemosensory proteins, odor receptors, honey bee venom peptides, and vitellogenin were downregulated in VI and VI+FPF bees.

CONCLUSIONS

Given the importance of these suppressed genes in honey bees' innate immunity, eicosanoid biosynthesis, and olfactory associative function, their inhibition because of the change in the mode of infection with BQCV and SBV to vector-mediated transmission (injection into haemocoel) could explain the high virulence observed in these viruses when they were experimentally injected into hosts. These changes may help explain why other viruses, such as DWV, represent such a threat to colony survival when transmitted by varroa mites.

Tropilaelaps mercedesae Infestation Is Correlated with Injury Numbers on the Brood and the Population Size of Honey Bee Apis mellifera

Tropilaelaps mercedesae, one of the most devastating parasitic mites of honey bee Apis mellifera hosts, is a major threat to honey products by causing severe damage to honey bee colonies. Here, we recorded injury numbers caused by T. mercedesae to different body parts of the larval, pupal, and crippled adult stages of honey bee A. mellifera. We evaluated the relationship between infestation rate and injury numbers per bee for both larvae and pupae. We also noted the total bee numbers per beehive and examined the relationship between the infestation rate and population size. T. mercedesae infested all developmental stages of honey bees, with the highest injury numbers in the abdomens of bee pupae and the antennas of crippled adult bees. Although larvae received more injury numbers than pupae, both infestation rate and injury numbers decreased as the larval stage progressed to the pupal stage. The infestation rate increased as the population size per beehive decreased. This study provided new perspectives to the understanding of changes in the effects of T. mercedesae infestations on different developmental stages of honey bees. It also showed useful baseline information for screening honey bee stock that might have high defensive behaviors against mite infestation.

Unraveling the Bombus terrestris Hemolymph, an Indicator of the Immune Response to Microbial Infections, through Complementary Mass Spectrometry Approaches

Pollinators, including Bombus terrestris, are crucial for maintaining biodiversity in ecosystems and for agriculture. Deciphering their immune response under stress conditions is a key issue for protecting these populations. To assess this metric, we analyzed the B. terrestris hemolymph as an indicator of their immune status. Hemolymph analysis was carried out using mass spectrometry, MALDI molecular mass fingerprinting was used for its effectiveness in assessing the immune status, and high-resolution mass spectrometry was used to measure the impact of experimental bacterial infections on the "hemoproteome". By infecting with three different types of bacteria, we observed that B. terrestris reacts in a specific way to bacterial attacks. Indeed, bacteria impact survival and stimulate an immune response in infected individuals, visible through changes in the molecular composition of their hemolymph. The characterization and label-free quantification of proteins involved in specific signaling pathways in bumble bees by bottom-up proteomics revealed differences in protein expression between the non-experimentally infected and the infected bees. Our results highlight the alteration of pathways involved in immune and defense reactions, stress, and energetic metabolism. Lastly, we developed molecular signatures reflecting the health status of B. terrestris to pave the way for diagnosis/prognosis tools in response to environmental stress.

Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole-genome sequencing data

Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole-genome queen genotypes, from a pool-seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well-defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F ROH) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F ROH (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F ROH (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F ROH of 2.33%. Within the breed-specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh-BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.

Effect of queen cell numbers on royal jelly production and quality

Royal jelly (RJ) is a popular functional food with a wealth of health-promoting effects. Over 90% of the global RJ is produced in China mainly by a high RJ-producing honeybee (RJB) strain that can accept and feed a great number of queen larvae for RJ production. To elucidate RJ changes due to queen cell numbers (QCNs), we compared the yield, larval acceptance rate, metabolic and proteomic profiles, and antioxidant activities of RJ from 1 to 5 strips of queen cells (64 per strip) in RJB colonies. As QCNs increased, the larval acceptance rate was not found to vary (p = 0.269) whereas the RJ weight per cell began to significantly decline in the 5-strip colonies (p < 0.05). Increased QCNs had a profound impact on RJ metabolic profiles and mainly reduced fatty acid levels. Remarkably, the 10-hydroxy-2-decenoic acid (10-HDA) content, a most important indicator of RJ quality, declined gradually from 2.01% in the 1-strip colonies to 1.52% in the 5-strip colonies (p < 0.001). RJ proteomic profiles were minimally altered and antioxidant activities were not significantly changed by QCNs. Collectively, the metabolomics and proteomics data and the antioxidant activity test represent a global evaluation of the quality of RJ produced with different QCNs. Our findings gain new insights into higher-quality RJ production using the high-yielding RJBs.

Changes in chemical composition and antioxidant activity of royal jelly produced at different floral periods during migratory beekeeping

Over 90% of global royal jelly (RJ), a functional food with various health benefits, is produced in China mainly by migratory beekeeping of a high RJ-producing honeybee (RJB) strain. To explore quality changes of RJ produced by migratory RJBs at different floral periods, we performed metabolomics and proteomics analysis and assessed RJ antioxidant activity. Overall, the RJ metabolic and proteomic profiles were observed to vary with floral periods. Minor sugars (raffinose, erlose, and sucrose) and major RJ protein 5 (MRJP5) were identified among the discriminating components mainly contributing to the altered profiles. Water, crude protein, and the trans-10-hydroxy-2-decenoic acid (10-HDA) content fulfill the requirements of the International Organization for Standardization regardless of floral periods. Notably, the 10-HDA content increased 11.05%-19.65% during tea blooming. Moreover, changes in antioxidants resulted in significant difference in RJ antioxidant activity. The integrated omics data provide a detailed view of chemical composition for RJ quality evaluation.

Genetic Drift and Purifying Selection Shaped Mitochondrial Genome Variation in the High Royal Jelly-Producing Honeybee Strain (Apis mellifera ligustica)

Mitochondrial genomes (mitogenomes) are involved in cellular energy metabolism and have been shown to undergo adaptive evolution in organisms with increased energy-consuming activities. The genetically selected high royal jelly-producing bees (RJBs, Apis mellifera ligustica) in China can produce 10 times more royal jelly, a highly nutritional and functional food, relative to unselected Italian bees (ITBs). To test for potential adaptive evolution of RJB mitochondrial genes, we sequenced mitogenomes from 100 RJBs and 30 ITBs. Haplotype network and phylogenetic analysis indicate that RJBs and ITBs are not reciprocally monophyletic but mainly divided into the RJB- and ITB-dominant sublineages. The RJB-dominant sublineage proportion is 6-fold higher in RJBs (84/100) than in ITBs (4/30), which is mainly attributable to genetic drift rather than positive selection. The RJB-dominant sublineage exhibits a low genetic diversity due to purifying selection. Moreover, mitogenome abundance is not significantly different between RJBs and ITBs, thereby rejecting the association between mitogenome copy number and royal jelly-producing performance. Our findings demonstrate low genetic diversity levels of RJB mitogenomes and reveal genetic drift and purifying selection as potential forces driving RJB mitogenome evolution.

High royal jelly production does not impact the gut microbiome of honey bees

Background

Honey bees are not only essential for pollination services, but are also economically important as a source of hive products (e.g., honey, royal jelly, pollen, wax, and propolis) that are used as foods, cosmetics, and alternative medicines. Royal jelly is a popular honey bee product with multiple potential medicinal properties. To boost royal jelly production, a long-term genetic selection program of Italian honey bees (ITBs) in China has been performed, resulting in honey bee stocks (here referred to as RJBs) that produce an order of magnitude more royal jelly than ITBs. Although multiple studies have investigated the molecular basis of increased royal jelly yields, one factor that has not been considered is the role of honey bee-associated gut microbes.

Results

Based on the behavioral, morphological, physiological, and neurological differences between RJBs and ITBs, we predicted that the gut microbiome composition of RJBs bees would differ from ITBs. To test this hypothesis, we investigated the bacterial composition of RJB and ITB workers from an urban location and RJBs from a rural location in China. Based on 16S rRNA gene profiling, we did not find any evidence that RJBs possess a unique bacterial gut community when compared to ITBs. However, we observed differences between honey bees from the urban versus rural sites.

Conclusions

Our results suggest that the environmental factors rather than stock differences are more important in shaping the bacterial composition in honey bee guts. Further studies are needed to investigate if the observed differences in relative abundance of taxa between the urban and rural bees correspond to distinct functional capabilities that impact honey bee health. Because the lifestyle, diet, and other environmental variables are different in rural and urban areas, controlled studies are needed to determine which of these factors are responsible for the observed differences in gut bacterial composition between urban and rural honeybees.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00124-1.

Metabolic profiling unravels the effects of enhanced output and harvesting time on royal jelly quality

Royal jelly (RJ), a natural secretion of nurse bees, is a functional food with health-promoting properties. To investigate RJ quality changes with enhanced output and harvesting time, we compared metabolic profiles of RJ harvested at 24 h, 48 h, and 72 h after larval grafting into queen cells from Italian bees (ITBs) and high RJ-producing bees (RJBs) selected from ITBs. A total of 77 high-abundance compounds were identified with reverse-phase liquid chromatography- and hydrophilic interaction liquid chromatography-high-resolution mass spectrometry. Metabolite abundance, water content, and antibacterial activity were found largely unchanged between both bee strains, indicating that a similar RJ quality was achieved. Notably, the observed 10-hydroxy-2-decenoic acid content in RJBs (1.98%) and ITBs (2.41%) meet the standard for premium RJ in China. RJ harvested at 48 h and 72 h had similar metabolic profiles, indicating little influence of the two harvesting times on RJ quality. The harvesting at 24 h could lead to a higher RJ quality with higher abundance levels of many bioactive compounds and lower water content. Our data based mainly on RJ metabolic composition constitute evidence for mass production of high-quality RJ.

Comparative analysis of the gut microbiota of Apis cerana in Yunnan using high-throughput sequencing

Gut microbes play an important role in host disease and health. The Asian honey bee Apis cerana is an important pollinator of agricultural crops in China. However, there are still few studies on the structure and composition of the microbiota in the intestine of A. cerana, especially A. cerana in Yunnan. To understand the species and composition of the microbiota in the intestine of A. cerana in Yunnan, we used high-throughput sequencing technology to carry out 16S rRNA sequencing on 50 samples from Kunming, Xishuangbanna and Mengzi. The results show that both from the phylum level and the genus level, the structure and abundance of the microbiota in the gut of A. cerana from the three regions tended to be the same. At the phylum level, the abundance of Proteobacteria, Firmicutes, Bacteroidetes, Actinobacteria, Acidobacteria and other species was high in A. cerana from different areas. At the genus level, the abundance of Lactobacillus, Gilliamella, Snodgrassella, Apibacter, Candidatus Schmidhempelia and other species was high in A. cerana from different areas. Due to its unique geographical environment and climatic conditions, at the genus level, the diversity of bacterial communities in Xishuangbanna was significantly lower than that in the other two regions, which was about 100 genera less. In conclusion, our results reveal the composition and structure of the intestinal microbiota of bees in Yunnan and deepen our understanding of the intestinal microbiota of bees.

The Emerging Proteomic Research Facilitates in-Depth Understanding of the Biology of Honeybees

Advances in instrumentation and computational analysis in proteomics have opened new doors for honeybee biological research at the molecular and biochemical levels. Proteomics has greatly expanded the understanding of honeybee biology since its introduction in 2005, through which key signaling pathways and proteins that drive honeybee development and behavioral physiology have been identified. This is critical for downstream mechanistic investigation by knocking a gene down/out or overexpressing it and being able to attribute a specific phenotype/biochemical change to that gene. Here, we review how emerging proteome research has contributed to the new understanding of honeybee biology. A systematic and comprehensive analysis of global scientific progress in honeybee proteome research is essential for a better understanding of research topics and trends, and is potentially useful for future research directions.

Varroa destructor parasitism has a greater effect on proteome changes than the deformed wing virus and activates TGF-β signaling pathways

Honeybee workers undergo metamorphosis in capped cells for approximately 13 days before adult emergence. During the same period, Varroa mites prick the defenseless host many times. We sought to identify proteome differences between emerging Varroa-parasitized and parasite-free honeybees showing the presence or absence of clinical signs of deformed wing virus (DWV) in the capped cells. A label-free proteomic analysis utilizing nanoLC coupled with an Orbitrap Fusion Tribrid mass spectrometer provided a quantitative comparison of 2316 protein hits. Redundancy analysis (RDA) showed that the combination of Varroa parasitism and DWV clinical signs caused proteome changes that occurred in the same direction as those of Varroa alone and were approximately two-fold higher. Furthermore, proteome changes associated with DWV signs alone were positioned above Varroa in the RDA. Multiple markers indicate that Varroa activates TGF-β-induced pathways to suppress wound healing and the immune response and that the collective action of stressors intensifies these effects. Furthermore, we indicate JAK/STAT hyperactivation, p53-BCL-6 feedback loop disruption, Wnt pathway activation, Wnt/Hippo crosstalk disruption, and NF-κB and JAK/STAT signaling conflict in the Varroa–honeybee–DWV interaction. These results illustrate the higher effect of Varroa than of DWV at the time of emergence. Markers for future research are provided.

pH-dependent stability of honey bee (Apis mellifera) major royal jelly proteins

Honey bee larval food jelly is a secretion of the hypopharyngeal and mandibular glands of young worker bees that take care of the growing brood in the hive. Food jelly is fed to all larvae (workers, drones and queens) and as royal jelly to the queen bee for her entire life. Up to 18% of the food jelly account for proteins the majority of which belongs to the major royal jelly protein (MRJP) family. These proteins are produced in the hypopharyngeal glands at a pH value of 7.0. Before being fed to the larvae, they are mixed with the fatty acids secreted by the mandibular glands of the worker bees resulting at a pH of 4.0 in the food jelly. Thus, MRJPs are exposed to a broad pH range from their site of synthesis to the actual secreted larval food. We therefore determined the pH-dependent stability of MRJP1, MRJP2 and MRJP3 purified from royal jelly using differential scanning fluorimetry. All MRJPs were much more stable at acidic pH values compared to neutral ones with all proteins showing highest stability at pH 4.0 or 4.5, the native pH of royal jelly.

In-depth Proteome of the Hypopharyngeal Glands of Honeybee Workers Reveals Highly Activated Protein and Energy Metabolism in Priming the Secretion of Royal Jelly

Royal jelly (RJ) is a secretion of the hypopharyngeal glands (HGs) of honeybee workers. High royal jelly producing bees (RJBs), a stock of honeybees selected from Italian bees (ITBs), have developed a stronger ability to produce RJ than ITBs. However, the mechanism underpinning the high RJ-producing performance in RJBs is still poorly understood. We have comprehensively characterized and compared the proteome across the life span of worker bees between the ITBs and RJBs. Our data uncover distinct molecular landscapes that regulate the gland ontogeny and activity corresponding with age-specific tasks. Nurse bees (NBs) have a well-developed acini morphology and cytoskeleton of secretory cells in HGs to prime the gland activities of RJ secretion. In RJB NBs, pathways involved in protein synthesis and energy metabolism are functionally induced to cement the enhanced RJ secretion compared with ITBs. In behavior-manipulated RJB NBs, the strongly expressed proteins implicated in protein synthesis and energy metabolism further demonstrate their critical roles in the regulation of RJ secretion. Our findings provide a novel understanding of the mechanism consolidating the high RJ-output in RJBs.

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.