RT-PCR analysis of Deformed wing virus in honeybees (Apis mellifera) and mites (Varroa destructor)

A model of infection in honeybee colonies with social immunity

Honeybees (Apis mellifera) play a significant role in the pollination of various food crops and plants. In the past decades, honeybee management has been challenged with increased pathogen and environmental pressure associating with increased beekeeping costs, having a marked economic impact on the beekeeping industry. Pathogens have been identified as a contributing cause of colony losses. Evidence suggested a possible route of pathogen transmission among bees via oral-oral contacts through trophallaxis. Here we propose a model that describes the transmission of an infection within a colony when bee members engage in the trophallactic activity to distribute nectar. In addition, we examine two important features of social immunity, defined as collective disease defenses organized by honeybee society. First, our model considers the social segregation of worker bees. The segregation limits foragers, which are highly exposed to pathogens during foraging outside the nest, from interacting with bees residing in the inner parts of the nest. Second, our model includes a hygienic response, by which healthy nurse bees exterminate infected bees to mitigate horizontal transmission of the infection to other bee members. We propose that the social segregation forms the first line of defense in reducing the uptake of pathogens into the colony. If the first line of defense fails, the hygienic behavior provides a second mechanism in preventing disease spread. Our study identifies the rate of egg-laying as a critical factor in maintaining the colony's health against an infection. We propose that winter conditions which cease or reduce the egg-laying activity combined with an infection in early spring can compromise the social immunity defenses and potentially cause colony losses.

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.

Transcriptomic Responses of the Honey Bee Brain to Infection with Deformed Wing Virus

Managed colonies of European honey bees (Apis mellifera) are under threat from Varroa destructor mite infestation and infection with viruses vectored by mites. In particular, deformed wing virus (DWV) is a common viral pathogen infecting honey bees worldwide that has been shown to induce behavioral changes including precocious foraging and reduced associative learning. We investigated how DWV infection of bees affects the transcriptomic response of the brain. The transcriptomes of individual brains were analyzed using RNA-Seq after experimental infection of newly emerged adult bees with DWV. Two analytical methods were used to identify differentially expressed genes from the ~15,000 genes in the Apis mellifera genome. The 269 genes that had increased expression in DWV infected brains included genes involved in innate immunity such as antimicrobial peptides (AMPs), Ago2, and Dicer. Single bee brain NMR metabolomics methodology was developed for this work and indicates that proline is strongly elevated in DWV infected brains, consistent with the increased presence of the AMPs abaecin and apidaecin. The 1361 genes with reduced expression levels includes genes involved in cellular communication including G-protein coupled, tyrosine kinase, and ion-channel regulated signaling pathways. The number and function of the downregulated genes suggest that DWV has a major impact on neuron signaling that could explain DWV related behavioral changes.

Direct Evidence for Infection of Varroa destructor Mites with the Bee-Pathogenic Deformed Wing Virus Variant B - but Not Variant A - via Fluorescence-in situ-Hybridization Analysis

Deformed wing virus (DWV) is a bee pathogenic, single- and positive-stranded RNA virus that has been involved in severe honey bee colony losses worldwide. DWV, when transmitted horizontally or vertically from bee to bee, causes mainly covert infections not associated with any visible symptoms or damage. Overt infections occur after vectorial transmission of DWV to the developing bee pupae through the ectoparasitic mite Varroa destructor Symptoms of overt infections are pupal death, bees emerging with deformed wings and shortened abdomens, or cognitive impairment due to brain infection. So far, three variants of DWV, DWV-A, DWV-B, and DWV-C, have been described. While it is widely accepted that V. destructor acts as vector of DWV, the question of whether the mite only functions as a mechanical vector or whether DWV can infect the mite thus using it as a biological vector is hotly debated, because in the literature data can be found that support both hypotheses. In order to settle this scientific dispute, we analyzed putatively DWV-infected mites with a newly established protocol for fluorescence-in situ-hybridization of mites and demonstrated DWV-specific signals inside mite cells. We provide compelling and direct evidence that DWV-B infects the intestinal epithelium and the salivary glands of V. destructor In contrast, no evidence for DWV-A infecting mite cells was found. Our data are key to understanding the pathobiology of DWV, the mite's role as a biological DWV vector and the quasispecies dynamics of this RNA virus when switching between insect and arachnid host species.IMPORTANCE Deformed wing virus (DWV) is a bee pathogenic, originally rather benign, single- and positive-stranded RNA virus. Only the vectorial transmission of this virus to honey bees by the ectoparasitic mite Varroa destructor leads to fatal or symptomatic infections of individuals, usually followed by collapse of the entire colony. Studies on whether the mite only acts as a mechanical virus vector or whether DWV can infect the mite and thus use it as a biological vector have led to disparate results. In our study using fluorescence-in situ-hybridization we provide compelling and direct evidence that at least the DWV-B variant infects the gut epithelium and the salivary glands of V. destructor Hence, the host range of DWV includes both, bees (Insecta) and mites (Arachnida). Our data contribute to a better understanding of the triangular relationship between honey bees, V. destructor and DWV and the evolution of virulence in this viral bee pathogen.

Deformed wing virus type a and b in managed honeybee colonies of Argentina

Apis mellifera is infected by more than 24 virus species worldwide, mainly positive-sense, single-stranded RNA viruses of the Dicistroviridae and Iflaviridae families. Among the viruses that infect honey bees, Deformed wing virus is the most prevalent and is present as three master variants DWV-A, B, and C. Given that the ectoparasitic mite Varroa destructor vectors these virus variants, recombination events between them are expected, and variants and their recombinants can co-exist in mites and honeybees at the same time. In this study, we detect, through RT-qPCR, the presence of DWV-A and B in the same samples of adult bees from colonies of Argentina. Total RNA was extracted from pools of ten adult bees from 45 apiaries distributed across the main beekeeping Provinces of Argentina (Buenos Aires, Santa Fe, Córdoba, Santiago del Estero, Río Negro, and Mendoza); then RT-qPCR reactions were performed to detect DWV-A and B, with specific primer pairs. After the amplifications, PCR products (204 and 660 bp amplicons for DWV-B, and ~250 bp for DWV-A) were purified and sequenced to verify that they corresponded to reported sequences, analyzing them using the Blast software. Of the 45 samples analyzed by RT-qPCR, over 90% were infected with DWV-A and 47% were also positive for DWV-B, where it was found in high prevalence specifically in colonies of A. mellifera of the Buenos Aires Province. Future studies will determine the impact of this type of the virus and its ability to recombine with the other DWV types in the apiaries of our country.

A Novel Iflavirus Was Discovered in Green Rice Leafhopper Nephotettix cincticeps and Its Proliferation Was Inhibited by Infection of Rice Dwarf Virus

The green rice leafhopper, Nephotettix cincticeps (Hemiptera: Cicadellidae), is a key insect vector transmitting rice dwarf virus (RDV) that causes rice dwarf disease. We discovered a novel iflavirus from the transcriptomes of N. cincticeps and named it as Nephotettix cincticeps positive-stranded RNA virus-1 (NcPSRV-1). The viral genome consists of 10,524 nucleotides excluding the poly(A) tail and contains one predicted open reading frame encoding a polyprotein of 3,192 amino acids, flanked by 5' and 3' untranslated regions. NcPSRV-1 has a typical iflavirus genome arrangement and is clustered with the members of the family Iflaviridae in the phylogenetic analysis. NcPSRV-1 was detected in all tested tissues and life stages of N. cincticeps and could be transmitted horizontally and vertically. Moreover, NcPSRV-1 had high prevalence in the laboratory populations and was widely spread in field populations of N. cincticeps. NcPSRV-1 could also infect the two-striped leafhopper, Nephotettix apicalis, at a 3.33% infection rate, but was absent in the zigzag leafhopper, Recilia dorsalis, and rice Oryza sativa variety TN1. The infection of RDV altered the viral load and infection rate of NcPSRV-1 in N. cincticeps, for which it seems that RDV has an antagonistic effect on NcPSRV-1 infection in the host.

