Identification of the proteome composition occurring during the course of embryonic development of bees (Apis mellifera)

Egg-size plasticity in Apis mellifera: Honey bee queens alter egg size in response to both genetic and environmental factors

Social evolution has led to distinct life-history patterns in social insects, but many colony-level and individual traits, such as egg size, are not sufficiently understood. Thus, a series of experiments was performed to study the effects of genotypes, colony size and colony nutrition on variation in egg size produced by honey bee (Apis mellifera) queens. Queens from different genetic stocks produced significantly different egg sizes under similar environmental conditions, indicating standing genetic variation for egg size that allows for adaptive evolutionary change. Further investigations revealed that eggs produced by queens in large colonies were consistently smaller than eggs produced in small colonies, and queens dynamically adjusted egg size in relation to colony size. Similarly, queens increased egg size in response to food deprivation. These results could not be solely explained by different numbers of eggs produced in the different circumstances but instead seem to reflect an active adjustment of resource allocation by the queen in response to colony conditions. As a result, larger eggs experienced higher subsequent survival than smaller eggs, suggesting that honey bee queens might increase egg size under unfavourable conditions to enhance brood survival and to minimize costly brood care of eggs that fail to successfully develop, and thus conserve energy at the colony level. The extensive plasticity and genetic variation of egg size in honey bees has important implications for understanding life-history evolution in a social context and implies this neglected life-history stage in honey bees may have trans-generational effects.

Mass Spectrometry-Based Analysis of Time-Resolved Proteome Quantification

The aspect of time is essential in biological processes and thus it is important to be able to monitor signaling molecules through time. Proteins are key players in cellular signaling and they respond to many stimuli and change their expression in many time-dependent processes. Mass spectrometry (MS) is an important tool for studying proteins, including their posttranslational modifications and their interaction partners-both in qualitative and quantitative ways. In order to distinguish the different trends over time, proteins, modification sites, and interacting proteins must be compared between different time points, and therefore relative quantification is preferred. In this review, the progress and challenges for MS-based analysis of time-resolved proteome dynamics are discussed. Further, aspects on model systems, technologies, sampling frequencies, and presentation of the dynamic data are discussed.

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.

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.

Proteomic analysis of trochophore and veliger larvae development in the small abalone Haliotis diversicolor

Background

Haliotis diversicolor is commercially important species. The trochophore and veliger are distinct larval stages in gastropod development. Their development involves complex morphological and physiological changes. We studied protein changes during the embryonic development of H. diversicolor using two dimensional electrophoresis (2-DE) and label-free methods, tandem mass spectrometry (MS/ MS), and Mascot for protein identification.

Results

A total of 150 2-DE gel spots were identified. Protein spots showed upregulation of 15 proteins and downregulation of 28 proteins as H. diversicolor developed from trochophore to veliger larvae. Trochophore and veliger larvae were compared using a label-free quantitative proteomic approach. A total of 526 proteins were identified from both samples, and 104 proteins were differentially expressed (> 1.5 fold). Compared with trochophore larvae, veliger larvae had 55 proteins upregulated and 49 proteins downregulated. These differentially expressed proteins were involved in shell formation, energy metabolism, cellular and stress response processes, protein synthesis and folding, cell cycle, and cell fate determination. Compared with the 5 protein (fructose-bisphosphate aldolase, 14–3-3ε, profilin, actin-depolymerizing factor (ADF)/cofilin) and calreticulin) expression patterns, the mRNA expression exhibited similar patterns except gene of fructose-bisphosphate aldolase.

Conclusion

Our results provide insight into novel aspects of protein function in shell formation, torsion, and nervous system development, and muscle system differentiation in H. diversicolor larvae. “Quality control” proteins were identified to be involved in abalone larval development.

Electronic supplementary material

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

Proteomic and metabolomic analysis reveals rapid and extensive nicotine detoxification ability in honey bee larvae

