Eggs sampling as an effective tool for identifying the incidence of viruses in honey bees involved in artificial queen rearing

The Carniolan honey bee (Apis mellifera carnica) plays an essential role in crop pollination, environment diversity, and the production of honey bee products. However, the health of individual honey bees and their colonies is under pressure due to multiple stressors, including viruses as a significant threat to bees. Monitoring various virus infections could be a crucial selection tool during queen rearing. In the present study, samples from all developmental stages (eggs, larvae, pupae, and queens) were screened for the incidence of seven viruses during queen rearing in Slovenia. The screening of a total of 108 samples from five queen breeders was performed by the RT-qPCR assays. The results showed that the highest incidence was observed for black queen cell virus (BQCV), Lake Sinai virus 3 (LSV3), deformed wing virus B (DWV-B), and sacbrood virus (SBV). The highest viral load was detected in queens (6.07 log10 copies/queen) and larvae (5.50 log10 copies/larva) for BQCV, followed by SBV in larvae (5.47 log10 copies/larva). When comparing all the honey bee developmental stages, the eggs exhibited general screening for virus incidence and load in queen mother colonies. The results suggest that analyzing eggs is a good indicator of resilience to virus infection during queen development.

Population-wide modelling reveals prospects of marker-assisted selection for parasitic mite resistance in honey bees

In 2019, a joint eight-variant model was published in which eight single nucleotide polymorphisms (SNPs) in seven Apis mellifera genes were associated with Varroa destructor drone brood resistance (DBR, i.e. mite non-reproduction in drone brood). As this model was derived from only one Darwinian Black Bee Box colony, it could not directly be applied on a population-overarching scale in the northern part of Belgium (Flanders), where beekeepers prefer the carnica subspecies. To determine whether these eight SNPs remained associated with the DBR trait on a Flemish colony-broad scope, we performed population-wide modelling through sampling of various A. mellifera carnica colonies, DBR scoring of Varroa-infested drone brood and variant genotyping. Novel eight-variant modelling was performed and the classification performance of the eight SNPs was evaluated. Besides, we built a reduced three-variant model retaining only three genetic variants and found that this model classified 76% of the phenotyped drones correctly. To examine the spread of beneficial alleles and predict the DBR probability distribution in Flanders, we determined the allelic frequencies of the three variants in 292 A. mellifera carnica queens. As such, this research reveals prospects of marker-assisted selection for Varroa drone brood resistance in honeybees.

Quantitative assessments of honeybee colony's response to an artificial vibrational pulse resulting in non-invasive measurements of colony's overall mobility and restfulness

In this work we aim to provide a quantitative method allowing the probing of the physiological status of honeybee colonies by providing them with a gentle, short, external artificial vibrational shockwave, and recording their response. The knock is provided by an external electromagnetic shaker attached to the outer wall of a hive, driven by a computer with a 0.1 s long, monochromatic vibration at 340Hz set to an amplitude that occasionally yields a mild response from the bees, recorded by an accelerometer placed in the middle of the central frame of the colony. To avoid habituation, the stimulus is supplied at randomised times, approximately every hour. The method is pioneered with a pilot study on a single colony hosted indoors, then extended onto eight outdoors colonies. The results show that we can quantitatively sense the colony's overall mobility, independently from another physiological aspect, which is phenomenologically explored. Using this, a colony that is queenless is easily discriminated from the others.

System Mapping of the Production and Value Chain to Explore Beekeeping Potential in Southwest Ethiopia

Ethiopia has a high potential for the production of honey and other apiary products due to its ideal agroecology. This potential is, however, not yet well utilized due to weak production and valorization systems. The study analyzed beekeeping systems and their honey value chain to detect the barriers and to explore ways to better exploit the existing potential. Descriptive statistics, a SWOT and PESTEL matrix, and system mapping were utilized for analysis. Ethiopian beekeeping is still dominated by traditional production systems, followed by modern and transitional systems, differing in types of beehives and the average amount of honey yield. The combined SWOT-PESTEL analysis revealed challenges like a limited supply and high cost of modern beehives, shortage of credit, absence of a honey marketing legal framework, pest and predator attacks, absconding, and uncontrolled application of agrochemicals. Opportunities include the globally increasing demand for honey, availability of good investment policy, conducive agroecology, and support from NGOs. The less productive techniques of smallholder beekeepers' crude honey production for local beverage making affected the good use of the potential and minimized its contribution to the local and national economy. On the contrary, strengthening private investors and cooperatives towards the production of fully and semi-processed honey impacted the utilization of the potential positively.

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

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

Shift in virus composition in honeybees (Apis mellifera) following worldwide invasion by the parasitic mite and virus vector Varroa destructor

Invasive vectors can induce dramatic changes in disease epidemiology. While viral emergence following geographical range expansion of a vector is well known, the influence a vector can have at the level of the host's pathobiome is less well understood. Taking advantage of the formerly heterogeneous spatial distribution of the ectoparasitic mite Varroa destructor that acts as potent virus vector among honeybees Apis mellifera, we investigated the impact of its recent global spread on the viral community of honeybees in a retrospective study of historical samples. We hypothesized that the vector has had an effect on the epidemiology of several bee viruses, potentially altering their transmissibility and/or virulence, and consequently their prevalence, abundance, or both. To test this, we quantified the prevalence and loads of 14 viruses from honeybee samples collected in mite-free and mite-infested populations in four independent geographical regions. The presence of the mite dramatically increased the prevalence and load of deformed wing virus, a cause of unsustainably high colony losses. In addition, several other viruses became more prevalent or were found at higher load in mite-infested areas, including viruses not known to be actively varroa-transmitted, but which may increase opportunistically in varroa-parasitized bees.

Convivina is a specialised core gut symbiont of the invasive hornet Vespa velutina

We provide a culturomics analysis of the cultivable bacterial communities of the crop, midgut and hindgut compartments, as well as the ovaries, of the invasive insect Vespa velutina, along with a cultivation-independent analysis of samples of the same nest through 16S rRNA amplicon sequencing. The Vespa velutina bacterial symbiont community was dominated by the genera Convivina, Fructobacillus, Lactiplantibacillus, Lactococcus, Sphingomonas and Spiroplasma. Lactococcus lactis and Lactiplantibacillus plantarum represented generalist core lactic acid bacteria (LAB) symbionts, while Convivina species and Fructobacillus fructosus represented highly specialised core LAB symbionts with strongly reduced genome sizes. Sphingomonas and Spiroplasma were the only non-LAB core symbionts but were not isolated. Convivina bacteria were particularly enriched in the hornet crop and included Convivina intestini, a species adapted towards amino acid metabolism, and Convivina praedatoris sp. nov. which was adapted towards carbohydrate metabolism.

Beekeeping in Europe facing climate change: A mixed methods study on perceived impacts and the need to adapt according to stakeholders and beekeepers