Deformed wing virus: using reverse genetics to tackle unanswered questions about the most important viral pathogen of honey bees

Deformed wing virus (DWV) is the most important viral pathogen of honey bees. It usually causes asymptomatic infections but, when vectored by the ectoparasitic mite Varroa destructor, it is responsible for the majority of overwintering colony losses globally. Although DWV was discovered four decades ago, research has been hampered by the absence of an in vitro cell culture system or the ability to culture pure stocks of the virus. The recent developments of reverse genetic systems for DWV go some way to addressing these limitations. They will allow the investigation of specific questions about strain variation, host tropism and pathogenesis to be answered, and are already being exploited to study tissue tropism and replication in Varroa and non-Apis pollinators. Three areas neatly illustrate the advances possible with reverse genetic approaches: (i) strain variation and recombination, in which reverse genetics has highlighted similarities rather than differences between virus strains; (ii) analysis of replication kinetics in both honey bees and Varroa, in studies that likely explain the near clonality of virus populations often reported; and (iii) pathogen spillover to non-Apis pollinators, using genetically tagged viruses to accurately monitor replication and infection.

Review of Mites Found on Various Animal Hosts and at Different Localities in Malaysia

Mite biodiversity and distribution in Malaysia is currently understudied. Most previous works on Malaysian Acari have focused on pest organisms of medical, veterinary, and agricultural concern, with a few recent studies centered on mites in forensic contexts. Previous literatures have targeted collection sites in forest reserves and/or mountains in either Peninsular or Malaysian Borneo, though the state of Sarawak had the least publications related to mite species descriptions despite having the highest number of nature parks of any state in the country. Most publications focused on the three states Selangor, Pahang and Sabah. Most of the mite species reported were from mammals (66.3%), with fewer species from birds (21.7%), arthropods (11.2%), and reptiles (0.8%). We believe that further work on the systematic documentation of mite species throughout Malaysia is necessary as it could generate useful tools, such as the use of mites as biogeographical markers or as forensic indicators. Therefore, this review catalogs mite species that have been documented in or on animal hosts in Malaysia and serves as a foundation for future work.

The Combined Effects of Varroa destructor Parasitism and Exposure to Neonicotinoids Affects Honey Bee (Apis mellifera L.) Memory and Gene Expression

Honey bees (Apis mellifera L.) are exposed biotic and abiotic stressors but little is known about their combined effect and impact on neural processes such as learning and memory, which could affect behaviours that are important for individual and colony survival. This study measured memory with the proboscis extension response (PER) assay as well as the expression of neural genes in bees chronically exposed to three different sublethal doses of the insecticide clothianidin and/or the parasitic mite Varroa destructor. The proportion of bees that positively responded to PER at 24 and 48 h post-training (hpt) was significantly reduced when exposed to clothianidin. V. destructor parasitism reduced the proportion of bees that responded to PER at 48 hpt. Combined effects between the lowest clothianidin dose and V. destructor for the proportion of bees that responded to PER were found at 24 hpt. Clothianidin, V. destructor and their combination differentially affected the expression of the neural-related genes, AmNrx-1 (neurexin), AmNlg-1 (neuroligin), and AmAChE-2 (acetylcholinesterase). Different doses of clothianidin down-regulated or up-regulated the genes, whereas V. destructor tended to have a down-regulatory effect. It appears that clothianidin and V. destructor affected neural processes in honey bees through different mechanisms.

Honey Bee Virus Transmission via Hive Products

The global trade of honey bee hive products has raised concern about pathogen transmission. However, the efficacy of hive products as virus vehicles is poorly understood. Here, we investigated the transmission capacity of hive products for Deformed wing virus genotype A (DWV-A) in a fully-crossed hoarding cage experiment and estimated the transmission risk by screening commercial products. Western honey bee workers were provided with honey, pollen and wax either contaminated with high (~2 × 10⁹), medium (~1.7 × 10⁸), low (~8 × 10⁶) or zero (control) DWV-A genome copies. For 10 days, mortality was monitored. Then, virus titers were quantified in bee heads and 38 commercial products using RT-qPCR. For honey and pollen, a positive association between DWV-A concentration and mortality was observed. High concentrations always resulted in infections, medium ones in 47% of cases and low ones in 20% of cases. No significant difference was observed between the tested products. In commercial honey and pollen, 7.7 × 10²-1.8 × 10⁵ and 1.4 × 10³-1.3 × 10⁴ DWV-A copies per gram were found, respectively. The results show that DWV-A transmission via hive products is feasible. The risk of introducing novel viruses and/or strains should be considered in trade regulations by including virus analyses for health certificates of hive products.

An updated genetic marker for detection of Lake Sinai Virus and metagenetic applications

Background

Lake Sinai Viruses (LSV) are common RNA viruses of honey bees (Apis mellifera) that frequently reach high abundance but are not linked to overt disease. LSVs are genetically heterogeneous and collectively widespread, but despite frequent detection in surveys, the ecological and geographic factors structuring their distribution in A. mellifera are not understood. Even less is known about their distribution in other species. Better understanding of LSV prevalence and ecology have been hampered by high sequence diversity within the LSV clade.

Methods

Here we report a new polymerase chain reaction (PCR) assay that is compatible with currently known lineages with minimal primer degeneracy, producing an expected 365 bp amplicon suitable for end-point PCR and metagenetic sequencing. Using the Illumina MiSeq platform, we performed pilot metagenetic assessments of three sample sets, each representing a distinct variable that might structure LSV diversity (geography, tissue, and species).

Results

The first sample set in our pilot assessment compared cDNA pools from managed A. mellifera hives in California (n = 8) and Maryland (n = 6) that had previously been evaluated for LSV2, confirming that the primers co-amplify divergent lineages in real-world samples. The second sample set included cDNA pools derived from different tissues (thorax vs. abdomen, n = 24 paired samples), collected from managed A. mellifera hives in North Dakota. End-point detection of LSV frequently differed between the two tissue types; LSV metagenetic composition was similar in one pair of sequenced samples but divergent in a second pair. Overall, LSV1 and intermediate lineages were common in these samples whereas variants clustering with LSV2 were rare. The third sample set included cDNA from individual pollinator specimens collected from diverse landscapes in the vicinity of Lincoln, Nebraska. We detected LSV in the bee Halictus ligatus (four of 63 specimens tested, 6.3%) at a similar rate as A. mellifera (nine of 115 specimens, 7.8%), but only one H. ligatus sequencing library yielded sufficient data for compositional analysis. Sequenced samples often contained multiple divergent LSV lineages, including individual specimens. While these studies were exploratory rather than statistically powerful tests of hypotheses, they illustrate the utility of high-throughput sequencing for understanding LSV transmission within and among species.

The Pathogen Profile of a Honey Bee Queen Does Not Reflect That of Her Workers

Throughout a honey bee queen's lifetime, she is tended to by her worker daughters, who feed and groom her. Such interactions provide possible horizontal transmission routes for pathogens from the workers to the queen, and as such a queen's pathogen profile may be representative of the workers within a colony. To explore this further, we investigated known honey bee pathogen co-occurrence, as well as pathogen transmission from workers to queens. Queens from 42 colonies were removed from their source hives and exchanged into a second, unrelated foster colony. Worker samples were taken from the source colony on the day of queen exchange and the queens were collected 24 days after introduction. All samples were screened for Nosema spp., Trypanosome spp., acute bee paralysis virus (ABPV), black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), Israeli acute paralysis virus (IAPV), Lake Sinai virus (LSV), and deformed wing virus master variants (DWV-A, B, and C) using RT-qPCR. The data show that LSV, Nosema, and DWV-B were the most abundant pathogens in colonies. All workers (n = 42) were LSV-positive, 88% were Nosema-positive, whilst pathogen loads were low (<1 × 106 genome equivalents per pooled worker sample). All queens (n = 39) were negative for both LSV and Nosema. We found no evidence of DWV transmission occurring from worker to queen when comparing queens to foster colonies, despite DWV being present in both queens and workers. Honey bee pathogen presence and diversity in queens cannot be revealed from screening workers, nor were pathogens successfully transmitted to the queen.

Interaction of Varroa destructor and Sublethal Clothianidin Doses during the Larval Stage on Subsequent Adult Honey Bee (Apis mellifera L.) Health, Cellular Immunity, Deformed Wing Virus Levels and Differential Gene Expression

Honeybees (Apis mellifera L.) are exposed to many parasites, but little is known about interactions with abiotic stressors on their health, particularly when affected as larvae. Larvae were exposed singly and in combination to the parasitic mite Varroa destructor and three sublethal doses of the neonicotinoid insecticide clothianidin to evaluate their effects on survivorship, weight, haemocyte counts, deformed wing virus (DWV) levels and gene expression of the adult bees that subsequently developed. Clothianidin significantly reduced bee weight at the highest dose and was associated with an increase in haemocyte counts at the lowest dose, whereas V. destructor parasitism increased DWV levels, reduced bee emergence, lowered weight and reduced haemocyte counts. An interaction between the two stressors was observed for weight at emergence. Among the differentially expressed genes (DEGs), V. destructor infestation resulted in broader down-regulatory effects related to immunity that was often shared with the combined stressors, while clothianidin resulted in a broader up-regulatory effect more related to central metabolic pathways that was often shared with the combined stressors. Parasites and abiotic stressors can have complex interactions, including additive effects on reduced weight, number of up-regulated DEGs and biological pathways associated with metabolism.

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.

Honey Bee Suppresses the Parasitic Mite Vitellogenin by Antimicrobial Peptide

The negative effects of honey bee parasitic mites and deformed wing virus (DWV) on honey bee and colony health have been well characterized. However, the relationship between DWV and mites, particularly viral replication inside the mites, remains unclear. Furthermore, the physiological outcomes of honey bee immune responses stimulated by DWV and the mite to the host (honey bee) and perhaps the pathogen/parasite (DWV/mite) are not yet understood. To answer these questions, we studied the tripartite interactions between the honey bee, Tropilaelaps mercedesae, and DWV as the model. T. mercedesae functioned as a vector for DWV without supporting active viral replication. Thus, DWV negligibly affected mite fitness. Mite infestation induced mRNA expression of antimicrobial peptides (AMPs), Defensin-1 and Hymenoptaecin, which correlated with DWV copy number in honey bee pupae and mite feeding, respectively. Feeding T. mercedesae with fruit fly S2 cells heterologously expressing honey bee Hymenoptaecin significantly downregulated mite Vitellogenin expression, indicating that the honey bee AMP manipulates mite reproduction upon feeding on bee. Our results provide insights into the mechanism of DWV transmission by the honey bee parasitic mite to the host, and the novel role of AMP in defending against mite infestation.

Effect of Oral Administration of 1,3-1,6 β-Glucans in DWV Naturally Infected Newly Emerged Bees (Apis mellifera L.)

Honeybee pathogens have an important role in honeybee colony mortality and colony losses; most of them are widely spread and necessitate worldwide solutions to contrast honeybee's decline. Possible accepted solutions to cope with the spread of honeybee's pathogens are focused on the study of experimental protocols to enhance the insect's immune defenses. Honeybee's artificial diet capable to stimulate the immune system is a promising field of investigation as ascertained by the introduction of 1,3-1,6 β-glucans as a dietary supplement. In this work, by collecting faecal samples of honeybees exposed to different dietary conditions of 1,3-1,6 β-glucans (0.5% and 2% w/w), it has been possible to investigate the Deformed wing virus (DWV) viral load kinetic without harming the insects. Virological data obtained by a one-step TaqMan RT-PCR highlighted the ability of 1,3-1,6 β-glucans to reduce the viral load at the 24th day of rearing. The results indicated that the diet supplemented with 1,3-1,6 β-glucans was associated with a dose-dependent activation of phenoloxidase. The control group showed a higher survival rate than the experimental groups. This research confirmed 1,3-1,6 β-glucans as molecules able to modulate honeybees' defense pathways, and this is the first report in which the kinetic of DWV infection in honeybee faeces has been monitored by a RT-qPCR.

The honeybee (Apis mellifera) developmental state shapes the genetic composition of the deformed wing virus-A quasispecies during serial transmission

The main biological threat to the western honeybee (Apis mellifera) is the parasitic mite Varroa destructor, largely because it vectors lethal epidemics of honeybee viruses that, in the absence of this mite, are relatively innocuous. The severe pathology is a direct consequence of excessive virus titres caused by this novel transmission route. However, little is known about how the virus adapts genetically during transmission and whether this influences the pathology. Here, we show that upon injection into honeybee pupae, the deformed wing virus type-A (DWV-A) quasispecies undergoes a rapid, extensive expansion of its sequence space, followed by strong negative selection towards a uniform, common shape by the time the pupae have completed their development, with no difference between symptomatic and asymptomatic adults in either DWV titre or genetic composition. This suggests that the physiological and molecular environment during pupal development has a strong, conservative influence on shaping the DWV-A quasispecies in emerging adults. There was furthermore no evidence of any progressive adaptation of the DWV-A quasispecies to serial intra-abdominal injection, simulating mite transmission, despite the generation of ample variation immediately following each transmission, suggesting that the virus either had already adapted to transmission by injection, or was unaffected by it.

Biological characterization of Euscelidius variegatus iflavirus 1

Virus-based biocontrol technologies represent sustainable alternatives to pesticides and insecticides. Phytoplasmas are prokaryotic plant pathogens causing severe losses to crops worldwide. Novel approaches are needed since insecticides against their insect vectors and rogueing of infected plants are the only available strategies to counteract phytoplasma diseases. A new iflavirus, named EVV-1, has been described in the leafhopper phytoplasma vector Euscelidius variegatus, raising the potential to use virus-based application strategies against phytoplasma disease. Here transmission routes of EVV-1 are characterized, and localization within the host reveals the mechanism of insect tolerance to virus infection. Both vertical and horizontal transmission of EVV-1 occur and vertical transmission was more efficient. The virus is systemic and occurs in all life-stages, with the highest loads measured in ovaries and first to third instar nymphs. The basic knowledge gained here on the biology of the virus is crucial for possible future application of iflaviruses as biocontrol agents.

Development of a Honey Bee RNA Virus Vector Based on the Genome of a Deformed Wing Virus

We developed a honey bee RNA-virus vector based on the genome of a picorna-like Deformed wing virus (DWV), the main viral pathogen of the honey bee (Apis mellifera). To test the potential of DWV to be utilized as a vector, the 717 nt sequence coding for the enhanced green fluorescent protein (eGFP), flanked by the peptides targeted by viral protease, was inserted into an infectious cDNA clone of DWV in-frame between the leader protein and the virus structural protein VP2 genes. The in vitro RNA transcripts from egfp-tagged DWV cDNA clones were infectious when injected into honey bee pupae. Stable DWV particles containing genomic RNA of the recovered DWV with egfp inserts were produced, as evidenced by cesium chloride density gradient centrifugation. These particles were infectious to honey bee pupae when injected intra-abdominally. Fluorescent microscopy showed GFP expression in the infected cells and Western blot analysis demonstrated accumulation of free eGFP rather than its fusions with DWV leader protein (LP) and/or viral protein (VP) 2. Analysis of the progeny egfp-tagged DWV showed gradual accumulation of genome deletions for egfp, providing estimates for the rate of loss of a non-essential gene an insect RNA virus genome during natural infection.

