Honey bee (Apis mellifera) colony health and pathogen composition in migratory beekeeping operations involved in California almond pollination

Towards a Stable Host-Parasite Relationship Between Honey Bees and Varroa Mites Through Innovative Beekeeping

Varroa destructor is a major factor in declining honey bee health worldwide. Conventional beekeeping involves multiple Varroa treatments, limiting bees' ability to adapt to the mite. To foster a stable host-parasite relationship, we tested an "innovative" beekeeping method with fewer Varroa treatments, focusing on its impact on honey bee health. We compared Varroa mite fall, immune responses, and seasonal dynamics of Deformed Wing Virus-B (DWV-B) in colonies managed under conventional and innovative practices. Viral loads of newly emerged honey bees and foragers were quantified three times during the season. Varroa mite fall was monitored and immune responses were assessed. In spring, bees managed with the innovative method had significantly lower haemocyte counts 48 h after emergence. DWV-B loads did not differ between groups in spring but were higher in summer in bees managed with the innovative method. After summer treatment, DWV-B loads and Varroa mite fall were similar between groups. Despite higher numbers in Varroa mite fall and DWV-B loads in summer, the innovative method reduced both by fall, ensuring healthy winter bee production and colony survival. These findings suggest that reducing Varroa treatments can support a stable host-parasite relationship while minimising negative effects on honey bee health.

Factors influencing the prevalence of acute bee paralysis virus in Apis mellifera and insights into its phylogenetic relationships

Acute bee paralysis virus (ABPV) is a notable pathogen frequently detected in managed honeybee (Apis mellifera) colonies. Although infections are often covert, they can lead to severe outcomes in the presence of Varroa destructor infestations. This study aims to evaluate the prevalence of ABPV and its correlation with a range of biotic and abiotic stressors in managed beehives located in the Central Anatolia and Mediterranean Regions of Türkiye during the spring-summer and autumn seasons of 2021. ABPV was identified in 38.6% of the samples (27/70) using real-time RT-PCR. The high prevalence observed was linked to Varroa destructor infestations, elevated temperatures and dry climatic conditions, migratory beekeeping practices, and disruptions in colony management during the COVID-19 pandemic. Phylogenetic relationships among ABPV strains were elucidated through partial sequencing of the capsid and RNA-dependent RNA polymerase protein coding genes, employing maximum likelihood tree construction with the Tamura 3-parameter model. Turkish ABPV strains clustered into a distinct subclade, sharing 98.4-99% nucleotide identity with European strains, indicative of a monophyletic origin and geographic segregation at the regional or continental level. These findings highlight the necessity for robust surveillance and research programs to monitor ABPV prevalence and mitigate its detrimental effects on colony health and productivity. Additionally, the phylogenetic insights provided by this study enhance our understanding of the geographic distribution and evolutionary dynamics of ABPV strains, offering critical information for future molecular epidemiological research and apicultural management strategies.

Vitellogenins Level as a Biomarker of the Honeybee Colony Strength in Urban and Rural Conditions

The study aimed to verify whether urban beekeeping affects the strength of the honeybee (Apis mellifera) colonies from urban apiaries and the variability of the crucial for their health and long-life protein-vitellogenins. For this purpose, honeybees were kept in two locations-in a city apiary on a roof in the city center and in agricultural areas. Each of the apiaries consisted of six colonies, with the sister queens artificially inseminated with semen from the same pool of drones. The bee colony strength and the variability of the vitellogenins in various tissues in foragers from both apiaries were analyzed from May to August. Here, we revealed that colonies from the urban apiary were more abundant than those from the rural apiary. We observed the compensation mechanism during periods of worker deficiency in the bee colony, which was expressed as a change in the Vgs level in the forager tissues. Using the vitellogenin level as a biomarker of the honeybee colony strength can predict the fate of colonies, especially those with low numbers. The high level of Vgs can be a candidate for bee colony depopulation biomarker.

Analysis of contaminant residues in honey bee hive matrices

Pollinators provide ecological services essential to maintaining our food supply and propagating natural habitats. Populations are in decline due to environmental stressors including pesticides, pathogens, and habitat loss. To better understand the impacts of pesticide exposures on colony health, a field survey in Ohio, USA was conducted to monitor the potential contamination of honey bee colonies by pesticides. Apiaries (n = 10) were situated across an agricultural gradient and samples were collected over a 4-week period encompassing corn planting. Dead bees from entrance traps (DBT), pollen, and in-hive (IH) matrices including bee bread, honey, larvae, and nurse bees were analyzed for a whole suite of pesticides. Out of 210 pesticides targeted, 68 residues were quantified across 306 samples. Neonicotinoids, miticides, and fungicides were the dominant pesticide classes identified throughout all the matrix types. Neonicotinoids were detected at higher concentrations and at higher frequencies compared to fungicides, specifically in field pollen samples. DBT also contained high concentrations of these two contaminant classes, although detection frequencies for neonicotinoids were typically lower. Overall, herbicides and non‑neonicotinoid insecticides were found with low frequency and at low concentrations. For most pesticide classes, trends for the mean concentrations were DBT > IH nurse bees > field pollen > IH larvae > IH honey. Pesticides were detected in 100 % of samples with concentrations ranging from 0.01 ppb (diphenylamine) to 2790 ppb (clothianidin). All samples were contaminated with at least two pesticide residues, while 19 samples presented over ten detects and maximum detections of 20 in DBT. Pesticide residues were positively correlated with agricultural gradients across sites and sampling periods. These findings reveal that foraging leads to the exposure of the entire colony to a wide range of pesticides. Moreover, residues determined in DBT serve as an effective proxy for monitoring hive matrices with significantly less disturbance to active hives.

Spray paint-derived microplastics and incorporated substances as ecotoxicological contaminants in the neotropical bumblebee Bombus atratus

While bumblebees may be exposed to microplastics (MPs), the effects on them are not well studied. Therefore, in this research, we assessed the cytotoxicity of pristine and photodegraded spray paint-derived MPs on the midgut, Malpighian tubules, and hepato-nephrocitic system cells of Bombus atratus workers exposed to 50 mg.L-1 MPs for 96 hours. Histological and histochemical analyses revealed that pristine MPs caused subtle cellular changes, while the exposure to photodegraded MPs led to significant vacuolization, nuclear condensation, and pyknosis. These effects are possibly linked to the release of potentially toxic elements (PTEs) like Copper, Manganese, and Iron from photodegraded MPs, which exceeded Brazil's CONAMA safety limits. Photodegraded MPs also reduced body weight, disrupting homeostasis and potentially decreasing bumblebee's fitness. These findings highlight the importance of studying the toxicity of environmentally realistic MPs, as plastic composition and weathering significantly influence their harmful effects.

Prevalence and Phylogenetic Network Analysis of Nosema apis and Nosema ceranae Isolates from Honeybee Colonies in Türkiye

PURPOSE

Nosemosis is a disease that infects both Western honeybees (Apis mellifera L.) and Asian honeybees (Apis cerana) and causes colony losses and low productivity worldwide. In order to control nosemosis, it is important to determine the distribution and prevalence of this disease agent in a particular region. For this purpose, a national study was conducted to assess the prevalence of Nosema ceranae and N. apis throughout Türkiye, to perform network analyses of the parasites, and to determine the presence of nosemosis.

METHODS

In this study which aimed to assess the prevalence of N. apis and N. ceranae in different colony types and regions where beekeeping is intensive in Türkiye, specimens were collected from hives with no clinical signs.

RESULTS

A total of 1194 Western honeybee colonies in 400 apiaries from 40 provinces of Türkiye were examined by microscopic and molecular techniques. Nosemosis was found in all of 40 provinces. The mean prevalence ratio was 64.3 ± 3.0, with 95% CI in apiaries and 40.5 ± 2.9, 95% CI in hives. Nosema ceranae DNA was detected in all of positive hives, while N. ceranae and N. apis co-infection was detected in only four colonies.

CONCLUSION

This study showed that nosemosis has spread to all provinces, and it is common in every region of Türkiye. All of the N. ceranae or N. apis samples examined were 100% identical within themselves. Network analysis showed that they were within largest haplotype reported worldwide.

Detection and characterization of Deformed Wing Virus (DWV) in apiaries with stationary and migratory management in the province of Entre Ríos, Argentina

