Honeybee Colony Thermoregulation – Regulatory Mechanisms and Contribution of Individuals in Dependence on Age, Location and Thermal Stress

Abstract

Honeybee larvae and pupae are extremely stenothermic, i.e. they strongly depend on accurate regulation of brood nest temperature for proper development (33–36°C). Here we study the mechanisms of social thermoregulation of honeybee colonies under changing environmental temperatures concerning the contribution of individuals to colony temperature homeostasis. Beside migration activity within the nest, the main active process is “endothermy on demand” of adults. An increase of cold stress (cooling of the colony) increases the intensity of heat production with thoracic flight muscles and the number of endothermic individuals, especially in the brood nest. As endothermy means hard work for bees, this eases much burden of nestmates which can stay ectothermic. Concerning the active reaction to cold stress by endothermy, age polyethism is reduced to only two physiologically predetermined task divisions, 0 to ∼2 days and older. Endothermic heat production is the job of bees older than about two days. They are all similarly engaged in active heat production both in intensity and frequency. Their active heat production has an important reinforcement effect on passive heat production of the many ectothermic bees and of the brood. Ectothermy is most frequent in young bees (<∼2 days) both outside and inside of brood nest cells. We suggest young bees visit warm brood nest cells not only to clean them but also to speed up flight muscle development for proper endothermy and foraging later in their life. Young bees inside brood nest cells mostly receive heat from the surrounding cell wall during cold stress, whereas older bees predominantly transfer heat from the thorax to the cell wall. Endothermic bees regulate brood comb temperature more accurately than local air temperature. They apply the heat as close to the brood as possible: workers heating cells from within have a higher probability of endothermy than those on the comb surface. The findings show that thermal homeostasis of honeybee colonies is achieved by a combination of active and passive processes. The differential individual endothermic and behavioral reactions sum up to an integrated action of the honeybee colony as a superorganism.

About the Authors:

Anton Stabentheiner

* E-mail: [email protected] (AS); [email protected] (HK)

Affiliation: Institut für Zoologie, Karl-Franzens-Universität Graz, Graz, Austria

Helmut Kovac

* E-mail: [email protected] (AS); [email protected] (HK)

Affiliation: Institut für Zoologie, Karl-Franzens-Universität Graz, Graz, Austria

Robert Brodschneider

Affiliation: Institut für Zoologie, Karl-Franzens-Universität Graz, Graz, Austria

Introduction

Among social insects, the cavity nesting honeybee species Apis mellifera and A. cerana display the most advanced regulation of the nest climate [1]. Of these two species, Apis mellifera is the more intensely investigated concerning thermal homeostasis of the colony. The endothermic heat production and the insulation of the breeding cavity allows the bees to regulate the brood nest temperature within the range of 32–36°C [e.g.], [ 2]–[10], or to survive cold winters [11], [12]. As the honeybee brood is extremely stenothermic, accurate temperature regulation is indispensable for its proper development [4], [13]. While eggs and larvae (in open brood cells) can tolerate lower temperatures for some time, the pupae (in sealed brood cells) are very sensitive to cooling. If they remain too long below 32°C there is a high incidence of shrivelled wings and legs, and malformations of the abdomen [4], and adults may suffer from neural and behavioral insufficiencies [14]–[17]. Accordingly, the accuracy of thermoregulation is high in the presence of brood [3], [4], [9], [10], and much more variable and generally lower in broodless colonies [8], [18].

The brood itself has a low metabolic rate at young ages [19]–[21]. As it lacks regulatory ability and does not provide enough heat on its own, it would not be able to achieve thermal constancy in a variable environment. Warming of the brood has therefore to be accomplished by the worker bees. This warming behavior is released by chemical stimuli and physical properties of the brood. Sealed cells are more attractive than open ones [13], [20]. If the hive is in danger of being overheated the bees cool it by fanning, and collected water is spread on the combs [20], [22]–[24].

Much research on the thermoregulation of breeding honeybee colonies concentrated on the colony as a whole “superorganism” [25], [26], describing aspects of the hive microclimate like the temperature of the brood nest and the combs, the air temperature and/or humidity between them, the temperature gradient across the colony at various ambient conditions, and diurnal rhythms of hive temperature and metabolic rate [e.g. 3]–[5], [8], [18], [20], [25], [27]. Investigations of this type, however, neglect the fact that thermal homeostasis of a honeybee colony is the result of the cooperation of thousands of individuals in a multi-factor regulatory system. On the individual level there is profound knowledge on the respiration and thus heat production of the eggs [28], the larvae [19], [21], and of resting and active bees of different ages [29]–[31]. Principles of the individual brood warming behavior by active heat production (endothermy) via the thoracic flight muscles were investigated in bees heating brood cells from outside [9], [32]–[34], and in bees heating the brood from within cells [10]. However, any comprehensive approach of the mechanisms of colony thermoregulation has to include the contribution of ectothermic bees to colonial heat production. In order to understand the relative importance of the regulatory mechanisms, parallel investigation of as many variables as possible is necessary. We here provide quantitative experimental data of honeybee body temperature in reaction to environmental temperature changes. This allows detailed analysis of the regulatory mechanisms of colony temperature homeostasis driving dynamics and relative frequency of endothermy (with active heat production) and ectothermy (no heat production beyond the resting level) of individuals in different hive locations. The use of infrared thermography made possible thermal investigations on a high number of bees without behavioral impairment.

Honeybee workers were reported to heat brood combs also from within cells [10], [35]. However, the importance of this behavior in relation to endothermic activity of bees on combs is unknown. We here not only analyze the frequency of cell visitation in dependence on changes of environmental temperature but also present a comparison of the frequency of endothermy between bees in cells and on combs. Special emphasis is given to the analysis of the direction of heat flow in cells visited by bees.

Temporal polyethism (age-related division of labor) within honeybee colonies usually addresses the visually observable activities like cell cleaning, nursing, food processing, guarding and foraging [see e.g. 36]–[43]. Harrison [32] was the first to investigate the contribution of individual workers of different age to another important task, colonial heat production. Applying a high cold stress to single combs, he was not able to find age-dependent differences in the intensity of endothermy of hive bees older than 2 days. Only 1-day-old bees displayed less intense endothermy. However, the age distribution of endothermic activity might well be different at less extreme temperatures. We here report extended measurements on this topic, including very young bees (age <24 h), with emphasis of the relative abundance of both endothermic and ectothermic bees in various age classes from emergence to the foraging age.

Materials and Methods

Colonies, Bee Treatment and Experimental Procedure

Experiments were carried out in observation hives with two vertically aligned combs which were covered by plastic films transparent to infrared radiation. The bees could leave and enter the hive at their will through an about 1.5 m long transparent plastic tube with 5 cm inner diameter. The cold stress of the test colonies was varied by placing them in an air stream at experimental temperatures (Texp) of 15, 20, 25, 30 or 34°C. Air temperature near the bees (Ta) was measured by a mesh of 20 thermocouples on each side of the colony, mounted at a height of 5–9 mm above the combs. Actual bee position on the combs was determined relative to a wire mesh with 3.1×3.4 cm rectangles mounted at a height of 10 mm above the combs. By dividing each mesh rectangle into 5 subsections (a-d and center) the bee position could be determined with a resolution of about ±10 mm. Together with laboratory temperature and humidity, and experimental air stream and outside temperature, thermocouple data were stored every 10 s on a laptop computer via a network of one ALMEMO 5590-2 (40-channels) and two ALMEMO 2290-8 data loggers (5-channels each; Ahlborn).

Observation colonies were stocked with about 4000–5000 workers and a queen from a standard colony (Apis mellifera carnica). Sealed brood combs from several standard colonies were transferred to an incubator and kept there at 34°C. Every second morning 150 bees that had emerged during the night before were marked with a color code for age identification by small paint dots (‘Edding 751 paint marker’) on the margin of abdomen and/or thorax to avoid disturbance of thermographic body temperature measurements, and added gently to the upper comb. After 4 weeks of adding bees, where the colonies had time to settle in and establish a brood nest, the first measurements were started. Freshly added bees were allowed to equilibrate and distribute for at least two hours before measurement. Data are from three colonies and 36 measurement sessions (years 2000–2002).

All marked bees of any age found on one side of the colony were measured by scanning the combs with an infrared (IR) thermographic camera (see below). The thermographed colony side was changed regularly between measurement days. As at low experimental temperatures (15 and 20°C) bees hang in several layers on the combs, we first measured all marked bees in the outermost layer. During the second scan the outer bees were gently pushed aside to get a view to the bees and the comb below. From these measurements the age of the youngest bee group was only known to be in the range of 3–27 h. To get a finer age resolution of very young bees, a second series of experiments was conducted (2 replicates at Texp = 15°C experimental temperature, and one at 20, 25, 30 and 34°C each). For each of these experiments, we marked 100 bees 0–1 h after emergence in the morning, and added them gently to the center of the upper comb. Starting after one hour of equilibration, the body surface temperatures just of the bees of this second series were measured at 10:00, 14:00 and 20:00 hours, on both sides of the colony. This way, at each Texp the youngest bees were always of an age of 1–3 h during measurement. Older bees which had been introduced on earlier days (at different Texp) were included. As especially the young bees were often seen to slip into open cells, a separate scan was made where all bees in cells were touched on the abdomen with a blunt needle until they left the cell. Their body temperatures, together with the temperatures of rim and base of the cells where they had been in, were measured immediately after their emergence (within 0–2 s after leaving the cell). In addition to marked bees, all naturally emerged unmarked bees coming out of cells were measured (unmarked bees and ‘nn’ in Figures).

Body and Comb Temperature Measurement

Dorsal honeybee body and comb surface temperature was measured thermographically with a ThermaCam SC2000 NTS equipped with a close-up lens (Fig. 1; FLIR, Inc.; 320×240 pixel sensor, thermal resolution <0.1°C). Thermography allowed measurement of body surface temperature without behavioral impairment. The infrared (IR) camera was calibrated for offset errors against an AGA1010 reference radiator visible in the IR picture via a highly reflective aluminium mirror, or against a custom-made peltier-driven miniature reference source of known temperature and emissivity. Non-uniformity within thermographic pictures was corrected by a VBA macro which took into account differential within-picture drift in dependence on the bees' position in the picture, between the camera's internal shutter calibrations. Attenuation of the IR radiation by the plastic films covering the colonies was compensated for by changing the atmospheric transmission value accordingly during evaluation. Using an infrared emissivity of 0.97 of the insect cuticle and of 0.95 of the comb wax, surface temperature was measured to the nearest 0.8°C by this procedure [44].

[Figure omitted. See PDF.]

Figure 1. Close-up thermogram of honeybees.

Left and right: ectothermic bees (light blue: Tthorax = 25.2°C and 24.9°C). Top and bottom: endothermic bees (white: Tthorax = 35.6°C; orange: Tthorax = 32.4°C). Ta = 25°C, Tcell rim ∼24°C. Measurement performed in a peripheral area of an observation hive during high cold stress (Texp = 20°C).

https://doi.org/10.1371/journal.pone.0008967.g001

Thermographic data were stored digitally with 14 bit resolution on a DOLCH FlexPac-400-XG portable computer at a rate of 2 frames s−1 for the bees on the combs, and at 10 frames s−1 for the second experimental series with very young bees and the bees from within the cells. In addition, thermographic scenes were stored in real time (25 frames s−1) on a SONY Hi8 Video Walkman, together with the spoken commentary on bee marking (age code), position in the wire mesh raster, behavior, and the prevailing type of comb cells at the bee position (open and closed honey or brood cells, pollen cells and empty cells). At the same time we pointed at the bee of interest with a needle which was visible in the video. After the measurements, once per day, the exact ranges of brood cells, honey and pollen stores (open and closed where applicable) and empty cells were drawn on a transparent plastic film laid over the side(s) of the hive where the bee temperatures had been measured.

Data Evaluation

The bees' behavior was judged from the thermographic video sequences and the spoken commentary on the tapes. Evaluation of the surface temperatures of head (Thead), thorax (Tthorax) and abdomen (Tabdomen), and of cell rim and cell base (Tcell rim, Tcell base) was done from the files after the measurements, with AGEMA Research software (FLIR) controlled by a custom programmed Excel VBA macro. Temperature data necessary for exact temperature calculation at the thermographic measurement points were automatically extracted from the logger files and interpolated temporally by this macro. Local air temperature at the actual positions of the bees (Ta) was calculated by triangular interpolation between adjacent thermocouples with a different Excel macro, which extracted the necessary temperature data from the logger files.

Statistics was done with the Statgraphics package (Statistical Graphics Corporation) or with self written Excel sheets (Microsoft Corporation) according to Sachs [45]. Correlations were calculated with Statgraphics or ORIGIN (OriginLab Corporation). For Chi2 statistics the significance level was adjusted according to the Bonferroni correction for multiple comparisons wherever applicable [45].

Results

Location

Fig. 1 shows a typical thermogram of ectothermic and endothermic bees. The dorsal thorax surface temperatures (Tthorax) of marked bees on different locations in the observation hives is shown in Fig. 2. The body surface temperature of the bees varied in a wide range. Both the lowest and the highest Tthorax were measured at an experimental temperature (Texp) of 15°C: minimum Tthorax = 17.5°C (Thead = 17.1°C, Tabdomen = 17.5°C, Ta = 19.8°C, Tcell rim = 18.9°C) and maximum Tthorax = 44.5°C (Thead = 41.4°C, Tabdomen = 39.5°C, Ta = 31.0°C, Tcell rim = 35.9°C). Two exceptionally hot bees on pollen stores at Texp = 34°C with a Tthorax of 46.9 and 47.3°C (asterisks in Fig. 2) were excluded from further analysis because they originated from bees heating their thorax up during intense examination by hivemates [46], [47]. In general, an increasing cold stress (i.e., a decreasing Texp) affected the bees' body surface temperature the least on sealed brood, and somewhat more on open brood. An increasingly higher dependence of Tthorax on Texp was observed in bees on empty and pollen cells, and even more on open and sealed honey cells (Fig. 2).

[Figure omitted. See PDF.]

Figure 2. Worker honeybee thorax temperature on different locations in the observation colonies in dependence on cold stress (Texp).

Box plots: median Tthorax with 1st and 3rd quartile, maximum and minimum. Notches which do not overlap indicate significant differences (P<0.05). Tcell rim and Ta: medians. Numbers  =  n (Tthorax). * are two bees during intense examination by hivemates [45], [46].

https://doi.org/10.1371/journal.pone.0008967.g002

At Texp = 34°C the median Tthorax was slightly higher than Tcell rim and clearly higher than the local ambient air temperature (Ta) on all locations (Fig. 2). At lower Texp, Tthorax was lower than Tcell rim and higher than Ta both on sealed and open brood. This was only partly the case on empty and pollen cells. At Texp <25°C Tthorax was mostly higher on sealed and open brood in comparison with non-brood areas. Bees on open and sealed honey cells were similar in their close decrease of Tthorax with Tcell rim and Ta. They differed from each other in so far as at Texp below 30°C the bees' median Tthorax was always above Ta and partly above Tcell rim on open honey, but always below Ta and Tcell rim on sealed honey (Fig. 2). The variability of body temperature increased considerably with cold stress. This means that at low temperatures an increasing amount of bees actively increased the body temperature noticeably above Ta. Active heat production was observed the least on open and sealed honey cells. A detailed analysis of data showed that ANOVA was only partly permissible due to not normally distributed data, considerable differences in variance and type of distributions to be compared, or missing linear relationships of parameters. A three factor ANOVA with brood status as the main effect and Texp and worker age as covariates showed that, if each factor was measured having removed the effects of all other factors, Texp had by far the highest influence on Tthorax-Tabdomen of bees on combs (in intercomb spaces) (mean square F-ratios1.76816: 1103.83 for Texp, 164.49 for age, and 14.99 for brood; P<0.0001). Median Worker Tthorax-Tabdomen was slightly but significantly higher on brood (0.3°C) than on non-brood areas (0.2°C) (P<<0.0001, Kruskal-Wallis test).

Age and Abundance

Marked bees were present both on the combs and inside cells. In general, considerably more bees were on the combs than inside cells. Concerning the bees on the combs the relative distribution of bees between brood and non-brood areas differed in dependence on age (Fig. 3A, Fig. S1). The young bees (0–2 d) clearly preferred the brood nest. With increasing age this preference disappeared, with an even distribution in the 13–17 d old and 18–22 d old bees. In bees of forager age (>22 d) the relation reversed, bees were more than twice as abundant outside the brood nest as on it. This trend was similarly visible at low and high cold stress (at all Texp; Figs 3A, S1).

[Figure omitted. See PDF.]

Figure 3. Frequency of honeybees of different ages on brood (br) and non-brood (n-br) areas.

Shading of bars  =  thermal stress (Texp). (A) Marked bees on combs; 100% = 12732 bees (measurements) of all age classes. Insert: detailed analysis of 0–2 d old bees. (B) Bees in cells; 100% = 350 bees of all age classes. Percentage of bees in brood nest cells per age class (young to old): 92.5 (184 of 199), 87.5 (70 of 80), 92.3 (36 of 39), 76.2 (16 of 21), 100.0 (5 of 5), 100.0 (6 of 6); no significant differences between age classes (P>0.05, Chi2 = 0.00–6.08). Insert: distribution of unmarked bees on brood and non-brood areas (‘nn’ in Figs 7, 8); n = bees = 100% (per Texp). (A, B) Significant differences between brood and non-brood: * P<0.05, *** P<0.001, Chi2 test.

https://doi.org/10.1371/journal.pone.0008967.g003

Most of the bees inside cells were seen in the brood nest area although many open cells were not occupied.. This was similar in the naturally emerged (unmarked) and the added (marked) bees, even at high temperature (Texp = 34°C; Fig. 3B). With increasing cold stress more bees entered empty cells in the brood nest area. Especially in young bees (<2d) this preference was very pronounced. The tendency to visit cells was much lower in older bees both on the brood nest and outside it (Fig. 3B). Unlike in the bees on combs (in intercomb spaces; Fig. 3A), however, the relation of bees in cells of brood and non-brood areas was similar in all age classes up to the forager age (see legend of Fig. 3B).

Age and Body Temperature

Fig. 4 shows the age dependence of the thorax surface temperature (Tthorax) and its relation to the temperature of the abdomen (Tabdomen), local ambient air (Ta) and comb (Tcell rim) on the brood nest. The medians of Tthorax of the investigated age classes were 35.2–36.0°C when Texp was 34°C (Fig. 4A). Median Tthorax decreased not below 31.4–33.6°C when a very high cold stress acted on the hive at Texp = 15°C. At all Texp conditions the bees younger than 12 h showed the lowest variation of Tthorax (Fig. 4A). In older bees the variation of Tthorax increased considerably, especially during high cold stress (low Texp), reaching the maximum at about 7 days. This was caused by a higher degree of endothermy (see below) and a higher variation of Ta in the locations where the older bees were seen. In the very young bees (<12 h) the median Tthorax decreased less in response to an increase of cold stress than in the older age classes, because they preferably stayed in the well regulated brood nest (Fig. 3).

[Figure omitted. See PDF.]

Figure 4. Body temperature and temperature excess of workers of different ages on brood in dependence on thermal stress (Texp).

Box plots: median with 1st and 3rd quartile, maximum and minimum of (A) Tthorax, (B) Tthorax-Tabdomen, (C) Tthorax-Ta, and (D) Tthorax-Tcell rim. n = measurements (bees) (range in parenthesis).

https://doi.org/10.1371/journal.pone.0008967.g004

