The Apis mellifera Filamentous Virus Genome

Ontogeny of immunity and natural viral infection in Apis mellifera drones and workers

The most common viral diseases affecting honey bees (Apis mellifera) in Israel include deformed wing viruses (DWV A and DWV-B) and acute paralysis viruses (ABPV and IAPV). These viruses are transmitted within and between colonies, both horizontally and vertically. All members of the colony contribute to this transmission, on the other hand individual and social immunity, particularly hygienic behaviour, may affect the outcome of the process. In this study, we evaluated the ontogeny of natural infections of DWV-A, DWV-B, ABPV and IAPV, their prevalence and loads, in workers and drones from high (H) and low (L) hygienic colonies. In parallel, we evaluated the expression of two immune genes: peptidoglycan recognition protein S2(PGRP-S2) and hymenoptaecin. The prevalence of DWV-B and IAPV increased with age and was higher in workers than in drones. ABPV was not detected in drones. The expression of both immune genes was significantly affected by age and sex. Drones from H colonies had higher expression of these genes. The increased expression of immune genes with drones' age, particularly in hygienic colonies, suggest additional value of honey bee breeding for hygienic behaviour for sustainable beekeeping.

Discovery and characterization of novel DNA viruses in Apis mellifera: expanding the honey bee virome through metagenomic analysis

To date, many viruses have been discovered to infect honey bees. In this study, we used high-throughput sequencing to expand the known virome of the honey bee, Apis mellifera, by identifying several novel DNA viruses. While the majority of previously identified bee viruses are RNA, our study reveals nine new genomes from the Parvoviridae family, tentatively named Bee densoviruses 1 to 9. In addition, we characterized a large DNA virus, Apis mellifera filamentous-like virus (AmFLV), which shares limited protein identities with the known Apis mellifera filamentous virus. The complete sequence of AmFLV, obtained by a combination of laboratory techniques and bioinformatics, spans 152,678 bp. Linear dsDNA genome encodes for 112 proteins, of which 49 are annotated. Another large virus we discovered is Apis mellifera nudivirus, which belongs to a group of Alphanudivirus. The virus has a length of 129,467 bp and a circular dsDNA genome, and has 106 protein encoding genes. The virus contains most of the core genes of the family Nudiviridae. This research demonstrates the effectiveness of viral binning in identifying viruses in honey bee virology, showcasing its initial application in this field.IMPORTANCEHoney bees contribute significantly to food security by providing pollination services. Understanding the virome of honey bees is crucial for the health and conservation of bee populations and also for the stability of the ecosystems and economies for which they are indispensable. This study unveils previously unknown DNA viruses in the honey bee virome, expanding our knowledge of potential threats to bee health. The use of the viral binning approach we employed in this study offers a promising method to uncovering and understanding the vast viral diversity in these essential pollinators.

Apis mellifera filamentous virus from a honey bee gut microbiome survey in Hungary

In Hungary, as part of a nationwide, climatically balanced survey for a next-generation sequencing-based study of the honey bee (Apis mellifera) gut microbiome, repeated sampling was carried out during the honey production season (March and May 2019). Among other findings, the presence of Apis mellifera filamentous virus (AmFV) was detected in all samples, some at very high levels. AmFV-derived reads were more abundant in the March samples than in the May samples. In March, a higher abundance of AmFV-originated reads was identified in samples collected from warmer areas compared to those collected from cooler areas. A lower proportion of AmFV-derived reads were identified in samples collected in March from the wetter areas than those collected from the drier areas. AmFV-read abundance in samples collected in May showed no significant differences between groups based on either environmental temperature or precipitation. The AmFV abundance correlated negatively with Bartonella apihabitans, Bartonella choladocola, and positively with Frischella perrara, Gilliamella apicola, Gilliamella sp. ESL0443, Lactobacillus apis, Lactobacillus kullabergensis, Lactobacillus sp. IBH004. De novo metagenome assembly of four samples resulted in almost the complete AmFV genome. According to phylogenetic analysis based on DNA polymerase, the Hungarian strains are closest to the strain CH-05 isolated in Switzerland.

A novel and diverse family of filamentous DNA viruses associated with parasitic wasps

Large dsDNA viruses from the Naldaviricetes class are currently composed of four viral families infecting insects and/or crustaceans. Since the 1970s, particles described as filamentous viruses (FVs) have been observed by electronic microscopy in several species of Hymenoptera parasitoids but until recently, no genomic data was available. This study provides the first comparative morphological and genomic analysis of these FVs. We analyzed the genomes of seven FVs, six of which were newly obtained, to gain a better understanding of their evolutionary history. We show that these FVs share all genomic features of the Naldaviricetes while encoding five specific core genes that distinguish them from their closest relatives, the Hytrosaviruses. By mining public databases, we show that FVs preferentially infect Hymenoptera with parasitoid lifestyle and that these viruses have been repeatedly integrated into the genome of many insects, particularly Hymenoptera parasitoids, overall suggesting a long-standing specialization of these viruses to parasitic wasps. Finally, we propose a taxonomical revision of the class Naldaviricetes in which FVs related to the Leptopilina boulardi FV constitute a fifth family. We propose to name this new family, Filamentoviridae.

Integrase-associated niche differentiation of endogenous large DNA viruses in crustaceans

IMPORTANCE

Crustacean genomes harbor sequences originating from a family of large DNA viruses called nimaviruses, but it is unclear why they are present. We show that endogenous nimaviruses selectively insert into repetitive sequences within the host genome, and this insertion specificity was correlated with different types of integrases, which are DNA recombination enzymes encoded by the nimaviruses themselves. This suggests that endogenous nimaviruses have colonized various genomic niches through the acquisition of integrases with different insertion specificities. Our results point to a novel survival strategy of endogenous large DNA viruses colonizing the host genomes. These findings may clarify the evolution and spread of nimaviruses in crustaceans and lead to measures to control and prevent the spread of pathogenic nimaviruses in aquaculture settings.

Virus replication in the honey bee parasite, Varroa destructor

IMPORTANCE

The parasitic mite Varroa destructor is a significant driver of worldwide colony losses of our most important commercial pollinator, the Western honey bee Apis mellifera. Declines in honey bee health are frequently attributed to the viruses that mites vector to honey bees, yet whether mites passively transmit viruses as a mechanical vector or actively participate in viral amplification and facilitate replication of honey bee viruses is debated. Our work investigating the antiviral RNA interference response in V. destructor demonstrates that key viruses associated with honey bee declines actively replicate in mites, indicating that they are biological vectors, and the host range of bee-associated viruses extends to their parasites, which could impact virus evolution, pathogenicity, and spread.

