REVIEW article

Front. Physiol., 17 September 2025

Sec. Invertebrate Physiology

Volume 16 - 2025 | https://doi.org/10.3389/fphys.2025.1623705

Climate change-induced stress in the honey bee Apis mellifera L.- a genetic review

  • 1. Laboratorio de Patología Apícola, Centro de Investigación Apícola y Agroambiental de Marchamalo, (CIAPA), Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal de Castilla-La Mancha (IRIAF), Marchamalo, Spain

  • 2. CREA Centro di Ricerca Agricoltura e Ambiente (CREA-AA), Bologna, Italy

  • 3. Centro de Investigação de Montanha (CIMO), LA SusTEC, Instituto Politécnico de Bragança, Bragança, Portugal

  • 4. Honey Bee Research Laboratory, Agricultural Research Service, United States Department of Agriculture, Beltsville, MD, United States

Article metrics

View details

2

Citations

5k

Views

498

Downloads

Abstract

Climate change is a powerful driver of stress, as it reinforces hotter and drier environments. For bees, the most concerning aspects of these new environmental conditions are the resistance and resilience of bees to changes in temperature, humidity and ultraviolet radiation, as well as the negative effect on diversity of food resources which can lead in nutritional stress. The climatic vulnerability of various bee species and subspecies varies worldwide, as they experience varying levels of stress and display distinct behaviors, weaknesses, and lifespans. To understand these differences, it is crucial to consider both the genetics and epigenetics of bees, as these factors play a key role in their response, resistance, and adaptation to new stressors. This review provides a guide of genetic and epigenetic markers involved in the cellular response of Apis mellifera to most common stressors derived from climate change. Understanding how the various molecular mechanisms interact to restore homeostasis during the stress response is essential for designing future studies based on molecular markers.

1 Introduction

Honey bees are social insects of great ecological and economic importance, but they have experienced substantial losses over the last years due to the action of several interacting biotic and abiotic stressors. These include invasive alien predators, parasites and pathogens, pesticides, and climate change [reviewed by Zhao et al. (2021), Even et al. (2012)]. One of the main impacts of climate change is the intensification of extreme weather events, such as heavy rainfall, prolonged droughts, and heat waves (IPCC-Intergovernmental Panel on Climate Change, 2007), which may severely impact honey bees’ wellbeing. On one hand, these events can indirectly affect honey bee nutrition by disrupting the availability of floral resources and reducing the quantity and quality of nectar and pollen [reviewed by Obeso and Herrera (2018)]. On the other hand, they can directly cause heat stress on honey bees, impairing their foraging efficiency and metabolic functions [reviewed by Zhao et al. (2021)]. Whether honey bees are able to adapt and survive changes in temperature, humidity, and ultraviolet (UV) radiation, and to a reduction in the quantity and diversity of food resources, is a question of utmost importance. However, the vulnerability of bees to climatic vagancies is not necessarily the same worldwide. Even in the same environment, different bee species and subspecies may suffer different stress levels and exhibit different behaviours, weaknesses and lifespans [reviewed by Zhao et al. (2021)]. This would mean that bee species and subspecies may differ in their ability to adapt, resist and be resilient to environmental changes.

Genetic variability is key to resistance and adaptation in stressful environments. Several candidate Single Nucleotide Polymorphisms (SNPs) mapped to genes involved in reproduction, immunity, olfaction, circadian rhythm, lipids biosynthesis and storage were linked to local adaptation in the western honey bee Apis mellifera (Chávez-Galarza et al., 2013; Wallberg et al., 2014; Chen et al., 2016; Henriques et al., 2018). However, beyond genetic determination, complex interactions between genetic and epigenetic factors are known to shape the diversity of organismal phenotypes (West-Eberhard, 2003). A phenotype can change in response to environmental signals without altering its genotype, resulting in modifications to an organism’s physiology and behaviour. This ability is called “phenotypic plasticity” and is remarkably developed in A. mellifera, which can respond to environmental cues to generate dramatically distinct phenotypes, such as queens or workers, through nutritional stimuli such as royal jelly (Winston, 1991; Corona et al., 2016). Interestingly, the A. mellifera genome has been described as structured with respect to plasticity, where stress-related genes are organised into clusters that show coordinated gene expression in response to environmental changes (Duncan et al., 2020). In addition to epigenetics, RNA processing mechanisms are also involved in phenotypic plasticity; splicing is the process by which introns are removed from a gene’s primary transcript, resulting either in a single functional protein (constitutive splicing) or in various structurally and functionally isoforms (alternative splicing). Alternative splicing (AS) can be triggered by physiological needs and environmental stimuli, often representing a primary source of phenotypic diversity within the proteome of eukaryotic cells [reviewed by Ast (2004), Blencowe (2006), Maniatis and Tasic (2002)]. It also plays a key role in cellular stress tolerance [reviewed by Biamonti and Caceres (2009)].

Genetic studies are a valuable tool to understand the effects of stressors on different molecular mechanisms. This review proposes genetic markers for Apis mellifera focusing mainly on climate change stressors such as temperature, humidity, UV exposure and food scarcity. These markers belong to different molecular mechanisms summarized in Figure 1, ranging from genetics to epigenetics, and affecting an organism´s phenotype. Thus, it can serve as a basis for designing new studies that consider the relationship with the stressor(s), as well as the type of gene expression (constitutive or inducible), splicing-related processes (key in PCR designs based on RNA) and epigenetic factors that might affect gene expression. In addition, the knowledge of the role of each marker in the corresponding cellular pathways will help future studies to select accurate molecular markers and correctly interpret the results.

FIGURE 1

Diagram illustrating the process of phenotypic plasticity. It shows how genotype leads to phenotype through genetic transcription, pre-mRNA, constitutive splicing, and mature mRNA, resulting in protein isoforms. Epigenetics and alternative splicing (AS) impact this process, allowing response to environmental cues. Changes in physiology and behavior result from this plasticity. The diagram highlights pathways like heat shock response, immunity, redox metabolism, and metabolic pathways, emphasizing adaptation through advantageous mutations.

Overview of the genetic and cellular processes that lead to phenotypic changes. Cellular biomarkers include genes related to stress resistance (genotypes and gene-expression levels), epigenetic factors influencing gene expression, detection of mRNA variants produced by alternative splicing (AS), and the study of the corresponding synthetized proteins and their role in the different stress-related pathways such as heat shock response, immunity, metabolism and RedOx homeostasis. All of these characteristics confer different phenotypic plasticity to each individual.

2 Methods

This review was based on a published bibliography of original articles, reviews, book chapters and web pages obtained from scientific sources such as PubMed (PubMed, 2025), Scopus (Scopus, 2025), Google Scholar (Google Scholar, 2025), and Web of Science (Web Of Science, 2025). Candidate search terms are shown in Supplementary Material (Supplementary Table S1).

Gene descriptions and nomenclature follow the conventions established by NCBI. Given that the Apis mellifera genome has been characterized as organized in clusters associated with phenotypic plasticity (Duncan et al., 2020), the genetic markers referenced throughout this study were visualized as loci mapped onto their respective chromosomes (Figure 2). This was accomplished using the R package chromoMap v0.4.1 (Anand and Rodriguez Lopez, 2022). Such spatial representation of genomic features may prove instrumental for future research, particularly when selecting genetic markers based on their chromosomal location.

FIGURE 2

Bar graph showing six panels labeled A to F, each with sixteen rows labeled LG1 to LG16. Panels depict colored segments on horizontal bars representing genomic elements across a 25-megabase span, with segments varying in color and position among panels.

Genomic map in Apis mellifera. Distribution of genes related to (A) heat shock response (Supplementary Table S2), (B) other heat stress-related pathways (Supplementary Table S3), (C) humidity (Supplementary Table S4), (D) Ultraviolet exposure (Supplementary Table S5), (E) nutrition (Supplementary Table S6), and (F) epigenetics (Supplementary Table S7). The figures were created using the R package chromoMap v0.4.1 (Anand and Rodriguez Lopez, 2022). The chromosomes were annotated using the honey bee reference genome Amel_HAv3.1.

3 What is stress for honey bees?

Although the term “stress” is commonly understood to have a negative connotation, stress is a natural state pervasive in all known biological systems. These systems maintain a complex dynamic equilibrium called homeostasis that is vulnerable to destabilization in the presence of stressors. After many years of controversy, stress is defined as a “state of homeostasis being challenged,” and stressors as “the factors with the potential to directly challenge homeostasis” [reviewed by Lu et al. (2021)]. An optimal stress level plays a key role in the health and adaptability of organisms. Depending on the level of exposure to the stressor(s), stress can be classified as: (i) distress (or bad stress), which occurs when the state of homeostasis is severely challenged by high levels of stress, inducing a severe response that impairs homeostasis and threatens health; (ii) sustress (or inadequate stress), which is the consequence of inappropriate stressor effects that do not challenge homeostasis but undermine its capacity, and threaten health, and (iii) eustress (or good stress), in which the system is mildly challenged by moderate levels of stressors, inducing a mild response, enhancing the buffering capacity of homeostasis and benefiting health [reviewed by Lu et al. (2021)]. Previous terminology defined “physiological” stress as encompassing environmental stress, intrinsic developmental stress and ageing [reviewed by Kagias et al. (2012)], initially understood as eustress. However, if an organism is unable to cope with one or more of these stressors, the level of stress increases (distress) and becomes pathological.

Regardless of its type or intensity, stress has various consequences at different levels. In honey bees, these consequences can be observed in their behavioral, physiological, and cellular responses. The first consists of the term “fight-or-flight.” For example, sting extension has been used to evaluate sensitivity to stressors in honey bees. This behaviour is widely considered to be indicative of stress as well as an aggressive response. Physiological responses include hormone and neurotransmitter levels. Finally, cellular responses consist of the activation of several molecular mechanisms to control the stressful state and restore the damage. Perhaps the best known and most studied cellular response is the production of stress-related proteins, such as Heat Shock Proteins (HSPs) and antioxidant enzymes, known as stress biomarkers (review by Even et al. (2012)]. At all these levels, cells and systems use different strategies to respond and adapt to environmental changes. Suppose these strategies are effective in coping with the stressful situation. In this case, the organism would undergo a process of adaptation, enabling it to survive and reproduce under new conditions.

4 Temperature

Rising temperature is one of the main effects of climate change (IPCC-Intergovernmental Panel on Climate Change, 2007) and is of concern because temperature is one of the most stressful abiotic stresses on living organisms. But in addition to the effects of temperature on individual organisms, the consequences for an entire biocenosis depend on how each organism changes to adapt. For example, changes in temperature are affecting the hibernation cycle of honey bees (which would affect pollination, and hence plant reproduction), but also the timing and amount of flowering (which changes the availability of food resources for bees) (Obeso and Herrera, 2018). This example is very important because nutritional stress is an important indirect stressor to honey bees from rising temperatures [reviewed by Quinlan and Grozinger (2023)] and is discussed in section 7 of this review.

Focusing only on honey bees, the effect of temperature must be analysed from two perspectives: as a super-organism (colony) and at the individual level (single bee). In both cases, honey bees have developed mechanisms to resist stress and restore homeostasis. As a super-organism, they can regulate the temperature inside the colony, maintaining it constant at approximately 35 °C. In a cold environment, all the bees cluster around the queen and brood to generate and maintain heat (Southwick, 1985; Southwick and Heldmaier, 1987). In a warm environment, the colony employs a combination of ventilation, achieved through coordinated wing fanning by honey bees at hive entrance (at a rate of up to 60 L/min), and evaporative cooling, achieved through the use of water from the bee´s bodies to create droplets that cool down the nest (Southwick and Heldmaier, 1987; Jones and Oldroyd, 2006). All these actions at the community level require also actions at the individual level. Here, cellular mechanisms are key to coping with temperature changes.

Heat stress negatively impacts key biological processes in honey bees, including physiological and behavioral development and immunocompetence (Bordier et al., 2017; Medina et al., 2018; Alqarni, 2020; Greenop et al., 2020). It is worth mentioning that the impact of heat stress on honey bees and their defence mechanisms vary across species and even subspecies [reviewed by Abou-Shaara et al. (2017), Zhao et al. (2021)]. In Apis mellifera, two of the 31 currently recognized subspecies (Ruttner, 1988; De la Rúa et al., 2009), A. mellifera ligustica (native to Italy) and A.m. carnica (native to the Balkan region), have been the focus of several comparative studies due to their widespread commercial use beyond their native ranges. A study on thermal limits and metabolic rates showed that A .m. ligustica foragers are more tolerant to high temperatures than A. m. carnica (Kovac et al., 2007), suggesting that the Italian honey bee is better fit for warmer climates. In Saudi Arabia, several studies have compared the two commercial subspecies with the native A. m. jemenitica, which evolved in semi-arid and desert environments. Interestingly, all the studies showed that A. m. ligustica and A. m. carnica had lower heat tolerance and survival rates than A. m. jemenitica (Abou-Shaara et al., 2012; Alattal et al., 2015; Alqarni et al., 2019). In another study comparing the native subspecies of Algeria, A. m. sahariensis and A. m. intermissa, a differential response to heat stress was also found. Interestingly, A.m. sahariensis, the subspecies adapted to the desert environment with extreme temperatures, reacted better to heat than A.m. intermissa, the subspecies native to the southern shore of the Mediterranean, where the average temperatures are comparatively milder (Khedidji et al., 2024). As for interspecific variability, a study in China showed a higher temperature adaptability in Apis cerana than in A. mellifera. However, the contrary was observed on survival assays; under constant heat and humidity, A. mellifera showed higher survival rates than A. cerana (Li et al., 2019). These findings suggest that variations within the genus Apis must be studied to understand the different thermotolerances and metabolic responses to heat stress, as well as the expected differential transcriptional regulation. However, there are still many honey bee species and subspecies for which no information is available.

4.1 Heat shock response

Heat is a highly negative stressor that alters important cellular structures and mechanisms, including actin filament organisation, protein aggregation, disruption of intracellular transport, fragmentation of the Golgi and endoplasmic reticulum, altering membrane-bound organelles, RNA splicing, ribosomal activity and translation, and ultimately leading to cell cycle arrest and stagnation of growth and proliferation [reviewed by Richter et al. (2010), Morimoto (2011)]. Due to its highly detrimental impact on cell integrity, heat triggers a quick and robust heat shock response (HSR) at the molecular level. This response involves the activation of heat shock genes encoding heat shock proteins (HSPs), antioxidant metabolism genes, and genes associated with alternative pathways related to heat stress (Figure 3). Activation of the HSR enhances cellular resistance to stress, a phenomenon known as thermotolerance [reviewed by Calderwood and Ciocca (2008)]. Thermotolerance is defined as the temperature range between the lowest temperature (critical thermal minimum, CTmin) and the highest temperature (critical thermal maximum, CTmax) at which an organism can maintain muscle control [reviewed by Perez and Aron (2020)]. Consequently, thermal tolerance and its degree of plasticity play key roles in determining the geographical distribution of species.

FIGURE 3

Diagram illustrating cellular stress response mechanisms. The process starts with external stressors causing reactive oxygen species (ROS) production, leading to unfolded proteins. Chaperones assist in refolding proteins. The heat shock response (HSR) involves transcription factors like HSF, influencing heat shock (HS) and antioxidant (AO) gene expression. In the cytoplasm, proteins like HSP90, HSP70, and HSP60, alongside small heat shock proteins (sHSPs), play roles in protein maintenance. The nucleus shows interactions between NF-Y, HS genes, AO genes, ZFP transcription factors, and AO enzymes, highlighting the cellular stress defense system.

General outline of heat stress in Apis mellifera showing related genes. One stressor, like heat, triggers the heat shock response (HSR) when proteins denaturation occurs. The HSR starts with the activation of Heat shock Factor (HSF), a transcription factor that induces the expression of heat shock (HS) genes. These genes encode various chaperones (HSP90, HSP70, HSP60 or chaperonines, HSP40 and small HSPs), which assist in refolding denatured proteins. Specially HSP40 (DnaJ) also regulates the expression of many antioxidant and HS genes. Stressors typically generate ROS, whose levels are controlled by antioxidant enzymes (AO). The genes encoding these enzymes can be activated by the presence of ROS, the action of DnaJ proteins (HSP40), the Nuclear Factor-Y (NF-Y), and Serine Proteases (SPs). Black, purple and yellow stars indicate alternative splicing (AS), involvement in immune pathways, and relationship with antioxidant metabolism, respectively.

However, the HSR is not triggered by heat per se but by the presence of unfolded proteins that result from a variety of stresses (Figure 3), including oxidative stress, heavy metals, ethanol or other toxic substances [reviewed by Richter et al. (2010)]. Therefore, some of the genes related to thermal stress will also be involved in another kind of resistance or adaptation associated with different stressors. On the other hand, the AS mechanism is also important in the HSR (Fujikake et al., 2005; Ruan et al., 2015) due to the generation of new and different isoforms relevant to function under new stressful conditions.

The temperatures experienced by honey bees during their normal activity result in a strong induction of HSR (McKinstry et al., 2017). The HSR usually starts with the activation of Hsp genes through transcriptional factors, the heat shock factor (HSF) family, which specifically bind to the heat shock element (HSE) in the promoter region of the Hsp genes, and regulate their activity in normal conditions [reviewed by Morimoto (2011), Garbuz (2017). Therefore, the HSF family is considered a transcriptional activator of all significant Hsp genes, and it is required in the regulation of various environmental stresses [reviewed by Zhang et al. (2011), Wu (1995), Zou et al. (1998)]. Four Hsf-related genes have been described in A. mellifera (Supplementary Table S2), and changes in their expression have been observed in honey bees suffering from heat stress (McKinstry et al., 2017; Al-Ghzawi et al., 2022).

4.1.1 Heat shock proteins

The key to HSR lies in the heat shock proteins (HSPs), which are proteins that protect or restore cellular structures or components that are prone to damage by heat shock. Not surprisingly, there is a broad functional classification of these stress-inducible proteins, such as molecular chaperones, components of the proteolytic system, RNA- and DNA-modifying enzymes, metabolic enzymes, regulatory proteins (like transcription factor or kinases), proteins involved in sustaining cellular structures such as cytoskeleton, and finally, transport, detoxifying, and membrane-modulation proteins (Richter et al., 2010). Under moderate stress, the synthesis of HSPs is the basis of resistance to such stress and may even provide “cross-protection” against other types of stress (Lindquist, 1986).

The most well-known HSPs are the molecular chaperones, both because they were the first to be discovered and because they are the most studied in relation to thermal stress. Their primary function is to assist in the folding and unfolding of proteins, and they are usually classified according to their molecular weight: HSP100 (78–104 kDa), HSP90 (82–96 kDa), HSP70 (68–78 kDa), HSP60 or chaperonins (60 kDa), HSP40 or DNAj domain proteins (40 kDa) and small HSPs (sHSPs) or α-crystallin proteins (10–30 kDa) [reviewed by Garbuz (2017)]. Another classification of HSPs is based on the pattern of gene expression of the corresponding genes, which may be inducible or constitutive. Inducible Hsp genes are typically expressed at extremely low levels, but under stress, their transcription intensity increases by a factor of several hundred or several thousand relative to baseline. Contrarily, constitutive Hsp genes are expressed at relatively high levels even under normal temperatures (they are supposed to possess a physiological function), and their transcription increases only several-fold under stress (Pardue et al., 1980). In honey bees, HSPs seem to serve a significant role in their responses to abiotic stressors (Table 1), like high and low temperature, dehydration, UV radiation, and starvation (Kim et al., 2019). The most important HSPs involved in HSR in A. mellifera are HSP90, HSP70, HSP60, HSP40, and small HSPs (sHSPs) [reviewed by Abou-Shaara (2024)].

TABLE 1

AbbreviationDescriptionStressorGene expression
TemperatureHumidityUV radiationNutritionImmunityPatternmRNA variants
HSFHeat shock factor proteinsxxId
HSP90Heat shock proteins 90 kDaxxCx
HSP70Heat shock proteins 70 kDaxxI
HSC70Heat shock proteins 70 kDa- cognate formsxxxCx
HSP60Heat shock protein 60 kDaxxId
HSP40Heat shock proteins 40 kDaxxxId
HSP10Heat shock protein 10 kDaxxId
sHSPSmall Heat shock proteinsxxxId
NF-YNuclear factor Y subunitsxxxIdx
ZFPZinc finger proteinsxxxIdx
SPSerine proteasesxxxxId
AcheAcetylcholinesterasexxxId
Crh-BPCorticotropin-releasing hormone-binding proteinxxId
AOAntioxidant enzymesxxxxxIdx
Pla2Phospholipase A2xxId
AfpAntifreeze proteinxxId
PPOProphenoloxidasexxxxId
FOXOForkhead box protein sub-group OxId
ILPInsulin-like peptidesxId
IGFInsulin-like growth factorxId
InRInsulin-like peptide receptorxId
ChicoInsulin receptor substratexId
PI3KPhosphatidylinositol 3-kinasexxId
TorSerine/threonine-protein kinase TorxxId
AKTAKT serine/threonine protein kinasexxId
AdcyAdenylyl cyclasexId
OAOctopamine receptorsxId
VgVitellogeninxxxId
HexHexamerinsxId

Genetic marker groups, their relationship to stress factors and immunity, and the characteristics of their gene expression. Grey letters indicate an indirect relationship through another stressor (humidity influences directly the heat stress degree, and hence, the corresponding molecular response). In the gene expression pattern, “C” and “I” stand for constitutive and inducible, respectively. “Id” indicates ‘indeterminate’ for those markers that are considered inducible or constitutive depending on the context.

Among the many HSP families, HSP90 is one of the most abundant and universally expressed stress proteins. Hsp90 gene is generally overexpressed during HSR (Brunt and Silver, 1991; Pratt and Toft, 2003; Sonoda et al., 2007), and has also been related to morphological evolution (Rutherford and Lindquist, 1998; Sollars et al., 2003; Specchia et al., 2010), reproduction, and brain development (Itoh et al., 1993; Izumoto and Herbert, 1993; Furay et al., 2006; Johnson and Brown, 2009). In A. mellifera, Hsp90 is caste- and age-specifically expressed in adult bees (Aamodt, 2008), suggesting that Hsp90 is expressed constitutively. In the genome of A. mellifera, two homologous cytoplasmic Hsp90 genes are described, Hsp83 and Hsp90, which are located in linkage group 1 (LG1) and LG7, respectively (Supplementary Table S2). These two genes produce two transcripts (A and B) by constitutive splicing (transcribed from the homolog in LG1 and LG7, respectively). Moreover, at least nine transcripts have been described from the LG7 gene by the AS under stress conditions (Xu et al., 2010).

HSP70s are a large family of ubiquitous molecular chaperones that protect cells from the damaging effects of many proteotoxic stresses [reviewed by Hendrick and Hartl (1993), Voos and Röttgers (2002), Richter et al. (2010), Walter and Ron (2011), Rosenzweig et al. (2019)]. Hsp70 genes are inducible, have no or relatively short introns and are preferentially translated, allowing HSP70 proteins to accumulate rapidly in response to adverse environmental stimuli (Gkouvitsas et al., 2009; Sørensen, 2010; Zhang and Denlinger, 2010). There are no relevant data about the genetic variability of Hsp70 in A. mellifera, and three genes are annotated in its genome: Hsp70Ab-like, Hsp70Cb, and the mitochondrial Trap-1 (Supplementary Table S2). Furthermore, the expression of Hsp70 genes is upregulated in response to many stressors (Hendrick and Hartl, 1993), and it seems to be a good indicator of colony stress in different A. mellifera subspecies (Alqarni et al., 2019; Morammazi and Shokrollahi, 2020).

Within the HSP70, there is another group of proteins called “cognate forms” or HSC70. The HSC70 protein family is structurally and functionally similar to HSP70, but its properties are different (Liu et al., 2012). Hsc70 genes are constitutive, and their corresponding proteins are involved in regulating the life cycle of various viruses, as it has been described in “Immunity” section. These genes contain more introns than those of Hsp70, and their number is conserved in vertebrates but it is variable in invertebrates (Chuang et al., 2007). In A. mellifera there are three annotated Hsc genes (Supplementary Table S2) and two of them have been linked to heat stress. HSC70-3 is a conserved endoplasmic reticulum chaperone (Pincus et al., 2010; Walter and Ron, 2011; Johnston et al., 2016; McKinstry et al., 2017) which is induced in honey bees under heat stress (i.e., 45 °C for 4 h) (McKinstry et al., 2017). Although viral infections alone induce Hsc70-3 expression, heat shock alone did not always lead to overexpression. However, honey bees that were both virus-infected and heat-shocked displayed greater Hsc70-3 expression (McMenamin et al., 2020). On the other hand, Hsc70-4 is a core heat shock response gene that is induced by exposing honey bees to heat stress (i.e., 42 °C and 45 °C for 4 h) (Elekonich, 2009; Mahat et al., 2016; Solís et al., 2016; McKinstry et al., 2017). There is some heterogeneity in the expression of this gene, and genetic differences between subspecies may be one of the causes (McMenamin et al., 2020). The antiviral effect of both HSC70-3 and HSC70-4 proteins is due to the relationship between heat shock and RNAi machinery (discussed in the “Immunity” section below).

Chaperones of the HSP60 family are one of the most important components of the protein folding system in the mitochondrial matrix (Martin, 1997). HSP60, or chaperonins, forms a large homo-oligomeric protein complex with an inner cavity that provides a protected environment for the ATP-dependent folding of unfolded or newly synthesized single proteins or protein domains (Voos and Röttgers, 2002). HSP60 works with the cochaperone HSP10 (Voos and Röttgers, 2002), which is supposed to coordinate the behaviour of the single HSP60 monomers and regulate the ATPase cycle (Martin et al., 1993). In A. m. ligustica, highly expressed Hsp60 has been described at 45°C (Alqarni et al., 2019), and Hsp60 and Hsp10 under heat and cold stress (Kim et al., 2019). In the genome of A. mellifera, there is one gene for each HSP60 and HSP10 protein (Supplementary Table S2).

The group of HSP40/J-domain-containing proteins are the largest class of HSP70 cofactors. They bind the nonnative protein and deliver it to HSP70. The J domains of these proteins interact with the ATPase domain of HSP70 and stimulate the hydrolysis of bound ATP (Kampinga and Craig, 2010). HSP40 is categorized into three subfamilies, namely, DnaJA, DnaJB, and DnaJC (Craig et al., 2006; Kampinga and Craig, 2010; Craig and Marszalek, 2017; Zhao et al., 2021). DnaJA1 regulates the expression of many antioxidant genes and heat shock genes (Figure 3), thereby improving the antioxidant ability of bees under heat stress (Li et al., 2020a). It is important to note that the expression levels of DnaJA1, DnaJB12 and DnaJC8 are upregulated under UV radiation, cold, and pesticide treatment in A. cerana cerana, and their silencing attenuates the resistance of this subspecies to λ-cyhalothrin stress (Li et al., 2018a). Within the genome of A. mellifera there are 18 genes of the DnaJ family (Protein lethal (2) essential for life-like) (Supplementary Table S2).

Finally, small HSPs (sHSPs) are a wide family of proteins involved in HSR but also have protective functions under different stresses, such as cold, drought, oxidation, hypertonic stress, UV, heavy metals, and even stress by high population density (Dasgupta et al., 1992; Wang et al., 2007; Waters et al., 2008; Shih et al., 2021). It is important to note that some sHSPs also have a chaperone function in development (Sun and MacRae, 2005). Thus, their corresponding genes appear to be a family of both constitutive and inducible genes. Unfortunately, at present, there is no formal classification of this gene family in A. mellifera according to the pattern of gene expression. Interspecific genetic variability has been described in the genes coding for small HSPs. Insect sHSPs are generally species-specific, suggesting that functions of most sHSPs may have diverged across species. This variability likely reflects the role of sHSPs in the adaptation of insects to diverse ecological niches (Li et al., 2009). The C-terminal sequences of these proteins harbours the conserved α-crystallin domain, while the N-terminal remain variable. This indicates that the conserved C-terminal has a significant part in sustaining the chaperone and other functions, whereas the N-terminal may be associated with diverse expressions, functions, and evolutionary patterns within sHSPs (Li et al., 2009). On the other hand, most sHSP genes that are located on a single chromosome are usually arranged in tandem. This arrangement may enable organisms to rapidly respond to changing environmental conditions due to regulatory advantages. Tandem sHsp genes may be a better way for insects to regulate gene expressions in diverse environments (Li et al., 2009). In the A. mellifera genome there are eight sHsp genes annotated as protein lethal (2) essential for life (Supplementary Table S2).

4.2 Other genes related to heat stress

4.2.1 Nuclear Factor-Y

The Nuclear Factor Y (NF-Y), also known as Heme Activator Protein (HAP) or CCAAT-Binding Factor (CBF), consists of three distinct subunits (NF-YA, NF-YB, and NF-YC) which are found in almost all organisms (Li et al., 2018b; Myers and Holt, 2018). Each animal NF-Y subunit is typically encoded by only one gene, the product of which can undergo different post-translational modifications and have various splicing forms (Li et al., 1992; Mantovani, 1999; Fujikake et al., 2005). In Drosophila melanogaster, NF-Y is essential for the growth and development of the thorax, eye, and R7 photoreceptors, as it regulates multiple signalling pathways, such as the extracellular signal-regulated kinases (ERK) and the c-Jun N-terminal kinases (JNK) pathways (Yamaguchi et al., 2017; Li et al., 2018b). In bees, the expression levels of NF-YA, NF-YB, and NF-YC are induced by long- and short-term heat stress in A. c. cerana and A. mellifera (Zhao et al., 2021). In A. c. cerana, the knockdown of NF-YB and NF-YC decreases the antioxidant capacity and increases the oxidative damage caused by heat (Figure 3). Upregulation of NF-Y may increase the heat resistance of bees under different heat stress conditions by reducing oxidative damage and enhancing antioxidant ability (Li et al., 2020b). The expression of NF-YA, NF-YB, and NF-YC in A. mellifera is not only induced by heat but also by several stress conditions (Table 1), including cold and UV light (Li et al., 2020b). There are four NF-Y genes annotated in the genome of A. mellifera (Supplementary Table S3).

4.2.2 Zinc finger proteins

Zinc finger proteins (ZFP) are among the most abundant proteins in eukaryotic genomes. Their functions are diverse and include DNA recognition, RNA packaging, transcriptional activation, regulation of apoptosis, protein folding and assembly, and lipid binding (Laity et al., 2001). The expression of some genes encoding ZFP has been linked to heat stress (Droll et al., 2013; Liu et al., 2015) (Figure 3), cold, pesticides, and UV exposure in A. cerana, suggesting an important role of these proteins in resistance to a variety of environmental stressors (Guo et al., 2021) (Table 1). In A. mellifera, five ZFP genes (ZFP271, ZFP37, ZFP239, ZFP776, and ZFP93) have been described as upregulated under high-temperature exposure, while two (ZFP62 and ZFP628) have been described as downregulated (Ma et al., 2019). There are several ZFP genes annotated along the genome of A. mellifera, so here we summarize the seven already cited (Supplementary Table S3). There are four loci annotated for ZFP271, two for ZFP37, two for ZFP239, and only one for zinc finger protein 776-like. The ZFP93 gene does not appear in the genome assembly of A. mellifera, but it encodes a ZFP of the KRAB (Kruppel-associated box) subfamily (Bellefroid et al., 1993), and there is a kruppel homolog 1 (Kr-h1) gene in A. mellifera. Finally, the two genes described downregulated are ZFP62 and ZFP628-like. Regarding AS, different isoforms from ZFP genes have been predicted and annotated in GenBank.

4.2.3 Serine proteases (SPs)

Serine proteases (SP) are endo- and exopeptidases involved in insect immunity and antioxidant systems (Ashida and Brey, 1997), as well as the heat stress response in different bee species (Table 1). For example, in A. cerana, the gene enconding for the Clip-domain serine protease1 (AccSp1) is upregulated by temperature (4, 24 °C and 44 °C), H2O2, heavy metals, UV-light, and pesticides, thus linking SPs to its defence against abiotic stresses (Gao et al., 2019). There is a Sp1 gene annotated in the genome of A. mellifera (Supplementary Table S3). Regarding biotic stressors, SPs are directly involved in immunity, participating in the prophenoloxidase (PPO) activation pathway, RNA interference, and SP proteolytic cascade in the Toll signalling, as described below in the “Immunity” section. Compared with D. melanogaster and Anopheles gambiae, A. mellifera has much smaller gene families of SP, SPH, serpin (Serine Protease Inhibitors), PPO and other immune proteins (Evans et al., 2006). A search of the A. mellifera genome yielded 57 sequences with significant similarity to the S1 protease family: 44 SP and 13 SPH genes (Zou et al., 2006). SPHs are similar in sequence to S1 proteases but lack one or more of the catalytic residues in SPs. In addition, seven annotated genes in the honey bee genome encode five serpins (serpin 1–5) and two serpin-like proteins (Zou et al., 2006). SP inhibitors of the serpin superfamily are present in insect haemolymph to remove excess proteases and maintain homeostasis (Kanost, 1999). Finally, genes of SP putative substrates prophenoloxidase (PPO) and spätzle are described in the A. mellifera genome (Zou et al., 2006) (Supplementary Table S3).

4.2.4 Antioxidant (AO) enzymes

Alterations in the oxidation states of intracellular metabolites and enzymes have historically been considered negative stressors, requiring strictly defensive responses. Cellular growth and survival require the coupling of electron-transfer reactions to the generation of ATP. These reactions depend on key cellular electron carriers and the stability of protein residues and cofactors. Redox enzymes are notoriously nonspecific, transferring electrons to any suitable acceptor they encounter. These molecular mechanisms represent a constant cellular stress balanced and maintained by redox homeostasis (Sporer et al., 2017). There are both primary and secondary antioxidant enzymes, which act directly or indirectly on reactive oxygen species (ROS) molecules. The first line of defence against ROS attack is provided by three different kinds of primary antioxidant enzymes, which act directly on ROS: superoxide dismutase, catalase and peroxidases. Superoxide dismutase (SOD) rearranges superoxide to oxygen and hydrogen peroxide, catalase prevents free hydroxyl radical formation by breaking down hydrogen peroxide into oxygen and water, and peroxidases (POD) catalyse an analogous reaction in which hydrogen peroxide is reduced to water by a reductant that acts as an electron donor, normally reduced thioredoxin (TRX) or glutathione (GSH). In addition, insects have three families of genes that encode antioxidant enzymes that act as peroxidases: TPXs, also known as peroxiredoxins (Radyuk et al., 2001), phospholipid-hydroperoxide GPX homologs with thioredoxin peroxidase activity (GTPX) (Missirlis et al., 2003), and glutathione S-transferases (GSTs) (Tang and Tu, 1994; Toba and Aigaki, 2000). Secondary antioxidant enzymes that act indirectly on ROS include thioredoxin (TrxR) and methionine sulphoxide reductases (MsrA and MsrB), which are involved in protein reparation by catalysing the TRX-dependent reduction of methionine sulphoxide to methionine (Moskovitz et al., 1996; Kumar et al., 2002).

In A. mellifera genome, 38 antioxidant genes were identified. In general, antioxidant genes encode small proteins with less than 250 amino acids, and most of them possess at least one intron (Corona and Robinson, 2006). Alternative splicing has been described as a common mechanism in the RNA processing of these genes, suggesting that the resulting different isoforms may play a role in stress resistance. About genetic variability, 52 and 29 alleles are described for SOD1 and SOD2 in D. melanogaster, as well as 34 alleles for catalase (FlyBase, 2025). On the other hand, there is a relationship between redox enzymes and HSR (Figure 3). High temperatures exponentially increase metabolic rates, which means higher oxygen consumption and therefore higher production of ROS and more oxidative damage to organisms (Finkel and Holbrook, 2000; Gillooly et al., 2001; Belhadj Slimen et al., 2014; González-Tokman et al., 2020). Under long-term stress, the expression of additional Hsp and other gene family members is upregulated, which in turn scavenges ROS, enhances the antioxidant defence system of bees, and increases their survival rate [reviewed by Zhao et al. (2021)]. As mentioned in the “Heat Shock Proteins” subsection, DnaJA1 (HSP40) can regulate the expression of many antioxidant genes, thereby improving the antioxidant ability of bees under heat stress (Li et al., 2020b). It is important to note that oxidative stress response mechanisms are activated by many types of factors, including not only heat stress but also nutritional stress, pesticides, pathogens, UV radiation, among others (Table 1). On the other hand, immunity molecular pathways can influence ROS levels by increasing the expression of antioxidant proteins (Zhang et al., 2016), as described in the “Immunity” section. Genes encoding antioxidative enzymes in A. mellifera are noted in Supplementary Material (Supplementary Table S3).

4.2.5 Acetylcholinesterase

Acetylcholinesterase (AChE) is a serine hydrolase that controls synaptic and neurohumoral cholinergic activity by hydrolyzing the neurotransmitter acetylcholine into acetic acid and choline [reviewed by Silman and Sussman (2008)]. A. mellifera has two genes, Ace1 and Ace2, that encode AChE1 and AChE2 proteins, respectively (Supplementary Table S3). In bee species belonging to Bombus and Apis, AChE2 acts as the synaptic enzyme, while AChE1 shows little catalytic activity. This suggests that AChE1 may have become specialized to play non-synaptic functions (Kim and Lee, 2013). In A. mellifera, the expression of Ace1 changes depending on the rearing stage of the colony (Kim et al., 2017), and it appears to be connected to HSR (Kim et al., 2019). In addition, the genetic variability of Ace1 linked to pesticide resistance has been described in several insect species (Nabeshima et al., 2003; 2004; Weill et al., 2004; Baek et al., 2005; Oh et al., 2006; Alon et al., 2008; Kozaki et al., 2008; Jiang et al., 2009; Ramphul et al., 2009; Wu et al., 2010). AChE is the main target of organophosphorus (OP) and carbamate (CB) insecticides, to which A. mellifera displays unique sensitivity profiles (Hardstone and Scott, 2010), perhaps because AChE2 is its main synaptic AChE (Kim and Lee, 2013).

4.2.6 Corticotropin-releasing hormone-binding protein (CRH-BP)

Throughout evolution, highly conserved signalling molecules have been utilized to integrate stress responses, highlighting their important roles in survival. The corticotropin-releasing hormone (CRH), also known as corticotropin-releasing factor (CRF), is a neuroendocrine peptide that regulates various physiological responses to stresses (Seasholtz et al., 2002). It is considered to act as a neurotransmitter, coordinating various autonomic, hormonal, and behavioural responses to stresses, and may be involved in developmental processes (Muglia et al., 1995; Cortright et al., 1997; Majzoub, 2006). The high homology of the CRH binding protein (CRH-BP) in honey bees to that in humans suggests that its function(s) has also been evolutionarily conserved (Ketchesin et al., 2017). In A. mellifera, CRH-BP is a 322-amino-acid soluble protein structurally unrelated to the CRH receptors. The CRH-BP gene (Supplementary Tables S3 and S5) is well-conserved and identifiable in insects. Although very few studies have focused on CRH-BP and stress in invertebrates, the Chinese honey bee A. cerana cerana subjected to UV light, heat, or cold exhibited increased CRH-BP mRNA in the head in a time-dependent manner (Liu et al., 2011). These studies highlight the upregulation of invertebrate CRH-BP mRNA in response to stress (Ketchesin et al., 2017). Therefore, CHR-BP seems to be a promising protein candidate as a potential element involved in the stress response, and its action would be worth investigating in future studies (Even et al., 2012).

5 Relative humidity

While temperature has been the main focus of climate change studies, changes in precipitation patterns are also key, not only because of water availability but also because of fluctuations in relative humidity (RH) (IPCC-Intergovernmental Panel on Climate Change, 2007). In A. mellifera colonies, a suitable RH of up to 75% is required for egg hatching (Ellis et al., 2008), and changes in RH can significantly affect larval development during incubation (Human et al., 2006) [reviewed by Abou-Shaara et al. (2017)]. In fact, when RH falls below 50%, there is a significant reduction in the number of normally hatched eggs (Li et al., 2016). Comparative studies of A. mellifera subspecies revealed that A.m. jemenitica eggs exhibit higher hatching rates than those of A.m. carnica at 50%–75% RH, while no eggs from either subspecies hatched at 30% RH (Al-Ghamdi et al., 2014).

The bees can regulate humidity within the colony by evaporating water from nectar and regurgitating liquid droplets, helping to restore favourable RH conditions while also aiding in thermoregulation (Heinrich, 1980; Human et al., 2006; Kovac et al., 2007). As for external conditions, no clear direct impact of RH in honey bees, including foraging activity, has been reported (Joshi and Joshi, 2010).

High humidity alone appears to have a limited impact on bees compared to high temperatures (Ma et al., 2019). However, rather than acting as an isolated stressor, humidity most relevant effect likely occurs in combination with temperature. These two variables are closely interconnected and, when combined, act as a new independent stressor. For example, low RH levels combined with high temperature can exacerbate heat stress, whereas high RH can reduce its severity. The highest rate of body water loss in A. mellifera occurs at 35 °C/0% RH, and the lowest occurs at 25 °C/75% RH and 30 °C/100% (Atmowidjojo et al., 1997). Studies on A. m. jemenitica and A. m. carnica showed that body water loss rates increased with higher temperature and lower RH, while higher RH generally improved workers’ survival (Abou-Shaara et al., 2012). In A. mellifera and A. cerana, high RH was shown to have a protective effect on bee survival under high temperatures, which may be due to the inhibition of body water loss (Li et al., 2019). Therefore, although the effect of RH on bees was previously considered negligible (Joshi and Joshi, 2010), it appears to become an important stress factor when it drops below 50% in warm environments, conditions increasingly common in arid regions under a climate change scenario.

The protective effect of high RH under high temperatures is also evident at the molecular level. Ma et al. (2019) demonstrated this in their study on A. mellifera when examining its responses to various combinations of temperatures and RH. The authors reported 434 differentially expressed genes (DEGs) under high-temperature treatments, 86 under high-humidity treatments, and 266 in combined high-temperature and humidity treatments. This suggests that high humidity reduces the expression of heat response genes by nearly half in a heat-stress environment. However, when focusing solely on molecular responses to humidity stress under optimal temperature conditions, the expression of genes such as Pla2 (phospholipase A2) and Afp (antifreeze protein) increases with increasing humidity (Table 1), starting at RH levels above 50% (Ma et al., 2019). Pla2 encodes a cell membrane enzyme (PLA2) that cleaves fatty acids. When its activity increases, the hydrolysis reaction is enhanced, producing a variety of fatty acids, disturbing cell membrane metabolism and nerve signal transmission (Mahalka et al., 2011). It is important to note that the PLA2 protein is involved in immunity through the arachidonic acid release in the eicosanoid pathway (McMenamin et al., 2018). There are 10 genes in the genome of A. mellifera referred to Pla2 (Supplementary Table S4). Afp encodes AFP, a proteinaceous compound with enhanced antifreeze properties, enabling it to bind to small ice crystals and inhibiting their growth and recrystallization. This mechanism helps minimize the damage caused by frozen water/ice to living organisms (Jia and Davies, 2002). There is one gene in the genome of A. mellifera referred to Afp (Supplementary Table S4).

6 Ultraviolet exposure

Climate change has been linked to stratospheric ozone depletion, which increases the amount of ultraviolet-B (UV-B) radiation that reaches the Earth’s surface (IPCC-Intergovernmental Panel on Climate Change, 2007). This affects how organisms and ecosystems respond to it (Bernhard et al., 2020). UV radiation is well-known for its detrimental effect on organisms. Yet, it is the universal source of non-ionizing radiation and is essential for life and its development on Earth. According to CIE (Commission Internationale de l’Eclairage, International Commission on Illumination), UV radiation can be divided into three ranges: UVA (320–400 nm), UVB (280–320 nm), and UVC (200–280 nm). Earth’s ozone layer blocks the majority of UVC, and a significant portion of UV-B, so the light near the Earth’s surface is enriched by UV-A, which has also been linked to oxidative types of mutation due to oxidative stress following irradiation [reviewed by Cockell and Knowland (1999)].

UV promotes photochemical reactions of ROS formation. In fact, one of the most susceptible biological targets of UV radiation is the DNA. When UV radiation enters a cell, it is absorbed by the aromatic rings of nucleotides and amino acids, leading to DNA and protein damage, respectively [reviewed by Cockell and Knowland (1999)]. However, it is important to note that repair mechanisms are not 100% efficient, so the lower the exposure to UV-radiation, the lower the damage and the greater the benefit to the organism [reviewed by Cockell and Knowland (1999)]. Hence, UV protection is key, particularly in a scenario of climate change with a prospective increase of UV radiation exposure.

Cells use a wide array of biological macromolecules to protect themselves against ROS by quenching oxygen free radicals. In insects, melanins are good examples of such macromolecules. Melanins can be classified into two groups: brown to black pigments termed eumelanin, and alkali-soluble yellow to reddish-brown pigments termed pheomelanin (Ito and Wakamatsu, 2008). They enhance the protective properties of the cuticle, acting as important barriers against environmental stressors such as UV radiation (Charlier et al., 2020). Melanins are produced by the enzymatic oxidation of tyrosine by tyrosinase followed by the conversion of dopa to 5,6-dihydroxy-indole. This is a phenolic and indolic compound and the basic building block of the eumelanin polymeric structure. This structure acts as the UV-absorbing chromophore [reviewed by Cockell and Knowland (1999)]. Recently, Charlier et al. (2020) used the electron paramagnetic resonance imaging (EPRI) technique to detect melanin in A. mellifera. They identified free radicals almost exclusively in the cuticle of the bee periphery of A. mellifera consistent with a eumelanin signal. This finding suggests that melanin–chitin complexes in the honey bee cuticle play a key role in UV defence (Charlier et al., 2020). Interestingly, these authors detected the presence of other free radicals in the centre of the honey bee head, suggesting the possible presence of neuromelanin in its brain, as in Drosophila (Barek et al., 2018).

At the molecular level, melanization is an immunological process that results from the combination of humoral and cellular processes that occur during encapsulation and healing to cope with non-pathogen-mediated and pathogen-mediated injuries. In the honey bee, and insects in general, melanization acts as an important cellular defence mechanism responsible for eliminating a large number of bacterial cells, parasites and xenobiotics (Eleftherianos et al., 2009; McMenamin et al., 2018). Simultaneously with the formation of melanin and its polymerization together with other proteins to encapsulate the invading agent, reactive oxygen and nitrogen intermediates are released, which collaborate with the destruction of the agent (Negri et al., 2016). Melanization is mediated by a humoral protein, prophenoloxidase (PPO) (González-Santoyo and Córdoba-Aguilar, 2012; Negri et al., 2019) and A. mellifera possesses only one PPO gene (Supplementary Table S5).

Several markers involved in temperature stress (see “Temperature” section) are also involved in UV stress (Table 1). For example, in A. cerana, the expression of DnaJA1, DnaJB12 and DnaJC8 are upregulated under exposure to both climatic stressors (Li et al., 2018a). Furthermore, as for temperature, sHSPs can develop the protection function under UV stress conditions not only in honey bees but also in other animals and even plants (Dasgupta et al., 1992; Waters et al., 2008). Although very few studies have focused on CRH-BP response to stress in invertebrates, it has been shown that A. cerana cerana subjected to UV radiation exhibited increased CRH-BP mRNA in the head in a time-dependent manner (Liu et al., 2011). On the other hand, it is worth noting that UV promotes photochemical reactions leading to ROS formation, so antioxidant enzymes are also involved in UV stress. Finally, the serine protease AccSp1 gene, which is directly involved in ROS metabolism, seems to play different roles in resistance to UV radiation in A. cerana (Gao et al., 2019).

7 Nutrition

Although the interactions between plants and their insect pollinators are often the result of a long history of co-evolution, climate change can quickly disrupt the timing of their life cycles by altering their phenology and distribution (Hughes, 2000; Scaven and Rafferty, 2013; Gérard et al., 2020). Changes in plant physiology under realistic climate change scenarios may alter flowering patterns and the duration and intensity of blooming (Scaven and Rafferty, 2013; Inouye, 2020). This disruption could lead to food shortages for pollinators and, in some cases, even contribute to their extinction (Memmott et al., 2007; Dalton et al., 2023).

The abundance and diversity of floral resources play a crucial role in honey bee health, as they significantly influence immune response (Alaux et al., 2010; Di Pasquale et al., 2013; Martelli et al., 2022; Corona et al., 2023). The survival of bee colonies relies heavily on the availability of pollen, which provides proteins, lipids, and micronutrients, and nectar, which supplies carbohydrates (Brodschneider and Crailsheim, 2010). In temperate climates, bees typically experience nutritional stress as winter approaches due to the diminishing availability of food resources (Mattila and Otis, 2007; Knoll et al., 2024). This issue has gained particular attention because many bees face low dietary diversity due to climate change (Goulson et al., 2015).

At the molecular level, bees that consume a protein-rich diet exhibit signs of anabolism, lipid metabolism, and an increased expression of genes encoding nursing-related proteins, such as vitellogenin and major royal jelly protein 1 (Bitondi and Simoes, 1996; Corona et al., 2007; 2023; Ament et al., 2008; Azzouz-Olden et al., 2018). In contrast, bees subjected to suboptimal nutrition, typically characterized by the intake of diets lacking proteins, display a reduction in nitrogen and lipid metabolism, reduction in antioxidant response, and different gene expression patterns (Martelli et al., 2022). These include altered expression in cuticle maturation genes, the over-expression of genes involved in the modulation of circadian rhythm and other genes related to the foraging behavioral state, such as those involved in the insulin/insulin-like growth signaling (IIS), juvenile hormone (JH) and Target of Rapamycin (TOR) (Figure 4) (Supplementary Table S6) (Corona et al., 2007; 2023; Ament et al., 2008; Azzouz-Olden et al., 2018; Martelli et al., 2022). These nutrient-sensing pathways include the Forkhead box O subfamily of gene (FOXO) [reviewed by Murtaza et al. (2017)], a transcription factor involved in multiple biological processes including the regulation of aging, nutrient levels response, and stress response (Hwangbo et al., 2004; Semaniuk et al., 2021a; Semaniuk et al., 2021b). FOXO proteins are known to regulate translation of environment-induced stimuli into gene expression and their antioxidant role [reviewed by Murtaza et al. (2017)]. In A. mellifera genome there is only one FOXO gene (Supplementary Table S6).

FIGURE 4

Diagram illustrating signaling pathways involved in cell processes. The pathways include IIS, Autophagy, and TOR. IIS shows interaction between InR, Chico, and Pi3K, leading to AKT/PKB activation and FOXO inhibition, promoting cell proliferation. Autophagy is regulated through Pi3K, AKT, and RAPTOR/TOR interaction. TOR pathway involves growth factors, Pi3K, AKT, GTP-Rheb, and GDP-Rheb, leading to TOR complex formation and control of cell growth, survival, and proliferation by inhibiting autophagy. Amino acids stimulate TOR activity, shown with blue circles.

Nutritional molecular pathways connected to immunity. In the insulin/insulin-like signalling (IIS) pathway, insulin-like peptides (ILPs) and insulin-like growth factors (IGFs) bind to the insulin receptor and CHICO receptor substrate. Phosfatidylinositol-3-OH kinase (Pi3K) activates the phosphatidylinositol (PiP2 to PiP3) reaction ultimately promoting cell differentiation, and organ growth and thereby an increase in lifespan (Pan and Finkel, 2017). The AGC family kinase AKT inhibits the Forkhead box O subfamily of gene (FOXO). The FOXO proteins play an antioxidant role and regulate the autophagy process. TOR pathway is stimulated by growth factors and free amino acids but inhibited by hypoxia and ATP depletion (Abraham, 2009). TOR is a serine/threonine protein kinase belonging to the phosphoinositide 3-kinase-related family, which is highly conserved among eukaryotes (Pan and Finkel, 2017). Pi3K plays a role in the TOR pathway when it is stimulated by growth factors, activating AKT. Yellow, purple and green stars indicate relationship with antioxidants, immune pathways, and nutrition, respectively (Abraham, 2009; Emlen et al., 2012; Koyama et al., 2013; McMenamin et al., 2018).

The connection between honey bee nutrition and behavioral development is further highlighted by the fact that nutritional stress is associated to precocious foraging (Fewell and Winston, 1992; Schulz et al., 1998), changes in colony demographics, and, ultimately, colony collapse (Perry et al., 2015). Research indicates that similarly to nutritional stress, thermal stress reduces the expression levels of nurse-associated genes (vg and mrjp1) and increases the expression of foraging-associated genes (e.g., juvenile hormone esterase) (Bordier et al., 2017; Corona et al., 2023) suggesting that thermal stress can induce the physiological changes linked to precocious foraging, potentially affecting the overall fitness of the colony. In support of this idea, it has been shown that bees raised under high temperatures exhibited an increased probability of dancing and foraging earlier in life (Becher et al., 2009).

7.1 Phenolamines

Unlike vertebrates, honey bees do not synthesize the catecholamines norepinephrine and epinephrine but use instead the phenolamines tyramine (TA) and octopamine (OA) to perform similar physiological functions (Roeder, 2020). Both biogenic amines affect the locomotor behaviour of adult worker honey bees (Fussnecker et al., 2006). Additionally, OA plays an important role in associative learning and memory (Menzel et al., 1996; Blenau and Baumann, 2001; Roeder et al., 2003) by mediating and modulating the reward in appetitive learning (Hammer and Menzel, 1995; Scheiner et al., 2006; Kim et al., 2013). The activation of receptors of tyramine and octopamine are closely related to adenylyl cyclase activity, which can be activated or inhibited in order of the receptor type or subtype (Blenau et al., 2000; Grohmann et al., 2003; Mustard et al., 2005; Beggs et al., 2011; Balfanz et al., 2014; Reim et al., 2017). Adenylyl cyclase (or adenylate cyclase) catalyzes the conversion of ATP to 3′,5′-cyclic AMP. Genes encoding adenylyl cyclase (Ac76E) and the octopamine receptor (OA) have been described as “starvation genes” in underfed honey bees (Azzouz-Olden et al., 2018). Six adenylyl cyclase genes and four octopamine receptor genes have been annotated in the A. mellifera genome (Supplementary Table S6).

7.2 Storage proteins

Insect storage proteins are an important source of amino acids, particularly during metamorphosis. In honey bees, the most known storage protein is the lipoprotein Vitellogenin (Vg), which possesses multiple functions such as royal jelly production (Amdam et al., 2003), promotion of longevity (Amdam and Omholt, 2002; Seehuus et al., 2006; Corona et al., 2007), and immunity (Amdam et al., 2004). Vg has also been proposed as a plausible candidate for a stress marker together with the juvenile hormone (JH) which is also considered to be related to stress responses. Both are involved in oxidative stress (Seehuus et al., 2006; Corona et al., 2007) and heat stress in A. mellifera (Bordier et al., 2017). There is only one Vg gene in A. mellifera genome (Supplementary Table S6), but three Vg-like genes have also been described with different evolutionary patterns and functions (Salmela et al., 2016).

Other important storage proteins in honey bees are the hexamerins. They are the most abundant proteins in larval haemolymph and essentially participate in the dynamics of amino acid storage and exploitation, which occurs during insect development (Telfer et al., 1991). There are four hexamerin genes in the honey bee genome corresponding to subunits 70A, 70B, 70C and 110 (Supplementary Table S6). The hex70a, hex70b and hex70c genes are arranged in tandem in chromosome 8, whereas hex110 is dispersed in chromosome 11, exhibiting unusual features throughout its sequence (Martins et al., 2010). The JH exerts a strong and positive influence on the expression of hex70b and hex70c, while its effect on the expression of hex70a and hex110 is comparatively weaker (Martins et al., 2010). As for nutritional effects, similar to Vg, hexamerins expression increases when honey bees are fed pollen (DeGrandi-Hoffman et al., 2021), thus being reliable indicators of good nutritional status.

8 Other cellular mechanisms involved in stress resistance

In addition to the pathways that are directly or indirectly involved in stress management, cells have several general mechanisms that are activated in stressful situations, such as changes in RNA processing and epigenetic factors. They all act independently of the organism’s genetic background and directly influence its phenotypic plasticity (Figure 1), and, therefore, its resistance and adaptability to new environments. This section reviews the major ones in A. mellifera and their role in stress.

8.1 Alternative splicing

Modulation of RNA processing is involved in stress tolerance in insects (Fujikake et al., 2005; Ding et al., 2014; Ruan et al., 2015; Pai et al., 2017). In RNA processing, the splicing machinery (spliceosome) recognises exons with high accuracy, removes the introns from the pre-mRNA molecule, and ligates the exons to form a mature mRNA. This process is known as “constitutive splicing.” On the other hand, alternative splicing (AS) is the process by which the exons of primary transcripts (pre-mRNAs) can be spliced in different arrangements to produce structurally and functionally distinct mRNA and protein variants, or isoforms (Figure 5). AS is involved in many physiological processes, including the response to biotic and abiotic stresses [reviewed by Biamonti and Caceres (2009)]. The production of proteins with diverse domain rearrangements from the same gene is the main AS mechanism for pathogen-resistance genes. These mRNA variants have been identified for many genes, particularly those involved in the regulation of stress responses, such as protein kinases, transcription factors, splicing regulators and pathogen-resistance genes (Mastrangelo et al., 2012). This occurs according to physiological needs and environmental stimuli, and often represents a primary source of phenotypic diversity within the proteome of eukaryotic cells (Ast, 2004; Blencowe, 2006; Neves-da-Rocha et al., 2019). Therefore, alternative-splicing isoforms produced by stress-related genes could be used as proxies of stress.

FIGURE 5

Diagram showing RNA splicing processes. On the left, constitutive splicing of pre-mRNA with exons and introns leads to a protein. On the right, stress-induced alternative splicing results in different isoforms. Variants include exon skipping, alternative splice sites, intron retention, and mutually exclusive exons.

Scheme of constitutive and alternative splicing. In constitutive splicing, total coding DNA is retained in the final mRNA product and translated into a protein, and all introns are fully spliced during posttranscriptional RNA processing. In alternative splicing, introns can be spliced in more than one way, resulting in different sets of exons in the mature mRNA, and hence, in different series of related proteins called isoforms (Ast, 2004).

8.2 Epigenetics

Epigenetic phenomena include all processes by which the expression of a gene can be altered (overexpressed or silenced), resulting in a phenotypic change while the genotype remains unchanged. The best-known epigenetic phenomena are the methylation of nucleotides (i.e., fixation of a methyl radical on a nucleotide) and changes in the configuration of histones. Both mechanisms are the major regulators of gene expression in all organisms (Glastad et al., 2019) and act as a major source of phenotypic plasticity (Berger et al., 2009). The genome of A. mellifera has been described as structured with respect to plasticity, with genes related to stress response being organized into clusters that are marked by histone modifications (Duncan et al., 2020). Epigenetic mechanisms can be triggered by environmental factors such as heavy metals or persistent organic pollutants, which can modulate epigenetic marks such as acetylation or methylation (Foret et al., 2009; Collotta et al., 2013). With regard to climate change, DNA methylation and histone/chromatin modifications have been linked to thermal stress responses and facilitate transgenerational epigenetic inheritance of thermal adaptation [reviewed by McCaw et al. (2020), de Carvalho (2023)]. These modifications also enable populations to adapt to local and global climate gradients.

8.2.1 DNA methylation

DNA methylation is a covalent modification that occurs by the addition of a methyl group to the fifth carbon of cytosines, mostly in CG dinucleotides (CpG) (Klose and Bird, 2006; Suzuki and Bird, 2008; Zemach et al., 2010; Moore et al., 2013) (Figure 6), although adenine methylation can also occur (Ratel et al., 2006). Moreover, DNA is not the only nucleic acid that can be modified; RNA can also undergo modifications that impact gene expression post-transcriptionally (Sieber et al., 2021), creating a new layer of dynamic gene regulation [reviewed by Kan et al. (2022)]. Recent studies suggest that the epitranscriptome of honeybees may play a role in stress responses (Bataglia et al., 2021; 2022).

FIGURE 6

Diagram illustrating DNA methylation and histone modification. Methylation of CpG sites by DNMT is shown affecting gene activation, with DNA wrapped around blue histones. Labels indicate "Inactive gene" and "Active gene" based on methylation status. Epigenetic factors influence histone tail modifications. A box details the acetylation of lysine residues, highlighting the roles of HAT and HDAT in histone modification.

Main epigenetic phenomena in honey bees. DNA methylation occurs by the addition of a methyl group to the fifth carbon of cytosines, mostly in CG dinucleotides (CpG). Histone modifications involve the disruption of histone-DNA interactions, causing nucleosomes to unwind and genes to become accessible to the transcriptional machinery, allowing subsequent gene activation.

In contrast to the heavily methylated genomes of mammals, the invertebrate genomes are sparsely methylated in a “mosaic” fashion, with most methylated CpG dinucleotides found across gene bodies (Wedd and Maleszka, 2016). The gene body methylation is frequently associated with active transcription (Wedd et al., 2022), and it seems to be an important feature of caste determination in social insects (Sieber et al., 2021). Some studies on the honey bee have shown strong links between gene body methylation and AS (Wedd et al., 2016). In addition, the existence of non-CpG methylation events in honey bee introns, potentially playing a role in the regulation of AS, has been also proposed (Cingolani et al., 2013).

DNA methylation is controlled by DNA methyltransferases (DNMTs). DNMT1 is responsible for maintaining methylation states across cell divisions, whereas DNMT3 is involved in de novo methylation (Klose and Bird, 2006), although these functions can overlap (Jeltsch and Jurkowska, 2014; Karemaker and Baubec, 2020). DNA methylation can be reversible, in particular through the action of the ten/eleven translocation (TET) family enzymes (Kohli and Zhang, 2013). Regarding maintained methylation states, allele-specific methylation is known to occur in honey bees (Remnant et al., 2016; Wedd et al., 2016). Interestingly, honey bee males seem to harbour individual-specific DNA methylation patterns in their semen and these patterns are often associated with genotypic variation. It means that genes that are variable at the epigenetic level appear to be more likely to be variable at the genetic level. This sequence polymorphism could be an important determinant of the DNA methylation state at many loci in honey bees, contributing both to the individual-specificity of epigenetic marks and to their retention across generations (Yagound et al., 2019). While there is evidence that gene expression patterns are sometimes heritable, additional experimental analyses are required to conclusively demonstrate that DNA methylation is the epigenetic cause of such heritable effect (Glastad et al., 2019).

8.2.2 Histone post-translational modifications (hPTMs)

Eukaryotic DNA is packaged in basic and repeating structural units (nucleosomes), where a segment of DNA wound around the histone cores (Figure 6). These histone cores are composed of several subunits (H2A, H2B, H3 and H4), and each one contains amino acid tails that are sites of post-translational regulation (Spotswood and Turner, 2002). Some modifications disrupt histone-DNA interactions, causing nucleosomes to unwind. In this open chromatin conformation, DNA becomes accessible to the transcriptional machinery, enabling subsequent gene activation. In contrast, modifications that reinforce histone-DNA interactions create a very compact chromatin where the transcriptional machinery cannot access the DNA, resulting in gene silencing.

Histone post-translational modifications (hPTMs) consist of a diverse set of epigenetic signals that can alter transcription either by the addition of a chemical group to a histone protein or by specific protein binding to histone tails (Glastad et al., 2019). Regarding the former mechanism, lysine acetylation is perhaps the most studied modification, as it was one of the first discovered to influence transcriptional regulation. Acetylation of lysine residues results in neutralization of histone charge, weakening the nucleosome structure and making DNA accessible for transcriptional factors binding, significantly increasing gene expression (Roth et al., 2001). Acetyl groups are added to lysine residues of histones H3 and H4 by histone acetyltransferases (HAT) and removed by deacetylases (HDAC). The action of these two types of enzymes results in opposite gene expression outcomes (Bernstein et al., 2007). Diverse compounds can inhibit deacetylases, and they are known as HDAC inhibitors (HDACis). They act by triggering histone tail acetylation and play an important role in epigenetic gene regulation (Marks et al., 2003). In insects, HDACis can accelerate growth, extend longevity, and help overcome injuries (Zhao et al., 2005; Mukherjee et al., 2012; Zhao et al., 2005; Mukherjee et al., 2012) and, when produced in high doses, they may arrest cell growth and induce apoptosis (Shao et al., 2004; Tabuchi et al., 2006). In A. mellifera, the HDACis activity has been linked to the regulation of immune and detoxification genes under stress from pesticides and Nosema (Hu et al., 2017), and, interestingly, to the epigenetic mechanism of royal jelly (Spannhoff et al., 2011). Importantly, similar to methylation patterns, emerging evidence suggests that information stored in nucleosomal hPTMs can be transmitted across cell divisions (Glastad et al., 2019).

8.2.3 Non-coding RNAs (ncRNAs)

Non-coding RNAs (ncRNAs) are a varied class of RNAs that are not translated into proteins. Some of ncRNA products may have no specific function, but others play a key role in regulating cellular processes (Glastad et al., 2019). There are four types of ncRNAs that have been suggested to have a potential epigenetic effect: PIWI-interacting RNAs (piRNAs), microRNAs (miRNAs), small interfering RNAs (siRNAs), and long noncoding RNAs (lncRNAs). The small regulatory RNAs (piRNAs, miRNAs and siRNAs) form the RNAi (RNA interference) pathway, which is responsible for RNA-based antiviral immunity (see “immunity” section). RNAi is a post-transcriptional sequence-specific gene silencing mechanism that is involved in regulating gene expression in most organisms [reviewed by Ding (2010)]. Among them, miRNAs and siRNAs are major post-transcriptional gene expression regulators (reviewed by Richard et al. (2021)]. On the other hand, there is evidence that piRNAs and lncRNAs are linked to epigenetic effects that are particularly strong in insects [reviewed by Chambeyron and Seitz (2014)].

The small regulatory piRNAs are highly variable, short (21–35 nucleotides), single-stranded ncRNAs that are the primary small RNAs mediating chromatin modifications within insect genomes [reviewed by Chambeyron and Seitz (2014)]. In various organisms they are associated with PIWI proteins, which are part of Argonaute proteins [reviewed by Luteijn and Ketting (2013)]. They act to silence transposable elements (TEs), which are ubiquitous in insect genomes and can cause damage through aberrant recombination events and deleterious mutations. Interestingly, some piRNAs are maternally transmitted to offspring (Chambeyron and Seitz, 2014; Glastad et al., 2019). Another class of small regulatory RNAs are miRNAs, which are one of the best-studied classes of ncRNAs, although molecular evidence that they influence chromatin or are passed on through cell divisions in a truly epigenetic manner is still largely lacking in insects (Glastad et al., 2019). Finally, the siRNAs interact with proteins from the same families of miRNAs. It is important to note that, in addition to their epigenetic properties, siRNAs constitute an important antiviral defence mechanism in plants, fungi, nematodes, and arthropods, and the role of RNAi in honey bee antiviral defence has been demonstrated in laboratory-based experiments (McMenamin et al., 2018). Lastly, lncRNAs are transcripts of more than 200 nucleotides that are present in almost all eukaryotic organisms (Mercer et al., 2009) and share the potential for regulating gene expression at both transcriptional and post-transcriptional levels (Engreitz et al., 2016). They have been found to function in developmental processes in honey bees (Choudhary et al., 2021).

9 Immunity

Immunity is one of the main physiological mechanisms that regulate the overall survival of an organism, and, therefore, immune stress requires a substantial nutritional and energetic cost [reviewed by Lochmiller and Deerenberg (2000)]. The immune system of honey bees consists of a complex set of individual immune mechanisms and a special type of behavioural adaptations that protect against biotic and abiotic stress factors. As social insects, at the colony level, honey bees exhibit a collective defence mechanism known as “social immunity,” while at the individual level, they rely on cellular and humoral immune reactions. Within humoral immunity, antimicrobial peptides (AMPs) play a key role (Danihlík et al., 2015). At the molecular level, several immune pathways have been described in A. mellifera. These are mainly Jak/STAT (Janus Kinase/Signal Transducer and Activator of Transcription), RNAi (RNA interference), Toll via NF-ĸΒ (Nuclear Factor ĸΒ/Dorsal), Imd (immune deficiency) via NF-ĸΒ/Relish, JNK (c-Jun N-terminal kinase), MAPK (Mitogen-Activated Protein Kinases), autophagy, eicosanoid biosynthesis, endocytosis, melanization, and Prophenoloxidase (PPO) (McMenamin et al., 2018) (Figure 7).

FIGURE 7

Diagram illustrating multiple cellular pathways and immune responses, including Jak/STAT, RNA interference, Toll, immune deficiency (Imd), autophagy, eicosanoid, and prophenoloxidase (PPO). Each pathway shows components and interactions like proteins and genetic elements, with icons indicating specific actions such as PAMP recognition, proteolytic cascades, and autophagy induction. Fungal and bacterial elements are depicted, showing recognition and response processes. Key molecules include Jak, STAT, RISC, Dicer-like siRNA, PGRPLC, NF-κB, PI3K, and PLA2, with additional references to arachidonic acid, defensins, and melanin.

Immunity pathways in the honey bee. The main immunity pathways are represented together with the biotic stressors that are able to activate them: Jak/STAT (Janus Kinase/Signal Transducer and Activator of Transcription), RNAi (RNA interference), Toll via NF-ĸΒ (Nuclear Factor ĸΒ/Dorsal), Imd (immune deficiency) via NF-ĸΒ/Relish, JNK (c-Jun N-terminal kinase), MAPK (Mitogen-Activated Protein Kinases), autophagy, eicosanoid biosynthesis, endocytosis, and melanization and Prophenoloxidase (PPO). Genes involved in abiotic stress induced by climate change are highlighted: red, green, white, purple, and yellow stars indicate an association with heat, nutrition, humidity, UV, and oxidative stress, respectively. Redrawn based on (McMenamin et al., 2018).

Most of the stressors described in this review, directly or indirectly can modify the expression of genes related to at least one of these pathways (Table 1). For example, the first step when pathogens enter honey bee is the host recognition pathogen-associated molecular patterns from the invading microorganisms (PAMP recognition), which rapidly promote the activation of an SP cascade (Figure 7) [reviewed by McMenamin et al. (2018)]. Serine proteases, which have been described related to heat stress (see section 4), are also directly involved in immunity, as they are key in melanisation, wound healing, and phagocytosis stimulation by participating in the prophenoloxidase (PPO) activation pathway, RNAi, and SP proteolytic cascade in the Toll signalling [reviewed by Lu et al. (2014), McMenamin et al. (2018)] (Figure 3). Three genes of SP putative substrates, PPO, spätzle-1, and spätzle-2, are described in the A. mellifera genome (Zou et al., 2006). However, Phenoloxidase has been also linked to heat stress, existing different levels of enzymatic activity in workers, queens and drones exposed to the same stressful situation during development (Medina et al., 2020).

Regarding heat, there is a close relationship between the iRNAs and HSPs in honey bees. The RNAi pathway is initiated by Dicer-2 cleavage of viral dsRNA into 21–22 bp siRNAs, which are then loaded into Argonaute-2 (Ago2), the catalytic component of the RNA Induced Silencing Complex (RISC). Within this route, a putative serine/threonine cyclin-dependent kinase (MF116383) acts in non-specific dsRNA-mediated antiviral responses [reviewed by McMenamin et al. (2018)] (Figure 7). McMenamin et al. (2020) conducted an in-depth study of A. mellifera, comparing the transcriptomic response to heat shock and viral infections. The mf116383 gene was the only immune-related gene of this pathway consistently induced by heat treatment alone and it has been suggested to serve as a point of cross-talk between the generalized antiviral immune response and the HSR in honey bees. On the other hand, the expression of Hsc70-3, Hsc70-4, and Hsp90 was found to be positively correlated with Dicer-like (Dcr-like) and Argonaute2 (Ago2), suggesting co-regulation of these genes (McMenamin et al., 2020). Cognate forms of Hsp70 are involved in regulating the life cycle of various viruses, such as mediating attachment and endocytosis (Iša et al., 2004), penetration and uncoating, transcription and replication (Du et al., 2011), assembly and budding (Prange et al., 1999), and modulating autophagy (Beere, 2004). HSC70-4 is an important chaperone for the assembly of the RISC in Drosophila S2 cells and other flies (Dorner et al., 2006; Iwasaki et al., 2010). While further studies are needed to determine the mechanisms leading to the co-regulation of immune genes and HSP-encoding genes, it may be beneficial to co-regulate HSPs and HSP client proteins (Iwasaki et al., 2010). Meanwhile, heat stress in honey bee colonies has been linked to a reduction in virus and parasite infections [reviewed by Zhao et al. (2021)]. Finally, RNAi is not only key to antiviral defence and closely related to HSR, but may also act as epigenetic factors. The RNAi pathways is constituted by piRNAs, miRNAs and siRNAs, and is involved in regulating gene expression in most organisms [reviewed by Ding (2010)] (see “epigenetics” subsection).

Within the Toll pathway (Figure 7), NF-ĸB transcriptional factors are crucial in immunity, inflammatory response, cellular adhesion, differentiation, proliferation, autophagy, senescence, and apoptosis (Bonizzi and Karin, 2004). In A. mellifera, NF-ĸB factors Dorsal are encoded by two dorsal homologues genes, Dorsal-1 and Dorsal-2 (Evans et al., 2006). Dorsal-1 produces two isoforms through alternative splicing (DORSAL-1A and DORSAL-1B) (Fan, 2001; Evans et al., 2006), and arguably this allows for a refinement of Toll-induced immune function. At this point, there is an important relationship between immunity and oxidative stress, consisting of both ROS influencing the activation of NF-κB pathway, and NF-κB pathway influencing the ROS levels by increasing the expression of antioxidant genes (Zhang et al., 2016). On the other hand, in adult bees, Dorsal genes are directly involved in the regulation of AMP genes together with Imd pathway (McMenamin et al., 2018), as Defensine1 (Lourenço et al., 2018) and Apidaecin (Prisco et al., 2016). The activation of both Toll and Imd pathways also governs the expression of Abaecin (McMenamin et al., 2018). Defensin genes have been described as predictive markers of honey bee colony health linked to viral infections (Barroso-Arévalo et al., 2019). Two Defensin genes are described in honey bees: Defensin1 and Defensin2. Sequences of these two genes revealed their different structure in a phylogenetic analysis, with Defensin1 forming a clade with two other hymenopteran defensins, whereas Defensin2 grouped with coleopteran defensins (Klaudiny et al., 2005). In contrast to other invertebrates, heat shock in A. mellifera represses the expression of certain antimicrobial peptide genes such as Hymenoptaecin, Defensin1 and Abaecin (McKinstry et al., 2017). This may be due to an energetic trade-off in the honeybee between cellular mechanisms for maintaining proteostasis and the immune response, as both are energetically costly (Powers and Balch, 2013).

The link between nutrition and immunity is well known. In fact, “nutritional immunology” is a discipline that studies the composition of the diet to optimise the immune response (Ponton et al., 2011). The interaction between nutritional and immune pathways is as wide as complex. For example, miRNAs have been described as acting on the insulin and TOR pathways (Shi et al., 2015), while PiP3, AKT and TOR inhibit autophagy (Figure 4) (McMenamin et al., 2018). On the other hand, inadequate nutrition stimulates adenylate kinase (Ak6) (involved in stress-induced pathways), NF-ĸB, and genotoxic/non-genotoxic stress (Hou et al., 2008; Barrio et al., 2014; Hasygar and Hietakangas, 2014; Corona et al., 2023). Furthermore, VG is involved in immunity in honey bees (Amdam et al., 2004). Recently, Azzouz-Olden et al. (2018) conducted an in-depth study comparing well- and poorly-fed honey bees between two groups comprising infected and non-infected individuals with the microsporidium Nosema apis. The authors showed that poorly fed and infected individuals underwent a significant upregulation of some genes related to the Toll pathway, such as Sph, serpine-1 (Nec), Tub, Cactus, PGRP-S3 and AMP genes (like Apisimin). On the other hand, during melanisation (see section 6), melanins are produced by the enzymatic oxidation of tyrosine by tyrosinase, which is the enzyme’s main activating system (PPO). Tyrosinase depends upon tyrosine, which derives from phenylalanine, an essential amino acid that can only be obtained from ingested food. In addition, the defensive compounds produced by the bees are rich in nitrogen, requiring a significant investment of their resources. Hence, the PO response involves a high energetic cost for each bee (González-Santoyo and Córdoba-Aguilar, 2012; Negri et al., 2019).

It is important to note that the PLA2 protein, involved in stress induced by humidity, is involved in immunity through the arachidonic acid release in the eicosanoid pathway (McMenamin et al., 2018).

10 Conclusion

Predicting how climate change will impact honey bees is complex and challenging. It will not only directly affect individual bees and colonies but also compound existing stressors, alter the ecosystems they inhabit, and introduce new and additional pressures. Higher temperatures, changes in humidity, and increased exposure to UV radiation will not only affect bee biology but also contribute to colony stress, alter plant flowering (and thus bee food resources), and influence disease dynamics. These factors can both benefit the pathogens under the new conditions and weaken the host organisms. In this review, we aimed to summarise a guide of molecular targets potentially linked to climate change stress, along with the cellular pathways directly involved or related. Genetic studies are a valuable tool for addressing these issues, and understanding the related mechanisms and pathways is essential for correctly interpreting molecular insights. Knowing which genes are involved in a specific pathway, what they are like, what their transcripts are like, how they express under stress, and how they affect other pathways, provides a basis for future studies aimed at measuring not only the degree of stress in honey bees but also its consequences for their cell biology.

Statements

Author contributions

SS: Conceptualization, Supervision, Writing – original draft, Writing – review and editing. GC: Supervision, Writing – original draft, Writing – review and editing, Conceptualization. DH: Writing – original draft, Supervision, Writing – review and editing. CY: Writing – review and editing, Writing – original draft. MC: Conceptualization, Writing – original draft, Writing – review and editing. MH: Supervision, Writing – review and editing, Writing – original draft. MP: Conceptualization, Writing – review and editing, Funding acquisition, Supervision, Writing – original draft. AN: Writing – original draft, Funding acquisition, Writing – review and editing, Supervision. RM-H: Supervision, Funding acquisition, Writing – review and editing, Conceptualization, Project administration, Writing – original draft.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. MEDIBEES project (Grant number 2011) funded by PRIMA (Partnership for research and innovation in the Mediterranean area) programme supported by the European Union. Fundação para a Ciência e a Tecnologia (FCT), Portugal, funded M.A.P., D.H. and C.Y.G. through national funds to CIMO, under UIDB/00690/2020 (DOI: 10.54499/UIDB/00690/2020) and UIDP/00690/2020 (DOI: 10.54499/UIDP/00690/2020), and SusTEC, under LA/P/0007/2020 (DOI: 10.54499/LA/P/0007/2020), via FCT/MCTES (PIDDAC). FCT also funded C.Y.G. through the individual PhD research grant 2021.06948.BD (DOI: 10.54499/2021.06948.BD).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2025.1623705/full#supplementary-material

Glossary

  • AChE

    Acetylcholinesterase

  • Ace

    Acetylcholinesterase gene

  • AMPs

    antimicrobial peptides

  • AO

    Antioxidant

  • AS

    Alternative splicing

  • CpG

    Cytosine-Guanine dinucleotide

  • CRH-BP

    Cortico-releasing hormone-binding protein

  • HDAC

    Histone Deacetylase

  • HDACis

    Histone Deacetylase inhibitors

  • DNMTs

    DNA methyltransferases

  • hPTMs

    Histone post-translational modifications

  • HS

    Heat Shock

  • HSC

    Heat Shock Protein cognate form

  • HSF

    Heat Shock Factor

  • HSE

    Heat Shock Element

  • HSR

    Heat Shock Response

  • HSP

    Heat Shock Proteins

  • FOXO

    Forkhead box O subfamily of genes

  • GSH

    glutathione

  • GST

    glutathione S-transferase

  • IIS

    Insulin/insulin-like growth factor signalling

  • RNAi

    RNA of interference

  • JH

    juvenile hormone

  • ncRNAs

    Non-coding RNAs

  • NF-Y

    Nuclear Factor Y

  • OA

    Octopamine

  • PPO

    prophenoloxidase

  • PO

    phenoloxidase

  • RH

    Relative humidity

  • ROS

    reactive oxygen species

  • sHSP

    small Heat Shock Proteins

  • SP

    Serine Protease

  • STK

    Serine/threonine kinase

  • SNP

    Single Nucleotide Polymorphism

  • SOD

    Superoxide dismutase

  • TA

    Tyramine

  • TEs

    Transposable elements

  • TOR

    Target of Rapamycin

  • TRX

    thioredoxin

  • UV

    Ultraviolet

  • VG

    Vitellogenin

  • ZFP

    Zinc Finger Protein

References

  • 1

    AamodtR. M. (2008). The caste- and age-specific expression signature of honeybee heat shock genes shows an alternative splicing-dependent regulation of Hsp90. Mech. Ageing Dev.129 (11), 632637. 10.1016/j.mad.2008.07.002

  • 2

    Abou-ShaaraH. F. (2024). The response of heat shock proteins in honey bees to abiotic and biotic stressors. J. Therm. Biol.119, 103784. 10.1016/j.jtherbio.2024.103784

  • 3

    Abou-ShaaraH. F.Al-GhamdiA. A.MohamedA. A. (2012). Tolerance of two honey bee races to various temperature and relative humidity gradients. Environ. Exp. Biol.10. 10.1371/journal.pone.0217921

  • 4

    Abou-ShaaraH. F.OwayssA. A.IbrahimY. Y.BasunyN. K. (2017). A review of impacts of temperature and relative humidity on various activities of honey bees. Insectes Sociaux. Birkhauser Verl. Ag.64, 455463. 10.1007/s00040-017-0573-8

  • 5

    AbrahamR. T. (2009). Regulation of the mTOR signaling pathway: from laboratory bench to bedside and back again. F1000 Biol. Rep.1, 8. 10.3410/B1-8

  • 6

    Al-GhamdiA. A.Abou-ShaaraH. F.MohamedA. A. (2014). Hatching rates and some characteristics of Yemeni and Carniolan honey bee eggs. ∼ 6 ∼. J. Entomol. Zoology Stud.2 (1), 610.

  • 7

    Al-GhzawiA. A. M. A.Al-ZghoulM. B.ZaitounS.Al-OmaryI. M.AlahmadN. A. (2022). Dynamics of heat shock proteins and heat shock factor expression during heat stress in daughter workers in pre-heat-treated (rapid heat hardening) Apis mellifera mother queens. J. Therm. Biol.104, 103194. 10.1016/j.jtherbio.2022.103194

  • 8

    AlattalY.Al-GhamdiA.AlghamdiA. (2015). Impact of temperature extremes on survival of indigenous and exotic honeybee subspecies, Apis mellifera, under desert and semiarid climates Impact of temperature extremes on survival of indigenous and exotic honey bee subspecies, Apis mellifera, under desert and semiarid climates. Bull. Insectol.68 (2), 219222.

  • 9

    AlauxC.DuclozF.CrauserD.Le ConteY. (2010). Diet effects on honeybee immunocompetence. Biol. Lett.6 (4), 562565. 10.1098/rsbl.2009.0986

  • 10

    AlonM.AlonF.NauenR.MorinS. (2008). ‘Organophosphates’ resistance in the B-biotype of Bemisia tabaci (Hemiptera: aleyrodidae) is associated with a point mutation in an ace1-type acetylcholinesterase and overexpression of carboxylesterase. Insect Biochem. Mol. Biol.38 (10), 940949. 10.1016/j.ibmb.2008.07.007

  • 11

    AlqarniA. S. (2020). Differential foraging of indigenous and exotic honeybee (Apis mellifera L.) races on nectar-rich flow in a subtropical ecosystem. Insects11 (4), 254. 10.3390/insects11040254

  • 12

    AlqarniA. S.AliH.IqbalJ.OwayssA. A.SmithB. H. (2019). Expression of heat shock proteins in adult honey bee (Apis mellifera L.) workers under hot-arid subtropical ecosystems. Saudi J. Biol. Sci.26 (7), 13721376. 10.1016/j.sjbs.2019.08.017

  • 13

    AmdamG. V.OmholtS. W. (2002). The regulatory anatomy of honeybee lifespan. J. Theor. Biol.216 (2), 209228. 10.1006/jtbi.2002.2545

  • 14

    AmdamG. V.NorbergK.HagenA.OmholtS. W. (2003). Social exploitation of vitellogenin. Proc. Natl. Acad. Sci.100 (4), 17991802. 10.1073/pnas.0333979100

  • 15

    AmdamG. V.SimõesZ. L. P.HagenA.NorbergK.SchrøderK.MikkelsenØ.et al (2004). Hormonal control of the yolk precursor vitellogenin regulates immune function and longevity in honeybees. Exp. Gerontol.39 (5), 767773. 10.1016/j.exger.2004.02.010

  • 16

    AmentS. A.CoronaM.PollockH. S.RobinsonG. E. (2008). Insulin signaling is involved in the regulation of worker division of labor in honey bee colonies. Proc. Natl. Acad. Sci. U. S. A.105 (11), 42264231. 10.1073/pnas.0800630105

  • 17

    AnandL.Rodriguez LopezC. M. (2022). ChromoMap: an R package for interactive visualization of multi-omics data and annotation of chromosomes. BMC Bioinforma.23 (1), 33. 10.1186/s12859-021-04556-z

  • 18

    AshidaM.BreyP. T. (1997). Recent advances in research on the insect prophenoloxidase cascade. Molecular mechanisms of immune responses in insects. Editor BreyP. T.HultmarkD. (London: Chapman and Hall). 135172.

  • 19

    AstG. (2004). How did alternative splicing evolve?Nat. Rev. Genet.5, 773782. 10.1038/nrg1451

  • 20

    AtmowidjojoA. H.WheelerD. E.EricksonE. H.CohenA. C. (1997). Temperature tolerance and water balance in feral and domestic honey bees, Apis mellifera L. Comparative Biochemistry and Physiology Part A: Physiology118 (4), 13991403. 10.1016/S0300-9629(97)00031-5

  • 21

    Azzouz-OldenF.HuntA.DeGrandi-HoffmanG. (2018). Transcriptional response of honey bee (Apis mellifera) to differential nutritional status and Nosema infection. BMC Genomics19 (1), 628. 10.1186/s12864-018-5007-0

  • 22

    BaekJ. H.KimJ. I.LeeD. W.ChungB. K.MiyataT.LeeS. H. (2005). Identification and characterization of ace1-type acetylcholinesterase likely associated with organophosphate resistance in Plutella xylostella. Pesticide Biochem. Physiol.81 (3), 164175. 10.1016/j.pestbp.2004.12.003

  • 23

    BalfanzS.JordanN.LangenstückT.BreuerJ.BergmeierV.BaumannA. (2014). Molecular, pharmacological, and signaling properties of octopamine receptors from honeybee (Apis mellifera) brain. J. Neurochem.129 (2), 284296. 10.1111/jnc.12619

  • 24

    BarekH.VeraksaA.SugumaranM. (2018). Drosophila melanogaster has the enzymatic machinery to make the melanic component of neuromelanin. Pigment Cell Melanoma Res.31 (6), 683692. 10.1111/pcmr.12709

  • 25

    BarrioL.DekantyA.MilánM. (2014). MicroRNA-mediated regulation of Dp53 in the Drosophila fat body contributes to metabolic adaptation to nutrient deprivation. Cell Rep.8 (2), 528541. 10.1016/j.celrep.2014.06.020

  • 26

    Barroso-ArévaloS.Vicente-RubianoM.PuertaF.MoleroF.Sánchez-VizcaínoJ. M. (2019). Immune related genes as markers for monitoring health status of honey bee colonies. BMC Veterinary Res.15 (1), 72. 10.1186/s12917-019-1823-y

  • 27

    BatagliaL.SimõesZ. L. P.NunesF. M. F. (2021). Active genic machinery for epigenetic RNA modifications in bees. Insect Mol. Biol.30 (6), 566579. 10.1111/imb.12726

  • 28

    BatagliaL.SimõesZ. L. P.NunesF. M. F. (2022). Transcriptional expression of m6A and m5C RNA methyltransferase genes in the brain and fat body of honey bee adult workers. Front. Cell Dev. Biol.10, 921503. 10.3389/fcell.2022.921503

  • 29

    BecherM. A.ScharpenbergH.MoritzR. F. A. (2009). Pupal developmental temperature and behavioral specialization of honeybee workers (Apis mellifera L.). J. Comp. Physiology A195, 673679. 10.1007/s00359-009-0442-7

  • 30

    BeereH. M. (2004). “The stress of dying”: the role of heat shock proteins in the regulation of apoptosis. J. Cell Sci.117 (13), 26412651. 10.1242/jcs.01284

  • 31

    BeggsK. T.TyndallJ. D. A.MercerA. R. (2011). Honey bee dopamine and octopamine receptors linked to intracellular calcium signaling have a close phylogenetic and pharmacological relationship. PLoS One6 (11), e26809. 10.1371/journal.pone.0026809

  • 32

    Belhadj SlimenI.NajarT.GhramA.DabbebiH.Ben MradM.AbdrabbahM. (2014). Reactive oxygen species, heat stress and oxidative-induced mitochondrial damage. A review. Int. J. Hyperth.30 (7), 513523. 10.3109/02656736.2014.971446

  • 33

    BellefroidE. J.MarineJ. C.RiedT.LecocqP. J.RivièreM.AmemiyaC.et al (1993). Clustered organization of homologous KRAB zinc-finger genes with enhanced expression in human T lymphoid cells. EMBO J.12 (4), 13631374. 10.1002/j.1460-2075.1993.tb05781.x

  • 34

    BergerS. L.KouzaridesT.ShiekhattarR.ShilatifardA. (2009). An operational definition of epigenetics: Figure 1. Genes Dev.23 (7), 781783. 10.1101/gad.1787609

  • 35

    BernhardG. H.NealeR. E.BarnesP. W.NealeP. J.ZeppR. G.WilsonS. R.et al (2020). Environmental effects of stratospheric ozone depletion, UV radiation and interactions with climate change: UNEP Environmental Effects Assessment Panel, update 2019. Photochem. Photobiological Sci.19 (5), 542584. 10.1039/d0pp90011g

  • 36

    BernsteinB. E.MeissnerA.LanderE. S. (2007). The mammalian epigenome. Cell128 (4), 669681. 10.1016/j.cell.2007.01.033

  • 37

    BiamontiG.CaceresJ. F. (2009). Cellular stress and RNA splicing. Trends Biochem. Sci.34 (3), 146153. 10.1016/j.tibs.2008.11.004

  • 38

    BitondiM. M. G.SimoesZ. L. P. (1996). The relationship between level of pollen in the diet, vitellogenin and juvenile hormone titres in Africanized Apis mellifera workers. J. Apic. Res.35 (1), 2736. 10.1080/00218839.1996.11100910

  • 39

    BlenauW.BaumannA. (2001). Molecular and pharmacological properties of insect biogenic amine receptors: lessons from Drosophila melanogaster and Apis mellifera. Archives Insect Biochem. Physiol.48 (1), 1338. 10.1002/arch.1055

  • 40

    BlenauW.BalfanzS.BaumannA. (2000). Amtyr1: characterization of a gene from honeybee (Apis mellifera) brain encoding a functional tyramine receptor. J. Neurochem.74 (3), 900908. 10.1046/j.1471-4159.2000.0740900.x

  • 41

    BlencoweB. J. (2006). Alternative splicing: new insights from global analyses. Cell126 (1), 3747. 10.1016/j.cell.2006.06.023

  • 42

    BonizziG.KarinM. (2004). The two NF-kappaB activation pathways and their role in innate and adaptive immunity. Trends Immunol.25 (6), 280288. 10.1016/j.it.2004.03.008

  • 43

    BordierC.SuchailS.PiozM.DevaudJ. M.ColletC.CharretonM.et al (2017). Stress response in honeybees is associated with changes in task-related physiology and energetic metabolism. J. Insect Physiol.98, 4754. 10.1016/j.jinsphys.2016.11.013

  • 44

    BrodschneiderR.CrailsheimK. (2010). Nutrition and health in honey bees. Apidologie41 (3), 278294. 10.1051/apido/2010012

  • 45

    BruntS. A.SilverJ. C. (1991). Molecular cloning and characterization of two distinct hsp 85 sequences from the steroid responsive fungus Achlya ambisexualis. Curr. Genet.19 (5), 383388. 10.1007/BF00309599

  • 46

    CalderwoodS. K.CioccaD. R. (2008). Heat shock proteins: stress proteins with Janus-like properties in cancer. Int. J. Hyperth.24 (1), 3139. 10.1080/02656730701858305

  • 47

    de CarvalhoC. F. (2023). Epigenetic effects of climate change on insects. Curr Opin Insect Sci.57, 101029. 10.1016/j.cois.2023.101029

  • 48

    ChambeyronS.SeitzH. (2014). Insect small non-coding RNA involved in epigenetic regulations. Current opinion in insect science1, 19. 10.1016/j.cois.2014.05.001

  • 49

    CharlierN.DesoilM.GossuinY.GillisP.GallezB. (2020). Electron paramagnetic resonance imaging of melanin in honey bee. Cell Biochem. Biophysics78 (2), 123126. 10.1007/s12013-020-00903-8

  • 50

    Chávez‐GalarzaJ.HenriquesD.JohnstonJ. S.AzevedoJ. C.PattonJ. C.MuñozI.et al (2013). Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms. Mol. Ecol.22 (23), 58905907. 10.1111/mec.12537

  • 51

    ChenC.LiuZ.PanQ.ChenX.WangH.GuoH.et al (2016). Genomic analyses reveal demographic history and temperate adaptation of the newly discovered honey bee subspecies Apis mellifera sinisxinyuan n. ssp. Mol. Biol. Evol.33 (5), 13371348. 10.1093/molbev/msw017

  • 52

    ChoudharyC.SharmaS.MeghwanshK. K.PatelS.MehtaP.ShuklaN.et al (2021). Long non-coding rnas in insects. Animals11 (4), 1118. 10.3390/ani11041118

  • 53

    ChuangK.-H.HoS.-H.SongY.-L. (2007). Cloning and expression analysis of heat shock cognate 70 gene promoter in tiger shrimp (Penaeus monodon). Gene405 (1–2), 1018. 10.1016/j.gene.2007.08.016

  • 54

    CingolaniP.CaoX.KhetaniR. S.ChenC. C.CoonM.SammakA.et al (2013). Intronic Non-CG DNA hydroxymethylation and alternative mRNA splicing in honey bees. BMC Genomics14 (1), 666. 10.1186/1471-2164-14-666

  • 55

    CockellC. S.KnowlandJ. (1999). Ultraviolet radiation screening compounds. Biol. Rev. Camb. Philosophical Soc.74 (3), 311345. 10.1017/S0006323199005356

  • 56

    CollottaM.BertazziP. A.BollatiV. (2013). Epigenetics and pesticides. Toxicology307, 3541. 10.1016/j.tox.2013.01.017

  • 57

    CoronaM.RobinsonG. E. (2006) Genes of the antioxidant system of the honey bee: annotation and phylogeny. Insect Mol Biol.15 (5), 687701. 10.1111/j.1365-2583.2006.00695.x

  • 58

    CoronaM.VelardeR. A.RemolinaS.Moran-LauterA.WangY.HughesK. A.et al (2007). Vitellogenin, juvenile hormone, insulin signaling, and queen honey bee longevity. Proc. Natl. Acad. Sci. U. S. A.104 (17), 71287133. 10.1073/pnas.0701909104

  • 59

    CoronaM.LibbrechtR.WheelerD. E. (2016). Molecular mechanisms of phenotypic plasticity in social insects. Current opinion in insect science13, 5560. 10.1016/j.cois.2015.12.003

  • 60

    CoronaM.BranchiccelaB.AlburakiM.Palmer-YoungE. C.MadellaS.ChenY.et al (2023). Decoupling the effects of nutrition, age, and behavioral caste on honey bee physiology, immunity, and colony health. Front. Physiol.14, 1149840. 10.3389/fphys.2023.1149840

  • 61

    CortrightD. N.GoosensK. A.LeshJ. S.SeasholtzA. F. (1997). Isolation and characterization of the rat corticotropin-releasing hormone (CRH)-binding protein gene: transcriptional regulation by cyclic adenosine monophosphate and CRH. Endocrinology138 (5), 20982108. 10.1210/endo.138.5.5128

  • 62

    CraigE. A.MarszalekJ. (2017). How do J-proteins get Hsp70 to do so many different things?Trends Biochem. Sci.42 (5), 355368. 10.1016/j.tibs.2017.02.007

  • 63

    CraigE. A.HuangP.AronR.AndrewA. (2006). “The diverse roles of J-proteins, the obligate Hsp70 co-chaperone,” in Reviews of physiology, biochemistry and pharmacology (Berlin, Heidelberg: Springer Berlin Heidelberg). 10.1007/s10254-005-0001-0

  • 64

    DaltonR. M.UnderwoodN. C.InouyeD. W.SouléM. E.InouyeB. D. (2023). Long‐term declines in insect abundance and biomass in a subalpine habitat. Ecosphere14 (8), e4620. 10.1002/ecs2.4620

  • 65

    DanihlíkJ.AronsteinK.PetřivalskýM. (2015). Antimicrobial peptides: a key component of honey bee innate immunity. J. Apic. Res.54, 123136. 10.1080/00218839.2015.1109919

  • 66

    DasguptaS.HohmanT. C.CarperD. (1992). Hypertonic stress induces αB-crystallin expression. Exp. Eye Res.54 (3), 461470. 10.1016/0014-4835(92)90058-Z

  • 67

    DeGrandi-HoffmanG.Corby-HarrisV.CarrollM.TothA. L.GageS.Watkins deJongE.et al (2021). The importance of time and place: nutrient composition and utilization of seasonal pollens by European honey bees (Apis mellifera l.). Insects12 (3), 235. 10.3390/insects12030235

  • 68

    DingS.-W. (2010). RNA-based antiviral immunity. Nat. Rev. Immunol.10 (9), 632644. 10.1038/nri2824

  • 69

    DingF.GilM. P.FranklinM.FerreiraJ.TatarM.HelfandS. L.et al (2014). Transcriptional response to dietary restriction in Drosophila melanogaster. J. Insect Physiol.69, 101106. 10.1016/j.jinsphys.2014.05.002

  • 70

    DornerS.LumL.KimM.ParoR.BeachyP. A.GreenR. (2006). A genomewide screen for components of the RNAi pathway in Drosophila cultured cells. Proc. Natl. Acad. Sci. U. S. A.103 (32), 1188011885. 10.1073/pnas.0605210103

  • 71

    DrollD.MiniaI.FaddaA.SinghA.StewartM.QueirozR.et al (2013). Post-transcriptional regulation of the trypanosome heat shock response by a zinc finger protein. PLoS Pathog.9 (4), e1003286. 10.1371/journal.ppat.1003286

  • 72

    DuN.-N.LiX.WangY. P.LiuF.LiuY. X.LiC. X.et al (2011). Synthesis, structure−activity relationship and biological evaluation of novel N-substituted matrinic acid derivatives as host heat-stress cognate 70 (Hsc70) down-regulators. Bioorg. Med. Chem. Lett.21 (16), 47324735. 10.1016/j.bmcl.2011.06.071

  • 73

    DuncanE. J.LeaskM. P.DeardenP. K. (2020). Genome architecture facilitates phenotypic plasticity in the honeybee (Apis mellifera). Mol. Biol. Evol.37 (7), 19641978. 10.1093/molbev/msaa057

  • 74

    EleftherianosI.FelföldiG.ffrench-ConstantR. H.ReynoldsS. E. (2009). Induced nitric oxide synthesis in the gut of Manduca sexta protects against oral infection by the bacterial pathogen Photorhabdus luminescens. Insect Mol. Biol.18 (4), 507516. 10.1111/j.1365-2583.2009.00899.x

  • 75

    ElekonichM. M. (2009). Extreme thermotolerance and behavioral induction of 70-kDa heat shock proteins and their encoding genes in honey bees. Cell Stress Chaperones14 (2), 219226. 10.1007/s12192-008-0063-z

  • 76

    EllisM. B.NicolsonS. W.CreweR. M.DietemannV. (2008). Hygropreference and brood care in the honeybee (Apis mellifera). J. Insect Physiol.54 (12), 15161521. 10.1016/j.jinsphys.2008.08.011

  • 77

    EmlenD. J.WarrenI. A.JohnsA.DworkinI.LavineL. C. (2012). A mechanism of extreme growth and reliable signaling in sexually selected ornaments and weapons. Science337 (6096), 860864. 10.1126/science.1224286

  • 78

    EngreitzJ.OllikainenN.GuttmanM. (2016). Long non-coding RNAs: spatial amplifiers that control nuclear structure and gene expression. Nat. Rev. Mol. Cell Biol.17, 756770. 10.1038/nrm.2016.126

  • 79

    EvansJ. D.AronsteinK.ChenY. P.HetruC.ImlerJ. L.JiangH.et al (2006). Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol. Biol.15, 645656. 10.1111/j.1365-2583.2006.00682.x

  • 80

    EvenN.DevaudJ. M.BarronA. B. (2012). General stress responses in the honey bee. Insects3, 12711298. 10.3390/insects3041271

  • 81

    FanX. (2001). Apis dorsal: gene structure, role in oogenesis and embryonic development, signal in the induction of an innate immune response. Tübingen: Erberhard-Karls Universität Tübingen.

  • 82

    FewellJ. H.WinstonM. L. (1992). Colony state and regulation of pollen foraging in the honey bee, Apis mellifera L. Behav. Ecol. Sociobiol.30, 387393. 10.1007/bf00176173

  • 83

    FinkelT.HolbrookN. J. (2000). Oxidants, oxidative stress and the biology of ageing. Nature408 (6809), 239247. 10.1038/35041687

  • 84

    FlyBase (2025) A database of Drosophila genes and genomes. Available online at: https://flybase.org/.

  • 85

    ForetS.KucharskiR.PittelkowY.LockettG. A.MaleszkaR. (2009). Epigenetic regulation of the honey bee transcriptome: unravelling the nature of methylated genes. BMC Genomics10 (1), 472. 10.1186/1471-2164-10-472

  • 86

    FujikakeN.NagaiY.PopielH. A.KanoH.YamaguchiM.TodaT. (2005). Alternative splicing regulates the transcriptional activity of Drosophila heat shock transcription factor in response to heat/cold stress. FEBS Lett.579 (17), 38423848. 10.1016/j.febslet.2005.05.074

  • 87

    FurayA. R.MurphyE. K.MattsonM. P.GuoZ.HermanJ. P. (2006). Region-specific regulation of glucocorticoid receptor/HSP90 expression and interaction in brain. J. Neurochem.98 (4), 11761184. 10.1111/j.1471-4159.2006.03953.x

  • 88

    FussneckerB. L.SmithB. H.MustardJ. A. (2006). Octopamine and tyramine influence the behavioral profile of locomotor activity in the honey bee (Apis mellifera). J. Insect Physiol.52 (10), 10831092. 10.1016/j.jinsphys.2006.07.008

  • 89

    GaoL.WangL.YangX.WangY.LiuZ.XuB.et al (2019). Role of a serine protease gene (AccSp1) from Apis cerana cerana in abiotic stress responses and innate immunity. Cell Stress Chaperones24 (1), 2943. 10.1007/s12192-018-0934-x

  • 90

    GarbuzD. G. (2017). Regulation of heat shock gene expression in response to stress. Mol. Biol.51, 352367. 10.1134/S0026893317020108

  • 91

    GérardM.VanderplanckM.WoodT.MichezD. (2020). Global warming and plant-pollinator mismatches. Emerg. Top. Life Sci.4 (1), 7786. 10.1042/ETLS20190139

  • 92

    GilloolyJ. F.BrownJ. H.WestG. B.SavageV. M.CharnovE. L. (2001). Effects of size and temperature on metabolic rate. Science293 (5538), 22482251. 10.1126/science.1061967

  • 93

    GkouvitsasT.KontogiannatosD.KourtiA. (2009). Cognate Hsp 70 gene is induced during deep larval diapause in the moth Sesamia nonagrioides. Insect Mol. Biol.18 (2), 253264. 10.1111/j.1365-2583.2009.00866.x

  • 94

    GlastadK. M.HuntB. G.GoodismanM. A. D. (2019). Epigenetics in insects: genome regulation and the generation of phenotypic diversity. Annu. Rev. Entomol.64, 185203. 10.1146/annurev-ento-011118-111914

  • 95

    González-SantoyoI.Córdoba-AguilarA. (2012). Phenoloxidase: a key component of the insect immune system. Entomologia Exp. Appl.142 (1), 116. 10.1111/j.1570-7458.2011.01187.x

  • 96

    González‐TokmanD.Córdoba-AguilarA.DáttiloW.Lira-NoriegaA.Sánchez-GuillénR. A.VillalobosF. (2020). Insect responses to heat: physiological mechanisms, evolution and ecological implications in a warming world. Biol. Rev.95 (3), 802821. 10.1111/brv.12588

  • 97

    Google Scholar (2025). Google scholar. Available online at: https://scholar.google.com/ (Accessed June 18, 2025).

  • 98

    GoulsonD.NichollsE.BotíasC.RotherayE. L. (2015). Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science347 (6229), 1255957. 10.1126/science.1255957

  • 99

    GreenopA.Mica-HawkyardN.WalkingtonS.WilbyA.CookS. M.PywellR. F.et al (2020). Equivocal evidence for colony level stress effects on bumble bee pollination services. Insects11 (3), 191. 10.3390/insects11030191

  • 100

    GrohmannL.BlenauW.ErberJ.EbertP. R.StrünkerT.BaumannA. (2003). Molecular and functional characterization of an octopamine receptor from honeybee (Apis mellifera) brain. J. Neurochem.86 (3), 725735. 10.1046/j.1471-4159.2003.01876.x

  • 101

    GuoH. juanWangL. J.WangC.GuoD. Z.XuB. H.GuoX. Q.et al (2021). Identification of an Apis cerana zinc finger protein 41 gene and its involvement in the oxidative stress response. Archives Insect Biochem. Physiol.108 (1), e21830. 10.1002/arch.21830

  • 102

    HammerM.MenzelR. (1995). Learning and memory in the honeybee. J. Neurosci.15 (3), 16171630. 10.1523/JNEUROSCI.15-03-01617.1995

  • 103

    HardstoneM. C.ScottJ. G. (2010). Is Apis mellifera more sensitive to insecticides than other insects?Pest Manag. Sci.66 (11), 11711180. 10.1002/ps.2001

  • 104

    HasygarK.HietakangasV. (2014). p53-and ERK7-dependent ribosome surveillance response regulates Drosophila insulin-like peptide secretion. PLoS Genet.10 (11), e1004764. 10.1371/journal.pgen.1004764

  • 105

    HeinrichB. (1980). Mechanisms of body-temperature regulation in honeybees, Apis mellifera. J. Exp. Biol.85 (1), 6172. 10.1242/jeb.85.1.61

  • 106

    HendrickJ. P.HartlF.-U. (1993). Molecular chaperone functions of heat-shock proteins. Annu. Rev. Biochem.62 (1), 349384. 10.1146/annurev.bi.62.070193.002025

  • 107

    HenriquesD.WallbergA.Chávez-GalarzaJ.JohnstonJ. S.WebsterM. T.PintoM. A. (2018). Whole genome SNP-associated signatures of local adaptation in honeybees of the Iberian Peninsula. Sci. Rep.8 (1), 11145. 10.1038/s41598-018-29469-5

  • 108

    HouY.-C. C.ChittaranjanS.BarbosaS. G.McCallK.GorskiS. M. (2008). Effector caspase Dcp-1 and IAP protein Bruce regulate starvation-induced autophagy during Drosophila melanogaster oogenesis. J. Cell Biol.182 (6), 11271139. 10.1083/jcb.200712091

  • 109

    HuY. T.WuT. C.YangE. C.WuP. C.LinP. T.WuY. L. (2017). Regulation of genes related to immune signaling and detoxification in Apis mellifera by an inhibitor of histone deacetylation. Sci. Rep.7, 41255. 10.1038/srep41255

  • 110

    HughesL. (2000). Biological consequences of global warming: is the signal already apparent?Trends Ecol. Evol.15 (2), 5661. 10.1016/S0169-5347(99)01764-4

  • 111

    HumanH.NicolsonS. W.DietemannV. (2006). Do honeybees, Apis mellifera scutellata, regulate humidity in their nest?Naturwissenschaften93 (8), 397401. 10.1007/s00114-006-0117-y

  • 112

    HwangboD. S.GershmanB.TuM. P.PalmerM.TatarM. (2004). Drosophila dFOXO controls lifespan and regulates insulin signalling in brain and fat body. Nature429 (6991), 562566. 10.1038/nature02549

  • 113

    InouyeD. W. (2020). Effects of climate change on alpine plants and their pollinators. Ann. N. Y. Acad. Sci.1469 (1), 2637. 10.1111/nyas.14104

  • 114

    IPCC-Intergovernmental Panel on Climate Change (2007) The intergovernmental panel on climate change, 52. World Meteorological Organization. Available online at: https://www.ipcc.ch/ (Accessed: June 18, 2025).

  • 115

    IšaP.RealpeM.RomeroP.LópezS.AriasC. F. (2004). Rotavirus RRV associates with lipid membrane microdomains during cell entry. Virology322 (2), 370381. 10.1016/j.virol.2004.02.018

  • 116

    ItoS.WakamatsuK. (2008). Chemistry of mixed melanogenesis—pivotal roles of dopaquinone. Photochem. Photobiol.84 (3), 582592. 10.1111/j.1751-1097.2007.00238.x

  • 117

    ItohH.TashimaY.EishiY.OkedaR. (1993). Localization of HSP90 in rat brain. Int. J. Biochem.25 (1), 9399. 10.1016/0020-711X(93)90494-Y

  • 118

    IwasakiS.KobayashiM.YodaM.SakaguchiY.KatsumaS.SuzukiT.et al (2010). Hsc70/Hsp90 chaperone machinery mediates ATP-dependent RISC loading of small RNA duplexes. Mol. Cell39 (2), 292299. 10.1016/j.molcel.2010.05.015

  • 119

    IzumotoS.HerbertJ. (1993). Widespread constitutive expression of HSP90 messenger RNA in rat brain. J. Neurosci. Res.35 (1), 2028. 10.1002/jnr.490350104

  • 120

    JeltschA.JurkowskaR. Z. (2014). New concepts in DNA methylation. Trends Biochem. Sci.39 (7), 310318. 10.1016/j.tibs.2014.05.002

  • 121

    JiaZ.DaviesP. L. (2002). Antifreeze proteins: an unusual receptor–ligand interaction. Trends Biochem. Sci.27 (2), 101106. 10.1016/S0968-0004(01)02028-X

  • 122

    JiangX.QuM.DenholmI.FangJ.JiangW.HanZ. (2009). Mutation in acetylcholinesterase1 associated with triazophos resistance in rice stem borer, Chilo suppressalis (Lepidoptera: Pyralidae). Biochem. Biophysical Res. Commun.378 (2), 269272. 10.1016/j.bbrc.2008.11.046

  • 123

    JohnsonJ. L.BrownC. (2009). Plasticity of the Hsp90 chaperone machine in divergent eukaryotic organisms. Cell Stress Chaperones14 (1), 8394. 10.1007/s12192-008-0058-9

  • 124

    JohnstonB. A.HooksK. B.McKinstryM.SnowJ. W. (2016). Divergent forms of endoplasmic reticulum stress trigger a robust unfolded protein response in honey bees. J. Insect Physiol.86, 110. 10.1016/j.jinsphys.2015.12.004

  • 125

    JonesJ. C.OldroydB. P. (2006). Nest thermoregulation in social insects. Adv. Insect Phys.153191. 10.1016/S0065-2806(06)33003-2

  • 126

    JoshiN. C.JoshiP. C. (2010). Foraging behaviour of Apis spp. on apple flowers in a subtropical environment. New York Science Journal3, 7176. Available online at: https://www.researchgate.net/publication/283046029.

  • 127

    KagiasK.NehammerC.PocockR. (2012). Neuronal responses to physiological stress. Front. Genet.3, 222. 10.3389/fgene.2012.00222

  • 128

    KampingaH. H.CraigE. A. (2010). The HSP70 chaperone machinery: J proteins as drivers of functional specificity. Nat. Rev. Mol. Cell Biol.11 (8), 579592. 10.1038/nrm2941

  • 129

    KanR. L.ChenJ.SallamT. (2022). Crosstalk between epitranscriptomic and epigenetic mechanisms in gene regulation. Trends Genet.38(2), 182193. 10.1016/j.tig.2021.06.014

  • 130

    KanostM. R. (1999). Serine proteinase inhibitors in arthropod immunity. Dev. Comp. Immunol.23 (4–5), 291301. 10.1016/S0145-305X(99)00012-9

  • 131

    KaremakerI. D.BaubecT. (2020). DNA methyltransferases hitchhiking on chromatin. Swiss Med. Wkly.150, w20329. 10.4414/smw.2020.20329

  • 132

    KetchesinK. D.StinnettG. S.SeasholtzA. F. (2017). Corticotropin-releasing hormone-binding protein and stress: from invertebrates to humans. Stress20 (5), 449464. 10.1080/10253890.2017.1322575

  • 133

    KhedidjiH.ChellaliM.CherreredN.GhaffarI.Ladjali-MohammediK.MohammediA. (2024). Effect of heat stress on Apis mellifera intermissa and Apis mellifera sahariensis in laboratory controlled conditions: a physiological approach. J. Apic. Res.64, 12431251. 10.1080/00218839.2024.2398903

  • 134

    KimY. H.LeeS. H. (2013). Which acetylcholinesterase functions as the main catalytic enzyme in the Class Insecta?Insect Biochem. Mol. Biol.43 (1), 4753. 10.1016/j.ibmb.2012.11.004

  • 135

    KimY.-C.LeeH. G.LimJ.HanK. A. (2013). Appetitive learning requires the alpha1-like octopamine receptor OAMB in the Drosophila mushroom body neurons. J. Neurosci.33 (4), 16721677. 10.1523/JNEUROSCI.3042-12.2013

  • 136

    KimY. H.KimJ. H.KimK.LeeS. H. (2017). Expression of acetylcholinesterase 1 is associated with brood rearing status in the honey bee, Apis mellifera. Sci. Rep.7 (1), 39864. 10.1038/srep39864

  • 137

    KimS.KimK.LeeJ. H.HanS. H.LeeS. H. (2019). Differential expression of acetylcholinesterase 1 in response to various stress factors in honey bee workers. Sci. Rep.9 (1), 10342. 10.1038/s41598-019-46842-0

  • 138

    KlaudinyJ.AlbertS.BachanováK.KopernickýJ.SimúthJ. (2005). Two structurally different defensin genes, one of them encoding a novel defensin isoform, are expressed in honeybee Apis mellifera. Insect Biochem. Mol. Biol.35 (1), 1122. 10.1016/j.ibmb.2004.09.007

  • 139

    KloseR. J.BirdA. P. (2006). Genomic DNA methylation: the mark and its mediators. Trends Biochem. Sci.31 (2), 8997. 10.1016/j.tibs.2005.12.008

  • 140

    KnollS.FaddaV.AhmedF.CappaiM. G. (2024). The nutritional year-cycle of Italian honey bees (Apis mellifera ligustica) in a southern temperate climate. Agriculture14 (5), 730. 10.3390/agriculture14050730

  • 141

    KohliR. M.ZhangY. (2013). TET enzymes, TDG and the dynamics of DNA demethylation. Nature502 (7472), 472479. 10.1038/nature12750

  • 142

    KovacH.StabentheinerA.HetzS. K.PetzM.CrailsheimK. (2007). Respiration of resting honeybees. J. Insect Physiol.53 (12), 12501261. 10.1016/j.jinsphys.2007.06.019

  • 143

    KoyamaT.MendesC. C.MirthC. K. (2013). Mechanisms regulating nutrition-dependent developmental plasticity through organ-specific effects in insects. Front. Physiol.4, 263. 10.3389/fphys.2013.00263

  • 144

    KozakiT.KimmelblattB. A.HammR. L.ScottJ. G. (2008). Comparison of two acetylcholinesterase gene cDNAs of the lesser mealworm,Alphitobius diaperinus, in insecticide susceptible and resistant strains. Archives Insect Biochem. Physiol.67 (3), 130138. 10.1002/arch.20229

  • 145

    KumarR. A.KocA.CernyR. L.GladyshevV. N. (2002). Reaction mechanism, evolutionary analysis, and role of zinc in Drosophila methionine-R-sulfoxide reductase. J. Biol. Chem.277 (40), 3752737535. 10.1074/jbc.M203496200

  • 146

    LaityJ. H.LeeB. M.WrightP. W. (2001). Zinc finger proteins. new insights into structural and functional diversity. Curr. Opin. Struct. Biol.11 (11), 3946. 10.1016/S0959-440X(00)00167-6

  • 147

    LiX. Y.Hooft van HuijsduijnenR.MantovaniR.BenoistC.MathisD. (1992). Intron-exon organization of the NF-Y genes. Tissue-specific splicing modifies an activation domain. J. Biol. Chem.267 (13), 89848990. 10.1016/S0021-9258(19)50377-5

  • 148

    LiZ. W.LiX.YuQ. Y.XiangZ. H.KishinoH.ZhangZ. (2009). The small heat shock protein (sHSP) genes in the silkworm, Bombyx mori, and comparative analysis with other insect sHSP genes. BMC Evol. Biol.9 (1), 215. 10.1186/1471-2148-9-215

  • 149

    LiZ.HuangZ. Y.SharmaD. B.XueY.WangZ.RenB. (2016). Drone and worker brood microclimates are regulated differentially in honey bees, Apis mellifera. PLoS One11 (2), e0148740. 10.1371/journal.pone.0148740

  • 150

    LiG.ZhaoH.ZhangX.ZhangY.ZhaoH.YangX.et al (2018a). Environmental stress responses of DnaJA1, DnaJB12 and DnaJC8 in Apis cerana cerana. Front. Genet.9, 445. 10.3389/fgene.2018.00445

  • 151

    LiG.ZhaoH.WangL.WangY.GuoX.XuB. (2018b). The animal nuclear factor Y: an enigmatic and important heterotrimeric transcription factor. Am. J. Cancer Res.8, 11061125. Available online at: www.ajcr.us/.

  • 152

    LiX.MaW.ShenJ.LongD.FengY.SuW.et al (2019). Tolerance and response of two honeybee species Apis cerana and Apis mellifera to high temperature and relative humidity. PLoS One14 (6), e0217921. 10.1371/journal.pone.0217921

  • 153

    LiG.ZhaoH.GuoH.WangY.CuiX.LiH.et al (2020a). Analyses of the function of DnaJ family proteins reveal an underlying regulatory mechanism of heat tolerance in honeybee. Sci. Total Environ.716, 137036. 10.1016/j.scitotenv.2020.137036

  • 154

    LiG.ZhaoH.GuoH.WangY.CuiX.XuB.et al (2020b). Functional and transcriptomic analyses of the NF-Y family provide insights into the defense mechanisms of honeybees under adverse circumstances. Cell. Mol. Life Sci.77 (23), 49774995. 10.1007/s00018-019-03447-0

  • 155

    LindquistS. (1986). The heat-shock response. Annu. Rev. Biochem.55 (1), 11511191. 10.1146/annurev.bi.55.070186.005443

  • 156

    LiuL.YuX.MengF.GuoX.XuB. (2011). Identification and characterization of a novel corticotropin-releasing hormone-binding protein (CRH-BP) gene from Chinese honeybee (Apis cerana cerana). Archives Insect Biochem. Physiol.78 (3), 161175. 10.1002/arch.20451

  • 157

    LiuT.DanielsC. K.CaoS. (2012). Comprehensive review on the HSC70 functions, interactions with related molecules and involvement in clinical diseases and therapeutic potential. Pharmacol. Ther.136 (3), 354374. 10.1016/j.pharmthera.2012.08.014

  • 158

    LiuP.WangY.DuX.YaoL.LiF.MengZ. (2015). Transcriptome analysis of thermal parthenogenesis of the domesticated silkworm. PLoS One10 (8), e0135215. 10.1371/journal.pone.0135215

  • 159

    LochmillerR. L.DeerenbergC. (2000). Trade‐offs in evolutionary immunology: just what is the cost of immunity?Oikos88 (1), 8798. 10.1034/j.1600-0706.2000.880110.x

  • 160

    LourençoA. P.FloreckiM. M.SimõesZ. L. P.EvansJ. D. (2018). Silencing of Apis mellifera dorsal genes reveals their role in expression of the antimicrobial peptide defensin-1. Insect Mol. Biol.27 (5), 577589. 10.1111/imb.12498

  • 161

    LuA.ZhangQ.ZhangJ.YangB.WuK.XieW.et al (2014). Insect prophenoloxidase: the view beyond immunity. Front. Res. Found.5, 252. 10.3389/fphys.2014.00252

  • 162

    LuS.WeiF.LiG. (2021). The evolution of the concept of stress and the framework of the stress system. Cell Stress26, 7685. 10.15698/CST2021.06.250

  • 163

    LuteijnM. J.KettingR. F. (2013). PIWI-interacting RNAs: from generation to transgenerational epigenetics. Nat. Rev. Genet.14 (8), 523534. 10.1038/nrg3495

  • 164

    MaW.LiX.ShenJ.DuY.XuK.JiangY. (2019). Transcriptomic analysis reveals Apis mellifera adaptations to high temperature and high humidity. Ecotoxicol. Environ. Saf.184, 109599. 10.1016/j.ecoenv.2019.109599

  • 165

    MahalkaA. K.CodeC.JahromiB. R.KirkegaardT.JäätteläM.KinnunenP. K. J. (2011). Activation of phospholipase A2 by Hsp70 in vitro. Biochimica Biophysica Acta (BBA) - Biomembr.1808 (10), 25692572. 10.1016/j.bbamem.2011.06.002

  • 166

    MahatD. B.SalamancaH. H.DuarteF. M.DankoC. G.LisJ. T. (2016). Mammalian heat shock response and mechanisms underlying its genome-wide transcriptional regulation. Mol. Cell62 (1), 6378. 10.1016/j.molcel.2016.02.025

  • 167

    MajzoubJ. A. (2006). Corticotropin-releasing hormone physiology. Eur. J. Endocrinol.155 (1), 7176. 10.1530/eje.1.02247

  • 168

    ManiatisT.TasicB. (2002). Alternative pre-mRNA splicing and proteome expansion in metazoans. Nature418 (6894), 236243. 10.1038/418236a

  • 169

    MantovaniR. (1999). The molecular biology of the CCAAT-binding factor NF-Y. Gene239 (1), 1527. 10.1016/S0378-1119(99)00368-6

  • 170

    MarksP. A.MillerT.RichonV. M. (2003). Histone deacetylases. Curr. Opin. Pharmacol.3 (4), 344351. 10.1016/S1471-4892(03)00084-5

  • 171

    MartelliF.FalconT.PinheiroD. G.SimõesZ. L. P.NunesF. M. F. (2022). Worker bees (Apis mellifera) deprived of pollen in the first week of adulthood exhibit signs of premature aging. Insect Biochem. Mol. Biol.146, 103774. 10.1016/j.ibmb.2022.103774

  • 172

    MartinJ. (1997). Molecular chaperones and mitochondrial protein folding. J. Bioenergetics Biomembr.29 (1), 3543. 10.1023/A:1022407705182

  • 173

    MartinJ.MayhewM.LangerT.HartlF. U. (1993). The reaction cycle of GroEL and GroES in chaperonin-assisted protein folding. Nature366 (6452), 228233. 10.1038/366228a0

  • 174

    MartinsJ. R.NunesF. M. F.CristinoA. S.SimõesZ. L. P.BitondiM. M. G. (2010). The four hexamerin genes in the honey bee: structure, molecular evolution and function deduced from expression patterns in queens, workers and drones. BMC Mol. Biol.11, 23. 10.1186/1471-2199-11-23

  • 175

    MastrangeloA. M.MaroneD.LaidòG.De LeonardisA. M.De VitaP. (2012). Alternative splicing: enhancing ability to cope with stress via transcriptome plasticity. Plant Sci.185–186, 4049. 10.1016/j.plantsci.2011.09.006

  • 176

    MattilaH. R.OtisG. W. (2007). Dwindling pollen resources trigger the transition to broodless populations of long‐lived honeybees each autumn. Ecol. Entomol.32 (5), 496505. 10.1111/j.1365-2311.2007.00904.x

  • 177

    McCawB. A.StevensonT. J.LancasterL. T. (2020). Epigenetic responses to temperature and climate. Integr. Comp. Biol.60 (6), 14691480. 10.1093/icb/icaa049

  • 178

    McKinstryM.ChungC.TruongH.JohnstonB. A.SnowJ. W. (2017). The heat shock response and humoral immune response are mutually antagonistic in honey bees. Sci. Rep.7 (1), 8850. 10.1038/s41598-017-09159-4

  • 179

    McMenaminA. J.DaughenbaughK. F.ParekhF.PizzornoM. C.FlennikenM. L. (2018). Honey bee and bumble bee antiviral defense. Viruses10, 395. 10.3390/v10080395

  • 180

    McMenaminA. J.DaughenbaughK. F.FlennikenM. L. (2020). The heat shock response in thewestern honey bee (Apis mellifera) is antiviral. Viruses12 (2), 245. 10.3390/v12020245

  • 181

    MedinaR. G.PaxtonR. J.De LunaE.Fleites-AyilF. A.Medina MedinaL. A.Quezada-EuánJ. J. G. (2018). Developmental stability, age at onset of foraging and longevity of Africanized honey bees (Apis mellifera L.) under heat stress (Hymenoptera: Apidae). J. Therm. Biol.74, 214225. 10.1016/j.jtherbio.2018.04.003

  • 182

    MedinaR. G.PaxtonR. J.Hernández-SotomayorS. M. T.Pech-JiménezC.Medina-MedinaL. A.Quezada-EuánJ. J. G. (2020). Heat stress during development affects immunocompetence in workers, queens and drones of Africanized honey bees (Apis mellifera L.) (Hymenoptera: Apidae). J. Therm. Biol.89, 102541. 10.1016/j.jtherbio.2020.102541

  • 183

    MemmottJ.CrazeP. G.WaserN. M.PriceM. V. (2007). Global warming and the disruption of plant–pollinator interactions. Ecol. Lett.10 (8), 710717. 10.1111/j.1461-0248.2007.01061.x

  • 184

    MenzelR.HammerM.MüllerU.RosenboomH. (1996). Behavioral, neural and cellular components underlying olfactory learning in the honeybee. J. Physiology-Paris90 (5–6), 395398. 10.1016/S0928-4257(97)87928-4

  • 185

    MercerT. R.DingerM. E.MattickJ. S. (2009). Long non-coding RNAs: insights into functions. Nat. Rev. Genet.10 (3), 155159. 10.1038/nrg2521

  • 186

    MissirlisF.RahlfsS.DimopoulosN.BauerH.BeckerK.HillikerA.et al (2003). A putative glutathione peroxidase of Drosophila encodes a thioredoxin peroxidase that provides resistance against oxidative stress but fails to complement a lack of catalase activity. Biol. Chem.384 (3), 463472. 10.1515/BC.2003.052

  • 187

    MooreL. D.LeT.FanG. (2013). DNA methylation and its basic function. Neuropsychopharmacology38, 2338. 10.1038/npp.2012.112

  • 188

    MorammaziS.ShokrollahiB. (2020). The pattern of HSP70 gene expression, flight activity and temperature in Apis mellifera meda colonies. J. Therm. Biol.91, 102647. 10.1016/j.jtherbio.2020.102647

  • 189

    MorimotoR. I. (2011). The heat shock response: systems biology of proteotoxic stress in aging and disease. Cold Spring Harb. Symposia Quantitative Biol.76 (0), 9199. 10.1101/sqb.2012.76.010637

  • 190

    MoskovitzJ.WeissbachH.BrotN. (1996). Cloning the expression of a mammalian gene involved in the reduction of methionine sulfoxide residues in proteins. Proc. Natl. Acad. Sci. U. S. A.93 (5), 20952099. 10.1073/pnas.93.5.2095

  • 191

    MugliaL.JacobsonL.DikkesP.MajzoubJ. A. (1995). Corticotropin-releasing hormone deficiency reveals major fetal but not adult glucocorticoid need. Nature373 (6513), 427432. 10.1038/373427a0

  • 192

    MukherjeeK.FischerR.VilcinskasA. (2012). Histone acetylation mediates epigenetic regulation of transcriptional reprogramming in insects during metamorphosis, wounding and infection. Front. Zool.9 (1), 25. 10.1186/1742-9994-9-25

  • 193

    MurtazaG.KhanA. K.RashidR.MuneerS.HasanS. M. F.ChenJ. (2017). FOXO transcriptional factors and long-term living. Oxidative Med. Cell. Longev.2017, 34942898. 10.1155/2017/3494289

  • 194

    MustardJ. A.KurshanP. T.HamiltonI. S.BlenauW.MercerA. R. (2005). Developmental expression of a tyramine receptor gene in the brain of the honey bee, Apis mellifera. J. Comp. Neurology483 (1), 6675. 10.1002/cne.20420

  • 195

    MyersZ. A.HoltB. F. (2018). NUCLEAR FACTOR-Y: still complex after all these years?Curr. Opin. Plant Biol.45, 96102. 10.1016/j.pbi.2018.05.015

  • 196

    NabeshimaT.KozakiT.TomitaT.KonoY. (2003). An amino acid substitution on the second acetylcholinesterase in the pirimicarb-resistant strains of the peach potato aphid, Myzus persicae. Biochem. Biophysical Res. Commun.307 (1), 1522. 10.1016/S0006-291X(03)01101-X

  • 197

    NabeshimaT.MoriA.KozakiT.IwataY.HidohO.HaradaS.et al (2004). An amino acid substitution attributable to insecticide-insensitivity of acetylcholinesterase in a Japanese encephalitis vector mosquito, Culex tritaeniorhynchus. Biochem. Biophysical Res. Commun.313 (3), 794801. 10.1016/j.bbrc.2003.11.141

  • 198

    NegriP.MaggiM.RamirezL.SzawarskiN.De FeudisL.LamattinaL.et al (2016). Cellular immunity in Apis mellifera: studying hemocytes brings light about bees skills to confront threats. Apidologie47 (3), 379388. 10.1007/s13592-015-0418-2

  • 199

    NegriP.VillalobosE.SzawarskiN.DamianiN.GendeL.GarridoM.et al (2019). Towards precision nutrition: a novel concept linking phytochemicals. Immune Response and Honey Bee Health Insects10 (11), 401. 10.3390/insects10110401

  • 200

    Neves-da-RochaJ.BitencourtT. A.OliveiraV. M. d.SanchesP. R.RossiA.Martinez-RossiN. M. (2019). Alternative splicing in heat shock protein transcripts as a mechanism of cell adaptation in trichophyton rubrum. Cells8 (10), 1206. 10.3390/cells8101206

  • 201

    ObesoJ. R.HerreraJ. M. (2018). Pollinators and climate change. Ecosistemas27 (2), 5259. 10.7818/ECOS.1371

  • 202

    OhS.-H.KozakiT.MizunoH.TomitaT.KonoY. (2006). Expression of Ace-paralogous acetylcholinesterase of Culex tritaeniorhynchus with an amino acid substitution conferring insecticide insensitivity in baculovirus-insect cell system. Pesticide Biochem. Physiol.85 (1), 4651. 10.1016/j.pestbp.2005.10.003

  • 203

    PaiA. A.HenriquesT.McCueK.BurkholderA.AdelmanK.BurgeC. B. (2017). The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture. eLife6, e32537. 10.7554/eLife.32537

  • 204

    PanH.FinkelT. (2017). Key proteins and pathways that regulate lifespan. J. Biol. Chem.292 (16), 64526460. 10.1074/jbc.R116.771915

  • 205

    PardueM. L.ScottM. P.StortiR. V.LengyelJ. A. (1980). “The heat shock response: a model system for the study of gene regulation in Drosophila,” in Development and neurobiology of Drosophila (Boston, MA: Springer US). 10.1007/978-1-4684-7968-3_5

  • 206

    Di PasqualeG.SalignonM.Le ConteY.BelzuncesL. P.DecourtyeA.KretzschmarA.et al (2013). Influence of pollen nutrition on honey bee health: do pollen quality and diversity matter?PloS One8 (8), e72016. 10.1371/journal.pone.0072016

  • 207

    PerezR.AronS. (2020). Adaptations to thermal stress in social insects: recent advances and future directions. Biol. Rev.95 (6), 15351553. 10.1111/brv.12628

  • 208

    PerryC. J.SøvikE.MyerscoughM. R.BarronA. B. (2015). Rapid behavioral maturation accelerates failure of stressed honey bee colonies. Proc. Natl. Acad. Sci. U. S. A.112 (11), 34273432. 10.1073/pnas.1422089112

  • 209

    PincusD.ChevalierM. W.AragónT.van AnkenE.VidalS. E.El-SamadH.et al (2010). BiP binding to the ER-stress sensor Ire1 tunes the homeostatic behavior of the unfolded protein response. PLoS Biol.8 (7), e1000415. 10.1371/journal.pbio.1000415

  • 210

    PontonF.WilsonK.CotterS. C.RaubenheimerD.SimpsonS. J. (2011). Nutritional immunology: a multi-dimensional approach. PLoS Pathog.7 (12), e1002223. 10.1371/journal.ppat.1002223

  • 211

    PowersE. T.BalchW. E. (2013). Diversity in the origins of proteostasis networks — a driver for protein function in evolution. Nat. Rev. Mol. Cell Biol.14 (4), 237248. 10.1038/nrm3542

  • 212

    PrangeR.WerrM.Löffler-MaryH. (1999). Chaperones involved in hepatitis B virus morphogenesis. Biol. Chem.380 (3), 305314. 10.1515/BC.1999.042

  • 213

    PrattW. B.ToftD. O. (2003). Regulation of signaling protein function and trafficking by the hsp90/hsp70-Based chaperone machinery. Exp. Biol. Med.228 (2), 111133. 10.1177/153537020322800201

  • 214

    PriscoG. D.AnnosciaD.MargiottaM.FerraraR.VarricchioP.ZanniV.et al (2016). A mutualistic symbiosis between a parasitic mite and a pathogenic virus undermines honey bee immunity and health. Proc. Natl. Acad. Sci. U. S. A.113 (12), 32033208. 10.1073/PNAS.1523515113

  • 215

    PubMed (2025). National library of medicine. National Center of Biotechnology Information. Available online at: https://pubmed.ncbi.nlm.nih.gov/ (Accessed June 18, 2025).

  • 216

    QuinlanG. M.GrozingerC. M. (2023). Honey bee nutritional ecology: from physiology to landscapes. Adv. Insect Physiol.64, 289345. 10.1016/bs.aiip.2023.01.003

  • 217

    RadyukS. N.KlichkoV. I.SpinolaB.SohalR. S.OrrW. C. (2001). The peroxiredoxin gene family in drosophila melanogaster. Free Radic. Biol. Med.31 (9), 10901100. 10.1016/S0891-5849(01)00692-X

  • 218

    RamphulU.BoaseT.BassC.OkediL. M.DonnellyM. J.MüllerP. (2009). Insecticide resistance and its association with target-site mutations in natural populations of Anopheles gambiae from eastern Uganda. Trans. R. Soc. Trop. Med. Hyg.103 (11), 11211126. 10.1016/j.trstmh.2009.02.014

  • 219

    RatelD.RavanatJ. L.BergerF.WionD. (2006). N6-methyladenine: the other methylated base of DNA. BioEssays28 (3), 309315. 10.1002/bies.20342

  • 220

    ReimT.BalfanzS.BaumannA.BlenauW.ThammM.ScheinerR. (2017). AmTAR2: functional characterization of a honeybee tyramine receptor stimulating adenylyl cyclase activity. Insect Biochem. Mol. Biol.80, 91100. 10.1016/j.ibmb.2016.12.004

  • 221

    RemnantE. J.AsheA.YoungP. E.BuchmannG.BeekmanM.AllsoppM. H.et al (2016). Parent-of-origin effects on genome-wide DNA methylation in the cape honey bee (Apis mellifera capensis) may be confounded by allele-specific methylation. BMC Genomics17 (1), 226. 10.1186/s12864-016-2506-8

  • 222

    RichardG.JaquiéryJ.Le TrionnaireG. (2021). Contribution of epigenetic mechanisms in the regulation of environmentally-induced polyphenism in insects. Ins. 12 (7), 649. 10.3390/insects12070649

  • 223

    RichterK.HaslbeckM.BuchnerJ. (2010). The heat shock response: life on the verge of death. Mol. Cell40, 253266. 10.1016/j.molcel.2010.10.006

  • 224

    RoederT. (2020). The control of metabolic traits by octopamine and tyramine in invertebrates. J. Exp. Biol.223 (7), jeb194282. 10.1242/jeb.194282

  • 225

    RoederT.SeifertM.KählerC.GeweckeM. (2003). Tyramine and octopamine: antagonistic modulators of behavior and metabolism. Archives Insect Biochem. Physiol.54 (1), 113. 10.1002/arch.10102

  • 226

    RosenzweigR.NillegodaN. B.MayerM. P.BukauB. (2019). The Hsp70 chaperone network. Nat. Rev. Mol. Cell Biol.20, 665680. 10.1038/s41580-019-0133-3

  • 227

    RothS. Y.DenuJ. M.AllisC. D. (2001). Histone acetyltransferases. Annu. Rev. Biochem.70 (1), 81120. 10.1146/annurev.biochem.70.1.81

  • 228

    De la RúaP.JafféR.Dall’OlioR.MuñozI.SerranoJ. (2009). Biodiversity, conservation and current threats to European honeybees. Apidologie40 (3), 263284. 10.1051/apido/2009027

  • 229

    RuanK.ZhuY.BrazillJ. M.ZhaiR. G. (2015). Alternative splicing of Drosophila Nmnat functions as a switch to enhance neuroprotection under stress. Nat. Commun.6 (1), 10057. 10.1038/ncomms10057

  • 230

    RutherfordS. L.LindquistS. (1998). Hsp90 as a capacitor for morphological evolution. Nature396 (6709), 336342. 10.1038/24550

  • 231

    RuttnerF. (1988). Biogeography and taxonomy of honeybees. Berlin, Heidelberg and New York: Springer-Verlag. 10.1007/978-3-642-72649-1

  • 232

    SalmelaH.StarkT.StuckiD.FuchsS.FreitakD.DeyA.et al (2016). Ancient duplications have led to functional divergence of vitellogenin-like genes potentially involved in inflammation and oxidative stress in honey bees. Genome Biol. Evol.8 (3), 495506. 10.1093/gbe/evw014

  • 233

    ScavenV. L.RaffertyN. E. (2013). Physiological effects of climate warming on flowering plants and insect pollinators and potential consequences for their interactions. Curr. Zool.59 (3), 418426. 10.1093/czoolo/59.3.418

  • 234

    ScheinerR.BaumannA.BlenauW. (2006). Aminergic control and modulation of honeybee behaviour. Curr. Neuropharmacol.4 (4), 259276. 10.2174/157015906778520791

  • 235

    SchulzD. J.HuangZ.-Y.RobinsonG. E. (1998). Effects of colony food shortage on behavioral development in honey bees. Behav. Ecol. Sociobiol.42, 295303. 10.1007/s002650050442

  • 236

    Scopus (2025). Scopus: a comprehensive abstract and citation database for impact makers. Available online at: https://www.scopus.com (Accessed June 18, 2025).

  • 237

    SeasholtzA.ValverdeR.DenverR. (2002). Corticotropin-releasing hormone-binding protein: biochemistry and function from fishes to mammals. J. Endocrinol.175 (1), 8997. 10.1677/joe.0.1750089

  • 238

    SeehuusS.-C.NorbergK.GimsaU.KreklingT.AmdamG. V. (2006). Reproductive protein protects functionally sterile honey bee workers from oxidative stress. Proc. Natl. Acad. Sci. U. S. A.103 (4), 962967. 10.1073/pnas.0502681103

  • 239

    SemaniukU.PiskovatskaV.StrilbytskaO.StrutynskaT.BurdyliukN.VaisermanA.et al (2021a). Drosophila insulin-like peptides: from expression to functions – a review. Entomologia Exp. Appl.169 (2), 195208. 10.1111/eea.12981

  • 240

    SemaniukU.StrilbytskaO.MalinovskaK.StoreyK. B.VaisermanA.LushchakV.et al (2021b). Factors that regulate expression patterns of insulin-like peptides and their association with physiological and metabolic traits in Drosophila. Insect Biochem. Mol. Biol.135, 103609. 10.1016/j.ibmb.2021.103609

  • 241

    ShaoY.GaoZ.MarksP. A.JiangX. (2004). Apoptotic and autophagic cell death induced by histone deacetylase inhibitors. Proc. Natl. Acad. Sci. U. S. A.101 (52), 1803018035. 10.1073/pnas.0408345102

  • 242

    ShiY.-Y.ZhengH. J.PanQ. Z.WangZ. L.ZengZ. J. (2015). Differentially expressed microRNAs between queen and worker larvae of the honey bee (Apis mellifera). Apidologie46 (1), 3545. 10.1007/s13592-014-0299-9

  • 243

    ShihS. R.BachD. M.RondeauN. C.SamJ.LovingerN. L.LopatkinA. J.et al (2021). Honey bee sHSP are responsive to diverse proteostatic stresses and potentially promising biomarkers of honey bee stress. Sci. Rep.11 (1), 22087. 10.1038/s41598-021-01547-1

  • 244

    SieberK. R.DormanT.NewellN.YanH. (2021). (Epi)genetic mechanisms underlying the evolutionary success of eusocial insects. Insects12 (6), 498. 10.3390/insects12060498

  • 245

    SilmanI.SussmanJ. L. (2008). Acetylcholinesterase: how is structure related to function?Chemico-Biological Interact.175 (1–3), 310. 10.1016/j.cbi.2008.05.035

  • 246

    SolísE. J.PandeyJ. P.ZhengX.JinD. X.GuptaP. B.AiroldiE. M.et al (2016). Defining the essential function of yeast Hsf1 reveals a compact transcriptional program for maintaining eukaryotic proteostasis. Mol. Cell63 (1), 6071. 10.1016/j.molcel.2016.05.014

  • 247

    SollarsV.LuX.XiaoL.WangX.GarfinkelM. D.RudenD. M. (2003). Evidence for an epigenetic mechanism by which Hsp90 acts as a capacitor for morphological evolution. Nat. Genet.33 (1), 7074. 10.1038/ng1067

  • 248

    SonodaS.AshfaqM.TsumukiH. (2007). A comparison of heat shock protein genes from cultured cells of the cabbage armyworm,Mamestra brassicae, in response to heavy metals. Archives Insect Biochem. Physiol.65 (4), 210222. 10.1002/arch.20178

  • 249

    SørensenJ. G. (2010). Application of heat shock protein expression for detecting natural adaptation and exposure to stress in natural populations. Curr. Zool.56 (6), 703713. 10.1093/czoolo/56.6.703

  • 250

    SouthwickE. E. (1985). Bee hair structure and the effect of hair on metabolism at low temperature. J. Apic. Res.24 (3), 144149. 10.1080/00218839.1985.11100664

  • 251

    SouthwickE. E.HeldmaierG. (1987). Temperature control in honey bee colonies. BioScience37 (6), 395399. 10.2307/1310562

  • 252

    SpannhoffA.KimY. K.RaynalN. J. M.GharibyanV.SuM. B.ZhouY. Y.et al (2011). Histone deacetylase inhibitor activity in royal jelly might facilitate caste switching in bees. EMBO Rep.12 (3), 238243. 10.1038/embor.2011.9

  • 253

    SpecchiaV.PiacentiniL.TrittoP.FantiL.D’AlessandroR.PalumboG.et al (2010). Hsp90 prevents phenotypic variation by suppressing the mutagenic activity of transposons. Nature463 (7281), 662665. 10.1038/nature08739

  • 254

    SporerA. J.KahlL. J.Price-WhelanA.DietrichL. E. P. (2017). Redox-based regulation of bacterial development and behavior. Annu. Rev. Biochem.86 (1), 777797. 10.1146/annurev-biochem-061516-044453

  • 255

    SpotswoodH. T.TurnerB. M. (2002). An increasingly complex code. J. Clin. Investigation110 (5), 577582. 10.1172/JCI16547

  • 256

    SunY.MacRaeT. H. (2005). Small heat shock proteins: molecular structure and chaperone function. Cell. Mol. Life Sci.62 (21), 24602476. 10.1007/s00018-005-5190-4

  • 257

    SuzukiM. M.BirdA. (2008). DNA methylation landscapes: provocative insights from epigenomics. Nat. Rev. Genet.9 (6), 465476. 10.1038/nrg2341

  • 258

    TabuchiY.TakasakiI.DoiT.IshiiY.SakaiH.KondoT. (2006). Genetic networks responsive to sodium butyrate in colonic epithelial cells. FEBS Lett.580 (13), 30353041. 10.1016/j.febslet.2006.04.048

  • 259

    TangA. H.TuC. P. (1994). Biochemical characterization of Drosophila glutathione S-transferases D1 and D21. J. Biol. Chem.269 (45), 2787627884. 10.1016/S0021-9258(18)46868-8

  • 260

    TelferW. H.PanM.LawJ. H. (1991). Lipophorin in developing adults of Hyalophora cecropia: support of yolk formation and preparation for flight. Insect Biochem.21 (6), 653663. 10.1016/0020-1790(91)90036-E

  • 261

    TobaG.AigakiT. (2000). Disruption of the Microsomal glutathione S-transferase-like gene reduces life span of Drosophila melanogaster. Gene253 (2), 179187. 10.1016/S0378-1119(00)00246-8

  • 262

    VoosW.RöttgersK. (2002). Molecular chaperones as essential mediators of mitochondrial biogenesis. Biochim. Biophys. Acta1592, 5162. 10.1016/s0167-4889(02)00264-1

  • 263

    WallbergA.HanF.WellhagenG.DahleB.KawataM.HaddadN.et al (2014). A worldwide survey of genome sequence variation provides insight into the evolutionary history of the honeybee Apis mellifera. Nat. Genet.46 (10), 10811088. 10.1038/ng.3077

  • 264

    WalterP.RonD. (2011). The unfolded protein response: from stress pathway to homeostatic regulation. Science334, 10811086. 10.1126/science.1209038

  • 265

    WangH.-S.WangX. H.ZhouC. S.HuangL. H.ZhangS. F.GuoW.et al (2007). cDNA cloning of heat shock proteins and their expression in the two phases of the migratory locust. Insect Mol. Biol.16 (2), 207219. 10.1111/j.1365-2583.2006.00715.x

  • 266

    WatersE. R.AevermannB. D.Sanders-ReedZ. (2008). Comparative analysis of the small heat shock proteins in three angiosperm genomes identifies new subfamilies and reveals diverse evolutionary patterns. Cell Stress Chaperones13 (2), 127142. 10.1007/s12192-008-0023-7

  • 267

    Web Of Science (2025). Web of science. Available online at: https://www.webofscience.com (Accessed June 18, 2025).

  • 268

    WeddL.MaleszkaR. (2016). DNA methylation and gene regulation in honeybees: from genome-wide analyses to obligatory epialleles. Adv. Exp. Med. Biol.945, 193211. 10.1007/978-3-319-43624-1_9

  • 269

    WeddL.KucharskiR.MaleszkaR. (2016). Differentially methylated obligatory epialleles modulate context-dependent LAM gene expression in the honeybee Apis mellifera. Epigenetics11 (1), 110. 10.1080/15592294.2015.1107695

  • 270

    WeddL.KucharskiR.MaleszkaR. (2022). DNA methylation in honey bees and the unresolved questions in insect methylomics. Adv. Exp. Med. Biol.1389, 159176. 10.1007/978-3-031-11454-0_7

  • 271

    WeillM.MalcolmC.ChandreF.MogensenK.BerthomieuA.MarquineM.et al (2004). The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol. Biol.13 (1), 17. 10.1111/j.1365-2583.2004.00452.x

  • 272

    West-EberhardM. J. (2003). Developmental plasticity and evolution. New York, NY: Oxford University Press. 10.1093/oso/9780195122343.001.0001

  • 273

    WinstonM. L. (1991). The biology of the honey bee. Cambridge (MA): Harvard University Press.

  • 274

    WuC. (1995). Heat shock transcription factors: structure and regulation. Annu. Rev. Cell Dev. Biol.11 (1), 441469. 10.1146/annurev.cb.11.110195.002301

  • 275

    WuS.LiM.TangP. A.FeltonG. W.WangJ. J. (2010). Cloning and characterization of acetylcholinesterase 1 genes from insecticide-resistant field populations of Liposcelis paeta Pearman (Psocoptera: liposcelididae). Insect Biochem. Mol. Biol.40 (5), 415424. 10.1016/j.ibmb.2010.04.001

  • 276

    XuP. J.XiaoJ. H.XiaQ. Y.MurphyB.HuangD. W. (2010). Apis mellifera has two isoforms of cytoplasmic HSP90. Insect Mol. Biol.19 (4), 593597. 10.1111/j.1365-2583.2010.01015.x

  • 277

    YagoundB.SmithN. M. A.BuchmannG.OldroydB. P.RemnantE. J. (2019). Unique DNA methylation profiles are associated with cis-variation in honey bees. Genome Biol. Evol.11 (9), 25172530. 10.1093/gbe/evz177

  • 278

    YamaguchiM.AliM. S.YoshiokaY.LyL. L.YoshidaH. (2017). NF-Y in invertebrates. Biochimica Biophysica Acta (BBA) - Gene Regul. Mech.1860 (5), 630635. 10.1016/j.bbagrm.2016.10.008

  • 279

    ZemachA.McDanielI. E.SilvaP.ZilbermanD. (2010). Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science328 (5980), 916919. 10.1126/science.1186366

  • 280

    ZhangQ.DenlingerD. L. (2010). Molecular characterization of heat shock protein 90, 70 and 70 cognate cDNAs and their expression patterns during thermal stress and pupal diapause in the corn earworm. J. Insect Physiol.56 (2), 138150. 10.1016/j.jinsphys.2009.09.013

  • 281

    ZhangY.ChouS. D.MurshidA.PrinceT. L.SchreinerS.StevensonM. A.et al (2011). The role of heat shock factors in stress-induced transcription. Methods Mol. Biol.787, 2132. 10.1007/978-1-61779-295-3_2

  • 282

    ZhangJ.WangX.VikashV.YeQ.WuD.LiuY.et al (2016). ROS and ROS-mediated cellular signaling. Oxidative Med. Cell. Longev.2016, 4350965. 10.1155/2016/4350965

  • 283

    ZhaoY.SunH.LuJ.LiX.ChenX.TaoD.et al (2005). Lifespan extension and elevated hsp gene expression in Drosophila caused by histone deacetylase inhibitors. J. Exp. Biol.208 (4), 697705. 10.1242/jeb.01439

  • 284

    ZhaoH.LiuH.WeiG.WangH.ZhuY.ZhangR.et al (2021). Comparative life cycle assessment of emergency disposal scenarios for medical waste during the COVID-19 pandemic in China. Apidologie126, 388399. 10.1016/j.wasman.2021.03.034

  • 285

    ZouJ.GuoY.GuettoucheT.SmithD. F.VoellmyR. (1998). Repression of heat shock transcription factor HSF1 activation by HSP90 (HSP90 complex) that forms a stress-sensitive complex with HSF1. Cell94 (4), 471480. 10.1016/S0092-8674(00)81588-3

  • 286

    ZouZ.LopezD. L.KanostM. R.EvansJ. D.JiangH. (2006). Comparative analysis of serine protease-related genes in the honey bee genome: possible involvement in embryonic development and innate immunity. Insect Mol. Biol.15 (5), 603614. 10.1111/j.1365-2583.2006.00684.x

Summary

Keywords

climate change, honey bee, Apis mellifera subspecies, stress, stress-related genes, HSPs, epigenetics, stress response

Citation

Sagastume S, Cilia G, Henriques D, Yadró C, Corona M, Higes M, Pinto MA, Nanetti A and Martín-Hernández R (2025) Climate change-induced stress in the honey bee Apis mellifera L.- a genetic review. Front. Physiol. 16:1623705. doi: 10.3389/fphys.2025.1623705

Received

07 May 2025

Accepted

27 August 2025

Published

17 September 2025

Volume

16 - 2025

Edited by

Makio Takeda, Kobe University, Japan

Reviewed by

Hamzeh Izadi, Vali-E-Asr University of Rafsanjan, Iran

Francis Morais Franco Nunes, Federal University of São Carlos, Brazil

Updates

Copyright

*Correspondence: Soledad Sagastume,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics