Honeybees share a mutual relationship with the plants in their habitats. Plants rely on the pollination services offered by honeybees for their reproduction and survival. Honeybees, in turn, accrue nutritional benefits from plants (de Sá-Otero et al., 2009). The bees’ nutrition is composed of plant-derived nectar and pollen (Roulston & Cane, 2000; Nicolson & Thornburg, 2007). Whereas nectar acts as the main carbohydrate source, pollen is the primary protein source for developing larvae and adult honeybees (Brodschneider & Crailsheim, 2010). The nutritional content of pollen varies with the host-plant species with some pollen species being nutritionally richer than others (Maurizio & Louveaux, 1965; Odoux et al., 2012), with noticeable differences in their protein, amino acid, lipid, micronutrient, and sugar content (Keller et al., 2005; Roulston & Cane, 2000).
The nutritional composition of pollen influences the foraging preferences of bees (de Sá-Otero et al., 2009). Honeybees are generalists collecting a diversity of polyfloral pollen constituting a balanced and optimal diet (Schmidt, 1984; Schmidt et al., 1987). Intensified agricultural practices and urbanization are contributing to the loss of pollen sources for honeybees. Consequently, rich habitats endowed with polyfloral pollen sources are slowly being altered towards habitats with unifloral (for monocultures) or few floral sources (Decourtye et al., 2010; Naug, 2009). This limits the quality and feeding choice of the honeybee, a trend also being witnessed in sub-Saharan Africa (Nkonya et al., 2016).
Multiple studies have already been conducted to investigate the effect of both a pollen and a no-pollen diet on the honeybee (Schmidt, 1984; Schmidt et al., 1987; Pernal & Currie 2000). Furthermore, there have been extensive studies on the effects of polyfloral versus monofloral pollen diet on the European honeybee subspecies in the northern Hemisphere (Schmidt, 1984; Schmidt et al., 1987; Danihlík et al., 2018; Dolezal et al., 2019), but to the best of our knowledge, this is little understood in African honeybees, including the widely distributed sub-species of
This study aimed to answer the following research questions: (i) What are the survival rates of honeybees fed on lowly diverse relative to highly diverse pollen diets?, (ii) How do different pollen diets affect the body weight and pollen consumption rates of honeybees?, (iii) Does variation in pollen diversity influence the immune gene responses in honeybees? To answer these questions, newly emerged honeybees were fed with different pollen regimens, and measured parameters related to the physiology as weight, immune response and survival were assessed. We hypothesized that the honeybees fed on highly diverse pollen would show increased survival and immunity. This study is important in understanding the dietary requirements of the African honeybee and its physiological responses to varying diversities of pollen diets.
This study was conducted at the International Centre of Insect Physiology and Ecology (
Age-matched honeybees were obtained using sealed brood combs containing late-stage pupae which were incubated overnight at 34±2°C and 50–70% relative humidity (Williams et al., 2013). The newly emerged workers (less than 24 h old) were transferred into cages and kept in an incubator at 34±2°C and 50–70% relative humidity. We used metal cages measuring 10 cm × 8.5 cm × 5 cm and circular plastic cages measuring 12 cm diameter × 8 cm height accompanied by their feeders to confine cohorts of five and 50 honeybees, respectively.
Adult honeybees were fed
We used commercially available corbicular pollen (Aurica Naturheilmittel und Naturwaren GmbH, Schwalbach-Elm, Germany, originally collected from Spain) because locally available pollen was insufficient for the experimental setups. The pollen loads were hand-sorted based on color and appearance. Further on, using microscopy, we examined the hand-sorted pollen loads to verify whether they were also morphologically distinct. We obtained a total of ten pollen morphotypes (Tab. S2). The pollen types were later constituted into two pollen diets: LD pollen consisting of a single pollen type and HD pollen comprising of ten pollen types at equal frequency. The pollen diets were immediately frozen at −20°C until use.
We analyzed the protein content of the pollen diets using the Bradford method with minor modifications (Bradford, 1976). Pollen was sampled from the daily food supplied to the bees kept in the metallic cages (n=5 bees). Each pollen sample was crushed into a fine powder using a pestle and mortar, transferred to a micro-centrifuge tube, and weighed on (Analytical Plus, Ohaus, Switzerland) a weighing balance (accuracy ±0.0001 g). Applying the method used by de Sá-Otero et al., (2009) with slight modifications (Ochungo et al., 2021), we extracted protein from the samples. Briefly, 5 ml of 30 mM TE-buffer was added to each sample, vortexed, and centrifuged for 10 min at 3,000 × g. We collected the supernatant in aliquots of 0.3 mL, added 4 mL of Bradford reagent and incubated the reaction mixture for 2 min at room temperature. This procedure was done in triplicates.
We used 0.3 mL of TE buffer as blank and Bovine Serum Albumin (BSA) (Sigma-Aldrich-Kobian, Kenya) as the standard (concentration 2 mg/mL). Then, to the various BSA dilutions prepared: 50 μg/mL, 150 μg/mL, 200 μg/mL, 250 μg/mL, and 400 μg/mL, we added 4 mL of Bradford reagent. Protein absorbances were measured at 595 nm using a Bio Spec-mini spectrophotometer (Shimadzu Corporation, Japan). We used absorbance values of the standards to generate a calibration curve and its corresponding linear formula. From this formula, we obtained crude protein values (μg protein/mL of pollen) and afterwards converted them to g protein/100 g pollen. We used a total of 19 and 16 replicates for HD and LD diets respectively.
We marked honeybees kept in small metallic cages in groups of five bees per cage with different paint marks on their thoraces and weighed them to record the initial weight of individual bees using a weighing balance (Analytical Plus, Ohaus, Switzerland) (accuracy ±0.0001 g). The honeybees were subjected to different diet regimens (LD and HD). We recorded the weight of the individual honeybees daily before supplying a fresh diet into the cages. To do this, we placed individual bees in a 1.5 mL Eppendorf tube and subtracted the weight of the empty tube from the total weight of the tube + weighed bee. Previous studies had measured the weight of whole cages (Harbo & Hoopingarner, 1997; Harbo & Harris, 1999), but we were interested in the individual body weight of honeybees. Since the bees were colour-marked, we recorded the weight of individual bees in all cages daily. We further estimated the weight of the daily pollen consumption per cage after collecting the remaining diet the next day. This experimental setup was replicated four times. To ascertain the effect of pollen diets on survival, we kept groups of 50 honeybees in plastic cages. Daily inspection of the bees was conducted to remove and score the number of dead honeybees as well as renew the food and water supply. We terminated the experimental setups when all honeybees in the cages died. We replicated the experiment four times.
For the immune gene expression assay, we set up cages each consisting of 50 freshly emerged honeybee workers. Subsequently, the bees received either of the two pollen treatments for six days and this experiment was conducted in triplicates. We changed the food and removed dead bees daily. On day six, we harvested ten bees per cage by freeze-killing at −80°C. These samples remained frozen until further processing.
Total RNA was isolated from individual bees (n=30 bees for each treatment group) in order to quantify the expression of the target immune gene. This involved grounding the honeybee abdominal tissues in liquid nitrogen for RNA isolation using TRIzol (Invitrogen, USA) as described by Evans et al. (2013). RNA samples were quantified based on absorbance at 260 nm using a NanodropTM 2000 spectrophotometer (Thermofisher Scientific, USA) and quality was inferred from 260/280 nm ratio as well as from 230/260 nm ratio. We used 10 μL total RNA to synthesize cDNA using the High-Capacity cDNA Reverse Transcription Kit that contained oligo-dT primers (Applied Biosystems, Lithuania) according to the manufacturer’s instruction. The thermal conditions for the reaction were 25°C for 10 min, 37°C for 120 min, and 85°C for 5 min on T100TM Thermal Cycler (Bio-Rad, Singapore). We pooled equal volumes of cDNA (10 μL) from two bees per treatment group as a template for qRT-PCR reactions.
We conducted quantitative real-time PCR (qRT-PCR) on 96-well microtitre plates using HOT FIREPol EvaGreen® qPCR Mix Plus with ROX (Solis Biodyne, Estonia) and gene-specific primers (Tab. S1) for
We tested all data for normality using the Shapiro-Wilks test. Since all our data were non-normally distributed, we used non-parametric tests for analysis. To compare differences in protein content between the LD and HD pollen diets, we used the Mann-Whitney-U-Test. We used the Cox proportional-hazard regression analysis as implemented in the
Tata Binding Protein (TBP) served as a housekeeping gene (HKG), and
We performed all our statistical analyses using R software 3.6.2 version (R Core Team 2020) with appropriate packages:
The protein content of the different pollen diets was similar (LD pollen: 5.3±4.4 mg/100 mg, HD pollen: 5.6±7.3 mg/100 mg; Mann-Whitney U-test: w=210, p=0.52; Fig. S1).
HD pollen diet significantly increased the lifespan of the bees relative to those supplied with the LD pollen diet (Cox hazard proportion test; Cox log-rank score=10.56; beta HR=0.34; df=1; p=0.001; Fig. 1).
HD pollen consumption per day ranged between 0.3 mg - 14 mg, while daily LD pollen consumption ranged between 0.4 mg-11.5 mg. HD-fed bees showed a higher pollen consumption relative to LD-fed bees across the days (Fig. S2a). The null model only including random effects (cage and replicate) was better than our full model (GLMM: LRT: 0.4532 p=0.32; AIC: null model: −338.74; full model: −335.57) implying that the fixed effects, feeding day and diet, and their interaction had no significant effect on daily pollen consumption.
The range in individual bodyweight of HD-fed bees was 59.6 mg - 101.7 mg and LD-fed bees 48.9 mg - 109.5 mg. LD and HD-fed bees shared a common trend from day 1–3. After day 3, the HD-reared bees experienced a decrease in weight persisting to day 4 before slightly increasing towards day 5 (Fig. S2b). The LD-reared bees recorded an onward sharp increase in weight past day 3 up to day 5 (Fig. S2b). Bee weight was affected by the feeding day (GLMM: df=7, t=−4.948, p<0.0001) and not by the type of diet (GLMM: df=7, t=0.972, p=0.33) or interaction between the two factors (GLMM; df=7, t=−1.720, p=0.09). HD pollen consumption was correlated to the bee weight gain as marginally significant (r=0.9, p=0.08), while this effect was less strong in LD (r=0.7, p=0.23) (Fig. 2).
LD pollen diet increased the expression levels of
Diet comparisons can be investigated at two levels (i) pollen vs. no-pollen diet, and (ii) HD vs. LD pollen diet. Previous studies have already investigated the effect of the 1st level of comparison on honeybee health (Schmidt, 1984; Schmidt et al., 1987; Pernal & Currie 2000). We instead focused on the influence of the 2nd level of comparison, the quality of the pollen diet, by comparing the HD vs. LD pollen diet on bee health. We attempted to answer whether pollen quality sourced from two environments, namely lowly diverse habitat and rich biodiverse habitat, influences honeybee development and wellbeing. Our study, therefore, highlighted the effects of these two different pollen diets on the honeybee-related parameters of survival, bee body weight change, pollen consumption and immune response.
The LD pollen consisted of one predominant pollen based on colour, while the HD pollen was a blend of ten equally frequent pollen morphotypes verified by microscopic examination (Tab. S2). Pollen types can share the same colour shade yet belong to different plant species owing to other morphological differences in shape, size, and weight (Almeida-Muradian et al., 2005; Komosinska-Vassev et al., 2015).
The protein values for both pollen types were, on average, within the acceptable range of 2.5%–61% (Roulston et al., 2000). However, we attribute the low amount of crude protein values to the stored, dried pollen that was used in this study. Drying pollen significantly reduces the crude protein content of pollen (Human & Nicolson, 2006). LD-fed bees received the same quantity of protein as HD-fed bees. Although the method used, Bradford assay, provides a general guideline on the protein quality, it is limited to determining the number of specific amino acids (Pernal & Currie, 2000). Therefore, we propose that any difference elicited in the phenotype of the bees is due to differing amounts of amino acids and fatty acids in the pollen. Previously, different pollen types had been reported to share a similar concentration of proteins yet differ in their amino acids and fatty acids concentration (DeGrandi-Hoffman et al., 2018). This could have drastic implications for either a sufficiency or scarcity of essential amino acids such as tryptophan and phenylalanine (Roulston & Cane, 2000).
The longer lifespan registered by the HD-fed bees could be attributed to the nutritional differences in the two pollen regimens apart from the protein content. Pollen also contains fatty acids and micronutrients in varying amounts and composition in different pollen types (Keller et al., 2005; Di Pasquale et al., 2013). Polyfloral pollen is composed of pollen obtained from different floral sources, which collectively influences honeybee nutrition (Donkersley et al., 2017). Our findings coincide with previous reports which demonstrated that caged bees fed on polyfloral pollen survived longer than their counterparts on monofloral pollen (Schmidt et al., 1987; Di Pasquale et al., 2013; Dolezal et al., 2019). A polyfloral diet mimics the natural pollen diet used by honeybees (Schmidt et al., 1987). Honeybees are generalists that collect diverse multifloral pollen from different plant species to constitute a balanced and diverse diet (Dimou & Thrasyvoulou, 2009). A high protein diet seems to boost the honeybee’s survival even under stressful conditions (Archer et al., 2014). In field studies, the quality and availability of pollen have been shown to influence colony survival and productivity. Both, low amounts of pollen food stores and pollen of low nutritional value makes brood-rearing difficult (Kleinschmidt & Kondos, 1976; Brodschneider & Crailsheim, 2010) and might even lead to absconding (Cheruiyot et al., 2020). A noticeable increase in pollen consumption in HD- and LD-fed bees in the first three days was consistent with other studies (Alqarni, 2006). Possibly, this increased appetite was to meet the nutritional demands of the honeybees during this early, critical phase of growth and development. Pollen diet is crucial for the first seven days of adult honeybee development (Brodschneider & Crailsheim, 2010). In addition, it promotes hypopharyngeal gland development in brood and young bees (Di Pasquale et al., 2013) and determines the lifespan of honeybee workers (Amdam & Omholt, 2002).
The honeybees fed under the LD pollen regimen consumed a lower amount of pollen. The pollen consumed is a reliable measure of the palatability of the diet (Pernal & Currie, 2000). We propose that this reduced palatability in the LD diet could be associated with its poor unbalanced nutritional state. LD-fed bees would be expected to ingest more pollen to compensate for any lacking nutrients. However, honeybees do not ingest more of a poor pollen diet to compensate for any nutritional deficiencies in their diet or counter stressful conditions (Knox et al., 1971; Archer et al., 2014). If they are forcefully fed, this results in lowered survival and compromises health status (Paoli et al., 2014). Under field conditions, whether bees evaluate which pollen to consume is presently unclear. While foragers have been shown to select nutritionally dense pollen (Hendriksma & Shafir, 2016), nurse bees do not evaluate their choice of pollen diet based on nutrition (Corby-Harris et al., 2018). Noteworthy is that besides the nutritional quality of pollen, the amount of pollen consumed is determined by physical and chemical cues (Schmidt & Johnson, 1984), including phagostimulants which occur in different combinations for different pollen types (Schmidt, 1985).
LD-fed bees increased in weight over time. However, there is a relatively weaker correlation of change in body weight to the pollen they consumed relative to the HD-fed bees. The higher weight recorded could not be traced only to pollen consumed. We suggest that the LD-fed bees reverted to ingesting more of the additional sugar solution supplied possibly to compensate for the nutritional deficiency in the LD diet. Previously, honeybees under laboratory conditions had been reported to regulate their nutrient intake by shifting to a carbohydrate-biased diet (Altaye et al., 2010). However, sugar-solution consumption was not quantified in this study.
We used a low number of bees (n=5) because individual weighing generally induces stress in bees (Sgolastra et al., 2017); a previous study had used a similar number of bees (Arien et al., 2018). We measured the pollen consumption and bee bodyweight up to day five as some cages had no surviving bees on day 6. Eusocial honeybees endure stress when kept in small numbers resulting in a shorter survival period (Rinderer & Baxter, 1978). Furthermore, in ants, isolated ones (n=1 or 2) were more hyperactive and increased their energy demand compared to grouped ants (n=10), thus leading to higher mortality (Koto et al., 2015). In bumblebees, isolated bees showed a reduced immune response compared to bees kept in groups (Richter et al., 2012).
In summary, although HD and LD pollen diets appear to be similar in their protein quantity, there were detectable differences in their influence on the survival and physiology of the caged honeybees. Our findings reveal that the HD pollen diet was consumed more than the LD pollen diet and promoted higher survivorship as we had hypothesized. Besides proteins, pollen is packed with fatty acids and lipids alongside micronutrients which vary with the pollen type.
Our HD pollen diet, polyfloral in its composition, possibly supplied macro- and micronutrients lacking in the LD pollen diet. Our pollen consumption results reveal that LD-fed bees did not ingest more pollen than the HD-fed bees to compensate for any nutritional deficiencies. Finally, we noted a stronger correlation of bee weight and pollen consumption in the HD-fed bees than in the LD-fed bees, which implies a highly diverse pollen diet has a stronger influence on honeybee physiology. Further investigation is required to determine the effect of individual pollen morphotypes on the survival and physiology of the African bee.
Our study conducted under laboratory conditions mimics two contrasting habitats and their possible effect on honeybee life-history traits and physiology. The findings support the concept that honeybees in highly degraded habitats are nutritionally deprived of an optimum diet unlike those in diverse and intact environments. Under malnourished conditions, honeybees shift their feeding patterns besides being more susceptible to other environmental stressors acting synergistically with malnutrition. Thus, existing diverse environments need to be conserved as a prerequisite for securing the health of the honeybee and other pollinators to ensure optimum delivery of the ecosystem service of pollination.