Infestation by Phorid Flies Disrupts Behavior and Immune Function in Honey Bees Monitored by Radio-frequency Identification
Gloria Ruiz-Guzmán, Oliverio Delgado-Carrillo, Francisco J. Balvino-Olvera, María de Jesús Aguilar-Aguilar, Violeta Patiño-Conde, Paulo de Souza, Ulises Olivares-Pinto, Mauricio Quesada

TL;DR
Phorid flies disrupt honey bee behavior and weaken their immune systems, potentially threatening colony health.
Contribution
This study is the first to use RFID to monitor phorid fly effects on honey bee behavior and immune function under natural conditions.
Findings
Phorid infestation reduces foraging trips and nocturnal activity in honey bees.
Immune function is disrupted with elevated proPO/PO activity and reduced hemocyte counts.
Physiological immune responses are consistent across adult and pupal stages.
Abstract
Honey bees (Apis mellifera Linnaeus) are important pollinators in natural and agricultural ecosystems and play a key role in global food production. However, colonies are increasingly affected by multiple stressors, including emerging parasites, which can compromise their health and performance. Among these, phorid flies (Diptera: Phoridae) have been observed in association with weak or collapsing honey bee colonies in various regions in Mexico, yet their effects on honey bee behavior and immune function under field conditions remain poorly understood. In this study, we investigated the impact of a phorid fly infestation on the behavior, survival, and immune defense of A. mellifera under natural conditions. Using radio-frequency identification (RFID) technology, we monitored the daily activity of individually tagged bees from a healthy hive and a naturally infested hive. Additionally,…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7- —PAPIIT
- —http://dx.doi.org/10.13039/501100003141Consejo Nacional de Ciencia y Tecnología
- —LANASE
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsInsect and Pesticide Research · Insect and Arachnid Ecology and Behavior · Neurobiology and Insect Physiology Research
Introduction
Over recent decades, bees and other pollinators have faced significant global declines, threatening pollinator-dependent ecosystems and food crops (Ashworth et al. 2009; Novais et al. 2016; Simmons et al. 2019; Hristov et al. 2020; Balvino-Olvera et al. 2024). Currently, managed bees continue to experience colony losses, which remain poorly understood in various parts of the world, including the tropics (vanEngelsdorp et al. 2009; Currie et al. 2010; Brodschneider et al. 2018; Gray et al. 2023; Medina-Flores et al. 2023; Lamas et al. 2024; Requier et al. 2024). These losses have been attributed to multiple factors including habitat degradation, land use change, pesticide exposure, extreme weather events, and parasitism, which are believed to synergistically alter bee performance (vanEngelsdorp et al. 2009; Perry et al. 2015; Potts et al. 2016; Roy et al. 2018).
Honey bees exhibit behaviors that significantly increase their vulnerability to parasites. Specifically, worker bees, the most numerous and active members of a colony, perform a wide range of age-related tasks essential for colony survival, including brood care, nest maintenance, defense, and foraging. Foraging bees, in particular, frequently interact with other colonies (Oliveira et al. 2021) and drone congregation areas (Ayup et al. 2021), thus increasing their exposure to parasites. This not only makes them susceptible to infection but also the likelihood that foragers act as vectors within the colony (Graystock et al. 2015; Nunes-Silva et al. 2019). Consequently, parasitic pressure can alter essential worker functions such as foraging activity, nest-related tasks, and defense-related behavior, and may compromise overall colony performance, particularly under conditions of limited floral resources (Polatto et al. 2014; Koch et al. 2017).
To counter parasites and pathogens, honey bees rely on a suite of behavioral and immune defenses. Behavioral defenses operate at the colony level, including social immunity mechanisms that mitigate pathogen spread and damage (Cremer 2019). In parallel, the immune response involves both cellular (e.g., hemocytes mediating phagocytosis and encapsulation) and humoral components (e.g., antimicrobial peptides, lysozymes, and the prophenoloxidase system) that target a wide range of invaders including viruses, bacteria, fungi, and parasitic arthropods (Gillespie et al. 1997; Kanost et al. 2004; Tang et al. 2006). Among these parasites, phorid flies have been recognized as a significant mortality factor in stingless bee colonies in Central America (Robroek et al. 2003) and are increasingly reported as emerging threats to honey bee health across several regions, including the USA, Italy, Saudi Arabia, Egypt, Algeria, and Central America (Core et al. 2012; Dutto and Ferrazzi 2014; Khattab and El-Hosseny 2014; Menail et al. 2016; Cham et al. 2018; Mohammed 2018; Dias de Freitas et al. 2023). Although most studies focus on the presence of these flies within honey bee colonies (see Table 1), little is known about their actual impact on colony health and defense mechanisms. Table 1. Reviewed studies (from 2010 to 2024) about the parasitic behavior of phorid flies observed in honey beesRegionPhorid fly speciesBee speciesParasitic behaviorEffect of infectionFly originDetection methodParasitism percentageCoinfectionBee defense traitReferenceUSA(California, South Dakota)Apocephalus borealis**Apis melliferaParasitoidBehavioral manipulation, disorientation, nocturnal abandonment, deathField/labMorphological, DNA barcode (COI), rRNAUp to 77% of sitesNosema ceranae, DWV, BQCV/Core et al. 2012Italy(Piedmont)Megaselia rufipes**Apis melliferaParasitoidLimited movementFieldMorphological23% of adultsDWV, Varroa destructor/Dutto and Ferrazzi 2014EgyptApocephalus borealis**Apis melliferaParasitoidAbnormal behavior, abandonmentFieldMorphologicalUp to 37% of colonies//Khattab and El-Hosseny 2014Algeria(Annaba)Megaselia scalaris**Apis mellifera intermissaParasitoid/FieldRT-qPCR/DWV in bees and flies/Menail et al. 2016Italy (central–south)Megaselia scalaris**Apis melliferaParasitoid/FieldMorphologicalUp to 80% of apiaries//Ricchiuti et al. 2016CameroonMegaselia scalaris**Apis melliferaParasitoid/FieldMorphological, DNA barcode, qPCR9.2% of colonies//Cham et al. 2018SlovakiaMegaselia spp.Apis melliferaParasitoid, beehive parasite (on eggs and larvae in capped cells)/FieldMorphological/Fungi/Sabo et al. 2020South KoreaApocephalus borealis**Apis mellifera//FieldRT-qPCR1.71% of adultsAspergillus flavus, ABPV, KBV/Truong et al. 2023The slash indicates missing information in each case
In Mexico, beekeepers have observed phorid flies during episodes of high colony mortality. In 2016, for example, the Comarca Lagunera region in northern Mexico reported the loss of approximately 3066 of 6256 hives. Similar events occurred in central (San Luis Potosí, Querétaro), southern (Oaxaca, Chiapas), and southeastern (Campeche) regions, based on personal reports from beekeepers. While the parasitic behavior of Apocephalus borealis is well documented (Core et al. 2012), other phorid species have also been observed infesting weak or decaying A. mellifera colonies (Disney 2008; Dutto and Ferrazzi 2014), and occasionally even in apparently healthy hives (Ricchiuti et al. 2016). However, no field studies have evaluated how these infestations affect bee behavior and immune function under natural conditions.
In this study, we evaluated the effect of phorid fly infestation on honey bee behavior and immune responses under field conditions. While conducting routine monitoring in ten apiaries in Mexico as part of the Global Initiative for Honey Bee Health (GIHH), led by the Australian Commonwealth Scientific and Industrial Research Organisation, we identified a case of phorid infestation in one hive (H2-infested) adjacent to a healthy reference hive (H1-healthy). We intensified monitoring of these two colonies by equipping individual worker bees with radio-frequency identification (RFID) tags to record their activity patterns (Susanto et al. 2018), capturing detailed records of incoming and outgoing movements. We evaluated the impact of phorid infestation by comparing survival, behavior, and immune response between hives. Immune parameters included humoral (phenoloxidase, prophenoloxidase, and lytic activity) and cellular (hemocyte counts) components, which reflect physiological condition of bees and their investment in defense (Hill 2011). We predicted that the H2-infested hive would show elevated immune activity and reduced survival, despite being managed under the same conditions as the H1-healthy hive. To control for confounding factors, we also screened for prevalent pathogens such as Ascosphaera apis, Melissococcus plutonius, and Vairimorpha (Nosema) ceranae in both colonies. To our knowledge, this is the first field-based study to explore the phorid fly infestation on behavior and defense in A. mellifera under a natural coinfection scenario.
Material and Methods
Study Site, Observations, and Sampling
This study was conducted in July and August 2016 at an apiary located in Los Tigrillos, Michoacán (19°41′58.68″ N, 101°0′2.46″ W) within a pine-oak forest at an elevation of 2023 m.a.s.l. Two hives were intensively monitored: H1 represented the healthy hive not infested by flies (and no other parasites); and H2, the hive colonized by phorid flies. Both hives were located nearby in the same apiary. The H1-healthy hive was monitored for eight consecutive weeks, while the H2-infested hive was monitored for five consecutive weeks, until it collapsed. In the H2-infested hive, phorid eggs, larvae, and pupae were directly collected from naturally parasitized bees at different developmental stages. Adult worker bees and pupae (aged 12–14 days) were sampled from both hives. The brood frames and the collected samples were transported to the laboratory and stored at − 70 °C until use.
Taxonomic and Molecular Classification of Phorid Flies
Phorid flies were classified based on morphological traits described by Disney (1994, 2008) and by molecular barcoding using the cytochrome oxidase I (COI) DNA sequence that distinguishes between Phoridae species (Boehme et al. 2010). The genomic DNA of 21 individual adult flies was extracted using Doyle’s (1991) protocol, modified by adding 1% PVP in the extraction buffer. Samples were incubated with 2.5 μL Proteinase K (20 mg/mL per 500 μL extraction buffer), for 3 h at 65 °C with constant agitation. The COI gene was amplified in a 12.5-μL reaction using the QIAGEN Multiplex PCR Kit with primers LCO1490 and HCO2198 (Folmer et al. 1994), and the PCR conditions adapted from Boehme et al. (2010), including an initial activation step of 15 min at 95 °C for the HotStarTaq DNA polymerase. All PCR reactions were performed on a Mastercycler® nexus gradient (Eppendorf). PCR products were analyzed using the QIAxcel Advanced system-(QIAGEN), where a ~ 740-bp fragment was determined using the Screening Kit (QIAxcel-QIAGEN) and method OM400 described in the QIAxcel DNA Handbook (sample uptake 10 s at 5 kV, separation 400 s at 6 kV). Amplicons were sequenced in both directions by Macrogen (Seoul, Korea), and assembled using Sequencher.
DNA Extraction and Detection of Parasites
The thorax of adult bees and entire pupae were homogenized in liquid nitrogen in pools of five individuals with adults and pupae processed separately. In total, we obtained six pools of adults and seven pools of pupae from the H1-healthy hive, and 13 pools of adults and 12 pools of pupae from the H2-infested hive. The DNA was extracted from each pool using the cetyltrimethylammonium bromide (CTAB) method with 1% PVP and 0.01% β-mercaptoethanol, followed by incubation with 2.5 μL Proteinase K for 1 h at 65 °C (600 rpm). The extracted DNA was quantified using an Eppendorf Biospectrometer® basic. The prevalence of parasites was assessed as the proportion of positive pools, including Vairimorpha ceranae in adults, Melissococcus plutonius, and Ascosphaera apis in pupae. Each pathogen was targeted with a single diagnostic primer (see Table 2 for references and PCR conditions). Table 2. Prevalence of the parasites Ascosphaera apis, Melissococcus plutonius, and Vairimorpha ceranae expressed as the percentage of positive pools in adult worker bees and pupae of the H1-healthy hive and the H2-infested hiveParasite prevalence (% of positive pools)HiveAgeNo. of pools (No. of bees)Ascosphaera apis^a^Melissococcus plutonius^a^Vairimorpha ceranae^b^H1-healthyPupae7 (35)00-Adults6 (30)--100H2-infestedPupae12 (60)00-Adults13 (65)--0
The dashes mean that those parasites were not tested at that developmental stage. Primer references: ^a^Garrido-Bailón et al. (2013); ^b^Martín-Hernández et al. (2007).
Bee Activity
As part of the Global Initiative for Bee Health (GIHH) implemented by the Australian Commonwealth Scientific and Industrial Research Organisation (CSIRO; https://www.csiro.au/en/news/All/News/2015/August/Honey-Bee-Health), the activity of the bees (incoming and outgoing records) was monitored using radio-frequency identification (RFID) tags that were glued to the thorax of adult bees with superglue (see: shellac has produced better results, Toppa et al. 2020). The tag’s size was 2.5 × 2.5 × 0.4 mm; each tag weighed 2.4 mg (Hitachi Chemical, Japan) and allowed the recording of a unique bee identification number in hexadecimal format (de Souza et al. 2018). One hundred tags per hive (H1-healthy and H2-infested) were tagged weekly over eight consecutive weeks. An RFID reader (USB Desktop RU-824) installed at each hive’s entrance recorded the daily tags’ signals, capturing data on each individual bee’s movements (Fig. 1).Fig. 1. Left: RFID reader installed at the entrance of each hive for data collection. Right: A worker honey bee with an RFID tag attached on the top of the thorax. As the bee passes through the only hive entrance, the tag is detected and the bee’s ID, date, and time are recorded in daily CSV files
The RFID system was autonomously controlled by an Edison Intel Atom Processor microcomputer (1 GB memory, 5 GB storage, Wi-Fi, and Bluetooth enabled), powered by solar cells. Data were wireless downloaded weekly. The system recorded each tagged bee’s unique ID upon passing through the reader. Activity was analyzed according to four behavior categories suggested by Susanto et al. (2018): (1) By the entry (BTE) or successive readings of the same bee within less than 3 min indicates that this particular bee is near the RFID reader; (2) short mission (SM) or the individual bee engaged in an activity of short duration (between 3 and 6 min), for example, making orientation or defecation flights and inspecting the hive surroundings; (3) foraging (FG) or the individual bee engaged in activities to search for or exploit food/water sources (6 min to more than 6 h), example bee roles undertaking such activities are scout, exploiter, recruit, and water carrier; and (4) departed bee (DB), the last detection of a bee during its lifetime. Daily activity was categorized according to the day and nighttime it was recorded: (1) sunrise, from 7 a.m. to 12 a.m.; (2) sunset, from 12 a.m. to 7 p.m. (Susanto et al. 2018); (3) night, from 7 p.m. to 12 p.m.; and (4) midnight, from 12 p.m. to 7 a.m. The activity duration was considered as the difference in time between one reading and the subsequent reading.
Immune Response and Health Parameters
A sample of the homogenized tissue of adults and pupae was used to assess immune system parameters including prophenoloxidase activity (proPO), phenoloxidase activity (PO), and lytic activity (LA). Each sample (30 mg) was homogenized in 400 μL of phosphate buffered saline (PBS, pH 7.4, Sigma), vortexed, and centrifuged (17,000 × g for 15 min). The supernatant was used for subsequent protein and immune analyses (Nicoletti et al. 2020). Total protein content was also measured as a proxy of nutritional condition (Contreras-Garduño et al. 2007; Lee et al. 2008; Nicoletti et al. 2020; Burciaga et al. 2023).
Prophenoloxidase (proPO) and Phenoloxidase (PO) Activities
The PO activity in adults and pupae from each hive was quantified by the oxidation of L-dihydroxyphenylalanine (L-DOPA; Sigma-Aldrich) to dopachrome (Contreras-Garduño et al. 2007). The proPO is the inactive zymogen precursor to PO, an important enzyme to innate immune function that results in the formation of melanin (Wilson-Rich et al. 2009), which leads to cuticle sclerotization, wound healing, and cellular defense responses against infectious agents (Ratcliffe et al. 1984; Söderhäll and Cerenius 1998; Chan et al. 2009; Laughton and Siva-Jothy 2011). For each reaction, 50 μL of sample, 40 μL PBS, and 10 μL of L-DOPA solution (4.0 mg/mL) were added to a 96-well plate. Absorbance was read every 3 min for 60 min at 490 nm in a Multiskan GO Microplate reader (Thermo Scientific). In each well of a 96-well microplate, 50 μL of the sample, 40 μL of PBS, and 10 μL of L-DOPA solution were added. In total, 50 μL of PBS was added to the control wells instead of the sample. Each plate was read at 490 nm every 3 min for 60 min in a Multiskan GO Microplate Spectrophotometer (Thermo Scientific) after 15 min of incubation at room temperature. Due to the cytotoxic nature of the products in the melanization cascade, PO is commonly stored in the inactive precursor proPO and activated naturally following the recognition of foreign compounds or artificially by using chymotrypsin that activates PO production (Ratcliffe et al. 1984; Söderhäll and Cerenius 1998; Laughton and Siva-Jothy 2011). Thus, proPO activity was activated using 2 μL of ɑ-chymotrypsin (0.25 mg/mL) as an artificial activator (Ratcliffe et al. 1984; Söderhäll and Cerenius 1998; Eleftherianos et al. 2008; Laughton and Siva-Jothy 2011), and 38 μL of PBS, instead of 40 μL. The PO and proPO activities were expressed as the slope of the absorbance curve during its linear phase (Nicoletti et al. 2020).
Lytic Activity
To quantify lytic activity for each sample, 30 μL of each sample was added to 200 μL of a bacterial suspension (72 mg of lyophilized Micrococcus lysodeikticus in 20 mL PBS). Absorbance at 490 nm was recorded every 3 min for 60 min following a 15-min incubation at room temperature. Activity was calculated as the rate of OD reduction over time (Kortet et al. 2007).
Total Hemocyte Counts
Hemolymph was extracted from anesthetized worker bees by puncturing between abdominal segments III and IV (Chan et al. 2006). Approximately 5 µL per bee was collected in a glass capillary and transferred to a PCR tube containing 5 µL of trypan blue 0.4% in a 1:1 ratio. For total hemocyte counts, samples were transferred to dual chambers of a counting slide. Cell counts and viability were determined using the TC10 automated cell counter (Bio-Rad) following manufacturer protocols.
Protein Content
The total protein content in adults and pupa bees was determined using the Pierce BCA Protein Assay Kit (Thermo Scientific). For each sample, 5 μL was mixed with 45 μL of PBS and 150 μL of BCA reagent (A:B were placed in a radius of 50:1). Absorbance at 562 nm was measured after 15-min incubation at room temperature using a Multiskan GO Microplate reader (Thermo Scientific). The total protein content (μg/5 μL) was estimated using a standard BSA calibration curve.
Statistical Analyses
Statistical analyses were performed in RStudio (version 2023.12.1 + 402). Chi-square tests were used to compare incoming vs. outgoing RFID records, behavioral categories, and daytime-scale frequencies, and T-tests were applied to compare activity duration between hives. Generalized linear models (GLMs) were used to analyze immune parameters. The unit of replication for protein and enzymatic assays (proPO, PO, lytic activity) was a pool of five bees, whereas hemocyte counts were obtained from individual bees. Because only two hives were studied, “Hive” (H1-healthy, H2-infested) was included as a fixed effect, with no random or blocking factors. Data transformations were applied where necessary according to model assumptions, in particular, proPO, PO, and lytic activity data were modeled with a Gamma distribution (link = log), and a constant of “1” was added to manage zeros in the dataset. We initially tested whether including “Age” (adults vs. pupae) or the interaction “Hive × Age” improved model fit, but models with “Hive” alone had the lowest AIC and were retained for parsimony. Parasite prevalence was analyzed by comparing the proportion of positive pools between the hives, and Levene’s test was used to assess variance homogeneity.
Results
Taxonomic and Molecular Classification of Phorid Flies
Adult phorid flies were classified based on morphological traits as outlined by Disney (1994, 2008). The specimens exhibited characteristic morphology of small flies (~ 2.5 mm), including a humpbacked thorax, laterally flattened hind femora, and distinctive wing venation indicative of the genus Megaselia (Fig. 2). Confirmation of their taxonomic classification was achieved through cytochrome oxidase I (COI) gene sequencing, which matched Megaselia sequences in the NCBI BLAST database (Online Resource 1).Fig. 2. Morphological characteristics of Megaselia sp. by SEM micrographs. (A) Entire body view of a third-instar larva. (B) Wing venation highlighting radial veins (R), radial sector (Rs), and subcostal (Sc). (C) Apex of the female terminalia. (D): Propleura with numerous scattered hairs (yellow arrow), compound eyes (e), and antennae (a). (E) Maxillary palpi (p). (F) Mouthparts of Megaselia sp., ventrolateral view of labellum showing large, sharply pointed teeth (yellow arrows)
Observation of Parasitic Behavior in Honey Bee Hives
In the field, Megaselia larvae were observed on both adult bees and brood. Eggs, larvae, pupae, and adult phorids were observed in the brood frames of the H2-infested hive, whereas only larvae were observed on the hive’s stored food (Fig. 3). The densities of Megaselia were estimated at 318 larvae/10 cm^2^ and 291 pupae/10 cm^2^, averaging four pupae per cell. The eggs were rare on the adult bees and were laid mainly on capped and uncapped brood cells. The developmental stages observed suggest that the oviposition occurred approximately 3 weeks prior, while the bees were still alive. No phorid flies were detected in the H1-healthy hive.Fig. 3. The H2-infested hive with Megaselia sp. eggs, larvae, and pupae. (A) Phorid fly larvae feeding on Apis mellifera food. (B) Phorid fly larvae of different ages feeding on Apis mellifera larvae. (C) Phorid fly larvae on capped brood cells. (D) Frame of the H2-infested hive showing the remains of Apis mellifera and Megaselia sp*.* after infestation (yellow arrows indicate some hatched eggs of phorid fly)
Bee Activity
The bees of the H2-infested hive recorded a significantly higher number of incoming and outgoing RFID detections (n = 77,380; 93.1%) compared to the H1-healthy hive (n = 5762; 6.9%) (χ^2^ = 53,344, df = 7, p < 0.0001; Fig. 4, Table 3). Weekly comparisons showed similar numbers of detected bees in week 1 (W1) (χ^2^ = 2.79, df = 1, p = 0.09), W3 (χ^2^ = 3.31, df = 1, p = 0.07), W4 (χ^2^ = 0.29, df = 1, p = 0.59), and W5 (χ^2^ = 0.33, df = 1, p = 0.56). However, during W2, the H1-healthy hive displayed a significantly higher activity level than the H2-infested hive (χ^2^ = 10.13, df = 1, p = 0.001). By W5, the H2-infested hive had collapsed following the onset of Megaselia sp. oviposition in W2, leading to a subsequent surge in larval activity and pupation in W3. Post-collapse, the H1-healthy hive showed increasing bee activity, whereas the H2-infested hive recorded no further activity.Fig. 4. Violin plot with error bar shows the honey bee IDs detected weekly during monitoring. Both H1-healthy and H2-infested hives had a similar number of detected honey bees, except for W2 when the H1-healthy hive (blue circles) had a higher number of active honey bees than the H2-infested hive (orange triangles)Table 3. Total filtered readings, weekly counts of RFID-detected honey bees, and behavioral category frequencies observed in H1-healthy and H2-infested hivesH1-healthyH2-infestedχ^2^dfp valueTotal filtered readings576277,38053,3447 < 0.0001Number of bees detected by week W110192.7910.09 W225710.1310.001 W319323.3110.07 W414170.2910.59 W5570.3310.56 W620020.001 < 0.0001 W726026.001 < 0.0001 W8909.001 < 0.0001Bees’ activity (behavior scale) BTE407176,65765,2651 < 0.0001 SM17329724.91 < 0.0001 FG1433373622.21 < 0.0001 DB85537.410.006Statistical results include Chi-square values, degrees of freedom, and p values
Bee behavior categories also differed between hives (χ^2^ = 16,454, df = 3, p < 0.0001). The H2-infested hive had more short activities “by the entry” (BTE; χ^2^ = 65,265, p < 0.0001) and “short mission” (SM; χ^2^ = 24.9, p < 0.0001) behaviors (76,657 and 297 detections for H2 versus 4071 and 173 detections for H1), while the H1-healthy hive showed more “foraging” (FG; χ^2^ = 622.2, p < 0.0001) and “departed bee” (DB; χ^2^ = 7.4, p = 0.006) behaviors (53 and 373 detections for H2-infested hive versus 85 and 1433 detections for H1-healthy hive, respectively) (Table 3). Importantly, when corrected by foraging effort (FG events), the probability of DB relative to FG was more than twice as high in the H2-infested hive (14.4%) compared with the H1-healthy hive (6.0%), indicating that bees from the infested colony were more likely to fail return from foraging.
There was also a significant difference in behavior duration between each category (F = 7.04, p < 0.0001). BTE and FG behaviors were significantly shorter in the H2-infested hive than in the H1-healthy hive (p < 0.0001, Fig. 5). Regarding the daily activity schedule, it was observed that the H1-healthy hive was active during the “sunrise” and “sunset” daytime categories, having a longer activity duration during the “sunset” (t = − 13.11, p < 0.0001). For hive H2-infested hive, two additional daytime categories were necessary, since it was also recorded during the “night” and “midnight.” At “midnight,” the activity duration was short (averaging 3.97 min) and was only recorded in W4 and W5, while the “night” included longer trips in W3 (t = − 35.54, p < 0.0001) and short trips in W4 (t = 3.73, p < 0.0001) and W5 (Fig. 6).Fig. 5. Behavioral activity over 8 weeks in the H1-healthy and H2-infested hives. The histogram shows the number of behaviors (log of count) and a density plot with each behavior’s duration (log of minutes): BTE, by the entry; DB, departed bee; FG, foraging; and SM, short mission. The H2-infested hive had more BTE and SM behaviors and fewer DB and FG behaviors than the H1-healthy hive. The BTE and FG behaviors had a shorter duration in the H2-infested hive than in the H1-healthy hive. DB behavior is excluded from the density plot as it represents unique detections without a duration. The log scale was required to display numerical data over a very wide range of values compactlyFig. 6Frequency of bee activity at different daytimes for both hives. The H1-healthy hive was active at sunrise and sunset, while the H2-infested hive was active at sunrise, sunset, night, and midnight
Immune Response and Health
Molecular analysis by end-point PCR using specific diagnostic primer showed that adult worker bees and pupa pools from the H2-infested hive had a 0% prevalence of the Ascosphaera apis, Melissococcus plutonius, and Vairimorpha ceranae. In contrast, the H1-healthy showed no detection of A. apis and M. plutonius in pupae, but it presented a 100% prevalence of V. ceranae in adult worker bee pools, indicating that at least one individual per pool was infected (Table 2). Bees from H1-healthy and H2-infested hives had a similar protein content (GLM, χ^2^ = 1.0, p = 3.17), suggesting that the physiological condition of bees at the time of sampling was comparable, regardless of whether the colony was infested with phorid flies.
Immune assays showed significantly higher proPO (GLM, χ^2^ = 13.98, p < 0.0001) and PO activity (GLM, χ^2^ = 22.05, p < 0.0001) in H2-infested hive compared to H1-healthy hive. However, lytic activity was higher in the H1-healthy hive (GLM, χ^2^ = 11.84, p < 0.001; Fig. 7). In all cases, the age of the bees did not significantly affect PO (χ^2^ = 0.008, p = 0.927), proPO (χ^2^ = 0.0001, p = 0.999), or lytic activity (χ^2^ = 1.71, p = 0.19). Total hemocyte counts (cells/mL) were significantly higher in the H1-healthy hive (χ^2^ = 145.22, p < 0.0001), whereas the number of viable hemocytes did not differ between hives (χ^2^ = 0.95, p = 0.33; Fig. 7).Fig. 7. Humoral and cellular immunity of adult honey bees from H1-healthy and H2-infested hives. Protein content did not differ significantly between hives. However, the prophenoloxidase (proPO) and phenoloxidase (PO) activities were significantly higher in bees from the H2-infested hive, whereas lytic activity was significantly higher in bees from the H1-healthy hive. Total hemocyte counts were also higher in the H1-healthy hive, but the number of viable hemocytes did not differ between hives
Discussion
Bees are generally more affected by macro- than microparasites. Boomsma et al. (2005) analyzed interactions between social hosts (bees, wasps, ants, and termites) and ten different parasites and diseases, showing that, after mites, dipterans are the most common. Phorid flies are a highly diverse dipteran group, often associated with social insects. Although their initial interactions are primarily detritivores or commensals (Hölldobler and Wilson 1990), increased population density or reduced food availability can trigger parasitic behavior (Wojcik 1994, Robroek et al. 2003). Several recent reports have documented the presence of phorid flies in A. mellifera colonies (Table 1). However, the physiological and behavioral impacts of these parasites on honey bees remain largely unexplored and are often underestimated. Although Cham et al. (2018) and Ricchiuti et al. (2016) provided preliminary evidence of Megaselia infestations in honey bees, no detailed field studies have evaluated their effects on bee behavior and defense mechanisms. Like other phorid species, Megaselia larvae can feed on bee tissues and hemolymph, potentially reducing the longevity of bees and impairing colony performance (Abou-Shaara and Darwish 2021). Our findings provide the first detailed description of this interaction, suggesting that phorid fly invasion may alter bee activity patterns and compromise colony condition.
The H2-infested hive showed increased short-duration activities (BTE and SM), shorter foraging trips (FG), and fewer departed bee events (DB) compared to the H1-healthy hive, possibly reflecting increased hive cleaning activity and short trips for essential functions such as excretion, resource exploration, and feeding in nearby points (Susanto et al. 2018). In addition, a decrease in honey bees may have altered the caste structure within the H2-infested hive. In this case, foraging bees, which normally make long trips in search of food, may have been forced to remain and defend the colony. A similar phenomenon was found by Mehdiabadi and Gilbert (2002) in the interaction between A. mellifera and the phorid fly Pseudacteon tricuspis, which decreased in half the worker population and reduced foraging rates for the social ant Solenopsis invicta (fire ants). Feener and Brown (1992) similarly observed that Pseudacteon spp*.* reduced foraging rates by up to 50% in the ant Solenopsis geminata. Bragança et al. (1998) also showed that relatively few leafcutter ants of Atta sexdens were attacked by the phorid fly Neodohrniphora sp., but these flies affected the behaviour of the outgoing foraging ants by making them defend the colony instead, this interaction caused a significant reduction in the number of foragers.
The Megaselia sp*.* invasion also appeared to affect the daily activity. H1-healthy hive had activityfrom sunrise to sunset, while the H2-infested hive was active at sunrise, sunset, night, and midnight. A similar behavior has been observed in the ant Pheidole titanis; the foraging workers stop foraging during the day because they defended the colony from the nocturnal oviposition attacks of parasitic phorid flies (Feener 1988). In this case, ants foraged at night when flies did not attacked them and took turns in their activity times to avoid parasitoids (Feener 1988). Core et al. (2012) associated the nocturnal activity of bees infested by phorids with hive abandonment, reporting that bees leaving the hive at night had higher parasitism than those leaving by day. They proposed that A. borealis could have manipulated the bees to the point of changing their circadian rhythms.
In our study, nocturnal activity in the H2-infested hive occurred mainly in weeks 4 and 5, with some evidence as early as week 3, during the period preceding colony collapse. The activity at “midnight” was short and restricted to weeks 4 and 5, while “night” activity included longer trips in week 3 and short trips in weeks 4 and 5. These events mostly consisted of short missions averaging 3.97 min, more consistent with host-driven stress or defensive responses (e.g., cleaning, guarding, agitation) than with sustained nocturnal foraging. This interpretation aligns with reports of nocturnal oviposition in Megaselia (Zulaikha and Zuha 2016), which may have triggered defensive activity, and with the strong hygienic behavior of Africanized bees (Aumeier et al. 2000), which often prioritize removal of infested individuals, eggs, or pupae at night. Such behaviors could explain repeated RFID detections near the reader as workers cleaned or defended the hive. At the same time, the nocturnal activity over several weeks does not rule out parasite-driven disruption of circadian rhythms. True manipulation should be expressed consistently across colonies and confer direct benefit to the parasite (Poulin 1995; Moore 2002), conditions not yet demonstrated for Megaselia. In our dataset, the proportion of “DB” (last detection) events relative to foraging trips “FG” was more than twice as high in H2-infested hive (14.4%) compared with H1-healthy hive (6.0%), suggesting that bees in the H2-infested hive were more likely to fail return, which could reflect stress-related losses or parasite induce effects. Future studies combining experimental infestations, dissections, or molecular screening of night active bees, and hygienic behavior assays will be necessary to determine whether nocturnal activity reflects host defenses or parasite manipulation or both. Alternative explanations such as RFID detection of dead bees removed by workers, or climatic effects, appear unlikely; meteorological data (Servicio Meteorológico Nacional (SMN), Comisión Nacional del Agua (CONAGUA) 2023) indicated stable conditions (9–26 °C) during the monitoring period.
Regarding the number of bees detected daily, both hives recorded similar bee IDs per week, except in W2, when fewer detected IDs were found in the H2-infested hive. This coincided with the estimated oviposition period of phorid flies in this hive. The lower number of detected IDs in W2 may reflect bees remaining inside the hive, as they spent their time cleaning cells and removing dead bodies and eggs, instead of foraging. The H2-infested hive may also have had weaker bees or a lower density, which could have favored the invasion and subsequent collapse.
Total protein concentration was similar between hives. Both hives were managed under identical conditions and had access to the same food sources. Given that detected IDs were similar between the two hives in the W1, for this reason, it is unlikely that the initial colony condition or the population density could predispose H2-infested hive to collapse, nor did it influence the bees’ ability to respond to the invasion or mount a defense against parasitoids. Although phorid flies are sometimes found on bees already dead, dying, or affected by other diseases (Dutto and Ferrazzi 2014), in our study, coinfections did not appear to predispose the H2-infested hive to collapse, with a zero prevalence of three of the most common parasites, M. plutonius, A. apis, and V. ceranae, which causes the European foulbrood, chalkbrood, and nosemosis diseases, respectively.
The strong proPO/PO activation observed in the H2-infested hive is consistent with defenses triggered by parasitoid oviposition, larval activity, and associated tissue damage (Pennacchio and Strand 2006). Although the ability of A. mellifera to melanize dipteran parasitoid eggs remains unclear, A. mellifera possesses a single phenoloxidase (AmPPO, Lourenço et al. 2005) which is activated through the proPO → PO cascade, typically induced against macroparasites such as nematodes or parasitoid eggs and larvae, although it can also act against some microparasites like bacteria or fungi. Our results suggest that enhanced PO activity reflects rapid activation of pre-existing proPO stores, allowing an immediate defense upon infection (Chan et al. 2009). In our assays, proPO and PO activities were measured as ΔAbs_490_/min in slope-based assays over 30 min, rather than Vmax, which typically yields higher absolute values (Wilson-Rich et al. 2008). In our data, proPO:PO ratios were close to 1 in both hives (H1-healthy hive: median 0.7, H2-infested hive: median 1.02), whereas Laughton et al. (2011) reported much lower ratios (~ 1:365 to ~ 1:50) in healthy bees, reflecting the predominance of inactive proPO in circulation. The ratios observed here suggest that a substantial fraction of the proPO pool had been converted into active PO, particularly in the H2-infested hive, consistent with a sustained immune activation. Melanization generally encapsulates large intruders, whereas against microparasites melanin tends to be more diffuse and acts through reactive intermediates (such as quinones and reactive radicals) with antimicrobial activity (Cerenius and Söderhäll 2004). Thus, the specific response of the honey bee immune system to Megaselia spp. oviposition remains an open question for further research.
In the H1-healthy hive, we detected a 100% prevalence of V. ceranae, but this was evaluated as presence/absence by endpoint PCR in pools. A positive pool means that at least one honey bee carried the parasite, but it does not indicate infection intensity. This point is important because parasite load (not measured here) could modulate both immune and behavioral responses. The higher lytic activity and hemocyte counts in the H1-healthy hive may reflect an antimicrobial response against V. ceranae since lysozymes and antimicrobial peptides are activated by Toll and IMD pathways to degrade the cellular layer of bacteria, fungi, and microsporidia (Keeling and Fast 2002; Lemaitre and Hoffmann 2007). Lytic activity is relatively stable under nutritional stress in other insects, even when other immune parameters vary (Adamo et al. 2016), supporting the idea that this pathway provides a robust line of defense. At high parasite loads, sustained lytic activity and cellular defense (granulocytes in larvae or plasmatocytes in adults) would be expected, although V. ceranae can suppress antimicrobial gene expression (suggesting that the microsporidium interferes with activation) and impair longevity and foraging (Antúnez et al. 2009; Higes et al. 2008). At low parasite loads, infections may appear with undetectable symptoms with limited impact on colony dynamics (Forsgren and Fries 2010), and in some cases even upregulate antimicrobial peptides (Lourenço et al. 2021). In fact, infected bees with V. ceranae have been shown to exhibit changes in their cuticular hydrocarbon (CHC) profiles but no significant differences in behavior or social interactions compared with healthy bees (McDonnell et al. 2013). We may speculate that the second scenario appears more consistent with our observations. Interestingly, although total hemocyte counts were significantly higher in the H1-healthy hive, the number of viable hemocytes did not differ between hives (H1, 566 cells/µL vs. H2, 599 cells/µL, see Fig. 7). Reported hemocyte concentrations in healthy honey bees average around 400–600 cells/µL (Wilson-Rich et al. 2008), which is comparable to the viable counts observed in both colonies. This suggests that the excess cells in H1-healthy hive consisted largely of non-viable hemocytes, possibly plasmatocytes undergoing rupture or apoptosis during melanization or coagulation (Terán-Murillo et al. 2025). Similar processes have been documented in Drosophila and Lepidoptera, where melanotic encapsulation involves hemocyte rupture and apoptotic pathways (Eleftherianos et al. 2021). Hemocyte death is recognized as a biologically relevant outcome of immune activation (Terán-Murillo et al. 2025). Hemocyte dynamics are also shaped by intrinsic and extrinsic factors; for example, counts decline with age but not task (Wilson-Rich et al. 2008), and experimental infections with Serratia marcescens reduce hemocyte numbers compared with asymptomatic bees (Burritt et al. 2016). In our study, the elevated but largely non-viable hemocyte count in the H1-healthy hive is consistent with immune activation and turnover under V. ceranae infection, in contrast to the reduced hemocytes in the H2-infested hive under phorid fly pressure.
In addition, trade-offs between lytic activity and proPO/PO pathways have been documented (Cotter et al. 2004; Schmid-Hempel 2005; Chan et al. 2009; Nicoletti et al. 2020). Our results support this trade-off because investment in lytic activity in the H1-healthy hive could have limited proPO/PO activation, whereas the opposite pattern prevailed in the H2-infested hive. However, there is no evidence of a trade-off between lytic activity and hemocyte counts in the H1-healthy hive. Moreover, the lower hemocyte counts in the H2-infested hive could reflect resource depletion or immune exhaustion after sustained activation of the proPO/PO pathway, as hemocytes are required both for both melanization and wound repair. Such patterns may reflect immune system reconfiguration rather than simple trade-offs (Adamo 2017), where some pathways (e.g., lytic activity) are maintained as robust front-line defenses while others (e.g., proPO/PO cascade) are differentially modulated depending on the stressor.
The coinfections in social insects are not yet well understood. Core et al. (2012) found that V. ceranae had a prevalence of 50% in adult phorids and 87.7% in larvae, while the deformed wing virus (DWV) had a prevalence of 25% in adult phorids and 75% in larvae. This suggests that phorids may selectively avoid hosts already harboring infections, possibly by detecting volatiles from infected bees, selecting hives free of other parasites before invading them. Richard et al. (2008) showed that experimentally infected bees with lipopolysaccharides significantly expressed different cuticular hydrocarbon profiles compared to healthy bees, emitting specific volatiles that other parasites might detect. In our case, although both hives were adjacent, it is possible that H1-healthy hive, already infected with V. ceranae, was less attractive, whereas phorids invaded the parasite-free H2-infested hive as a higher-quality resource to develop their offspring.
It is important to note that this study was based on the comparison of a single infested hive and a single healthy hive. Consequently, our results should be regarded as a case study rather than a generalized pattern. Although the behavioral and immunological changes we observed are consistent with beekeeper reports of colony abandonment and collapse following phorid infestations, replicated studies including more colonies and sites will be necessary to confirm these trends and to clarify the mechanisms involved.
Phorid species can be considered an emerging threat that may also affect native bees (Core et al. 2012) that could contribute to colony weakening and finally collapse or abandonment, reinforcing the need for improved management strategies. More attention should be given to a potential increase of infestations by phorid flies in honey bee colonies. For beekeepers, early detection of parasites is crucial for maintaining a healthy and productive apiary, and the infected hives are constantly monitored and treated until collapse or clearance; however, the healthy hives are monitored quickly or less frequently.
Future research should test whether the behavioral and immunological patterns reported here are consistent across multiple colonies and landscapes. A key avenue will be to evaluate the hypothesis that phorids might selectively target uninfected or healthier hives (without or with a low prevalence of microparasites), possibly because they represent higher-quality resources for fly larval development. Therefore, we emphasize the importance of giving special attention to the management of healthy hives and not taking for granted their good condition over time. Confirming or rejecting this idea under replicated and controlled conditions will be essential to clarify the dynamics of phorid infestations and their consequences for honey bee health.
Supplementary Information
Below is the link to the electronic supplementary material.Supplementary file1 (DOCX 2.27 MB)
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Mc Donnell CM, Alaux C, Parrinello H, Desvignes JP, Crauser D, Durbesson E, ... Le Conte Y (2013) Ecto-and endoparasite induce similar chemical and brain neurogenomic responses in the honey bee (Apis mellifera). BMC Ecol 13(1):25. 10.1186/1472-6785-13-2510.1186/1472-6785-13-25PMC 372516223866001 · doi ↗ · pubmed ↗
- 2Robroek BJ, de Jong H, Sommeijer MJ (2003) The behaviour of the kleptoparasite, Pseudohypocera kerteszi (Diptera, Phoridae), in hives of stingless bees (Hymenoptera, Apidae) in Central America. In: Proceedings Of The Section Experimental And Applied Entomology-Netherlands Entomological Society, pp 65–70
- 3Servicio Meteorológico Nacional (SMN), Comisión Nacional del Agua (CONAGUA) (2023) Información estadística climatológica. Available at: https://smn.conagua.gob.mx/es/climatologia/informacion-climatologica/informacion-estadistica-climatologica . Accessed 14 Aug 2023
