The Nutritional Year-Cycle of Italian Honey Bees ( Apis mellifera ligustica ) in a Southern Temperate Climate

This research work is the part of Ph.D. thesis of Evaluation of the nutritional status of worker honey bees (Apis mellifera ligustica S., 1806) across temporal patterns through morphological analysis. Simple Summary: In this research, the nutritional status of Italian bees from Sassari was monitored over a yearly cycle leading to the first report of how honey bee nutrition varies over time and according to external factors in a warm temperate Mediterranean climate. During spring and summer, the nutritional status of sampled bees changed in parallel with the availability of feed resources: when flowering plants were plentiful bees were in a good nutritional state, and their nutritional state declined when flowers disappeared. During this period, rainfall was of great importance, with summer droughts representing a particularly challenging period for bees in the study area. In fall and winter, honey bee nutrition was in opposition to the availability of feed resources as deteriorating environmental factors and the disappearance of flowering plants caused honey bees to transform into their winter state (called winter bees) with increased individual nutrient storage. Nevertheless, winter bees were only present for a limited time, which was accredited to high winter temperatures and continuous (but limited) availability of flowering plants. These results provide valuable insights into the nutritional dynamics of Italian bees in the Mediterranean that could support management decisions to improve overwintering success and prevent unnecessary losses. Abstract: Nutrition is a key aspect influencing honey bee health and overwintering. Since honey bee seasonality in southern temperate climates represents a significant research gap, this study conducted long-term monitoring of honey bees in the Mediterranean (Sassari, Italy). Specifically, individual weight, fat body, and size measurements (head, thorax, abdomen, and total body) were recorded monthly so to detect changes in the nutrient storage of worker bees during an annual cycle. Data were analysed according to sampling date, climate (temperature, precipitation, and daylength), and flower diversity and were conducted for nurse and forager bees separately. The nutritional honey bee year-cycle generally followed the nectar flow and showed two critical timepoints: summer and winter dearth. A short cessation of activities in late fall/early winter coupled with an increase in nutrient storage indicated the presence of winter bees. Precipitation was found to play an important role in honey bee nutrition in the study area through its impacts on colony demography and plants in particular illustrating how climate change could pose a threat to European honey bee populations in the future. These results provide valuable insights into the nutritional dynamics of Apis mellifera ligustica in the Mediterranean that could support management decisions to improve overwintering success and prevent unnecessary colony losses.


