Bridging the Gap between Field Experiments and Machine Learning: The EC H2020 B-GOOD Project as a Case Study towards Automated Predictive Health Monitoring of Honey Bee Colonies.

Honey bee colonies have great societal and economic importance. The main challenge that beekeepers face is keeping bee colonies healthy under ever-changing environmental conditions. In the past two decades, beekeepers that manage colonies of Western honey bees (Apis mellifera) have become increasingly concerned by the presence of parasites and pathogens affecting the bees, the reduction in pollen and nectar availability, and the colonies' exposure to pesticides, among others. Hence, beekeepers need to know the health condition of their colonies and how to keep them alive and thriving, which creates a need for a new holistic data collection method to harmonize the flow of information from various sources that can be linked at the colony level for different health determinants, such as bee colony, environmental, socioeconomic, and genetic statuses. For this purpose, we have developed and implemented the B-GOOD (Giving Beekeeping Guidance by computational-assisted Decision Making) project as a case study to categorize the colony's health condition and find a Health Status Index (HSI). Using a 3-tier setup guided by work plans and standardized protocols, we have collected data from inside the colonies (amount of brood, disease load, honey harvest, etc.) and from their environment (floral resource availability). Most of the project's data was automatically collected by the BEEP Base Sensor System. This continuous stream of data served as the basis to determine and validate an algorithm to calculate the HSI using machine learning. In this article, we share our insights on this holistic methodology and also highlight the importance of using a standardized data language to increase the compatibility between different current and future studies. We argue that the combined management of big data will be an essential building block in the development of targeted guidance for beekeepers and for the future of sustainable beekeeping.

Diverse pollen nutrition can improve the development of solitary bees but does not mitigate negative pesticide impacts.

Floral resource loss and pesticide exposure are major threats to bees in intensively managed agroecosystems, but interactions among these drivers remain poorly understood. Altered composition and lowered diversity of pollen nutrition may reinforce negative pesticide impacts on bees. Here we investigated the development and survival of the solitary bee Osmia bicornis provisioned with three different pollen types, as well as a mixture of these types representing a higher pollen diversity. We exposed bees of each nutritional treatment to five pesticides at different concentrations in the laboratory. Two field-realistic concentrations of three nicotinic acetylcholine receptor (nAChR) modulating insecticides (thiacloprid, sulfoxaflor and flupyradifurone), as well as of two fungicides (azoxystrobin and tebuconazole) were examined. We further measured the expression of two detoxification genes (CYP9BU1, CYP9BU2) under exposure to thiacloprid across different nutrition treatments as a potential mechanistic pathway driving pesticide-nutrition interactions. We found that more diverse pollen nutrition reduced development time, enhanced pollen efficacy (cocoon weight divided by consumed pollen weight) and pollen consumption, and increased weight of O. bicornis after larval development (cocoon weight). Contrary to fungicides, high field-realistic concentrations of all three insecticides negatively affected O. bicornis by extending development times. Moreover, sulfoxaflor and flupyradifurone also reduced pollen efficacy and cocoon weight, and sulfoxaflor reduced pollen consumption and increased mortality. The expression of detoxification genes differed across pollen nutrition types, but was not enhanced after exposure to thiacloprid. Our findings highlight that lowered diversity of pollen nutrition and high field-realistic exposure to nAChR modulating insecticides negatively affected the development of O. bicornis, but we found no mitigation of negative pesticide impacts through increased pollen diversity. These results have important implications for risk assessment for bee pollinators, indicating that negative effects of nAChR modulating insecticides to developing solitary bees are currently underestimated.

Real-time monitoring of honeybee colony daily activity and bee loss rates can highlight the risk posed by a pesticide.