Despite potential links between pesticides and bee declines, toxicology information on honey bee larvae (Apis mellifera) is scarce and detoxification mechanisms in this development stage are virtually unknown. Larvae are exposed to natural and synthetic toxins present in pollen and nectar through consumption of brood food. Due to the characteristic intensive brood care displayed by honey bees, which includes progressive feeding throughout larval development, it is generally assumed that larvae rely on adults to detoxify for them and exhibit a diminished detoxification ability. We found the opposite. We examined the proteomic and metabolomic responses of in vitro reared larvae fed nicotine (an alkaloid found in nectar and pollen) to understand how larvae cope on a metabolic level with dietary toxins. Larvae were able to effectively detoxify nicotine through an inducible detoxification mechanism. A coordinated stress response complemented the detoxification processes, and we detected significant enrichment of proteins functioning in energy and carbohydrate metabolism, as well as in development pathways, suggesting that nicotine may promote larval growth. Further exploration of the metabolic fate of nicotine using targeted mass spectrometry analysis demonstrated that, as in adult bees, formation of 4-hydroxy-4-(3-pyridyl) butanoic acid, the result of 2'C-oxidation of nicotine, is quantitatively the most significant pathway of nicotine metabolism. We provide conclusive evidence that larvae are capable of effectively catabolising a dietary toxin, suggesting that increased larval sensitivity to specific toxins is not due to diminished detoxification abilities. These findings broaden the current understanding of detoxification biochemistry at different organizational levels in the colony, bringing us closer to understanding the capacity of the colony as a superorganism to tolerate and resist toxic compounds, including pesticides, in the environment.

Proteome Analysis Unravels Mechanism Underling the Embryogenesis of the Honeybee Drone and Its Divergence with the Worker (Apis mellifera lingustica)

The worker and drone bees each contain a separate diploid and haploid genetic makeup, respectively. Mechanisms regulating the embryogenesis of the drone and its mechanistic difference with the worker are still poorly understood. The proteomes of the two embryos at three time-points throughout development were analyzed by applying mass spectrometry-based proteomics. We identified 2788 and 2840 proteins in the worker and drone embryos, respectively. The age-dependent proteome driving the drone embryogenesis generally follows the worker's. The two embryos however evolve a distinct proteome setting to prime their respective embryogenesis. The strongly expressed proteins and pathways related to transcriptional-translational machinery and morphogenesis at 24 h drone embryo relative to the worker, illustrating the earlier occurrence of morphogenesis in the drone than worker. These morphogenesis differences remain through to the middle-late stage in the two embryos. The two embryos employ distinct antioxidant mechanisms coinciding with the temporal-difference organogenesis. The drone embryo's strongly expressed cytoskeletal proteins signify key roles to match its large body size. The RNAi induced knockdown of the ribosomal protein offers evidence for the functional investigation of gene regulating of honeybee embryogenesis. The data significantly expand novel regulatory mechanisms governing the embryogenesis, which is potentially important for honeybee and other insects.

Proteomic analysis of Apis cerana and Apis mellifera larvae fed with heterospecific royal jelly and by CSBV challenge

Chinese honeybee Apis cerana (Ac) is one of the major Asian honeybee species for local apiculture. However, Ac is frequently damaged by Chinese sacbrood virus (CSBV), whereas Apis mellifera (Am) is usually resistant to it. Heterospecific royal jelly (RJ) breeding in two honeybee species may result in morphological and genetic modification. Nevertheless, knowledge on the resistant mechanism of Am to this deadly disease is still unknown. In the present study, heterospecific RJ breeding was conducted to determine the effects of food change on the larval mortality after CSBV infection at early larval stage. 2-DE and MALDI-TOF/TOF MS proteomic technology was employed to unravel the molecular event of the bees under heterospecific RJ breeding and CSBV challenge. The change of Ac larval food from RJC to RJM could enhance the bee resistance to CSBV. The mortality rate of Ac larvae after CSBV infection was much higher when the larvae were fed with RJC compared with the larvae fed with RJM. There were 101 proteins with altered expressions after heterospecific RJ breeding and viral infection. In Ac larvae, 6 differential expression proteins were identified from heterospecific RJ breeding only, 21 differential expression proteins from CSBV challenge only and 7 differential expression proteins from heterospecific RJ breeding plus CSBV challenge. In Am larvae, 17 differential expression proteins were identified from heterospecific RJ breeding only, 26 differential expression proteins from CSBV challenge only and 24 differential expression proteins from heterospecific RJ breeding plus CSBV challenge. The RJM may protect Ac larvae from CSBV infection, probably by activating the genes in energy metabolism pathways, antioxidation and ubiquitin-proteasome system. The present results, for the first time, comprehensively descript the molecular events of the viral infection of Ac and Am after heterospecific RJ breeding and are potentially useful for establishing CSBV resistant populations of Ac for apiculture.

In-depth proteomics characterization of embryogenesis of the honey bee worker (Apis mellifera ligustica)