The beekeeping sector is suffering from the detrimental effects of climate change, both directly and indirectly. Despite numerous studies conducted on this subject, large-scale research incorporating stakeholders' and beekeepers' perspectives has remained elusive. This study aims to bridge this gap by assessing the extent to which stakeholders involved in the European beekeeping sector and European beekeepers perceive and experience the impacts of climate change on their operations, and whether they had to adapt their practices accordingly. To this end, a mixed-methods study including in-depth stakeholder interviews (n = 41) and a pan-European beekeeper survey (n = 844) was completed within the frame of the EU-funded H2020-project B-GOOD. The development of the beekeeper survey was informed by insights from literature and the stakeholder interviews. The results highlighted significant regional disparities in the perceived impacts of climate change, with beekeepers in Southern European regions expressing more negative outlooks, while Northern European beekeepers reported more favourable experiences. Furthermore, survey analysis revealed beekeepers who were classified as 'heavily impacted' by climate change. These beekeepers reported lower average honey yields, higher colony winter loss rates and a stronger perceived contribution of honey bees to pollination and biodiversity, underscoring climate change's detrimental impacts on the beekeeping sector. Multinomial logistic regression revealed determinants of the likelihood of beekeepers being classified as 'heavily impacted' by climate change. This analysis indicates that Southern European beekeepers experienced a 10-fold likelihood of being classified as heavily impacted by climate change compared to Northern European beekeepers. Other significant factors distinguishing 'winners' and 'losers' were self-reported level of professionalism as a beekeeper (ranging from pure hobbyist to fully professional, Odds Ratio (OR) = 1.31), number of years active in beekeeping (OR = 1.02), availability of floral resources throughout the bee season (OR = 0.78), beehives located in a forested environment (OR = 1.34), and the presence of local policy measures addressing climate change-related challenges (OR = 0.76).

Expression of Molecular Markers of Resilience against Varroa destructor and Bee Viruses in Ethiopian Honey Bees (Apis mellifera simensis) Focussing on Olfactory Sensing and the RNA Interference Machinery

Varroa destructor mites and the viruses it vectors are two major factors leading to high losses of honey bees (Apis mellifera) colonies worldwide. However, honey bees in some African countries show resilience to varroa infestation and/or virus infections, although little is known about the mechanisms underlying this resilience. In this study, we investigated the expression profiles of some key molecular markers involved in olfactory sensing and RNA interference, as these processes may contribute to the bees' resilience to varroa infestation and virus infection, respectively. We found significantly higher gene expression of the odorant binding protein, OBP14, in the antennae of Ethiopian bees compared to Belgian bees. This result suggests the potential of OBP14 as a molecular marker of resilience to mite infestation. Scanning electron microscopy showed no significant differences in the antennal sensilla occurrence and distribution, suggesting that resilience arises from molecular processes rather than morphological adaptations. In addition, seven RNAi genes were upregulated in the Ethiopian honey bees and three of them-Dicer-Drosha, Argonaute 2, and TRBP2-were positively correlated with the viral load. We can conclude that the antiviral immune response was triggered when bees were experiencing severe viral infection and that this might contribute to the bees' resilience to viruses.

Cold case: The disappearance of Egypt bee virus, a fourth distinct master strain of deformed wing virus linked to honeybee mortality in 1970's Egypt

In 1977, a sample of diseased adult honeybees (Apis mellifera) from Egypt was found to contain large amounts of a previously unknown virus, Egypt bee virus, which was subsequently shown to be serologically related to deformed wing virus (DWV). By sequencing the original isolate, we demonstrate that Egypt bee virus is in fact a fourth unique, major variant of DWV (DWV-D): more closely related to DWV-C than to either DWV-A or DWV-B. DWV-A and DWV-B are the most common DWV variants worldwide due to their close relationship and transmission by Varroa destructor. However, we could not find any trace of DWV-D in several hundred RNA sequencing libraries from a worldwide selection of honeybee, varroa and bumblebee samples. This means that DWV-D has either become extinct, been replaced by other DWV variants better adapted to varroa-mediated transmission, or persists only in a narrow geographic or host range, isolated from common bee and beekeeping trade routes.

Risk and protective indicators of beekeeping management practices

Explaining the reasons for the high honey bee (Apis mellifera) colony loss rate in recent years has become a top global research priority in apicultural and agricultural sciences. Although there are indications of the role played by beekeeping management practices on honey bee health, very little information is currently available. Our study aimed to characterize the beekeeping management practices carried out in Belgium, and to determine the relationship between beekeeping management practices and colony losses. Variables obtained from face-to-face questioning of a representative randomized and stratified sample of Belgian beekeepers (n = 186) were integrated into a logistic regression model (univariate and multivariate) and correlated to the declared colony loss rates to identify risk and protective indicators. We used a classification tree analysis to validate the results. We present evidence of a relationship between poor beekeeping management practices and colony losses. The main factors protecting honey bee colonies are the aptitude of the beekeeper to change his management practices, the hive type, the equipment origin and hygiene, wintering in proper conditions (the use of divider boards, i.e. board blocks or space fillers off part of the hive body), the colony strength estimation before wintering, winter monitoring, and last but not least, appropriate integrated pest management. Proper estimation of the Varroa infestation level should be performed prior to treatment. The consequences of poor beekeeping practices on honey bee health can be addressed by proper training of beekeepers. An online tool was developed and published for beekeepers allowing them to evaluate the effect of their management practices on colony health.

Pyrethroid target-site resistance mutations in populations of the honey bee parasite Varroa destructor (Acari: Varroidae) from Flanders, Belgium

The honey bee ectoparasite Varroa destructor is considered the major threat to apiculture, as untreated colonies of Apis mellifera usually collapse within a few years. In order to control this mite, many beekeepers rely on a limited number of approved synthetic acaricides, including the pyrethroids tau-fluvalinate and flumethrin. Due to the intensive use of these products, resistance is now commonplace in many beekeeping regions across the world. In the present study, the occurrence of amino acid substitutions at residue L925 of the voltage-gate sodium channel-the pyrethroid target site-was studied in Varroa populations collected throughout Flanders, Belgium. Dose-response bioassays supported the involvement of the frequently observed L925V substitution in flumethrin resistance, resulting in a 12.64-fold increase of the LC₅₀ in a Varroa population mostly consisting of homozygous 925 V/V mites. With the presence of L925 substitutions in about four out of 10 screened apiaries, the use of pyrethroid-based varroacides in Flanders, including the recently released PolyVar® Yellow, should be carefully considered.

Honey bee exposure scenarios to selected residues through contaminated beeswax

Twenty-two pesticides and veterinary drugs of which residues were detected in beeswax in Europe were selected according to different criteria. The risk to honey bee health posed by the presence of these residues in wax was assessed based on three exposure scenarios. The first one corresponds to the exposure of larvae following their close contact with wax constituting the cells in which they develop. The second one corresponds to the exposure of larvae following consumption of the larval food that was contaminated from contact with contaminated wax. The third one corresponds to the exposure of adult honey bees following wax chewing when building cells and based on a theoretical worst-case scenario (= intake of contaminants from wax). Following these three scenarios, maximum concentrations which should not be exceeded in beeswax in order to protect honey bee health were calculated for each selected substance. Based on these values, provisional action limits were proposed. Beeswax exceeding these limits should not be put on the market.