Information on honeybee foraging performance and especially bee loss rates at the colony level are crucial for evaluating the magnitude of effects due to pesticide exposure, thereby ensuring that protection goals for honeybee colonies are met (i.e. threshold of acceptable effects). However, current methods for monitoring honeybee foraging activity and mortality are very approximate (visual records) or are time-limited and mostly based on single cohort analysis. We therefore assess the potential of bee counters, that enable a colony-level and continuous monitoring of bee flight activity and mortality, in pesticide risk assessment. After assessing the background activity and bee loss rates, we exposed colonies to two concentrations of sulfoxaflor (a neurotoxic insecticide) in sugar syrup: a concentration that was considered to be field realistic (0.59 µg/ml) and a higher concentration (2.36 µg/ml) representing a worst-case exposure scenario. We did not find any effect of the field-realistic concentration on flight activity and bee loss rates. However, a two-fold decrease in daily flight activity and a 10-fold increase in daily bee losses were detected in colonies exposed to the highest sulfoxaflor concentration as compared to before exposure. When compared to the theoretical trigger values associated with the specific protection goal of 7 % colony-size reduction, the observed fold changes in daily bee losses were often found to be at risk for colonies. In conclusion, the real-time and colony-level monitoring of bee loss rates, combined with threshold values indicating at which levels bee loss rates threaten the colony, have great potential for improving regulatory pesticide risk assessments for honeybees under field conditions.

Using physiology to better support wild bee conservation.

There is accumulating evidence that wild bees are experiencing a decline in terms of species diversity, abundance or distribution, which leads to major concerns about the sustainability of both pollination services and intrinsic biodiversity. There is therefore an urgent need to better understand the drivers of their decline, as well as design conservation strategies. In this context, the current approach consists of linking observed occurrence and distribution data of species to environmental features. While useful, a highly complementary approach would be the use of new biological metrics that can link individual bee responses to environmental alteration with population-level responses, which could communicate the actual bee sensitivity to environmental changes and act as early warning signals of bee population decline or sustainability. We discuss here through several examples how the measurement of bee physiological traits or performance can play this role not only in better assessing the impact of anthropogenic pressures on bees, but also in guiding conservation practices with the help of the documentation of species' physiological needs. Last but not least, because physiological changes generally occur well in advance of demographic changes, we argue that physiological traits can help in predicting and anticipating future population trends, which would represent a more proactive approach to conservation. In conclusion, we believe that future efforts to combine physiological, ecological and population-level knowledge will provide meaningful contributions to wild bee conservation-based research.

Nutritional stress exacerbates impact of a novel insecticide on solitary bees' behaviour, reproduction and survival.

Pesticide exposure and food stress are major threats to bees, but their potential synergistic impacts under field-realistic conditions remain poorly understood and are not considered in current pesticide risk assessments. We conducted a semi-field experiment to examine the single and interactive effects of the novel insecticide flupyradifurone (FPF) and nutritional stress on fitness proxies in the solitary bee Osmia bicornis. Individually marked bees were released into flight cages with monocultures of buckwheat, wild mustard or purple tansy, which were assigned to an insecticide treatment (FPF or control) in a crossed design. Nutritional stress, which was high in bees foraging on buckwheat, intermediate on wild mustard and low on purple tansy, modulated the impact of insecticide exposure. Within the first day after application of FPF, mortality of bees feeding on buckwheat was 29 times higher compared with control treatments, while mortality of FPF exposed and control bees was similar in the other two plant species. Moreover, we found negative synergistic impacts of FPF and nutritional stress on offspring production, flight activity, flight duration and flower visitation frequency. These results reveal that environmental policies and risk assessment schemes that ignore interactions among anthropogenic stressors will fail to adequately protect bees and the pollination services they provide.

Pesticide risk assessment: honeybee workers are not all equal regarding the risk posed by exposure to pesticides.

Toxicological studies in honeybees have long shown that a single pesticide dose or concentration does not necessarily induce a single response. Inter-individual differences in pesticide sensitivity and/or the level of exposure (e.g., ingestion of pesticide-contaminated matrices) may explain this variability in risk posed by a pesticide. Therefore, to better inform pesticide risk assessment for honeybees, we studied the risk posed by pesticides to two behavioral castes, nurse, and forager bees, which are largely represented within colonies and which exhibit large differences in their physiological backgrounds. For that purpose, we determined the sensitivity of nurses and foragers to azoxystrobin (fungicide) and sulfoxaflor (insecticide) upon acute or chronic exposure. Azoxystrobin was found to be weakly toxic to both types of bees. However, foragers were more sensitive to sulfoxaflor than nurses upon acute and chronic exposure. This phenomenon was not explained by better sulfoxaflor metabolization in nurses, but rather by differences in body weight (nurses being 1.6 times heavier than foragers). Foragers consistently consumed more sugar syrup than nurses, and this increased consumption was even more pronounced with pesticide-contaminated syrup (at specific concentrations). Altogether, the stronger susceptibility and exposure of foragers to sulfoxaflor contributed to increases of 2 and tenfold for the acute and chronic risk quotients, respectively, compared to nurses. In conclusion, to increase the safety margin and avoid an under-estimation of the risk posed by insecticides to honeybees, we recommend systematically including forager bees in regulatory tests.

