Next Article in Journal
Nematodes in the Pine Forests of Northern and Central Greece
Previous Article in Journal
Organophosphate Insecticides Resistance in Field Populations of House Flies, Musca domestica L.: Levels of Resistance and Acetylcholinesterase Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

1H NMR Profiling of Honey Bee Bodies Revealed Metabolic Differences between Summer and Winter Bees

1
Department of Food Science, Faculty of Agrobiology, Food and Natural Resources, Czech University of Life Sciences Prague, Kamycka 129, 165 00 Prague, Czech Republic
2
Campus de Santiago, University of Aveiro (CICECO), 3810-19 Aveiro, Portugal
3
Department of Zoology and Fisheries, Faculty of Agrobiology, Food and Natural Resources, Czech University of Life Sciences Prague, Kamycka 129, 165 00 Prague, Czech Republic
4
Honeybee Research Institute, Dol 94, 252 66 Maslovice, Czech Republic
5
Department of Experimental Biology, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
6
Department of Biochemistry, Faculty of Science, Palacky University Olomouc, Slechtitelu 27, 783 71 Olomouc, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 15 December 2021 / Revised: 27 January 2022 / Accepted: 8 February 2022 / Published: 12 February 2022
(This article belongs to the Section Insect Societies and Sociality)

Abstract

:

Simple Summary

The European honey bee, Apis mellifera, is well-known to have two distinct populations in temperate climate zone: short-living summer bees and long-living winter bees. Several biological factors related to the different lifespans of the two populations have been studied. However, the link between the metabolic changes and basic physiological features in the bodies of summer bees and winter bees is limited. This study aimed to identify the metabolic fingerprints that characterize summer and winter bees using proton nuclear magnetic resonance (1H NMR) spectroscopy. In total, we found 28 significantly changed metabolites between the two populations. The results suggest that the metabolites detected in honey bee bodies can distinguish the summer and winter bees. Changes in carbohydrates, amino acids, choline-containing compounds, and an unknown compound were noticeable during the transition from summer bees to winter bees. The results from this study give us a broad perspective on honey bee metabolism that will support future research related to honey bee lifespan and overwintering management.

Abstract

In temperate climates, honey bee workers of the species Apis mellifera have different lifespans depending on the seasonal phenotype: summer bees (short lifespan) and winter bees (long lifespan). Many studies have revealed the biochemical parameters involved in the lifespan differentiation of summer and winter bees. However, comprehensive information regarding the metabolic changes occurring in their bodies between the two is limited. This study used proton nuclear magnetic resonance (1H NMR) spectroscopy to analyze the metabolic differences between summer and winter bees of the same age. The multivariate analysis showed that summer and winter bees could be distinguished based on their metabolic profiles. Among the 36 metabolites found, 28 metabolites have displayed significant changes from summer to winter bees. Compared to summer bees, trehalose in winter bees showed 1.9 times higher concentration, and all amino acids except for proline and alanine showed decreased patterns. We have also detected an unknown compound, with a CH3 singlet at 2.83 ppm, which is a potential biomarker that is about 13 times higher in summer bees. Our results show that the metabolites in summer and winter bees have distinctive characteristics; this information could provide new insights and support further studies on honey bee longevity and overwintering.

1. Introduction

The European honey bee (Apis mellifera) is one of the most important pollinators for agriculture and helps support ecological diversity [1]. However, the steadily reported decline in honey bee populations’ ability to overwinter in various regions, including North America and Europe, is a worrisome phenomenon [2,3,4,5]. The elevated overwintering mortality as observed in the last few decades [6] is probably caused by several biotic and abiotic stressors, such as Varroa destructor [7], Nosema sp. [8], viruses [9], pesticides [10], diet and nutrition [11], poor winter conditions [12], and climate change [13]. As winter is a critical period for honey bee colonies, understanding the positive and negative environmental factors that influence colony mortality will be a crucial task to saving honey bees.
In temperate climates, European honey bee workers emerge into roles as either short-living summer bee or long-living winter bee populations, depending on the emergence period and labor division within the colony. Summer bees mainly take care of their brood, collect nectar, pollen, and water, and have a short life span which lasts an average of 15 to 38 days [14]. When brood rearing begins to cease in the autumn, the winter bees begin to emerge [15,16]. Their main tasks are to keep the colony warm throughout the winter and rear the first summer generation around late winter or in early spring [17], and live approximately 140 and up to 320 days [14].
The differences in longevity are linked to several physiological and biochemical characteristics. The protein vitellogenin, juvenile hormone, and their interconnections are one of the most well-known mechanisms that influence the aging process in honey bees. Winter bees have low levels of juvenile hormone along with high levels of vitellogenin [18]; this is similar to queen bees, who have the longest lifespan among the colony members [19]. However, summer bees show the opposite pattern, high levels of juvenile hormone and low levels of vitellogenin [18]. Vitellogenin appears to play a role that affects the aging and lifespan of honey bees in beneficial ways such as its antioxidant function, its role in immunity, and controlling inflammation [20].
Apart from the differences in vitellogenin and juvenile hormone, winter bees have a greater dry weight and a higher amount of total proteins, glucose, glycogen, and lipids compared to summer bees [21,22]. Additionally, they have higher antibacterial activity [22] and different gut microbiota profiles, which appear to be affected by dietary differences in comparison to summer bees [23]. Moreover, winter bees showed a negligible level of carbonylation damage in the brain [24] and negligible cognitive senescence [25] when compared to summer foragers.
Although many studies have explained the characteristics of winter bees, their biochemical pathways and associated mechanisms have not been fully elucidated. Hence, screening the biomarkers that link the physiological factors that underlie aging and lifespan in honey bees is important considering that the emergence of the long-living population in late summer is crucial for successful overwintering. That could infer the diagnostic biomarkers that are potentially influenced by both positive and negative environmental factors.
Metabolomics entails the comprehensive profiling of metabolites in biofluids, cells, and tissues, which potentially enables biomarker discovery and can provide new insights into the mechanisms that determine various physiological conditions [26]. Nuclear magnetic resonance (NMR) offers an exquisitely informative way to obtain a holistic quantitative picture of small molecules [27]. This approach has been employed in the field of entomology, analyzing fruit flies [28], locusts [29], parasitoid wasps [30], aphids [31], and recently in honey bees [32].
The aim of this study is to apply NMR spectroscopy-based metabolomics to identify the metabolic fingerprints that characterize summer and winter bees of the same age in the temperate zone of the Northern Hemisphere. The period of sample collecting, late June and late August for summer and winter bees, respectively, ensures a consistency in the age of bees. This allows for an unbiased comparison of the transition of the metabolites from short-living populations to long-living populations. The results are intended to provide additional information that links previous research regarding the lifespan of honey bees and to propose potential indicator metabolites necessary for successful overwintering in honey bee colonies.

2. Materials and Methods

2.1. Sample Collection