In Argentina, migratory activity in search of floral diversity has become a common approach to maximizing honey production. The Entre Ríos province possesses a floral diversity that allows beekeepers to perform migratory or stationary management. Beyond the impact caused by transhumance, migratory colonies in this province start and end the season in monoculture areas. To study the effect of these practices on viral infection, we assayed for the presence, abundance and genetic characterization of the Deformed Wing Virus (DWV) in honey bees from apiaries with both types of management. In migratory apiaries, DWV was detectable in 86.2% of the colonies at the beginning of the season (September 2018), and 66% at the end of the season (March 2019). On the other hand, DWV was detected in 44.11% and 53.12% of stationary samples, at the beginning and the end of the season, respectively. Sequence analysis from migratory and stationary colonies revealed that all samples belonged to DWV-A type. The highest viral loads were detected in migratory samples collected in September. Higher DWV presence and abundance were associated with migratory management and the sampling time. Based on our findings we propose that the benefit of migration to wild flowering areas can be dissipated when the bee colonies end the season with monoculture.

Varroa destructor shapes the unique viral landscapes of the honey bee populations of the Azores archipelago

The worldwide dispersal of the ectoparasitic mite Varroa destructor from its Asian origins has fundamentally transformed the relationship of the honey bee (Apis mellifera) with several of its viruses, via changes in transmission and/or host immunosuppression. The extent to which honey bee-virus relationships change after Varroa invasion is poorly understood for most viruses, in part because there are few places in the world with several geographically close but completely isolated honey bee populations that either have, or have not, been exposed long-term to Varroa, allowing for separate ecological, epidemiological, and adaptive relationships to develop between honey bees and their viruses, in relation to the mite's presence or absence. The Azores is one such place, as it contains islands with and without the mite. Here, we combined qPCR with meta-amplicon deep sequencing to uncover the relationship between Varroa presence, and the prevalence, load, diversity, and phylogeographic structure of eight honey bee viruses screened across the archipelago. Four viruses were not detected on any island (ABPV-Acute bee paralysis virus, KBV-Kashmir bee virus, IAPV-Israeli acute bee paralysis virus, BeeMLV-Bee macula-like virus); one (SBV-Sacbrood virus) was detected only on mite-infested islands; one (CBPV-Chronic bee paralysis virus) occurred on some islands, and two (BQCV-Black queen cell virus, LSV-Lake Sinai virus,) were present on every single island. This multi-virus screening builds upon a parallel survey of Deformed wing virus (DWV) strains that uncovered a remarkably heterogeneous viral landscape featuring Varroa-infested islands dominated by DWV-A and -B, Varroa-free islands naïve to DWV, and a refuge of the rare DWV-C dominating the easternmost Varroa-free islands. While all four detected viruses investigated here were affected by Varroa for one or two parameters (usually prevalence and/or the Richness component of ASV diversity), the strongest effect was observed for the multi-strain LSV. Varroa unambiguously led to elevated prevalence, load, and diversity (Richness and Shannon Index) of LSV, with these results largely shaped by LSV-2, a major LSV strain. Unprecedented insights into the mite-virus relationship were further gained from implementing a phylogeographic approach. In addition to enabling the identification of a novel LSV strain that dominated the unique viral landscape of the easternmost islands, this approach, in combination with the recovered diversity patterns, strongly suggests that Varroa is driving the evolutionary change of LSV in the Azores. This study greatly advances the current understanding of the effect of Varroa on the epidemiology and adaptive evolution of these less-studied viruses, whose relationship with Varroa has thus far been poorly defined.

Seasonal Variability in the Prevalence of DWV Strains in Individual Colonies of European Honeybees in Hawaii

The most prevalent viral pathogen of honeybees is Deformed Wing Virus (DWV) and its two most widely studied and common master-variants are DWV-A and DWV-B. The prevalence of DWV variants in the UK and in the US is changing, with the prevalence of the DWV-A strain declining and DWV-B increasing over time. In 2012, only DWV-A was detected on the Hawaiian Islands of Oahu. In this study we focused on a colony-level survey of DWV strains in a single apiary and examined the prevalence of DWV variants over the course of two years. In 2018 and 2019, a total of 16 colonies underwent viral testing in January, May, and September. Of those 16 colonies, four were monitored in both 2018 and 2019. Individual colonies showed variability of DWV master variants throughout the sampling period. DWV-A was consistently detected; however, the detection of DWV-B was variable across time in individual colonies. Ultimately, this study demonstrated a seasonal variation in both viral prevalence and load for DWV-B, providing a perspective on the dynamic nature of DWV master variants emerging in Hawaii.

Prevalence and genome features of lake sinai virus isolated from Apis mellifera in the Republic of Korea

Lake Sinai Virus (LSV) is an emerging pathogen known to affect the honeybee (Apis mellifera). However, its prevalence and genomic characteristics in the Republic of Korea (ROK) remain unexplored. This study aimed to assess the prevalence of and analyze the LSVs by examining 266 honeybee samples from the ROK. Our findings revealed that LSV exhibited the highest infection rate among the pathogens observed in Korean apiaries, particularly during the reported period of severe winter loss (SWL) in A. mellifera apiaries in 2022. Three LSV genotypes- 2, 3, and 4 -were identified using RNA-dependent RNA polymerase gene analysis. Importantly, the infection rates of LSV2 (65.2%) and LSV3 (73.3%) were significantly higher in colonies experiencing SWL than in those experiencing normal winter loss (NWL) (p < 0.03). Furthermore, this study provides the first near-complete genome sequences of the Korean LSV2, LSV3, and LSV4 strains, comprising 5,759, 6,040, and 5,985 nt, respectively. Phylogenetic analysis based on these near-complete genome sequences demonstrated a close relationship between LSVs in the ROK and China. The high LSV infection rate in colonies experiencing a heightened mortality rate during winter suggests that this pathogen might contribute to SWL in ROK. Moreover, the genomic characteristic information on LSVs in this study holds immense potential for epidemiological information and the selection of specific genes suitable for preventing and treating LSV, including the promising utilization of RNA interference medicine in the future.

The Varroa paradox: infestation levels and hygienic behavior in feral scutellata-hybrid and managed Apis mellifera ligustica honey bees

The Varroa destructor mite is a parasitic threat to managed and feral honey bee colonies around the world. Beekeepers use miticides to eliminate Varroa in commercial hives, but these chemicals can diminish bee health and increase miticide resistance. In contrast, feral honey bees have developed multiple ways to counteract mites without chemical treatment. We compared mite levels, grooming habits, and mite-biting behavior between feral Africanized honey bees (genomically verified Apis mellifera scutellata hybrids) and managed Italian honey bees (A. mellifera ligustica). Surprisingly, there was no difference in mite infestation levels between scutellata-hybrids and managed bees over one year despite the regular use of miticides in managed colonies. We also found no differences in the social immunity responses of the two groups, as measured by their hygienic habits (through worker brood pin-kill assays), self-grooming, and mite-biting behavior. However, we provide the first report that both scutellata-hybrids and managed honey bees bite off mite chemosensory forelegs, which the mites use to locate brood cells for reproduction, to a significantly greater degree than other legs (a twofold greater reduction in foreleg length relative to the most anterior legs). Such biting may impair mite reproduction.

Promising Algerian essential oils as natural acaricides against the honey bee mite Varroa destructor (Acari: Varroidae)

Varroosis induced by Varroa destructor Anderson and Trueman represents the most pathogenic and destructive disease affecting the western honey bee, Apis mellifera. In this study, we investigated the acaricidal activity against the Varroa mite using essential oils (EOs) from the aerial parts of four autochthonous Algerian herbal species, namely Artemisia herba alba, Artemisia campestris, Artemisia judaica and Ruta montana. EOs were obtained by means of hydrodistillation and their composition was characterized by gas chromatography-mass spectrometry. The toxicity of the selected EOs toward V. destructor and A. mellifera adult honey bees was evaluated using the complete exposure method. The results indicate the predominance of davanone (66.9%) in A. herba alba, β-pinene (19.5%) in A. campestris, piperitone (68.7%) in A. judaica and 2-undecanone (70.1%) in R. montana EOs. Interestingly, the LC50 values coupled to bee mortality rates revealed that all tested oils exhibited significant acaricidal efficiency with selectivity ratio (SR) values of 10.77, 8.78, 5.62 and 3.73 for A. campestris, A. judaica, A. herba alba, and R. montana, respectively. These values were better than that of thymol (SR = 3.65), the positive control. These findings suggest that these EOs could be used as plant-derived veterinary acaricides to control varroosis in field conditions.

Comparison of APIStrip passive sampling with conventional sample techniques for the control of acaricide residues in honey bee hives