Data mining reveals tissue-specific expression and host lineage-associated forms of Apis mellifera filamentous virus

Background

Apis mellifera filamentous virus (AmFV) is a large double-stranded DNA virus of uncertain phylogenetic position that infects honey bees (Apis mellifera). Little is known about AmFV evolution or molecular aspects of infection. Accurate annotation of open-reading frames (ORFs) is challenged by weak homology to other known viruses. This study was undertaken to evaluate ORFs (including coding-frame conservation, codon bias, and purifying selection), quantify genetic variation within AmFV, identify host characteristics that covary with infection rate, and examine viral expression patterns in different tissues.

Methods

Short-read data were accessed from the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI). Sequence reads were downloaded from accessions meeting search criteria and scanned for kmers representative of AmFV genomic sequence. Samples with kmer counts above specified thresholds were downloaded in full for mapping to reference sequences and de novo assembly.

Results

At least three distinct evolutionary lineages of AmFV exist. Clade 1 predominates in Europe but in the Americas and Africa it is replaced by the other clades as infection level increases in hosts. Only clade 3 was found at high relative abundance in hosts with African ancestry, whereas all clades achieved high relative abundance in bees of non-African ancestry. In Europe and Africa, clade 2 was generally detected only in low-level infections but was locally dominant in some North American samples. The geographic distribution of clade 3 was consistent with an introduction to the Americas with 'Africanized' honey bees in the 1950s. Localized genomic regions of very high nucleotide divergence in individual isolates suggest recombination with additional, as-yet unidentified AmFV lineages. A set of 155 high-confidence ORFs was annotated based on evolutionary conservation in six AmFV genome sequences representative of the three clades. Pairwise protein-level identity averaged 94.6% across ORFs (range 77.1-100%), which generally exhibited low evolutionary rates and moderate to strong codon bias. However, no robust example of positive diversifying selection on coding sequence was found in these alignments. Most of the genome was detected in RNA short-read alignments. Transcriptome assembly often yielded contigs in excess of 50 kb and containing ORFs in both orientations, and the termini of long transcripts were associated with tandem repeats. Lower levels of AmFV RNA were detected in brain tissue compared to abdominal tissue, and a distinct set of ORFs had minimal to no detectable expression in brain tissue. A scan of DNA accessions from the parasitic mite Varroa destructor was inconclusive with respect to replication in that species.

Discussion

Collectively, these results expand our understanding of this enigmatic virus, revealing transcriptional complexity and co-evolutionary associations with host lineage.

Developments in the classification and nomenclature of arthropod-infecting large DNA viruses that contain pif genes

Viruses of four families of arthropod-specific, large dsDNA viruses (the nuclear arthropod large DNA viruses, or NALDVs) possess homologs of genes encoding conserved components involved in the baculovirus primary infection mechanism. The presence of such homologs encoding per os infectivity factors (pif genes), along with their absence from other viruses and the occurrence of other shared characteristics, suggests a common origin for the viruses of these families. Therefore, the class Naldaviricetes was recently established, accommodating these four families. In addition, within this class, the ICTV approved the creation of the order Lefavirales for three of these families, whose members carry homologs of the baculovirus genes that code for components of the viral RNA polymerase, which is responsible for late gene expression. We further established a system for the binomial naming of all virus species in the order Lefavirales, in accordance with a decision by the ICTV in 2019 to move towards a standardized nomenclature for all virus species. The binomial species names for members of the order Lefavirales consist of the name of the genus to which the species belongs (e.g., Alphabaculovirus), followed by a single epithet that refers to the host species from which the virus was originally isolated. The common names of viruses and the abbreviations thereof will not change, as the format of virus names lies outside the remit of the ICTV.

Metagenomic analysis of viromes in honey bee colonies (Apis mellifera; Hymenoptera: Apidae) after mass disappearance in Korea

After the nationwide, massive winter losses of honey bees in Korea during the winter of 2021, samplings were conducted from live honey bees in colonies and dead honey bees nearby colonies in the same bee-farms in six regions in Korea. Each sample was subjected to virome analysis using high-throughput sequencing technology. The number of viral reads was the lowest in the live honey bee group sample with 370,503 reads and the highest in the dead honey bee group sample with 42,659,622 reads. Viral contigs were matched with the viral genomes of the black queen cell virus, deformed wing virus, Israeli acute paralysis virus, and sacbrood virus, all of which have been previously reported in Korea. However, Apis rhabdovirus 5, bee macula-like virus, Varroa orthomyxovirus-1, Hubei partiti-like virus 34, Lake Sinai virus 2, 3, and 4, and the Ditton virus, were also discovered in this study, which are the first records in Korea. Plant viral sequences resembling those of Arabidopsis latent virus 1, and a novel viral sequence was also discovered. In the present study 55 complete viral genome sequences were identified. This study is the first virome analysis of domestic honey bees and provides the latest information on the diversity of honey bee viruses in Korea.

Emerging Risk of Cross-Species Transmission of Honey Bee Viruses in the Presence of Invasive Vespid Species

The increase in invasive alien species is a concern for the environment. The establishment of some of these species may be changing the balance between pathogenicity and host factors, which could alter the defense strategies of native host species. Vespid species are among the most successful invasive animals, such as the genera Vespa, Vespula and Polistes. Bee viruses have been extensively studied as an important cause of honey bee population losses. However, knowledge about the transmission of honey bee viruses in Vespids is a relevant and under-researched aspect. The role of some mites such as Varroa in the transmission of honey bee viruses is clearer than in the case of Vespidae. This type of transmission by vectors has not yet been clarified in Vespidae, with interspecific relationships being the main hypotheses accepted for the transmission of bee viruses. A majority of studies describe the presence of viruses or their replicability, but aspects such as the symptomatology in Vespids or the ability to infect other hosts from Vespids are scarcely discussed. Highlighting the case of Vespa velutina as an invader, which is causing huge losses in European beekeeping, is of special interest. The pressure caused by V. velutina leads to weakened hives that become susceptible to pathogens. Gathering this information is necessary to promote further research on the spread of bee viruses in ecosystems invaded by invasive species of Vespids, as well as to prevent the decline of bee populations due to bee viruses.

Functional Peroral Infectivity Complex of White Spot Syndrome Virus of Shrimp