Breeding for Virus Resistance and Its Effects on Deformed Wing Virus Infection Patterns in Honey Bee Queens

Viruses, and in particular the deformed wing virus (DWV), are considered as one of the main antagonists of honey bee health. The 'suppressed in ovo virus infection' trait (SOV) described for the first time that control of a virus infection can be achieved from genetically inherited traits and that the virus state of the eggs is indicative for this. This research aims to explore the effect of the SOV trait on DWV infections in queens descending from both SOV-positive (QDS+) and SOV-negative (QDS-) queens. Twenty QDS+ and QDS- were reared from each time four queens in the same starter-finisher colony. From each queen the head, thorax, ovaries, spermatheca, guts and eviscerated abdomen were dissected and screened for the presence of the DWV-A and DWV-B genotype using qRT-PCR. Queens descending from SOV-positive queens showed significant lower infection loads for DWV-A and DWV-B as well as a lower number of infected tissues for DWV-A. Surprisingly, differences were less expressed in the reproductive tissues, the ovaries and spermatheca. These results confirm that selection on the SOV trait is associated with increased virus resistance across viral genotypes and that this selection drives DWV towards an increased tissue specificity for the reproductive tissues. Further research is needed to explore the mechanisms underlying the interaction between the antiviral response and DWV.

qPCR assays with dual-labeled probes for genotyping honey bee variants associated with varroa resistance

Background

The varroa mite is one of the main causes of honey bee mortality. An important mechanism by which honey bees increase their resistance against this mite is the expression of suppressed mite reproduction. This trait describes the physiological inability of mites to produce viable offspring and was found associated with eight genomic variants in previous research.

Results

This paper presents the development and validation of high-throughput qPCR assays with dual-labeled probes for discriminating these eight single-nucleotide variants. Amplicon sequences used for assay validation revealed additional variants in the primer/probe binding sites in four out of the eight assays. As for two of these the additional variants interfered with the genotyping outcome supplementary primers and/or probes were developed. Inclusion of these primers and probes in the assay mixes allowed for the correct genotyping of all eight variants of interest within our bee population.

Conclusion

These outcomes underline the importance of checking for interfering variants in designing qPCR assays. Ultimately, the availability of this assay allows genotyping for the suppressed mite reproduction trait and paves the way for marker assisted selection in breeding programs.

Reduced nest development of reared Bombus terrestris within apiary dense human-modified landscapes

Wild bees are in decline on a local to global scale. The presence of managed honey bees can lead to competition for resources with wild bee species, which has not been investigated so far for human-modified landscapes. In this study we assess if managed honey bee hive density influence nest development (biomass) of bumble bees, an important trait affecting fitness. We hypothesize that domesticated honey bees can negatively affect Bombus terrestris nest development in human-modified landscapes. In Flanders, Belgium, where such landscapes are dominantly present, we selected 11 locations with landscape metrics ranging from urban to agricultural. The bee hive locations were mapped and each location contained one apiary dense (AD) and one apiary sparse (AS) study site (mean density of 7.6 ± 5.7 managed honey bee hives per km² in AD sites). We assessed the effect of apiary density on the reproduction of reared B. terrestris nests. Reared B. terrestris nests had more biomass increase over 8 weeks in apiary sparse (AS) sites compared to nests located in apiary dense (AD) sites. This effect was mainly visible in urban locations, where nest in AS sites have 99.25 ± 60.99 g more biomass increase compared to nest in urban AD sites. Additionally, we found that managed bumble bee nests had higher biomass increase in urban locations. We conclude that the density of bee hives is a factor to consider in regard to interspecific competition between domesticated honey bees and bumble bees.

Beekeepers perception of risks affecting colony loss: A pilot survey

Understanding amateur beekeepers' perception of risks affecting bee health and mortality is essential to analyse the reasons for adopting or rejecting good management practices. A perception survey on how beekeepers perceive and manage factors related to climate change, Varroa infestation, management practices, and pesticide exposure was designed and launched online. This unpreceded sociological survey involved 355 beekeepers spread all over Belgium. A two-sample t test with unequal variances comparing beekeepers with colony loss rates below or exceeding the acceptable level, that is <10% and ≥10%, indicates that beekeepers (N = 213) with colony loss rates <10% generally have greater average levels of perceived risks and the benefits of action that lead to increased motivation to act in better ways. The results of this survey highlight the importance of looking beyond socio-economic determinants in any risk mitigation strategy associated with bee mortality when dealing with amateur beekeepers.

Pesticide and veterinary drug residues in Belgian beeswax: Occurrence, toxicity, and risk to honey bees

Pesticide and veterinary drug residues are one of the stress factors affecting bee health and mortality. To investigate the occurrence, the concentration and the toxicity risk to bees of pesticide residues in four different types of beeswax (brood comb wax, recycled comb wax, honey comb wax, and cappings wax), 182 samples were collected from apiaries located all over the Belgian territories, during spring 2016 and analysed by LC-MS/MS and GC-MS/MS for the presence of 294 chemical residues. The toxicity risk to bees expressed as the Hazard Quotient (HQ) was calculated for each wax sample, according to two scenarios with different tau-fluvalinate LD₅₀ values. Residues showing the highest prevalence were correlated to bee mortality in a multivariate logistic regression model and a risk-based model was used to predict colony bee mortality. Altogether, 54 different pesticide and veterinary drug residues were found in the four types of beeswax. The residues with a higher likelihood to be retained in beeswax are applied in-hive or with a high lipophilic nature. The multivariate logistic regression model showed a statistically significant influence of chlorfenvinphos on bee mortality. All our results indicated that cappings wax was substantially less contaminated. This national survey on beeswax contamination provides guidelines on the re-use of beeswax by beekeepers and shows the necessity to introduce maximum residue levels for global trade in beeswax. An online tool was developed to enable beekeepers and wax traders to estimate the risk to honey bee health associated with contaminated wax.

Nationwide Screening for Bee Viruses and Parasites in Belgian Honey Bees

The health of honey bees is threatened by multiple factors, including viruses and parasites. We screened 557 honey bee (Apis mellifera) colonies from 155 beekeepers distributed all over Belgium to determine the prevalence of seven widespread viruses and two parasites (Varroa sp. and Nosema sp.). Deformed wing virus B (DWV-B), black queen cell virus (BQCV), and sacbrood virus (SBV) were highly prevalent and detected by real-time RT-PCR in more than 95% of the colonies. Acute bee paralysis virus (ABPV), chronic bee paralysis virus (CBPV) and deformed wing virus A (DWV-A) were prevalent to a lower extent (between 18 and 29%). Most viruses were only present at low or moderate viral loads. Nevertheless, about 50% of the colonies harbored at least one virus at high viral load (>10⁷ genome copies/bee). Varroa mites and Nosema sp. were found in 81.5% and 59.7% of the honey bee colonies, respectively, and all Nosema were identified as Nosema ceranae by real time PCR. Interestingly, we found a significant correlation between the number of Varroa mites and DWV-B viral load. To determine the combined effect of these and other factors on honey bee health in Belgium, a follow up of colonies over multiple years is necessary.