Identifying proteome changes of honey bee embryogenesis is of prime importance for unraveling the molecular mechanisms that they underlie. However, many proteomic changes during the embryonic period are not well characterized. We analyzed the proteomic alterations over the complete time course of honey bee worker embryogenesis at 24, 48, and 72 h of age, using mass spectrometry-based proteomics, label-free quantitation, and bioinformatics. Of the 1460 proteins identified the embryo of all three ages, the core proteome (proteins shared by the embryos of all three ages, accounting for 40%) was mainly involved in protein synthesis, metabolic energy, development, and molecular transporter, which indicates their centrality in driving embryogenesis. However, embryos at different developmental stages have their own specific proteome and pathway signatures to coordinate and modulate developmental events. The young embryos (<24 h) stronger expression of proteins related to nutrition storage and nucleic acid metabolism may correlate with the cell proliferation occurring at this stage. The middle aged embryos (24-48 h) enhanced expression of proteins associated with cell cycle control, transporters, antioxidant activity, and the cytoskeleton suggest their roles to support rudimentary organogenesis. Among these proteins, the biological pathways of aminoacyl-tRNA biosynthesis, β-alanine metabolism, and protein export are intensively activated in the embryos of middle age. The old embryos (48-72 h) elevated expression of proteins implicated in fatty acid metabolism and morphogenesis indicate their functionality for the formation and development of organs and dorsal closure, in which the biological pathways of fatty acid metabolism and RNA transport are highly activated. These findings add novel understanding to the molecular details of honey bee embryogenesis, in which the programmed activation of the proteome matches with the physiological transition observed during embryogenesis. The identified biological pathways and key node proteins allow for further functional analysis and genetic manipulation for both the honey bee embryos and other eusocial insects.

The identification and oxidative stress response of a zeta class glutathione S-transferase (GSTZ1) gene from Apis cerana cerana

Glutathione-S-transferases (GSTs) play an important role in protecting organisms against the toxicity of reactive oxygen species (ROS). However, no information is available for GSTs in the Chinese honey bee (Apis cerana cerana). In this study, we isolated and characterized a zeta class GST gene (AccGSTZ1) from the Chinese honey bee. This gene is present in a single copy and harbors five exons. The deduced amino acid sequence of AccGSTZ1 shared high sequence identity with homologous proteins and contained the highly conserved features of this gene family. The temporal and spatial expression profiles of AccGSTZ1 showed that AccGSTZ1 was highly expressed in fourth instar larvae during development, and the mRNA level of AccGSTZ1 was higher in the epidermis than that in other tissues. The expression pattern under oxidative stress revealed that AccGSTZ1 transcription was significantly upregulated by external factors, such as temperature challenges and H(2)O(2) treatment. The characterization of the purified protein revealed that AccGSTZ1 had low glutathione-conjugating activity, but the recombinant AccGSTZ1 protein displayed high antioxidant activity under oxidative stress. These data suggest that AccGSTZ1 is an oxidative stress-inducible antioxidant enzyme that plays an important role in the protection against oxidative stress and may be of critical importance for the survival of the honey bees.

AccERK2, a map kinase gene from Apis cerana cerana, plays roles in stress responses, developmental processes, and the nervous system

Extracellular signal-regulated kinase (ERK), a mitogen-activated protein kinase (MAPK), plays roles in a variety of cellular responses. However, limited information is available on the relationship between ERKs and environmental stresses. In this report, an ERK gene, AccERK2, was cloned and characterized from Apis cerana cerana. Polypeptide sequence alignment revealed that the single-copied AccERK2 shares high identity with other known ERKs and contains the typical conserved Thr-Glu-Tyr (TEY) motif in its activation loop. Genomic sequence analysis revealed that the seven exons of AccERK2 are interrupted by six introns, and the seventh intron is located in the 3' untranslated region. Semi-quantitative reverse transcription (RT-PCR) indicated that AccERK2 was expressed at higher levels in the larval and pupal stages than in the adult stage. AccERK2 was also most highly expressed in the brain. The expression of AccERK2 was induced by abiotic stresses, including heat, ion irradiation, oxidative stress, and heavy metal ions. Based on these results, it appears that AccERK2 in A. cerana cerana participates in developmental processes, the nervous system, and responses to environmental stressors.

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.

Honeybee (Apis mellifera ligustica) drone embryo proteomes