Haemolymph removal by Varroa mite destabilizes the dynamical interaction between immune effectors and virus in bees, as predicted by Volterra's model

The association between the deformed wing virus and the parasitic mite Varroa destructor has been identified as a major cause of worldwide honeybee colony losses. The mite acts as a vector of the viral pathogen and can trigger its replication in infected bees. However, the mechanistic details underlying this tripartite interaction are still poorly defined, and, particularly, the causes of viral proliferation in mite-infested bees. Here, we develop and test a novel hypothesis that mite feeding destabilizes viral immune control through the removal of both virus and immune effectors, triggering uncontrolled viral replication. Our hypothesis is grounded on the predator-prey theory developed by Volterra, which predicts prey proliferation when both predators and preys are constantly removed from the system. Consistent with this hypothesis, we show that the experimental removal of increasing volumes of haemolymph from individual bees results in increasing viral densities. By contrast, we do not find consistent support for alternative proposed mechanisms of viral expansion via mite immune suppression or within-host viral evolution. Our results suggest that haemolymph removal plays an important role in the enhanced pathogen virulence observed in the presence of feeding Varroa mites. Overall, these results provide a new model for the mechanisms driving pathogen-parasite interactions in bees, which ultimately underpin honeybee health decline and colony losses.

Genetic Strain Diversity of Multi-Host RNA Viruses that Infect a Wide Range of Pollinators and Associates is Shaped by Geographic Origins

Emerging viruses have caused concerns about pollinator population declines, as multi-host RNA viruses may pose a health threat to pollinators and associated arthropods. In order to understand the ecology and impact these viruses have, we studied their host range and determined to what extent host and spatial variation affect strain diversity. Firstly, we used RT-PCR to screen pollinators and associates, including honey bees (Apis mellifera) and invasive Argentine ants (Linepithema humile), for virus presence and replication. We tested for the black queen cell virus (BQCV), deformed wing virus (DWV), and Kashmir bee virus (KBV) that were initially detected in bees, and the two recently discovered Linepithema humile bunya-like virus 1 (LhuBLV1) and Moku virus (MKV). DWV, KBV, and MKV were detected and replicated in a wide range of hosts and commonly co-infected hymenopterans. Secondly, we placed KBV and DWV in a global phylogeny with sequences from various countries and hosts to determine the association of geographic origin and host with shared ancestry. Both phylogenies showed strong geographic rather than host-specific clustering, suggesting frequent inter-species virus transmission. Transmission routes between hosts are largely unknown. Nonetheless, avoiding the introduction of non-native species and diseased pollinators appears important to limit spill overs and disease emergence.

Foodborne Transmission and Clinical Symptoms of Honey Bee Viruses in Ants Lasius spp

Emerging infectious diseases are often the products of host shifts, where a pathogen jumps from its original host to a novel species. Viruses in particular cross species barriers frequently. Acute bee paralysis virus (ABPV) and deformed wing virus (DWV) are viruses described in honey bees (Apis mellifera) with broad host ranges. Ants scavenging on dead honey bees may get infected with these viruses via foodborne transmission. However, the role of black garden ants, Lasius niger and Lasius platythorax, as alternative hosts of ABPV and DWV is not known and potential impacts of these viruses have not been addressed yet. In a laboratory feeding experiment, we show that L. niger can carry DWV and ABPV. However, negative-sense strand RNA, a token of virus replication, was only detected for ABPV. Therefore, additional L. niger colonies were tested for clinical symptoms of ABPV infections. Symptoms were detected at colony (fewer emerging workers) and individual level (impaired locomotion and movement speed). In a field survey, all L. platythorax samples carried ABPV, DWV-A and -B, as well as the negative-sense strand RNA of ABPV. These results show that L. niger and L. platythorax are alternative hosts of ABPV, possibly acting as a biological vector of ABPV and as a mechanical one for DWV. This is the first study showing the impact of honey bee viruses on ants. The common virus infections of ants in the field support possible negative consequences for ecosystem functioning due to host shifts.

The detection of honey bee (Apis mellifera)-associated viruses in ants

Interspecies virus transmission involving economically important pollinators, including honey bees (Apis mellifera), has recently sparked research interests regarding pollinator health. Given that ants are common pests within apiaries in the southern U.S., the goals of this study were to (1) survey ants found within or near managed honey bee colonies, (2) document what interactions are occurring between ant pests and managed honey bees, and 3) determine if any of six commonly occurring honey bee-associated viruses were present in ants collected from within or far from apiaries. Ants belonging to 14 genera were observed interacting with managed honey bee colonies in multiple ways, most commonly by robbing sugar resources from within hives. We detected at least one virus in 89% of the ant samples collected from apiary sites (n = 57) and in 15% of ant samples collected at non-apiary sites (n = 20). We found that none of these ant samples tested positive for the replication of Deformed wing virus, Black queen cell virus, or Israeli acute paralysis virus, however. Future studies looking at possible virus transmission between ants and bees could determine whether ants can be considered mechanical vectors of honey bee-associated viruses, making them a potential threat to pollinator health.

Foodborne Transmission of Deformed Wing Virus to Ants (Myrmica rubra)

Virus host shifts occur frequently, but the whole range of host species and the actual transmission pathways are often poorly understood. Deformed wing virus (DWV), an RNA virus described from honeybees (Apis mellifera), has been shown to have a broad host range. Since ants are often scavenging on dead honeybees, foodborne transmission of these viruses may occur. However, the role of the ant Myrmica rubra as an alternative host is not known and foodborne transmission to ants has not been experimentally addressed yet. Here, we show with a 16-week feeding experiment that foodborne transmission enables DWV type-A and -B to infect M. rubra and that these ants may serve as a virus reservoir. However, the titers of both plus- and minus-sense viral RNA strands decreased over time. Since the ants were fed with highly virus-saturated honeybee pupae, this probably resulted in initial viral peaks, then approaching lower equilibrium titers in infected individuals later. Since DWV infections were also found in untreated field-collected M. rubra colonies, our results support the wide host range of DWV and further suggest foodborne transmission as a so far underestimated spread mechanism.

Evidence of presence and replication of honey bee viruses among wild bee pollinators in subtropical environments

We determined the presence of six viruses in different bee species collected in subtropical environments. Deformed wing virus (DWV) and black queen cell virus (BQCV) were detected in >90% of honey bee samples and in 50-100% of four stingless bee, two bumble bee and one solitary bee species. Additionally, minus DWV and BQCV RNA strands were detected, indicating that the viruses replicate in several hosts. This is the first report of honey bee viruses replicating in six wild bee species in the tropics. If pathogenic to them, viral infections could result in negative impacts in agricultural and unmanaged ecosystems.

Viral infections alter antennal epithelium ultrastructure in honey bees

Varroa destructor and its associated viruses, in particular deformed wing virus (DWV), have been identified as probable causes of honey bee (Apis mellif era L.) colony losses. Evidence suggests that elevated DWV titres in bees could compromise sensory and communication abilities resulting in negative consequences for hygienic behaviour. As antennae play a central role in this behaviour, we compared antennal ultrastructure in DWV-symptomatic and asymptomatic bees. The results show that virus capsids accumulate in the basal regions of the antennal epithelium, close to the haemolymph. No virus particles were detected at the level of sensory sensilla, such as pore plates, nor within the sensory cell dendrites associated with these sensilla. However, membranous structures appeared to be more prevalent in supporting cells surrounding the dendrites of DWV-symptomatic bees. Para-crystalline arrays containing large numbers of virus particles were detected in the antennae of DWV-symptomatic bees but not in asymptomatic bees.

Feeding by Tropilaelaps mercedesae on pre- and post-capped brood increases damage to Apis mellifera colonies

