Metagenomic Approach with the NetoVIR Enrichment Protocol Reveals Virus Diversity within Ethiopian Honey Bees (Apis mellifera simensis).

Metagenomics studies have accelerated the discovery of novel or divergent viruses of the honey bee. However, most of these studies predominantly focused on RNA viruses, and many suffer from the relatively low abundance of viral nucleic acids in the samples (i.e., compared to that of the host). Here, we explored the virome of the Ethiopian honey bee, Apis mellifera simensis, using an unbiased metagenomic approach in which the next-generation sequencing step was preceded by an enrichment protocol for viral particles. Our study revealed five well-known bee viruses and 25 atypical virus species, most of which have never been found in A. mellifera before. The viruses belong to Iflaviridae, Dicistroviridae, Secoviridae, Partitiviridae, Parvoviridae, Potyviridae, and taxonomically unclassified families. Fifteen of these atypical viruses were most likely plant-specific, and the remaining ten were presumed to be insect-specific. Apis mellifera filamentous virus (AmFV) was found in one sampling site out of 10. Two samples contained high read counts of a virus similar to Diatraea saccharales densovirus (DsDNV), which is a virus that causes high mortality in the sugarcane borer. AmFV and the DsDNV-like virus were the only DNA viruses found. Three viruses that primarily infect Drosophila spp. were also discovered: La Jolla virus (LJV), Kilifi virus (KiV), and Thika virus. Our study suggests that phoretic varroa mites are involved in the transmission of LJV and KiV and that both viruses replicate in mites and adult bees. We also found an overwhelming dominance of the deformed wing virus type B variant, which fits well with the apparently harmless infestation by Varroa destructor. It was suggested that Ethiopian bees have developed tolerance against virus infections as the result of natural selection.

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.

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.

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.

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.