Western honey bees are very sensitive bioindicators for studying environmental conditions, hence frequently included in many investigations. However, it is very common in both research studies and health surveillance programs to sample different components of the colony, including adult bees, brood and their food reserves. These practices are undoubtedly aggressive for the colony as a whole, and may affect its normal functioning and even compromise its viability. APIStrip-based passive sampling allows long-term monitoring of residues without affecting the colony in any way. In this study, we compared the effectiveness in the control of acaricide residues by using passive and conventional sampling, where the residue levels of the acaricides coumaphos, tau-fluvalinate and amitraz were evaluated. Conventional and APIStrip-based sampling differ in methods for evaluating bee exposure to residues. APIStrip is less invasive than conventional sampling, offers more efficient measurement of environmental contaminants, and can be stored at room temperature, saving costs and minimizing operator error.

Virus Prevalence in Egg Samples Collected from Naturally Selected and Traditionally Managed Honey Bee Colonies across Europe

Monitoring virus infections can be an important selection tool in honey bee breeding. A recent study pointed towards an association between the virus-free status of eggs and an increased virus resistance to deformed wing virus (DWV) at the colony level. In this study, eggs from both naturally surviving and traditionally managed colonies from across Europe were screened for the prevalence of different viruses. Screenings were performed using the phenotyping protocol of the 'suppressed in ovo virus infection' trait but with qPCR instead of end-point PCR and a primer set that covers all DWV genotypes. Of the 213 screened samples, 109 were infected with DWV, 54 were infected with black queen cell virus (BQCV), 3 were infected with the sacbrood virus, and 2 were infected with the acute bee paralyses virus. It was demonstrated that incidences of the vertical transmission of DWV were more frequent in naturally surviving than in traditionally managed colonies, although the virus loads in the eggs remained the same. When comparing virus infections with queen age, older queens showed significantly lower infection loads of DWV in both traditionally managed and naturally surviving colonies, as well as reduced DWV infection frequencies in traditionally managed colonies. We determined that the detection frequencies of DWV and BQCV in honey bee eggs were lower in samples obtained in the spring than in those collected in the summer, indicating that vertical transmission may be lower in spring. Together, these patterns in vertical transmission show that honey bee queens have the potential to reduce the degree of vertical transmission over time.

Isolation, Assessments of Risk Factors, and Antimicrobial Susceptibility Test of Klebsiella from Gut of Bee in and around Haramaya University Bee Farm, East Hararghe, Oromia Regional State, Ethiopia

A cross-sectional study was employed from March 2021 to October 2021 to isolate and identify Klebsiella species found in the gut of honey bees collected from worker of honey bee (Apis mellifera) from hives in Haramaya University bee farm, Damota and Finqile's, managed under traditional and modern beekeeping apiculture. From the selected farm, a total of 60 samples of live adult honey bees were collected purposively. The live adult worker of the honey bee was individually surface-sterilized and complete alimentary canals of the worker bee were dissected and processed for Klebsiella isolation. Descriptive statistics were used to describe the occurrence of Klebsiella species and the proportion of Klebsiella found in the gut was analyzed for the association with study variables by the Pearson chi-square test. The overall prevalence of Klebsiella spp. was 50% from samples. The prevalence of Klebsiella pneumoniae was 26.7% and that of Klebsiella oxytoca was 23.3% from isolated using bacteriological examined samples. The isolates were characterized for the antimicrobial susceptibility test using the disc diffusion method. Among the isolated colonies, Klebsiella pneumoniae had the highest resistance to ampicillin (84.2%) and showed less resistance to gentamycin and trimethoprim sulfamethoxazole (26.3%). Klebsiella oxytoca was highly resistant to ampicillin (54.5%) and erythromycin (54.5%) and showed low and equal resistance to gentamycin and amoxicillin (18.2%). Molecular characterization should be conducted to identify Klebsiella spp. from honey bees. Monitoring antimicrobial effectiveness is recommended to tackle the existing problem in apiculture farms, and its possible public health threat should be noted for community by public health professionals.

Bee foraging preferences, microbiota and pathogens revealed by direct shotgun metagenomics of honey

Honeybees (Apis mellifera) continue to succumb to human and environmental pressures despite their crucial role in providing essential ecosystem services. Owing to their foraging and honey production activities, honeybees form complex relationships with species across all domains, such as plants, viruses, bacteria and other hive pests, making honey a valuable biomonitoring tool for assessing their ecological niche. Thus, the application of honey shotgun metagenomics (SM) has paved the way for a detailed description of the species honeybees interact with. Nevertheless, SM bioinformatics tools and DNA extraction methods rely on resources not necessarily optimized for honey. In this study, we compared five widely used taxonomic classifiers using simulated species communities commonly found in honey. We found that Kraken 2 with a threshold of 0.5 performs best in assessing species distribution. We also optimized a simple NaOH-based honey DNA extraction methodology (Direct-SM), which profiled species seasonal variability similarly to an established column-based DNA extraction approach (SM). Both approaches produce results consistent with melissopalinology analysis describing the botanical landscape surrounding the apiary. Interestingly, we detected a strong stability of the bacteria constituting the core and noncore gut microbiome across seasons, pointing to the potential utility of honey for noninvasive assessment of bee microbiota. Finally, the Direct-SM approach to detect Varroa correlates well with the biomonitoring of mite infestation observed in hives. These observations suggest that Direct-SM methodology has the potential to comprehensively describe honeybee ecological niches and can be tested as a building block for large-scale studies to assess bee health in the field.

Molecular Detection and Differentiation of Arthropod, Fungal, Protozoan, Bacterial and Viral Pathogens of Honeybees

The honeybee Apis mellifera is highly appreciated worldwide because of its products, but also as it is a pollinator of crops and wild plants. The beehive is vulnerable to infections due to arthropods, fungi, protozoa, bacteria and/or viruses that manage to by-pass the individual and social immune mechanisms of bees. Due to the close proximity of bees in the beehive and their foraging habits, infections easily spread within and between beehives. Moreover, international trade of bees has caused the global spread of infections, several of which result in significant losses for apiculture. Only in a few cases can infections be diagnosed with the naked eye, by direct observation of the pathogen in the case of some arthropods, or by pathogen-associated distinctive traits. Development of molecular methods based on the amplification and analysis of one or more genes or genomic segments has brought significant progress to the study of bee pathogens, allowing for: (i) the precise and sensitive identification of the infectious agent; (ii) the analysis of co-infections; (iii) the description of novel species; (iv) associations between geno- and pheno-types and (v) population structure studies. Sequencing of bee pathogen genomes has allowed for the identification of new molecular targets and the development of specific genotypification strategies.

Phenomic analysis of the honey bee pathogen-web and its dynamics on colony productivity, health and social immunity behaviors

Many pathogens and parasites have evolved to overwhelm and suppress their host's immune system. Nevertheless, the interactive effects of these agents on colony productivity and wintering success have been relatively unexplored, particularly in large-scale phenomic studies. As a defense mechanism, honey bees have evolved remarkable social behaviors to defend against pathogen and parasite challenges, which reduce the impact of disease and improve colony health. To investigate the complex role of pathogens, parasites and social immunity behaviors in relation to colony productivity and outcomes, we extensively studied colonies at several locations across Canada for two years. In 2016 and 2017, colonies founded with 1-year-old queens of diverse genetic origin were evaluated, which represented a generalized subset of the Canadian bee population. During each experimental year (May through April), we collected phenotypic data and sampled colonies for pathogen analysis in a standardized manner. Measures included: colony size and productivity (colony weight, cluster size, honey production, and sealed brood population), social immunity traits (hygienic behavior, instantaneous mite population growth rate, and grooming behavior), as well as quantification of gut parasites (Nosema spp., and Lotmaria passim), viruses (DWV-A, DWV-B, BQCV and SBV) and external parasites (Varroa destructor). Our goal was to examine: 1) correlations between pathogens and colony phenotypes; 2) the dynamics of pathogens and parasites on colony phenotypes and productivity traits; and 3) the effects of social immunity behaviors on colony pathogen load. Our results show that colonies expressing high levels of some social immunity behaviors were associated with low levels of pathogens/parasites, including viruses, Nosema spp., and V. destructor. In addition, we determined that elevated viral and Nosema spp. levels were associated with low levels of colony productivity, and that five out of six pathogenic factors measured were negatively associated with colony size and weight in both fall and spring periods. Finally, this study also provides information about the incidence and abundance of pathogens, colony phenotypes, and further disentangles their inter-correlation, so as to better understand drivers of honey bee colony health and productivity.

Impact of Honey Bee Migratory Management on Pathogen Loads and Immune Gene Expression is Affected by Complex Interactions With Environment, Worker Life History, and Season