Tropilaelaps mercedesae parasitism can cause Apis mellifera colony mortality in Asia. Here, we report for the first time that tropilaelaps mites feed on both pre- and post-capped stages of honey bees. Feeding on pre-capped brood may extend their survival outside capped brood cells, especially in areas where brood production is year-round. In this study, we examined the types of injury inflicted by tropilaelaps mites on different stages of honey bees, the survival of adult honey bees, and level of honey bee viruses in 4th instar larvae and prepupae. The injuries inflicted on different developing honey bee stages were visualised by staining with trypan blue. Among pre-capped stages, 4th instar larvae sustained the highest number of wounds (4.6 ± 0.5/larva) while 2nd-3rd larval instars had at least two wounds. Consequently, wounds were evident on uninfested capped brood (5th-6th instar larvae = 3.91 ± 0.64 wounds; prepupae = 5.25 ± 0.73 wounds). Tropilaelaps mite infestations resulted in 3.4- and 6-fold increases in the number of wounds in 5th-6th instar larvae and prepupae as compared to uninfested capped brood, respectively. When wound-inflicted prepupae metamorphosed to white-eyed pupae, all wound scars disappeared with the exuviae. This healing of wounds contributed to the reduction of the number of wounds (≤10) observed on the different pupal stages. Transmission of mite-borne virus such as Deformed Wing Virus (DWV) was also enhanced by mites feeding on early larval stages. DWV and Black Queen Cell Virus (BQCV) were detected in all 4th instar larvae and prepupae analysed. However, viral levels were more pronounced in scarred 4th instar larvae and infested prepupae. The remarkably high numbers of wounds and viral load on scarred or infested developing honey bees may have caused significant weight loss and extensive injuries observed on the abdomen, wings, legs, proboscis and antennae of adult honey bees. Together, the survival of infested honey bees was significantly compromised. This study demonstrates the ability of tropilaelaps mites to inflict profound damage on A. mellifera hosts. Effective management approaches need to be developed to mitigate tropilaelaps mite problems.

Deformed wing virus type A, a major honey bee pathogen, is vectored by the mite Varroa destructor in a non-propagative manner

Honey bees, the primary managed insect pollinator, suffer considerable losses due to Deformed wing virus (DWV), an RNA virus vectored by the mite Varroa destructor. Mite vectoring has resulted in the emergence of virulent DWV variants. The basis for such changes in DWV is poorly understood. Most importantly, it remains unclear whether replication of DWV occurs in the mite. In this study, we exposed Varroa mites to DWV type A via feeding on artificially infected honey bees. A significant, 357-fold increase in DWV load was observed in these mites after 2 days. However, after 8 additional days of passage on honey bee pupae with low viral loads, the DWV load dropped by 29-fold. This decrease significantly reduced the mites’ ability to transmit DWV to honey bees. Notably, negative-strand DWV RNA, which could indicate viral replication, was detected only in mites collected from pupae with high DWV levels but not in the passaged mites. We also found that Varroa mites contain honey bee mRNAs, consistent with the acquisition of honey bee cells which would additionally contain DWV replication complexes with negative-strand DWV RNA. We propose that transmission of DWV type A by Varroa mites occurs in a non-propagative manner.

Transcriptomic responses to diet quality and viral infection in Apis mellifera

BACKGROUND

Parts of Europe and the United States have witnessed dramatic losses in commercially managed honey bees over the past decade to what is considered an unsustainable extent. The large-scale loss of bees has considerable implications for the agricultural economy because bees are one of the leading pollinators of numerous crops. Bee declines have been associated with several interactive factors. Recent studies suggest nutritional and pathogen stress can interactively contribute to bee physiological declines, but the molecular mechanisms underlying interactive effects remain unknown. In this study, we provide insight into this question by using RNA-sequencing to examine how monofloral diets and Israeli acute paralysis virus inoculation influence gene expression patterns in bees.

RESULTS

We found a considerable nutritional response, with almost 2000 transcripts changing with diet quality. The majority of these genes were over-represented for nutrient signaling (insulin resistance) and immune response (Notch signaling and JaK-STAT pathways). In our experimental conditions, the transcriptomic response to viral infection was fairly limited. We only found 43 transcripts to be differentially expressed, some with known immune functions (argonaute-2), transcriptional regulation, and muscle contraction. We created contrasts to explore whether protective mechanisms of good diet were due to direct effects on immune function (resistance) or indirect effects on energy availability (tolerance). A similar number of resistance and tolerance candidate differentially expressed genes were found, suggesting both processes may play significant roles in dietary buffering from pathogen infection.

CONCLUSIONS

Through transcriptional contrasts and functional enrichment analysis, we contribute to our understanding of the mechanisms underlying feedbacks between nutrition and disease in bees. We also show that comparing results derived from combined analyses across multiple RNA-seq studies may allow researchers to identify transcriptomic patterns in bees that are concurrently less artificial and less noisy. This work underlines the merits of using data visualization techniques and multiple datasets to interpret RNA-sequencing studies.

DWV-A Lethal to Honey Bees (Apis mellifera): A Colony Level Survey of DWV Variants (A, B, and C) in England, Wales, and 32 States across the US

The strong association between Varroa destructor, deformed wing virus (DWV), and high overwintering colony losses (OCL) of honey bees is well established. Three DWV master variants (DWV-A, -B, and -C) have been described, and their role in colony mortality remains an open question. Therefore, the aim of this study is to investigate the seasonal prevalence, viral load, and changing distribution of the three DWV master variants within honey bee colonies from England, Wales, and 32 states across the United States. Here, we report that in 2016, DWV-B was prevalent (100%, n = 249) and dominant (95%) in England and Wales, compared to the US. (56%, n = 217 and 23%, respectively), where DWV-A was prevalent (83%, n = 217) and dominant (63%). DWV-C was regularly detected in low viral loads (<1 × 10⁷ genome equivalents per bee) and at lower prevalence (58% in England and Wales, n = 203, and 14% across the United States, n = 124) compared to DWV-A and -B. DWV-B prevalence and dominance in England and Wales coincided with low OCL (6%). Meanwhile, a 60% loss was reported by participating U.S. beekeepers. In the United States, DWV-A prevalence (89%, n = 18) and viral load were significantly (p = 0.002) higher (1 × 10 ⁸–1 × 10¹¹) in colonies that died when compared to the surviving colonies (49% (n = 27), 1 × 10⁶–1 × 10¹⁰). DWV-B had low prevalence (56%, n = 18) in the colonies that died with viral loads of <1 × 10¹⁰. However, DWV-B was routinely detected in high viral loads (>1 × 10¹⁰) in surviving colonies from all sample locations, providing further supporting evidence of DWV-A exhibiting increased virulence over DWV-B at the colony level.

Honey Bee Parasitic Mite Contains the Sensilla-Rich Sensory Organ on the Foreleg Tarsus Expressing Ionotropic Receptors With Conserved Functions

Honey bee parasitic mites (Tropilaelaps mercedesae and Varroa destructor) detect temperature, humidity, and odor but the underlying sensory mechanisms are poorly understood. To uncover how T. mercedesae responds to environmental stimuli inside a hive, we first identified the sensilla-rich sensory organ on the foreleg tarsus. The organ appeared to correspond to Haller's organ in ticks and contained four types of sensilla, which may respond to different stimuli based on their morphology. We searched for differentially expressed genes between the forelegs and hindlegs to identify mRNAs potentially associated with the sensory organ. The forelegs were enriched with mRNAs encoding sensory proteins such as ionotropic receptors (IRs) and gustatory receptors, as well as proteins involved in ciliary transport. We also found that T. mercedesae IR25a and IR93a were capable of rescuing temperature and humidity preference defects in Drosophila melanogaster IR25a and IR93a mutants. These results demonstrate that the structures and physiological functions of ancient IRs have been conserved during arthropod evolution. Our study provides insight into the sensory mechanisms of honey bee parasitic mites, as well as potential targets for methods to control the most serious honey bee pest.

