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.

ICTV Virus Taxonomy Profile: Solinviviridae.

Solinviviridae is a family of picorna/calici-like viruses with non-segmented, linear, positive-sense RNA genomes of approximately 10-11 kb. Unusually, their capsid proteins are encoded towards the 3'-end of the genome where they can be expressed both from a subgenomic RNA and as an extension of the replication (picorna-like helicase-protease-polymerase) polyprotein. Members of two species within the family infect ants, but related unclassified virus sequences derive from a large variety of insects and other arthropods. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the Solinviviridae, which is available at www.ictv.global/report/solinviviridae.

ICTV Virus Taxonomy Profile: Polycipiviridae.

Polycipiviridae is a family of picorna-like viruses with non-segmented, linear, positive-sense RNA genomes of approximately 10-12 kb. Unusually for viruses within the order Picornavirales, their genomes are polycistronic, with four (or more) consecutive 5'-proximal open reading frames (ORFs) encoding structural (and possibly other) proteins and a long 3' ORF encoding the replication polyprotein. Members of species within the family have all been detected in ants or via arthropod transcriptomic datasets. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the Polycipiviridae, which is available at www.ictv.global/report/polycipiviridae.

A Genetic Network for Systemic RNA Silencing in Plants.

Non-cell autonomous RNA silencing can spread from cell to cell and over long distances in animals and plants. However, the genetic requirements and signals involved in plant mobile gene silencing are poorly understood. Here, we identified a DICER-LIKE2 (DCL2)-dependent mechanism for systemic spread of posttranscriptional RNA silencing, also known as posttranscriptional gene silencing (PTGS), in Nicotiana benthamiana Using a suite of transgenic DCL RNAi lines coupled with a GFP reporter, we demonstrated that N. benthamiana DCL1, DCL2, DCL3, and DCL4 are required to produce microRNAs and 22, 24, and 21nt small interfering RNAs (siRNAs), respectively. All investigated siRNAs produced in local incipient cells were present at low levels in distal tissues. Inhibition of DCL2 expression reduced the spread of gene silencing, while suppression of DCL3 or DCL4 expression enhanced systemic PTGS. In contrast to DCL4 RNAi lines, DCL2-DCL4 double-RNAi lines developed systemic PTGS similar to that observed in DCL2 RNAi. We further showed that the 21 or 24 nt local siRNAs produced by DCL4 or DCL3 were not involved in long-distance gene silencing. Grafting experiments demonstrated that DCL2 was required in the scion to respond to the signal, but not in the rootstock to produce/send the signal. These results suggest a coordinated DCL genetic pathway in which DCL2 plays an essential role in systemic PTGS in N. benthamiana, while both DCL4 and DCL3 attenuate systemic PTGS. We discuss the potential role of 21, 22, and 24 nt siRNAs in systemic PTGS.

Roles of Dicer-Like Proteins 2 and 4 in Intra- and Intercellular Antiviral Silencing.

RNA silencing is an innate antiviral mechanism conserved in organisms across kingdoms. Such a cellular defense involves DICER or DICER-LIKEs (DCLs) that process plant virus RNAs into viral small interfering RNAs (vsiRNAs). Plants encode four DCLs that play diverse roles in cell-autonomous intracellular virus-induced RNA silencing (known as VIGS) against viral invasion. VIGS can spread between cells. However, the genetic basis and involvement of vsiRNAs in non-cell-autonomous intercellular VIGS remains poorly understood. Using GFP as a reporter gene together with a suite of DCL RNAi transgenic lines, here we show that despite the well-established activities of DCLs in intracellular VIGS and vsiRNA biogenesis, DCL4 acts to inhibit intercellular VIGS whereas DCL2 is required (likely along with DCL2-processed/dependent vsiRNAs and their precursor RNAs) for efficient intercellular VIGS trafficking from epidermal to adjacent cells. DCL4 imposed an epistatic effect on DCL2 to impede cell-to-cell spread of VIGS. Our results reveal previously unknown functions for DCL2 and DCL4 that may form a dual defensive frontline for intra- and intercellular silencing to double-protect cells from virus infection in Nicotiana benthamiana.

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.