Honey bees (Apis mellifera carnica) were collected in 2017 over two days, 27th June and 29th August, from three hives at the Honeybee Research Institute, Kyvalka, Czech Republic (49°11′28.7″ N 16°26′58.4″ E). The honey bees were kept under typical bee management conditions and fed sucrose solution (3:2 w/v sucrose:water) for overwintering on August 8th, 19th, and 24th; no other supplements were added to the hives. The colonies did not show any clinical symptoms of diseases or Varroa infestation at the time of sampling. Samples were collected as a subset of another study published earlier [33]. The first sample collection in June represented normal brood rearing activities, whereas the second sample collection in August showed reduced brood activities, which is a sign that the colony is entering the winter bee stage.
To prepare the same age of honey bees, the method from [34] was applied. Three isolated frames (39 × 24 cm) containing capped brood that will emerge in one to three days were placed in frame cages and kept in three hives. Newly emerged bees were collected from each hive, marked with a color on their thoraxes, then returned to their respective hives. After 10 days of living in the hives, the marked bees were re-collected. The re-collected honey bees were transferred to a laboratory and kept in cages at 30 °C for 24 h on 60% sucrose solution to remove the bias of short-term dietary influences. Thus, the total age of the bees was 11 days. Upon collection, whole honey bees were stored at −80 °C until the time of the analysis. Previously, we have demonstrated on the samples from the same sampling that honey bees collected in August showed significantly higher concentrations of vitellogenin and vitellogenin gene expression (p < 0.005) compared with the honey bees collected in June [33]. The color-marked honey bees from August continued to be found in the hive until February. Therefore, this study considers these samples as winter bee population.

2.2. Sample Preparation for NMR Analysis

Only the body (head, thorax with legs and wings; and abdomen) was used for NMR analysis. Venom sacks and gastrointestinal tracts of the bees were removed with tweezers and excluded. Honey bee bodies were individually homogenized in a 2.0 mL Eppendorf tube by freezing with liquid nitrogen and then ground to a fine powder by a Dremel MultiTool with an attached burr (Dremel, Mount Prospect, WI, USA). The total number of samples included in the analysis was 120. June samples were 75 (25 per hive) and August samples were 45 (15 per hive). Homogenized bodies were extracted with 1.2 mL of methanol in an ultrasonic bath for five minutes, centrifuged at 14,000× g for five minutes, and evaporated in a centrifugal vacuum concentrator for two hours (MiVac Duo, Genevac, Ipswitch, UK). After evaporation, the sample was re-suspended in 540 µL of D2O. Subsequently, 60 µL of NMR solution (1.5 M K2HPO4/NaH2PO4 pH 7.4 in D2O, 5 mM Trimethylsilylpropanoic acid, and 0.2% NaN3) was added and vortexed, and the sample was transferred into NMR tubes (5 mm, 7”, High-Throughput, SP Wilmad-Labglass, Wineland, NJ, USA).

2.3. 1H NMR Analysis

Standard 1D 1H NMR spectra (‘noesypr1d’ Bruker pulse program) of honey bee methanol extracts reconstituted in D2O were acquired on a Bruker Avance III spectrometer (Bruker, Billerica, MA, USA) operating at 500.23 MHz for 1H observation, at 298 K, and using a 5 mm Broadband obseve probe. The spectra were acquired with 4.9 s acquisition time, 1 s relaxation delay (D1), 0.1 s mixing time, 6510.42 Hz spectral width, 64 K data points, and 64 scans.

2.4. Data Processing and Statistical Analysis

All spectra were zero-filled to 128 K, exponentially multiplied (line broadening 0.3 Hz), manually phased, and baseline-corrected using Topspin 3.5 (Bruker, Billerica, MA, US). For metabolite identification and quantification, we used Chenomx NMR Suite 8.6 (Chenomx, Edmonton, AB, Canada) comprising the Reference Library 10 and our in-house library of approx. 80 compounds. The additional process for identifying metabolites was confirmed through the standard compound analysis (Sigma-Aldrich, St. Louis, MO, USA) and VWR International, Stribrna Skalice, Czech Republic). A set of 36 metabolites was fitted to the spectra and data were manually curated. Compounds below the limit of detection were replaced with one-third of the minimum concentration detected in the dataset. Metabolite concentrations were expressed in µg per bee. The average weight of a bee was 157 mg, of which 48 mg was a digestive tract that was not used for analysis.
This study used R ver.4.1.2 (RStudio, Boston, MA, USA) for the multivariate analysis. Data were processed using probabilistic quotient normalization, log transformation, and Pareto scaling by the package ‘MetaboAnalystR’ [35]. Orthogonal partial least squares-discriminant analysis (OPLS-DA) with permutation test validation was applied in MetaboAnalyst 4.0 [36]. The packages ‘FactoMineR’ [37] and ‘ggplot2′ [38] were used for the principal component analysis and its visualization.
An independent t-test was used to compare the concentration of metabolites in summer bees and winter bees using SPSS statistics ver. 26 (IBM, Armonk, NY, USA), and those metabolites with p < 0.05 after Benjamini–Hochberg false discovery rate correction were considered statistically significant. Comparative visualizations for summer and winter bees were performed using the package ‘ggplot2′ [38] in R.

3. Results

3.1. NMR Spectrum of Honey Bee’s Body

A total of 36 metabolites were identified in the NMR spectra of the honey bee methanolic extracts: 13 amino acids, 4 carboxylic acids, 5 amines, 3 choline-containing compounds, 3 nucleosides, 2 nucleotides, 2 monosaccharides, 2 disaccharides, 1 nucleobase, and 1 unknown compound.
Figure 1 shows a typical 1H NMR spectrum of honey bees’ body with the expansion of the signal between the 6.0 and 9.0 ppm region. The region from 0.5 to 3.5 ppm contains peaks mainly from amino acids, organic acids, and lipids. The region between 3.5 and 5.5 ppm contains resonances of mainly carbohydrates, and the region beyond 6.0 ppm contain mainly aromatic compounds.

3.2. Multivariate Analysis Explaining the Discrimination between Summer (June) and Winter Bees (Late August)

This study posed a major question: Can summer and winter bees be discriminated by metabolites? If so, what are the key metabolites that could separate these two groups? To answer the question, we used sample set of worker bees comprised of 75 summer bees (June) and 45 winter bees (late August) both of equal age of 11 days.
In this study, the principal component analysis (PCA) showed the separation between summer and winter bees, and the correlated variables that contribute to each of the principal components (Figure 2). The first two principal components explain 54.7% of the total variation in the data (PC1 42.6%, PC2 12.1%) and have an eigenvalue > 1 in PC1. Summer and winter bees showed a notable separation where summer bees formed a cluster to about half of the PC1 whereas, winter bees more tightly clustered over the PC2 than summer bees (Figure 2a).
The significant variables which separate the summer and winter bees are presented in the variable correlation plot (Figure 2b). The arrow represents the amount of variance positively or negatively correlated on the principal component. The color and length refers the contribution of variables on the principal component in percentage [37]. Among the top five distinguished variables on the PC1, an unknown compound had the highest contribution at 11.6%, followed by sucrose, phenylalanine, trimethylamine N-oxide, and leucine (10.7%, 7.2%, 6.9%, and 6.3%, respectively). Fructose, sucrose, and glucose were the majority of variables on the PC2 which contributed 23.2%, 20.7%, and 17.9%, respectively (Figure 2b, Table S1). Taken together, the unknown compound, sucrose, fructose, glucose, phenylalanine, trimethylamine N-oxide, and leucine were the meaningful variables which contributed to the separation between summer and winter bees. In addition, the permutation test from the orthogonal partial least squares-discriminant analysis (OPLS-DA) confirmed the high validity at Q2 = 0.908, R2Y = 0.946 when calculated from the 1000 different model permutations.

3.3. Univariate Analysis Confirmed Statistically Significant Metabolites

After evaluating the link between the honey bees in summer and winter with certain metabolites by PCA (Figure 2), an independent t-test with Benjamini–Hochberg correction was applied to demonstrate the changes in concentrations of individual metabolites between summer and winter bees. From the 36 metabolites, 28 were significantly different between these two groups (Figure 3, Table 1).
The concentration of all the carbohydrates (fructose p = 0.0001, glucose p > 0.05, sucrose p < 0.0001, trehalose p < 0.0001, Figure 3a) increased in winter bees compared to summer bees. For example, trehalose in winter bees (801.8 ± 212.0 μg/bee, Table 1) showed 1.9 times higher concentration than in summer bees (416.0 ± 181.1 μg/bee, Table 1).
In the amino acids group, all amino acids (Asparagine p < 0.0001, glutamine p = 0.0001, isoleucine p < 0.0001, leucine p < 0.0001, lysine p < 0.0001, phenylalanine p < 0.0001, threonine p < 0.0001, and valine p < 0.0001, Figure 3b) except for alanine and proline showed a noticeable decrease in winter bees. Phenylalanine showed the largest difference, with levels 3.4 times higher in summer bees (4.39 ± 2.53 μg/bee, Table 1) than in winter bees (1.31 ± 0.55 μg/bee, Table 1). Alanine (p = 0.0001, Figure 3b) and proline (p < 0.0001, Figure 3b) were significantly increased in winter bees compared to summer bees.
For the choline-containing compounds, O-phosphocholine and choline presented opposite trends (Figure 3c). The level of O-phosphocholine showed a statistically significant increase in winter bees (p < 0.0001, Figure 3c), while the levels of choline showed a statistically significant decrease in winter bees (p < 0.0001, Figure 3c).
The group of nucleobases, nucleotides, and nucleosides are presented in Figure 3d. AMP, adenosine, and NAD+ showed a statistically significant increase in winter bees compared to summer bees (all p-value < 0.0001, Figure 3d) while adenine, inosine, and uridine decreased (Adenine: p < 0.0001, inosine: p < 0.0001, uridine: p = 0.0499, Figure 3d). Adenine showed the largest differences between the groups at 2.8 times higher in summer bees (1.29 ± 0.82 μg/bee, Table 1) than in winter bees (0.46 ± 0.25 μg/bee, Table 1).
The rest of the metabolites, such as carboxylic acids and amines, are presented as a group in Figure 3e. Trimethylamine N-oxide showed notable differences between the two groups, where the concentration in winter bees (2.97 ± 1.54 μg/bee, Table 1) was 3.5 times higher than in summer bees (0.86 ± 0.63 μg/bee, Table 1). Most noteworthy is that the unknown compound, a methyl singlet at 2.83 ppm, showed the greatest change among all metabolites (p < 0.0001, Figure 3e), which was about 13 times more concentrated in summer bees (19.61 ± 7.62 μg/bee, Table 1) compared to winter bees (1.53 ± 0.55 μg/bee, Table 1).

4. Discussion

This study utilized NMR metabolomics to investigate whether summer bees and winter bees could be distinguished by specific metabolites and identified that 28 metabolites significantly changed during the transition from summer to winter bees. Among the 28 metabolites distinguishing summer and winter bees, the unknown compound is by far the most noticeable. As confirmed in the PCA (Figure 2, Table S1), the unknown compound, showing a singlet peak in the spectrum at 2.83 ppm is an important indicator of the distinction between the two groups. This metabolite has a particularly high concentration in the body of summer bees, while the presence is negligible in winter bees (Figure 3e, Table 1). A preliminary elucidation suggests that the signal belongs to N-methyl group of a non-proteinogenic amino acid and appears to be a powerful biomarker of seasonality or longevity. Structural elucidation is currently ongoing.
Concerning carbohydrates, winter bees had a significantly higher concentration of fructose, sucrose, and trehalose than in summer bees (Figure 3a, Table 1). Carbohydrates are a critical source of energy for bees, and are the exclusive fuel used for flight [39,40]. They can be converted into glycogen and stored in the fat body [41]. Higher carbohydrate intake was translated into higher glycogen levels and better diapause survival in reared bumble bees [42]. Trehalose is the primary blood sugar in most insects and is involved in multiple physiological roles [43]. One of which is a role as a stress protector. Studies of several insects have shown that stress conditions, such as drought, cold stress, and diapause, lead to trehalose accumulation in the body [44,45,46]. We observed increased levels of trehalose in the winter bees’ bodies (Figure 3a), this could be attributed to osmoregulation or enhancing cold stress tolerance in preparation for winter. Supporting this idea, an overwintering study of the mountain pine beetle (Dendroctonus ponderosae) using 1H NMR showed that the level of trehalose begins to increase from September until December as temperatures declined [47], a similar pattern has been shown in European spruce bark beetles (Ips typographus) [44], as well as in honey bees [48].
The standard deviation of sucrose, glucose, and fructose was higher than the mean value (Table 1). This is possibly originating from the crop leak during the sample disection or individual differences in sucrose feeding before the bees were harvested for this experiement.
We observed that most free amino acids except for alanine and proline were significantly lower in winter compared to summer bees (Figure 3b). Although winter bees are characterized by generally higher protein content, vitellogenin concentration, and more antibacterial peptides than summer bees [16,21,22,33], their protein turnover is lower [49]. In contrast, the June samples were nurse bees, producing royal jelly proteins during the active brood rearing period [50,51], requiring a high protein synthesis and demand for free amino acids.
Proline is the predominant amino acid in the hemolymph of honey bees [52,53] and is known as a cryoprotectant of insects, such as trehalose mentioned above [54,55]. In freeze tolerance studies in fruit flies (Drosophila melanogaster), proline and trehalose were the most abundant metabolites that accumulated during the cold assimilation [55,56,57]. We speculate that increased levels of trehalose and proline in the winter bees act as cryoprotectants that prevent membrane structure damages occuring during the cold winter, maintain water levels, or serve as an energy pool.
Environmental stressors such as pathogen exposure, starvation, and temperature fluctuation affect changes in insect metabolites. Winter bees emerge in autumn when brood rearing slows in the hive due to a drop in available pollen resources [15,58]. Such environmental changes could affect the metabolites in honey bees; this may be a part of their longevity mechanism. Low protein intake is known to be linked with lifespan extension [59]. A study found that honey bees kept on a diet low in essential amino acids showed longer lifespans through the Sir2 expression than those consuming a high essential amino acid diet [60]. Consistent with this result, Gomez-Moracho et al. [61] reported that insects fed a low-protein, high-carbohydrate diet showed the highest parasite prevalence but had the highest longevity. The decrease in amino acids in winter bees can be a contributing factor to the longevity mechanism in the long-living winter bees. Thus, future research focused on the metabolites correlation with the pathways that are known to regulate lifespan would expand the knowledge of honey bees’ longevity.
Vitellogenin is one of the most robust biochemical indicators in honey bees regarding their aging and longevity [18,20,62]. It is a phospholipoglycoprotein, a major precursor of egg-yolk protein synthesized by the fat body in female insects [20]. Vitellogenin contains phospholipids, mainly phosphatidylcholine and lysophosphatidylcholine [63,64]. In this study, we found significant changes in detected levels of choline and O-Phosphocholine in winter bees compared with summer bees (Figure 3c). Considering these two metabolites are the intermediates of phosphatidylcholine synthesis [65], we hypothesize that the changes in the choline-containing compounds from summer bees to winter bees could be related with the phosphatidylcholine synthesis as a vitellogenin component.
The levels of adenosine monophosphate (AMP) and oxidized nicotinamide adenine dinucleotide (NAD+) have significantly increased in winter bees compared to summer bees (Figure 3d). Changes in AMP and NAD+ concentrations indicate cellular energy status and are meaningful for the study of longevity due to their correlation with the activation of AMP-activated protein kinase (AMPK) and sirtuins [66,67]. Sirtuins are NAD+-dependent protein deacylases. Increased NAD+ levels can enhance activated sirtuins that improve metabolic function and longevity [67]. AMP is one of the AMPK activators, an important energy sensor that is activated by increasing cellular AMP coupled with falling ATP [66]. Overexpression of AMPK has been shown to extend lifespan in Caenorhabditis elegans [68] and Drosophila melanogaster [69]. Previous honey bee studies have shown a lower concentration of NAD+ and decreased activity of AMPK and sirtuin in the abdomen of older bees compared to that of younger bees [70,71]. Therefore, we hypothesize that the significant increase in AMP and NAD+ in winter bees is one of the signals of switching to the long-living generation. Thus, further research is needed on this subject to more accruately link AMP and NAD+ with the regulation of longevity pathways.
For successful overwintering practices, we consider it is important to examine the metabolic changes from short-living summer bees to long-living winter bees. The molecular mechanisms and the detailed metabolic pathways regarding the changes in metabolites was not discerned during this experiment and require further research. The highlight of this study is that summer and winter bees have distinctive metabolic characteristics, and some metabolites could be useful biomarkers to distinguish these two groups.

5. Conclusions

This study observed and analyzed the metabolic differences between summer and winter bees of the same age. We identified and quantified a total of 36 metabolites, 28 metabolites which showed significant differences between the two groups.
Our results have the potential to provide a new analytical approach to confirm the presence of a long-living population in honey bee colonies. Furthermore, knowing robust biomarkers between summer and winter bees will guide future research on honey bee longevity and overwintering, and help to develop a practical overwintering strategy for the apiculture industry.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/insects13020193/s1, Table S1: Contributions of the variables to the principal components in percentage (%).

Author Contributions

Conceptualization, J.H. and P.H.; methodology, J.H., P.H., P.D., M.K. (Martin Kunc), and S.L.; software, S.L. and I.F.D.; validation, J.H., P.H., P.D. and M.K. (Martin Kunc); formal analysis, S.L, and P.M.; investigation, S.L, P.M., A.P. and F.K.; resources, D.T. and M.K. (Martin Kamler); data curation, S.L. and J.H.; writing—original draft preparation, S.L.; writing—review and editing, S.L., F.K, P.H., P.D., M.K. (Martin Kunc), I.F.D., J.D. and J.H.; visualization, S.L.; supervision, J.H.; project administration, J.H. and P.H.; funding acquisition, P.H., I.F.D. and J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the project of the Ministry of Agriculture of the Czech Republic (QJ1610248 and QK1910286), by the European Regional Development Fund-Project No. CZ.02.1.01/0.0/0.0/16_019/0000845 and by a METROFOOD-CZ research infrastructure project (MEYS Grant No: LM2018100) including access to its facilities. IFD acknowledges the project CICECO-Aveiro Institute of Materials (UIDB/50011/2020 and UIDP/50011/2020) financed by national funds through the Foundation for Science and Technology/MCTES, Portugal.

Institutional Review Board Statement

Not applicable.

Acknowledgments