Impact of Varroa destructor and deformed wing virus on emergence, cellular immunity, wing integrity and survivorship of Africanized honey bees in Mexico

The ectoparasitic mite Varroa destructor is the primary health problem of honey bees (Apis mellifera) worldwide. Africanized honey bees in Brazil have demonstrated tolerance to the mite, but there is controversy about the degree of mite tolerance of Africanized bees in other countries. This study was conducted to quantify the effect of V. destructor parasitism on emergence, hemocyte concentration, wing integrity and longevity of Africanized honey bees in Mexico. Africanized bee brood were artificially infested with V. destructor mites and held in an incubator until emergence as adults and compared to non-infested controls. Deformed wing virus (DWV) presence was determined in the mites used to infest the bees. After emergence, the bees were maintained in an incubator to determine survivorship. The percentage of worker bees that emerged from parasitized cells (69%) was significantly lower than that of bees emerged from non-infested cells (96%). Newly-emerged parasitized bees had a significantly lower concentration of hemocytes in the hemolymph than non-parasitized bees. Additionally, the proportion of bees with deformed wings that emerged from V. destructor-parasitized cells was significantly higher (54%) than that of the control group (0%). The mean survival time of bees that emerged from infested and non-infested cells was 8.5 ± 0.3 and 14.4 ± 0.4 days, respectively, and the difference was significant. We conclude that V. destructor parasitism and DWV infections kill, cause deformities and inhibit cellular immunity in developing Africanized honey bees, and significantly reduce the lifespan of adult bees in Mexico. These results suggest that the tolerance of Africanized bees to V. destructor is related to adult bee mechanisms.

Disentangling host-parasite-pathogen interactions in a varroa-resistant honeybee population reveals virus tolerance as an independent, naturally adapted survival mechanism

The ectoparasitic mite, Varroa destructor, is unarguably the leading cause of honeybee (Apis mellifera) mortality worldwide through its role as a vector for lethal viruses, in particular, strains of the Deformed wing virus (DWV) and Acute bee paralysis virus (ABPV) complexes. This multi-level system of host-parasite-pathogen interactions makes it difficult to investigate effects of either the mite or the virus on natural host survival. The aim of this study was to remove confounding effects of varroa to examine the role of virus susceptibility in the enhanced survival of a naturally adapted Swedish mite-resistant (MR) honeybee population, relative to mite-susceptible (MS) honeybees. Caged adult bees and laboratory reared larvae, from varroa-free colonies, were inoculated with DWV and ABPV in a series of feeding infection experiments, while control groups received virus-free food. Virus infections were monitored using RT-qPCR assays in individuals sampled over a time course. In both adults and larvae the DWV and ABPV infection dynamics were nearly identical between MR and MS groups, but MS adults suffered significantly higher mortality than MR adults. Results suggest virus tolerance, rather than reduced susceptibility or virus resistance, is an important component of the natural survival of this honeybee population.

Deformed Wing Virus Infection in Honey Bees (Apis mellifera L.)

Deformed wing virus (DWV) is a single-stranded RNA virus of honey bees (Apis mellifera L.) transmitted by the parasitic mite Varroa destructor. Although DWV represents a major threat to honey bee health worldwide, the pathological basis of DWV infection is not well documented. The objective of this study was to investigate clinicopathological and histological aspects of natural DWV infection in honey bee workers. Emergence of worker honey bees was observed in 5 colonies that were clinically affected with DWV and the newly emerged bees were collected for histopathology. DWV-affected bees were 2 times slower to emerge and had 30% higher mortality compared to clinically normal bees. Hypopharyngeal glands in bees with DWV were hypoplastic, with fewer intracytoplasmic secretory vesicles; cells affected by apoptosis were observed more frequently. Mandibular glands were hypoplastic and were lined by cuboidal epithelium in severely affected bees compared to tall columnar epithelium in nonaffected bees. The DWV load was on average 1.7 × 10⁶ times higher (P < .001) in the severely affected workers compared to aged-matched sister honey bee workers that were not affected by deformed wing disease based on gross examination. Thus, DWV infection is associated with prolonged emergence, increased mortality during emergence, and hypoplasia of hypopharyngeal and mandibular glands in newly emerged worker honey bees in addition to previously reported deformed wing abnormalities.

The Two Prevalent Genotypes of an Emerging Infectious Disease, Deformed Wing Virus, Cause Equally Low Pupal Mortality and Equally High Wing Deformities in Host Honey Bees

Deformed wing virus (DWV) is an emerging infectious disease of the honey bee (Apis mellifera) that is considered a major cause of elevated losses of honey bee colonies. DWV comprises two widespread genotypes: the originally described genotype A, and genotype B. In adult honey bees, DWV-B has been shown to be more virulent than DWV-A. However, their comparative effects on earlier host developmental stages are unknown. Here, we experimentally inoculated honey bee pupae and tested for the relative impact of DWV-A versus DWV-B on mortality and wing deformities in eclosing adults. DWV-A and DWV-B caused similar, and only slightly elevated, pupal mortality (mean 18% greater mortality than control). Both genotypes caused similarly high wing deformities in eclosing adults (mean 60% greater wing deformities than control). Viral titer was high in all of the experimentally inoculated eclosing adults, and was independent of wing deformities, suggesting that the phenotype 'deformed wings' is not directly related to viral titer or viral genotype. These viral traits favor the emergence of both genotypes of DWV by not limiting the reproduction of its vector, the ectoparasitic Varroa destructor mite, in infected pupae, and thereby facilitating the spread of DWV in honey bees infested by the mite.

Pathogen shifts in a honeybee predator following the arrival of the Varroa mite

Emerging infectious diseases (EIDs) are a global threat to honeybees, and spillover from managed bees threaten wider insect populations. Deformed wing virus (DWV), a widespread virus that has become emergent in conjunction with the spread of the mite Varroa destructor, is thought to be partly responsible for global colony losses. The arrival of Varroa in honeybee populations causes a dramatic loss of viral genotypic diversity, favouring a few virulent strains. Here, we investigate DWV spillover in an invasive Hawaiian population of the wasp, Vespula pensylvanica, a honeybee predator and honey-raider. We show that Vespula underwent a parallel loss in DWV variant diversity upon the arrival of Varroa, despite the mite being a honeybee specialist. The observed shift in Vespula DWV and the variant-sharing between Vespula and Apis suggest that these wasps can acquire DWV directly or indirectly from honeybees. Apis prey items collected from Vespula foragers were positive for DWV, indicating predation is a possible route of transmission. We also sought cascading effects of DWV shifts in a broader Vespula pathogen community. We identified concurrent changes in a suite of additional pathogens, as well as shifts in the associations between these pathogens in Vespula. These findings reveal how hidden effects of the Varroa mite can, via spillover, transform the composition of pathogens in interacting species, with potential knock-on effects for entire pathogen communities.

Putative Drone Copulation Factors Regulating Honey Bee (Apis mellifera) Queen Reproduction and Health: A Review

Honey bees are major pollinators of agricultural and non-agricultural landscapes. In recent years, honey bee colonies have exhibited high annual losses and commercial beekeepers frequently report poor queen quality and queen failure as the primary causes. Honey bee colonies are highly vulnerable to compromised queen fertility, as each hive is headed by one reproductive queen. Queens mate with multiple drones (male bees) during a single mating period early in life in which they obtain enough spermatozoa to fertilize their eggs for the rest of their reproductive life span. The process of mating initiates numerous behavioral, physiological, and molecular changes that shape the fertility of the queen and her influence on the colony. For example, receipt of drone semen can modulate queen ovary activation, pheromone production, and subsequent worker retinue behavior. In addition, seminal fluid is a major component of semen that is primarily derived from drone accessory glands. It also contains a complex mixture of proteins such as proteases, antioxidants, and antimicrobial proteins. Seminal fluid proteins are essential for inducing post-mating changes in other insects such as Drosophila and thus they may also impact honey bee queen fertility and health. However, the specific molecules in semen and seminal fluid that initiate post-mating changes in queens are still unidentified. Herein, we summarize the mating biology of honey bees, the changes queens undergo during and after copulation, and the role of drone semen and seminal fluid in post-mating changes in queens. We then review the effects of seminal fluid proteins in insect reproduction and potential roles for honey bee drone seminal fluid proteins in queen reproduction and health. We finish by proposing future avenues of research. Further elucidating the role of drone fertility in queen reproductive health may contribute towards reducing colony losses and advancing honey bee stock development.

