Virus replication in the honey bee parasite, Varroa destructor.

IMPORTANCE: The parasitic mite Varroa destructor is a significant driver of worldwide colony losses of our most important commercial pollinator, the Western honey bee Apis mellifera. Declines in honey bee health are frequently attributed to the viruses that mites vector to honey bees, yet whether mites passively transmit viruses as a mechanical vector or actively participate in viral amplification and facilitate replication of honey bee viruses is debated. Our work investigating the antiviral RNA interference response in V. destructor demonstrates that key viruses associated with honey bee declines actively replicate in mites, indicating that they are biological vectors, and the host range of bee-associated viruses extends to their parasites, which could impact virus evolution, pathogenicity, and spread.

The final frontier: ecological and evolutionary dynamics of a global parasite invasion.

Studying rapid biological changes accompanying the introduction of alien organisms into native ecosystems can provide insights into fundamental ecological and evolutionary theory. While powerful, this quasi-experimental approach is difficult to implement because the timing of invasions and their consequences are hard to predict, meaning that baseline pre-invasion data are often missing. Exceptionally, the eventual arrival of Varroa destructor (hereafter Varroa) in Australia has been predicted for decades. Varroa is a major driver of honeybee declines worldwide, particularly as vectors of diverse RNA viruses. The detection of Varroa in 2022 at over a hundred sites poses a risk of further spread across the continent. At the same time, careful study of Varroa's spread, if it does become established, can provide a wealth of information that can fill knowledge gaps about its effects worldwide. This includes how Varroa affects honeybee populations and pollination. Even more generally, Varroa invasion can serve as a model for evolution, virology and ecological interactions between the parasite, the host and other organisms.

The role of honey in the ecology of the hive: Nutrition, detoxification, longevity, and protection against hive pathogens.

Honey is the source of energy for the European honey bee, Apis mellifera. Beyond simple nutrition and a hedge against the seasonal, geographic, and chemical unpredictability of nectar, honey has properties that protect the hive against various stresses. Enzyme-mediated detoxification during honey ripening neutralizes potentially toxic phytochemicals, and bees that consume honey have enhanced tolerance to other ingested toxins. Catalase and antioxidant phenolics protect honey bees from oxidative damage caused by reactive oxygen species, promoting their longevity. Phytochemical components of honey and microRNAs have the potential to influence developmental pathways, with diet playing a large role in honey bee caste determination. Components of honey mediate stress response and promote cold tolerance during overwintering. Honey has a suite of antimicrobial mechanisms including osmotic pressure, low water activity, low pH, hydrogen peroxide, and plant-, honey bee-, and microbiota-derived compounds such as phytochemicals and antimicrobial peptides. Certain types of honey, particularly polyfloral honeys, have been shown to inhibit important honey bee pathogens including the bacteria responsible for American and European Foulbrood, the microsporidian Nosema ceranae, and the fungi responsible for Stonebrood. Understanding the diverse functional properties of honey has far-ranging implications for honey bee and hive health and management by beekeepers.

Viral communities in the parasite Varroa destructor and in colonies of their honey bee host (Apis mellifera) in New Zealand.

The parasitic mite Varroa destructor is a leading cause of mortality for Western honey bee (Apis mellifera) colonies around the globe. We sought to confirm the presence and likely introduction of only one V. destructor haplotype in New Zealand, and describe the viral community within both V. destructor mites and the bees that they parasitise. A 1232 bp fragment from mitochondrial gene regions suggests the likely introduction of only one V. destructor haplotype to New Zealand. Seventeen viruses were found in bees. The most prevalent and abundant was the Deformed wing virus A (DWV-A) strain, which explained 95.0% of the variation in the viral community of bees. Black queen cell virus, Sacbrood virus, and Varroa destructor virus 2 (VDV-2) played secondary roles. DWV-B and the Israeli acute paralysis virus appeared absent from New Zealand. Ten viruses were observed in V. destructor, with > 99.9% of viral reads from DWV-A and VDV-2. Substantially more variation in viral loads was observed in bees compared to mites. Where high levels of VDV-2 occurred in mites, reduced DWV-A occurred in both the mites and the bees co-occurring within the same hive. Where there were high loads of DWV-A in mites, there were typically high viral loads in bees.

Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera.

BACKGROUND: Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development. RESULTS: Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination. CONCLUSIONS: We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

Cold case: The disappearance of Egypt bee virus, a fourth distinct master strain of deformed wing virus linked to honeybee mortality in 1970's Egypt.

