First detection and complete genome sequence of Deformed wing virus in Chilean honeybees

Complete genome sequence of deformed wing virus and black queen cell virus isolated from honeybees (Apis mellifera) in Argentina

We report the complete genome sequence of deformed wing virus and black queen cell virus isolated from Argentinean's honeybees. These sequence data will be valuable for future research on the viral variants present in the country and the development of strategies to control the spread of these viruses in apiaries.

Genome analysis and phylogenetic characterization of two deformed wing virus strains from Apis cerana in Vietnam

Background

Deformed wing virus (DWV) is a virulent virus that causes honeybee disease. DWV can exist as a latent infection in honeybees, outbreak into epidemics, and cause serious damage to beekeeping cross the world, including Vietnam.

Methods

The two DWV strains circulating in Vietnamese honeybee, Apis cerana, were first isolated from adult honeybees in North Vietnam (DWV-NVN) and South Vietnam (DWV-SVN). Their complete nucleotide sequences were determined, aligned, and compared with other DWV strains.

Results

The two Vietnamese DWV strains comprised 10,113 bp and contained a large single open reading frame (ORF) of 2,893 amino acids, initiating at nucleotide 1,130 and terminating at nucleotide 9,812. Multiple nucleotide sequence alignment between these two DWV-VN strains and DWV strains in A. mellifera was performed. The DWV-VN strains showed a low genetic identity (from 91.4% to 92.0%) with almost of these strains, but lower identities (89.2% and 89.4%) with UK2 and (89.6%) with the China2 strain. Low identities (91.7% and 91.9%) were also observed between the China3 strain (in A. cerana) and the DWV-VN strains, respectively. The deduced amino acid sequence alignment showed high genetic similarities (97.0%–97.9%) when the USA1, Chile, Italy1, France, UK1, UK2, Japan, Korea2, China1, China2 and China3 strains were compared to the DWV-VN strains. This ratio was 96.7% and 96.8% when the Korea1 strain was compared to the DWV-SVN and DWV-NVN strains, respectively. Numerous amino acid substitutions were identified in the L, VP3, and RdRp sequences. Notably, we observed six substitutions positioned at amino acids 27 (E > I), 98 (S > T), 120 (A > V), 153 (M > T), 170 (D > F), and 174 (Y > F) in the L protein, two amino acid changes at positions 980 (S > A) and 1032 (E > T) in VP3, and one amino acid change at position 2627 (R > C) unique to the DWV-VN strains. Phylogenetic analysis based on complete genome sequences, RdRp sequences and Simplot analysis indicated that there was a significant difference between DWV-VN strains in A. cerana and DWV strains in A. mellifera. The results suggested that the genetic variations of the DWV-VN strains in A. cerana help them to adapt geographical conditions and may lead to change the viral pathogenicity of DWV-VN strains.

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.

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.

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.

Prevalence of honey bee (Apis mellifera) viruses in temperate and subtropical regions from Argentina

In Argentina, bee virus studies are still incipient, and there are no studies regarding the climatic effect. The aim of this study was to assess and compare the presence of honeybee viruses in different climatic regions from Argentina. A total of 385 colonies distributed in five Argentinean eco-regions were examined to evaluate the percentage of infestation with Varroa destructor and the presence of seven virus species (Deformed wing virus, DWV; Acute bee paralysis virus, ABPV; Chronic bee paralysis virus, CBPV; Black queen cell virus, BQCV; Kashmer bee virus, KBV; Israeli acute bee paralysis virus, IAPV; and Sacbrood bee virus, SBV) after honey yield. Two viruses, KBV and IAPV, were not detected. The other five viruses were found in different prevalences: DWV (35%), ABPV (21.5%), BQCV (8.0%), CBPV (2.2%), and SBV (1.1%). We found double and triple viral associations in approximately 25% of the sampled colonies. The mean V. destructor infestation in the colonies prior to the acaricide treatment was 7.12%±8.7%. The knowledge of the prevalence of these viruses in the region and their relation with the mite and other possible influencing factors is important for preventing colony losses. Further studies are necessary to identify the risk factors associated with virus presence and its relationship with other pathogens such as V. destructor.

Infectivity of DWV associated to flower pollen: experimental evidence of a horizontal transmission route

Deformed wing virus (DWV) is a honeybee pathogen whose presence is generally associated with infestation of the colony by the mite Varroa destructor, leading to the onset of infections responsible for the collapse of the bee colony. DWV contaminates bee products such as royal jelly, bee-bread and honey stored within the infected hive. Outside the hive, DWV has been found in pollen loads collected directly from infected as well as uninfected forager bees. It has been shown that the introduction of virus-contaminated pollen into a DWV-free hive results in the production of virus-contaminated food, whose role in the development of infected bees from virus-free eggs has been experimentally demonstrated. The aim of this study was twofold: (i) to ascertain the presence of DWV on pollen collected directly from flowers visited by honeybees and then quantify the viral load and (ii) determine whether the virus associated with pollen is infective. The results of our investigation provide evidence that DWV is present on pollen sampled directly from visited flowers and that, following injection in individuals belonging to the pollinator species Apis mellifera, it is able to establish an active infection, as indicated by the presence of replicating virus in the head of the injected bees. We also provide the first indication that the pollinator species Osmia cornuta is susceptible to DWV infection.

Sequence recombination and conservation of Varroa destructor virus-1 and deformed wing virus in field collected honey bees (Apis mellifera)

We sequenced small (s) RNAs from field collected honeybees (Apis mellifera) and bumblebees (Bombuspascuorum) using the Illumina technology. The sRNA reads were assembled and resulting contigs were used to search for virus homologues in GenBank. Matches with Varroadestructor virus-1 (VDV1) and Deformed wing virus (DWV) genomic sequences were obtained for A. mellifera but not B. pascuorum. Further analyses suggested that the prevalent virus population was composed of VDV-1 and a chimera of 5’-DWV-VDV1-DWV-3’. The recombination junctions in the chimera genomes were confirmed by using RT-PCR, cDNA cloning and Sanger sequencing. We then focused on conserved short fragments (CSF, size > 25 nt) in the virus genomes by using GenBank sequences and the deep sequencing data obtained in this study. The majority of CSF sites confirmed conservation at both between-species (GenBank sequences) and within-population (dataset of this study) levels. However, conserved nucleotide positions in the GenBank sequences might be variable at the within-population level. High mutation rates (Pi>10%) were observed at a number of sites using the deep sequencing data, suggesting that sequence conservation might not always be maintained at the population level. Virus-host interactions and strategies for developing RNAi treatments against VDV1/DWV infections are discussed.