Promiscuous feeding on multiple adult honey bee hosts amplifies the vectorial capacity of Varroa destructor.

Varroa destructor is a cosmopolitan pest and leading cause of colony loss of the European honey bee. Historically described as a competent vector of honey bee viruses, this arthropod vector is the cause of a global pandemic of Deformed wing virus, now endemic in honeybee populations in all Varroa-infested regions. Our work shows that viral spread is driven by Varroa actively switching from one adult bee to another as they feed. Assays using fluorescent microspheres were used to indicate the movement of fluids in both directions between host and vector when Varroa feed. Therefore, Varroa could be in either an infectious or naïve state dependent upon the disease status of their host. We tested this and confirmed that the relative risk of a Varroa feeding depended on their previous host's infectiousness. Varroa exhibit remarkable heterogeneity in their host-switching behavior, with some Varroa infrequently switching while others switch at least daily. As a result, relatively few of the most active Varroa parasitize the majority of bees. This multiple-feeding behavior has analogs in vectorial capacity models of other systems, where promiscuous feeding by individual vectors is a leading driver of vectorial capacity. We propose that the honeybee-Varroa relationship offers a unique opportunity to apply principles of vectorial capacity to a social organism, as virus transmission is both vectored and occurs through multiple host-to-host routes common to a crowded society.

Susceptible and infectious states for both vector and host in a dynamic pathogen-vector-host system.

Deformed wing virus (DWV) is a resurgent insect pathogen of honeybees that is efficiently transmitted by vectors and through host social contact. Continual transmission of DWV between hosts and vectors is required to maintain the pathogen within the population, and this vector-host-pathogen system offers unique disease transmission dynamics for pathogen maintenance between vectors and a social host. In a series of experiments, we measured vector-vector, host-host and host-vector transmission routes and show how these maintain DWV in honeybee populations. We found co-infestations on shared hosts allowed for movement of DWV from mite to mite. Additionally, two social behaviours of the honeybee, trophallaxis and cannibalization of pupae, provide routes for horizontal transmission from bee to bee. Circulation of the virus solely among hosts through communicable modes provides a reservoir of DWV for naïve Varroa to acquire and subsequently vector the pathogen. Our findings illustrate the importance of community transmission between hosts and vector transmission. We use these results to highlight the key avenues used by DWV during maintenance and infection and point to similarities with a handful of other infectious diseases of zoonotic and medical importance.

The Temporal and Geographical Dynamics of Potato Virus Y Diversity in Russia.

Potato virus Y, an important viral pathogen of potato, has several genetic variants and geographic distributions which could be affected by environmental factors, aphid vectors, and reservoir plants. PVY is transmitted to virus-free potato plants by aphids and passed on to the next vegetative generations through tubers, but the effects of tuber transmission in PVY is largely unknown. By using high-throughput sequencing, we investigated PVY populations transmitted to potato plants by aphids in different climate zones of Russia, namely the Moscow and Astrakhan regions. We analyzed sprouts from the tubers produced by field-infected plants to investigate the impact of tuber transmission on PVY genetics. We found a significantly higher diversity of PVY isolates in the Astrakhan region, where winters are shorter and milder and summers are warmer compared to the Moscow region. While five PVY types, NTNa, NTNb, N:O, N-Wi, and SYR-I, were present in both regions, SYRI-II, SYRI-III, and 261-4 were only found in the Astrakhan region. All these recombinants were composed of the genome sections derived from PVY types O and N, but no full-length sequences of such types were present. The composition of the PVY variants in the tuber sprouts was not always the same as in their parental plants, suggesting that tuber transmission impacts PVY genetics.

Honeybee intestines retain low yeast titers, but no bacterial mutualists, at emergence.