Comparison between Apicystis cryptica sp. n. and Apicystis bombi (Arthrogregarida, Apicomplexa): Gregarine parasites that cause fat body hypertrophism in bees

The molecular divergence, morphology and pathology of a cryptic gregarine that is related to the bee parasite Apicystis bombi Lipa and Triggiani, 1996 is described. The 18S ribosomal DNA gene sequence of the new gregarine was equally dissimilar to that of A. bombi and the closest related genus Mattesia Naville, 1930, although phylogenetic analysis supported a closer relation to A. bombi. Pronounced divergence with A. bombi was found in the ITS1 sequence (69.6% similarity) and seven protein-coding genes (nucleotide 78.05% and protein 90.2% similarity). The new gregarine was isolated from a Bombus pascuorum Scopoli, 1763 female and caused heavy hypertrophism of the fat body tissue in its host. In addition, infected cells of the hypopharyngeal gland tissue, an important excretory organ of the host, were observed. Mature oocysts were navicular in shape and contained four sporozoites, similar to A. bombi oocysts. Given these characteristics, we proposed the name Apicystis cryptica sp. n. Detections so far indicated that distribution and host species occupation of Apicystis spp. overlap at least in Europe, and that historical detections could not discriminate between them. Specific molecular assays were developed that can be implemented in future pathogen screens that aim to discriminate Apicystis spp. in bees.

Honeybee and consumer's exposure and risk characterisation to glyphosate-based herbicide (GBH) and its degradation product (AMPA): Residues in beebread, wax, and honey

In order to assess bee and human exposure to residues of glyphosate-based herbicide (GBH) and its main degradation products aminomethylphosphonic acid (AMPA) and to characterise the risk posed by these substances, we analysed 3 different bee matrices; beebread (N = 81), wax (N = 100) and 10-paired samples of wax/honey collected in 2016/2017 from 379 Belgian apiaries. A high-performance liquid chromatography-electrospray ionisation tandem mass spectrometry (HPLC-ESI-MS-MS) was used as analytical method. Limit of quantification and detection (LOQ and LOD) for GBH residues and AMPA in the 3 matrices was respectively of 10 ng g-1 and 1 ng g-1. In beebread, 81.5% of the samples showed a residue concentration > LOQ and 9.9% of the samples a residue concentration < LOQ (detection without quantification); no significant difference in detection rate was found between the north and the south of the country. Glyphosate was detected in beeswax less frequently than in beebread (i.e. 26% >LOQ versus 81.5% >LOQ). The maximum GBH residues and AMPA concentration found in beebread (respectively 700 ng g-1 and 250 ng g-1) led to sub-lethal exposure to bees. The Hazard Quotient (HQ) for beebread and beeswax (7 and 3.2, respectively) were far below the "safety" oral and contact thresholds for bees. For human health, the highest exposure to GBH residues in pollen corresponded to 0.312% and 0.187% of the ADI and of the ARfD respectively and, to 0.002% and to 0.001% for beeswax. No transfer of glyphosate from wax to honey was detected. Considering our results and the available regulatory data on the glyphosate molecule considered solely, not including the adjuvants in GBH formulation, the consumption of these three contaminated matrices would not be a food safety issue. Nonetheless, caution should be taken in the interpretation of the results as new studies indicate possible glyphosate/GBH residues toxicity below regulatory limits and at chronic sub-lethal doses.

Belgian case study on flumethrin residues in beeswax: Possible impact on honeybee and prediction of the maximum daily intake for consumers

To assess the health risk posed by flumethrin residues in beeswax to honeybees and honey consumers, 124 wax samples randomly distributed in Belgium were analysed for flumethrin residues using liquid chromatography/tandem mass spectrometry. The risk posed by flumethrin residues in beeswax to honeybee health was assessed through the calculation of a non-pondered and a pondered Hazard Quotient by the prevalence rate of flumethrin considering an oral or topical exposure. No statistical difference was found when comparing both the average flumethrin residues concentrations and contact and oral pondered hazard quotients between apiaries with lower and equal or higher than 10% of colony loss. Flumethrin residues estimated daily intake by Belgian consumers through honey and wax ingestion was estimated via a deterministic (worst-case scenario) and a probabilistic approach. The probabilistic approach was not possible for beeswax consumption due to the lack of individual consumption data. The highest estimated exposure was <0.1% of the theoretical maximum daily intake for both approaches, meaning no risk for human health.

Characterization of the honeybee venom proteins C1q-like protein and PVF1 and their allergenic potential

Honeybee (Apis mellifera) venom (HBV) represents an ideal model to study the role of particular venom components in allergic reactions in sensitized individuals as well as in the eusociality of Hymenoptera species. The aim of this study was to further characterize the HBV components C1q-like protein (C1q) and PDGF/VEGF-like factor 1 (PVF1). C1q and PVF1 were produced as recombinant proteins in insect cells. Their allergenic properties were examined by determining the level of specific IgE antibodies in the sera of HBV-allergic patients (n = 26) as well as by their capacity to activate patients' basophils (n = 11). Moreover, the transcript heterogeneity of PVF1 was analyzed. It could be demonstrated that at least three PVF1 variants are present in the venom gland, which all result from alternative splicing of one transcript. Additionally, recombinant C1q and PVF1 from Spodoptera frugiperda insect cells exhibited specific IgE reactivity with approximately 38.5% of sera of HBV-allergic patients. Interestingly, both proteins were unable to activate basophils of the patients, questioning their role in the context of clinically relevant sensitization. Recombinant C1q and PVF1 can build the basis for a deeper understanding of the molecular mechanisms of Hymenoptera venoms. Moreover, the conflicting results between IgE sensitization and lack of basophil activation, might in the future contribute to the identification of factors that determine the allergenic potential of proteins.

Study of the Metatranscriptome of Eight Social and Solitary Wild Bee Species Reveals Novel Viruses and Bee Parasites

Bees are associated with a remarkable diversity of microorganisms, including unicellular parasites, bacteria, fungi, and viruses. The application of next-generation sequencing approaches enables the identification of this rich species composition as well as the discovery of previously unknown associations. Using high-throughput polyadenylated ribonucleic acid (RNA) sequencing, we investigated the metatranscriptome of eight wild bee species (Andrena cineraria, Andrena fulva, Andrena haemorrhoa, Bombus terrestris, Bombus cryptarum, Bombus pascuorum, Osmia bicornis, and Osmia cornuta) sampled from four different localities in Belgium. Across the RNA sequencing libraries, 88–99% of the taxonomically informative reads were of the host transcriptome. Four viruses with homology to insect pathogens were found including two RNA viruses (belonging to the families Iflaviridae and Tymoviridae that harbor already viruses of honey bees), a double stranded DNA virus (family Nudiviridae) and a single stranded DNA virus (family Parvoviridae). In addition, we found genomic sequences of 11 unclassified arthropod viruses (related to negeviruses, sobemoviruses, totiviruses, rhabdoviruses, and mononegaviruses), seven plant pathogenic viruses, and one fungal virus. Interestingly, nege-like viruses appear to be widespread, host-specific, and capable of attaining high copy numbers inside bees. Next to viruses, three novel parasite associations were discovered in wild bees, including Crithidia pragensis and a tubulinosematid and a neogregarine parasite. Yeasts of the genus Metschnikowia were identified in solitary bees. This study gives a glimpse of the microorganisms and viruses associated with social and solitary wild bees and demonstrates that their diversity exceeds by far the subset of species first discovered in honey bees.

Environmental contaminants of honeybee products in Uganda detected using LC-MS/MS and GC-ECD

Pollinator services and the development of beekeeping as a poverty alleviating tool have gained considerable focus in recent years in sub-Saharan Africa. An improved understanding of the pervasive environmental extent of agro-chemical contaminants is critical to the success of beekeeping development and the production of clean hive products. This study developed and validated a multi-residue method for screening 36 pesticides in honeybees, honey and beeswax using LC-MS/MS and GC-ECD. Of the 36 screened pesticides, 20 were detected. The highest frequencies occurred in beeswax and in samples from apiaries located in the proximity of citrus and tobacco farms. Fungicides were the most prevalent chemical class. Detected insecticides included neonicotinoids, organophosphates, carbamates, organophosphorus, tetrazines and diacylhydrazines. All detected pesticide levels were below maximum residue limits (according to EU regulations) and the lethal doses known for honeybees. However, future risk assessment is needed to determine the health effects on the African genotype of honeybees by these pesticide classes and combinations of these. In conclusion, our data present a significant challenge to the burgeoning organic honey sector in Uganda, but to achieve this, there is an urgent need to regulate the contact routes of pesticides into the beehive products. Interestingly, the "zero" detection rate of pesticides in the Mid-Northern zone is a significant indicator of the large potential to promote Ugandan organic honey for the export market.

The buzz about bees and poverty alleviation: Identifying drivers and barriers of beekeeping in sub-Saharan Africa

The potential of beekeeping to mitigate the exposure of rural sub-Sahara African farmers to economic stochasticity has been widely promoted by an array of development agencies. Robust outcome indicators of the success of beekeeping to improve household well-being are unfortunately lacking. This study aimed to identify the key drivers and barriers of beekeeping adoption at the household level, and quantified the associated income contribution in three agro-ecological zones in Uganda. Beekeepers were generally the most economically disadvantaged people in the study areas and tended to adopt beekeeping following contact with non-government organisations and access to training. Whilst incomes were not statistically lower than their non-beekeeping counterparts; their mean household well-being scores were significantly lower than non-beekeeping households. The inability of beekeeping to significantly improve well-being status can in part be attributed to a lack of both training in bee husbandry and protective equipment provision such as suits, gloves and smokers. These are critical tools for beekeepers as they provide the necessary confidence to manage honey bees. Rather than focussing solely on the socio-economic conditions of farmers to effectively adopt beekeeping, future research should also attempt to evaluate the effectiveness of development agencies’ provision to the beekeeping sector.

Stress indicator gene expression profiles, colony dynamics and tissue development of honey bees exposed to sub-lethal doses of imidacloprid in laboratory and field experiments

In this study, different context-dependent effects of imidacloprid exposure on the honey bee response were studied. Honey bees were exposed to different concentrations of imidacloprid during a time period of 40 days. Next to these variables, a laboratory-field comparison was conducted. The influence of the chronic exposure on gene expression levels was determined using an in-house developed microarray targeting different immunity-related and detoxification genes to determine stress-related gene expression changes. Increased levels of the detoxification genes encoding, CYP9Q3 and CYT P450, were detected in imidacloprid-exposed honey bees. The different context-dependent effects of imidacloprid exposure on honey bees were confirmed physiologically by decreased hypopharyngeal gland sizes. Honey bees exposed to imidacloprid in laboratory cages showed a general immunosuppression and no detoxification mechanisms were triggered significantly, while honey bees in-field showed a resilient response with an immune stimulation at later time points. However, the treated colonies had a brood and population decline tendency after the first brood cycle in the field. In conclusion, this study highlighted the different context-dependent effects of imidacloprid exposure on the honey bee response. These findings warn for possible pitfalls concerning the generalization of results based on specific experiments with short exposure times. The increased levels of CYT P450 and CYP9Q3 combined with an immune response reaction can be used as markers for bees which are exposed to pesticides in the field.

Covert deformed wing virus infections have long-term deleterious effects on honeybee foraging and survival

Several studies have suggested that covert stressors can contribute to bee colony declines. Here we provide a novel case study and show using radiofrequency identification tracking technology that covert deformed wing virus (DWV) infections in adult honeybee workers seriously impact long-term foraging and survival under natural foraging conditions. In particular, our experiments show that adult workers injected with low doses of DWV experienced increased mortality rates, that DWV caused workers to start foraging at a premature age, and that the virus reduced the workers' total activity span as foragers. Altogether, these results demonstrate that covert DWV infections have strongly deleterious effects on honeybee foraging and survival. These results are consistent with previous studies that suggested DWV to be an important contributor to the ongoing bee declines in Europe and the USA. Overall, our study underlines the strong impact that covert pathogen infections can have on individual and group-level performance in bees.

Unbiased RNA Shotgun Metagenomics in Social and Solitary Wild Bees Detects Associations with Eukaryote Parasites and New Viruses

The diversity of eukaryote organisms and viruses associated with wild bees remains poorly characterized in contrast to the well-documented pathosphere of the western honey bee, Apis mellifera. Using a deliberate RNA shotgun metagenomic sequencing strategy in combination with a dedicated bioinformatics workflow, we identified the (micro-)organisms and viruses associated with two bumble bee hosts, Bombus terrestris and Bombus pascuorum, and two solitary bee hosts, Osmia cornuta and Andrena vaga. Ion Torrent semiconductor sequencing generated approximately 3.8 million high quality reads. The most significant eukaryote associations were two protozoan, Apicystis bombi and Crithidia bombi, and one nematode parasite Sphaerularia bombi in bumble bees. The trypanosome protozoan C. bombi was also found in the solitary bee O. cornuta. Next to the identification of three honey bee viruses Black queen cell virus, Sacbrood virus and Varroa destructor virus-1 and four plant viruses, we describe two novel RNA viruses Scaldis River bee virus (SRBV) and Ganda bee virus (GABV) based on their partial genomic sequences. The novel viruses belong to the class of negative-sense RNA viruses, SRBV is related to the order Mononegavirales whereas GABV is related to the family Bunyaviridae. The potential biological role of both viruses in bees is discussed in the context of recent advances in the field of arthropod viruses. Further, fragmentary sequence evidence for other undescribed viruses is presented, among which a nudivirus in O. cornuta and an unclassified virus related to Chronic bee paralysis virus in B. terrestris. Our findings extend the current knowledge of wild bee parasites in general and addsto the growing evidence of unexplored arthropod viruses in valuable insects.