Warmer winters are associated with lower levels of the cryoprotectant glycerol, a slower decrease in vitellogenin expression and reduced virus infections in winter honeybees.

Within the context of climate change, winter temperatures at high latitudes are predicted to rise faster than summer temperatures. This phenomenon is expected to negatively affect the diapause performance and survival of insects, since they largely rely on low temperatures to lower their metabolism and preserve energy. However, some insects like honeybees, remain relatively active during the winter and elevate their metabolic rate to produce endothermic heat when temperatures drop. Warming winters are thus expected to improve overwintering performance of honeybees. In order to verify this hypothesis, for two consecutive years, we exposed honeybee colonies to either a mild or cold winter. We then monitored the influence of wintering conditions on several parameters of honeybee overwintering physiology, such as levels of the cryoprotectant glycerol, expression levels of immune and antioxidant genes, and genes encoding multifunctional proteins, including vitellogenin, which promotes bee longevity. Winter conditions had no effect on the expression of antioxidant genes, and genes related to immunity were not consistently affected. However, mild winters were consistently associated with a lower investment in glycerol synthesis and a higher expression of fat body genes, especially apidaecin and vitellogenin. Finally, while we found that viral loads generally decreased through the winter, this trend was more pronounced under mild winter conditions. In conclusion, and without considering how warming temperatures might affect other aspects of honeybee biology before overwintering, our data suggest that warming temperatures will likely benefit honeybee vitality by notably reducing their viral loads over the winter.

Delayed effects of a single dose of a neurotoxic pesticide (sulfoxaflor) on honeybee foraging activity.

Pesticide risk-assessment guidelines for honeybees (Apis mellifera) generally require determining the acute toxicity of a chemical over the short-term through fix-duration tests. However, potential long-lasting or delayed effects resulting from an acute exposure (e.g. a single dose) are often overlooked, although the modification of a developmental process may have life-long consequences. To investigate this question, we exposed young honeybee workers to a single sublethal field-realistic dose of a neurotoxic pesticide, sulfoxaflor, at one of two amounts (16 or 60 ng), at the moment when they initiated orientation flights (preceding foraging activity). We then tracked in the field their flight activity and lifespan with automated life-long monitoring devices. Both amounts of sulfoxaflor administered reduced the total number of flights but did not affect bee survival and flight duration. When looking at the time series of flight activity, effects were not immediate but delayed until foraging activity with a decrease in the daily number of foraging flights and consequently in their total number (24 and 33% less for the 16 and 60 ng doses, respectively). The results of our study therefore blur the general assumption in honeybee toxicology that acute exposure results in immediate and rapid effects and call for long-term recording and/or time-to-effect measurements, even upon exposure to a single dose of pesticide.

Pollen nutrition fosters honeybee tolerance to pesticides.

A reduction in floral resource abundance and diversity is generally observed in agro-ecosystems, along with widespread exposure to pesticides. Therefore, a better understanding on how the availability and quality of pollen diets can modulate honeybee sensitivity to pesticides is required. For that purpose, we evaluated the toxicity of acute exposure and chronic exposures to field realistic and higher concentrations of azoxystrobin (fungicide) and sulfoxaflor (insecticide) in honeybees provided with pollen diets of differing qualities (named S and BQ pollens). We found that pollen intake reduced the toxicity of the acute doses of pesticides. Contrary to azoxystrobin, chronic exposures to sulfoxaflor increased by 1.5- to 12-fold bee mortality, which was reduced by pollen intake. Most importantly, the risk of death upon exposure to a high concentration of sulfoxaflor was significantly lower for the S pollen diet when compared with the BQ pollen diet. This reduced pesticide toxicity was associated with a higher gene expression of vitellogenin, a glycoprotein that promotes bee longevity, a faster sulfoxaflor metabolization and a lower concentration of the phytochemical p-coumaric acid, known to upregulate detoxification enzymes. Thus, our study revealed that pollen quality can influence the ability of bees to metabolize pesticides and withstand their detrimental effects, providing another strong argument for the restoration of suitable foraging habitat.