Honeybee symbionts, predominantly bacteria, play important roles in honeybee health, nutrition, and pathogen protection, thereby supporting colony health. On the other hand, fungi are often considered indicators of poor bee health, and honeybee microbiome studies generally exclude fungi and yeasts. We hypothesized that yeasts may be an important aspect of early honeybee biology, and if yeasts provide a mutual benefit to their hosts, then honeybees could provide a refuge during metamorphosis to ensure the presence of yeasts at emergence. We surveyed for yeast and fungi during pupal development and metamorphosis in worker bees using fungal-specific quantitative polymerase chain reaction (qPCR), next-generation sequencing, and standard microbiological culturing. On the basis of yeast presence in three distinct apiaries and multiple developmental stages, we conclude that yeasts can survive through metamorphosis and in naïve worker bees, albeit at relatively low levels. In comparison, known bacterial mutualists, like Gilliamella and Snodgrassella, were generally not found in pre-eclosed adult bees. Whether yeasts are actively retained as an important part of the bee microbiota or are passively propagating in the colony remains unknown. Our demonstration of the constancy of yeasts throughout development provides a framework to further understand the honeybee microbiota.

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.

Dynamic evolution in the key honey bee pathogen deformed wing virus: Novel insights into virulence and competition using reverse genetics.

The impacts of invertebrate RNA virus population dynamics on virulence and infection outcomes are poorly understood. Deformed wing virus (DWV), the main viral pathogen of honey bees, negatively impacts bee health, which can lead to colony death. Despite previous reports on the reduction of DWV diversity following the arrival of the parasitic mite Varroa destructor, the key DWV vector, we found high genetic diversity of DWV in infested United States honey bee colonies. Phylogenetic analysis showed that divergent US DWV genotypes are of monophyletic origin and were likely generated as a result of diversification after a genetic bottleneck. To investigate the population dynamics of this divergent DWV, we designed a series of novel infectious cDNA clones corresponding to coexisting DWV genotypes, thereby devising a reverse-genetics system for an invertebrate RNA virus quasispecies. Equal replication rates were observed for all clone-derived DWV variants in single infections. Surprisingly, individual clones replicated to the same high levels as their mixtures and even the parental highly diverse natural DWV population, suggesting that complementation between genotypes was not required to replicate to high levels. Mixed clone-derived infections showed a lack of strong competitive exclusion, suggesting that the DWV genotypes were adapted to coexist. Mutational and recombination events were observed across clone progeny, providing new insights into the forces that drive and constrain virus diversification. Accordingly, our results suggest that Varroa influences DWV dynamics by causing an initial selective sweep, which is followed by virus diversification fueled by negative frequency-dependent selection for new genotypes. We suggest that this selection might reflect the ability of rare lineages to evade host defenses, specifically antiviral RNA interference (RNAi). In support of this hypothesis, we show that RNAi induced against one DWV strain is less effective against an alternate strain from the same population.

Invertebrate RNA virus diversity from a taxonomic point of view.

Invertebrates are hosts to diverse RNA viruses that have all possible types of encapsidated genomes (positive, negative and ambisense single stranded RNA genomes, or a double stranded RNA genome). These viruses also differ markedly in virion morphology and genome structure. Invertebrate RNA viruses are present in three out of four currently recognized orders of RNA viruses: Mononegavirales, Nidovirales, and Picornavirales, and 10 out of 37 RNA virus families that have yet to be assigned to an order. This mini-review describes general properties of the taxonomic groups, which include invertebrate RNA viruses on the basis of their current classification by the International Committee on Taxonomy of Viruses (ICTV).

A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission.