White spot syndrome virus (WSSV) is a major cause of disease in shrimp cultures worldwide. The infection process of this large circular double-stranded DNA virus has been well studied, but its entry mechanism remains controversial. The major virion envelope protein VP28 has been implicated in oral and systemic viral infection in shrimp. However, genetic analysis of viral DNA has shown the presence of a few genes related to proteins of per os infectivity factor (PIF) complex in baculoviruses. This complex is essential for the entry of baculoviruses, large terrestrial circular DNA viruses, into the midgut epithelial cells of insect larvae. In this study, we aimed to determine whether a PIF complex exists in WSSV, the components of this complex, whether it functions as an oral infectivity complex in shrimp, and the biochemical properties that contribute to its function in a marine environment. The results revealed a WSSV PIF complex (~720 kDa) comprising at least eight proteins, four of which were not identified as PIF homologs: WSV134, VP124 (WSV216), WSSV021, and WSV136. WSV134 is suggested to be a PIF4 homolog due to predicted structural similarity and amino acid sequence identity. The WSSV PIF complex is resistant to alkali, proteolysis, and high salt, properties that are important for maintaining infectivity in aquatic environments. Oral infection can be neutralized by PIF-specific antibodies but not by VP28-specific antibodies. These results indicate that the WSSV PIF complex is critical for WSSV entry into shrimp; the complex's evolutionary significance is also discussed. IMPORTANCE White spot disease, caused by the white spot syndrome virus (WSSV), is a major scourge in cultured shrimp production facilities worldwide. This disease is only effectively controlled by sanitation. Intervention strategies are urgently needed but are limited by a lack of appropriate targets. Our identification of a per os infectivity factor (PIF) complex, which is pivotal for the entry of WSSV into shrimp, could provide new targets for antibody- or dsRNA-based intervention strategies. In addition, the presence of a PIF complex with at least eight components in WSSV, which is ancestrally related to the PIF complex of invertebrate baculoviruses, suggests that this complex is structurally and functionally conserved in disparate virus taxa.

Reconstructing the ecosystem context of a species: Honey-borne DNA reveals the roles of the honeybee

To assess a species’ impact on its environment–and the environment’s impact upon a species–we need to pinpoint its links to surrounding taxa. The honeybee (Apis mellifera) provides a promising model system for such an exercise. While pollination is an important ecosystem service, recent studies suggest that honeybees can also provide disservices. Developing a comprehensive understanding of the full suite of services and disservices that honeybees provide is a key priority for such a ubiquitous species. In this perspective paper, we propose that the DNA contents of honey can be used to establish the honeybee’s functional niche, as reflected by ecosystem services and disservices. Drawing upon previously published genomic data, we analysed the DNA found within 43 honey samples from Northern Europe. Based on metagenomic analysis, we find that the taxonomic composition of DNA is dominated by a low pathogenicity bee virus with 40.2% of the reads, followed by bacteria (16.7%), plants (9.4%) and only 1.1% from fungi. In terms of ecological roles of taxa associated with the bees or taxa in their environment, bee gut microbes dominate the honey DNA, with plants as the second most abundant group. A range of pathogens associated with plants, bees and other animals occur frequently, but with lower relative read abundance, across the samples. The associations found here reflect a versatile the honeybee’s role in the North-European ecosystem. Feeding on nectar and pollen, the honeybee interacts with plants–in particular with cultivated crops. In doing so, the honeybee appears to disperse common pathogens of plants, pollinators and other animals, but also microbes potentially protective of these pathogens. Thus, honey-borne DNA helps us define the honeybee’s functional niche, offering directions to expound the benefits and drawbacks of the associations to the honeybee itself and its interacting organisms.

Viral communities in the parasite Varroa destructor and in colonies of their honey bee host (Apis mellifera) in New Zealand

The parasitic mite Varroa destructor is a leading cause of mortality for Western honey bee (Apis mellifera) colonies around the globe. We sought to confirm the presence and likely introduction of only one V. destructor haplotype in New Zealand, and describe the viral community within both V. destructor mites and the bees that they parasitise. A 1232 bp fragment from mitochondrial gene regions suggests the likely introduction of only one V. destructor haplotype to New Zealand. Seventeen viruses were found in bees. The most prevalent and abundant was the Deformed wing virus A (DWV-A) strain, which explained 95.0% of the variation in the viral community of bees. Black queen cell virus, Sacbrood virus, and Varroa destructor virus 2 (VDV-2) played secondary roles. DWV-B and the Israeli acute paralysis virus appeared absent from New Zealand. Ten viruses were observed in V. destructor, with > 99.9% of viral reads from DWV-A and VDV-2. Substantially more variation in viral loads was observed in bees compared to mites. Where high levels of VDV-2 occurred in mites, reduced DWV-A occurred in both the mites and the bees co-occurring within the same hive. Where there were high loads of DWV-A in mites, there were typically high viral loads in bees.

Genomics and Proteomics of Apis mellifera Filamentous Virus Isolated from Honeybees in China

Apis mellifera filamentous virus (AmFV) is a large DNA virus that is endemic in honeybee colonies. The genome sequence of the AmFV Swiss isolate (AmFV CH–C05) has been reported, but so far very few molecular studies have been conducted on this virus. In this study, we isolated and purified AmFV (AmFV CN) from Chinese honeybee (Apis mellifera) colonies and elucidated its genomics and proteomics. Electron microscopy showed ovoid purified virions with dimensions of 300–500 ​× ​210–285 ​nm, wrapping a 3165 ​× ​40 ​nm filamentous nucleocapsid in three figure-eight loops. Unlike AmFV CH–C05, which was reported to have a circular genome, our data suggest that AmFV CN has a linear genome of approximately 493 ​kb. A total of 197 ORFs were identified, among which 36 putative genes including 18 baculoviral homologs were annotated. The overall nucleotide similarity between the CN and CH–C05 isolates was 96.9%. Several ORFs were newly annotated in AmFV CN, including homologs of per os infectivity factor 4 (PIF4) and a putative integrase. Phylogenomic analysis placed AmFVs on a separate branch within the newly proposed virus class Naldaviricetes. Proteomic analysis revealed 47 AmFV virion-associated proteins, of which 14 had over 50% sequence coverage, suggesting that they are likely to be main structural proteins. In addition, all six of the annotated PIFs (PIF-0–5) were identified by proteomics, suggesting that they may function as entry factors in AmFV infection. This study provides fundamental information regarding the molecular biology of AmFV.

•

The AmFV CN contains a 493 ​kb linear genome encoding 197 ORFs.

•

