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.

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.

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.

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.

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.

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.

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.

The dynamic DNA methylation cycle from egg to sperm in the honey bee Apis mellifera.

In honey bees (Apis mellifera), the epigenetic mark of DNA methylation is central to the developmental regulation of caste differentiation, but may also be involved in additional biological functions. In this study, we examine the whole genome methylation profiles of three stages of the haploid honey bee genome: unfertilised eggs, the adult drones that develop from these eggs and the sperm produced by these drones. These methylomes reveal distinct patterns of methylation. Eggs and sperm show 381 genes with significantly different CpG methylation patterns, with the vast majority being more methylated in eggs. Adult drones show greatly reduced levels of methylation across the genome when compared with both gamete samples. This suggests a dynamic cycle of methylation loss and gain through the development of the drone and during spermatogenesis. Although fluxes in methylation during embryogenesis may account for some of the differentially methylated sites, the distinct methylation patterns at some genes suggest parent-specific epigenetic marking in the gametes. Extensive germ line methylation of some genes possibly explains the lower-than-expected frequency of CpG sites in these genes. We discuss the potential developmental and evolutionary implications of methylation in eggs and sperm in this eusocial insect species.

Parent-of-origin effects on genome-wide DNA methylation in the Cape honey bee (Apis mellifera capensis) may be confounded by allele-specific methylation.

BACKGROUND: Intersexual genomic conflict sometimes leads to unequal expression of paternal and maternal alleles in offspring, resulting in parent-of-origin effects. In honey bees reciprocal crosses can show strong parent-of-origin effects, supporting theoretical predictions that genomic imprinting occurs in this species. Mechanisms behind imprinting in honey bees are unclear but differential DNA methylation in eggs and sperm suggests that DNA methylation could be involved. Nonetheless, because DNA methylation is multifunctional, it is difficult to separate imprinting from other roles of methylation. Here we use a novel approach to investigate parent-of-origin DNA methylation in honey bees. In the subspecies Apis mellifera capensis, reproduction of females occurs either sexually by fertilization of eggs with sperm, or via thelytokous parthenogenesis, producing female embryos derived from two maternal genomes. RESULTS: We compared genome-wide methylation patterns of sexually-produced, diploid embryos laid by a queen, with parthenogenetically-produced diploid embryos laid by her daughters. Thelytokous embryos inheriting two maternal genomes had fewer hypermethylated genes compared to fertilized embryos, supporting the prediction that fertilized embryos have increased methylation due to inheritance of a paternal genome. However, bisulfite PCR and sequencing of a differentially methylated gene, Stan (GB18207) showed strong allele-specific methylation that was maintained in both fertilized and thelytokous embryos. For this gene, methylation was associated with haplotype, not parent of origin. CONCLUSIONS: The results of our study are consistent with predictions from the kin theory of genomic imprinting. However, our demonstration of allele-specific methylation based on sequence shows that genome-wide differential methylation studies can potentially confound imprinting and allele-specific methylation. It further suggests that methylation patterns are heritable or that specific sequence motifs are targets for methylation in some genes.

Gene duplication in the major insecticide target site, Rdl, in Drosophila melanogaster.

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.

A parent-of-origin effect on honeybee worker ovary size.

Apis mellifera capensis is unique among honeybees in that unmated workers can produce pseudo-clonal female offspring via thelytokous parthenogenesis. Workers use this ability to compete among themselves and with their queen to be the mother of new queens. Males could therefore enhance their reproductive success by imprinting genes that enhance fertility in their daughter workers. This possibility sets the scene for intragenomic conflict between queens and drones over worker reproductive traits. Here, we show a strong parent-of-origin effect for ovary size (number of ovarioles) in reciprocal crosses between two honeybee subspecies, A. m. capensis and Apis mellifera scutellata. In this cross, workers with an A. m. capensis father had 30% more ovarioles than genotypically matched workers with an A. m. scutellata father. Other traits we measured (worker weight at emergence and the presence/absence of a spermatheca) are influenced more by rearing conditions than by parent-of-origin effects. Our study is the first to show a strong epigenetic (or, less likely, cytoplasmic maternal) effect for a reproductive trait in the honeybee and suggests that a search for parent-of-origin effects in other social insects may be fruitful.

Reproductive interference between honeybee species in artificial sympatry.

Reproductive isolation between closely related species is often incomplete. The Western honeybee, Apis mellifera, and the Eastern hive bee, Apis cerana, have been allopatric for millions of years, but are nonetheless similar in morphology and behaviour. During the last century, the two species were brought into contact anthropogenically, providing potential opportunities for interspecific matings. Hybrids between A. mellifera and A. cerana are inviable, so natural interspecific matings are of concern because they may reduce the viability of A. cerana and A. mellifera populations - two of the world's most important pollinators. We examined the mating behaviour of A. mellifera and A. cerana queens and drones from Caoba Basin, China and Cairns, Australia. Drone mating flight times overlap in both areas. Analysis of the spermathecal contents of queens with species-specific genetic markers indicated that in Caoba Basin, 14% of A. mellifera queens mated with at least one A. cerana male, but we detected no A. cerana queens that had mated with A. mellifera males. Similarly, in Cairns, no A. cerana queens carried A. mellifera sperm, but one-third of A. mellifera queens had mated with at least one A. cerana male. No hybrid embryos were detected in eggs laid by interspecifically mated A. mellifera queens in either location. However, A. mellifera queens artificially inseminated with A. cerana sperm produced inviable hybrid eggs or unfertilized drones. This suggests that reproductive interference will impact the viability of honeybee populations wherever A. cerana and A. mellifera are in contact.

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.

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.

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.

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.