The globally distributed ectoparasite Varroa destructor is a vector for viral pathogens of the Western honeybee (Apis mellifera), in particular the Iflavirus Deformed Wing Virus (DWV). In the absence of Varroa low levels DWV occur, generally causing asymptomatic infections. Conversely, Varroa-infested colonies show markedly elevated virus levels, increased overwintering colony losses, with impairment of pupal development and symptomatic workers. To determine whether changes in the virus population were due Varroa amplifying and introducing virulent virus strains and/or suppressing the host immune responses, we exposed Varroa-naïve larvae to oral and Varroa-transmitted DWV. We monitored virus levels and diversity in developing pupae and associated Varroa, the resulting RNAi response and transcriptome changes in the host. Exposed pupae were stratified by Varroa association (presence/absence) and virus levels (low/high) into three groups. Varroa-free pupae all exhibited low levels of a highly diverse DWV population, with those exposed per os (group NV) exhibiting changes in the population composition. Varroa-associated pupae exhibited either low levels of a diverse DWV population (group VL) or high levels of a near-clonal virulent variant of DWV (group VH). These groups and unexposed controls (C) could be also discriminated by principal component analysis of the transcriptome changes observed, which included several genes involved in development and the immune response. All Varroa tested contained a diverse replicating DWV population implying the virulent variant present in group VH, and predominating in RNA-seq analysis of temporally and geographically separate Varroa-infested colonies, was selected upon transmission from Varroa, a conclusion supported by direct injection of pupae in vitro with mixed virus populations. Identification of a virulent variant of DWV, the role of Varroa in its transmission and the resulting host transcriptome changes furthers our understanding of this important viral pathogen of honeybees.

MosaicSolver: a tool for determining recombinants of viral genomes from pileup data.

Viral recombination is a key evolutionary mechanism, aiding escape from host immunity, contributing to changes in tropism and possibly assisting transmission across species barriers. The ability to determine whether recombination has occurred and to locate associated specific recombination junctions is thus of major importance in understanding emerging diseases and pathogenesis. This paper describes a method for determining recombinant mosaics (and their proportions) originating from two parent genomes, using high-throughput sequence data. The method involves setting the problem geometrically and the use of appropriately constrained quadratic programming. Recombinants of the honeybee deformed wing virus and the Varroa destructor virus-1 are inferred to illustrate the method from both siRNAs and reads sampling the viral genome population (cDNA library); our results are confirmed experimentally. Matlab software (MosaicSolver) is available.

A strong immune response in young adult honeybees masks their increased susceptibility to infection compared to older bees.

Honeybees, Apis mellifera, show age-related division of labor in which young adults perform maintenance ("housekeeping") tasks inside the colony before switching to outside foraging at approximately 23 days old. Disease resistance is an important feature of honeybee biology, but little is known about the interaction of pathogens and age-related division of labor. We tested a hypothesis that older forager bees and younger "house" bees differ in susceptibility to infection. We coupled an infection bioassay with a functional analysis of gene expression in individual bees using a whole genome microarray. Forager bees treated with the entomopathogenic fungus Metarhizium anisopliae s.l. survived for significantly longer than house bees. This was concomitant with substantial differences in gene expression including genes associated with immune function. In house bees, infection was associated with differential expression of 35 candidate immune genes contrasted with differential expression of only two candidate immune genes in forager bees. For control bees (i.e. not treated with M. anisopliae) the development from the house to the forager stage was associated with differential expression of 49 candidate immune genes, including up-regulation of the antimicrobial peptide gene abaecin, plus major components of the Toll pathway, serine proteases, and serpins. We infer that reduced pathogen susceptibility in forager bees was associated with age-related activation of specific immune system pathways. Our findings contrast with the view that the immunocompetence in social insects declines with the onset of foraging as a result of a trade-off in the allocation of resources for foraging. The up-regulation of immune-related genes in young adult bees in response to M. anisopliae infection was an indicator of disease susceptibility; this also challenges previous research in social insects, in which an elevated immune status has been used as a marker of increased disease resistance and fitness without considering the effects of age-related development.

Host-driven temperature dependence of Deformed wing virus infection in honey bee pupae.

The temperature dependence of infection reflects changes in performance of parasites and hosts. High temperatures often mitigate infection by favoring heat-tolerant hosts over heat-sensitive parasites. Honey bees exhibit endothermic thermoregulation-rare among insects-that can favor resistance to parasites. However, viruses are heavily host-dependent, suggesting that viral infection could be supported-not threatened-by optimum host function. To understand how temperature-driven changes in performance of viruses and hosts shape infection, we compared the temperature dependence of isolated viral enzyme activity, three honey bee traits, and infection of honey bee pupae. Viral enzyme activity varied <2-fold over a > 30 °C interval spanning temperatures typical of ectothermic insects and honey bees. In contrast, honey bee performance peaked at high (≥ 35 °C) temperatures and was highly temperature-sensitive. Although these results suggested that increasing temperature would favor hosts over viruses, the temperature dependence of pupal infection matched that of pupal development, falling only near pupae's upper thermal limits. Our results reflect the host-dependent nature of viruses, suggesting that infection is accelerated-not curtailed-by optimum host function, contradicting predictions based on relative performance of parasites and hosts, and suggesting tradeoffs between infection resistance and host survival that limit the viability of bee 'fever'.

Deformed wing virus of honey bees is inactivated by cold plasma ionized hydrogen peroxide.

Deformed wing virus (DWV) is a widespread pathogen of Apis mellifera honey bees, and is considered a major causative factor for the collapse of infected honey bee colonies. DWV can be horizontally transmitted among bees through various oral routes, including via food sharing and by interactions of bees with viral-contaminated solid hive substrates. Cold plasma ionized hydrogen peroxide (iHP) is used extensively by the food production, processing and medical industries to clean surfaces of microbial contaminants. In this study, we investigated the use of iHP to inactivate DWV particles in situ on a solid substrate. iHP-treated DWV sources were ~105-fold less infectious when injected into naïve honey bee pupae compared to DWV receiving no iHP treatment, matching injected controls containing no DWV. iHP treatment also greatly reduced the incidence of overt DWV infections (i.e., pupae having >109 copies of DWV). The level of DWV inactivation achieved with iHP treatment was higher than other means of viral inactivation such as gamma irradiation, and iHP treatment is likely simpler and safer. Treatment of DWV contaminated hive substrates with iHP, even with honey bees present, may be an effective way to decrease the impacts of DWV infection on honey bees.

Apis mellifera Solinvivirus-1, a Novel Honey Bee Virus That Remained Undetected for over a Decade, Is Widespread in the USA.

A metagenomic analysis of the virome of honey bees (Apis mellifera) from an apiary with high rates of unexplained colony losses identified a novel RNA virus. The virus, which was named Apis mellifera solinvivirus 1 (AmSV1), contains a 10.6 kb positive-strand genomic RNA with a single ORF coding for a polyprotein with the protease, helicase, and RNA-dependent RNA polymerase domains, as well as a single jelly-roll structural protein domain, showing highest similarity with viruses in the family Solinviviridae. The injection of honey bee pupae with AmSV1 preparation showed an increase in virus titer and the accumulation of the negative-strand of AmSV1 RNA 3 days after injection, indicating the replication of AmSV1. In the infected worker bees, AmSV1 was present in heads, thoraxes, and abdomens, indicating that this virus causes systemic infection. An analysis of the geographic and historic distribution of AmSV1, using over 900 apiary samples collected across the United States, showed AmSV1 presence since at least 2010. In the year 2021, AmSV1 was detected in 10.45% of apiaries (95%CI: 8.41-12.79%), mostly sampled in June and July in Northwestern and Northeastern United States. The diagnostic methods and information on the AmSV1 distribution will be used to investigate the connection of AmSV1 to honey bee colony losses.

Densovirus induces winged morphs in asexual clones of the rosy apple aphid, Dysaphis plantaginea.

