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

The evolutionary significance of epigenetic inheritance is controversial. While epigenetic marks such as DNA methylation can affect gene function and change in response to environmental conditions, their role as carriers of heritable information is often considered anecdotal. Indeed, near-complete DNA methylation reprogramming, as occurs during mammalian embryogenesis, is a major hindrance for the transmission of nongenetic information between generations. Yet it remains unclear how general DNA methylation reprogramming is across the tree of life. Here we investigate the existence of epigenetic inheritance in the honey bee. We studied whether fathers can transfer epigenetic information to their daughters through DNA methylation. We performed instrumental inseminations of queens, each with four different males, retaining half of each male’s semen for whole genome bisulfite sequencing. We then compared the methylation profile of each father’s somatic tissue and semen with the methylation profile of his daughters. We found that DNA methylation patterns were highly conserved between tissues and generations. There was a much greater similarity of methylomes within patrilines (i.e., father-daughter subfamilies) than between patrilines in each colony. Indeed, the samples’ methylomes consistently clustered by patriline within colony. Samples from the same patriline had twice as many shared methylated sites and four times fewer differentially methylated regions compared to samples from different patrilines. Our findings indicate that there is no DNA methylation reprogramming in bees and, consequently, that DNA methylation marks are stably transferred between generations. This points to a greater evolutionary potential of the epigenome in invertebrates than there is in mammals.

The Resistance to Dieldrin gene, Rdl, encodes a GABA-gated chloride channel subunit that is targeted by cyclodiene and phenylpyrazole insecticides. The gene was first characterized in Drosophila melanogaster by genetic mapping of resistance to the cyclodiene dieldrin. The 4,000-fold resistance observed was due to a single amino acid replacement, Ala 301 to Ser. The equivalent change was subsequently identified in Rdl orthologs of a large range of resistant insect species. Here, we report identification of a duplication at the Rdl locus in D. melanogaster. The 113-kb duplication contains one WT copy of Rdl and a second copy with two point mutations: an Ala 301 to Ser resistance mutation and Met 360 to Ile replacement. Individuals with this duplication exhibit intermediate dieldrin resistance compared with single copy Ser 301 homozygotes, reduced temperature sensitivity, and altered RNA editing associated with the resistant allele. Ectopic recombination between Roo transposable elements is involved in generating this genomic rearrangement. The duplication phenotypes were confirmed by construction of a transgenic, artificial duplication integrating the 55.7-kb Rdl locus with a Ser 301 change into an Ala 301 background. Gene duplications can contribute significantly to the evolution of insecticide resistance, most commonly by increasing the amount of gene product produced. Here however, duplication of the Rdl target site creates permanent heterozygosity, providing unique potential for adaptive mutations to accrue in one copy, without abolishing the endogenous role of an essential gene.T he single point mutation in the Resistance to dieldrin (Rdl) gene represents one of the most significant cases of target site resistance to an insecticide yet observed. Cyclodiene resistance was reported in 62% of insecticide resistant species in the 1980s, following widespread use of cyclodiene insecticides, including dieldrin, which started in the 1950s (1). The nature of the genetic target, Rdl, was discovered after dieldrin was discontinued because of the widespread evolution of resistance in many species. Rdl was first discovered in Drosophila melanogaster using a positional cloning approach. High homology to human GABA receptors confirmed it was the first insect ligand-gated chloride channel subunit identified (2-4). A point mutation in the chloride channel pore-lining domain, replacing alanine 301 with serine, was present in all resistant D. melanogaster strains (5). This mutation provided 4,000-fold resistance when homozygous and lower levels of resistance in heterozygotes (2, 6). The homologous mutation was subsequently found in a large number of cyclodiene-resistant species from many insect orders (7-9), as well as a glycine replacement at the homologous site in some resistant strains of Drosophila simulans and other species (5, 10, 11).Characterization of deficiency lines and inversions in D. melanogaster showed that Rdl is an essential gene (3). Thus, the Ala 301 to Ser or Gly mutation in Rdl exhibits unique properties, p...

The arrival of the ectoparasitic mite Varroa destructor on the western honeybee Apis mellifera saw a change in the diversity and prevalence of honeybee RNA viruses. One virus in particular, deformed wing virus (DWV) has become closely associated with V. destructor , leading many to conclude that V. destructor has affected viral virulence by changing the mode of transmission. While DWV is normally transmitted via feeding and faeces, V. destructor transmits viruses by direct injection. This change could have resulted in higher viral prevalence causing increased damage to the bees. Here we test the effect of a change in the mode of transmission on the composition and levels of honeybee RNA viruses in the absence of V. destructor . We find a rapid increase in levels of two viruses, sacbrood virus (SBV) and black queen cell virus (BQCV) after direct injection of viral extracts into honeybee pupae. In pupae injected with high levels of DWV extracted from symptomatic adult bees, DWV levels rapidly decline in the presence of SBV and BQCV. Further, we observe high mortality in honeybee pupae when injected with SBV and BQCV, whereas injecting pupae with high levels of DWV results in near 100% survival. Our results suggest a different explanation for the observed association between V. destructor and DWV. Instead of V. destructor causing an increase in DWV virulence, we hypothesize that direct virus inoculation, such as that mediated by a vector, quickly eliminates the most virulent honeybee viruses resulting in an association with less virulent viruses such as DWV.