The effects of honey bee management, such as intensive migratory beekeeping, are part of the ongoing debate concerning causes of colony health problems. Even though comparisons of disease and pathogen loads among differently managed colonies indicate some effects, the direct impact of migratory practices on honey bee pathogens is poorly understood. To test long- and short-term impacts of managed migration on pathogen loads and immunity, experimental honey bee colonies were maintained with or without migratory movement. Individuals that experienced migration as juveniles (e.g., larval and pupal development), as adults, or both were compared to control colonies that remained stationary and therefore did not experience migratory relocation. Samples at different ages and life-history stages (hive bees or foragers), taken at the beginning and end of the active season, were analyzed for pathogen loads and physiological markers of health. Bees exposed to migratory management during adulthood had increased levels of the AKI virus complex (Acute bee paralysis, Kashmir bee, and Israeli acute bee paralysis viruses) and decreased levels of antiviral gene expression (dicer-like). However, those in stationary management as adults had elevated gut parasites (i.e. trypanosomes). Effects of environment during juvenile development were more complex and interacted with life-history stage and season. Age at collection, life-history stage, and season all influenced numerous factors from viral load to immune gene expression. Although the factors that we examined are not independent, the results illuminate potential factors in both migratory and nonmigratory beekeeping that are likely to contribute to colony stress, and also indicate potential mitigation measures.

Global similarity, and some key differences, in the metagenomes of Swedish varroa-surviving and varroa-susceptible honeybees

There is increasing evidence that honeybees (Apis mellifera L.) can adapt naturally to survive Varroa destructor, the primary cause of colony mortality world-wide. Most of the adaptive traits of naturally varroa-surviving honeybees concern varroa reproduction. Here we investigate whether factors in the honeybee metagenome also contribute to this survival. The quantitative and qualitative composition of the bacterial and viral metagenome fluctuated greatly during the active season, but with little overall difference between varroa-surviving and varroa-susceptible colonies. The main exceptions were Bartonella apis and sacbrood virus, particularly during early spring and autumn. Bombella apis was also strongly associated with early and late season, though equally for all colonies. All three affect colony protein management and metabolism. Lake Sinai virus was more abundant in varroa-surviving colonies during the summer. Lake Sinai virus and deformed wing virus also showed a tendency towards seasonal genetic change, but without any distinction between varroa-surviving and varroa-susceptible colonies. Whether the changes in these taxa contribute to survival or reflect demographic differences between the colonies (or both) remains unclear.

Chemical Stimulants and Stressors Impact the Outcome of Virus Infection and Immune Gene Expression in Honey Bees (Apis mellifera)

Western honey bees (Apis mellifera) are ecologically, agriculturally, and economically important plant pollinators. High average annual losses of honey bee colonies in the US have been partially attributed to agrochemical exposure and virus infections. To examine the potential negative synergistic impacts of agrochemical exposure and virus infection, as well as the potential promise of phytochemicals to ameliorate the impact of pathogenic infections on honey bees, we infected bees with a panel of viruses (i.e., Flock House virus, deformed wing virus, or Sindbis virus) and exposed to one of three chemical compounds. Specifically, honey bees were fed sucrose syrup containing: (1) thyme oil, a phytochemical and putative immune stimulant, (2) fumagillin, a beekeeper applied fungicide, or (3) clothianidin, a grower-applied insecticide. We determined that virus abundance was lower in honey bees fed 0.16 ppb thyme oil augmented sucrose syrup, compared to bees fed sucrose syrup alone. Parallel analysis of honey bee gene expression revealed that honey bees fed thyme oil augmented sucrose syrup had higher expression of key RNAi genes (argonaute-2 and dicer-like), antimicrobial peptide expressing genes (abaecin and hymenoptaecin), and vitellogenin, a putative honey bee health and age indicator, compared to bees fed only sucrose syrup. Virus abundance was higher in bees fed fumagillin (25 ppm or 75 ppm) or 1 ppb clothianidin containing sucrose syrup relative to levels in bees fed only sucrose syrup. Whereas, honey bees fed 10 ppb clothianidin had lower virus levels, likely because consuming a near lethal dose of insecticide made them poor hosts for virus infection. The negative impact of fumagillin and clothianidin on honey bee health was indicated by the lower expression of argonaute-2, dicer-like, abaecin, and hymenoptaecin, and vitellogenin. Together, these results indicate that chemical stimulants and stressors impact the outcome of virus infection and immune gene expression in honey bees.

Pollen Treated with a Combination of Agrochemicals Commonly Applied During Almond Bloom Reduces the Emergence Rate and Longevity of Honey Bee (Hymenoptera: Apidae) Queens

Honey bee (Apis mellifera L.) colonies that pollinate California's almond orchards are often exposed to mixtures of agrochemicals. Although agrochemicals applied during almond bloom are typically considered bee-safe when applied alone, their combined effects to honey bees are largely untested. In recent years, beekeepers providing pollination services to California's almond orchards have reported reductions in queen quality during and immediately after bloom, raising concerns that pesticide exposure may be involved. Previous research identified a synergistic effect between the insecticide active ingredient chlorantraniliprole and the fungicide active ingredient propiconazole to lab-reared worker brood, but their effects to developing queens are unknown. To test the individual and combined effects of these pesticides on the survival and emergence of developing queens, we fed worker honey bees in closed queen rearing boxes with pollen artificially contaminated with formulated pesticides containing these active ingredients as well as the spray adjuvant Dyne-Amic, which contains both organosilicone and alkyphenol ethoxylate. The translocation of pesticides from pesticide-treated pollen into the royal jelly secretions of nurse bees was also measured. Despite consistently low levels of all pesticide active ingredients in royal jelly, the survival of queens from pupation to 7 d post-emergence were reduced in queens reared by worker bees fed pollen containing a combination of formulated chlorantraniliprole (Altacor), propiconazole (Tilt), and Dyne-Amic, as well as the toxic standard, diflubenzuron (Dimilin 2L), applied in isolation. These results support recommendations to protect honey bee health by avoiding application of pesticide tank-mixes containing insecticides and adjuvants during almond bloom.

Rapidly quantitative detection of Nosema ceranae in honeybees using ultra-rapid real-time quantitative PCR

BACKGROUND

The microsporidian parasite Nosema ceranae is a global problem in honeybee populations and is known to cause winter mortality. A sensitive and rapid tool for stable quantitative detection is necessary to establish further research related to the diagnosis, prevention, and treatment of this pathogen.

OBJECTIVES

The present study aimed to develop a quantitative method that incorporates ultra-rapid real-time quantitative polymerase chain reaction (UR-qPCR) for the rapid enumeration of N. ceranae in infected bees.

METHODS

A procedure for UR-qPCR detection of N. ceranae was developed, and the advantages of molecular detection were evaluated in comparison with microscopic enumeration.

RESULTS

UR-qPCR was more sensitive than microscopic enumeration for detecting two copies of N. ceranae DNA and 24 spores per bee. Meanwhile, the limit of detection by microscopy was 2.40 × 10⁴ spores/bee, and the stable detection level was ≥ 2.40 × 10⁵ spores/bee. The results of N. ceranae calculations from the infected honeybees and purified spores by UR-qPCR showed that the DNA copy number was approximately 8-fold higher than the spore count. Additionally, honeybees infected with N. ceranae with 2.74 × 10⁴ copies of N. ceranae DNA were incapable of detection by microscopy. The results of quantitative analysis using UR-qPCR were accomplished within 20 min.

CONCLUSIONS

UR-qPCR is expected to be the most rapid molecular method for Nosema detection and has been developed for diagnosing nosemosis at low levels of infection.

Microplastic Polystyrene Ingestion Promotes the Susceptibility of Honeybee to Viral Infection

Microplastics (MPs) are an emerging threat to ecological conservation and biodiversity; however, little is known of the types and possible impacts of MPs in pollinators. To examine whether MPs were present in honeybees, we analyzed the honeybee samples collected in fields from six provinces in China. Four types MPs were identified in honeybee including polystyrene (PS) by Raman spectroscopic analysis, and these plastic polymers were detected in 66.7% bee samples. Then, we assessed the physical and biological impacts of PS of three sizes (0.5, 5, and 50 μm) on bees for 21 days. Next, we measured how the presence of PS affected the Israeli acute paralysis virus proliferation, a small RNA virus associated with bee colony decline. Experimental evidence showed that a large mass of PS was ingested and accumulated within the midgut and enhanced the susceptibility of bees to viral infection. Not only histological analysis showed that PS, especially 0.5 μm PS, damaged the midgut tissue and was subsequently transferred to the hemolymph, trachea, and Malpighian tubules, but also qPCR and transcriptomic results indicated that genes correlated with membrane lipid metabolism, immune response, detoxification, and the respiratory system were significantly regulated after PS ingestion. Our results highlight neglected MP contamination to the bees, a pollination ecosystem stressed by the anthropogenic pollution, and have implications for human health via ingestion of bee products.