Introduction
Honey bee colony losses represent a grave and yet relatively poorly understood issue in modern apiculture [1][2][3][4][5][6].Regardless of climate, most losses occur in winter, which is a particularly challenging period for these social insects as there is little to no natural forage available [7][8][9].
Historically, as honey bees (Apis mellifera Linnaeus, 1758) spread from tropical /subtropical regions to temperate climates of the northern hemisphere they evolved unique adaptations allowing colonies to bridge harsh winter conditions without entering a dormant state [10][11][12][13].Specifically, honey bees synchronized their activities with plant phenology, greatly reduced brood rearing in winter, and assumed the formation of a thermoregulating cluster during the coldest months [7,8,14].Furthermore, A. mellifera adopted significant seasonal changes in individual lifespan within its yearly cycle.This has led to the description of two temporally distinct worker bee types; while the honey bee workforce is made up of classical short-lived "summer bees" during most of the year, in winter, these bees are replaced by long-lived winter or diutinus bees [15].
Besides assuring colony survival through thermoregulation [7,16] winter bees effectively function as a "nutrient storage caste" [12,17,18].These bees store large amounts of fat and protein within their bodies (through the accumulation of vitellogenin: Vg) which are conserved throughout winter and subsequently utilized to reinitiate brood rearing when the return of favourable environmental conditions is anticipated [8,19].Moreover, it is this same Vg that grants winter bees their longevity [17,[20][21][22], illustrating the fundamental role of nutrition for the survival of cold-adapted honey bees.Other typical features of diutinus bees (hypertrophied hypopharyngeal glands, enlarged fat bodies, and elevated hemolymph protein content) are also related to nutrition [7,8,14,17,[23][24][25].
Extensive research has allowed for the description of an elegant system showing how honey bees in temperate zones have adapted mechanisms of age division of labour into a bimodal, biannual worker caste system governed by a multitude of internal and external factors with varying sensitivity (reviewed in [15]).In brief, deteriorating environmental conditions and the disappearance of nutrient resources (nectar and pollen) likely cause a drastic reduction in brood rearing, triggering the transition of newly emerging bees into diutinus bees.It is noteworthy that this seasonal shift is mainly linked to the dwindling pollen availability in fall rather than to fluctuations in meteorological factors offering temporal plasticity and adaptability in a changing climate [23].
Whereas overwintering of honey bees in northern regions have been well studied [7,8,15], much less is known regarding the seasonal dynamics of these insects in southern temperate climates.At these latitudes, warm summers and soft winters generally allow for a long foraging season and only a short cessation of activities in winter [26][27][28][29].While this seems advantageous, relatively high winter temperatures can lead to unsustainable brood rearing causing exhaustion of worker bees towards spring [16,26].Moreover, extended periods of foraging resource dearth (e.g., during summer droughts) can put nutritional stress on a colony, hampering preparations for winter [28,[30][31][32].
Lastly, research efforts in relation to seasonal adaptations of honey bees have mainly focused on northern subspecies (e.g., Apis mellifera mellifera) [13].Since a higher survival rate of locally adapted subspecies [52] as well as adaptation to specific climatic conditions [53][54][55] has been shown, knowledge of southern honey bee populations is of increasing interest, especially in the face of accelerated climate change [30,56].
Against this background, this study aimed to conduct long-term monitoring of the nutritional status of locally adapted Italian honey bees (Apis mellifera ligustica, Spinola 1806) in a Mediterranean climate (Sassari, Italy).The goal of this research was to provide a better understanding of the activity and nutritional status of worker (both nurse and forager) bees in southern temperate climates and to generate new insights on the dynamics of the summer and winter bee transition in correlation with seasonal changes in environmental factors and feed resource availability.In addition, authors aimed to provide novel knowledge regarding the seasonal dynamics of Italian honey bees specifically, and the possible challenges these bees face in a changing climate.

Study Site and Apiary
Monitoring was conducted between February 2022 and January 2023 (12 months).A total of n. 5 colonies of Italian honey bees (Apis mellifera ligustica Spinola, 1806), located in a private apiary in the province of Sassari (Sardinia, Italy; 40 • 37 ′ 14.5 ′′ N 8 • 20 ′ 43.1 ′′ E), were studied.The southern temperate Mediterranean climate of the study area, with hot dry summers and mild wet winters, typically allows for a long foraging season and only a short cessation of activities in winter.The initiation of the study was planned according to the seasonal pattern of Italian honey bees in the region, coinciding with the start of the foraging season.
Meteorological data over the course of the study, including mean monthly temperature, precipitation, the number of days with precipitation, relative humidity, windspeed, and daylength (hours of daylight), as retrieved from the weather station of the meteorological services of the Military Airforce of ENAV (Ente Nazionale Assistenza al Volo) located approximately 15 km from the study area are summarized in Table 1.Hives were located in a semi-natural agricultural area, surrounded by managed and unmanaged fields, vineyards, olive groves, and small-scale mixed agriculture (vegetable gardens).The botanical composition of spontaneous flora and the phenological state of plants, with particular regard to pollen availability, in the direct vicinity of the apiary, was monitored throughout the sampling period.Specifically, in order to assess the diversity of flowering plants, 3 100 m × 2 m transects were defined prior to the initiation of the study.At each sampling date (concomitantly with bee monitoring, as described below), the 3 transects were walked by a single observer and the various species of flowering plants known to be visited by honey bees recorded.
Honey bee samples were collected from 5 individual hives selected by the responsible apiarist based on overall health and uniformity.Selected colonies were separated from the rest of the apiary by a distance of 25 m before the initiation of the study.Colonies received standardized care during the study period and were inspected weekly insuring good health.No clinical signs of disease were noted during the course of this study.Treatment against the ectoparasitic mite Varroa destructor was applied in March, August, and December using Amitraz and oxalic acid.Colonies were fed a homemade sucrose solution (3:2 sucrose/water) in spring and fall.No other nutritional supplements were provided.
Apis m. ligustica queens with a nucleus were acquired from a commercial queen breeder and introduced to each respective colony the year before.Queens remained during the whole duration of the study.Colonies were housed in wooden Dadant-type hives with 10 commercial brood frames with a cell diameter of 5.4 mm.
Brood rearing patterns were consistent with that of Italian bees in a southern temperate climate (exhibiting a "Mediterranean pattern" as has been described for bees on the neighbouring island of Corsica [57]); showing a steady increase from spring until peaking in June and subsequently decreasing during hot summer months.A second minor peak was seen in early fall.Three out of the five hives showed a cessation of brood rearing (for approximately 2 weeks) in early December and brood rearing remained relatively low until spring.No foraging stop was observed for any of the hives.

Sample Collection
Ten forager and ten nurse bees from each hive were collected separately on the last week of each month (100 individuals; 50 foragers/50 nurse bees per month).(1) Foragers: bees returning to the hive were captured from the flight deck using a horsehair brush.
(2) Nurse bees: young bees from the centre of the brood nest were collected.Captured bees were stored in 250 cl glass containers with breathable fabric lids and transported to the laboratory of animal production and nutrition of the university of Sassari (UNISS) in a cooler box with icepacks.Individual bees were weighed using a digital scale (OHAUS Pioneer Corp. Las Vegas (Nevada) Pioneer USA corporation, LA, mod.PA512C; precision of 0.01 g) before being frozen (−18 • C) and stored in 1.5 mL microcentrifuge tubes until further analysis.Any pollen or visible attachments were removed manually prior to weighing.

Morphological Analysis
Sampled bees were analysed in their entirety and within a frozen state insuring correct proportional morphological retention.Using a digital calliper (precision 0.01 mm) under a stereomicroscope (Leica ® EZ4 HD), six size measurements were taken for each individual bee; (1) Head width (HW), (2) thoracal width (TW), (3) thoracal length (TL), (4) abdominal width (AW), (5) abdominal length (AL), and (6) total body length (T).Width measurements of each respective body part were taken at the widest point.Length measurements of the thorax and abdomen were taken from the anterior end of the protergum to the caudal end of the first abdominal tergum (T 1 -IT; T 1 includes the scutum and scutellum) and the anterior end of the second abdominal tergum to the caudal end of seventh abdominal tergum (IIT-VIIT) not including the stinger, respectively.All size measurements were taken in duplicate and averaged creating a single observation.

Fat Body Quantification
Ether extraction was performed to estimate the weight and relative size of the fat body of bees according to Wilson-Rich et al. [58].Briefly, the abdomen of each bee was severed using surgical scissors and placed into separate holding cups to dry at 25 • C for 3 days.Next, abdomens were placed in individual 1.5 mL microcentrifuge tubes to which 500 µL of diethyl ether was added.Abdomens were removed after 24 h and dried again for 3 days (same conditions).A Binder ED 53 drying oven was used to insure continuity of drying conditions over the duration of the study.Dried abdomens were weighed before and after ether extraction using a ORMA BCA200 electric laboratory balance with a precision 0.0001 g.The fat body weight (FBW) was calculated as the difference between the weight of each abdomen before and after washing with diethyl ether.The relative size of the fat body (FB%) was calculated as the proportional weight of the fat body relative to the weight of the dried abdomens prior to ether extraction [58,59].

Data Analysis
All procedures were carried out using a software package (Minitab statistical software package, Minitab © , New York, NY, USA).Statistical significance was set at p-value < 0.05 and Tukey test was used for the post hoc pairwise comparison of means.

Worker Bee Type
Analysis of variance (ANOVA) was performed to detect significant differences in monitored metrics between the two types of sampled worker bees.A balanced linear model with interaction was used as follows: where y is the dependent variable (n = 9; Weight, HW, TW, TL, AW, AL, T, FBW, FB%), µ is the overall mean, W is the fixed factor representing worker type (2 levels; Forager, Nurse), H is the fixed factor of hive (n = 5; H1, H2, H3, H4, H5), W * H is the interaction term, and ε is the random error.
Further analysis for any dependent variable significantly affected by W was conducted for forager and nurse bees separately.Unaffected variables were analysed using the whole dataset.

Effect of Sampling Date
Analysis of variance (ANOVA) was performed to detect significant differences in monitored metrics between the sampling months.A balanced linear model with interaction was used as follows: where y is the dependent variable (n = 9; Weight, HW, TW, TL, AW, AL, T, FBW, FB%), µ is the overall mean, M is the fixed factor of sampling month (n = 12; January-December), H is the fixed factor of hive (n = 5; H1, H2, H3, H4, H5), M * H is the interaction term, and ε is the random error.If a significant effect of H was found, the dataset was split accordingly, and the effect of M analysed separately.

Effect of Environmental Factors and Flower Diversity
Analysis of variance (ANOVA) was performed to detect any changes in dependent variables of sampled bees according to environmental factors and flower diversity.All data were analysed following a general linear model procedure with interaction as follows: where y is the dependent variable (n = 9; Weight, HW, TW, TL, AW, AL, T, FBW, FB%), µ is the overall mean, T is the fixed factor of monthly average environmental temperature (monthly average ambient temperature was divided into 4 levels; Ta, Tb, Tc, Td; range: 10-15, >15-20, >20-25, >25 • C; Table 1), P is the fixed factor representing mean monthly precipitation (mean monthly precipitation was divided into 3 levels; Pa, Pb, Pc; range: 0-50, >50-100, >100 mm; see Table 1), O is the fixed factor representing monthly average daylength (monthly average daylength was divided into 4 levels; Oa, Ob, Oc, Od; range: 550-650, >650-750, >750-850, >950 h of daylight; see Table 1), F is the fixed factor effect of the monthly flower diversity (monthly flower diversity was dived into 3 levels; Fa, Fb, Fc; range: <5, 5-10, >10 species of flowering plants see Table 2), H is the fixed factor of hive (n = 5; H1, H2, H3, H4, H5), T * H is the interaction term between temperature and hive, P * H between precipitation and hive, O * H between daylight and hive, F * H between plant diversity and hive, and ε is the random error.If a significant effect of H was found, the dataset was split accordingly, and the effect of T/P/F/O analysed separately.

Results
The diversity of flowering plants steeply increased during spring before peaking in early summer.This peak was followed by a drastic decrease over the course of the summer, bottoming in August.Summer dearth was followed by a mild restoration in fall.Limited flower diversity was noted in early winter which increased in February marking the onset of the foraging season.The various species of flowering plants encountered during the study period are reported per sampling month in Table 2.
Post hoc analysis showed H5 to be significantly different from other hives for TW and TL, while no decisive pattern was revealed for T. Mean No significant difference in TW and TL was found over the months for any of the hives.
Results of the analysis of variance on the effect of hive and sampling month for nurse and forager bees are reported in Table 4. Figures 1 and 2 show boxplots of the different variables significantly affected by month for nurse and forager bees respectively.The effect of month on the Weight and T of nurse bees according to hives is depicted in Figures 3 and 4. The effect of month on the AL of forager bees for the different hives is shown in Figure 5. Significantly different months are indicated by a red "*".No significant interaction effect was found.Analysis of variance revealed no significant effect of environmental factors nor flower availability on TW for any of the hives.A significant effect of temperature (F (3,229) = 2.66, p = 0.049) and flower diversity (F (3,229) = 3.44, p = 0.034) on TL was found for H3, and a significant effect of temperature (F (3,229) = 2.66, p = 0.049) for H2.However, post hoc analysis showed no difference in TL between the groupings of various factors.
Results of the analysis of variance and post hoc analysis on the effect of environmental factors (mean monthly temperature, precipitation, daylength) and flower diversity for the dependent variables of nurse and forager bees are reported in Table 5.A significant effect of hive was found for Weight and T in nurse bees and are therefore reported here.Specifically, a significant effect of temperature was found on Weight and T for all hives (Weight: H1: F            Results of the correlation analysis between various metrics are reported for nurse and forager bees in Table 6.

Discussion
Nutrition is a key aspect influencing honey bee health and overwintering success [38,49,61,62].Nevertheless, there are relatively few studies that explore honey bee seasonal activity in southern temperate climates [28,29].In this research, the nutritional status of the Italian bee (A. m. ligustica), a subspecies well adapted to the warm temperate climate of the Mediterranean, was studied [26,55,57].Specifically, individual weight, fat body, and size measurements (head, thorax, abdomen, and total body) were recorded on a monthly basis in order to detect temporal changes in the nutrient storage of worker bees during a complete annual cycle (2022-2023).Recorded parameters were analysed according to climatological factors and the availability of feed recourses (flower diversity) in order to get a better understanding of the annual bimodal dynamics of the honey bee workforce in a southern temperate Mediterranean climate.
Besides following seasonal variations in honey bee nutrition, novel data regarding two distinct worker bee types with varying biological age; in-hive (nurse bees) Vs. out-hive (forager bees) is presented.Given the consistent and fundamental behavioural, physiological, and nutritional differences between these two worker bee types [17,[63][64][65][66][67][68][69] authors hypothesised nutrition-related size metrics to vary significantly between them.Furthermore, as nurse and forager bees have different responses to similar conditions [64], the analysis of fixed factors (sampling time, environmental factors, and feed resource availability) was conducted separately for both cohorts.
While body size is a known indicator of nutritional stress reflecting the quantity and quality of food available during development in honey bees [44,47,[70][71][72][73], to the best of our knowledge, no empirical evidence has so far been produced showing size variations between worker honey bees to be related to age division of labour.In fact, body size variations of worker bees within a single A. mellifera colony are believed to be negligible [74][75][76].Here we show significant differences in nutrition-related size measurements between individual forager and nurse bees.With the exception of thoracal dimensions, all measured metrics differed between both worker bee types.Correspondingly, individual size measurements were strongly or mildly positively correlated to known biological markers of honey bee nutrition (body and fat body weight [3,25,44,77]) (Table 6).These findings are in accordance with worker physiology, showing nurse bees to have substantially larger nutrient stores as compared to foragers [17,[62][63][64][65][67][68][69].The weak correlation between Weight and FBW Vs.HW in foragers is explained by the fact that these bees have hypotrophied hypopharyngeal glands [7,17].In contrast, these glands, which serve for the production of brood food, are well-developed in nurse bees [7,17,[78][79][80][81].Since brood food is produced from Vg [18], it is logical HW to be correlated to nutritional markers in this cohort.
The long-term monitoring of selected metrics allowed us to paint a detailed picture of the annual cycle of Italian bees in the study area from a nutritional point of view.In accordance with the seasonal adaptations of worker honey bees (summer Vs. winter bees) [7,8,15,25,29,82], a functional bimodal division of the honey bee cycle is followed.
The "summer-bee portion" of the nutritional cycle, running from mid-winter (end of December) to early fall (September) in this study, closely followed the nectar flow.When feed resources were abundant, individual honey bee nutrition increased and the opposite was seen during resource dearth [28,70,[83][84][85].Contrarily, a general increase in HW and TW was seen over the course of the foraging season.This corresponds to previous findings describing an increase in worker bee size within a yearly cycle [26,86].
Present data reflects a controversial increase in W, AL, T, FBW, and FB% of forager bees during the summer dearth period (peaking in July; Figures 2 and 5) which has been accredited to an explicit sampling error.Specifically, ambient temperatures at the time of sampling were so high that a large portion of bees had exited their hives and were found clustered around their respective hive entrances.This common strategy to prevent overheating [87,88] likely resulted in the sampling of a mixed population of worker bees rather than bees of a single biological age.Our deduction of this finding to be a sampling error is supported by the overlapping ranges of measured metrics between both worker bee types for the month of July, as well as the increased variability seen for that sampling date.
Consistent with our present understanding of honey bee physiology in southern temperate climates of the northern hemisphere [26,28,29,57], the "winter-bee portion" of the nutritional cycle in this study was short and restricted to late fall/early winter (November-December).October can be considered a transition month as the shift from summer to winter bee-state is known to occur gradually within a colony and thus a balanced number of both castes is most likely present at this time [7,8,23].This portion of the honey bee year-cycle was characterized by a steep increase in nutrient storage in opposition to the overall diversity of feed resources (Table 2) with a subsequent decrease over the course of the winter period [17].This correlates well with current knowledge of honey bee seasonality with the arrival of winter bees primarily related to the disappearance of flowering plants in fall [7,8,15].As indicated by the significant difference in average monthly FBW, FB%, and AL (Figure 1), nutrient storage of nurse bees peaked in November, showing the presence of winter bees [13,25,70,81].Whereas the overwintering state of honey bees in warmer climates differs from northern regions (e.g., sustained foraging and brood rearing activities), accumulation of fat and protein (Vg) is believed to be universal for overwintering honey bees in temperate zones [28].Recent research monitoring Vg levels in the fat body of worker bees over a yearly cycle in the Czech Republic revealed a strikingly similar pattern even though nutrient storage in said research peaked in December [25].Nevertheless, because sampling in the present study was conducted at the end of each month, nutrient storage of nurse bees could have peaked early to mid-December (as the noted brood rearing patterns would suggest) rather than in November.It is also necessary to stress on the fact that homemade sucrose solution was offered in negligible amounts, out of the long-term monitoring, because strictly necessary to colony survival and for a very limited period of days (like reported above).
The enlarged Weight, FBW, HW, AW, and T of forager bees in November, indicate that monitored hives exited their winter state somewhere between November and December.Indeed, increased morphological dimensions of forager bees during the winter dearth are indicative of a winter-bee-like state and can be considered remnants of the nutrient accumulation that occurred during in-hive activities [89].Analogously, previous research has identified forager bees with increased morphological dimensions in early spring likely to be winter bees hatched the year before [86].Authors expected to see a delay in the detection of winter-bee-like foragers as compared to nurse bees.Nevertheless, the cessation of activities of Italian bees in this research was shorter than the sampling frequency likely resulting in the absence of a notable temporal divergence in seasonal transition between both worker bee types.
Overall changes in recorded metrics corresponded to the variation in environmental factors observed within the study period known to influence seasonal honey bee colony activity [7,8,15,16,25].
For both nurse and forager bees, monthly average ambient temperatures between 15-20 • C were correlated to the highest degree of nutrient storage (Table 5).These temperatures coincide with peak honey bee activity during the nectar flow in spring as well as the appearance of winter bees in fall.The fact that honey bees show two distinct physiological states within similar temperature ranges illustrates it is unlikely mean temperature alone influences the seasonal transition of honey bee colonies.Alternatively, interaction of temperature with other factors (e.g., photoperiod and feed resource availability), or the direction of temperature change in combination with reaching a threshold value could serve as a possible seasonal trigger [7,14,16,23,81].
Current understanding of honey bee behaviour in temperate climates describes the formation of a thermoregulating cluster when ambient temperatures drop below 10 • C [7,90].While average temperatures were well above this mark in November in the present research (and remained so for the whole duration of the study), minimum ambient temperatures did dip below 10 • C in November.More significantly, temperatures below this threshold were first recorded the month before (October).Given winter bees start appearing during this transition month, the first cold nights in fall could signal colonies to prepare for winter.The physiological mechanisms of how dropping ambient temperatures allow worker bees to accumulate Vg has previously been described [22,25,91,92].
In accordance with previous research efforts [25,62,63,93], an association between decreasing daylength and the accumulation of nutrients in the fat body of in-hive bees was noted.These results strengthen the hypothesis that decreasing photoperiod is involved in the seasonal appearance of winter bees [7,15,16] although present morphometrical results did not reveal further insights into the possible influence of daylength on the nutritional cycle of Italian honey bees.
Lastly, high average monthly precipitation (>100 mm) was consistently associated with an elevated nutrient status in both nurse and forager bees and coincided with the presence of winter bees.With the exception of AW in foragers, all measured metrics were significantly higher during months with high precipitation (Table 5).This finding is intriguing and could point towards weather conditions to be of particular importance in the seasonal dynamics of honey bees in southern temperate climates.Indeed, impaired meteorological conditions ("bad weather") are known to influence honey bee demography by affecting the pheromone balance of a colony resulting in the active suppression of the biological maturation of young bees and the appearance of winter bees in fall [7,15,92,[94][95][96].
The flower diversity surveys conducted in the honey bee flight area provided a significant contribution to the study.We detected substantial variations in forage diversity (with specific regard to pollen availability) over the course of the monitoring period with an explicit pattern matching that of brood rearing.Despite honey bee colonies do not generally prefer to store large amounts of pollen [97] (and when they do, they is mix it with nectar and seal the compounds with wax) the availability of this resource (the main nutrient supply for brood rearing) is chiefly correlated to the brood rearing activity [23,27,31,38,76,98].
This pattern together with rest of the present data allowed us to identify two critical periods for honey bee health and nutrition in southern temperate climates, i.e., summer and winter dearth.While seasonal fluctuation in pollen availability showing one or two distinct peaks is not unusual [26,27], large temporal variations in feed resource availability (nutritional irregularity) are known to affect honey bee health and longevity [84,99].Indeed, poor foraging conditions and related malnutrition are believed to be key factors in global colony losses [36,38,49,85,100], especially in warm temperate climates [29].Sugars from nectar, on one side, and amino-acids and sterols from pollen, on the other side, are differently involved in metabolic patterns of honeybees, in which energy storage is limited in foragers and while being physiologically higher in nurse bees (depending on the period of the year and according to feeding source availability, like we observed).
High winter temperatures [16] together with the prolonged availability of pollen [23] offer a viable explanation for the late appearance of diutinus bees in this research as well as the sustained brood rearing observed for two out of the five hives [29,101].Although this might seem beneficial, continuous brood rearing during periods of limited pollen availability can cause premature exhaustion of fat and protein nutrient stores leaving colonies in a vulnerable state [16,28,29,31,67].In effect, in-hive colony reserves and reserves within bees themselves are rapidly depleted in times of pollen dearth [38].Besides, flower diversity has been shown to be an important factor in honey bee nutrition since different pollen and nectar sources vary significantly in their nutritive value, e.g., protein and mineral contents [35,41,85,99,102].Hence, even though nutritional resources would be available during winter months, the limited variety of flowering plants during this time might not provide adequate nutrition in order to support brood rearing or honey bee colonies in general [28,38,40,41,60,102,103].
A prominent finding of the present study is that the nutritional state of A. m. ligustica workers was significantly negatively affected during periods of high ambient temperatures (>25 • C) and low precipitation (0-50 mm) (Table 5).With the exception of T and HW, all nurse metrics were lowest during the summer drought period (June-August) (Figures 1, 3 and 4).
The precipitation pattern during the study period coincided with that of plant diversity during the "summer-bee portion" of the year which can be considered an illustration of the bottleneck effect of precipitation on plant growth in warm and dry Mediterranean climates [104,105].The noted influence of weather on honey bee nutrition in summer therefore likely stems from an indirect effect on plants resulting in an overall resource dearth [8,30,90,[106][107][108].Moreover, hot and dry conditions have been shown to reduce nectar and pollen production and the overall nutritional quality of these resources as well [30,90,105].For these reasons, in addition to the winter dearth, summer food shortages could be of serious concern for honey bee colonies in southern temperate climates.This could be especially true in the face of accelerated climate change [5,30,31,45,[109][110][111] as conditions in the Mediterranean head towards a similar scenario seen in particularly arid climates such as in the Middle East [111,112] where summer droughts are a key factor in colony losses since many plants suffer from heat stress leading to feed shortage for honey bees [32].
Temporal mismatches with possible nutritional consequences were pointed out at Mediterranean latitudes [16,104,105,113,114]. Indeed, a particularly early initiation of the foraging season, well before the start of the nectar flow, was noted followed by a sharp decrease in nutrient storage over the course of winter (Figure 1).

Conclusions
The present research contributes to our understanding of the seasonal dynamics of honey bees in a southern temperate climate showing a short cessation of activities in late fall/early winter coupled with an increase in nutrient storage of in-hive bees.While the fall decrease in feed resources appears to be the main factor governing honey bee seasonality, a combination of changing environmental factors seems to be required for the arrival of winter bees.The continuous but limited availability of flowering plants and forgiving ambient temperatures during winter likely allowed for the observed brood rearing pattern and consequential sharp decrease in nutrient storage over the winter dearth period.In addition, a first description of the annual nutritional honey bee cycle in a southern temperate climate is presented showing two critical timepoints.Overall, our results contradict the common assumption that warm climates are more suited for honey bees as besides winter, the Mediterranean summer, which is characterised by droughts and high temperatures, was identified as a second critical timepoint.It seems precipitation plays a particularly important role in southern latitudes, influencing nutrition in both the summer-and winter-bee portion of the honey bee year-cycle.Finally, present data seem to support the notion that the shift of environmental conditions have significant effects on honey bees in temperate Europe through a pronounced impact on melliferous plants, indirectly affecting health and nutrition of honeybees.Our results provide valuable insights into the seasonal and nutritional dynamics of locally adapted A. m. ligustica populations that could aid beekeepers to make management decisions in relation to environmental factors and availability of flowering plants with the ultimate goal of improving overwintering success and preventing undesirable colony losses.
January, Feb = February, Mar = March, Apr = April, Jun = June, Ju l= July, Aug = August, Sep = September, Oct = October, Nov = November, and Dec = December.** Data for the month of January were collected in 2023 while the rest of the months regard the year 2022.*** Categories of environmental factors for statistical analysis: monthly flower diversity is dived into 3 levels; Fa, Fb, Fc; range: <5, 5-10, >10 species of flowering plants.

Figure 1 .
Figure 1.Box plots of Head width (mm), Abdominal width (mm), Abdominal length (mm), Fat body weight (mg), and Fat body size (%) of Apis mellifera ligustica nurse bees according to sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represent outliers and the black line represents the mean connect line.Red "*" indicates months significantly different from unmarked months.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.

Figure 1 .
Figure 1.Box plots of Head width (mm), Abdominal width (mm), Abdominal length (mm), Fat body weight (mg), and Fat body size (%) of Apis mellifera ligustica nurse bees according to sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represent outliers and the black line represents the mean connect line.Red "*" indicates months significantly different from unmarked months.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.

Figure 2 .
Figure 2. Box plots of Weight (g), Head width (mm), Abdominal width (mm), Total body length (mm), Fat body weight (mg), and fat body size (%) of Apis mellifera ligustica forager bees according to sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represents outliers and the black line represents the mean connect line.Red "*" indicates months significantly different from unmarked months.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.

Figure 2 .
Figure 2. Box plots of Weight (g), Head width (mm), Abdominal width (mm), Total body length (mm), Fat body weight (mg), and fat body size (%) of Apis mellifera ligustica forager bees according to sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represents outliers and the black line represents the mean connect line.Red "*" indicates months significantly different from unmarked months.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.Agriculture 2024, 14, x FOR PEER REVIEW 10 of 22

Figure 3 .
Figure 3. Box plots of the Weight (g) of Apis mellifera ligustica nurse bees according to hive and sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represent outliers and the black line represents the mean connect line.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.

Figure 3 .
Figure 3. Box plots of the Weight (g) of Apis mellifera ligustica nurse bees according to hive and sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represent outliers and the black line represents the mean connect line.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.

Figure 3 .
Figure 3. Box plots of the Weight (g) of Apis mellifera ligustica nurse bees according to hive and sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represent outliers and the black line represents the mean connect line.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.

Figure 4 .
Figure 4. Box plots of the Total body length (mm) of Apis mellifera ligustica nurse bees according to hive and sampling months.Red "*" indicates months significantly different from unmarked months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represents outliers and the black line represents the mean connect line.Red "*" indicates months significantly different from unmarked months.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.

Figure 4 .
Figure 4. Box plots of the Total body length (mm) of Apis mellifera ligustica nurse bees according to hive and sampling months.Red "*" indicates months significantly different from unmarked months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represents outliers and the black line represents the mean connect line.Red "*" indicates months significantly different from unmarked months.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.Agriculture 2024, 14, x FOR PEER REVIEW 11 of 22

Figure 5 .
Figure 5. Box plots of the Abdominal length (mm) of Apis mellifera ligustica forager bees according to hive and sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represents outliers and the black line represents the mean connect line.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov = November, Dec = December, and Jan = January.Analysis of variance revealed no significant effect of environmental factors nor flower availability on TW for any of the hives.A significant effect of temperature (F(3,229) = 2.66, p = 0.049) and flower diversity (F(3,229) = 3.44, p = 0.034) on TL was found for H3, and a significant effect of temperature (F(3,229) = 2.66, p = 0.049) for H2.However, post hoc analysis showed no difference in TL between the groupings of various factors.

Figure 5 .
Figure 5. Box plots of the Abdominal length (mm) of Apis mellifera ligustica forager bees according to hive and sampling months.The boxplot represents the interquartile range (IQR = Q3 − Q1) and bars represent first (Q1, top) and third quartiles (Q3, bottom) of metric values.Black "*" represents outliers and the black line represents the mean connect line.Feb = February, Mar = March, Apr = April, Jun = June, Jul = July, Aug = August, Sep = September, Oct = October, Nov, =; November, Dec = December, and Jan = January.

Table 1 .
Summarizing meteorological data for the province of Sassari (Italy), February 2022-January 2023.

Table 2 .
Summarizing table of the monthly diversity of flowering plants in the honey bee flight area over the study period (February 2022-January 2023) in Sassari (Italy).
* Overall mean values for the whole dataset (nurse and forager bees together).** Standard Deviation.

Table 4 .
Results of the analysis of variance for the effect of hive and sampling month on various metrics of Apis mellifera ligustica nurse and forager honey bees.
a Statistical significance set at p < 0.005.* Degrees of freedom.

Table 5 .
Results of the analysis of variance for the effects of environmental factors (mean monthly temperature, precipitation, daylength, and flower diversity) on various metrics of Apis mellifera ligustica nurse and forager honey bees.
a Statistical significance set at p < 0.005.* Degrees of freedom.** Categories that were shown to be different through post hoc analysis (Tukey test) are separated by ">", while groupings that are not different from each other are separated by ",".

Table 6 .
Results of Pearson correlation analysis between morphologic metrics of nurse and forager bees.