In 1977, a sample of diseased adult honeybees (Apis mellifera) from Egypt was found to contain large amounts of a previously unknown virus, Egypt bee virus, which was subsequently shown to be serologically related to deformed wing virus (DWV). By sequencing the original isolate, we demonstrate that Egypt bee virus is in fact a fourth unique, major variant of DWV (DWV-D): more closely related to DWV-C than to either DWV-A or DWV-B. DWV-A and DWV-B are the most common DWV variants worldwide due to their close relationship and transmission by Varroa destructor. However, we could not find any trace of DWV-D in several hundred RNA sequencing libraries from a worldwide selection of honeybee, varroa and bumblebee samples. This means that DWV-D has either become extinct, been replaced by other DWV variants better adapted to varroa-mediated transmission, or persists only in a narrow geographic or host range, isolated from common bee and beekeeping trade routes.

A Diverse Viral Community from Predatory Wasps in Their Native and Invaded Range, with a New Virus Infectious to Honey Bees.

Wasps of the genus Vespula are social insects that have become major pests and predators in their introduced range. Viruses present in these wasps have been studied in the context of spillover from honey bees, yet we lack an understanding of the endogenous virome of wasps as potential reservoirs of novel emerging infectious diseases. We describe the characterization of 68 novel and nine previously identified virus sequences found in transcriptomes of Vespula vulgaris in colonies sampled from their native range (Belgium) and an invasive range (New Zealand). Many viruses present in the samples were from the Picorna-like virus family (38%). We identified one Luteo-like virus, Vespula vulgaris Luteo-like virus 1, present in the three life stages examined in all colonies from both locations, suggesting this virus is a highly prevalent and persistent infection in wasp colonies. Additionally, we identified a novel Iflavirus with similarity to a recently identified Moku virus, a known wasp and honey bee pathogen. Experimental infection of honey bees with this novel Vespula vulgaris Moku-like virus resulted in an active infection. The high viral diversity present in these invasive wasps is a likely indication that their polyphagous diet is a rich source of viral infections.

DNA methylation is not a driver of gene expression reprogramming in young honey bee workers.

The presence of DNA methylation marks within genic intervals, also called gene body methylation, is an evolutionarily-conserved epigenetic hallmark of animal and plant methylomes. In social insects, gene body methylation is thought to contribute to behavioural plasticity, for example between foragers and nurse workers, by modulating gene expression. However, recent studies have suggested that the majority of DNA methylation is sequence-specific, and therefore cannot act as a flexible mediator between environmental cues and gene expression. To address this paradox, we examined whole-genome methylation patterns in the brains and ovaries of young honey bee workers that had been subjected to divergent social contexts: the presence or absence of the queen. Although these social contexts are known to bring about extreme changes in behavioral and reproductive traits through differential gene expression, we found no significant differences between the methylomes of workers from queenright and queenless colonies. In contrast, thousands of regions were differentially methylated between colonies, and these differences were not associated with differential gene expression in the subset of genes examined. Methylation patterns were highly similar between brain and ovary tissues and only differed in nine regions. These results strongly indicate that DNA methylation is not a driver of differential gene expression between tissues or behavioral morphs. Finally, despite the lack of difference in methylation patterns, queen presence affected the expression of all four DNA methyltransferase genes, suggesting that these enzymes have roles beyond DNA methylation. Therefore, the functional role of DNA methylation in social insect genomes remains an open question.

Adaptation to vector-based transmission in a honeybee virus.

Global pollinator declines as a result of emerging infectious diseases are of major concern. Managed honeybees Apis mellifera are susceptible to numerous parasites and pathogens, many of which appear to be transmissible to sympatric non-Apis taxa. The ectoparasitic mite Varroa destructor is considered to be the most significant threat to honeybees due to its role in vectoring RNA viruses, particularly Deformed wing virus (DWV). Vector transmission of DWV has resulted in the accumulation of high viral loads in honeybees and is often associated with colony death. DWV has two main genotypes, A and B. DWV-A was more prevalent during the initial phase of V. destructor establishment. In recent years, the global prevalence of DWV-B has increased, suggesting that DWV-B is better adapted to vector transmission than DWV-A. We aimed to determine the role vector transmission plays in DWV genotype prevalence at a colony level. We experimentally increased or decreased the number of V. destructor mites in honeybee colonies, and tracked DWV-A and DWV-B loads over a period of 10 months. Our results show that the two DWV genotypes differ in their response to mite numbers. DWV-A accumulation in honeybees was positively correlated with mite numbers yet DWV-A was largely undetected in the absence of the mite. In contrast, colonies had high loads of DWV-B even when mite numbers were low. DWV-B loads persisted in miticide-treated colonies, indicating that this genotype has a competitive advantage over DWV-A irrespective of mite numbers. Our findings suggest that the global increase in DWV-B prevalence is not driven by selective pressure by the vector. Rather, DWV-B is able to persist in colonies at higher viral loads relative to DWV-A in the presence and absence of V. destructor. The interplay between V. destructor and DWV genotypes within honeybee colonies may have broad consequences upon viral diversity in sympatric taxa as a result of spillover.