Residues in Beeswax: A Health Risk for the Consumer of Honey and Beeswax?

A scenario analysis in regard to the risk of chronic exposure of consumers to residues through the consumption of contaminated honey and beeswax was conducted. Twenty-two plant protection products and veterinary substances of which residues have already been detected in beeswax in Europe were selected. The potential chronic exposure was assessed by applying a worst-case scenario based on the addition of a "maximum" daily intake through the consumption of honey and beeswax to the theoretical maximum daily intake through other foodstuffs. For each residue, the total exposure was finally compared to the acceptable daily intake. It is concluded that the food consumption of honey and beeswax contaminated with these residues considered separately does not compromise the consumer's health, provided proposed action limits are met. In regard to residues of flumethrin in honey and in beeswax, "zero tolerance" should be applied.

Multiple Locus Variable number of tandem repeat Analysis: A molecular genotyping tool for Paenibacillus larvae

American Foulbrood, caused by Paenibacillus larvae, is the most severe bacterial disease of honey bees (Apis mellifera). To perform genotyping of P. larvae in an epidemiological context, there is a need of a fast and cheap method with a high resolution. Here, we propose Multiple Locus Variable number of tandem repeat Analysis (MLVA). MLVA has been used for typing a collection of 209 P. larvae strains from which 23 different MLVA types could be identified. Moreover, the developed methodology not only permits the identification of the four Enterobacterial Repetitive Intergenic Consensus (ERIC) genotypes, but allows also a discriminatory subdivision of the most dominant ERIC type I and ERIC type II genotypes. A biogeographical study has been conducted showing a significant correlation between MLVA genotype and the geographical region where it was isolated.

First detection of Paenibacillus larvae the causative agent of American Foulbrood in a Ugandan honeybee colony

Paenibacillus larvae is a highly contagious and often lethal widely distributed pathogen of honeybees, Apis mellifera but has not been reported in eastern Africa to date. We investigated the presence of P. larvae in the eastern and western highland agro-ecological zones of Uganda by collecting brood and honey samples from 67 honeybee colonies in two sampling occasions and cultivated them for P. larvae. Also, 8 honeys imported and locally retailed in Uganda were sampled and cultivated for P. larvae. Our aim was to establish the presence and distribution of P. larvae in honeybee populations in the two highland agro-ecological zones of Uganda and to determine if honeys that were locally retailed contained this lethal pathogen. One honeybee colony without clinical symptoms for P. larvae in an apiary located in a protected area of the western highlands of Uganda was found positive for P. larvae. The strain of this P. larvae was genotyped and found to be ERIC I. In order to compare its virulence with P. larvae reference strains, in vitro infection experiments were conducted with carniolan honeybee larvae from the research laboratory at Ghent University, Belgium. The results show that the virulence of the P. larvae strain found in Uganda was at least equally high. The epidemiological implication of the presence of P. larvae in a protected area is discussed.

Factors influencing the prevalence and infestation levels of Varroa destructor in honeybee colonies in two highland agro-ecological zones of Uganda

Varroa mites are ecto-parasites of honeybees and are a threat to the beekeeping industry. We identified the haplotype of Varroa mites and evaluated potential factors that influence their prevalence and infestation levels in the eastern and western highland agro-ecological zones of Uganda. This was done by collecting samples of adult worker bees between December 2014 and September 2015 in two sampling moments. Samples of bees were screened for Varroa using the ethanol wash method and the mites were identified by molecular techniques. All DNA sequences obtained from sampled mite populations in the two zones were 100 % identical to the Korean Haplotype (AF106899). Mean mite prevalence in the apiaries was 40 and 53 % for the western and eastern zones, respectively, during the first sampling. Over the second sampling, mean mite prevalence increased considerably in the western (59 %) but not in the eastern (51 %) zone. Factors that were associated with Varroa mite infestation levels include altitude, nature of apiary slope and apiary management practices during the first sampling. Our results further showed that Varroa mites were spreading from lower to higher elevations. Feral colonies were also infested with Varroa mites at infestation levels not significantly different from those in managed colonies. Colony productivity and strength were not correlated to mite infestation levels. We recommend a long-term Varroa mite monitoring strategy in areas of varying landscape and land use factors for a clear understanding of possible changes in mite infestation levels among African honeybees for informed decision making.

Honeybee venom proteome profile of queens and winter bees as determined by a mass spectrometric approach

Venoms of invertebrates contain an enormous diversity of proteins, peptides, and other classes of substances. Insect venoms are characterized by a large interspecific variation resulting in extended lists of venom compounds. The venom composition of several hymenopterans also shows different intraspecific variation. For instance, venom from different honeybee castes, more specifically queens and workers, shows quantitative and qualitative variation, while the environment, like seasonal changes, also proves to be an important factor. The present study aimed at an in-depth analysis of the intraspecific variation in the honeybee venom proteome. In summer workers, the recent list of venom proteins resulted from merging combinatorial peptide ligand library sample pretreatment and targeted tandem mass spectrometry realized with a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS/MS). Now, the same technique was used to determine the venom proteome of queens and winter bees, enabling us to compare it with that of summer bees. In total, 34 putative venom toxins were found, of which two were never described in honeybee venoms before. Venom from winter workers did not contain toxins that were not present in queens or summer workers, while winter worker venom lacked the allergen Api m 12, also known as vitellogenin. Venom from queen bees, on the other hand, was lacking six of the 34 venom toxins compared to worker bees, while it contained two new venom toxins, in particularly serine proteinase stubble and antithrombin-III. Although people are hardly stung by honeybees during winter or by queen bees, these newly identified toxins should be taken into account in the characterization of a putative allergic response against Apis mellifera stings.

Differential diagnosis of the honey bee trypanosomatids Crithidia mellificae and Lotmaria passim

Trypanosomatids infecting honey bees have been poorly studied with molecular methods until recently. After the description of Crithidia mellificae (Langridge and McGhee, 1967) it took about forty years until molecular data for honey bee trypanosomatids became available and were used to identify and describe a new trypanosomatid species from honey bees, Lotmaria passim (Evans and Schwarz, 2014). However, an easy method to distinguish them without sequencing is not yet available. Research on the related bumble bee parasites Crithidia bombi and Crithidia expoeki revealed a fragment length polymorphism in the internal transcribed spacer 1 (ITS1), which enabled species discrimination. In search of fragment length polymorphisms for differential diagnostics in honey bee trypanosomatids, we studied honey bee trypanosomatid cell cultures of C. mellificae and L. passim. This research resulted in the identification of fragment length polymorphisms in ITS1 and ITS1-2 markers, which enabled us to develop a diagnostic method to differentiate both honey bee trypanosomatid species without the need for sequencing. However, the amplification success of the ITS1 marker depends probably on the trypanosomatid infection level. Further investigation confirmed that L. passim is the dominant species in Belgium, Japan and Switzerland. We found C. mellificae only rarely in Belgian honey bee samples, but not in honey bee samples from other countries. C. mellificae was also detected in mason bees (Osmia bicornis and Osmia cornuta) besides in honey bees. Further, the characterization and comparison of additional markers from L. passim strain SF (published as C. mellificae strain SF) and a Belgian honey bee sample revealed very low divergence in the 18S rRNA, ITS1-2, 28S rRNA and cytochrome b sequences. Nevertheless, a variable stretch was observed in the gp63 virulence factor.