Pesticide risk assessment in honeybees: Toward the use of behavioral and reproductive performances as assessment endpoints.

The growing gap between new evidence of pesticide toxicity in honeybees and conventional toxicological assays recommended by regulatory test guidelines emphasizes the need to complement current lethal endpoints with sublethal endpoints. In this context, behavioral and reproductive performances have received growing interest since the 2000s, likely due to their ecological relevance and/or the emergence of new technologies. We review the biological interests and methodological measurements of these predominantly studied endpoints and discuss their possible use in the pesticide risk assessment procedure based on their standardization level, simplicity and ecological relevance. It appears that homing flights and reproduction have great potential for pesticide risk assessment, mainly due to their ecological relevance. If exploratory research studies in ecotoxicology have paved the way toward a better understanding of pesticide toxicity in honeybees, the next objective will then be to translate the most relevant behavioral and reproductive endpoints into regulatory test methods. This will require more comparative studies and improving their ecological relevance. This latter goal may be facilitated by the use of population dynamics models for scaling up the consequences of adverse behavioral and reproductive effects from individuals to colonies.

Interactions Between Thiamethoxam and Deformed Wing Virus Can Drastically Impair Flight Behavior of Honey Bees.

Exposure to multiple stress factors is believed to contribute to honey bee colony decline. However, little is known about how co-exposure to stress factors can alter the survival and behavior of free-living honey bees in colony conditions. We therefore studied the potential interaction between a neonicotinoid pesticide, thiamethoxam, and a highly prevalent honey bee pathogen, Deformed wing virus (DWV). For this purpose, tagged bees were exposed to DWV by feeding or injection, and/or to field-relevant doses of thiamethoxam, then left in colonies equipped with optical bee counters to monitor flight activity. DWV loads and the expression of immune genes were quantified. A reduction in vitellogenin expression level was observed in DWV-injected bees and was associated with precocious onset of foraging. Combined exposure to DWV and thiamethoxam did not result in higher DWV loads compared to bees only exposed to DWV, but induced precocious foraging, increased the risk of not returning to the hive after the first flight, and decreased survival when compared to single stress exposures. We therefore provided the first evidence for deleterious interactions between DWV and thiamethoxam in natural conditions.

Honey bee survival mechanisms against the parasite Varroa destructor: a systematic review of phenotypic and genomic research efforts.

The ectoparasitic mite Varroa destructor is the most significant pathological threat to the western honey bee, Apis mellifera, leading to the death of most colonies if left untreated. An alternative approach to chemical treatments is to selectively enhance heritable honey bee traits of resistance or tolerance to the mite through breeding programs, or select for naturally surviving untreated colonies. We conducted a literature review of all studies documenting traits of A. mellifera populations either selectively bred or naturally selected for resistance and tolerance to mite parasitism. This allowed us to conduct an analysis of the diversity, distribution and importance of the traits in different honey bee populations that can survive V. destructor globally. In a second analysis, we investigated the genetic bases of these different phenotypes by comparing 'omics studies (genomics, transcriptomics, and proteomics) of A. mellifera resistance and tolerance to the parasite. Altogether, this review provides a detailed overview of the current state of the research projects and breeding efforts against the most devastating parasite of A. mellifera. By highlighting the most promising traits of Varroa-surviving bees and our current knowledge on their genetic bases, this work will help direct future research efforts and selection programs to control this pest. Additionally, by comparing the diverse populations of honey bees that exhibit those traits, this review highlights the consequences of anthropogenic and natural selection in the interactions between hosts and parasites.

Honeybee lifespan: the critical role of pre-foraging stage.

Assessing the various anthropogenic pressures imposed on honeybees requires characterizing the patterns and drivers of natural mortality. Using automated lifelong individual monitoring devices, we monitored worker bees in different geographical, seasonal and colony contexts creating a broad range of hive conditions. We measured their life-history traits and notably assessed whether lifespan is influenced by pre-foraging flight experience. Our results show that the age at the first flight and onset of foraging are critical factors that determine, to a large extent, lifespan. Most importantly, our results indicate that a large proportion (40%) of the bees die during pre-foraging stage, and for those surviving, the elapsed time and flight experience between the first flight and the onset of foraging is of paramount importance to maximize the number of days spent foraging. Once in the foraging stage, individuals experience a constant mortality risk of 9% and 36% per hour of foraging and per foraging day, respectively. In conclusion, the pre-foraging stage during which bees perform orientation flights is a critical driver of bee lifespan. We believe these data on the natural mortality risks in honeybee workers will help assess the impact of anthropogenic pressures on bees.

