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.

Honey-bee-associated prokaryotic viral communities reveal wide viral diversity and a profound metabolic coding potential.

Honey bees (Apis mellifera) produce an enormous economic value through their pollination activities and play a central role in the biodiversity of entire ecosystems. Recent efforts have revealed the substantial influence that the gut microbiota exert on bee development, food digestion, and homeostasis in general. In this study, deep sequencing was used to characterize prokaryotic viral communities associated with honey bees, which was a blind spot in research up until now. The vast majority of the prokaryotic viral populations are novel at the genus level, and most of the encoded proteins comprise unknown functions. Nevertheless, genomes of bacteriophages were predicted to infect nearly every major bee-gut bacterium, and functional annotation and auxiliary metabolic gene discovery imply the potential to influence microbial metabolism. Furthermore, undiscovered genes involved in the synthesis of secondary metabolic biosynthetic gene clusters reflect a wealth of previously untapped enzymatic resources hidden in the bee bacteriophage community.

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.

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.

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 LC50 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.

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.

Varroa destructor Macula-like virus, Lake Sinai virus and other new RNA viruses in wild bumblebee hosts (Bombus pascuorum, Bombus lapidarius and Bombus pratorum).

Pollinators such as bumblebees (Bombus spp.) are in decline worldwide which poses a threat not only for ecosystem biodiversity but also to human crop production services. One main cause of pollinator decline may be the infection and transmission of diseases including RNA viruses. Recently, new viruses have been discovered in honeybees, but information on the presence of these in wild bumblebees is largely not available. In this study, we investigated the prevalence of new RNA viruses in Bombus species, and can report for the first time Varroa destructor Macula-like virus (VdMLV) and Lake Sinai virus (LSV) infection in multiple wild bumblebee hosts of Bombus pascuorum, Bombus lapidarius and Bombus pratorum. We sampled in 4 locations in Flanders, Belgium. Besides, we confirmed Slow bee paralysis virus (SBPV) in wild bumblebees, but no positive samples were obtained for Big Sioux river virus (BSRV). Secondly, we screened for the influence of apiaries on the prevalence of these viruses. Our results indicated a location effect for the prevalence of VdMLV in Bombus species, with a higher prevalence in the proximity of honeybee apiaries mainly observed in one location. For LSV, the prevalence was not different in the proximity or at a 1.5 km-distance of apiaries, but we reported a different isolate with similarities to LSV-2 and "LSV-clade A" as described by Ravoet et al. (2015), which was detected both in Apis mellifera and Bombus species. In general, our results indicate the existence of a disease pool of new viruses that seems to be associated to a broad range of Apoidae hosts, including multiple Bombus species.

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.

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.

Bee pathogens found in Bombus atratus from Colombia: A case study.

Bombus atratus bumblebees from Colombia that were caught in the wild and from breeding programs were screened for a broad set of bee pathogens. We discovered for the first time Lake Sinai Virus and confirmed the infection by other common viruses. The prevalence of Apicystis bombi, Crithidia bombi and Nosema ceranae was remarkably high. According to other studies the former two could have been co-introduced in South America with exotic bumble bees as Bombus terrestris or Bombus ruderatus. Given the fact that none of these species occur in Colombia, our data puts a new light on the spread of these pathogens over the South American continent.

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.

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.

Finding the missing honey bee genes: lessons learned from a genome upgrade.

BACKGROUND: The first generation of genome sequence assemblies and annotations have had a significant impact upon our understanding of the biology of the sequenced species, the phylogenetic relationships among species, the study of populations within and across species, and have informed the biology of humans. As only a few Metazoan genomes are approaching finished quality (human, mouse, fly and worm), there is room for improvement of most genome assemblies. The honey bee (Apis mellifera) genome, published in 2006, was noted for its bimodal GC content distribution that affected the quality of the assembly in some regions and for fewer genes in the initial gene set (OGSv1.0) compared to what would be expected based on other sequenced insect genomes. RESULTS: Here, we report an improved honey bee genome assembly (Amel_4.5) with a new gene annotation set (OGSv3.2), and show that the honey bee genome contains a number of genes similar to that of other insect genomes, contrary to what was suggested in OGSv1.0. The new genome assembly is more contiguous and complete and the new gene set includes ~5000 more protein-coding genes, 50% more than previously reported. About 1/6 of the additional genes were due to improvements to the assembly, and the remaining were inferred based on new RNAseq and protein data. CONCLUSIONS: Lessons learned from this genome upgrade have important implications for future genome sequencing projects. Furthermore, the improvements significantly enhance genomic resources for the honey bee, a key model for social behavior and essential to global ecology through pollination.

Reproduction of honeybee workers is regulated by epidermal growth factor receptor signaling.

Eusocial insect societies display a remarkable reproductive division of labor between a single fertile queen and thousands of largely sterile workers. In most species, however, the workers retain the capacity to reproduce, particularly in queenless colonies where typically many workers lay eggs. As yet, the molecular determinants that initiate this shift in worker fertility are still poorly documented. By using RNA interference we here demonstrate that the knockdown of epidermal growth factor receptor, a gene which was previously shown to be involved in queen-worker caste differentiation, also induces reproduction in worker honeybees (Apis mellifera). These data show that worker fertility and queen-worker caste determination partly rely on the same gene regulatory networks, thereby providing a major breakthrough in our understanding of the molecular determinants of the social insects' spectacular reproductive division of labor.

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.

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.

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.

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.

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.