ANOVA indicated a dependence of the temperature excess Tthorax-Tabdomen on age both on brood and non-brood areas if the effects of Ta and Tcell rim were removed (mean square F-ratios1,47078: 6.57 for brood and 5.33 for non-brood; P<0.0001). On the brood nest, median Tthorax-Tabdomen was quite similar and at a low level in all age classes at Texp = 34°C (0.10–0.18°C, Fig. 4B). An increase of cold stress influenced Tthorax-Tabdomen not much in the very young bees (≤12 h) but considerably more in the older ones (maximum effect reached in the 2–7 d old bees; P<0.0001 at Texp <30°C, Mann-Whitney U test). Especially in the age classes >12 h at least 75% of the bees had a thorax which was warmer than the abdomen (see 1st quartiles in Fig. 4B). The Tthorax was in most bees higher than the Ta at Texp = 34°C, with no visible tendency to any age dependency (Fig. 4C). During high cold stress (Texp = 15°C) the median Tthorax-Ta increased in general, with mostly more than 75% of an age class showing a Tthorax higher than Ta (see 1st quartiles in Fig. 4C). The variation of Tthorax-Ta increased with cold stress in bees older than 12 h. The Tthorax was quite similar to the comb surface (Tcell rim) at Texp = 34°C (Fig. 4D). Variation of data was again slightly higher in bees older than 12 h. With increasing cold stress the thorax became increasingly cooler than Tcell rim, with a median difference of up to −2.8°C in the >7–12 days old bees at Texp = 15°C. In parallel, in the bees older than 12 h also variation increased. Nevertheless in these age classes a considerable amount of bees (in most cases more than 25%, (see 3rd quartiles in Fig. 4D) remained warmer than the comb surface even during very high cold stress.

Outside the brood nest (Fig. S2) the Tthorax was similar to that on the brood nest in the general trend (especially at high Texp) but different in detail. As bees younger than 54 h were not very abundant outside the brood nest (Fig. 3A) a comparison with the brood nest and between age classes was not possible at all Texp because of low sample sizes (despite a total of 12732 measurements). In general, Texp had a stronger effect on body temperature outside than on the brood nest (Fig. 2). At Texp = 15°C medians were 22.1–33.6°C in comparison to 31.4–33.6°C on the brood nest. Especially in the bees older than 54 h Tthorax decreased stronger with increasing cold stress than on the brood nest (median Tthorax at Texp = 15°C: 22.1–25.7°C outside the brood nest, and 31.4–33.0°C on the brood nest). Tthorax-Tabdomen of bees outside the brood nest, on food and empty cells, was similar to that of bees on the brood nest for all ages and all Texp conditions (medians at Texp = 15°C: 0.32–1.95 versus 0.12–1.32°C, respectively). Tthorax-Ta was on average somewhat smaller outside the brood nest than on it (medians at Texp = 15°C: −0.02–5.9 versus 1.44–4.34°C). Tthorax-Tcell rim was higher outside the brood nest than on it at high Texp (34 and 30°C), especially in bees older than 54 h. At a low Texp of 15°C Tthorax was less frequently below Tcell rim outside the brood nest than on it (medians at Texp = 15°C: −4.07–+0.47 outside versus −2.82–−0.94°C on the brood nest).

Endo- and Ectothermy

In order to get an overview of the number of reliably endothermic bees on the combs we filtered those bees out where Tthorax was not only higher than Tabdomen but also at least 0.2°C higher than Thead: Thead+0.2°C ≤ Tthorax ≥ Tabdomen+dT (dT = +0.2 or +1.0°C; Fig. 5). The percentage of probably endothermic bees on the combs increased with cold stress in all age classes and at all investigated levels of endothermy on the brood nest as well as outside it (Fig. 5A-B). On the brood nest, at Texp = 15°C 18.2% (46 of 253) of the young bees (≤2 d) and 47.6% (654 of 1374) of the bees older than 2 days were endothermic at the dT >1.0°C level. On the other extreme of experimental temperatures (Texp = 34°C), this percentage was reduced to 0.0% (0 of 0) in the young and to 0.5% (4 of 758 bees) in the older bees (statistics in Fig. 5). Young bees (0–2 d) were less often endothermic than the older ones at Texp = 15 and 25°C, and showed a lower tendency at Texp = 20°C. This tendency seemed to be more pronounced outside the brood nest (inserts in Fig. 5A-B) than on it.

[Figure omitted. See PDF.]

Figure 5. Percentage of endothermic bees in dependence on thermal stress (Texp).

Shading of bars  =  age. (A-B) Percentage of actively heating workers on brood (main graphs) and non-brood areas (inserts); (C-D) percentage of actively heating workers in open brood nest cells. Parts (A, C) and (B, D) display different levels of endothermy (Tthorax-Tabdomen). Number of investigated workers of each age class ( = 100%) written above each column in (C) and (D). Number of measurements for non-brood areas (inserts) with increasing age: Texp = 15°C: 14, 47, 140, 971; Texp = 20°C: 12, 127, 157, 651; Texp = 25°C: 85, 227, 125, 698; Texp = 30°C: 72, 178, 172, 562; Texp = 34°C: 68, 189, 163, 776. Statistics: *, differences significant (P<0.05, chi23.841 >3.94) between 0–2 d old and >2 d old bees at a certain Texp (at Texp = 20°C, P<0.1, chi22.706 = 3.484). #, differences not significant (P>0.05) between Texp = 15°C and Texp >15°C in each age class; otherwise P<0.05 (chi26.239 >6.94 in A, B and chi25.731 >7.39 in C, D).

https://doi.org/10.1371/journal.pone.0008967.g005

Inside brood nest cells (Fig. 5C-D), at Texp = 15°C, 66.7% (30 of 45) of the bees older than 2 days were endothermic at the dT >1.0°C level. At Texp = 34°C, this percentage was 0.0% (0 of 8 bees). Young bees in cells (0–2 d) were endothermic at the dT >1°C level at a rate of 47.1% (8 of 17) and 0.0% (0 of 14) at Texp = 15 and 34°C, respectively (statistics in Fig. 5).

A comparison of endothermic activity between bees inside cells and on combs on the basis of the above thermal classification has to consider that, if Tthorax-Tabdomen is used as a correlate of endothermic activity, local thermal gradients might be higher for bees in cells and lead to an overestimation of their endothermic activity. To compare these local ambient gradients, reliably ectothermic bees which are unable of heating the thorax (0–24 h old) can be used as gradient thermometers. A correlate of the ambient gradient along the body is the Thead-Tabdomen. In the brood nest, even at very high thermal stress (Texp = 15°C), it was not significantly different between bees in cells and on combs despite a visible trend (1.12±1.05°C in cells vs. 0.49±0.84°C on combs; n = 8 vs. 79; P = 0.142, t test). At moderate thermal stress (Texp = 25°C) gradients of Thead-Tabdomen were nearly identical (0.11±0.35°C in cells vs. 0.09±0.35°C on combs; n = 43 vs. 274; P = 0.732).

The rate of bees classified as endothermic at the dT >1°C level was not different between bees in brood nest cells (Fig. 5B) and on combs (Fig. 5D) at Texp = 34°C (chi23.841 = 0–0.02). Heating at this Texp was extremely rare both inside and outside cells in all age classes (0–2 d, 3–7 d, ≥7 d tested). At lower Texp the 3–7 and >7 days old bees were always more often classified as endothermic in the cells than on the combs (Fig. 5). Differences were significant at a Texp of 30 and 25°C (P<0.05; chi23.841 = 4.76–21.1). At a Texp of 15°C the difference was weakly beyond significance in the ≥7 days old bees (P<0.1; chi23.841 = 3.47), and not significant in the 3–7 days old bees (chi23.841 = 2.52). The young bees (0–2 d) were more often endothermic in the cells than on the combs only at Texp = 15°C (chi23.841 = 9.16) but not at higher Texp.

Local Thermal Gradients

We found highly significant correlations between Ta and Texp, and between Tcell base and Tcell rim (Table 1, Fig. S3). A comparison between brood and non-brood areas showed that these gradients were influenced by the bees' endothermic activity (Fig. S3). The effect of endothermy on thermal gradients was especially pronounced for the relation of Tcell rim and Ta, which could not be described by a simple linear regression (Fig. 6). The temperature gradient in cells (median Tcell base-Tcell rim) increased with increasing cold stress. In empty cells not visited by bees it was −0.1, 0.4, 0.5, and 1.3°C at Texp = 34, 30, 25 and 15°C, respectively. In cells visited by unmarked bees it was quite similar, amounting to −0.2, 0.7, 0.6 and 1.4°C (Fig. 7C), respectively. There was no statistical difference between not visited and visited cells (P≥0.079, Mann Whitney U test) except at Texp = 30°C (P<0.001). The median cell temperature gradients were −0.25, 0.5, 0.4 and 1.7°C in cells visited by 0–2 days old bees, −0.1 0.5, 0.6 and 2.4°C in cells visited by >3–7 days old bees, and −0.4, 0.8, 0.8 and 1.3°C in cells visited by ≥7 days old bees, respectively (Fig. 7C).

[Figure omitted. See PDF.]

Figure 6. Intensity of endothermy (Tthorax-Tabdomen) of bees on combs in relation to individual bees' local ambient temperatures (Tcell rim and Ta).

See Fig. S1 for correlations between Tcell base and Tcell rim, and between Ta and Texp.

https://doi.org/10.1371/journal.pone.0008967.g006

[Figure omitted. See PDF.]

Figure 7. Thermal relations in open brood nest cells visited by bees.

Box plots: median with 1st and 3rd quartile, maximum and minimum. Notches which do not overlap indicate significant differences (P<0.05). nn  =  unmarked workers; red crosses  =  medians of all age groups of marked bees (not significantly different from nn bees, except Tthorax-Tabdomen at Texp = 25 and 30°C, P<0.05, U test).Tcell 30% is the cell wall temperature at the estimated mean contact position of the thorax with the cell wall at 30% cell depth as measured from the base, calculated by linear interpolation between Tcell base and Tcell rim. Numbers  =  measurements (bees).

https://doi.org/10.1371/journal.pone.0008967.g007

[Figure omitted. See PDF.]

Table 1. Regressions between hive temperatures and experimental temperature.

https://doi.org/10.1371/journal.pone.0008967.t001

In unmarked bees inside brood nest cells, which had emerged naturally in the test colonies (‘nn’ in Fig. 7A), the median gradient of Tthorax-Tabdomen increased from just 0.1°C at Texp = 34°C to 2.0°C at Texp = 15°C. In parallel, the variation of Tthorax-Tabdomen increased considerably. The range of central 50% (1st - 3rd quartile), for example, increased from 0.2 at Texp = 34°C to 1.5°C at Texp = 15°C, the total range from 1.0 at Texp = 34°C to 7.2°C at Texp = 15°C. Fig. 7A also shows that the increase with cold stress was smaller in young bees (0–2 d: 0.1 to 1.6°C) in comparison to the older ones (>2 d: 0.1 to 1.9–2.3°C; comparison of 0–2d old with >2d old bees: P<0.05 at Texp = 15°C, P<<0.0001 at Texp = 25 and 30°C, P = 0.67 at Texp = 34°C; Mann Whitney U test).

Heat Transfer between Bees in Cells and Combs

A considerable amount of the difference of Tthorax-Tabdomen was caused by the temperature gradient in the cells (compare Fig. 7A with Fig. 7C). The center of honeybee endothermy is the thorax. We estimated the position of contact between the thorax (the heat source during endothermy) and the cell wall to be at about 30% of the cell depth as measured from its base on average. In order to judge the direction of the heat flow between bees and combs, we calculated the difference between thorax and its cell contact position (Tcell 30%) by interpolating the temperature between cell base and cell rim linearly (Fig. 7B). If the test colonies were not exposed to cold stress (Texp = 34°C) the majority of the bees in brood nest cells was warmer than the cell wall (Fig. 8A), and there were no differences between the age groups (P>0.05; Chi23 comparisons = 0.35–1.75<5.731). At this Texp, however, the presence inside cells was generally low (Fig. 3B). The unmarked bees in the cells had a median Tthorax of 0.65°C above Tcell 30% at Texp = 34°C but of 0.34°C below Tcell 30% at Texp = 15°C (‘nn’ in Fig. 7B). Variation of central 50% increased from 0.5 to 3.4°C, and the total range from 1.8 to 6.0°C. The age analysis showed that predominantly the young bees (0–2 d) became cooler (−0.44°C) than Tcell 30% at the lowest Texp of 15°C. In older bees (>2–7 d) less individuals had a Tthorax below Tcell 30%. The majority of bees older than 7 days had a Tthorax higher than Tcell 30% at Texp = 15°C (Fig. 7B).

[Figure omitted. See PDF.]

Figure 8. Heat transfer between bees in cells and brood combs.

Percentage of workers of different ages in brood nest cells with (A) higher Tthorax than Tcell 30%, (B) approximately equal Tthorax and Tcell 30%, and (C) lower Tthorax than Tcell 30%, in dependence on thermal stress (Texp  =  shading of bars). Sample numbers in (C) represent 100%. nn  =  naturally emerged, unmarked bees. nn bees did not differ from the averaged age classes of marked bees (P>0.1; Chi22.706 = 0.00–2.143) with exception of bees at Texp = 30°C in (B) (P<0.05; Chi23.841 = 4.163).

https://doi.org/10.1371/journal.pone.0008967.g008

The rate of bees inside cells being warmer than the cell wall (Fig. 8A) was not different between the >2–7 days old and the >7 days old bees (P>0.05; Chi23 comparisons = 0.64–2.19<5.731; all Texp). At Texp = 25 and 30°C, 0–2 days old bees transferred less often heat to the comb than the combined older age classes (P<0.01; Chi2 = 10.14 and 8.63>6.635). At Texp = 15°C the visible trend could not be proved statistically (P<0.1; Chi2 = 3.46<3.841). The unmarked bees in cells (‘nn’), which had emerged naturally in the observation hives, serve as an independent corroboration for the reliability and unbiased nature of these results. They showed the same dependence on Texp as the marked bees both concerning body temperature relations (Fig. 7) and probability of heat transfer to the combs (Fig. 8).

Discussion

Mechanisms of Thermal Homeostasis

The honeybees' endothermic activity is one of the main factors to establish thermal homeostasis in a colony. In insects capable of endothermy, the flight muscles are the source of active heat production [e.g. 1], [44], [48], [49]. Since endothermy increases the thorax temperature in relation to the abdomen, the temperature difference between thorax and abdomen (Tthorax-Tabdomen) is a good correlate of endothermic activity. In our experiments Tthorax-Tabdomen of individuals was strongly affected by Texp. As the thorax temperature excess Tthorax-Tabdomen is always a mixture of local temperature gradient and degree of endothermic heat production, more clarity on endothermic activity can be achieved by filtering those bees out where the thorax was the warmest body part. Such analysis shows that the frequency of endothermic bees is higher on the brood nest than in non-brood areas, especially at extreme cold stress (Fig. 5, 6). However, there is always a considerable number of endothermic bees outside the brood nest. These bees contribute indirectly to thermal homeostasis of the brood because a minimized temperature gradient means a reduced heat flow out of the brood nest. Foragers, which are always endothermic during their returns to the colony [50]–[53], also contribute to this reduction of temperature gradients though they preferably stay in peripheral comb areas near the entrance during daytime [54], [55]. During night-time, by contrast, foragers were observed to cluster in peripheral hive areas. Though the majority of them is at rest [56] and not endothermic in a standard colony (our own unpublished observations), their high number nevertheless provides a considerable amount of resting heat production [31]. In addition, their increased number at nighttime contributes to colony isolation.

It has to be kept in mind that a Texp of 15°C in our 2-comb observation hives represents a much higher thermal stress for the bees than in standard colonies. The heat loss of standard colonies is essentially smaller because of the parallel arrangement of several combs, better isolation, a higher number of bees, and (usually) less external convection than in our experiments. We estimated an environmental temperature (Texp) of about 0–10°C to initiate an extent of endothermy in standard colonies comparable to our observation hives at Texp = 15°C (our own unpublished experiments; [49]).

Active compensation of heat loss.

Our measurements showed that from low to high cold stress (Texp = 30 to 15°C) a considerable number (mostly noticeably more than 50%) of the bees on the combs were cooler than the comb surface (Tcell rim; Fig. 4D). In addition, the surface of open cells was cooler than their base (Fig. 7; see also [10]). This means that in many locations the heat flow is from the combs to the bee ways (intercomb spaces) and into the ectothermic bees there. To achieve thermal constancy of the brood this heat loss has to be compensated and minimized. Compensation is done by the endothermic bees on the comb surface and inside open cells (Figs 4–8; [9], [10]), and for the smaller part also by drones [57]. It might seem strange that average thorax temperatures of bees on combs (in intercomb spaces) are below the cell rim temperatures at low environmental temperatures, as the bees are the main source of the heat. However, thermal gradients may still persist as long as the heat production of the larvae and pupae themselves suffices to compensate for the heat loss. The duty of the bees is to minimize but not completely remove thermal gradients. Gradients have just to remain stable. If larvae are young and thus small, compensation by the bees has to be increased. Even more, heat production of the brood necessitates maintenance of (small) thermal gradients. If environmental temperature (Texp) is 34°C, where gradients in cells already reverse (Fig. 7), bees already spread water on combs to get rid of excess larval heat [24], [49].

The intensity of endothermy, however, varies in a wide range, and many bees are only weakly endothermic (Fig. 4). Our observations on endothermic reactions represent a mean of general trends. Actually, individual bees do not heat continuously but show heating bouts lasting several minutes, alternating with ectothermic phases [10]. The differential decisions of individual bees to heat or not is suggested to be determined by temperature threshold differences between individuals, according to the threshold model of task allocation in polyandrous social insect colonies [e.g. 58]. A higher genetic diversity (with a higher probability of differing thresholds) was reported to improve thermal homeostasis of honeybee colonies [27]. We suggest that individual endothermy is kept up until an individual optimum of local temperature (probably Tcomb, see below) is achieved. The sum of temporally and spatially distributed individual phases of endothermy enables the superorganism honeybee colony to achieve thermal homeostasis.

Active (endothermic) compensation of heat loss is done by relatively few bees. On the basis of an estimation of the relation between thorax temperature excess above ambient air (Tthorax-Ta) and oxygen consumption in endothermic bees [59], [60], and measurements of the resting metabolism [30], [31] a rough estimation is possible how much heat an endothermic bee dissipates in comparison to ectothermic bees. An endothermic bee with a temperature excess (Tthorax-Ta) of 1°C above the resting level (at Ta ∼34°C) produces roughly 2–4 times the heat of resting older bees (>2 d), and approximately 4–6 times the heat of freshly hatched bees (<12 h). A temperature excess of 5°C compensates for the resting heat production of approximately 11–16 bees capable of endothermy (>2 d), and of 16–23 young bees (<12 h) which are not able of endothermic heat production, respectively. So heat production of a few endothermic bees eases much burden of many nestmates. Their heat surely has a considerable effect on local ambient temperature [compare 9]. This affects local comb temperature most if bees heat cells from within [10].

The thermal activity of endothermic bees does not only apply heat directly to the surrounding air and the combs but has also an important reinforcement effect on passive heat production of the ectothermic bees: it shifts them to a higher basal metabolic rate along their curve of resting metabolism [29]–[31], [61]. The same holds for the larvae [21]. In Fig. 6 this is represented by the fact that the distribution of endothermic bees bases upon a large ‘carpet’ of bees with no or just weak active heat production.

Passive minimisation of heat loss.

Heat flow minimisation is achieved by the presence of bees in the bee ways and inside of part of the open cells. They not only reduce the convective but also the diffusive heat loss by reduction of the thermal gradient. The resting metabolism is, because of the high number of ectothermic bees, an important parameter for heat flow minimisation and compensation.

Heat transfer in cells.

Usually, many of the open brood cells are not visited by bees for longer periods of time. A high presence inside cells was only observed after increased emergence the day before. Most of these bees are young (Fig. 3B) and therefore not capable of much endothermy [62]. Of the older bees it is mainly the nurses (3–12 d old) which visit open brood cells, for controlling and feeding the larvae, and for active heat production (Figs 5, 7, 8).

The temperature gradient in brood nest cells (Tcell base-Tcell rim) increased significantly with increasing cold stress (decreasing Texp; Fig. 7C), but we did not find a statistical difference between empty and visited cells. This, however, would not justify the conclusion that visitation has no effect on cell temperature. Rather, visitation may have compensated for temperature differences. Direct thermographic observations demonstrate this [10]. At Texp = 34°C, the average (median) heat flow was always directed inwards, from the bee ways to the comb (Tcell base–Tcell rim = −0.25°C; Fig. 7C). At lower Texp the average heat flow was always outwards. We suggest the equilibrium environmental temperature (Texp) of our observation hives concerning heat flow in visited open brood cells to have been at about 32.5–33°C (see Fig. 7C). In standard colonies we suggest it to be below this level because of a better volume to surface ratio, better insulation and less internal convection because of the parallel arrangement of the combs.

If it is to judge the heat flow inside visited cells in more detail, the difference between the thorax (the source of heat production) and the cell wall is of particular interest. Inside cells rather weak endothermy may reverse the direction of heat flow. At high cold stress the median Tthorax especially of the 0–2 days old bees was below the cell temperature at the point of contact with the cell wall (Tcell 30%, Fig. 7C). In these cases the heat flow was from the cell wall to the thorax (Fig. 8). As the older bees (>2 d) were more frequently endothermic, the direction of heat flow was in more bees directed from the thorax to the cell wall. It has to be kept in mind, however, that these considerations refer to a value of Tcell 30% interpolated between Tcell base and Tcell rim. This works best in ectothermic bees. In the case of endothermic bees in cells, the gradient between cell base and cell rim is probably not linear. Somewhat higher temperatures than calculated by interpolation have to be expected at the contact site of a heated thorax with the cell wall [10]. Nevertheless, as the heat applied to the contact site distributes within the cells in time, our estimation describes the average heat flow inside visited cells. Bees inside cells, however, contribute to comb temperature stabilisation not only in the endothermic but also in the ectothermic state. Ectothermic bees are not just passive insulators which reduce convective heat loss due to their presence. In the narrowness of a cell, their rather high resting metabolism (and thus heat production) [31] at brood comb temperature additionally reduces the thermal gradient and thus the heat flow out of the cells. Heat transport from the cells to the bee ways via respiration is suggested to be of minor importance (though not negligible) in ectothermic bees, because at Ta = 30 and 34°C (as is usual on a brood nest), resting bees breath only every 37 and 28 seconds on average, respectively [31], and because the temperature gradient inside cells is rather small anyway at low to moderate temperature stress (Fig. 7A).

Age Polyethism

Behavioral analysis of temporal polyethism within honeybee colonies usually addresses task allocation in terms of the visually observable activity like cell cleaning (days 1–3), nursing (days 4–12), food processing and nest maintenance (days 13–20), and foraging (days >20) [38]–[43]. Harrison [32] was the first to investigate the temporal polyethism of workers concerning active heat production of individual bees for the maintenance of thermal homeostasis. He reported that, under very high cold stress (one single comb exposed to 15°C), he was not able to find age-dependent differences in the degree of endothermy in bees older than two days. Only 1-day-old bees showed a smaller increase of Tthorax above Ta. In our experiments this was best visible with the temperature excess of Tthorax-Tabdomen (Fig. 4B). Especially the youngest bees (<24 h) showed no or just a very weak endothermic activity both in terms of intensity and percentage of endothermic workers on the combs (Figs 4, 5). We confirm the finding of Harrison [32] that the nurses do not differ from middle aged bees and foragers. This finding, however, might have resulted from the extreme cold stress in those experiments, forcing more age divisions to take part in active heat production as usual. Our graduated manipulations of cold stress did not result in different age distributions of endothermic activity. The same holds for drones, though their presence on the brood nest is considerably lower in comparison to the workers [57]. In addition, our experiments show for the first time that this not only applies for workers on the combs but also for bees in cells. The rate of bees inside cells being warmer than the estimated cell wall temperature (Tcell 30%), i.e. the number of bees transferring heat to the combs, also was not different between the age groups older than two days. Only the youngest bees (0–2 d) showed a lower probability of heat transfer to the combs (Fig. 8A).

Concerning active heat production, therefore, the four classical divisions of task allocation reduce to only two: 0 to ∼2days and older. The shift between these two task divisions is gradual as is common in honeybee polyethism, starting at an age of about 24 hours and being finished at about 5 days (Fig. 4). This resembles the age dependent development of endothermic ability [62]. Task allocation concerning endothermy, therefore, directly depends on physiological ability. This has the consequence that the passive effect of cold stress is highest in very young bees, whereas the bees' physiological (endothermic) reaction to cold stress is highest in older bees. Young bees compensate for their inability of active heat production by behavioral mechanisms, i.e.by their increased preference of the warm brood nest ambience (Fig. 3A, [63]) and their tendency to enter empty cells (Fig. 3B). As endothermic capacity increases in the nurse age (3–13 d, [62]) the brood nest preference decreases (Fig. 3A, Fig. S1). Bees engaged in hive maintenance and food storing (>13–22 d) are already as abundant outside the brood nest as on it. All ages of honeybees can be seen as ‘active isolators’ with respect to environmental changes. The young ones preferably in terms of migration activity into or out of brood nest cells (Fig. 3B, insert), and the older bees both by migration activity to or from the brood nest (Fig. S1; [49]) and additionally in terms of active reduction of thermal gradients by endothermy.

In conclusion it can be said that thermal homeostasis in honeybee nests is the result of the cooperation of all colony members and all temporal castes of adults. This includes heat production by the larvae and pupae [19], [21], passive worker and drone heat production via the resting metabolism [30], [31], [64], active (endothermic) heat production by hive bees and drones on the combs [9], [57, this paper] and inside cells [10, this paper], active heat production by foragers [48], [50]–[53], migratory movements of workers on the combs (changes of bee density; this paper and [49]), cooling of combs by foragers [22], [24] and via social control of colony ventilation [23], and heat shielding [65], [66].

Individual Development of Endothermy

Our measurements showed that especially the very young bees of 0–2 days age seek the warm area of the brood nest (Fig. 3A). This preference is less pronounced in 3–7 days old bees, and disappears completely in older workers. If empty cells are available it is predominantly the very youngest bees which enter them (Fig. 3B). These (young) bees are usually thought to enter the cells to clean them [e.g. 37], [38]. However, in search for the causes that govern the young bees' preference for the brood nest and the empty cells, one has to look beyond arguments of social behavior, and additionally consider physiological necessities. From a physiological point of view bees are not adults after emergence. They are not yet capable of proper activation of their flight muscles both for flight and endothermic heat production. Capability of endothermy develops only within the first days after emergence [44], [62]. Morphological and enzymatic make-ups are completely developed only at an age of about 8–9 days [67]–[71]. As freshly emerged bees are poikilothermic [44], their respiratory turnover, which is the basis for development, is strongly temperature dependent [30]. We suggest the young bees to seek the high brood nest and cell temperature for self-warming to guarantee proper development.

Stimuli of Endothermy

We suggest the individual decision of bees to heat or not to depend on their local environment, i.e. Tcomb and/or Ta, but most probably not on the outside temperature (Texp). We do not know, however, whether the ‘local environment’ is the actual environment, or whether bees make a temporal integration on their way across the combs. Recent work has shown that especially the older age classes (‘middle aged bees’, 13–20 d) show considerable migration activity, presumably also for information sampling, beside search for work [43].

Only suggestions exist about what bees measure, comb or air temperature. Besides measuring comb temperature they might just regulate Ta at a level which guarantees a minimum heat flow away from the combs. From experiments with combs cooled with a heat exchange plate, Kronenberg and Heller [20] suggested the comb temperature to be the stimulus. Bujok et al. [9] reported intense antennal contact with the comb surface of endothermic workers pressing their body on the comb surface to facilitate heat transfer. An analysis searching for the stimulus temperatures for the release of active heating which can be perceived by bees, has to take into consideration that there exists a chain of interdependence between Tcell base, Tcell rim, Ta, and Texp. While we found direct correlations between Tcell base and Tcell rim and between Ta and Texp (Table 1, Fig. S3), the interrelation of Tcell rim and Ta at the bees' positions was more complex, especially on the brood nest (Fig. 6). This ‘complexity’ is caused by the bees' endothermic counteraction to stabilize brood temperature as a reaction to changes of environmental temperature (our manipulations of cold stress). On the brood nest, there was a clear concentration of strongly endothermic bees (Tthorax–Tabdomen >3°C) around a Tcell rim of about 35°C, whereas along the Ta axis endothermy was more evenly distributed. We conclude that the attempts to stabilize Tcell rim directly were higher than to stabilize Ta, which is supported by observations of direct heat transfer to the comb surface [9] and to the cell interior [10]. Our data show that bees inside brood nest cells are more often endothermic than bees on the surface of brood combs (Fig 5). This supports the hypothesis that bees undertake special efforts to stabilize comb temperature by applying heat as near to the brood as possible. Therefore, we suggest the bees to react to comb temperature first of all and to air temperature in the second place.

Comparison of endothermic activity between brood and non-brood areas (Fig. 6) and literature reports on the relation of brood and non-brood temperatures [13], [18], [20] suggest that accurate regulation of brood temperature requires a trigger of higher order to coordinate endothermic activity and behavioral reactions between colony members. In addition to temperature, therefore, brood pheromones are suggested to play an important role in brood temperature homeostasis [13].

Conclusion

Thermal homeostasis of honeybee colonies is achieved by a combination of active regulatory processes and passive effects. The central parameter of active processes is “endothermy on demand” of adult bees (older than about 2 days), which has an important reinforcement effect on passive heat production of the many ectothermic adult (>∼2 days) and young (<∼2 days) bees, as well as of larvae and pupae. The second foundation of thermal homeostasis is behavior. This includes migration activity within the colony (changes of bee density), evaporative cooling with water, and regulation of internal heat transport via convection (fanning).

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Distribution of workers of different ages (A-J) between brood (open and sealed) and non-brood areas in dependence on age and thermal stress (Texp). Mixed brood additionally classified in young bees (2–54 h old) (A - E). Numbers  =  n = 100% of columns.

https://doi.org/10.1371/journal.pone.0008967.s001

(7.51 MB EPS)

Figure S2.

Body temperature of workers of different ages on non-brood areas in dependence on thermal stress (Texp). Box plots: median with 1st and 3rd quartile, maximum and minimum of (A) Tthorax, (B) Tthorax-Tabdomen, (C) Tthorax-Ta, and (D) Tthorax-Tcell rim. n  =  measurements (bees), low sample sizes of 1 or 2 are indicated.

https://doi.org/10.1371/journal.pone.0008967.s002

(1.82 MB EPS)

Figure S3.

Correlations of temperatures in observation hives. (A, B) Correlation of cell temperatures. (B, C) Correlation of individual bees' local ambient temperature (Ta) and environmental temperature (Texp). See Fig. 6 for relation of Tcell rim and Ta.

https://doi.org/10.1371/journal.pone.0008967.s003

(1.06 MB TIF)

Acknowledgments

We thank Gabriel Stabentheiner and Stefan K. Hetz for software and electronics support. Thanks are also due to Manuela Ablasser, Bettina Germann, Lisbeth Postl and Jutta Vollmann for help with data evaluation, and two anonymous reviewers for helpful comments.

Author Contributions

Conceived and designed the experiments: AS HK. Performed the experiments: AS HK. Analyzed the data: AS HK RB. Contributed reagents/materials/analysis tools: AS HK. Wrote the paper: AS HK RB.

Citation: Stabentheiner A, Kovac H, Brodschneider R (2010) Honeybee Colony Thermoregulation – Regulatory Mechanisms and Contribution of Individuals in Dependence on Age, Location and Thermal Stress. PLoS ONE5(1): e8967. https://doi.org/10.1371/journal.pone.0008967

References

1. Heinrich B (1993) The Hot-Blooded Insects. Berlin, Heidelberg, London, Paris: Springer. B. Heinrich1993The Hot-Blooded Insects.Berlin, Heidelberg, London, ParisSpringer601 p

2. Gates V (1914) The temperature of the bee-colony. Bull US Dept Agric 96: 1–19.V. Gates1914The temperature of the bee-colony.Bull US Dept Agric96119

3. Hess WR (1926) Die Temperaturregulation im Bienenvolk. Z Vgl Physiol 4: 465–487.WR Hess1926Die Temperaturregulation im Bienenvolk.Z Vgl Physiol4465487

4. Himmer A (1932) Die Temperaturverhältnisse bei den sozialen Hymenopteren. Biol Rev 7: 224–253.A. Himmer1932Die Temperaturverhältnisse bei den sozialen Hymenopteren.Biol Rev7224253

5. Büdel A (1955) Schwankungen der Lufttemperatur in der Wabengasse eines brütenden Bienenvolkes. Z Bienenforsch 3: 88–92.A. Büdel1955Schwankungen der Lufttemperatur in der Wabengasse eines brütenden Bienenvolkes.Z Bienenforsch38892

6. Büdel A (1960) Bienenphysik. In: Büdel A and Herold E, eds. Biene und Bienenzucht. Munich: Ehrenwirth. A. Büdel1960Bienenphysik. In: Büdel A and Herold E, eds. Biene und Bienenzucht.Munich: Ehrenwirth435 p

7. Simpson J (1961) Nest climate regulation in honey bee colonies. Science 133: 1327–1333.J. Simpson1961Nest climate regulation in honey bee colonies.Science13313271333

8. Fahrenholz L, Lamprecht I, Schricker B (1992) Calorimetric investigations of different castes of honey bees, Apis mellifera carnica. J Comp Physiol B 162: 119–130.L. FahrenholzI. LamprechtB. Schricker1992Calorimetric investigations of different castes of honey bees, Apis mellifera carnica.J Comp Physiol B162119130

9. Bujok B, Kleinhenz M, Fuchs S, Tautz J (2002) Hot spots in the bee hive. Naturwissenschaften 89: 299–301.B. BujokM. KleinhenzS. FuchsJ. Tautz2002Hot spots in the bee hive.Naturwissenschaften89299301

10. Kleinhenz M, Bujok B, Fuchs S, Tautz J (2003) Hot bees in empty broodnest cells: heating from within. J Exp Biol 206: 4217–4231.M. KleinhenzB. BujokS. FuchsJ. Tautz2003Hot bees in empty broodnest cells: heating from within.J Exp Biol20642174231

11. Southwick EE (1985) Allometric relations, metabolism and heat conductance in clusters of honey bees at cool temperatures. J Comp Physiol B 156: 143–149.EE Southwick1985Allometric relations, metabolism and heat conductance in clusters of honey bees at cool temperatures.J Comp Physiol B156143149

12. Stabentheiner A, Pressl H, Papst T, Hrassnigg N, Crailsheim K (2003) Endothermic heat production in honeybee winter clusters. J Exp Biol 206: 353–358.A. StabentheinerH. PresslT. PapstN. HrassniggK. Crailsheim2003Endothermic heat production in honeybee winter clusters.J Exp Biol206353358

13. Koeniger N (1978) Das Wärmen der Brut bei der Honigbiene (Apis mellifera L.). Apidologie 9: 305–320.N. Koeniger1978Das Wärmen der Brut bei der Honigbiene (Apis mellifera L.).Apidologie9305320

14. Tautz J, Maier S, Groh C, Roessler W, Brockmann A (2003) Behavioral performance in adult honey bees is influenced by the temperature experienced during their pupal development. Proc Natl Acad Sci USA 100: 7343–7347.J. TautzS. MaierC. GrohW. RoesslerA. Brockmann2003Behavioral performance in adult honey bees is influenced by the temperature experienced during their pupal development.Proc Natl Acad Sci USA10073437347

15. Groh C, Tautz J, Roessler W (2004) Synaptic organization in the adult honey-bee brain is influenced by brood-temperature control during pupal development. Proc Natl Acad Sci USA 101: 4268–4273.C. GrohJ. TautzW. Roessler2004Synaptic organization in the adult honey-bee brain is influenced by brood-temperature control during pupal development.Proc Natl Acad Sci USA10142684273

16. Jones CJ, Helliwell P, Beekman M, Maleszka R, Oldroyd BP (2005) The effects of rearing temperature on developmental stability and learning and memory in the honeybee, Apis mellifera. J Comp Physiol A 191: 1121–1129.CJ JonesP. HelliwellM. BeekmanR. MaleszkaBP Oldroyd2005The effects of rearing temperature on developmental stability and learning and memory in the honeybee, Apis mellifera.J Comp Physiol A19111211129

17. Becher MA· Scharpenberg H, · Moritz RFA (2009) Pupal developmental temperature and behavioral specialization of honeybee workers (Apis mellifera L.). J Comp Physiol A 195: 673–679.H. Becher MA· ScharpenbergRFA · Moritz2009Pupal developmental temperature and behavioral specialization of honeybee workers (Apis mellifera L.).J Comp Physiol A195673679

18. Ritter W (1982) Experimenteller Beitrag zur Thermoregulation des Bienenvolkes (Apis mellifera L.). Apidologie 13: 169–195.W. Ritter1982Experimenteller Beitrag zur Thermoregulation des Bienenvolkes (Apis mellifera L.).Apidologie13169195

19. Melampy RM, Willis ER (1939) Respiratory metabolism during larval and pupal development of the female honeybee (Apis mellifica L.). Physiol Zool 12: 302–311.RM MelampyER Willis1939Respiratory metabolism during larval and pupal development of the female honeybee (Apis mellifica L.).Physiol Zool12302311

20. Kronenberg F, Heller C (1982) Colonial Thermoregulation in Honey Bees (Apis mellifera). J Comp Physiol 148: 65–76.F. KronenbergC. Heller1982Colonial Thermoregulation in Honey Bees (Apis mellifera).J Comp Physiol1486576

21. Petz M, Stabentheiner A, Crailsheim K (2004) Respiration of individual honeybee larvae in relation to age and ambient temperature. J Comp Physiol B 174: 511–518.M. PetzA. StabentheinerK. Crailsheim2004Respiration of individual honeybee larvae in relation to age and ambient temperature.J Comp Physiol B174511518

22. Lindauer M (1954) Temperaturregulierung und Wasserhaushalt im Bienenstaat. Z Vgl Physiol 36: 391–432.M. Lindauer1954Temperaturregulierung und Wasserhaushalt im Bienenstaat.Z Vgl Physiol36391432

23. Southwick EE, Moritz RFA (1987) Social control of air ventilation in colonies of honey bees, Apis mellifera. J Insect Physiol 33: 623–626.EE SouthwickRFA Moritz1987Social control of air ventilation in colonies of honey bees, Apis mellifera.J Insect Physiol33623626

24. Mandl M, Stabentheiner A, Kovac H (2004) Cooling by water droplets in a breeding honeybee colony. M. MandlA. StabentheinerH. Kovac2004Cooling by water droplets in a breeding honeybee colony.In: Proceedings of the First European Conference of Apidology EurBee (Bernardinelle, I & Milani, N, eds.), Arti Grafiche Friulane SpA, Udine. 59 p.

25. Southwick EE (1991) The colony as a thermoregulating superorganism. In: Goodman LJ, Fisher RC, editors. The Behaviour and Physiology of Bees. Wallingford UK: CAB International. pp. 28–47.EE Southwick1991The colony as a thermoregulating superorganism.LJ GoodmanRC FisherThe Behaviour and Physiology of BeesWallingford UKCAB International2847

26. Moritz RFA, Southwick EE (1992) Bees as superorganisms. Berlin, Heidelberg New York: Springer. RFA MoritzEE Southwick1992Bees as superorganisms.Berlin, Heidelberg New YorkSpringer601 p

27. Jones JC, Myerscough MR, Graham S, Oldroyd BP (2004) Honey bee nest thermoregulation: diversity promotes stability. Science 305: 402–404.JC JonesMR MyerscoughS. GrahamBP Oldroyd2004Honey bee nest thermoregulation: diversity promotes stability.Science305402404

28. Mackasmiel LAM, Fell RD (2000) Respiration rates in eggs of the honey bee, Apis mellifera. J Apic Res 39: 125–135.LAM MackasmielRD Fell2000Respiration rates in eggs of the honey bee, Apis mellifera.J Apic Res39125135

29. Schmolz E, Hoffmeister D, Lamprecht I (2002) Calorimetric investigations on metabolic rates and thermoregulation of sleeping honeybees (Apis mellifera carnica). Thermochim Acta 382: 221–227.E. SchmolzD. HoffmeisterI. Lamprecht2002Calorimetric investigations on metabolic rates and thermoregulation of sleeping honeybees (Apis mellifera carnica).Thermochim Acta382221227

30. Stabentheiner A, Vollmann J, Kovac H, Crailsheim K (2003) Oxygen consumption and body temperature of active and resting honeybees. J Insect Physiol 49: 881–889.A. StabentheinerJ. VollmannH. KovacK. Crailsheim2003Oxygen consumption and body temperature of active and resting honeybees.J Insect Physiol49881889

31. Kovac H, Stabentheiner A, Hetz SK, Petz M, Crailsheim K (2007) Respiration of resting honeybees. J Insect Physiol 53: 1250–1261.H. KovacA. StabentheinerSK HetzM. PetzK. Crailsheim2007Respiration of resting honeybees.J Insect Physiol5312501261

32. Harrison JM (1987) Roles of individual honeybee workers and drones in colonial thermogenesis. J Exp Biol 129: 53–61.JM Harrison1987Roles of individual honeybee workers and drones in colonial thermogenesis.J Exp Biol1295361

33. Stabentheiner A, Kovac H (2002) Beitrag unterschiedlich alter Arbeiterinnen zur Wärmeproduktion im Brutnest von Bienenvölkern. Apidologie 33: 499–500.A. StabentheinerH. Kovac2002Beitrag unterschiedlich alter Arbeiterinnen zur Wärmeproduktion im Brutnest von Bienenvölkern.Apidologie33499500

34. Basile R, Pirk CWW, Tautz J (2008) Trophallactic activities in the honeybee brood nest - Heaters get supplied with high performance fuel. Zoology 111: 433–441.R. BasileCWW PirkJ. Tautz2008Trophallactic activities in the honeybee brood nest - Heaters get supplied with high performance fuel.Zoology111433441

35. Fehler M, Kleinhenz M, Klügl F, Puppe F, Tautz J (2007) Caps and gaps: a computer model for studies on brood incubation strategies in honeybees (Apis mellifera carnica). Naturwissenschaften 94: 675–680.M. FehlerM. KleinhenzF. KlüglF. PuppeJ. Tautz2007Caps and gaps: a computer model for studies on brood incubation strategies in honeybees (Apis mellifera carnica).Naturwissenschaften94675680

36. Lindauer M (1952) Ein Beitrag zur Frage der Arbeitsteilung im Bienenstaat. Z Vgl Physiol 34: 299–345.M. Lindauer1952Ein Beitrag zur Frage der Arbeitsteilung im Bienenstaat.Z Vgl Physiol34299345

37. Seeley TD (1982) Adaptive significance of the age polyethism schedule in honeybee colonies. Behav Ecol Sociobiol 11: 287–293.TD Seeley1982Adaptive significance of the age polyethism schedule in honeybee colonies.Behav Ecol Sociobiol11287293

38. Seeley TD (1995) The wisdom of the hive: the social physiology of honeybee colonies. Cambridge: Havard University Press. TD Seeley1995The wisdom of the hive: the social physiology of honeybee colonies.CambridgeHavard University Press295 p

39. Seeley TD, Kolmes SA (1991) Age polyethism for hive duties in honey bees - illusion or reality? Ethology 87: 284–297.TD SeeleySA Kolmes1991Age polyethism for hive duties in honey bees - illusion or reality?Ethology87284297

40. Page RE, Robinson GE (1991) The genetics of division of labor in honey bee colonies. Adv Insect Physiol 23: 117–169.RE PageGE Robinson1991The genetics of division of labor in honey bee colonies.Adv Insect Physiol23117169

41. Robinson GE (1992) Regulation of division of labor in insect societies. Ann Rev Entomol 37: 637–665.GE Robinson1992Regulation of division of labor in insect societies.Ann Rev Entomol37637665

42. Johnson BR (2002) Reallocation of labor in honeybee colonies during heat stress: the relative roles of task switching and the activation of reserve labor. Behav Ecol Sociobiol 51: 188–196.BR Johnson2002Reallocation of labor in honeybee colonies during heat stress: the relative roles of task switching and the activation of reserve labor.Behav Ecol Sociobiol51188196

43. Johnson BR (2008) Within-nest temporal polyethism in the honey bee. Behav Ecol Sociobiol 62: 777–784.BR Johnson2008Within-nest temporal polyethism in the honey bee.Behav Ecol Sociobiol62777784

44. Stabentheiner A, Schmaranzer S (1987) Thermographic Determination of Body Temperatures in Honey Bees and Hornets: Calibration and Applications. Thermology 2: 563–572.A. StabentheinerS. Schmaranzer1987Thermographic Determination of Body Temperatures in Honey Bees and Hornets: Calibration and Applications.Thermology2563572

45. Sachs L (1997) Angewandte Statistik. Berlin, Heidelberg, New York: Springer. L. Sachs1997Angewandte Statistik.Berlin, Heidelberg, New YorkSpringer889 p

46. Stabentheiner A, Kovac H, Schmaranzer S (2002) Honeybee nestmate recognition: the thermal behaviour of guards and their examinees. J Exp Biol 205: 2637–2642.A. StabentheinerH. KovacS. Schmaranzer2002Honeybee nestmate recognition: the thermal behaviour of guards and their examinees.J Exp Biol20526372642

47. Stabentheiner A, Kovac H, Schmaranzer S (2007) Thermal behaviour of honeybees during aggressive interactions. Ethology 113: 995–1006.A. StabentheinerH. KovacS. Schmaranzer2007Thermal behaviour of honeybees during aggressive interactions.Ethology1139951006

48. Esch H (1960) Über die Körpertemperaturen und den Wärmehaushalt von Apis mellifica. Z Vgl Physiol 43: 305–335.H. Esch1960Über die Körpertemperaturen und den Wärmehaushalt von Apis mellifica.Z Vgl Physiol43305335

49. Mandl M (2005) Thermoregulation im brütenden Bienenvolk. Diploma thesis, Karl-Franzens-Universität Graz, Austria. M. Mandl2005Thermoregulation im brütenden Bienenvolk.Diploma thesis, Karl-Franzens-Universität Graz, Austria

50. Stabentheiner A, Hagmüller K (1991) Sweet food means ‘hot dancing’ in honey bees. Naturwissenschaften 78: 471–473.A. StabentheinerK. Hagmüller1991Sweet food means ‘hot dancing’ in honey bees.Naturwissenschaften78471473

51. Stabentheiner A, Kovac H, Hagmüller K (1995) Thermal behavior of round and wagtail dancing honeybees. J Comp Physiol B 165: 433–444.A. StabentheinerH. KovacK. Hagmüller1995Thermal behavior of round and wagtail dancing honeybees.J Comp Physiol B165433444

52. Stabentheiner A (1996) Effect of foraging distance on the thermal behavior of honeybees during dancing, walking and trophallaxis. Ethology 102: 360–370.A. Stabentheiner1996Effect of foraging distance on the thermal behavior of honeybees during dancing, walking and trophallaxis.Ethology102360370

53. Stabentheiner A (2001) Thermoregulation of dancing bees: thoracic temperature of pollen and nectar foragers in relation to profitability of foraging and colony need. J Insect Physiol 47: 385–392.A. Stabentheiner2001Thermoregulation of dancing bees: thoracic temperature of pollen and nectar foragers in relation to profitability of foraging and colony need.J Insect Physiol47385392

54. Seeley TD (1992) The tremble dance of the honey bee: message and meanings. Behav Ecol Sociobiol 31: 375–383.TD Seeley1992The tremble dance of the honey bee: message and meanings.Behav Ecol Sociobiol31375383

55. Seeley TD (1994) Honey bee foragers as sensory units of their colonies. Behav Ecol Sociobiol 34: 51–62.TD Seeley1994Honey bee foragers as sensory units of their colonies.Behav Ecol Sociobiol345162

56. Klein BA, Olzsowsky KM, Klein A, Saunders KM, Seeley TD (2008) Caste-dependent sleep of worker honey bees. J Exp Biol 211: 3028–3040.BA KleinKM OlzsowskyA. KleinKM SaundersTD Seeley2008Caste-dependent sleep of worker honey bees.J Exp Biol21130283040

57. Kovac H, Stabentheiner A, Brodschneider R (2009) Contribution of honeybee drones of different age to colonial thermoregulation. Apidologie 40: 82–95.H. KovacA. StabentheinerR. Brodschneider2009Contribution of honeybee drones of different age to colonial thermoregulation.Apidologie408295

58. Graham S, Myerscough MR, Jones JC, Oldroyd BP (2006) Modelling the role of intracolonial genetic diversity on regulation of brood temperature in honey bee (Apis mellifera L.) colonies. Insect Soc 53: 226–232.S. GrahamMR MyerscoughJC JonesBP Oldroyd2006Modelling the role of intracolonial genetic diversity on regulation of brood temperature in honey bee (Apis mellifera L.) colonies.Insect Soc53226232

59. Stabentheiner A, Crailsheim K (1999) The effect of activity level and ambient temperature on thermoregulation in isolated honeybees (Hymenoptera: Apidae). Entomol Gener 24: 13–21.A. StabentheinerK. Crailsheim1999The effect of activity level and ambient temperature on thermoregulation in isolated honeybees (Hymenoptera: Apidae).Entomol Gener241321

60. Crailsheim K, Stabentheiner A, Hrassnigg N, Leonhard B (1999) Oxygen consumption at different activity levels and ambient temperatures in isolated honeybees (Hymenoptera: Apidae). Entomol Gener 24: 1–12.K. CrailsheimA. StabentheinerN. HrassniggB. Leonhard1999Oxygen consumption at different activity levels and ambient temperatures in isolated honeybees (Hymenoptera: Apidae).Entomol Gener24112

61. Rothe U, Nachtigall W (1989) Flight of the honeybee, IV. Respiratory quotients and metabolic rates during sitting, walking and flying. J Comp Physiol B 158: 739–749.U. RotheW. Nachtigall1989Flight of the honeybee, IV. Respiratory quotients and metabolic rates during sitting, walking and flying.J Comp Physiol B158739749

62. Vollmann J, Stabentheiner A, Kovac H (2004) Die Entwicklung der Endothermie bei Honigbienen (Apis mellifera carnica Pollm.). Mitt dtsch Ges allg angew Ent 14: 467–470.J. VollmannA. StabentheinerH. Kovac2004Die Entwicklung der Endothermie bei Honigbienen (Apis mellifera carnica Pollm.).Mitt dtsch Ges allg angew Ent14467470

63. Heran H (1952) Untersuchungen über den Temperatursinn der Honigbiene unter besonderer Berücksichtigung der Wahrnehmung strahlender Wärme. Z Vgl Physiol 34: 179–206.H. Heran1952Untersuchungen über den Temperatursinn der Honigbiene unter besonderer Berücksichtigung der Wahrnehmung strahlender Wärme.Z Vgl Physiol34179206

64. Cahill K, Lustick S (1976) Oxygen consumption and thermoregulation in Apis mellifera workers and drones. Comp Biochem Physiol 55A: 355–357.K. CahillS. Lustick1976Oxygen consumption and thermoregulation in Apis mellifera workers and drones.Comp Biochem Physiol55A355357

65. Starks PT, Gilley DC (1999) Heat shielding: a novel method of colonial thermoregulation in honey bees. Naturwissenschaften 86: 438–440.PT StarksDC Gilley1999Heat shielding: a novel method of colonial thermoregulation in honey bees.Naturwissenschaften86438440

66. Starks PT, Johnson RN, Siegel AJ, Decelle MM (2005) Heat shielding: a task for youngsters. Behav Ecol 16: 128–132.PT StarksRN JohnsonAJ SiegelMM Decelle2005Heat shielding: a task for youngsters.Behav Ecol16128132

67. Sakagami SF, Maruyama K (1956) Aktivitätsänderungen der Adenosintriposphatase im Flügelmuskel der Bienenarbeiterinnen während ihres ersten Imaginallebens. Z Vgl Physiol 39: 21–24.SF SakagamiK. Maruyama1956Aktivitätsänderungen der Adenosintriposphatase im Flügelmuskel der Bienenarbeiterinnen während ihres ersten Imaginallebens.Z Vgl Physiol392124

68. Herold RC, Borei H (1963) Cytochrome changes during honeybee flight muscle development. Dev Biol 8: 67–79.RC HeroldH. Borei1963Cytochrome changes during honeybee flight muscle development.Dev Biol86779

69. Herold RC (1965) Development and Ultrastructural Changes of Sarcosomes during Honey Bee Flight Muscle Development. Dev Biol 12: 269–286.RC Herold1965Development and Ultrastructural Changes of Sarcosomes during Honey Bee Flight Muscle Development.Dev Biol12269286

70. Hersch MI, Crewe RM, Hepburn HR, Thompson P, Savage N (1978) Sequential development of glycolytic competence in the muscles of worker honeybees. Comp Biochem Physiol 61B: 427–431.MI HerschRM CreweHR HepburnP. ThompsonN. Savage1978Sequential development of glycolytic competence in the muscles of worker honeybees.Comp Biochem Physiol61B427431

71. Moritz RFA (1988) Biochemical changes during honey bee flight muscle development. Biona-report 6: 51–64.RFA Moritz1988Biochemical changes during honey bee flight muscle development.Biona-report65164