Proteomics revealed 14 putative major structural proteins.

•

AmFV belongs to the class Naldaviricetes but not the order Lefavirales.

Transitions in symbiosis: evidence for environmental acquisition and social transmission within a clade of heritable symbionts

A dynamic continuum exists from free-living environmental microbes to strict host-associated symbionts that are vertically inherited. However, knowledge of the forces that drive transitions in symbiotic lifestyle and transmission mode is lacking. Arsenophonus is a diverse clade of bacterial symbionts, comprising reproductive parasites to coevolving obligate mutualists, in which the predominant mode of transmission is vertical. We describe a symbiosis between a member of the genus Arsenophonus and the Western honey bee. The symbiont shares common genomic and predicted metabolic properties with the male-killing symbiont Arsenophonus nasoniae, however we present multiple lines of evidence that the bee Arsenophonus deviates from a heritable model of transmission. Field sampling uncovered spatial and seasonal dynamics in symbiont prevalence, and rapid infection loss events were observed in field colonies and laboratory individuals. Fluorescent in situ hybridisation showed Arsenophonus localised in the gut, and detection was rare in screens of early honey bee life stages. We directly show horizontal transmission of Arsenophonus between bees under varying social conditions. We conclude that honey bees acquire Arsenophonus through a combination of environmental exposure and social contacts. These findings uncover a key link in the Arsenophonus clades trajectory from free-living ancestral life to obligate mutualism, and provide a foundation for studying transitions in symbiotic lifestyle.

Pathogens Spillover from Honey Bees to Other Arthropods

Honey bees, and pollinators in general, play a major role in the health of ecosystems. There is a consensus about the steady decrease in pollinator populations, which raises global ecological concern. Several drivers are implicated in this threat. Among them, honey bee pathogens are transmitted to other arthropods populations, including wild and managed pollinators. The western honey bee, Apis mellifera, is quasi-globally spread. This successful species acted as and, in some cases, became a maintenance host for pathogens. This systematic review collects and summarizes spillover cases having in common Apis mellifera as the mainteinance host and some of its pathogens. The reports are grouped by final host species and condition, year, and geographic area of detection and the co-occurrence in the same host. A total of eighty-one articles in the time frame 1960-2021 were included. The reported spillover cases cover a wide range of hymenopteran host species, generally living in close contact with or sharing the same environmental resources as the honey bees. They also involve non-hymenopteran arthropods, like spiders and roaches, which are either likely or unlikely to live in close proximity to honey bees. Specific studies should consider host-dependent pathogen modifications and effects on involved host species. Both the plasticity of bee pathogens and the ecological consequences of spillover suggest a holistic approach to bee health and the implementation of a One Health approach.

An Investigation of Honey Bee Viruses Prevalence in Managed Honey Bees (Apis mellifera and Apis cerana) Undergone Colony Decline

Objective:

In the absence of known clinical symptoms, viruses were considered to be the most probable key pathogens of honey bee. Therefore, the aim of this study was to investigate the prevalence and distribution of honey bee viruses in managed Apis mellifera and Apis cerana in China.

Methods:

We conducted a screening of 8 honey bee viruses on A. mellifera and A. cerana samples collected from 54 apiaries from 13 provinces in China using RT-PCR.

Results:

We found that the types and numbers of viral species significantly differed between A. mellifera and A. cerana. Black Queen Cell Virus (BQCV), Chronic Bee Paralysis Virus (CBPV), Apis mellifera filamentous virus (AmFV), and Kakugo virus (DWV-A/KV) were the primary viruses found in A. mellifera colonies, whereas Chinese Sacbrood Bee Virus (CSBV) and Sacbrood Bee Virus (SBV) were the primary viruses found in A. cerana. The percentage infection of BQCV and CSBV were 84.6% and 61.6% in all detected samples. We first detected the occurrences of Varroa destructor virus-1 (VDV-1 or DWV-B) and DWV-A/KV in China but not ABPV in both A. mellifera and A. cerana.

Conclusion:

This study showed that BQCV and CSBV are the major threat to investigated A. mellifera and A. cerana colonies.

Pathogens Detection in the Small Hive Beetle (Aethina tumida (Coleoptera: Nitidulidae))

Aethina tumida Murray is currently a worldwide emergent pest of Apis mellifera L. hives. Although the damaging effect on the colony stores and brood is well known, the possible role of these beetles as a disease carrier is not clear. This is the first report of DNA presence of the trypanosome honeybee parasite Lotmaria passim and Crithidia bombi, and the Apis mellifera filamentous virus (AmFV) in A. tumida. Further studies will be needed to determine if A. tumida is indeed a mechanical or biological vector of these pathogens.

Viruses that affect Argentinian honey bees (Apis mellifera)

Beekeeping is a widespread activity in Argentina, mainly producing honey that has gained both national and international recognition. There are more than 3,000,000 hives in the country, mainly concentrated in Buenos Aires Province (approximately 1,000,000 hives). In recent decades, worrying rates of hive loss have been observed in many countries around the world. In Latin America, the estimated loss of hives is between 13% (Peru and Ecuador) and 53% (Chile). Argentina had annual losses of 34% for the period of October 1, 2016 to October 1, 2017. The causes of these losses are not clear but probably involve multiple stressors that can act simultaneously. One of the main causes of loss of bee colonies worldwide is infestation by the ectoparasitic mite Varroa destructor in combination with viral infections. To date, 10 viruses have been detected that affect honey bees (Apis mellifera) in Argentina. Of these, deformed wing virus, sacbrood virus, acute bee paralysis virus, chronic bee paralysis virus, and Israeli acute bee paralysis can be transmitted by mites. Deformed wing virus and the AIK complex are the viruses most often associated with loss of hives worldwide. Considering that bee viruses have been detected in Argentina in several hymenopteran and non-hymenopteran insects, these hosts could act as important natural reservoirs for viruses and play an important role in their dispersal in the environment. Further studies to investigate the different mechanisms by which viruses spread in the environment will enable us to develop various strategies for the control of infected colonies and the spread of viruses in the habitat where they are found.

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

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

Shotgun sequencing of honey DNA can describe honey bee derived environmental signatures and the honey bee hologenome complexity