Viruses in the Invasive Hornet Vespa velutina.

The Asian yellow-legged hornet Vespa velutina nigrithorax, a major predator of honeybees, is spreading in Europe in part due to a lack of efficient control methods. In this study, as a first step to identify biological control agents, we characterized viral RNA sequences present in asymptomatic or symptomatic hornets. Among 19 detected viruses, the honey bee virus Deformed wing virus-B was predominant in all the samples, particularly in muscles from the symptomatic hornet, suggesting a putative cause of the deformed wing symptom. Interestingly, two new viruses closely related to Acyrthosiphon pisumvirus and Himetobi Pvirus and viruses typically associated with honey bees, Acute bee paralysis virus and Black queen cell virus, were detected in the brain and muscles, and may correspond to the circulation and possible replication forms of these viruses in the hornet. Aphid lethal paralysis virus, Bee Macula-like virus, and Moku virus, which are known to infect honey bees, were also identified in the gut virus metagenome of hornets. Therefore, our study underlined the urgent need to study the host range of these newly discovered viruses in hornets to determine whether they represent a new threat for honey bees or a hope for the biocontrol of V. velutina.

MALDI-MS Profiling to Address Honey Bee Health Status under Bacterial Challenge through Computational Modeling.

Honey bees play a critical role in the maintenance of plant biodiversity and sustainability of food webs. In the past few decades, bees have been subjected to biotic and abiotic threats causing various colony disorders. Therefore, monitoring solutions to help beekeepers to improve bee health are necessary. Matrix-assisted laser desorption ionization-mass spectrometry (MALDI-MS) profiling has emerged within this decade as a powerful tool to identify in routine micro-organisms and is currently used in real-time clinical diagnosis. MALDI BeeTyping is developed to monitor significant hemolymph molecular changes in honey bees upon infection with a series of entomopathogenic Gram-positive and -negative bacteria. A Serratia marcescens strain isolated from one naturally infected honey bee collected from the field is also considered. A series of hemolymph molecular mass fingerprints is individually recorded and to the authors' knowledge, the first computational model harboring a predictive score of 97.92% and made of nine molecular signatures that discriminate and classify the honey bees' systemic response to the bacteria is built. Hence, the model is challenged by classifying a training set of hemolymphs and an overall recognition of 91.93% is obtained. Through this work, a novel, time and cost saving high-throughput strategy that addresses honey bee health on an individual scale is introduced.

Toward the protection of bees and pollination under global change: present and future perspectives in a challenging applied science.

Over the past 30 years (1987-2016), bibliometric data have shown a drastic change in the scientific investigation of threats to bee populations. Bee research efforts committed to studying bioagressors of honeybees (mainly Varroa sp.) were predominant, but now appear to be shifting from bioagressors to global change in the published literature. This rise of global change science reveals prevailing topics, for current and future years: climate change, landscape alteration, agricultural intensification and invasive species. We argue that with increased investment in applied research and development, the scientific, beekeeping and agricultural communities will be able to find management strategies for productive agrosystems and enhanced resilience of pollination and beekeeping. This implies the need for restoring and improving food resources and shelters of bees by ecological intensification of diversified farming systems, and also reconciling sustainable beekeeping with wild pollinator conservation.

Influence of chronic exposure to thiamethoxam and chronic bee paralysis virus on winter honey bees.

Co-exposure to pesticides and viruses is likely to occur in honey bee colonies. Pesticides can be present in pollen, nectar, and persist in stored food (honey and bee bread), and viruses can be highly prevalent in honey bee colonies. Therefore, the present study describes the influence of chronic co-exposure to thiamethoxam and Chronic bee paralysis virus (CBPV) on bee survival, virus loads, expression level of immune and detoxication genes, and pesticide metabolism Experiments were performed on honey bees collected from a winter apiary with reduced viral contaminations. No synergistic effect of co-exposure was observed on bee survival, nor on the ability of bees to metabolise the pesticide into clothianidin. However, we found that co-exposure caused an increase in CBPV loads that reached the viral levels usually found in overt infections. The effect of co-exposure on CBPV replication was associated with down-regulation of vitellogenin and dorsal-1a gene transcription. Nevertheless, the observed effects might be different to those occurring in spring or summer bees, which are more likelyco-exposed to thiamethoxam and CBPV and exhibit a different physiology.