Winged morphs of aphids are essential for their dispersal and survival. We discovered that the production of the winged morph in asexual clones of the rosy apple aphid, Dysaphis plantaginea, is dependent on their infection with a DNA virus, Dysaphis plantaginea densovirus (DplDNV). Virus-free clones of the rosy apple aphid, or clones infected singly with an RNA virus, rosy apple aphid virus (RAAV), did not produce the winged morph in response to crowding and poor plant quality. DplDNV infection results in a significant reduction in aphid reproduction rate, but such aphids can produce the winged morph, even at low insect density, which can fly and colonize neighboring plants. Aphids infected with DplDNV produce a proportion of virus-free aphids, which enables production of virus-free clonal lines after colonization of a new plant. Our data suggest that a mutualistic relationship exists between the rosy apple aphid and its viruses. Despite the negative impact of DplDNV on rosy apple aphid reproduction, this virus contributes to their survival by inducing wing development and promoting dispersal.

Beeporter: Tools for high-throughput analyses of pollinator-virus infections.

Pollinators are in decline thanks to the combined stresses of disease, pesticides, habitat loss, and climate. Honey bees face numerous pests and pathogens but arguably none are as devastating as Deformed wing virus (DWV). Understanding host-pathogen interactions and virulence of DWV in honey bees is slowed by the lack of cost-effective high-throughput screening methods for viral infection. Currently, analysis of virus infection in bees and their colonies is tedious, requiring a well-equipped molecular biology laboratory and the use of hazardous chemicals. Here we describe virus clones tagged with green fluorescent protein (GFP) or nanoluciferase (nLuc) that provide high-throughput detection and quantification of virus infections. GFP fluorescence is measured noninvasively in living bees via commonly available long-wave UV light sources and a smartphone camera, or a standard ultraviolet transilluminator gel imaging system. Nonlethal monitoring with GFP allows continuous screening of virus growth and serves as a direct breeding tool for identifying honey bee parents with increased antiviral resistance. Expression using the nLuc reporter strongly correlates with virus infection levels and is especially sensitive. Using multiple reporters, it is also possible to visualize competition, differential virulence, and host tissue targeting by co-occuring pathogens. Finally, it is possible to directly assess the risk of cross-species "spillover" from honey bees to other pollinators and vice versa.

The vectoring competence of the mite Varroa destructor for deformed wing virus of honey bees is dynamic and affects survival of the mite.

The ectoparasitic mite, Varroa destructor and the viruses it vectors, including types A and B of Deformed wing virus (DWV), pose a major threat to honey bees, Apis mellifera. Analysis of 256 mites collected from the same set of field colonies on five occasions from May to October 2021 showed that less than a half of them, 39.8% (95% confidence interval (CI): 34.0 - 46.0%), were able to induce a high (overt) level DWV infection with more than 109 viral genomes per bee in the pupa after 6 days of feeding, with both DWV-A and DWV-B being vectored at similar rates. To investigate the effect of the phoretic (or dispersal) stage on adult bees on the mites' ability to vector DWV, the mites from two collection events were divided into two groups, one of which was tested immediately for their infectiveness, and the other was kept with adult worker bees in cages for 12 days prior to testing their infectiveness. We found that while 39.2% (95% CI: 30.0 - 49.1%) of the immediately tested mites induced overt-level infections, 12-day passage on adult bees significantly increased the infectiousness to 89.8% (95% CI: 79.2 - 95.6%). It is likely that Varroa mites that survive brood interruptions in field colonies are increasingly infectious. The mite lifespan was affected by the DWV type it transmitted to pupae. The mites, which induced high DWV-B but not DWV-A infection had an average lifespan of 15.5 days (95% CI: 11.8 - 19.2 days), which was significantly shorter than those of the mites which induced high DWV-A but not DWV-B infection, with an average lifespan of 24.3 days (95% CI: 20.2 - 28.5), or the mites which did not induce high levels of DWV-A or DWV-B, with an average survival of 21.2 days (95% CI: 19.0 - 23.5 days). The mites which transmitted high levels of both DWV-A and DWV-B had an intermediate average survival of 20.5 days (95% CI: 15.1 - 25.9 days). The negative impact of DWV-B on mite survival could be a consequence of the ability of DWV-B, but not DWV-A to replicate in Varroa.

Recombinants between Deformed wing virus and Varroa destructor virus-1 may prevail in Varroa destructor-infested honeybee colonies.

We have used high-throughput Illumina sequencing to identify novel recombinants between Deformed wing virus (DWV) and Varroa destructor virus-1 (VDV-1), which accumulate to higher levels than DWV in both honeybees and Varroa destructor mites. The recombinants, VDV-1(VVD) and VDV-1(DVD), exhibit crossovers between the 5'-UTR and the regions encoding the structural (capsid) and non-structural viral proteins. This implies that the genomes are modular and that each region may evolve independently, as demonstrated in human enteroviruses. Individual honeybee pupae were infected with a mixture of observed recombinants and DWV. A strong correlation was observed between VDV-1(DVD) levels in honeybee pupae and associated mites, suggesting that this recombinant, with a DWV-derived 5'-UTR and non-structural protein region flanking a VDV-1-derived capsid-encoding region, is better adapted to transmission between V. destructor and honeybees than the parental DWV or a recombinant bearing the VDV-1-derived 5'-UTR (VDV-1(VVD)).

Special Issue "Evolution and Diversity of Insect Viruses".

Insects are crucial for ecosystem functions and services and directly influence human well-being and health [...].

Pupal cannibalism by worker honey bees contributes to the spread of deformed wing virus.

Transmission routes impact pathogen virulence and genetics, therefore comprehensive knowledge of these routes and their contribution to pathogen circulation is essential for understanding host-pathogen interactions and designing control strategies. Deformed wing virus (DWV), a principal viral pathogen of honey bees associated with increased honey bee mortality and colony losses, became highly virulent with the spread of its vector, the ectoparasitic mite Varroa destructor. Reproduction of Varroa mites occurs in capped brood cells and mite-infested pupae from these cells usually have high levels of DWV. The removal of mite-infested pupae by worker bees, Varroa Sensitive Hygiene (VSH), leads to cannibalization of pupae with high DWV loads, thereby offering an alternative route for virus transmission. We used genetically tagged DWV to investigate virus transmission to and between worker bees following pupal cannibalisation under experimental conditions. We demonstrated that cannibalization of DWV-infected pupae resulted in high levels of this virus in worker bees and that the acquired virus was then transmitted between bees via trophallaxis, allowing circulation of Varroa-vectored DWV variants without the mites. Despite the known benefits of hygienic behaviour, it is possible that higher levels of VSH activity may result in increased transmission of DWV via cannibalism and trophallaxis.

Interaction of a plant virus-encoded protein with the major nucleolar protein fibrillarin is required for systemic virus infection.

The nucleolus and specific nucleolar proteins are involved in the life cycles of some plant and animal viruses, but the functions of these proteins and of nucleolar trafficking in virus infections are largely unknown. The ORF3 protein of the plant virus, groundnut rosette virus (an umbravirus), has been shown to cycle through the nucleus, passing through Cajal bodies to the nucleolus and then exiting back into the cytoplasm. This journey is absolutely required for the formation of viral ribonucleoprotein particles (RNPs) that, themselves, are essential for the spread of the virus to noninoculated leaves of the shoot tip. Here, we show that these processes rely on the interaction of the ORF3 protein with fibrillarin, a major nucleolar protein. Silencing of the fibrillarin gene prevents long-distance movement of groundnut rosette virus but does not affect viral replication or cell-to-cell movement. Repressing fibrillarin production also localizes the ORF3 protein to multiple Cajal body-like aggregates that fail to fuse with the nucleolus. Umbraviral ORF3 protein and fibrillarin interact in vitro and, when mixed with umbravirus RNA, form an RNP complex. This complex has a filamentous structure with some regular helical features, resembling the RNP complex formed in vivo during umbravirus infection. The filaments formed in vitro are infectious when inoculated to plants, and their infectivity is resistant to RNase. These results demonstrate previously undescribed functions for fibrillarin as an essential component of translocatable viral RNPs and may have implications for other plant and animal viruses that interact with the nucleolus.