Honey bees are large-scale monitoring tools due to their extensive environmental exploration. In their activities and from the hive ecosystem complex, they get in close contact with many organisms whose traces can be transferred into the honey, which can represent an interesting reservoir of environmental DNA (eDNA) signatures and information useful to analyse the honey bee hologenome complexity. In this study, we tested a deep shotgun sequencing approach of honey DNA coupled with a specifically adapted bioinformatic pipeline. This methodology was applied to a few honey samples pointing out DNA sequences from 191 organisms spanning different kingdoms or phyla (viruses, bacteria, plants, fungi, protozoans, arthropods, mammals). Bacteria included the largest number of species. These multi-kingdom signatures listed common hive and honey bee gut microorganisms, honey bee pathogens, parasites and pests, which resembled a complex interplay that might provide a general picture of the honey bee pathosphere. Based on the Apis mellifera filamentous virus genome diversity (the most abundant detected DNA source) we obtained information that could define the origin of the honey at the apiary level. Mining Apis mellifera sequences made it possible to identify the honey bee subspecies both at the mitochondrial and nuclear genome levels.

Bee Viruses: Routes of Infection in Hymenoptera

Numerous studies have recently reported on the discovery of bee viruses in different arthropod species and their possible transmission routes, vastly increasing our understanding of these viruses and their distribution. Here, we review the current literature on the recent advances in understanding the transmission of viruses, both on the presence of bee viruses in Apis and non-Apis bee species and on the discovery of previously unknown bee viruses. The natural transmission of bee viruses will be discussed among different bee species and other insects. Finally, the research potential of in vivo (host organisms) and in vitro (cell lines) serial passages of bee viruses is discussed, from the perspective of the host-virus landscape changes and potential transmission routes for emerging bee virus infections.

Contrasting impacts of a novel specialist vector on multihost viral pathogen epidemiology in wild and managed bees

Typically, pathogens infect multiple host species. Such multihost pathogens can show considerable variation in their degree of infection and transmission specificity, which has important implications for potential disease emergence. Transmission of multihost pathogens can be driven by key host species and changes in such transmission networks can lead to disease emergence. We study two viruses that show contrasting patterns of prevalence and specificity in managed honeybees and wild bumblebees, black queen cell virus (BQCV) and slow bee paralysis virus (SBPV), in the context of the novel transmission route provided by the virus‐vectoring Varroa destructor. Our key result is that viral communities and RNA virus genetic variation are structured by location, not host species or V. destructor presence. Interspecific transmission is pervasive with the same viral variants circulating between pollinator hosts in each location; yet, we found virus‐specific host differences in prevalence and viral load. Importantly, V. destructor presence increases the prevalence in honeybees and, indirectly, in wild bumblebees, but in contrast to its impact on deformed wing virus (DWV), BQCV and SBPV viral loads are not increased by Varroa presence, and do not show genetic evidence of recent emergence. Effective control of Varroa in managed honeybee colonies is necessary to mitigate further disease emergence, and alleviate disease pressure on our vital wild bee populations. More generally, our results highlight the over‐riding importance of geographical location to the epidemiological outcome despite the complexity of multihost‐parasite interactions.

Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids

Since the discovery that honey bee viruses play a role in colony decline, researchers have made major breakthroughs in understanding viral pathology and infection processes in honey bees. Work on virus transmission patterns and virus vectors, such as the mite Varroa destructor, has prompted intense efforts to manage honey bee health. However, little is known about the occurrence of honey bee viruses in bee predators, such as vespids. In this study, we characterized the occurrence of 11 honey bee viruses in five vespid species and one wasp from four provinces in China and two vespid species from four locations in France. The results showed that all the species from China carried certain honey bee viruses, notably Apis mellifera filamentous virus (AmFV), Deformed wing virus (DWV), and Israeli acute paralysis virus (IAPV); furthermore, in some vespid colonies, more than three different viruses were identified. In France, DWV was the most common virus; Sacbrood virus (SBV) and Black queen cell virus (BQCV) were observed in one and two samples, respectively. Phylogenetic analyses of IAPV and BQCV sequences indicated that most of the IAPV sequences belonged to a single group, while the BQCV sequences belonged to several groups. Additionally, our study is the first to detect Lake Sinai virus (LSV) in a hornet from China. Our findings can guide further research into the origin and transmission of honey bee viruses in Vespidae, a taxon of ecological, and potentially epidemiological, relevance.

Assessing the Diversity of Endogenous Viruses Throughout Ant Genomes

Endogenous viral elements (EVEs) can play a significant role in the evolution of their hosts and have been identified in animals, plants, and fungi. Additionally, EVEs potentially provide an important snapshot of the evolutionary frequency of viral infection. The purpose of this study is to take a comparative host-centered approach to EVE discovery in ant genomes to better understand the relationship of EVEs to their ant hosts. Using a comprehensive bioinformatic pipeline, we screened all nineteen published ant genomes for EVEs. Once the EVEs were identified, we assessed their phylogenetic relationships to other closely related exogenous viruses. A diverse group of EVEs were discovered in all screened ant host genomes and in many cases are similar to previously identified exogenous viruses. EVEs similar to ssRNA viral proteins are the most common viral lineage throughout the ant hosts, which is potentially due to more chronic infection or more effective endogenization of certain ssRNA viruses in ants. In addition, both EVEs similar to viral glycoproteins and retrovirus-derived proteins are also abundant throughout ant genomes, suggesting their tendency to endogenize. Several of these newly discovered EVEs are found to be potentially functional within the genome. The discovery and analysis of EVEs is essential in beginning to understand viral–ant interactions over evolutionary time.

Multiple Virus Infections in Western Honeybee (Apis mellifera L.) Ejaculate Used for Instrumental Insemination

Instrumental insemination of Apis mellifera L. queens is a widely employed technique used in honeybee breeding that enables the effective control of mating. However, drone semen represents a potential source of honeybee viruses. In this study, 43 semen doses collected from apparently healthy drones, and consequently used in instrumental insemination, were analysed using PCR or RT-PCR to detect the presence of viral genome of 11 honeybee viruses. In 91% of samples, viral infection was detected. The survey revealed genomes of five viruses, namely Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Black queen cell virus (BQCV), Sacbrood virus (SBV), and A. mellifera filamentous virus (AmFV) in 84%, 19%, 14%, 2%, and 67% of samples, respectively. Single infection (30% of samples) as well as multiple infection (61% of samples) of two, three or four pathogens were also evaluated. To the best of our knowledge, this is the first study describing the presence of the BQCV and SBV genome sequence in drone ejaculate. Phylogenetic analysis of BQCV partial helicase gene sequence revealed the high similarity of nucleotide sequence of described Czech strains, which varied from 91.4% to 99.6%. The findings of our study indicate the possibility of venereal transmission of BQCV and SBV.