Exposure to pollen-bound pesticide mixtures induces longer-lived but less efficient honey bees.

Due to the widespread use of pesticides and their persistence in the environment, non-target organisms are chronically exposed to mixtures of toxic residues. Fungicides, herbicides and insecticides are all found at low doses in the diet of pollinators such as honey bees, but due to the lack of data on the toxicological effects of these mixtures, determining their risk is difficult to assess. We therefore developed a study combining the identification of common pollen-bound pesticide mixtures associated with poor colony development and tested their effects on bee behavior and physiology. We exposed bees to the identified pesticide mixtures during the first days of their adult life, a crucial period for physiological development. Using optic bee counters we recorded the behavior of bees throughout their lives and identified two pesticide mixtures that delay the onset of foraging and slow-down foraging activity. Furthermore, one of these mixtures hampers pollen foraging. As bee longevity is strongly influenced by the time spent foraging, bees exposed to these pesticide mixtures outlived control bees. Physiological analysis revealed that perturbations of the energetic metabolism preceded the altered behavior. In conclusion, we found that early-life exposure to low doses of pesticide mixtures can have long-term effects that translate into longer-lived but slower and less efficient bees. These surprising findings contrast with the commonly reported increase in bee mortality upon pesticide exposure, and demonstrate that exposure that may seem harmless (e.g., very low doses, pesticides not intended to kill insects) can have undesirable effects on non-target organisms.

Temperature-driven changes in viral loads in the honey bee Apis mellifera.

Many of the physiological traits in insects are shaped by environmental temperatures, which can influence their interactions with pathogens. Therefore, quantifying the thermal responses of the host-pathogen system is crucial for better understanding and predicting their dynamics due to environmental changes. This is particularly important in honey bees, which are experiencing severe colony losses around the world, notably due to infection with the Deformed wing virus (DWV). To investigate the influence of temperature on the honey bee/DWV relationship we exposed adult bees to low or high temperatures and determined the effects on viral titers and bee survival. Emerging bees naturally infected with DWV were reared in vitro at different temperatures ranging from 15 °C to 37 °C. In addition, some bees reared at 37 °C were exposed daily to acute heat treatments (40 and 43 °C). High temperatures significantly decreased DWV titers close to the initial viral load at emergence but increased bee mortality. The lowest temperature resulted in higher mortality, but virus load was not significantly impacted. In conclusion, our results indicate that temperature could contribute to seasonal variations in viral loads but do not suggest temperature to be used as a tool to eliminate viruses, even given that high temperatures limit viral multiplication.

Measuring biological age to assess colony demographics in honeybees.

Honeybee colonies are increasingly exposed to environmental stress factors, which can lead to their decline or failure. However, there are major gaps in stressor risk assessment due to the difficulty of assessing the honeybee colony state and detecting abnormal events. Since stress factors usually induce a demographic disturbance in the colony (e.g. loss of foragers, early transition from nurse to forager state), we suggest that disturbances could be revealed indirectly by measuring the age- and task-related physiological state of bees, which can be referred to as biological age (an indicator of the changes in physiological state that occur throughout an individual lifespan). We therefore estimated the biological age of bees from the relationship between age and biomarkers of task specialization (vitellogenin and the adipokinetic hormone receptor). This relationship was determined from a calibrated sample set of known-age bees and mathematically modelled for biological age prediction. Then, we determined throughout the foraging season the evolution of the biological age of bees from colonies with low (conventional apiary) or high Varroa destructor infestation rates (organic apiary). We found that the biological age of bees from the conventional apiary progressively decreased from the spring (17 days) to the fall (6 days). However, in colonies from the organic apiary, the population aged from spring (13 days) to summer (18.5 days) and then rejuvenated in the fall (13 days) after Varroa treatment. Biological age was positively correlated with the amount of brood (open and closed cells) in the apiary with low Varroa pressure, and negatively correlated with Varroa infestation level in the apiary with high Varroa pressure. Altogether, these results show that the estimation of biological age is a useful and effective method for assessing colony demographic state and likely detrimental effects of stress factors.