Investigating Virus-Host Interactions in Cultured Primary Honey Bee Cells

Honey bee (Apis mellifera) health is impacted by viral infections at the colony, individual bee, and cellular levels. To investigate honey bee antiviral defense mechanisms at the cellular level we further developed the use of cultured primary cells, derived from either larvae or pupae, and demonstrated that these cells could be infected with a panel of viruses, including common honey bee infecting viruses (i.e., sacbrood virus (SBV) and deformed wing virus (DWV)) and an insect model virus, Flock House virus (FHV). Virus abundances were quantified over the course of infection. The production of infectious virions in cultured honey bee pupal cells was demonstrated by determining that naïve cells became infected after the transfer of deformed wing virus or Flock House virus from infected cell cultures. Initial characterization of the honey bee antiviral immune responses at the cellular level indicated that there were virus-specific responses, which included increased expression of bee antiviral protein-1 (GenBank: MF116383) in SBV-infected pupal cells and increased expression of argonaute-2 and dicer-like in FHV-infected hemocytes and pupal cells. Additional studies are required to further elucidate virus-specific honey bee antiviral defense mechanisms. The continued use of cultured primary honey bee cells for studies that involve multiple viruses will address this knowledge gap.

Virus Prevalence and Genetic Diversity Across a Wild Bumblebee Community

Viruses are key population regulators, but we have limited knowledge of the diversity and ecology of viruses. This is even the case in wild host populations that provide ecosystem services, where small fitness effects may have major ecological impacts in aggregate. One such group of hosts are the bumblebees, which have a major role in the pollination of food crops and have suffered population declines and range contractions in recent decades. In this study, we investigate the diversity of four recently discovered bumblebee viruses (Mayfield virus 1, Mayfield virus 2, River Liunaeg virus, and Loch Morlich virus), and two previously known viruses that infect both wild bumblebees and managed honeybees (Acute bee paralysis virus and Slow bee paralysis virus) from isolates in Scotland. We investigate the ecological and environmental factors that determine viral presence and absence. We show that the recently discovered bumblebee viruses were more genetically diverse than the viruses shared with honeybees. Coinfection is potentially important in shaping prevalence: we found a strong positive association between River Liunaeg virus and Loch Morlich virus presence after controlling for host species, location and other relevant ecological variables. We tested for a relationship between environmental variables (temperature, UV radiation, wind speed, and prevalence), but as we had few sampling sites, and thus low power for site-level analyses, we could not conclude anything regarding these variables. We also describe the relationship between the bumblebee communities at our sampling sites. This study represents a first step in the description of predictors of bumblebee infection in the wild.

Characterization of Apis mellifera Gastrointestinal Microbiota and Lactic Acid Bacteria for Honeybee Protection-A Review

Numerous honeybee (Apis mellifera) products, such as honey, propolis, and bee venom, are used in traditional medicine to prevent illness and promote healing. Therefore, this insect has a huge impact on humans’ way of life and the environment. While the population of A. mellifera is large, there is concern that widespread commercialization of beekeeping, combined with environmental pollution and the action of bee pathogens, has caused significant problems for the health of honeybee populations. One of the strategies to preserve the welfare of honeybees is to better understand and protect their natural microbiota. This paper provides a unique overview of the latest research on the features and functioning of A. mellifera. Honeybee microbiome analysis focuses on both the function and numerous factors affecting it. In addition, we present the characteristics of lactic acid bacteria (LAB) as an important part of the gut community and their special beneficial activities for honeybee health. The idea of probiotics for honeybees as a promising tool to improve their health is widely discussed. Knowledge of the natural gut microbiota provides an opportunity to create a broad strategy for honeybee vitality, including the development of modern probiotic preparations to use instead of conventional antibiotics, environmentally friendly biocides, and biological control agents.

Metatranscriptome Analysis of Sympatric Bee Species Identifies Bee Virus Variants and a New Virus, Andrena-Associated Bee Virus-1

Bees are important plant pollinators in agricultural and natural ecosystems. High average annual losses of honey bee (Apis mellifera) colonies in some parts of the world, and regional population declines of some mining bee species (Andrena spp.), are attributed to multiple factors including habitat loss, lack of quality forage, insecticide exposure, and pathogens, including viruses. While research has primarily focused on viruses in honey bees, many of these viruses have a broad host range. It is therefore important to apply a community level approach in studying the epidemiology of bee viruses. We utilized high-throughput sequencing to evaluate viral diversity and viral sharing in sympatric, co-foraging bees in the context of habitat type. Variants of four common viruses (i.e., black queen cell virus, deformed wing virus, Lake Sinai virus 2, and Lake Sinai virus NE) were identified in honey bee and mining bee samples, and the high degree of nucleotide identity in the virus consensus sequences obtained from both taxa indicates virus sharing. We discovered a unique bipartite + ssRNA Tombo-like virus, Andrena-associated bee virus-1 (AnBV-1). AnBV-1 infects mining bees, honey bees, and primary honey bee pupal cells maintained in culture. AnBV-1 prevalence and abundance was greater in mining bees than in honey bees. Statistical modeling that examined the roles of ecological factors, including floral diversity and abundance, indicated that AnBV-1 infection prevalence in honey bees was greater in habitats with low floral diversity and abundance, and that interspecific virus transmission is strongly modulated by the floral community in the habitat. These results suggest that land management strategies that aim to enhance floral diversity and abundance may reduce AnBV-1 spread between co-foraging bees.

Wide diversity of parasites in Bombus terrestris (Linnaeus, 1758) revealed by a high-throughput sequencing approach

Assessing the extent of parasite diversity requires the application of appropriate molecular tools, especially given the growing evidence of multiple parasite co-occurrence. Here, we compared the performance of a next-generation sequencing technology (Ion PGM ™ System) in 12 Bombus terrestris specimens that were PCR-identified as positive for trypanosomatids (Leishmaniinae) in a previous study. These bumblebees were also screened for the occurrence of Nosematidae and Neogregarinorida parasites using both classical protocols (either specific PCR amplification or amplification with broad-range primers plus Sanger sequencing) and Ion PGM sequencing. The latter revealed higher parasite diversity within individuals, especially among Leishmaniinae (which were present as a combination of Lotmaria passim, Crithidia mellificae and Crithidia bombi), and the occurrence of taxa never reported in these hosts: Crithidia acanthocephali and a novel neogregarinorida species. Furthermore, the complementary results produced by the different sets of primers highlighted the convenience of using multiple markers to minimize the chance of some target organisms going unnoticed. Altogether, the deep sequencing methodology offered a more comprehensive way to investigate parasite diversity than the usual identification methods and provided new insights whose importance for bumblebee health should be further analysed.

Colony-Level Effects of Amygdalin on Honeybees and Their Microbes

Amygdalin, a cyanogenic glycoside, is found in the nectar and pollen of almond trees, as well as in a variety of other crops, such as cherries, nectarines, apples and others. It is inevitable that western honeybees (Apis mellifera) consistently consume amygdalin during almond pollination season because almond crops are almost exclusively pollinated by honeybees. This study tests the effects of a field-relevant concentration of amygdalin on honeybee microbes and the activities of key honeybee genes. We executed a two-month field trial providing sucrose solutions with or without amygdalin ad libitum to free-flying honeybee colonies. We collected adult worker bees at four time points and used RNA sequencing technology and our HoloBee database to assess global changes in microbes and honeybee transcripts. Our hypothesis was that amygdalin will negatively affect bee microbes and possibly immune gene regulation. Using a log2 fold-change cutoff at two and intraday comparisons, we show no large change of bacterial counts, fungal counts or key bee immune gene transcripts, due to amygdalin treatment in relation to the control. However, relatively large titer decreases in the amygdalin treatment relative to the control were found for several viruses. Chronic bee paralysis virus levels had a sharp decrease (-14.4) with titers then remaining less than the control, Black queen cell virus titers were lower at three time points (<-2) and Deformed wing virus titers were lower at two time points (<-6) in amygdalin-fed compared to sucrose-fed colonies. Titers of Lotmaria passim were lower in the treatment group at three of the four dates (<-4). In contrast, Sacbrood virus had two dates with relative increases in its titers (>2). Overall, viral titers appeared to fluctuate more so than bacteria, as observed by highly inconstant patterns between treatment and control and throughout the season. Our results suggest that amygdalin consumption may reduce several honeybee viruses without affecting other microbes or colony-level expression of immune genes.