Diverse single-stranded DNA viruses associated with honey bees (Apis mellifera)

Honey bees (Apis mellifera) research has increased in light of their progressive global decline over the last decade and the important role they play in pollination. One expanding area of honey bee research is analysis of their microbial community including viruses. Several RNA viruses have been characterized but little is known about DNA viruses associated with bees. Here, using a metagenomics based approach, we reveal the presence of a broad range of novel single-stranded DNA viruses from the hemolymph and brain of nurse and forager (worker divisions of labour) bees belonging to two honey bees subspecies, Italian (Apis mellifera linguistica) and New World Carniolan (Apis mellifera carnica). Genomes of 100 diverse viruses were identified, designated into three groupings; genomoviruses (family Genomoviridae) (n = 4), unclassified replication associated protein encoding single-stranded DNA viruses (n = 28), and microviruses (family Microviridae; subfamily Gokushovirinae) (n = 70). Amongst the viruses identified, it appears that nurses harbour a higher diversity of these viruses comparative to the foragers. Between subspecies, the most striking outcome was the extremely high number of diverse microviruses identified in the Italian bees comparative to the New World Carniolan, likely indicating an association to the diversity of the bacterial community associated with these subspecies.

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

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

Baculovirus Per Os Infectivity Factor Complex: Components and Assembly

Baculovirus entry into insect midgut cells is dependent on a multiprotein complex of per os infectivity factors (PIFs) on the envelopes of occlusion-derived virions (ODVs). The structure and assembly of the PIF complex are largely unknown. To reveal the complete members of the complex, a combination of blue native polyacrylamide gel electrophoresis, liquid chromatography-tandem mass spectrometry, and Western blotting was conducted on three different baculoviruses. The results showed that the PIF complex has a molecular mass of ∼500 kDa and consists of nine PIFs, including a newly discovered member (PIF9). To decipher the assembly process, each pif gene was knocked out from the Autographa californica multicapsid nucleopolyhedrovirus (AcMNPV) genome individually by use of synthetic baculovirus technology, and the impact on PIF complex formation was investigated. Deletion of pif8 resulted in the formation of an ∼400-kDa subcomplex. Deletion of pif0, -4, -6, -7, or -9 resulted in a subcomplex of ∼230 kDa, but deletion of pif1, -2, or -3 abolished formation of any complex. Taken together, our data identified a core complex of ∼230 kDa, consisting of PIF1, -2, and -3. This revised the previous knowledge that the core complex was about 170 kDa and contained PIF1 to -4. Analysis of the PIF complex in cellular fractions suggested that it is assembled in the cytoplasm before being transported to the nucleus and subsequently incorporated into the envelopes of ODVs. Only the full complex, not the subcomplex, is resistant to proteolytic attack, indicating the essentiality of correct complex assembly for oral infection.IMPORTANCE Entry of baculovirus into host insects is mediated by a per os infectivity factor (PIF) complex on the envelopes of occlusion-derived viruses (ODVs). Knowledge of the composition and structure of the PIF complex is fundamental to understanding its mode of action. By using multiple approaches, we determined the complete list of proteins (nine) in the PIF complex. In contrast to previous knowledge in the field, the core complex is revised to ∼230 kDa and consists of PIF1 to -3 but not PIF4. Interestingly, our results suggest that the PIF complex is formed in the cytoplasm prior to its transport to the nucleus and subsequent incorporation into ODVs. Only the full complex is resistant to proteolytic degradation in the insect midgut, implying the critical role of the entire complex. These findings provide the baseline for future studies on the ODV entry mechanism mediated by the multiprotein complex.

First molecular detection of Apis mellifera filamentous virus in honey bees (Apis mellifera) in Hungary

Western honey bees (Apis mellifera) are important pollinators in the ecosystem and also play a crucial economic role in the honey industry. During the last decades, a continuous decay was registered in honey bee populations worldwide, including Hungary. In our study, we used metagenomic approaches and conventional PCR screening on healthy and winter mortality affected colonies from multiple sites in Hungary. The major goal was to discover presumed bee pathogens with viral metagenomic experiments and gain prevalence and distribution data by targeted PCR screening. We examined 664 honey bee samples that had been collected during winter mortality from three seemingly healthy colonies and from one colony infested heavily by the parasitic mite Varroa destructor in 2016 and 2017. The subsequent PCR screening of honey bee samples revealed the abundant presence of Apis mellifera filamentous virus (AmFV) for the first time in Central Europe. Based on phylogeny reconstruction, the newly-detected virus was found to be most closely related to a Chinese AmFV strain. More sequence data from multiple countries would be needed for studying the detailed phylogeographical patterns and worldwide spreading process of AmFV. Here we report the prevalent presence of this virus in Hungarian honey bee colonies.

New Viruses from the Ectoparasite Mite Varroa destructor Infesting Apis mellifera and Apis cerana

Varroa destructor is an ectoparasitic mite of Asian or Eastern honeybees Apis cerana(A. cerana) which has become a serious threat to European subspecies of Western honeybees Apis mellifera (A. mellifera) within the last century. V.destructor and its vectored honeybee viruses became serious threats for colony survival. This is a short period for pathogen- and host-populations to adapt. To look for possible variation in the composition of viral populations we performed RNA metagenomic analysis of the Western honeybee subspecies A. m. ligustica, A. m.syriaca, A. m. intermissa, and A. cerana and their respective V. destructor mites. The analysis revealed two novel viruses: Varroa orthomyxovirus-1 (VOV-1) in A. mellifera and V. destructor and a Hubei like-virga virus-14 homolog in V. destructor. VOV-1 was more prevalent in V. destructor than in A. mellifera and we found evidence for viral replication in both hosts. Interestingly, we found differences in viral loads of A. cerana and their V. destructor, A. m. intermissa, and its V. destructor showed partial similarity, while A. m.ligustica and A. m.syriaca and their varroa where very similar. Deformed wing virus exhibited 82.20%, 99.20%, 97.90%, and 0.76% of total viral reads in A. m. ligustica, A. m. syriaca, A. m. intermissa, and A. cerana, respectively. This is the first report of a complete segmented-single-stranded negative-sense RNA virus genome in honeybees and V. destructor mites.

Crustacean Genome Exploration Reveals the Evolutionary Origin of White Spot Syndrome Virus