Little attention has been paid to the drone honeybee (Apis mellifera ligustica) which is a haploid individual carrying only the set of alleles that it inherits from its mother. Molecular mechanisms underlying drone embryogenesis are poorly understood. This study evaluated protein expression profiles of drone embryogenesis at embryonic ages of 24, 48 and 72h. More than 100 reproducible proteins were analyzed by mass spectrometry on 2D electrophoresis gels. Sixty-two proteins were significantly changed at the selected three experimental age points. Expression of the metabolic energy requirement-related protein peaked at the embryonic age of 48h, whereas development and metabolizing amino acid-related proteins expressed optimally at 72h. Cytoskeleton, protein folding and antioxidant-related proteins were highly expressed at 48 and 72h. Protein networks of the identified proteins were constructed and protein expressions were validated at the transcription level. This first proteomic study of drone embryogenesis in the honeybee may provide geneticists an exact timetable and candidate protein outline for further manipulations of drone stem cells.

Proteome comparison of hypopharyngeal gland development between Italian and royal jelly producing worker honeybees (Apis mellifera L.)

The hypopharyngeal gland (HG) of the honeybee (Apis mellifera L.) produces royal jelly (RJ) that is essential to feed and raise broods and queens. A strain of bees (high royal jelly producing bee, RJb) has been selected for its high RJ production, but the mechanisms of its higher yield are not understood. In this study, we compared HG acini size, RJ production, and protein differential expressions between the RJb and nonselected honeybee (Italian bee, ITb) using proteomics in combination with an electron microscopy, Western blot, and quantitative real-time PCR (qRT-PCR). Generally, the HG of both bees showed age-dependent changes in acini sizes and protein expression as worker behaviors changed from brood nursing to nectar ripening, foraging, and storage activities. The electron microscopic analysis revealed that the HG acini diameter of the RJb strain was large and produced 5 times more RJ than the ITb, demonstrating a positive correlation between the yield and HG acini size. In addition, the proteomic analysis showed that RJb significantly upregulated a large group of proteins involved in carbohydrate metabolism and energy production, those involved in protein biosynthesis, development, amino acid metabolism, nucleotide and fatty acid, transporter, protein folding, cytoskeleton, and antioxidation, which coincides with the fact that the HGs of the RJb strain produce more RJ than the ITb strain that is owing to selection pressure. We also observed age-dependent major royal jelly proteins (MRJPs) changing both in form and expressional intensity concurrent with task-switching. In addition to MRJPs, the RJb overexpressed proteins such as enolase and transitional endoplasmic reticulum ATPase, protein biosynthesis, and development proteins compared to the ITb strain to support its large HG growth and RJ secretion. Because of selection pressure, RJb pursued a different strategy of increased RJ production by involving additional proteins compared to its original counterpart ITb. To our knowledge, this morphological and proteomic comparison study on the HG of the two strains of worker honeybees associated with their age-dependent division of labor is the first of its kind. The study provided not only the quantity and quality differences in the HG from the RJb and the ITb, but also addressed the cellular and behavioral biology development question of how the RJb strain can produce RJ more efficiently than its wild type strain (ITb).

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.

Proteomic analysis of honeybee worker (Apis mellifera) hypopharyngeal gland development

Background

Hypopharyngeal glands (HG) of honeybee workers play an important role in honeybee nutrition and caste differentiation. Previous research mainly focused on age-dependent morphological, physiological, biochemical and genomic characters of the HG. Here proteomics and biochemical network analysis were used to follow protein changes during the HG development.

Results

A total of 87, 76, 85, 74, 71, and 55 proteins were unambiguously identified on day 1, 3, 6, 12, 15 and 20, respectively. These proteins were major royal jelly proteins (MRJPs), metabolism of carbohydrates, lipids and proteins, cytoskeleton, development regulation, antioxidant, molecule transporter, regulation of transcription/translation, proteins with folding functions. The most interesting is that MRJP's that have been detected in the HG of the newly emerged worker bees. The MRJP's expression is at peak level from 6-12 days, was validated by western blot analysis of MRJP1, 2 and 3. Moreover, 35 key node proteins were found in the biochemical networks of the HG.

Conclusions

HG secretes RJ at peak level within 6-12 days, but the worker bee can secrete royal jelly (RJ) since birth, which is a new finding. Several key node proteins play an important role in the biochemical networks of the developing HG. This provides us some target proteins when genetically manipulating honeybees.