Intergenerational transfer of DNA methylation marks in the honey bee.

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.

High-Quality Assemblies for Three Invasive Social Wasps from the Vespula Genus.

Social wasps of the genus Vespula have spread to nearly all landmasses worldwide and have become significant pests in their introduced ranges, affecting economies and biodiversity. Comprehensive genome assemblies and annotations for these species are required to develop the next generation of control strategies and monitor existing chemical control. We sequenced and annotated the genomes of the common wasp (Vespula vulgaris), German wasp (Vespula germanica), and the western yellowjacket (Vespula pensylvanica). Our chromosome-level Vespula assemblies each contain 176-179 Mb of total sequence assembled into 25 scaffolds, with 10-200 unanchored scaffolds, and 16,566-18,948 genes. We annotated gene sets relevant to the applied management of invasive wasp populations, including genes associated with spermatogenesis and development, pesticide resistance, olfactory receptors, immunity and venom. These genomes provide evidence for active DNA methylation in Vespidae and tandem duplications of venom genes. Our genomic resources will contribute to the development of next-generation control strategies, and monitoring potential resistance to chemical control.

A Single Gene Causes Thelytokous Parthenogenesis, the Defining Feature of the Cape Honeybee Apis mellifera capensis.

In honeybees, the ability of workers to produce daughters asexually, i.e., thelytokous parthenogenesis, is restricted to a single subspecies inhabiting the Cape region of South Africa, Apis mellifera capensis. Thelytoky has unleashed new selective pressures and the evolution of traits such as social parasitism, invasiveness, and social cancer. Thelytoky arises from an abnormal meiosis that results in the fusion of two maternal pronuclei, restoring diploidy in newly laid eggs. The genetic basis underlying thelytoky is disputed. To resolve this controversy, we generated a backcross between thelytokous A. m. capensis and non-thelytokous A. m. scutellata from the neighboring population and looked for evidence of genetic markers that co-segregated with thelytokous reproduction in 49 backcross females. We found that markers associated with the gene GB45239 on chromosome 11, including non-synonymous variants, showed consistent co-segregation with thelytoky, whereas no other region did so. Alleles associated with thelytoky were present in all A. m. capensis genomes examined but were absent from all other honeybees worldwide including A. m. scutellata. GB45239 is derived in A. m. capensis and has a putative role in chromosome segregation. It is expressed in ovaries and is downregulated in thelytokous bees, likely because of polymorphisms in the promoter region. Our study reveals how mutations affecting the sequence and/or expression of a single gene can change the reproductive mode of a population.

Accumulation and Competition Amongst Deformed Wing Virus Genotypes in Naïve Australian Honeybees Provides Insight Into the Increasing Global Prevalence of Genotype B.

Honeybee colony deaths are often attributed to the ectoparasitic mite Varroa destructor and deformed wing virus (DWV), vectored by the mite. In the presence of V. destructor both main genotypes (DWV-A and DWV-B) have been correlated with colony loss. Studies show that DWV-B is the most prevalent genotype in the United Kingdom and Europe. More recently DWV-B has increased in prevalence in the United States. The increasing prevalence of DWV-B at the expense of DWV-A suggests that competition exists between the genotypes. Competition may be due to disparities in virulence between genotypes, differences in fitness, such as rate of replication, or a combination of factors. In this study we investigated if DWV genotypes differ in their rate of accumulation in Australian honeybees naïve to both V. destructor and DWV, and if viral load was associated with mortality in honeybee pupae. We singly and co-infected pupae with DWV-A, DWV-B, and a recombinant strain isolated from a V. destructor tolerant bee population. We monitored viral accumulation throughout pupation, up to 192 h post-injection. We found significant differences in accumulation, where DWV-A accumulated to significantly lower loads than DWV-B and the DWV-recombinant. We also found evidence of competition, where DWV-B loads were significantly reduced in the presence of DWV-A, but still accumulated to the highest loads overall. In contrast to previous studies, we found significant differences in virulence between pupae injected with DWV-A and DWV-B. The average mortality associated with DWV-B (0.4% ± 0.33 SE) and DWV-recombinant (2.2% ± 0.83 SE) injection were significantly less than observed for DWV-A (11% ± 1.2 SE). Our results suggest that a higher proportion of DWV-B infected pupae will emerge into adults, compared to DWV-A. Overall, our data suggest that low mortality in pupae and the ability of DWV-B to accumulate to higher loads relative to DWV-A even during co-infection may favor vector transmission by V. destructor, and may thus be contributing factors to the increasing prevalence of DWV-B globally.

Paternally-biased gene expression follows kin-selected predictions in female honey bee embryos.

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent-of-origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA-seq to identify transcripts that are over- or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele.

Unique DNA Methylation Profiles Are Associated with cis-Variation in Honey Bees.

DNA methylation is an important epigenetic modification that mediates diverse processes such as cellular differentiation, phenotypic plasticity, and genomic imprinting. Mounting evidence suggests that local DNA sequence variation can be associated with particular DNA methylation states, indicating that the interplay between genetic and epigenetic factors may contribute synergistically to the phenotypic complexity of organisms. Social insects such as ants, bees, and wasps have extensive phenotypic plasticity manifested in their different castes, and this plasticity has been associated with variation in DNA methylation. Yet, the influence of genetic variation on DNA methylation state remains mostly unknown. Here we examine the importance of sequence-specific methylation at the genome-wide level, using whole-genome bisulfite sequencing of the semen of individual honey bee males. We find that individual males harbor unique DNA methylation patterns in their semen, and that genes that are more variable at the epigenetic level are also more likely to be variable at the genetic level. DNA sequence variation can affect DNA methylation by modifying CG sites directly, but can also be associated with local variation in cis that is not CG-site specific. We show that covariation in sequence polymorphism and DNA methylation state contributes to the individual-specificity of epigenetic marks in social insects, which likely promotes their retention across generations, and their capacity to influence evolutionary adaptation.

Direct transmission by injection affects competition among RNA viruses in honeybees.

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.

Viable Triploid Honey Bees (Apis mellifera capensis) Are Reliably Produced in the Progeny of CO2 Narcotised Queens.

The haplodiploid system of sex determination of Hymenoptera acts as an exaptation for species to evolve novel forms of asexual reproduction including thelytoky (clonal offspring of the mother). During normal reproduction in Hymenoptera, three of the four products of meiosis that are present in newly-laid eggs are lost as polar bodies, while the remaining pronucleus either develops as a haploid male or fuses with a sperm nucleus to produce a diploid zygote. In contrast, in thelytokous reproduction, which is uncommon but taxonomically widespread, two of the four products of meiosis fuse, as if one acted as a sperm. Queenless workers of Apis mellifera capensis, a subspecies of honey bee from South Africa, routinely reproduce thelytokously. Unmated A. m. capensis queens can also be induced to lay thelytokously by narcosis with carbon dioxide, but mated queens are never thelytokous. We artificially inseminated A. m. capensis queens using CO2 narcosis. Up to 1/3 of offspring workers carried two maternal alleles and an allele of one father whereas no three-allele progeny were seen in control queens of the arrhenotokous (unfertilized eggs result in males) subspecies A. m. scutellata Flow cytometry of three-allele individuals revealed that they were triploid and arose from the fertilization of a thelytokous fusion nucleus. We then reared six queens from a narcotized A. m. capensis queen and determined the ploidy of the offspring queens based on microsatellites. One of the five daughters was triploid. Following artificial insemination, this queen produced unfertilized thelytokous diploid eggs at high frequency, and unfertilized triploid eggs at much lower frequency. If fertilized, thelytokous diploid eggs were non-viable, even though triploidy in itself does not impede normal development. In contrast, when the rarer triploid eggs were fertilized, a proportion developed into viable tetraploids. Our study highlights the extraordinary developmental flexibility of haplo-diploid systems.

A Diverse Range of Novel RNA Viruses in Geographically Distinct Honey Bee Populations.

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.