White spot syndrome virus (WSSV) is a crustacean-infecting, double-stranded DNA virus and is the most serious viral pathogen in the global shrimp industry. WSSV is the sole recognized member of the family Nimaviridae, and the lack of genomic data on other nimaviruses has obscured the evolutionary history of WSSV. Here, we investigated the evolutionary history of WSSV by characterizing WSSV relatives hidden in host genomic data. We surveyed 14 host crustacean genomes and identified five novel nimaviral genomes. Comparative genomic analysis of Nimaviridae identified 28 "core genes" that are ubiquitously conserved in Nimaviridae; unexpected conservation of 13 uncharacterized proteins highlighted yet-unknown essential functions underlying the nimavirus replication cycle. The ancestral Nimaviridae gene set contained five baculoviral per os infectivity factor homologs and a sulfhydryl oxidase homolog, suggesting a shared phylogenetic origin of Nimaviridae and insect-associated double-stranded DNA viruses. Moreover, we show that novel gene acquisition and subsequent amplification reinforced the unique accessory gene repertoire of WSSV. Expansion of unique envelope protein and nonstructural virulence-associated genes may have been the key genomic event that made WSSV such a deadly pathogen.IMPORTANCE WSSV is the deadliest viral pathogen threatening global shrimp aquaculture. The evolutionary history of WSSV has remained a mystery, because few WSSV relatives, or nimaviruses, had been reported. Our aim was to trace the history of WSSV using the genomes of novel nimaviruses hidden in host genome data. We demonstrate that WSSV emerged from a diverse family of crustacean-infecting large DNA viruses. By comparing the genomes of WSSV and its relatives, we show that WSSV possesses an expanded set of unique host-virus interaction-related genes. This extensive gene gain may have been the key genomic event that made WSSV such a deadly pathogen. Moreover, conservation of insect-infecting virus protein homologs suggests a common phylogenetic origin of crustacean-infecting Nimaviridae and other insect-infecting DNA viruses. Our work redefines the previously poorly characterized crustacean virus family and reveals the ancient genomic events that preordained the emergence of a devastating shrimp pathogen.

Bee Viruses: Ecology, Pathogenicity, and Impacts

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

A prokaryotic-eukaryotic relation in the fat body of Bombus terrestris

The interaction between the insect host and its microbiota plays a central role in insect health and is mostly studied in relation to the digestive system. Nonetheless, there are numerous microorganisms occupying multiple habitats in and on insects. We studied microbial communities in the gut and fat body of bumblebees (Bombus terrestris) using the V4 region of the 16S rRNA gene on the Illumina MiSeq platform. In one of the two study locations, the fat body microbial composition was marked by the dominant presence of Arsenophonus sp. and Phyllobacterium sp. Bumblebees infected with Apicystis bombi, a eukaryotic parasite multiplying in the fat body, had a significant higher relative abundance of Arsenophonus sp. compared with the non-infected individuals. In general, the infection of A. bombi correlated with a more interlinked microbial association network, as we observed an increase of significant associations between the relative abundance of bacteria present in the gut and fat body of infected bumblebees. The causality within this potential prokaryotic-eukaryotic relation is important when assessing the health impact on bees.

Characterisation of the British honey bee metagenome

The European honey bee (Apis mellifera) plays a major role in pollination and food production. Honey bee health is a complex product of the environment, host genetics and associated microbes (commensal, opportunistic and pathogenic). Improved understanding of these factors will help manage modern challenges to bee health. Here we used DNA sequencing to characterise the genomes and metagenomes of 19 honey bee colonies from across Britain. Low heterozygosity was observed in many Scottish colonies which had high similarity to the native dark bee. Colonies exhibited high diversity in composition and relative abundance of individual microbiome taxa. Most non-bee sequences were derived from known honey bee commensal bacteria or pathogens. However, DNA was also detected from additional fungal, protozoan and metazoan species. To classify cobionts lacking genomic information, we developed a novel network analysis approach for clustering orphan DNA contigs. Our analyses shed light on microbial communities associated with honey bees and demonstrate the power of high-throughput, directed metagenomics for identifying novel biological threats in agroecosystems.

Shotgun metagenomics of honey DNA: Evaluation of a methodological approach to describe a multi-kingdom honey bee derived environmental DNA signature

Honey bees are considered large-scale monitoring tools due to their environmental exploration and foraging activities. Traces of these activities can be recovered in the honey that also may reflect the hive ecological micro-conditions in which it has been produced. This study applied a next generation sequencing platform (Ion Torrent) for shotgun metagenomic analysis of honey environmental DNA (eDNA). The study tested a methodological framework to interpret DNA sequence information useful to describe the complex ecosystems of the honey bee colony superorganism, its pathosphere and the heterogeneity of the agroecological environments and environmental sources that left DNA marks in the honey. Analysis of two honeys reported sequence reads from five main organism groups (kingdoms or phyla): arthropods (that mainly included reads from Apis mellifera, several other members of the Hymenotpera, in addition to members of the Diptera, Coleoptera and Lepidoptera, as well as aphids and mites), plants (that clearly confirmed the botanical origin of the two honeys, i.e. orange tree blossom and eucalyptus tree blossom honeys), fungi and bacteria (including common hive and honey bee gut microorganisms, honey bee pathogens and plant pathogens), and viruses (which accounted for the largest number of reads in both honeys, mainly assigned to Apis mellifera filamentous virus). The shotgun metagenomic approach that was used in this study can be applied in large scale experiments that might have multiple objectives according to the multi-kingdom derived eDNA that is contained in the honey.

Honey Bee and Bumble Bee Antiviral Defense

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

Investigating the viral ecology of global bee communities with high-throughput metagenomics

Bee viral ecology is a fascinating emerging area of research: viruses exert a range of effects on their hosts, exacerbate impacts of other environmental stressors, and, importantly, are readily shared across multiple bee species in a community. However, our understanding of bee viral communities is limited, as it is primarily derived from studies of North American and European Apis mellifera populations. Here, we examined viruses in populations of A. mellifera and 11 other bee species from 9 countries, across 4 continents and Oceania. We developed a novel pipeline to rapidly and inexpensively screen for bee viruses. This pipeline includes purification of encapsulated RNA/DNA viruses, sequence-independent amplification, high throughput sequencing, integrated assembly of contigs, and filtering to identify contigs specifically corresponding to viral sequences. We identified sequences for (+)ssRNA, (−)ssRNA, dsRNA, and ssDNA viruses. Overall, we found 127 contigs corresponding to novel viruses (i.e. previously not observed in bees), with 27 represented by >0.1% of the reads in a given sample, and 7 contained an RdRp or replicase sequence which could be used for robust phylogenetic analysis. This study provides a sequence-independent pipeline for viral metagenomics analysis, and greatly expands our understanding of the diversity of viruses found in bee communities.

Recently identified bee viruses and their impact on bee pollinators

Bees are agriculturally and ecologically important plant pollinators. Recent high annual losses of honey bee colonies, and reduced populations of native and wild bees in some geographic locations, may impact the availability of affordable food crops and the diversity and abundance of native and wild plant species. Multiple factors including viral infections affect pollinator health. The majority of well-characterized bee viruses are picorna-like RNA viruses, which may be maintained as covert infections or cause symptomatic infections or death. Next generation sequencing technologies have been utilized to identify additional bee-infecting viruses including the Lake Sinai viruses and Rhabdoviruses. In addition, sequence data is instrumental for defining specific viral strains and characterizing associated pathogenicity, such as the recent characterization of Deformed wing virus master variants (DWV-A, DWV-B, and DWV-C) and their impact on bee health.

Genome Sequencing of the Behavior Manipulating Virus LbFV Reveals a Possible New Virus Family

Parasites are sometimes able to manipulate the behavior of their hosts. However, the molecular cues underlying this phenomenon are poorly documented. We previously reported that the parasitoid wasp Leptopilina boulardi which develops from Drosophila larvae is often infected by an inherited DNA virus. In addition to being maternally transmitted, the virus benefits from horizontal transmission in superparasitized larvae (Drosophila that have been parasitized several times). Interestingly, the virus forces infected females to lay eggs in already parasitized larvae, thus increasing the chance of being horizontally transmitted. In a first step towards the identification of virus genes responsible for the behavioral manipulation, we present here the genome sequence of the virus, called LbFV. The sequencing revealed that its genome contains an homologous repeat sequence (hrs) found in eight regions in the genome. The presence of this hrs may explain the genomic plasticity that we observed for this genome. The genome of LbFV encodes 108 ORFs, most of them having no homologs in public databases. The virus is however related to Hytrosaviridae, although distantly. LbFV may thus represent a member of a new virus family. Several genes of LbFV were captured from eukaryotes, including two anti-apoptotic genes. More surprisingly, we found that LbFV captured from an ancestral wasp a protein with a Jumonji domain. This gene was afterwards duplicated in the virus genome. We hypothesized that this gene may be involved in manipulating the expression of wasp genes, and possibly in manipulating its behavior.

A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations

Understanding the diversity and consequences of viruses present in honey bees is critical for maintaining pollinator health and managing the spread of disease. The viral landscape of honey bees (Apis mellifera) has changed dramatically since the emergence of the parasitic mite Varroa destructor, which increased the spread of virulent variants of viruses such as deformed wing virus. Previous genomic studies have focused on colonies suffering from infections by Varroa and virulent viruses, which could mask other viral species present in honey bees, resulting in a distorted view of viral diversity. To capture the viral diversity within colonies that are exposed to mites but do not suffer the ultimate consequences of the infestation, we examined populations of honey bees that have evolved naturally or have been selected for resistance to Varroa This analysis revealed seven novel viruses isolated from honey bees sampled globally, including the first identification of negative-sense RNA viruses in honey bees. Notably, two rhabdoviruses were present in three geographically diverse locations and were also present in Varroa mites parasitizing the bees. To characterize the antiviral response, we performed deep sequencing of small RNA populations in honey bees and mites. This provided evidence of a Dicer-mediated immune response in honey bees, while the viral small RNA profile in Varroa mites was novel and distinct from the response observed in bees. Overall, we show that viral diversity in honey bee colonies is greater than previously thought, which encourages additional studies of the bee virome on a global scale and which may ultimately improve disease management.IMPORTANCE Honey bee populations have become increasingly susceptible to colony losses due to pathogenic viruses spread by parasitic Varroa mites. To date, 24 viruses have been described in honey bees, with most belonging to the order Picornavirales Collapsing Varroa-infected colonies are often overwhelmed with high levels of picornaviruses. To examine the underlying viral diversity in honey bees, we employed viral metatranscriptomics analyses on three geographically diverse Varroa-resistant populations from Europe, Africa, and the Pacific. We describe seven novel viruses from a range of diverse viral families, including two viruses that are present in all three locations. In honey bees, small RNA sequences indicate that these viruses are processed by Dicer and the RNA interference pathway, whereas Varroa mites produce strikingly novel small RNA patterns. This work increases the number and diversity of known honey bee viruses and will ultimately contribute to improved disease management in our most important agricultural pollinator.

Per os infectivity factors: a complicated and evolutionarily conserved entry machinery of baculovirus

Baculoviruses are a family of arthropod-specific large DNA viruses that infect insect species belonging to the orders Lepidoptera, Hymenoptera and Diptera. In nature, occlusion-derived viruses (ODVs) initiate baculovirus primary infection in the midgut epithelium of insect hosts, and this process is largely dependent on a number of ODV envelope proteins designated as per os infectivity factors (PIFs). Interestingly, PIF homologs are also present in other invertebrate large DNA viruses, which is indicative that per os infection is an ancient and phylogenetically conserved entry mechanism shared by these viruses. Here, we review the advances in the knowledge of the functions of individual PIFs and recent discoveries about the PIF complex, and discuss the evolutionary implications of PIF homologs in invertebrate DNA viruses. Furthermore, future research highlights on the per os infection mechanism are also prospected.

Advances in the application of high-throughput sequencing in invertebrate virology

Over the last decade, advances in high-throughput sequencing technologies have revolutionised biological research, making it possible for DNA/RNA sequencing of any organism of interest to be undertaken. Sequencing approaches are now routinely used in the detection and characterisation of (novel) viruses, investigation of host-pathogen interactions, and effective development of disease treatment strategies. For the sequencing and identification of viruses of interest, metagenomics approaches using infected host tissue are frequently used, as it is not always possible to culture and isolate these pathogens. High-throughput sequencing can also be used to investigate host-pathogen interactions by investigating (temporal) transcriptomic responses of both the host and virus, potentially leading to the discovery of novel opportunities for treatment and drug targets. In addition, viruses in environmental samples (e.g. water or soil samples) can be identified using eDNA/metagenomics approaches. The promise that recent developments in sequencing brings to the field of invertebrate virology are not devoid of technical challenges, including the need for better laboratory and bioinformatics strategies to sequence and assemble virus genomes within complex tissue or environmental samples, and the difficulties associated with the annotation of the large number of novel viruses being discovered.