Unraveling the venom proteome of the bumblebee (Bombus terrestris) by integrating a combinatorial peptide ligand library approach with FT-ICR MS

Within the Apidae, the largest family of bees with over 5600 described species, the honeybee is the sole species with a well studied venom proteome. So far, only little research has focused on bumblebee venom. Recently, the genome sequence of the European large earth bumblebee (Bombus terrestris) became available and this allowed the first in-depth proteomic analysis of its venom composition. We identified 57 compounds, with 52 of them never described in bumblebee venom. Remarkably, 72% of the detected compounds were found to have a honeybee venom homolog, which reflects the similar defensive function of both venoms and the high degree of homology between both genomes. However, both venoms contain a selection of species-specific toxins, revealing distinct damaging effects that may have evolved in response to species-specific attackers. Further, this study extends the list of potential venom allergens. The availability of both the honeybee and bumblebee venom proteome may help to develop a strategy that solves the current issue of false double sensitivity in allergy diagnosis, which is caused by cross-reactivity between both venoms. A correct diagnosis is important as it is recommended to perform an immunotherapy with venom of the culprit species.

Pesticides for apicultural and/or agricultural application found in Belgian honey bee wax combs

In a Belgian pilot study honey bee wax combs from ten hives were analyzed on the presence of almost 300 organochlorine and organophosphorous compounds by LC-MS/MS and GC-MS/MS. Traces of 18 pesticides were found and not a single sample was free of residues. The number of residues found per sample ranged from 3 to 13, and the pesticides found could be categorized as (1) pesticides for solely apicultural (veterinary) application, (2) pesticides for solely agricultural (crop protection) application, (3) pesticides for mixed agricultural and apicultural (veterinary) application. The frequencies and quantities of some environmental pollutants bear us high concerns. Most alarming was the detection of lindane (gamma-HCH) and dichlorodiphenyltrichloroethane (including its breakdown product dichlorodiphenyldichloroethylene), two insecticides that are banned in Europe. The present comprehensive residue analysis, however, also reveals residues of pesticides never found in beeswax before, i.e. DEET, propargite and bromophos.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.

Characterization of Two Species of Trypanosomatidae from the Honey Bee Apis mellifera: Crithidia mellificae Langridge and McGhee, and Lotmaria passim n. gen., n. sp

Trypanosomatids are increasingly recognized as prevalent in European honey bees (Apis mellifera) and by default are attributed to one recognized species, Crithidia mellificae Langridge and McGhee, 1967. We provide reference genetic and ultrastructural data for type isolates of C. mellificae (ATCC 30254 and 30862) in comparison with two recent isolates from A. mellifera (BRL and SF). Phylogenetics unambiguously identify strains BRL/SF as a novel taxonomic unit distinct from C. mellificae strains 30254/30862 and assign all four strains as lineages of a novel clade within the subfamily Leishmaniinae. In vivo analyses show strains BRL/SF preferably colonize the hindgut, lining the lumen as adherent spheroids in a manner identical to previous descriptions from C. mellificae. Microscopy images show motile forms of C. mellificae are distinct from strains BRL/SF. We propose the binomial Lotmaria passim n. gen., n. sp. for this previously undescribed taxon. Analyses of new and previously accessioned genetic data show C. mellificae is still extant in bee populations, however, L. passim n. gen., n. sp. is currently the predominant trypanosomatid in A. mellifera globally. Our findings require that previous reports of C. mellificae be reconsidered and that subsequent trypanosomatid species designations from Hymenoptera provide genetic support.

Vertical transmission of honey bee viruses in a Belgian queen breeding program

Background

The Member States of European Union are encouraged to improve the general conditions for the production and marketing of apicultural products. In Belgium, programmes on the restocking of honey bee hives have run for many years. Overall, the success ratio of this queen breeding programme has been only around 50%. To tackle this low efficacy, we organized sanitary controls of the breeding queens in 2012 and 2014.

Results

We found a high quantity of viruses, with more than 75% of the egg samples being infected with at least one virus. The most abundant viruses were Deformed Wing Virus and Sacbrood Virus (≥40%), although Lake Sinai Virus and Acute Bee Paralysis Virus were also occasionally detected (between 10-30%). In addition, Aphid Lethal Paralysis Virus strain Brookings, Black Queen Cell Virus, Chronic Bee Paralysis Virus and Varroa destructor Macula-like Virus occurred at very low prevalences (≤5%). Remarkably, we found Apis mellifera carnica bees to be less infected with Deformed Wing Virus than Buckfast bees (p < 0.01), and also found them to have a lower average total number of infecting viruses (p < 0.001). This is a significant finding, given that Deformed Wing Virus has earlier been shown to be a contributory factor to winter mortality and Colony Collapse Disorder. Moreover, negative-strand detection of Sacbrood Virus in eggs was demonstrated for the first time.

Conclusions

High pathogen loads were observed in this sanitary control program. We documented for the first time vertical transmission of some viruses, as well as significant differences between two honey bee races in being affected by Deformed Wing Virus. Nevertheless, we could not demonstrate a correlation between the presence of viruses and queen breeding efficacies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-015-0386-9) contains supplementary material, which is available to authorized users.

Genome sequence heterogeneity of Lake Sinai Virus found in honey bees and Orf1/RdRP-based polymorphisms in a single host

Honey bees (Apis mellifera) are susceptible to a wide range of pathogens, including a broad set of viruses. Recently, next-generation sequencing has expanded the list of viruses with, for instance, two strains of Lake Sinai Virus. Soon after its discovery in the USA, LSV was also discovered in other countries and in other hosts. In the present study, we assemble four almost complete LSV genomes, and show that there is remarkable sequence heterogeneity based on the Orf1, RNA-dependent RNA polymerase and capsid protein sequences in comparison to the previously identified LSV 1 and 2 strains. Phylogenetic analyses of LSV sequences obtained from single honey bee specimens further revealed that up to three distinctive clades could be present in a single bee. Such superinfections have not previously been identified for other honey bee viruses. In a search for the putative routes of LSV transmission, we were able to demonstrate the presence of LSV in pollen pellets and in Varroa destructor mites. However, negative-strand analyses demonstrated that the virus only actively replicates in honey bees and mason bees (Osmia cornuta) and not in Varroa mites.

Exploring the hidden honeybee (Apis mellifera) venom proteome by integrating a combinatorial peptide ligand library approach with FTMS

At present, 30 compounds have been described in the venom of the honeybee, and 16 of them were confirmed by mass spectrometry. Previous studies typically combined 2-D PAGE with MALDI-TOF/TOF MS, a technology which now appears to lack sensitivity to detect additional venom compounds. Here, we report an in-depth study of the honeybee venom proteome using a combinatorial peptide ligand library sample pretreatment to enrich for minor components followed by shotgun LC-FT-ICR MS analysis. This strategy revealed an unexpectedly rich venom composition: in total 102 proteins and peptides were found, with 83 of them never described in bee venom samples before. Based on their predicted function and subcellular location, the proteins could be divided into two groups. A group of 33 putative toxins is proposed to contribute to venom activity by exerting toxic functions or by playing a role in social immunity. The other group, considered as venom trace molecules, appears to be secreted for their functions in the extracellular space, or is unintentionally secreted by the venom gland cells due to insufficient protein recycling or co-secretion with other compounds. In conclusion, our approach allowed to explore the hidden honeybee venom proteome and extended the list of potential venom allergens.

BIOLOGICAL SIGNIFICANCE

This study dug deeper into the complex honeybee venom proteome than ever before by applying a highly performing sample pretreatment and mass spectrometric technology. We present putative biological functions for all identified compounds, largely extending our knowledge of the venom toxicity. In addition, this study offers a long list of potential new venom allergens.

IgE recognition of chimeric isoforms of the honeybee (Apis mellifera) venom allergen Api m 10 evaluated by protein array technology

Api m 10 has recently been established as novel major allergen that is recognized by more than 60% of honeybee venom (HBV) allergic patients. Previous studies suggest Api m 10 protein heterogeneity which may have implications for diagnosis and immunotherapy of HBV allergy. In the present study, RT-PCR revealed the expression of at least nine additional Api m 10 transcript isoforms by the venom glands. Two distinct mechanisms are responsible for the generation of these isoforms: while the previously known variant 2 is produced by an alternative splicing event, novel identified isoforms are intragenic chimeric transcripts. To the best of our knowledge, this is the first report of the identification of chimeric transcripts generated by the honeybee. By a retrospective proteomic analysis we found evidence for the presence of several of these isoforms in the venom proteome. Additionally, we analyzed IgE reactivity to different isoforms by protein array technology using sera from HBV allergic patients, which revealed that IgE recognition of Api m 10 is both isoform- and patient-specific. While it was previously demonstrated that the majority of HBV allergic patients display IgE reactivity to variant 2, our study also shows that some patients lacking IgE antibodies for variant 2 display IgE reactivity to two of the novel identified Api m 10 variants, i.e. variants 3 and 4.

Holistic screening of collapsing honey bee colonies in Spain: a case study

Background

Here we present a holistic screening of collapsing colonies from three professional apiaries in Spain. Colonies with typical honey bee depopulation symptoms were selected for multiple possible factors to reveal the causes of collapse.

Results

Omnipresent were Nosema ceranae and Lake Sinai Virus. Moderate prevalences were found for Black Queen Cell Virus and trypanosomatids, whereas Deformed Wing Virus, Aphid Lethal Paralysis Virus strain Brookings and neogregarines were rarely detected. Other viruses, Nosema apis, Acarapis woodi and Varroa destructor were not detected. Palinologic study of pollen demonstrated that all colonies were foraging on wild vegetation. Consequently, the pesticide residue analysis was negative for neonicotinoids. The genetic analysis of trypanosomatids GAPDH gene, showed that there is a large genetic distance between Crithidia mellificae ATCC30254, an authenticated cell strain since 1974, and the rest of the presumed C. mellificae sequences obtained in our study or published. This means that the latter group corresponds to a highly differentiated taxon that should be renamed accordingly.

Conclusion

The results of this study demonstrate that the drivers of colony collapse may differ between geographic regions with different environmental conditions, or with different beekeeping and agricultural practices. The role of other pathogens in colony collapse has to bee studied in future, especially trypanosomatids and neogregarines. Beside their pathological effect on honey bees, classification and taxonomy of these protozoan parasites should also be clarified.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-649) contains supplementary material, which is available to authorized users.

Widespread occurrence of honey bee pathogens in solitary bees

Solitary bees and honey bees from a neighbouring apiary were screened for a broad set of putative pathogens including protists, fungi, spiroplasmas and viruses. Most sampled bees appeared to be infected with multiple parasites. Interestingly, viruses exclusively known from honey bees such as Apis mellifera Filamentous Virus and Varroa destructor Macula-like Virus were also discovered in solitary bees. A microsporidium found in Andrena vaga showed most resemblance to Nosema thomsoni. Our results suggest that bee hives represent a putative source of pathogens for other pollinators. Similarly, solitary bees may act as a reservoir of honey bee pathogens.

Gamma irradiation of pollen and eradication of Israeli acute paralysis virus

Honeybees and bumblebees are the most important pollinators of agricultural crops. For this purpose honeybees and bumblebees are reared and transported. A pathogen-free status of bees in general, is crucial. Indeed anthropogenic transports of hosts carrying parasites could alter the natural host/pathogen association, inducing an extra pathogenic stress. Therefore the creation of a pathogen-free rearing environment is needed. For bumblebees this is possible, as these species are reared in a closed environment. Although, a link remains between reared bumblebees and the outside bee community, as honeybee-collected pollen is essential food for bumblebee mass rearing. Here we evaluated if gamma irradiation can minimize the risk of this potential route of exposure and can inactivate viral particles present in honeybee-collected pollen. We show that 16.9kGy gamma irradiation induced a 100-1000 fold reduction on the ability of IAPV to cause mortality after injections. This result opens avenues toward rearing pathogen-free bumblebees and towards eliminating the risks of pathogen spillover to native wild bee species.

Effect of oral infection with Kashmir bee virus and Israeli acute paralysis virus on bumblebee (Bombus terrestris) reproductive success

Israeli acute paralysis virus (IAPV) together with Acute bee paralysis virus (ABPV) and Kashmir bee virus (KBV) constitute a complex of closely related dicistroviruses. They are infamous for their high mortality after injection in honeybees. These viruses have also been reported in non-Apis hymenopteran pollinators such as bumblebees, which got infected with IAPV when placed in the same greenhouse with IAPV infected honeybee hives. Here we orally infected Bombus terrestris workers with different doses of either IAPV or KBV viral particles. The success of the infection was established by analysis of the bumblebees after the impact studies: 50days after infection. Doses of 0.5×10(7) and 1×10(7) virus particles per bee were infectious over this period, for IAPV and KBV respectively, while a dose of 0.5×10(6) IAPV particles per bee was not infectious. The impact of virus infection was studied in micro-colonies consisting of 5 bumblebees, one of which becomes a pseudo-queen which proceeds to lay unfertilized (drone) eggs. The impact parameters studied were: the establishment of a laying pseudo-queen, the timing of egg-laying, the number of drones produced, the weight of these drones and worker mortality. In this setup KBV infection resulted in a significant slower colony startup and offspring production, while only the latter can be reported for IAPV. Neither virus increased worker mortality, at the oral doses used. We recommend further studies on how these viruses transmit between different pollinator species. It is also vital to understand how viral prevalence can affect wild bee populations because disturbance of the natural host-virus association may deteriorate the already critically endangered status of many bumblebee species.