Factors Associated with Honey Bee Colony Losses: A Mini-Review

The Western honey bee (Apis mellifera L., Hymenoptera: Apidae) is a species of crucial economic, agricultural and environmental importance. In the last ten years, some regions of the world have suffered from a significant reduction of honey bee colonies. In fact, honey bee losses are not an unusual phenomenon, but in many countries worldwide there has been a notable decrease in honey bee colonies. The cases in the USA, in many European countries, and in the Middle East have received considerable attention, mostly due to the absence of an easily identifiable cause. It has been difficult to determine the main factors leading to colony losses because of honey bees’ diverse social behavior. Moreover, in their daily routine, they make contact with many agents of the environment and are exposed to a plethora of human activities and their consequences. Nevertheless, various factors have been considered to be contributing to honey bee losses, and recent investigations have established some of the most important ones, in particular, pests and diseases, bee management, including bee keeping practices and breeding, the change in climatic conditions, agricultural practices, and the use of pesticides. The global picture highlights the ectoparasitic mite Varroa destructor as a major factor in colony loss. Last but not least, microsporidian parasites, mainly Nosema ceranae, also contribute to the problem. Thus, it is obvious that there are many factors affecting honey bee colony losses globally. Increased monitoring and scientific research should throw new light on the factors involved in recent honey bee colony losses. The present review focuses on the main factors which have been found to have an impact on the increase in honey bee colony losses.

Longitudinal monitoring of honey bee colonies reveals dynamic nature of virus abundance and indicates a negative impact of Lake Sinai virus 2 on colony health

Honey bees (Apis mellifera) are important pollinators of plants, including those that produce nut, fruit, and vegetable crops. Therefore, high annual losses of managed honey bee colonies in the United States and many other countries threaten global agriculture. Honey bee colony deaths have been associated with multiple abiotic and biotic factors, including pathogens, but the impact of virus infections on honey bee colony population size and survival are not well understood. To further investigate seasonal patterns of pathogen presence and abundance and the impact of viruses on honey bee colony health, commercially managed colonies involved in the 2016 California almond pollination event were monitored for one year. At each sample date, colony health and pathogen burden were assessed. Data from this 50-colony cohort study illustrate the dynamic nature of honey bee colony health and the temporal patterns of virus infection. Black queen cell virus, deformed wing virus, sacbrood virus, and the Lake Sinai viruses were the most readily detected viruses in honey bee samples obtained throughout the year. Analyses of virus prevalence and abundance revealed pathogen-specific trends including the overall increase in deformed wing virus abundance from summer to fall, while the levels of Lake Sinai virus 2 (LSV2) decreased over the same time period. Though virus prevalence and abundance varied in individual colonies, analyses of the overall trends reveal correlation with sample date. Total virus abundance increased from November 2015 (post-honey harvest) to the end of the almond pollination event in March 2016, which coincides with spring increase in colony population size. Peak total virus abundance occurred in late fall (August and October 2016), which correlated with the time period when the majority of colonies died. Honey bee colonies with larger populations harbored less LSV2 than weaker colonies with smaller populations, suggesting an inverse relationship between colony health and LSV2 abundance. Together, data from this and other longitudinal studies at the colony level are forming a better understanding of the impact of viruses on honey bee colony losses.

An updated genetic marker for detection of Lake Sinai Virus and metagenetic applications

Background

Lake Sinai Viruses (LSV) are common RNA viruses of honey bees (Apis mellifera) that frequently reach high abundance but are not linked to overt disease. LSVs are genetically heterogeneous and collectively widespread, but despite frequent detection in surveys, the ecological and geographic factors structuring their distribution in A. mellifera are not understood. Even less is known about their distribution in other species. Better understanding of LSV prevalence and ecology have been hampered by high sequence diversity within the LSV clade.

Methods

Here we report a new polymerase chain reaction (PCR) assay that is compatible with currently known lineages with minimal primer degeneracy, producing an expected 365 bp amplicon suitable for end-point PCR and metagenetic sequencing. Using the Illumina MiSeq platform, we performed pilot metagenetic assessments of three sample sets, each representing a distinct variable that might structure LSV diversity (geography, tissue, and species).

Results

The first sample set in our pilot assessment compared cDNA pools from managed A. mellifera hives in California (n = 8) and Maryland (n = 6) that had previously been evaluated for LSV2, confirming that the primers co-amplify divergent lineages in real-world samples. The second sample set included cDNA pools derived from different tissues (thorax vs. abdomen, n = 24 paired samples), collected from managed A. mellifera hives in North Dakota. End-point detection of LSV frequently differed between the two tissue types; LSV metagenetic composition was similar in one pair of sequenced samples but divergent in a second pair. Overall, LSV1 and intermediate lineages were common in these samples whereas variants clustering with LSV2 were rare. The third sample set included cDNA from individual pollinator specimens collected from diverse landscapes in the vicinity of Lincoln, Nebraska. We detected LSV in the bee Halictus ligatus (four of 63 specimens tested, 6.3%) at a similar rate as A. mellifera (nine of 115 specimens, 7.8%), but only one H. ligatus sequencing library yielded sufficient data for compositional analysis. Sequenced samples often contained multiple divergent LSV lineages, including individual specimens. While these studies were exploratory rather than statistically powerful tests of hypotheses, they illustrate the utility of high-throughput sequencing for understanding LSV transmission within and among species.

Longitudinal analysis on parasite diversity in honeybee colonies: new taxa, high frequency of mixed infections and seasonal patterns of variation

To evaluate the influence that parasites have on the losses of Apis mellifera it is essential to monitor their presence in the colonies over time. Here we analysed the occurrence of nosematids, trypanosomatids and neogregarines in five homogeneous colonies for up to 21 months until they collapsed. The study, which combined the use of several molecular markers with the application of a massive parallel sequencing technology, provided valuable insights into the epidemiology of these parasites: (I) it enabled the detection of parasite species rarely reported in honeybees (Nosema thomsoni, Crithidia bombi, Crithidia acanthocephali) and the identification of two novel taxa; (II) it revealed the existence of a high rate of co-infections (80% of the samples harboured more than one parasite species); (III) it uncovered an identical pattern of seasonal variation for nosematids and trypanosomatids, that was different from that of neogregarines; (IV) it showed that there were no significant differences in the fraction of positive samples, nor in the levels of species diversity, between interior and exterior bees; and (V) it unveiled that the variation in the number of parasite species was not directly linked with the failure of the colonies.

The Pathogen Profile of a Honey Bee Queen Does Not Reflect That of Her Workers

Throughout a honey bee queen's lifetime, she is tended to by her worker daughters, who feed and groom her. Such interactions provide possible horizontal transmission routes for pathogens from the workers to the queen, and as such a queen's pathogen profile may be representative of the workers within a colony. To explore this further, we investigated known honey bee pathogen co-occurrence, as well as pathogen transmission from workers to queens. Queens from 42 colonies were removed from their source hives and exchanged into a second, unrelated foster colony. Worker samples were taken from the source colony on the day of queen exchange and the queens were collected 24 days after introduction. All samples were screened for Nosema spp., Trypanosome spp., acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), Israeli acute paralysis virus (IAPV), Lake Sinai virus (LSV), and deformed wing virus master variants (DWV-A, B, and C) using RT-qPCR. The data show that LSV, Nosema, and DWV-B were the most abundant pathogens in colonies. All workers (n = 42) were LSV-positive, 88% were Nosema-positive, whilst pathogen loads were low (<1 × 106 genome equivalents per pooled worker sample). All queens (n = 39) were negative for both LSV and Nosema. We found no evidence of DWV transmission occurring from worker to queen when comparing queens to foster colonies, despite DWV being present in both queens and workers. Honey bee pathogen presence and diversity in queens cannot be revealed from screening workers, nor were pathogens successfully transmitted to the queen.

Comparative analysis of viruses in four bee species collected from agricultural, urban, and natural landscapes

Managed honey bees (Apis mellifera L.) and wild bees provide critical ecological services that shape and sustain natural, agricultural, and urban landscapes. In recent years, declines in bee populations have highlighted the importance of the pollination services they provide and the need for more research into the reasons for global bee losses. Several stressors cause declining populations of managed and wild bee species such as habitat degradation, pesticide exposure, and pathogens. Viruses, which have been implicated as a key stressor, are able to infect a wide range of species and can be transmitted both intra- and inter-specifically from infected bee species to uninfected bee species via vertical (from parent to offspring) and/or horizontal (between individuals via direct or indirect contact) transmission. To explore how viruses spread both intra- and inter-specifically within a community, we examined the impact of management, landscape type, and bee species on the transmission of four common viruses in Nebraska: Deformed wing virus (DWV), Israeli acute paralysis virus (IAPV), Black queen cell virus (BQCV), and Sacbrood virus (SBV). Results indicated the prevalence of viruses is significantly affected (P < 0.005) by bee species, virus type, and season, but not by landscape or year (P = 0.290 and 0.065 respectively). The higher prevalence of DWV detected across bee species (10.4% on Apis mellifera, 5.3% on Bombus impatiens, 6.1% on Bombus griseocollis, and 22.44% on Halictus ligatus) and seasons (10.8% in early-mid summer and 11.4% in late summer) may indicate a higher risk of interspecific transmission of DWV. However, IAPV was predominately detected in Halictus ligatus (20.7%) and in late season collections (28.1%), which may suggest species-specific susceptibility and seasonal trends in infection rates associated with different virus types. However, there were limited detections of SBV and BQCV in bees collected during both sampling periods, indicating SBV and BQCV may be less prevalent among bee communities in this area.

First Complete Genome of Lake Sinai Virus Lineage 3 and Genetic Diversity of Lake Sinai Virus Strains From Honey Bees and Bumble Bees

The complete genome of Lake Sinai virus 3 (LSV3) was sequenced by the Ion Torrent next-generation sequencing (NGS) technology from an archive sample of honey bees collected in 2010. This strain M92/2010 is the first complete genome sequence of LSV lineage 3. From October 2016 to December 2017, 56 honey bee samples from 32 different locations and 41 bumble bee samples from five different locations were collected. These samples were tested using a specific reverse transcriptase-polymerase chain reaction (RT-PCR) method; 75.92% of honey bee samples and 17.07% of bumble bee samples were LSV-positive with the RT-PCR method. Phylogenetic comparison of 557-base pair-long RNA-dependent RNA polymerase (RdRp) genome region of selected 23 positive samples of honey bees and three positive bumble bee samples identified three different LSV lineages: LSV1, LSV2, and LSV3. The LSV3 lineage was confirmed for the first time in Slovenia in 2010, and the same strain was later detected in several locations within the country. The LSV strains detected in bumble bees are from 98.6 to 99.4% identical to LSV strains detected among honey bees in the same territory.

The impact of lowbush blueberry (Vaccinium angustifolium Ait.) and cranberry (Vaccinium macrocarpon Ait.) pollination on honey bee (Apis mellifera L.) colony health status

Commercial lowbush blueberry (Vaccinium angustifolium Ait.) and cranberry (Vaccinium macrocarpon Ait.) crops benefit from the presence of honey bee (Apis mellifera L.) for pollination. Unfortunately, beekeepers are observing negative impacts of pollination services on honey bee colonies. In this study, we investigated three beekeeping management strategies (MS) and measured their impact on honey bee colony health and development. Experimental groups (five colonies/MS) were: A) Control farmland honey producing MS (control MS); B) Blueberry pollination MS (blueberry MS); C) Cranberry pollination MS (cranberry MS) and D) Double pollination MS, blueberry followed by cranberry (double MS). Our goals were to 1) compare floral abundance and attractiveness of foraging areas to honey bees between apiaries using a Geographic Information System, and 2) compare honey bee colony health status and population development between MS during a complete beekeeping season. Our results show significantly lower floral abundance and honey bee attractiveness of foraging areas during cranberry pollination compared to the other environments. The blueberry pollination site seemed to significantly reduce brood population in the colonies who provided those services (blueberry MS and double MS). The cranberry pollination site seemed to significantly reduce colony weight gain (cranberry MS and double MS) and induce a significantly higher winter mortality rate (cranberry MS). We also measured significantly higher levels of Black queen cell virus and Sacbrood virus in the MS providing cranberry pollination (cranberry MS and double MS).

Defining Pollinator Health: A Holistic Approach Based on Ecological, Genetic, and Physiological Factors

Evidence for global bee population declines has catalyzed a rapidly evolving area of research that aims to identify the causal factors and to effectively assess the status of pollinator populations. The term pollinator health emerged through efforts to understand causes of bee decline and colony losses, but it lacks a formal definition. In this review, we propose a definition for pollinator health and synthesize the available literature on the application of standardized biomarkers to assess health at the individual, colony, and population levels. We focus on biomarkers in honey bees, a model species, but extrapolate the potential application of these approaches to monitor the health status of wild bee populations. Biomarker-guided health measures can inform beekeeper management decisions, wild bee conservation efforts, and environmental policies. We conclude by addressing challenges to pollinator health from a One Health perspective that emphasizes the interplay between environmental quality and human, animal, and bee health.

DWV-A Lethal to Honey Bees (Apis mellifera): A Colony Level Survey of DWV Variants (A, B, and C) in England, Wales, and 32 States across the US

The strong association between Varroa destructor, deformed wing virus (DWV), and high overwintering colony losses (OCL) of honey bees is well established. Three DWV master variants (DWV-A, -B, and -C) have been described, and their role in colony mortality remains an open question. Therefore, the aim of this study is to investigate the seasonal prevalence, viral load, and changing distribution of the three DWV master variants within honey bee colonies from England, Wales, and 32 states across the United States. Here, we report that in 2016, DWV-B was prevalent (100%, n = 249) and dominant (95%) in England and Wales, compared to the US. (56%, n = 217 and 23%, respectively), where DWV-A was prevalent (83%, n = 217) and dominant (63%). DWV-C was regularly detected in low viral loads (<1 × 10⁷ genome equivalents per bee) and at lower prevalence (58% in England and Wales, n = 203, and 14% across the United States, n = 124) compared to DWV-A and -B. DWV-B prevalence and dominance in England and Wales coincided with low OCL (6%). Meanwhile, a 60% loss was reported by participating U.S. beekeepers. In the United States, DWV-A prevalence (89%, n = 18) and viral load were significantly (p = 0.002) higher (1 × 10 ⁸–1 × 10¹¹) in colonies that died when compared to the surviving colonies (49% (n = 27), 1 × 10⁶–1 × 10¹⁰). DWV-B had low prevalence (56%, n = 18) in the colonies that died with viral loads of <1 × 10¹⁰. However, DWV-B was routinely detected in high viral loads (>1 × 10¹⁰) in surviving colonies from all sample locations, providing further supporting evidence of DWV-A exhibiting increased virulence over DWV-B at the colony level.

Relative abundance and molecular evolution of Lake Sinai Virus (Sinaivirus) clades

Lake Sinai Viruses (Sinaivirus) are commonly detected in honey bees (Apis mellifera) but no disease phenotypes or fitness consequences have yet been demonstrated. This viral group is genetically diverse, lacks obvious geographic structure, and multiple lineages can co-infect individual bees. While phylogenetic analyses have been performed, the molecular evolution of LSV has not been studied extensively. Here, I use LSV isolates from GenBank as well as contigs assembled from honey bee Sequence Read Archive (SRA) accessions to better understand the evolutionary history of these viruses. For each ORF, substitution rate variation, codon usage, and tests of positive selection were evaluated. Outlier regions of high or low diversity were sought with sliding window analysis and the role of recombination in creating LSV diversity was explored. Phylogenetic analysis consistently identified two large clusters of sequences that correspond to the current LSV1 and LSV2 nomenclature, however lineages sister to LSV1 were the most frequently detected in honey bee SRA accessions. Different expression levels among ORFs suggested the occurrence of subgenomic transcripts. ORF1 and RNA-dependent RNA polymerase had higher evolutionary rates than the capsid and ORF4. A hypervariable region of the ORF1 protein-coding sequence was identified that had reduced selective constraint, but a site-based model of positive selection was not significantly more likely than a neutral model for any ORF. The only significant recombination signals detected between LSV1 and LSV2 initiated within this hypervariable region, but assumptions of the test (single-frame coding and independence of substitution rate by site) were violated. LSV codon usage differed strikingly from that of honey bees and other common honey-bee viruses, suggesting LSV is not strongly co-evolved with that host. LSV codon usage was significantly correlated with that of Varroa destructor, however, despite the relatively weak codon bias exhibited by the latter. While codon usage between the LSV1 and LSV2 clusters was similar for three ORFs, ORF4 codon usage was uncorrelated between these clades, implying rapid divergence of codon use for this ORF only. Phylogenetic placement and relative abundance of LSV isolates reconstructed from SRA accessions suggest that detection biases may be over-representing LSV1 and LSV2 in public databases relative to their sister lineages.