The authors thank Zuzana Hroncova for technical assistance. We also thank the reviewers for the valuable critical assessments and helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Potts, S.G.; Biesmeijer, J.C.; Kremen, C.; Neumann, P.; Schweiger, O.; Kunin, W.E. Global pollinator declines: Trends, impacts and drivers. Trends Ecol. Evol. 2010, 25, 345–353. [Google Scholar] [CrossRef]
  2. Brodschneider, R.; Gray, A.; Adjlane, N.; Ballis, A.; Brusbardis, V.; Charrière, J.D.; Chlebo, R.; Coffey, M.F.; Dahle, B.; de Graaf, D.C.; et al. Multi-country loss rates of honey bee colonies during winter 2016/2017 from the COLOSS survey. J. Apic. Res. 2018, 57, 452–457. [Google Scholar] [CrossRef] [Green Version]
  3. Vanengelsdorp, D.; Caron, D.; Hayes, J.; Underwood, R.; Henson, M.; Rennich, K.; Spleen, A.; Andree, M.; Snyder, R.; Lee, K.; et al. A national survey of managed honey bee 2010–11 winter colony losses in the USA: Results from the Bee Informed Partnership. J. Apic. Res. 2012, 51, 115–124. [Google Scholar] [CrossRef] [Green Version]
  4. Gray, A.; Brodschneider, R.; Adjlane, N.; Ballis, A.; Brusbardis, V.; Charrière, J.-D.; Chlebo, R.; Coffey, M.F.; Cornelissen, B.; da Costa, C.A.; et al. Loss rates of honey bee colonies during winter 2017/18 in 36 countries participating in the COLOSS survey including effects of forage sources. J. Apic. Res. 2019, 58, 479–485. [Google Scholar] [CrossRef] [Green Version]
  5. Steinhauer, N.A.; Rennich, K.; Wilson, M.E.; Caron, D.M.; Lengerich, E.J.; Pettis, J.S.; Rose, R.; Skinner, J.A.; Tarpy, D.R.; Wilkes, J.T.; et al. A national survey of managed honey bee 2012–2013 annual colony losses in the USA: Results from the Bee Informed Partnership. J. Apic. Res. 2014, 53, 1–18. [Google Scholar] [CrossRef]
  6. Van Engelsdorp, D.; Evans, J.D.; Saegerman, C.; Mullin, C.; Haubruge, E.; Nguyen, B.K.; Frazier, M.; Frazier, J.; Cox-Foster, D.; Chen, Y.; et al. Colony collapse disorder: A descriptive study. PLoS ONE 2009, 4, e6481. [Google Scholar] [CrossRef]
  7. Van Dooremalen, C.; Gerritsen, L.; Cornelissen, B.; van der Steen, J.J.M.; van Langevelde, F.; Blacquière, T. Winter survival of individual honey bees and honey bee colonies depends on level of varroa destructor infestation. PLoS ONE 2012, 7, e36285. [Google Scholar] [CrossRef] [Green Version]
  8. Higes, M.; Martín-Hernández, R.; Botías, C.; Bailón, E.G.; González-Porto, A.V.; Barrios, L.; Del Nozal, M.J.; Bernal, J.L.; Jiménez, J.J.; Palencia, P.G.; et al. How natural infection by Nosema ceranae causes honeybee colony collapse. Environ. Microbiol. 2008, 10, 2659–2669. [Google Scholar] [CrossRef] [PubMed]
  9. Highfield, A.C.; El Nagar, A.; Mackinder, L.C.M.; Noël, L.M.L.J.; Hall, M.J.; Martin, S.J.; Schroeder, D.C. Deformed wing virus implicated in overwintering honeybee colony losses. Appl. Environ. Microbiol. 2009, 75, 7212–7220. [Google Scholar] [CrossRef] [Green Version]
  10. Simon-Delso, N.; Martin, G.S.; Bruneau, E.; Minsart, L.A.; Mouret, C.; Hautier, L. Honeybee colony disorder in crop areas: The role of pesticides and viruses. PLoS ONE 2014, 9, e103073. [Google Scholar] [CrossRef]
  11. Requier, F.; Odoux, J.F.; Henry, M.; Bretagnolle, V. The carry-over effects of pollen shortage decrease the survival of honeybee colonies in farmlands. J. Appl. Ecol. 2017, 54, 1161–1170. [Google Scholar] [CrossRef]
  12. Seitz, N.; Traynor, K.S.; Steinhauer, N.; Rennich, K.; Wilson, M.E.; Ellis, J.D.; Rose, R.; Tarpy, D.R.; Sagili, R.R.; Caron, D.M.; et al. A national survey of managed honey bee 2014–2015 annual colony losses in the USA Nicola. J. Apic. Res. 2015, 54, 292–304. [Google Scholar] [CrossRef]
  13. Le Conte, Y.; Navajas, M. Climate change: Impact on honey bee populations and diseases. Rev. Sci. Tech. 2008, 27, 485–497, 499–510. [Google Scholar] [CrossRef]
  14. Page, R.E.; Peng, C.Y.S. Aging and development in social insects with emphasis on the honey bee, Apis mellifera L. Exp. Gerontol. 2001, 36, 695–711. [Google Scholar] [CrossRef]
  15. Mattila, H.R.; Harris, J.L.; Otis, G.W. Timing of production of winter bees in honey bee (Apis mellifera) colonies. Insectes Soc. 2001, 48, 88–93. [Google Scholar] [CrossRef]
  16. Fluri, P.; Lüscher, M.; Wille, H.; Gerig, L. Changes in weight of the pharyngeal gland and haemolymph titres of juvenile hormone, protein and vitellogenin in worker honey bees. J. Insect Physiol. 1982, 28, 61–68. [Google Scholar] [CrossRef]
  17. Seeley, T.D. The Lives of Bees; Princeton University Press: Princeton, NJ, USA, 2019; ISBN 978-0-691-16676-6. [Google Scholar]
  18. Omholt, S.W.; Amdam, G.V. Epigenetic regulation of aging in honeybee workers. Sci. Aging Knowl. Environ. 2004, 2004, pe28. [Google Scholar] [CrossRef]
  19. Corona, M.; Velarde, R.A.; Remolina, S.; Moran-lauter, A.; Wang, Y.; Hughes, K.A.; Robinson, G.E. Vitellogenin, juvenile hormone, insulin signaling, and queen honey bee longevity. Pnas 2007, 104, 7128–7133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Harwood, G.P.; Ihle, K.E.; Salmela, H.; Amdam, G.V. Regulation of Honeybee Worker (Apis mellifera) Life Histories by Vitellogenin, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2017; Volume 2, ISBN 9780128036082. [Google Scholar]
  21. Kunert, K.; Crailsheim, K. Seasonal changes in carbohydrate, lipid and protein content in emerging worker honeybees and their mortality. J. Apic. Res. 1988, 27, 13–21. [Google Scholar] [CrossRef]
  22. Kunc, M.; Dobeš, P.; Hurychová, J.; Vojtek, L.; Poiani, S.B.; Danihlík, J.; Havlík, J.; Titěra, D.; Hyršl, P. The year of the honey bee (Apis mellifera L.) with respect to its physiology and immunity: A search for biochemical markers of longevity. Insects 2019, 10, 244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kešnerová, L.; Emery, O.; Troilo, M.; Liberti, J.; Erkosar, B.; Engel, P. Gut microbiota structure differs between honeybees in winter and summer. ISME J. 2020, 14, 801–814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Seehuus, S.C.; Krekling, T.; Amdam, G.V. Cellular senescence in honey bee brain is largely independent of chronological age. Exp. Gerontol. 2006, 41, 1117–1125. [Google Scholar] [CrossRef] [Green Version]
  25. Behrends, A.; Scheiner, R. Learning at old age: A study on winter bees. Front. Behav. Neurosci. 2010, 4, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Johnson, C.H.; Ivanisevic, J.; Siuzdak, G. Metabolomics: Beyond biomarkers and towards mechanisms. Nat. Rev. Mol. Cell Biol. 2016, 17, 451–459. [Google Scholar] [CrossRef] [Green Version]
  27. Fan, T.W.M.; Lane, A.N. Applications of NMR spectroscopy to systems biochemistry. Prog. Nucl. Magn. Reson. Spectrosc. 2016, 92–93, 18–53. [Google Scholar] [CrossRef] [Green Version]
  28. Cox, J.E.; Thummel, C.S.; Tennessen, J.M. Metabolomic studies in Drosophila. Genetics 2017, 206, 1169–1185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Lenz, E.M.; Hägele, B.F.; Wilson, I.D.; Simpson, S.J. High resolution 1H NMR spectroscopic studies of the composition of the haemolymph of crowd- and solitary-reared nymphs of the desert locust, Schistocerca gregaria. Insect Biochem. Mol. Biol. 2001, 32, 51–56. [Google Scholar] [CrossRef]
  30. Kapranas, A.; Snart, C.J.P.; Williams, H.; Hardy, I.C.W.; Barrett, D.A. Metabolomics of aging assessed in individual parasitoid wasps. Sci. Rep. 2016, 6, 34848. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, Y.; Carolan, J.C.; Hao, F.; Nicholson, J.K.; Wilkinson, T.L.; Douglas, A.E. Integrated metabonomic-proteomic analysis of an insect-bacterial symbiotic system. J. Proteome Res. 2010, 9, 1257–1267. [Google Scholar] [CrossRef] [Green Version]
  32. Jousse, C.; Dalle, C.; Abila, A.; Traikia, M.; Diogon, M.; Lyan, B.; El Alaoui, H.; Vidau, C.; Delbac, F. A combined LC-MS and NMR approach to reveal metabolic changes in the hemolymph of honeybees infected by the gut parasite Nosema ceranae. J. Invertebr. Pathol. 2020, 176, 1–8. [Google Scholar] [CrossRef]
  33. Dostálková, S.; Dobeš, P.; Kunc, M.; Hurychová, J.; Škrabišová, M.; Petřivalský, M.; Titěra, D.; Havlík, J.; Hyršl, P.; Danihlík, J. Winter honeybee (Apis mellifera) populations show greater potential to induce immune responses than summer populations after immune stimuli. J. Exp. Biol. 2021, 224. [Google Scholar] [CrossRef] [PubMed]
  34. Williams, G.R.; Alaux, C.; Costa, C.; Csáki, T.; Doublet, V.; Eisenhardt, D.; Fries, I.; Kuhn, R.; McMahon, D.P.; Medrzycki, P.; et al. Standard methods for maintaining adult Apis mellifera in cages under in vitro laboratory conditions. J. Apic. Res. 2013, 52, 1–36. [Google Scholar] [CrossRef] [Green Version]
  35. Pang, Z.; Chong, J.; Li, S.; Xia, J. Metaboanalystr 3.0: Toward an optimized workflow for global metabolomics. Metabolites 2020, 10, 186. [Google Scholar] [CrossRef] [PubMed]
  36. Chong, J.; Wishart, D.S.; Xia, J. Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr. Protoc. Bioinforma. 2019, 68, 1–128. [Google Scholar] [CrossRef]
  37. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef] [Green Version]
  38. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016; ISBN 978-3-319-24277-4. [Google Scholar]
  39. Blatt, J.; Roces, F. Haemolymph sugar levels in foraging honeybees (Apis mellifera carnica): Dependence on metabolic rate and in vivo measurement of maximal rates of trehalose synthesis. J. Exp. Biol. 2001, 204, 2709–2716. [Google Scholar] [CrossRef]
  40. Suarez, R.K.; Darveau, C.A.; Welch, K.C.; O’Brien, D.M.; Roubik, D.W.; Hochachka, P.W. Energy metabolism in orchid bee flight muscles: Carbohydrate fuels all. J. Exp. Biol. 2005, 208, 3573–3579. [Google Scholar] [CrossRef] [Green Version]
  41. Arrese, E.L.; Soulages, J.L. Insect Fat Body: Energy, Metabolism, and Regulation. Annu. Rev. Entomol. 2010, 55, 207–225. [Google Scholar] [CrossRef] [Green Version]
  42. Hollis Woodard, S.; Duennes, M.A.; Watrous, K.M.; Jha, S. Diet and nutritional status during early adult life have immediate and persistent effects on queen bumble bees. Conserv. Physiol. 2019, 7, coz048. [Google Scholar] [CrossRef]
  43. Thompson, S.N. Trehalose—The Insect “Blood” Sugar. Adv. Insect Phys. 2003, 31, 205–285. [Google Scholar] [CrossRef]
  44. Koštál, V.; Zahradníčková, H.; Šimek, P.; Zelený, J. Multiple component system of sugars and polyols in the overwintering spruce bark beetle, Ips typographus. J. Insect Physiol. 2007, 53, 580–586. [Google Scholar] [CrossRef] [PubMed]
  45. Overgaard, J.; Malmendal, A.; Sørensen, J.G.; Bundy, J.G.; Loeschcke, V.; Nielsen, N.C.; Holmstrup, M. Metabolomic profiling of rapid cold hardening and cold shock in Drosophila melanogaster. J. Insect Physiol. 2007, 53, 1218–1232. [Google Scholar] [CrossRef]
  46. Watanabe, M.; Kikawada, T.; Minagawa, N.; Yukuhiro, F.; Okuda, T. Mechanism allowing an insect to survive complete dehydration and extreme temperatures. J. Exp. Biol. 2002, 205, 2799–2802. [Google Scholar] [CrossRef] [PubMed]
  47. Thompson, K.M.; Huber, D.P.W.; Murray, B.W. Autumn shifts in cold tolerance metabolites in overwintering adult mountain pine beetles. PLoS ONE 2020, 15, e0227203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Qin, M.; Wang, H.; Liu, Z.; Wang, Y.; Zhang, W.; Xu, B. Changes in cold tolerance during the overwintering period in Apis mellifera ligustica. J. Apic. Res. 2019, 58, 702–709. [Google Scholar] [CrossRef]
  49. Crailsheim, K. Dependence of protein metabolism on age and season in the honeybee (Apis mellifica carnica Pollm). J. Insect Physiol. 1986, 32, 629–634. [Google Scholar] [CrossRef]
  50. Haydak, M.H. Honey Bee Nutrition. Annu. Rev. Entomol. 1970, 15, 143–156. [Google Scholar] [CrossRef]
  51. Crailsheim, K. The protein balance of the honey bee worker. Apidologie 1990, 21, 417–429. [Google Scholar] [CrossRef]
  52. Crailsheim, K.; Leonhard, B. Amino acids in honeybee worker haemolymph. Amino Acids 1997, 13, 141–153. [Google Scholar] [CrossRef]
  53. Hrassnigg, N.; Leonhard, B.; Crailsheim, K. Free amino acids in the haemolymph of honey bee queens (Apis mellifera L.). Amino Acids 2003, 24, 205–212. [Google Scholar] [CrossRef]
  54. Storey, K.B.; Storey, J.M. Natural freezing survival in animals. Annu. Rev. Ecol. Syst. 1996, 27, 365–386. [Google Scholar] [CrossRef] [Green Version]
  55. Koštál, V.; Šimek, P.; Zahradníčková, H.; Cimlová, J.; Štětina, T. Conversion of the chill susceptible fruit fly larva (Drosophila melanogaster) to a freeze tolerant organism. Proc. Natl. Acad. Sci. USA 2012, 109, 3270–3274. [Google Scholar] [CrossRef] [Green Version]
  56. Koštál, V.; Korbelová, J.; Rozsypal, J.; Zahradníčková, H.; Cimlová, J.; Tomčala, A.; Šimek, P. Long-term cold acclimation extends survival time at 0°C and modifies the metabolomic profiles of the larvae of the fruit fly drosophila melanogaster. PLoS ONE 2011, 6, e25025. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Williams, C.M.; McCue, M.D.; Sunny, N.E.; Szejner-Sigal, A.; Morgan, T.J.; Allison, D.B.; Hahn, D.A. Cold adaptation increases rates of nutrient flow and metabolic plasticity during cold exposure in Drosophila melanogaster. Proc. R. Soc. B Biol. Sci. 2016, 283. [Google Scholar] [CrossRef] [Green Version]
  58. Mattila, H.R.; Otis, G.W. Dwindling pollen resources trigger the transition to broodless populations of long-lived honeybees each autumn. Ecol. Entomol. 2007, 32, 496–505. [Google Scholar] [CrossRef]
  59. Le Couteur, D.G.; Solon-Biet, S.; Cogger, V.C.; Mitchell, S.J.; Senior, A.; De Cabo, R.; Raubenheimer, D.; Simpson, S.J. The impact of low-protein high-carbohydrate diets on aging and lifespan. Cell. Mol. Life Sci. 2016, 73, 1237–1252. [Google Scholar] [CrossRef] [PubMed]
  60. Paoli, P.P.; Wakeling, L.A.; Wright, G.A.; Ford, D. The dietary proportion of essential amino acids and Sir2 influence lifespan in the honeybee. Age 2014, 36, 1239–1247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Gómez-Moracho, T.; Durand, T.; Pasquaretta, C.; Heeb, P.; Lihoreau, M. Artificial diets modulate infection rates by nosema ceranae in bumblebees. Microorganisms 2021, 9, 158. [Google Scholar] [CrossRef] [PubMed]
  62. Amdam, G.V.; Simões, Z.L.P.; Hagen, A.; Norberg, K.; Schrøder, K.; Mikkelsen, Ø.; Kirkwood, T.B.L.; Omholt, S.W. Hormonal control of the yolk precursor vitellogenin regulates immune function and longevity in honeybees. Exp. Gerontol. 2004, 39, 767–773. [Google Scholar] [CrossRef]
  63. Imboden, H.; König, R.; Ott, P.; Lustig, A.; Kämpfer, U.; Lanzrein, B. Characterization of the native vitellogenin and vitellin of the cockroach, Nauphoeta cinerea, and comparison with other species. Insect Biochem. 1987, 17, 353–365. [Google Scholar] [CrossRef]
  64. Wheeler, D.E.; Kawooya, J.K. Purification and characterization of honey bee vitellogenin. Arch. Insect Biochem. Physiol. 1990, 14, 253–267. [Google Scholar] [CrossRef] [PubMed]
  65. Moessinger, C.; Klizaite, K.; Steinhagen, A.; Philippou-Massier, J.; Shevchenko, A.; Hoch, M.; Ejsing, C.S.; Thiele, C. Two different pathways of phosphatidylcholine synthesis, the Kennedy Pathway and the Lands Cycle, differentially regulate cellular triacylglycerol storage. BMC Cell Biol. 2014, 15, 43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Hardie, D.G. AMP-activated/SNF1 protein kinases: Conserved guardians of cellular energy. Nat. Rev. Mol. Cell Biol. 2007, 8, 774–785. [Google Scholar] [CrossRef] [PubMed]
  67. Bonkowski, M.S.; Sinclair, D.A. Slowing ageing by design: The rise of NAD+ and sirtuin-activating compounds. Nat. Rev. Mol. Cell Biol. 2016, 17, 679–690. [Google Scholar] [CrossRef]
  68. Apfeld, J.; O’Connor, G.; McDonagh, T.; DiStefano, P.S.; Curtis, R. The AMP-activated protein kinase AAK-2 links energy levels and insulin-like signals to lifespan in C. elegans. Genes Dev. 2004, 18, 3004–3009. [Google Scholar] [CrossRef] [Green Version]
  69. Funakoshi, M.; Tsuda, M.; Muramatsu, K.; Hatsuda, H.; Morishita, S.; Aigaki, T. A gain-of-function screen identifies wdb and lkb1 as lifespan-extending genes in Drosophila. Biochem. Biophys. Res. Commun. 2011, 405, 667–672. [Google Scholar] [CrossRef] [PubMed]
  70. Chuang, Y.L.; Hsu, C.Y. Changes in mitochondrial energy utilization in young and old worker honeybees (Apis mellifera). Age 2013, 35, 1867–1879. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Hsu, C.Y.; Chuang, Y.L. Changes in energy-regulated molecules in the trophocytes and fat cells of young and old worker honeybees (Apis mellifera). J. Gerontol. Ser. A Biol. Sci. Med. Sci. 2014, 69, 955–964. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Typical proton nuclear magnetic resonance (1H NMR) spectrum of honey bee worker body (gut and venom sac excluded). The numbers correspond to these quantified compounds: 1. Adenosine monophosphate (AMP), 2. Oxidized nicotinamide adenine dinucleotide (NAD+), 3. Inosine, 4. Adenine, 5. Histidine, 6. Uridine, 7. Phenylalanine, 8. Sucrose, 9. Glucose, 10. Trehalose, 11. sn-glycero-3-phosphocholine, 12. Fructose, 13. Proline, 14. O-phosphocholine, 15. Taurine, 16. Trimethylamine N-oxide, 17. Choline, 18. β-alanine, 19. Dimethyl sulfone, 20. Malonate, 21. Putrescine, 22. Creatine, 23. Lysine, 24. Asparagine, 25. Unknown, 26. Sarcosine, 27. Glutamine, 28. Succinate, 29. Suberate, 30. Acetate, 31. Alanine, 32. Threonine, 33. Valine, 34. Isoleucine, 35. Leucine, 36. Glycine, and 37. Methanol.