Temporal changes in the viromes of Swedish Varroa-resistant and Varroa-susceptible honeybee populations

The parasitic mite, Varroa destructor, in combination with the viruses it vectors, is the main cause for global colony losses of the European honeybee, Apis mellifera. However, an isolated honeybee population established in 1999 on the Island of Gotland, Sweden has naturally acquired resistance to the mite, and has survived without mite control treatment for more than 18 years. A recent study has shown that this mite resistant (MR) population also appears to be resistant to Black queen cell virus (BQCV) and Sacbrood virus (SBV) and tolerant to Deformed wing virus (DWV), relative to nearby mite susceptible (MS) honeybee populations. In this study, RNA sequencing was employed to corroborate these previous findings and identify other viral factors that may play a role in the enhanced survival of this mite resistant honeybee population. Two additional honeybee-infecting viruses, Apis rhabdovirus-1 (ARV-1) and Lake Sinai virus (LSV), were identified and near-complete genomes of these two viruses were obtained. Phylogenetic analyses of the assembled virus sequences revealed consistent separation between the MR and MS honeybee populations, although it is unclear whether this is due to pre-existing differences between the viruses in the two populations when they were established, and isolated, or due to virus genetic adaptation towards reduced virulence in the MR population, to promote colony survival. Reverse transcription quantitative polymerase chain reaction(RT-qPCR) analyses show higher ARV and LSV titres in MS colonies compared to MR colonies, gradually increasing from summer to autumn 2009, and reaching maximum titres in the following spring 2010. While the DWV and BQCV titres in MR colonies increased between autumn 2009 and spring 2010, the SBV practically disappeared entirely by spring 2010. Possible explanations for the apparent virus tolerance or resistance in the Gotland mite-resistant honeybee population are discussed.

In vivo evolution of viral virulence: switching of deformed wing virus between hosts results in virulence changes and sequence shifts

The health of the Western honey bee is threatened by a global epidemic of deformed wing virus (DWV) infections driven by the ectoparasitic mite Varroa destructor acting as mechanical and biological virus vector. Three different variants of DWV, DWV-A, -B and -C exist. Virulence differences between these variants and their relation to V. destructor are still controversially discussed. We performed laboratory experiments to analyze the virulence of DWV directly isolated from crippled bees (DWV_P0 ) or after one additional passage in bee pupae (DWV_P1 ). We demonstrated that DWV_P0 was more virulent than DWV_P1 for pupae, when pupal mortality was taken as virulence marker, and for adult bees, when neurotropism and cognitive impairment were taken as virulence markers. Phylogenetic analysis supported that DWV exists as quasispecies and showed that DWV_P0 clustered with DWV-B and DWV_P1 with DWV-A when the phylogeny was based on the master sequences of the RNA-dependent RNA polymerase but not so when it was based on the VP3 region master sequences. We propose that switching of DWV between the bee and the mite host is accompanied by changes in viral sequence, tissue tropism and virulence and that the RNA-dependent RNA polymerase is involved in determining host range and virulence.

Emerging Viruses in Bees: From Molecules to Ecology

Emerging infectious diseases arise as a result of novel interactions between populations of hosts and pathogens, and can threaten the health and wellbeing of the entire spectrum of biodiversity. Bees and their viruses are a case in point. However, detailed knowledge of the ecological factors and evolutionary forces that drive disease emergence in bees and other host-pathogen communities is surprisingly lacking. In this review, we build on the fundamental insight that viruses evolve and adapt over timescales that overlap with host ecology. At the same time, we integrate the role of host community ecology, including community structure and composition, biodiversity loss, and human-driven disturbance, all of which represent significant factors in bee virus ecology. Both of these evolutionary and ecological perspectives represent major advances but, in most cases, it remains unclear how evolutionary forces actually operate across different biological scales (e.g., from cell to ecosystem). We present a molecule-to-ecology framework to help address these issues, emphasizing the role of molecular mechanisms as key bottom-up drivers of change at higher ecological scales. We consider the bee-virus system to be an ideal one in which to apply this framework. Unlike many other animal models, bees constitute a well characterized and accessible multispecies assemblage, whose populations and interspecific interactions can be experimentally manipulated and monitored in high resolution across space and time to provide robust tests of prevailing theory.

Genetic characterization of a novel picorna-like virus in Culex spp. mosquitoes from Mozambique

BACKGROUND

Mosquitoes are the potential vectors for a variety of viruses that can cause diseases in the human and animal populations. Viruses in the order Picornavirales infect a broad range of hosts, including mosquitoes. In this study, we aimed to characterize a novel picorna-like virus from the Culex spp. of mosquitoes from the Zambezi Valley of Mozambique.

METHODS

The extracted RNA from mosquito pools was pre-amplified with the sequence independent single primer amplification (SISPA) method and subjected to high-throughput sequencing using the Ion Torrent platform. Reads that are classified as Iflaviridae, Picornaviridae and Dicistroviridae were assembled by CodonCode Aligner and SPAdes. Gaps between the viral contigs were sequenced by PCR. The genomic ends were analyzed by 5' and 3' RACE PCRs. The ORF was predicted with the NCBI ORF finder. The conserved domains were identified with ClustalW multiple sequence alignment, and a phylogenetic tree was built with MEGA. The presence of the virus in individual mosquito pools was detected by RT-PCR assay.

RESULTS

A near full-length viral genome (9740 nt) was obtained in Culex mosquitoes that encoded a complete ORF (3112 aa), named Culex picorna-like virus (CuPV-1). The predicted ORF had 38% similarity to the Hubei picorna-like virus 35. The sequence of the conserved domains, Helicase-Protease-RNA-dependent RNA polymerase, were identified by multiple sequence alignment and found to be at the 3' end, similar to iflaviruses. Phylogenetic analysis of the putative RdRP amino acid sequences indicated that the virus clustered with members of the Iflaviridae family. CuPV-1 was detected in both Culex and Mansonia individual pools with low infection rates.

CONCLUSIONS

The study reported a highly divergent, near full-length picorna-like virus genome from Culex spp. mosquitoes from Mozambique. The discovery and characterization of novel viruses in mosquitoes is an initial step, which will provide insights into mosquito-virus interaction mechanisms, genetic diversity and evolution.

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.

Quantitative patterns of vertical transmission of deformed wing virus in honey bees

Deformed wing virus (DWV) is an important pathogen in a broad range of insects, including honey bees. Concordant with the spread of Varroa, DWV is present in the majority of honey bee colonies and can result in either low-level infections with asymptomatic bees that nonetheless exhibit increased colony loss under stress, or high-level infections with acute effects on bee health and viability. DWV can be transmitted vertically or horizontally and evidence suggests that horizontal transmission via Varroa is associated with acute symptomatic infections. Vertical transmission also occurs and is presumably important for the maintenance of DWV in honey bee populations. To further our understanding the vertical transmission of DWV through queens, we performed three experiments: we studied the quantitative effectiveness of vertical transmission, surveyed the prevalence of successful egg infection under commercial conditions, and distinguished among three possible mechanisms of transmission. We find that queen-infection level predicts the DWV titers in their eggs, although the transmission is not very efficient. Our quantitative assessment of DWV demonstrates that eggs in 1/3 of the colonies are infected with DWV and highly infected eggs are rare in newly-installed spring colonies. Additionally, our results indicate that DWV transmission occurs predominantly by virus adhering to the surface of eggs (transovum) rather than intracellularly. Our combined results suggest that the queens' DWV vectoring capacity in practice is not as high as its theoretical potential. Thus, DWV transmission by honey bee queens is part of the DWV epidemic with relevant practical implications, which should be further studied.

A real-time PCR method for quantification of the total and major variant strains of the deformed wing virus

European honey bees (Apis mellifera) are critically important to global food production by virtue of their pollination services but are severely threatened by deformed wing virus (DWV) especially in the presence of the external parasite Varroa destructor. DWV exists as many viral strains with the two major variants (DWV-A and DWV-B) varying in virulence. A single plasmid standard was constructed containing three sections for the specific determination of DWV-A (VP2 capsid region), DWV-B (IRES) and a conserved region suitable for total DWV (helicase region). The assays were confirmed as specific and discriminatory with limits of detections of 25, 25 and 50 genome equivalents for DWV-A, DWV-B and total-DWV, respectively. The methods were successfully tested on Apis mellifera and V. destructor samples with varying DWV profiles. The new method determined a more accurate total DWV titre in samples with substantial DWV-B than the method currently described in the COLOSS Beebook. The proposed assays could be utilized for the screening of large quantities of bee material for both a total DWV load overview along with more detailed investigations into DWV-A and DWV-B profiles.

Presence of Apis Rhabdovirus-1 in Populations of Pollinators and Their Parasites from Two Continents

The viral ecology of bee communities is complex, where viruses are readily shared among co-foraging bee species. Additionally, in honey bees (Apis mellifera), many viruses are transmitted - and their impacts exacerbated - by the parasitic Varroa destructor mite. Thus far, the viruses found to be shared across bee species and transmitted by V. destructor mites are positive-sense single-stranded RNA viruses. Recently, a negative-sense RNA enveloped virus, Apis rhabdovirus-1 (ARV-1), was found in A. mellifera honey bees in Africa, Europe, and islands in the Pacific. Here, we describe the identification - using a metagenomics approach - of ARV-1 in two bee species (A. mellifera and Bombus impatiens) and in V. destructor mites from populations collected in the United States and Israel. We confirmed the presence of ARV-1 in pools of A. mellifera, B. impatiens, and V. destructor from Israeli and U.S. populations by RT-PCR and found that it can reach high titers in individual honey bees and mites (107-108 viral genomic copies per individual). To estimate the prevalence of ARV-1 in honey bee populations, we screened 104 honey bee colonies across Israel, with 21 testing ARV-1-positive. Tagged-primer-mediated RT-PCR analysis detected the presence of the positive-sense ARV-1 RNA in A. mellifera and V. destructor, indicating that ARV-1 replicates in both hosts. This is the first report of the presence of ARV-1 in B. impatiens and of the replication of a rhabdovirus in A. mellifera and V. destructor. Our data suggest that Varroa mites could act as an ARV-1 vector; however, the presence of ARV-1 in B. impatiens (which are not parasitized by Varroa) suggests that it may not require the mite for transmission and ARV-1 may be shared among co-foraging bee species. Given that ARV-1 is found in non-Apis bee species, and because "ARV" is used for the Adelaide River virus, we propose that this virus should be called bee rhabdovirus 1 and abbreviated BRV-1. These results greatly expand our understanding of the diversity of viruses that can infect bee communities, though further analysis is required to determine how infection with this virus impacts these different hosts.

Protein nutrition governs within-host race of honey bee pathogens

Multiple infections are common in honey bees, Apis mellifera, but the possible role of nutrition in this regard is poorly understood. Microsporidian infections, which are promoted by protein-fed, can negatively correlate with virus infections, but the role of protein nutrition for the microsporidian-virus interface is unknown. Here, we challenged naturally deformed wing virus - B (DWV-B) infected adult honey bee workers fed with or without pollen ( = protein) in hoarding cages, with the microsporidian Nosema ceranae. Bee mortality was recorded for 14 days and N. ceranae spore loads and DWV-B titers were quantified. Amongst the groups inoculated with N. ceranae, more spores were counted in protein-fed bees. However, N. ceranae infected bees without protein-diet had reduced longevity compared to all other groups. N. ceranae infection had no effect on protein-fed bee's longevity, whereas bees supplied only with sugar-water showed reduced survival. Our data also support that protein-feeding can have a significant negative impact on virus infections in insects. The negative correlation between N. ceranae spore loads and DWV-B titers was stronger expressed in protein-fed hosts. Proteins not only enhance survival of infected hosts, but also significantly shape the microsporidian-virus interface, probably due to increased spore production and enhanced host immunity.

Immune related gene expression in worker honey bee (Apis mellifera carnica) pupae exposed to neonicotinoid thiamethoxam and Varroa mites (Varroa destructor)

Varroa destructor is one of the most common parasites of honey bee colonies and is considered as a possible co-factor for honey bee decline. At the same time, the use of pesticides in intensive agriculture is still the most effective method of pest control. There is limited information about the effects of pesticide exposure on parasitized honey bees. Larval ingestion of certain pesticides could have effects on honey bee immune defense mechanisms, development and metabolic pathways. Europe and America face the disturbing phenomenon of the disappearance of honey bee colonies, termed Colony Collapse Disorder (CCD). One reason discussed is the possible suppression of honey bee immune system as a consequence of prolonged exposure to chemicals. In this study, the effects of the neonicotinoid thiamethoxam on honey bee, Apis mellifera carnica, pupae infested with Varroa destructor mites were analyzed at the molecular level. Varroa-infested and non-infested honey bee colonies received protein cakes with or without thiamethoxam. Nurse bees used these cakes as a feed for developing larvae. Samples of white-eyed and brown-eyed pupae were collected. Expression of 17 immune-related genes was analyzed by real-time PCR. Relative gene expression in samples exposed only to Varroa or to thiamethoxam or simultaneously to both Varroa and thiamethoxam was compared. The impact from the consumption of thiamethoxam during the larval stage on honey bee immune related gene expression in Varroa-infested white-eyed pupae was reflected as down-regulation of spaetzle, AMPs abaecin and defensin-1 and up-regulation of lysozyme-2. In brown-eyed pupae up-regulation of PPOact, spaetzle, hopscotch and basket genes was detected. Moreover, we observed a major difference in immune response to Varroa infestation between white-eyed pupae and brown-eyed pupae. The majority of tested immune-related genes were upregulated only in brown-eyed pupae, while in white-eyed pupae they were downregulated.

Characterization of the Copy Number and Variants of Deformed Wing Virus (DWV) in the Pairs of Honey Bee Pupa and Infesting Varroa destructor or Tropilaelaps mercedesae

Recent honey bee colony losses, particularly during the winter, have been shown to be associated with the presence of both ectoparasitic mites and Deformed Wing Virus (DWV). Whilst the role of Varroa destructor mites as a viral vector is well established, the role of Tropilaelaps mercedesae mites in viral transmission has not been fully investigated. In this study, we tested the effects that V. destructor and T. mercedesae infestation have on fluctuation of the DWV copy number and alteration of the virus variants in honey bees by characterizing individual pupae and their infesting mites. We observed that both mite species were associated with increased viral copy number in honey bee pupae. We found a positive correlation between DWV copy number in pupae and copy number in infesting mites, and the same DWV type A variant was present in either low or high copy number in both honey bee pupae and infesting V. destructor. These data also suggest that variant diversity is similar between honey bee pupae and the mites that infest them. These results support a previously proposed hypothesis that DWV suppresses the honey bee immune system when virus copy number reaches a specific threshold, promoting greater replication.