Putative Drone Copulation Factors Regulating Honey Bee (Apis mellifera) Queen Reproduction and Health: A Review

Honey bees are major pollinators of agricultural and non-agricultural landscapes. In recent years, honey bee colonies have exhibited high annual losses and commercial beekeepers frequently report poor queen quality and queen failure as the primary causes. Honey bee colonies are highly vulnerable to compromised queen fertility, as each hive is headed by one reproductive queen. Queens mate with multiple drones (male bees) during a single mating period early in life in which they obtain enough spermatozoa to fertilize their eggs for the rest of their reproductive life span. The process of mating initiates numerous behavioral, physiological, and molecular changes that shape the fertility of the queen and her influence on the colony. For example, receipt of drone semen can modulate queen ovary activation, pheromone production, and subsequent worker retinue behavior. In addition, seminal fluid is a major component of semen that is primarily derived from drone accessory glands. It also contains a complex mixture of proteins such as proteases, antioxidants, and antimicrobial proteins. Seminal fluid proteins are essential for inducing post-mating changes in other insects such as Drosophila and thus they may also impact honey bee queen fertility and health. However, the specific molecules in semen and seminal fluid that initiate post-mating changes in queens are still unidentified. Herein, we summarize the mating biology of honey bees, the changes queens undergo during and after copulation, and the role of drone semen and seminal fluid in post-mating changes in queens. We then review the effects of seminal fluid proteins in insect reproduction and potential roles for honey bee drone seminal fluid proteins in queen reproduction and health. We finish by proposing future avenues of research. Further elucidating the role of drone fertility in queen reproductive health may contribute towards reducing colony losses and advancing honey bee stock development.

Bee Viruses: Ecology, Pathogenicity, and Impacts

Bees-including solitary, social, wild, and managed species-are key pollinators of flowering plant species, including nearly three-quarters of global food crops. Their ecological importance, coupled with increased annual losses of managed honey bees and declines in populations of key wild species, has focused attention on the factors that adversely affect bee health, including viral pathogens. Genomic approaches have dramatically expanded understanding of the diversity of viruses that infect bees, the complexity of their transmission routes-including intergenus transmission-and the diversity of strategies bees have evolved to combat virus infections, with RNA-mediated responses playing a prominent role. Moreover, the impacts of viruses on their hosts are exacerbated by the other major stressors bee populations face, including parasites, poor nutrition, and exposure to chemicals. Unraveling the complex relationships between viruses and their bee hosts will lead to improved understanding of viral ecology and management strategies that support better bee health.

Honey bee viruses in Serbian colonies of different strength

Protection of honey bees is of great economic importance because of their role in pollination. Crucial steps towards this goal are epidemiological surveys of pathogens connected with honey bee losses. In this study deformed wing virus (DWV), chronic bee paralysis virus (CBPV), acute bee paralysis virus (ABPV) and sacbrood virus (SBV) were investigated in colonies of different strength located in five regions of Serbia. The relationship between colony strength and virus occurrence/infection intensity were assessed as well as the genetic relationship between virus sequences from Serbia and worldwide. Real-time RT-PCR analyses detected at least one virus in 87.33% of colonies. Single infection was found in 28.67% colonies (21.33%, 4.00%, 2.67% and 0.67% in cases of DWV, ABPV, SBV and CBPV, respectively). In the majority of colonies (58.66%) more than one virus was found. The most prevalent was DWV (74%), followed by ABPV, SBV and CBPV (49.30%, 24.00% and 6.70%, respectively). Except for DWV, the prevalence of the remaining three viruses significantly varied between the regions. No significant differences were found between colony strength and either (i) the prevalence of DWV, ABPV, SBV, CBPV and their combinations, or (ii) DWV infection levels. The sequences of honey bee viruses obtained from bees in Serbia were 93-99% identical with those deposited in GenBank.

Home sick: impacts of migratory beekeeping on honey bee (Apis mellifera) pests, pathogens, and colony size

Honey bees are important pollinators of agricultural crops and the dramatic losses of honey bee colonies have risen to a level of international concern. Potential contributors to such losses include pesticide exposure, lack of floral resources and parasites and pathogens. The damaging effects of all of these may be exacerbated by apicultural practices. To meet the pollination demand of US crops, bees are transported to areas of high pollination demand throughout the year. Compared to stationary colonies, risk of parasitism and infectious disease may be greater for migratory bees than those that remain in a single location, although this has not been experimentally established. Here, we conducted a manipulative experiment to test whether viral pathogen and parasite loads increase as a result of colonies being transported for pollination of a major US crop, California almonds. We also tested if they subsequently transmit those diseases to stationary colonies upon return to their home apiaries. Colonies started with equivalent numbers of bees, however migratory colonies returned with fewer bees compared to stationary colonies and this difference remained one month later. Migratory colonies returned with higher black queen cell virus loads than stationary colonies, but loads were similar between groups one month later. Colonies exposed to migratory bees experienced a greater increase of deformed wing virus prevalence and load compared to the isolated group. The three groups had similar infestations of Varroa mites upon return of the migratory colonies. However, one month later, mite loads in migratory colonies were significantly lower compared to the other groups, possibly because of lower number of host bees. Our study demonstrates that migratory pollination practices has varying health effects for honey bee colonies. Further research is necessary to clarify how migratory pollination practices influence the disease dynamics of honey bee diseases we describe here.

Honey bees as models for gut microbiota research

The gut microbiota of the honey bee (Apis mellifera) offers several advantages as an experimental system for addressing how gut communities affect their hosts and for exploring the processes that determine gut community composition and dynamics. A small number of bacterial species dominate the honey bee gut community. These species are restricted to bee guts and can be grown axenically and genetically manipulated. Large numbers of microbiota-free hosts can be economically reared and then inoculated with single isolates or defined communities to examine colonization patterns and effects on host phenotypes. Honey bees have been studied extensively, due to their importance as agricultural pollinators and as models for sociality. Because of this history of bee research, the physiology, development, and behavior of honey bees is relatively well understood, and established behavioral and phenotypic assays are available. To date, studies on the honey bee gut microbiota show that it affects host nutrition, weight gain, endocrine signaling, immune function, and pathogen resistance, while perturbation of the microbiota can lead to reduced host fitness. As in humans, the microbiota is concentrated in the distal part of the gut, where it contributes to digestion and fermentation of plant cell wall components. Much like the human gut microbiota, many bee gut bacteria are specific to the bee gut and can be directly transmitted between individuals through social interaction. Although simpler than the human gut microbiota, the bee gut community presents opportunities to understand the processes that govern the assembly of specialized gut communities as well as the routes through which gut communities impact host biology.

Honey Bee and Bumble Bee Antiviral Defense

Bees are important plant pollinators in both natural and agricultural ecosystems. Managed and wild bees have experienced high average annual colony losses, population declines, and local extinctions in many geographic regions. Multiple factors, including virus infections, impact bee health and longevity. The majority of bee-infecting viruses are positive-sense single-stranded RNA viruses. Bee-infecting viruses often cause asymptomatic infections but may also cause paralysis, deformity or death. The severity of infection is governed by bee host immune responses and influenced by additional biotic and abiotic factors. Herein, we highlight studies that have contributed to the current understanding of antiviral defense in bees, including the Western honey bee (Apis mellifera), the Eastern honey bee (Apis cerana) and bumble bee species (Bombus spp.). Bee antiviral defense mechanisms include RNA interference (RNAi), endocytosis, melanization, encapsulation, autophagy and conserved immune pathways including Jak/STAT (Janus kinase/signal transducer and activator of transcription), JNK (c-Jun N-terminal kinase), MAPK (mitogen-activated protein kinases) and the NF-κB mediated Toll and Imd (immune deficiency) pathways. Studies in Dipteran insects, including the model organism Drosophila melanogaster and pathogen-transmitting mosquitos, provide the framework for understanding bee antiviral defense. However, there are notable differences such as the more prominent role of a non-sequence specific, dsRNA-triggered, virus limiting response in honey bees and bumble bees. This virus-limiting response in bees is akin to pathways in a range of organisms including other invertebrates (i.e., oysters, shrimp and sand flies), as well as the mammalian interferon response. Current and future research aimed at elucidating bee antiviral defense mechanisms may lead to development of strategies that mitigate bee losses, while expanding our understanding of insect antiviral defense and the potential evolutionary relationship between sociality and immune function.