Figure 1. Typical proton nuclear magnetic resonance (1H NMR) spectrum of honey bee worker body (gut and venom sac excluded). The numbers correspond to these quantified compounds: 1. Adenosine monophosphate (AMP), 2. Oxidized nicotinamide adenine dinucleotide (NAD+), 3. Inosine, 4. Adenine, 5. Histidine, 6. Uridine, 7. Phenylalanine, 8. Sucrose, 9. Glucose, 10. Trehalose, 11. sn-glycero-3-phosphocholine, 12. Fructose, 13. Proline, 14. O-phosphocholine, 15. Taurine, 16. Trimethylamine N-oxide, 17. Choline, 18. β-alanine, 19. Dimethyl sulfone, 20. Malonate, 21. Putrescine, 22. Creatine, 23. Lysine, 24. Asparagine, 25. Unknown, 26. Sarcosine, 27. Glutamine, 28. Succinate, 29. Suberate, 30. Acetate, 31. Alanine, 32. Threonine, 33. Valine, 34. Isoleucine, 35. Leucine, 36. Glycine, and 37. Methanol.
Insects 13 00193 g001
Figure 2. Principal component analysis of summer and winter bees (a) Score plot. Dim1(PC1) on the x-axis explains 42.6%, Dim2 (PC2) on the y-axis explains 12.1% of variability between the summer bees and winter bees; (b) Variable correlation plot. The length and colors of the arrows explain how each variable expresses the differences between the two compared groups (summer and winter). The longer the arrow, the stronger its influences in discrimination. The red and orange colors show a high contribution on each principal component, whereas the green and blue show a relatively low contribution.
Figure 2. Principal component analysis of summer and winter bees (a) Score plot. Dim1(PC1) on the x-axis explains 42.6%, Dim2 (PC2) on the y-axis explains 12.1% of variability between the summer bees and winter bees; (b) Variable correlation plot. The length and colors of the arrows explain how each variable expresses the differences between the two compared groups (summer and winter). The longer the arrow, the stronger its influences in discrimination. The red and orange colors show a high contribution on each principal component, whereas the green and blue show a relatively low contribution.
Insects 13 00193 g002
Figure 3. Comparison of 36 metabolites between summer bees and winter bees: (a) Carbohydrates; (b) Amino acids; (c) Choline-containing compounds; (d) Nucleoside, nucleotides, and nucleobases; (e) Others. X-axis: the groups summer (red) and winter (green), y-axis: metabolite concentration (μg/bee), black bar: mean. The P-values were obtained from the independent t-test. Asterisk indicates significant differences * p < 0.05, ** p < 0.005, *** p < 0.0001, ns = not significant. Benjamini–Hochberg’s procedure (false discovery rate, FDR 0.05) is applied after the independent t-test. Details of the independent t-test and Benjamini–Hochberg’s procedure informed in Table 1.
Figure 3. Comparison of 36 metabolites between summer bees and winter bees: (a) Carbohydrates; (b) Amino acids; (c) Choline-containing compounds; (d) Nucleoside, nucleotides, and nucleobases; (e) Others. X-axis: the groups summer (red) and winter (green), y-axis: metabolite concentration (μg/bee), black bar: mean. The P-values were obtained from the independent t-test. Asterisk indicates significant differences * p < 0.05, ** p < 0.005, *** p < 0.0001, ns = not significant. Benjamini–Hochberg’s procedure (false discovery rate, FDR 0.05) is applied after the independent t-test. Details of the independent t-test and Benjamini–Hochberg’s procedure informed in Table 1.
Insects 13 00193 g003
Table 1. Mean (μg/bee), standard deviation (SD), independent t-test, and Benjamini–Hochberg’s procedure (FDR 0.05).
Table 1. Mean (μg/bee), standard deviation (SD), independent t-test, and Benjamini–Hochberg’s procedure (FDR 0.05).
Summber Bees
(n = 75)
Winter Bees
(n = 45)
Independent
t-Test
Benjamini–Hochberg
MetaboliteMeanSDMeanSDDF 1p-Valuep-Value
AMP27.566.6036.127.551180.0019<0.0001
Acetate11.298.427.712.69118<0.00010.0027
Adenine1.290.820.460.25118<0.0001<0.0001
Adenosine1.250.682.761.581180.0001<0.0001
Alanine23.388.9028.566.72118 [115.6]<0.00010.0001
Asparagine10.453.555.522.09118<0.0001<0.0001
Choline27.359.6612.573.50118<0.0001<0.0001
Creatine2.720.713.080.721180.01200.0166
Fructose509.4611.61109.11008.3118 [84.2]0.00010.0001
Glucose655.8442.2915.0715.6118 [68.8]0.16640.1997
Glutamine74.2518.2960.9314.421180.00010.0002
Glycine11.373.1213.5211.001180.36080.4059
Histidine7.821.857.751.731180.84490.8449
Inosine14.734.769.172.36118<0.0001<0.0001
Isoleucine9.824.416.061.77118 [112.8]<0.0001<0.0001
Leucine10.356.023.541.31118<0.0001<0.0001
Lysine10.834.456.672.59118<0.0001<0.0001
Malonate6.680.936.460.971180.23470.2726
NAD+4.541.788.662.48118 [116.8]<0.0001<0.0001
O-Phosphocholine68.3817.7588.0115.99118 [117.8]<0.0001<0.0001
Phenylalanine4.392.531.310.55118<0.0001<0.0001
Proline99.2429.01123.4729.05118<0.0001<0.0001
Putrescine7.412.627.281.75118 [116.6]0.75550.7771
Sarcosine4.001.804.550.97118 [117.5]0.00080.0012
Suberate5.102.904.582.191180.59160.7204
Succinate3.861.535.492.471180.6804<0.0001
Sucrose284.3301.91690.32135.0118<0.0001<0.0001
Taurine87.5716.4779.6218.12118<0.00010.0167
Threonine10.194.736.002.971180.0125<0.0001
Trehalose416.0181.1801.8212.0118 [115.7]<0.0001<0.0001
Trimethylamine
N-oxide
0.860.632.971.54118<0.0001<0.0001
Unknown19.617.621.530.55118<0.0001<0.0001
Uridine2.380.852.010.59118<0.00010.0499
Valine8.833.603.611.521180.0388<0.0001
sn-Glycero-3-phosphocholine61.2821.2863.4218.68118<0.00010.6454
β-Alanine33.816.8031.567.311180.07610.0945
1 Correction to the degree of freedom (DF) has applied when the equal variance is not assumed and the corrected value is shown in brackets.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, S.; Kalcic, F.; Duarte, I.F.; Titera, D.; Kamler, M.; Mrna, P.; Hyrsl, P.; Danihlik, J.; Dobes, P.; Kunc, M.; et al. 1H NMR Profiling of Honey Bee Bodies Revealed Metabolic Differences between Summer and Winter Bees. Insects 2022, 13, 193. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13020193

AMA Style

Lee S, Kalcic F, Duarte IF, Titera D, Kamler M, Mrna P, Hyrsl P, Danihlik J, Dobes P, Kunc M, et al. 1H NMR Profiling of Honey Bee Bodies Revealed Metabolic Differences between Summer and Winter Bees. Insects. 2022; 13(2):193. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13020193

Chicago/Turabian Style

Lee, Saetbyeol, Filip Kalcic, Iola F. Duarte, Dalibor Titera, Martin Kamler, Pavel Mrna, Pavel Hyrsl, Jiri Danihlik, Pavel Dobes, Martin Kunc, and et al. 2022. "1H NMR Profiling of Honey Bee Bodies Revealed Metabolic Differences between Summer and Winter Bees" Insects 13, no. 2: 193. https://0-doi-org.brum.beds.ac.uk/10.3390/insects13020193

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop