Immunogene and viral transcript dynamics during parasitic Varroa destructor mite infection of developing honey bee (Apis mellifera) pupae

Functional redundancy and niche specialization in honeybee and Varroa microbiomes

The honeybee (Apis mellifera) is a key pollinator critical to global agriculture, facing threats from various stressors, including the ectoparasitic Varroa mite (Varroa destructor). Previous studies have identified shared bacteria between Varroa mites and honeybees, yet it remains unclear if these bacteria assemble similarly in both species. This study builds on existing knowledge by investigating co-occurrence patterns in the microbiomes of both Varroa mites and honeybees, shedding light on potential interactions. Leveraging 16S rRNA datasets, we conducted co-occurrence network analyses, explored Core Association Networks (CAN) and assess network robustness. Comparative network analyses revealed structural differences between honeybee and mite microbiomes, along with shared core features and microbial motifs. The mite network exhibited lower robustness, suggesting less resistance to taxa extension compared to honeybees. Furthermore, analyses of predicted functional profiling and taxa contribution revealed that common central pathways in the metabolic networks have different taxa contributing to Varroa mites and honeybee microbiomes. The results show that while both microbial systems exhibit functional redundancy, in which different taxa contribute to the functional stability and resilience of the ecosystem, there is evidence for niche specialization resulting in unique contributions to specific pathways in each part of this host-parasite system. The specificity of taxa contribution to key pathways offers targeted approaches to Varroa microbiome management and preserving honeybee microbiome. Our findings provide valuable insights into microbial interactions, aiding farmers and beekeepers in maintaining healthy and resilient bee colonies amid increasing Varroa mite infestations.

Deformed wing virus genotypes A and B do not elicit immunologically different responses in naïve honey bee hosts

Iflavirus aladeformis (Picornavirales: Iflaviridae), commonly known as deformed wing virus(DWV), in association with Varroa destructor Anderson and Trueman (Mesostigmata: Varroidae), is a leading factor associated with honey bee (Apis mellifera L. [Hymenoptera: Apidae]) deaths. The virus and mite have a near global distribution, making it difficult to separate the effect of one from the other. The prevalence of two main DWV genotypes (DWV-A and DWV-B) has changed over time, leading to the possibility that the two strains elicit a different immune response by the host. Here, we use a honey bee population naïve to both the mite and the virus to investigate if honey bees show a different immunological response to DWV genotypes. We examined the expression of 19 immune genes by reverse transcription quantitative PCR (RT-qPCR) and analysed small RNA after experimental injection with DWV-A and DWV-B. We found no evidence that DWV-A and DWV-B elicit different immune responses in honey bees. RNA interference genes were up-regulated during DWV infection, and small interfering RNA (siRNA) responses were proportional to viral loads yet did not inhibit DWV accumulation. The siRNA response towards DWV was weaker than the response to another honey bee pathogen, Triatovirus nigereginacellulae (Picornavirales: Dicistroviridae; black queen cell virus), suggesting that DWV is comparatively better at evading host antiviral defences. There was no evidence for the production of virus-derived Piwi-interacting RNAs (piRNAs) in response to DWV. In contrast to previous studies, and in the absence of V. destructor, we found no evidence that DWV has an immunosuppressive effect. Overall, our results advance our understanding of the immunological effect that DWV in isolation elicits in honey bees.

Physiological trade-offs in male social insects: Interactions among infection, immunity, fertility, size, and age in honey bee drones

Female social insects represent a dramatic exception of the evolutionarily conserved physiological trade-off between reproduction and life span, where aging is positively correlated with reproduction. However, whether this facet of life history also pertains to male social insects, remains largely unknown. Male honey bees (drones) die in the act of copulation, placing them under opposing selective pressures. At the individual level, there is inter-male competition for a single successful mating attempt, leading to selective pressure that favors an increase in male fitness. Honey bee drones are haploid individuals and lack the allelic variation in their genome compared to diploid females. We hypothesized that this genetic limitation may result in trade-offs between pathological stress and fitness traits in honey bee males. In our study, we observed differences in size and fertility measures in old and young drones along with stressors of several endemic viruses and the transcriptional immune response. We found that infection does not appear to decrease fertility in old drones, despite evidence for a shift in immune expression away from established mechanisms. Contrary to our expectations, drones additionally do not appear to exhibit a physiological trade-off between size and fertility. These findings demonstrate that drones of different size are likely of different mating quality and that higher quality drones likely favor retaining reproductive output over immune function.

Modeling seasonal immune dynamics of honey bee (Apis mellifera L.) response to injection of heat-killed Serratia marcescens

The honey bee, Apis mellifera L., is one of the main pollinators worldwide. In a temperate climate, seasonality affects the life span, behavior, physiology, and immunity of honey bees. In consequence, it impacts their interaction with pathogens and parasites. In this study, we used Bayesian statistics and modeling to examine the immune response dynamics of summer and winter honey bee workers after injection with the heat-killed bacteria Serratia marcescens, an opportunistic honey bee pathogen. We investigated the humoral and cellular immune response at the transcriptional and functional levels using qPCR of selected immune genes, antimicrobial activity assay, and flow cytometric analysis of hemocyte concentration. Our data demonstrate increased antimicrobial activity at transcriptional and functional levels in summer and winter workers after injection, with a stronger immune response in winter bees. On the other hand, an increase in hemocyte concentration was observed only in the summer bee population. Our results indicate that the summer population mounts a cellular response when challenged with heat-killed S. marcescens, while winter honey bees predominantly rely on humoral immune reactions. We created a model describing the honey bee immune response dynamics to bacteria-derived components by applying Bayesian statistics to our data. This model can be employed in further research and facilitate the investigating of the honey bee immune system and its response to pathogens.

Varroa destructor parasitism and Deformed wing virus infection in honey bees are linked to peroxisome-induced pathways

The ectoparasitic mite Varroa destructor transmits and triggers viral infections that have deleterious effects on honey bee colonies worldwide. We performed a manipulative experiment in which worker bees collected at emergence were exposed to Varroa for 72 h, and their proteomes were compared with those of untreated control bees. Label-free quantitative proteomics identified 77 differentially expressed A. mellifera proteins (DEPs). In addition, viral proteins were identified by orthogonal analysis, and most importantly, Deformed wing virus (DWV) was found at high levels/intensity in Varroa-exposed bees. Pathway enrichment analysis suggested that the main pathways affected included peroxisomal metabolism, cyto-/exoskeleton reorganization, and cuticular proteins. Detailed examination of individual DEPs revealed that additional changes in DEPs were associated with peroxisomal function. In addition, the proteome data support the importance of TGF-β signaling in Varroa-DWV interaction and the involvement of the mTORC1 and Hippo pathways. These results suggest that the effect of DWV on bees associated with Varroa feeding results in aberrant autophagy. In particular, autophagy is selectively modulated by peroxisomes, to which the observed proteome changes strongly corresponded. This study complements previous research with different study designs and suggests the importance of the peroxisome, which plays a key role in viral infections.

Poor Air Quality Is Linked to Stress in Honeybees and Can Be Compounded by the Presence of Disease

Climate change-related extreme weather events have manifested in the western United States as warmer and drier conditions with an increased risk of wildfires. Honeybees, essential for crop pollination in California, are at the center of these extreme weather events. We associated the maximum daily temperature and air quality index values with the performance of colonies placed in wildfire-prone areas and determined the impact of these abiotic stressors on gene expression and histopathology. Our results indicate that poor air quality was associated with higher maximum daily temperatures and a lower gene expression level of Prophenoloxidase (ProPO), which is tied to immune system strength; however, a higher gene expression level of Vitellogenin (Vg) is tied to oxidative stress. There was a positive relationship between Varroa mites and N. ceranae pathogen loads, and a negative correlation between Varroa mites and Heat Shock Protein 70 (HSP70) gene expression, suggesting the limited ability of mite-infested colonies to buffer against extreme temperatures. Histological analyses did not reveal overt signs of interaction between pathology and abiotic stressors, but N. ceranae infections were evident. Our study provides insights into interactions between abiotic stressors, their relation to common biotic stressors, and the expression of genes related to immunity and oxidative stress in bees.

Viral Co-Infections and Antiviral Immunity in Honey Bees

Over the past few decades, honey bees have been facing an increasing number of stressors. Beyond individual stress factors, the synergies between them have been identified as a key factor in the observed increase in colony mortality. However, these interactions are numerous and complex and call for further research. Here, in line with our need for a systemic understanding of the threats that they pose to bee health, we review the interactions between honey bee viruses. As viruses are obligate parasites, the interactions between them not only depend on the viruses themselves but also on the immune responses of honey bees. Thus, we first summarise our current knowledge of the antiviral immunity of honey bees. We then review the interactions between specific pathogenic viruses and their interactions with their host. Finally, we draw hypotheses from the current literature and suggest directions for future research.

Omics-based analysis of honey bee (Apis mellifera) response to Varroa sp. parasitisation and associated factors reveals changes impairing winter bee generation

The extensive annual loss of honey bees (Apis mellifera L.) represents a global problem affecting agriculture and biodiversity. The parasitic mite Varroa destructor, associated with viral co-infections, plays a key role in this loss. Despite years of intensive research, the complex mechanisms of Varroa - honey bee interaction are still not fully defined. Therefore, this study employed a unique combination of transcriptomic, proteomic, metabolomic, and functional analyses to reveal new details about the effect of Varroa mites and naturally associated factors, including viruses, on honey bees. We focused on the differences between Varroa parasitised and unparasitised ten-day-old worker bees collected before overwintering from the same set of colonies reared without anti-mite treatment. Supplementary comparison to honey bees collected from colonies with standard anti-Varroa treatment can provide further insights into the effect of a pyrethroid flumethrin. Analysis of the honey bees exposed to mite parasitisation revealed alterations in the transcriptome and proteome related to immunity, oxidative stress, olfactory recognition, metabolism of sphingolipids, and RNA regulatory mechanisms. The immune response and sphingolipid metabolism were strongly activated, whereas olfactory recognition and oxidative stress pathways were inhibited in Varroa parasitised honey bees compared to unparasitised ones. Moreover, metabolomic analysis confirmed the depletion of nutrients and energy stores, resulting in a generally disrupted metabolism in the parasitised workers. The combined omics-based analysis conducted on strictly parasitised bees revealed the key molecular components and mechanisms underlying the detrimental effects of Varroa sp. and its associated pathogens. This study provides the theoretical basis and interlinked datasets for further research on honey bee response to biological threats and the development of efficient control strategies against Varroa mites.

Honey Bee Genetic Stock Determines Deformed Wing Virus Symptom Severity but not Viral Load or Dissemination Following Pupal Exposure

Honey bees exposed to Varroa mites incur substantial physical damage in addition to potential exposure to vectored viruses such as Deformed wing virus (DWV) that exists as three master variants (DWV-A, DWV-B, and DWV-C) and recombinants. Although mite-resistant bees have been primarily bred to mitigate the impacts of Varroa mites, mite resistance may be associated with increased tolerance or resistance to the vectored viruses. The goal of our study is to determine if five honey bee stocks (Carniolan, Italian, Pol-Line, Russian, and Saskatraz) differ in their resistance or tolerance to DWV based on prior breeding for mite resistance. We injected white-eyed pupae with a sublethal dose (105) of DWV or exposed them to mites and then evaluated DWV levels and dissemination and morphological symptoms upon adult emergence. While we found no evidence of DWV resistance across stocks (i.e., similar rates of viral replication and dissemination), we observed that some stocks exhibited reduced symptom severity suggestive of differential tolerance. However, DWV tolerance was not consistent across mite-resistant stocks as Russian bees were most tolerant, while Pol-Line exhibited the most severe symptoms. DWV variants A and B exhibited differential dissemination patterns that interacted significantly with the treatment group but not bee stock. Furthermore, elevated DWV-B levels reduced adult emergence time, while both DWV variants were associated with symptom likelihood and severity. These data indicate that the genetic differences underlying bee resistance to Varroa mites are not necessarily correlated with DWV tolerance and may interact differentially with DWV variants, highlighting the need for further work on mechanisms of tolerance and bee stock-specific physiological interactions with pathogen variants.

The antiviral role of NF-κB-mediated immune responses and their antagonism by viruses in insects

The antiviral role of innate immune responses mediated by the NF-κB family of transcription factors is well established in vertebrates but was for a long time less clear in insects. Insects encode two canonical NF-κB pathways, the Toll and Imd ('immunodeficiency') pathways, which are best characterised for their role in antibacterial and antifungal defence. An increasing body of evidence has also implicated NF-κB-mediated innate immunity in antiviral responses against some, but not all, viruses. Specific pattern recognition receptors (PRRs) and molecular events leading to NF-κB activation by viral pathogen-associated molecular patterns (PAMPs) have been elucidated for a number of viruses and insect species. Particularly interesting are recent findings indicating that the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway detects viral RNA to activate NF-κB-regulated gene expression. We summarise the literature on virus-NF-κB pathway interactions across the class Insecta, with a focus on the dipterans Drosophila melanogaster and Aedes aegypti. We discuss potential reasons for differences observed between different virus-host combinations, and highlight similarities and differences between cGAS-STING signalling in insects versus vertebrates. Finally, we summarise the increasing number of known molecular mechanisms by which viruses antagonise NF-κB responses, which suggest that NF-κB-mediated immunity exerts strong evolutionary pressures on viruses. These developments in our understanding of insect antiviral immunity have relevance to the large number of insect species that impact on humans through their transmission of human, livestock and plant diseases, exploitation as biotechnology platforms, and role as parasites, pollinators, livestock and pests.

Impact of Honey Bee Migratory Management on Pathogen Loads and Immune Gene Expression is Affected by Complex Interactions With Environment, Worker Life History, and Season

The effects of honey bee management, such as intensive migratory beekeeping, are part of the ongoing debate concerning causes of colony health problems. Even though comparisons of disease and pathogen loads among differently managed colonies indicate some effects, the direct impact of migratory practices on honey bee pathogens is poorly understood. To test long- and short-term impacts of managed migration on pathogen loads and immunity, experimental honey bee colonies were maintained with or without migratory movement. Individuals that experienced migration as juveniles (e.g., larval and pupal development), as adults, or both were compared to control colonies that remained stationary and therefore did not experience migratory relocation. Samples at different ages and life-history stages (hive bees or foragers), taken at the beginning and end of the active season, were analyzed for pathogen loads and physiological markers of health. Bees exposed to migratory management during adulthood had increased levels of the AKI virus complex (Acute bee paralysis, Kashmir bee, and Israeli acute bee paralysis viruses) and decreased levels of antiviral gene expression (dicer-like). However, those in stationary management as adults had elevated gut parasites (i.e. trypanosomes). Effects of environment during juvenile development were more complex and interacted with life-history stage and season. Age at collection, life-history stage, and season all influenced numerous factors from viral load to immune gene expression. Although the factors that we examined are not independent, the results illuminate potential factors in both migratory and nonmigratory beekeeping that are likely to contribute to colony stress, and also indicate potential mitigation measures.

A SNP assay for assessing diversity in immune genes in the honey bee (Apis mellifera L.)

With a growing number of parasites and pathogens experiencing large-scale range expansions, monitoring diversity in immune genes of host populations has never been so important because it can inform on the adaptive potential to resist the invaders. Population surveys of immune genes are becoming common in many organisms, yet they are missing in the honey bee (Apis mellifera L.), a key managed pollinator species that has been severely affected by biological invasions. To fill the gap, here we identified single nucleotide polymorphisms (SNPs) in a wide range of honey bee immune genes and developed a medium-density assay targeting a subset of these genes. Using a discovery panel of 123 whole-genomes, representing seven A. mellifera subspecies and three evolutionary lineages, 180 immune genes were scanned for SNPs in exons, introns (< 4 bp from exons), 3' and 5´UTR, and < 1 kb upstream of the transcription start site. After application of multiple filtering criteria and validation, the final medium-density assay combines 91 quality-proved functional SNPs marking 89 innate immune genes and these can be readily typed using the high-sample-throughput iPLEX MassARRAY system. This medium-density-SNP assay was applied to 156 samples from four countries and the admixture analysis clustered the samples according to their lineage and subspecies, suggesting that honey bee ancestry can be delineated from functional variation. In addition to allowing analysis of immunogenetic variation, this newly-developed SNP assay can be used for inferring genetic structure and admixture in the honey bee.

Tropilaelaps mercedesae parasitism changes behavior and gene expression in honey bee workers

Tropilaelaps mercedesae is one of the most problematic honey bee parasites and has become more threatening to the beekeeping industry. Tropilaelaps can easily parasitize immature honey bees (larvae and pupae) and have both lethal and sublethal effects on the individual worker bees. Our study for the first time experimentally assessed the effects of T. mercedesae on olfactory learning, flight ability, homing ability as well as transcriptional changes in parasitized adult honey bees. T. mercedesae infestation had negative impacts on olfactory associated function, flight ability, and homing rate. The volume of the mushroom body significantly increased in infested honey bees, which may be correlated to the lower sucrose responsiveness as well as lower learning ability in the infested bees. The gene expression involved in immune systems and carbohydrate transport and metabolism were significantly different between infested bees and non-infested bees. Moreover, genes function in cell adhesion play an essential role in olfactory sensory in honey bees. Our findings provide a comprehensive understanding of European honey bees in response to T. mercedesae infestation, and could be used to further investigate the complex molecular mechanisms in honey bees under parasitic stress.

In recent decades, there has been serious concern about the decline of honey bees in the world. One of the most serious factors contributing to bee population declines is mite parasitism. Although Varroa destructor is the most widespread globally, Tropilaelaps mercedesae displays greater threat to bee colonies due to its smaller size, shorter phoretic phase, more rapid locomotion, as well as faster reproductive rate. Tropilaelaps mites, originally parasite of the giant Asian honey bees, now becoming an emerging threat of European honey bees (Apis mellifera) in Asian area. This work aimed to investigate the influence of T. mercedesae infestation on behavior and gene expression in A. mellifera. Our results highlight the T. mercedesae infestation induced negative effects of olfactory learning, flight ability, homing ability of honey bee workers. Moreover, we found that T. mercedesae infestation caused the up-regulation of genes involved in immune systems and carbohydrate mechanism which were correlated to the different olfactory learning performance in infested honeybee. In addition, genes function in cell adhesion play an essential role in olfactory sensory in honey bees. Our results increase the knowledge of proximate mechanisms in honey bee responding to parasitic stress.

Transcriptome Profiling Reveals a Novel Mechanism of Antiviral Immunity Upon Sacbrood Virus Infection in Honey Bee Larvae (Apis cerana)

The honey bee is one of the most important pollinators in the agricultural system and is responsible for pollinating a third of all food we eat. Sacbrood virus (SBV) is a member of the virus family Iflaviridae and affects honey bee larvae and causes particularly devastating disease in the Asian honey bees, Apis cerana. Chinese Sacbrood virus (CSBV) is a geographic strain of SBV identified in China and has resulted in mass death of honey bees in China in recent years. However, the molecular mechanism underlying SBV infection in the Asian honey bee has remained unelucidated. In this present study, we employed high throughput next-generation sequencing technology to study the host transcriptional responses to CSBV infection in A. cerana larvae, and were able to identify genome-wide differentially expressed genes associated with the viral infection. Our study identified 2,534 differentially expressed genes (DEGs) involved in host innate immunity including Toll and immune deficiency (IMD) pathways, RNA interference (RNAi) pathway, endocytosis, etc. Notably, the expression of genes encoding antimicrobial peptides (abaecin, apidaecin, hymenoptaecin, and defensin) and core components of RNAi such as Dicer-like and Ago2 were found to be significantly upregulated in CSBV infected larvae. Most importantly, the expression of Sirtuin target genes, a family of signaling proteins involved in metabolic regulation, apoptosis, and intracellular signaling was found to be changed, providing the first evidence of the involvement of Sirtuin signaling pathway in insects’ immune response to a virus infection. The results obtained from this study provide novel insights into the molecular mechanism and immune responses involved in CSBV infection, which in turn will contribute to the development of diagnostics and treatment for the diseases in honey bees.

Effects of Deformed Wing Virus Infection on Expressions of Immune- and Apoptosis-Related Genes in Western Honeybees (Apis mellifera)

Honeybees are globally threatened by several pathogens, especially deformed wing virus (DWV), as the presence of DWV in western honeybees is indicative of colony loss. The high mortality rate is further exacerbated by the lack of effective treatment, and therefore understanding the immune and apoptosis responses could pave an avenue for the treatment method. In this study, DWV was directly injected into the white-eyed pupae stage of western honeybees (Apis mellifera). The DWV loads and selected gene responses were monitored using the real-time PCR technique. The results showed that honeybee pupae that were injected with the highest concentration of viral loads showed a significantly higher mortality rate than the control groups. Deformed wings could be observed in newly emerged adult bees when the infected bees harbored high levels of viral loads. However, the numbers of viral loads in both normal and crippled wing groups were not significantly different. DWV-injected honeybee pupae with 104 and 107 copy numbers per bee groups showed similar viral loads after 48 h until newly emerged adult bees. Levels of gene expression including immune genes (defensin, abaecin, and hymenoptaecin) and apoptosis genes (buffy, p53, Apaf1, caspase3-like, caspase8-like, and caspase9-like) were analyzed after DWV infection. The expressions of immune and apoptosis genes were significantly different in infected bees compared to those of the control groups. In the pupae stage, the immune genes were activated by injecting DWV (defensin and hymenoptaecin) or Escherichia coli (defensin, abaecin, and hymenoptaecin), a positive control. On the contrary, the expression of apoptosis-related genes (buffy, caspase3-like, caspase8-like, and caspase9-like genes) was suppressed at 96 h post-infection. In DWV-infected newly emerged adult bees, abaecin, hymenoptaecin, Apaf1, and caspase8-like genes were upregulated. However, these genes were not significantly different between the normal and crippled wing bees. Our results suggested that DWV could activate the humoral immunity in honeybees and that honeybee hosts may be able to protect themselves from the virus infection through immune responses. Apoptosis gene expressions were upregulated in newly emerged adult bees by the virus, however, they were downregulated during the initial phase of viral infection.

The Role of Pathogen Dynamics and Immune Gene Expression in the Survival of Feral Honey Bees

Studies of the ecoimmunology of feral organisms can provide valuable insight into how host–pathogen dynamics change as organisms transition from human-managed conditions back into the wild. Honey bees (Apis mellifera Linnaeus) offer an ideal system to investigate these questions as colonies of these social insects often escape management and establish in the wild. While managed honey bee colonies have low probability of survival in the absence of disease treatments, feral colonies commonly survive in the wild, where pathogen pressures are expected to be higher due to the absence of disease treatments. Here, we investigate the role of pathogen infections [Deformed wing virus (DWV), Black queen cell virus (BQCV), and Nosema ceranae] and immune gene expression (defensin-1, hymenoptaecin, pgrp-lc, pgrp-s2, argonaute-2, vago) in the survival of feral and managed honey bee colonies. We surveyed a total of 25 pairs of feral and managed colonies over a 2-year period (2017–2018), recorded overwintering survival, and measured pathogen levels and immune gene expression using quantitative polymerase chain reaction (qPCR). Our results showed that feral colonies had higher levels of DWV but it was variable over time compared to managed colonies. Higher pathogen levels were associated with increased immune gene expression, with feral colonies showing higher expression in five out of the six examined immune genes for at least one sampling period. Further analysis revealed that differential expression of the genes hymenoptaecin and vago increased the odds of overwintering survival in managed and feral colonies. Our results revealed that feral colonies express immune genes at higher levels in response to high pathogen burdens, providing evidence for the role of feralization in altering pathogen landscapes and host immune responses.

Comparing Survival of Israeli Acute Paralysis Virus Infection among Stocks of U.S. Honey Bees

Among numerous viruses that infect honey bees (Apis mellifera), Israeli acute paralysis virus (IAPV) can be linked to severe honey bee health problems. Breeding for virus resistance may improve honey bee health. To evaluate the potential for this approach, we compared the survival of IAPV infection among stocks from the U.S. We complemented the survival analysis with a survey of existing viruses in these stocks and assessing constitutive and induced expression of immune genes. Worker offspring from selected queens in a common apiary were inoculated with IAPV by topical applications after emergence to assess subsequent survival. Differences among stocks were small compared to variation within stocks, indicating the potential for improving honey bee survival of virus infections in all stocks. A positive relation between worker survival and virus load among stocks further suggested that honey bees may be able to adapt to better cope with viruses, while our molecular studies indicate that toll-6 may be related to survival differences among virus-infected worker bees. Together, these findings highlight the importance of viruses in queen breeding operations and provide a promising starting point for the quest to improve honey bee health by selectively breeding stock to be better able to survive virus infections.

Transcriptomic and Epigenomic Dynamics of Honey Bees in Response to Lethal Viral Infection

Honey bees (Apis mellifera L.) suffer from many brood pathogens, including viruses. Despite considerable research, the molecular responses and dynamics of honey bee pupae to viral pathogens remain poorly understood. Israeli Acute Paralysis Virus (IAPV) is emerging as a model virus since its association with severe colony losses. Using worker pupae, we studied the transcriptomic and methylomic consequences of IAPV infection over three distinct time points after inoculation. Contrasts of gene expression and 5 mC DNA methylation profiles between IAPV-infected and control individuals at these time points - corresponding to the pre-replicative (5 h), replicative (20 h), and terminal (48 h) phase of infection - indicate that profound immune responses and distinct manipulation of host molecular processes accompany the lethal progression of this virus. We identify the temporal dynamics of the transcriptomic response to with more genes differentially expressed in the replicative and terminal phases than in the pre-replicative phase. However, the number of differentially methylated regions decreased dramatically from the pre-replicative to the replicative and terminal phase. Several cellular pathways experienced hyper- and hypo-methylation in the pre-replicative phase and later dramatically increased in gene expression at the terminal phase, including the MAPK, Jak-STAT, Hippo, mTOR, TGF-beta signaling pathways, ubiquitin mediated proteolysis, and spliceosome. These affected biological functions suggest that adaptive host responses to combat the virus are mixed with viral manipulations of the host to increase its own reproduction, all of which are involved in anti-viral immune response, cell growth, and proliferation. Comparative genomic analyses with other studies of viral infections of honey bees and fruit flies indicated that similar immune pathways are shared. Our results further suggest that dynamic DNA methylation responds to viral infections quickly, regulating subsequent gene activities. Our study provides new insights of molecular mechanisms involved in epigenetic that can serve as foundation for the long-term goal to develop anti-viral strategies for honey bees, the most important commercial pollinator.

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.

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.

Varroa destructor parasitism has a greater effect on proteome changes than the deformed wing virus and activates TGF-β signaling pathways

Honeybee workers undergo metamorphosis in capped cells for approximately 13 days before adult emergence. During the same period, Varroa mites prick the defenseless host many times. We sought to identify proteome differences between emerging Varroa-parasitized and parasite-free honeybees showing the presence or absence of clinical signs of deformed wing virus (DWV) in the capped cells. A label-free proteomic analysis utilizing nanoLC coupled with an Orbitrap Fusion Tribrid mass spectrometer provided a quantitative comparison of 2316 protein hits. Redundancy analysis (RDA) showed that the combination of Varroa parasitism and DWV clinical signs caused proteome changes that occurred in the same direction as those of Varroa alone and were approximately two-fold higher. Furthermore, proteome changes associated with DWV signs alone were positioned above Varroa in the RDA. Multiple markers indicate that Varroa activates TGF-β-induced pathways to suppress wound healing and the immune response and that the collective action of stressors intensifies these effects. Furthermore, we indicate JAK/STAT hyperactivation, p53-BCL-6 feedback loop disruption, Wnt pathway activation, Wnt/Hippo crosstalk disruption, and NF-κB and JAK/STAT signaling conflict in the Varroa–honeybee–DWV interaction. These results illustrate the higher effect of Varroa than of DWV at the time of emergence. Markers for future research are provided.

Stock-specific chemical brood signals are induced by Varroa and Deformed Wing Virus, and elicit hygienic response in the honey bee

The health of the honey bee Apis mellifera is challenged by the ectoparasitic mite Varroa destructor, and the numerous harmful pathogens it vectors. Existing pesticide-based Varroa controls are not sustainable. In contrast, one promising approach for improved honey bee health is the breeding of hygienic bees, capable of detecting and removing brood that is parasitized or diseased. In three experiments we find evidence to support the hypothesis that stock-specific chemical brood signals are induced by Varroa and Deformed Wing Virus, and elicit hygienic response in the honey bee. By collecting, analyzing, and running bioassays involving mite-infested and control brood extracts from three honey bee breeding stocks we: 1) found evidence that a transferrable chemical signal for hygienic behavior is present in Varroa-infested brood extracts, 2) identified ten stock-specific hydrocarbons as candidates of hygienic signaling, and 3) found that two of these hydrocarbons linked to Varroa and DWV were also elevated in brood targeted for hygienic behavior. These findings expand our understanding of honey bee chemical communication, and facilitate the development of improved hygienic selection tools to breed honey bees with greater resistance to Varroa and associated pathogens.

Knock-on community impacts of a novel vector: spillover of emerging DWV-B from Varroa-infested honeybees to wild bumblebees

Novel transmission routes can directly impact the evolutionary ecology of infectious diseases, with potentially dramatic effect on host populations and knock‐on effects on the wider host community. The invasion of Varroa destructor, an ectoparasitic viral vector in Western honeybees, provides a unique opportunity to examine how a novel vector affects disease epidemiology in a host community. This specialist honeybee mite vectors deformed wing virus (DWV), an important re‐emerging honeybee pathogen that also infects wild bumblebees. Comparing island honeybee and wild bumblebee populations with and without V. destructor, we show that V. destructor drives DWV prevalence and titre in honeybees and sympatric bumblebees. Viral genotypes are shared across hosts, with the potentially more virulent DWV‐B overtaking DWV‐A in prevalence in a current epidemic. This demonstrates disease emergence across a host community driven by the acquisition of a specialist novel transmission route in one host, with dramatic community level knock‐on effects.

Immune related genes as markers for monitoring health status of honey bee colonies

Background

Honey bee population decline threatens the beekeeping sector, agriculture and global biodiversity. Early detection of colony mortality may facilitate rapid interventions to contain and prevent mortality spread. Among others, deformed wing virus (DWV) is capable of inducing colony losses, especially when combined with Varroa destructor mite. Since the bee immune system plays a crucial role in ensuring that bees are able to face these pathogens, we explored whether expression of immune genes could serve as biomarkers of colony health.

Results

Herein, we describe a preliminary immunological marker composed of two immune genes (relish and defensin), which provide insight on honey bee antiviral defense mechanism. Of the tested genes, relish expression correlated with the presence of DWV-Varroa complex, while decreased defensin expression correlated with poor resistance to this complex.

Conclusions

The monitoring of these genes may help us to better understand the complex physiology of honey bees’s immune system and to develop new approaches for managing the health impacts of DWV infection and varroa infestation in the field.

Varroa destructor feeds primarily on honey bee fat body tissue and not hemolymph

The parasitic mite Varroa destructor is the greatest single driver of the global honey bee health decline. Better understanding of the association of this parasite and its host is critical to developing sustainable management practices. Our work shows that this parasite is not consuming hemolymph, as has been the accepted view, but damages host bees by consuming fat body, a tissue roughly analogous to the mammalian liver. Both hemolymph and fat body in honey bees were marked with fluorescent biostains. The fluorescence profile in the guts of mites allowed to feed on these bees was very different from that of the hemolymph of the host bee but consistently matched the fluorescence profile unique to the fat body. Via transmission electron microscopy, we observed externally digested fat body tissue in the wounds of parasitized bees. Mites in their reproductive phase were then fed a diet composed of one or both tissues. Mites fed hemolymph showed fitness metrics no different from the starved control. Mites fed fat body survived longer and produced more eggs than those fed hemolymph, suggesting that fat body is integral to their diet when feeding on brood as well. Collectively, these findings strongly suggest that Varroa are exploiting the fat body as their primary source of sustenance: a tissue integral to proper immune function, pesticide detoxification, overwinter survival, and several other essential processes in healthy bees. These findings underscore a need to revisit our understanding of this parasite and its impacts, both direct and indirect, on honey bee health.

The Dynamics of Deformed Wing Virus Concentration and Host Defensive Gene Expression after Varroa Mite Parasitism in Honey Bees, Apis mellifera

The synergistic interactions between the ectoparasitic mite Varroa destructor and Deformed wing virus (DWV) lead to the reduction in lifespan of the European honey bee Apis mellifera and often have been implicated in colony losses worldwide. However, to date, the underlying processes and mechanisms that form the multipartite interaction between the bee, mite, and virus have not been fully explained. To gain a better understanding of honey bees’ defense response to Varroa mite infestation and DWV infection, the DWV titers and transcription profiles of genes originating from RNAi, immunity, wound response, and homeostatic signaling pathways were monitored over a period of eight days. With respect to DWV, we observed low viral titers at early timepoints that coincided with high levels of Toll pathway transcription factor Dorsal, and its downstream immune effector molecules Hymenoptaecin, Apidaecin, Abaecin, and Defensin 1. However, we observed a striking increase in viral titers beginning after two days that coincided with a decrease in Dorsal levels and its corresponding immune effector molecules, and the small ubiquitin-like modifier (SUMO) ligase repressor of Dorsal, PIAS3. We observed a similar expression pattern for genes expressing transcripts for the RNA interference (Dicer/Argonaute), wound/homeostatic (Janus Kinase), and tissue growth (Map kinase/Wnt) pathways. Our results demonstrate that on a whole, honey bees are able to mount an immediate, albeit, temporally limited, immune and homeostatic response to Varroa and DWV infections, after which downregulation of these pathways leaves the bee vulnerable to expansive viral replication. The critical insights into the defense response upon Varroa and DWV challenges generated in this study may serve as a solid base for future research on the development of effective and efficient disease management strategies in honey bees.

Israeli Acute Paralysis Virus: Honey Bee Queen⁻Worker Interaction and Potential Virus Transmission Pathways

Queen loss or failure is an important cause of honey bee colony loss. A functional queen is essential to a colony, and the queen is predicted to be well protected by worker bees and other mechanisms of social immunity. Nevertheless, several honey bee pathogens (including viruses) can infect queens. Here, we report a series of experiments to test how virus infection influences queen⁻worker interactions and the consequences for virus transmission. We used Israeli acute paralysis virus (IAPV) as an experimental pathogen because it is relevant to bee health but is not omnipresent. Queens were observed spending 50% of their time with healthy workers, 32% with infected workers, and 18% without interaction. However, the overall bias toward healthy workers was not statistically significant, and there was considerable individual to individual variability. We found that physical contact between infected workers and queens leads to high queen infection in some cases, suggesting that IAPV infections also spread through close bodily contact. Across experiments, queens exhibited lower IAPV titers than surrounding workers. Thus, our results indicate that honey bee queens are better protected by individual and social immunity, but this protection is insufficient to prevent IAPV infections completely.

Bee Viruses: Ecology, Pathogenicity, and Impacts

Bees-including solitary, social, wild, and managed species-are key pollinators of flowering plant species, including nearly three-quarters of global food crops. Their ecological importance, coupled with increased annual losses of managed honey bees and declines in populations of key wild species, has focused attention on the factors that adversely affect bee health, including viral pathogens. Genomic approaches have dramatically expanded understanding of the diversity of viruses that infect bees, the complexity of their transmission routes-including intergenus transmission-and the diversity of strategies bees have evolved to combat virus infections, with RNA-mediated responses playing a prominent role. Moreover, the impacts of viruses on their hosts are exacerbated by the other major stressors bee populations face, including parasites, poor nutrition, and exposure to chemicals. Unraveling the complex relationships between viruses and their bee hosts will lead to improved understanding of viral ecology and management strategies that support better bee health.

Brood Affects Hygienic Behavior in the Honey Bee (Hymenoptera: Apidae)

Despite receiving much attention, the ectoparasitic mite Varroa destructor (Anderson and Trueman) and the pathogens it vectors remain critical threats to the health of the honey bee Apis mellifera (Linnaeus) (Hymenoptera: Apidae). One promising intervention approach is the breeding of hygienic honey bees, which have an improved ability to detect and remove unhealthy brood from the colony, and are thus more resistant to Varroa. While much hygienic behavior-related research has focused on enhanced adult honey bee olfaction, less attention has been paid to the olfactory signals that originate inside the brood cell, triggering hygienic removal. Here, we hypothesized that selection for hygienic behavior in honey bees has influenced brood signaling, predicting that: 1) in a common social environment, removal rates differ among brood with different selective breeding histories, and 2) the removal rates of brood positively correlate to the hygiene level of the brood's colony of origin. To test these predictions, we cross-fostered brood subjected to control, wound, or Varroa treatment in unselected (UNS), Minnesota Hygienic (HYG), and Varroa-Sensitive Hygienic (VSH) colonies, and monitored individual brood cells for hygienic removal. Results confirmed both predictions, as brood from hygienic colonies was more likely to be removed than brood from UNS colonies, regardless of where the brood was fostered. These findings suggest that hygiene-related brood signals complement previously identified characteristics of hygienic adults, constituting an important mechanism of social immunity in honey bees. Thus, selective breeding for honey bee hygienic behavior may be improved through the utilization of field assays containing compounds related to larval signaling.

Pathogenicity of Serratia marcescens Strains in Honey Bees

Recently, it has become apparent that multiple factors are responsible for honey bee decline, including climate change, pests and pathogens, pesticides, and loss of foraging habitat. Of the large number of pathogens known to infect honey bees, very few are bacteria. Because adult workers abandon hives when diseased, many of their pathogens may go unnoticed. Here we characterized the virulence of Serratia marcescens strains isolated from honey bee guts and hemolymph. Our results indicate that S. marcescens, an opportunistic pathogen of many plants and animals, including humans, is a virulent opportunistic pathogen of honey bees, which could contribute to bee decline. Aside from the implications for honey bee health, the discovery of pathogenic S. marcescens strains in honey bees presents an opportunity to better understand how opportunistic pathogens infect and invade hosts.

Although few honey bee diseases are known to be caused by bacteria, pathogens of adult worker bees may be underrecognized due to social immunity mechanisms. Specifically, infected adult bees typically abandon the hive or are removed by guards. Serratia marcescens, an opportunistic pathogen of many plants and animals, is often present at low abundance in the guts of honey bee workers and has recently been isolated from Varroa mites and from the hemolymph of dead and dying honey bees. However, the severity and prevalence of S. marcescens pathogenicity in honey bees have not been fully investigated. Here we characterized three S. marcescens strains isolated from the guts of honey bees and one previously isolated from hemolymph. In vivo tests confirmed that S. marcescens is pathogenic in workers. All strains caused mortality when a few cells were injected into the hemocoel, and the gut-isolated strains caused mortality when administered orally. In vitro assays and comparative genomics identified possible mechanisms of virulence of gut-associated strains. Expression of antimicrobial peptide and phenoloxidase genes was not elevated following infection, suggesting that these S. marcescens strains derived from honey bees can evade the immune response in their hosts. Finally, surveys from four locations in the United States indicated the presence of S. marcescens in the guts of over 60% of the worker bees evaluated. Taken together, these results suggest that S. marcescens is a widespread opportunistic pathogen of adult honey bees and that it may be highly virulent under some conditions such as perturbation of the normal gut microbiota or the presence of Varroa mites that puncture the integument, thereby enabling entry of bacterial cells.

Honey Bee and Bumble Bee Antiviral Defense

Bees are important plant pollinators in both natural and agricultural ecosystems. Managed and wild bees have experienced high average annual colony losses, population declines, and local extinctions in many geographic regions. Multiple factors, including virus infections, impact bee health and longevity. The majority of bee-infecting viruses are positive-sense single-stranded RNA viruses. Bee-infecting viruses often cause asymptomatic infections but may also cause paralysis, deformity or death. The severity of infection is governed by bee host immune responses and influenced by additional biotic and abiotic factors. Herein, we highlight studies that have contributed to the current understanding of antiviral defense in bees, including the Western honey bee (Apis mellifera), the Eastern honey bee (Apis cerana) and bumble bee species (Bombus spp.). Bee antiviral defense mechanisms include RNA interference (RNAi), endocytosis, melanization, encapsulation, autophagy and conserved immune pathways including Jak/STAT (Janus kinase/signal transducer and activator of transcription), JNK (c-Jun N-terminal kinase), MAPK (mitogen-activated protein kinases) and the NF-κB mediated Toll and Imd (immune deficiency) pathways. Studies in Dipteran insects, including the model organism Drosophila melanogaster and pathogen-transmitting mosquitos, provide the framework for understanding bee antiviral defense. However, there are notable differences such as the more prominent role of a non-sequence specific, dsRNA-triggered, virus limiting response in honey bees and bumble bees. This virus-limiting response in bees is akin to pathways in a range of organisms including other invertebrates (i.e., oysters, shrimp and sand flies), as well as the mammalian interferon response. Current and future research aimed at elucidating bee antiviral defense mechanisms may lead to development of strategies that mitigate bee losses, while expanding our understanding of insect antiviral defense and the potential evolutionary relationship between sociality and immune function.

Transcriptome profiling of the honeybee parasite Varroa destructor provides new biological insights into the mite adult life cycle

Background

The parasite Varroa destructor represents a significant threat to honeybee colonies. Indeed, development of Varroa infestation within colonies, if left untreated, often leads to the death of the colony. Although its impact on bees has been extensively studied, less is known about its biology and the functional processes governing its adult life cycle and adaptation to its host. We therefore developed a full life cycle transcriptomic catalogue in adult Varroa females and included pairwise comparisons with males, artificially-reared and non-reproducing females (10 life cycle stages and conditions in total).

Results

Extensive remodeling of the Varroa transcriptome was observed, with an upregulation of energetic and chitin metabolic processes during the initial and final phases of the life cycle (e.g. phoretic and post-oviposition stages), whereas during reproductive stages in brood cells genes showing functions related to transcriptional regulation were overexpressed. Several neurotransmitter and neuropeptide receptors involved in behavioural regulation, as well as active compounds of salivary glands, were also expressed at a higher level outside the reproductive stages. No difference was detected between artificially-reared phoretic females and their counterparts in colonies, or between females who failed to reproduce and females who successfully reproduced, indicating that phoretic individuals can be reared outside host colonies without impacting their physiology and that mechanisms underlying reproductive failure occur before oogenesis.

Conclusions

We discuss how these new findings reveal the remarkable adaptation of Varroa to its host biology and notably to the switch from living on adults to reproducing in sealed brood cells. By spanning the entire adult life cycle, our work captures the dynamic changes in the parasite gene expression and serves as a unique resource for deciphering Varroa biology and identifying new targets for mite control.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4668-z) contains supplementary material, which is available to authorized users.

Quantitative proteomics reveals divergent responses in Apis mellifera worker and drone pupae to parasitization by Varroa destructor

Varroa destructor is a haemophagous ectoparasite of honeybees and is considered a major causal agent of colony losses in Europe and North America. Although originating in Eastern Asia where it parasitizes Apis cerana, it has shifted hosts to the western honeybee Apis mellifera on which it has a greater deleterious effect on the individual and colony level. To investigate this important host-parasite interaction and to determine whether Varroa causes different effects on different castes we conducted a label free quantitative proteomic analysis of Varroa-parasitized and non-parasitized drone and worker Apis mellifera pupae. 1195 proteins were identified in total, of which 202 and 250 were differentially abundant in parasitized drone and worker pupae, respectively. Both parasitized drone and worker pupae displayed reduced abundance in proteins associated with the cuticle, lipid transport and innate immunity. Proteins involved in metabolic processes were more abundant in both parasitized castes although the response in workers was more pronounced. A number of caste specific responses were observed including differential abundance of numerous cytoskeletal and muscle proteins, which were of higher abundance in parasitized drones in comparison to parasitized workers. Proteins involved in fatty acid and carbohydrate metabolism were more abundant in parasitized workers as were a large number of ribosomal proteins highlighting either potentially divergent responses to Varroa or a different strategy by the mite when parasitizing the different castes. This data improves our understanding of this interaction and may provide a basis for future studies into improvements to therapy and control of Varroasis.

Varroa destructor induces changes in the expression of immunity-related genes during the development of Apis mellifera worker and drone broods

The ectoparasitic mite Varroa destructor has emerged as the major pest of honeybees. Despite extensive research efforts, the pathogenesis of varroosis has not been fully explained. Earlier studies suggested that V. destructor infestation leads to the suppression of the host's immune system. The aim of this study was to analyze the immune responses of 14 genes in the Toll signal transduction pathways, including effector genes of antimicrobial peptides (AMPs), in developing Apis mellifera workers and drones infested with V. destructor. Four developmental stages (L5 larvae, prepupae, and 2 pupal stages) and newly emerged imagines were analyzed. In workers, the most significant changes were observed in L5 larvae in the initial stages of infestation. A significant increase in the relative expression of 10 of the 14 analyzed genes, including defensin-1 and defensin-2, was observed in infested bees relative to non-infested individuals. The immune response in drones developed at a slower rate. The expression of genes regulating cytoplasmic signal transduction increased in prepupae, whereas the expression of defensin-1 and defensin-2 effector genes increased in P3 pupae with red eyes. The expression of many immunity-related genes was silenced in successive life stages and in imagines, and it was more profound in workers than in drones. The results indicate that V. destructor significantly influences immune responses regulated by the Toll signal transduction pathway in bees. In infested bees, the observed changes in Toll pathway genes varied between life stages and the sexes.

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.

Transcriptional signatures of parasitization and markers of colony decline in Varroa-infested honey bees (Apis mellifera)

Extensive annual losses of honey bee colonies (Apis mellifera L.) reported in the northern hemisphere represent a global problem for agriculture and biodiversity. The parasitic mite Varroa destructor, in association with deformed wing virus (DWV), plays a key role in this phenomenon, but the underlying mechanisms are still unclear. To elucidate these mechanisms, we analyzed the gene expression profile of uninfested and mite infested bees, under laboratory and field conditions, highlighting the effects of parasitization on the bee's transcriptome under a variety of conditions and scenarios. Parasitization was significantly correlated with higher viral loads. Honey bees exposed to mite infestation exhibited an altered expression of genes related to stress response, immunity, nervous system function, metabolism and behavioural maturation. Additionally, mite infested young bees showed a gene expression profile resembling that of forager bees. To identify potential molecular markers of colony decline, the expression of genes that were commonly regulated across the experiments were subsequently assessed in colonies experiencing increasing mite infestation levels. These studies suggest that PGRP-2, hymenoptaecin, a glucan recognition protein, UNC93 and a p450 cytocrome maybe suitable general biomarkers of Varroa-induced colony decline. Furthermore, the reliability of vitellogenin, a yolk protein previously identified as a good marker of colony survival, was confirmed here.

Virus and dsRNA-triggered transcriptional responses reveal key components of honey bee antiviral defense

Recent high annual losses of honey bee colonies are associated with many factors, including RNA virus infections. Honey bee antiviral responses include RNA interference and immune pathway activation, but their relative roles in antiviral defense are not well understood. To better characterize the mechanism(s) of honey bee antiviral defense, bees were infected with a model virus in the presence or absence of dsRNA, a virus associated molecular pattern. Regardless of sequence specificity, dsRNA reduced virus abundance. We utilized next generation sequencing to examine transcriptional responses triggered by virus and dsRNA at three time-points post-infection. Hundreds of genes exhibited differential expression in response to co-treatment of dsRNA and virus. Virus-infected bees had greater expression of genes involved in RNAi, Toll, Imd, and JAK-STAT pathways, but the majority of differentially expressed genes are not well characterized. To confirm the virus limiting role of two genes, including the well-characterized gene, dicer, and a probable uncharacterized cyclin dependent kinase in honey bees, we utilized RNAi to reduce their expression in vivo and determined that virus abundance increased, supporting their involvement in antiviral defense. Together, these results further our understanding of honey bee antiviral defense, particularly the role of a non-sequence specific dsRNA-mediated antiviral pathway.

Role of Human Action in the Spread of Honey Bee (Hymenoptera: Apidae) Pathogens

The increased annual losses in European honey bee (Apis mellifera) colonies in North America and some other countries is usually attributed to a range of factors including pathogens, poor nutrition, and insecticides. In this essay, I will argue that the global trade in honey bees and migratory beekeeping practices within countries has enabled pathogens to spread quickly. Beekeepers' management strategies have also contributed to the spread of pathogens as well as the development of resistance to miticides and antibiotics, and exacerbated by hobby beekeepers. The opportunities for arresting honey bee declines rest as strongly with individual beekeepers as they do with the dynamics of disease.

Effects of Imidacloprid and Varroa destructor on survival and health of European honey bees, Apis mellifera

There has been growing concern over declines in populations of honey bees and other pollinators which are a vital part to our food security. It is imperative to identify factors responsible for accelerated declines in bee populations and develop solutions for reversing bee losses. While exact causes of colony losses remain elusive, risk factors thought to play key roles are ectoparasitic mites Varroa destructor and neonicotinoid pesticides. The present study aims to investigate effects of a neonicotinoid pesticide Imidacloprid and Varroa mites individually on survivorship, growth, physiology, virus dynamics and immunity of honey bee workers. Our study provides clear evidence that the exposure to sublethal doses of Imidacloprid could exert a significantly negative effect on health and survival of honey bees. We observed a significant reduction in the titer of vitellogenin (Vg), an egg yolk precursor that regulates the honey bees development and behavior and often are linked to energy homeostasis, in bees exposed to Imidacloprid. This result indicates that sublethal exposure to neonicotinoid could lead to increased energy usage in honey bees as detoxification is a energy-consuming metabolic process and suggests that Vg could be a useful biomarker for measuring levels of energy stress and sublethal effects of pesticides on honey bees. Measurement of the quantitative effects of different levels of Varroa mite infestation on the replication dynamic of Deformed wing virus (DWV), an RNA virus associated with Varroa infestation, and expression level of immune genes yields unique insights into how honey bees respond to stressors under laboratory conditions.

Oldest Varroa tolerant honey bee population provides insight into the origins of the global decline of honey bees

The ecto-parasitic mite Varroa destructor has transformed the previously inconsequential Deformed Wing Virus (DWV) into the most important honey bee viral pathogen responsible for the death of millions of colonies worldwide. Naturally, DWV persists as a low level covert infection transmitted between nest-mates. It has long been speculated that Varroa via immunosuppression of the bees, activate a covert infection into an overt one. Here we show that despite Varroa feeding on a population of 20-40 colonies for over 30 years on the remote island of Fernando de Noronha, Brazil no such activation has occurred and DWV loads have remained at borderline levels of detection. This supports the alternative theory that for a new vector borne viral transmission cycle to start, an outbreak of an overt infection must first occur within the host. Therefore, we predict that this honey bee population is a ticking time-bomb, protected by its isolated position and small population size. This unique association between mite and bee persists due to the evolution of low Varroa reproduction rates. So the population is not adapted to tolerate Varroa and DWV, rather the viral quasispecies has simply not yet evolved the necessary mutations to produce a virulent variant.

Unity in defence: honeybee workers exhibit conserved molecular responses to diverse pathogens

BACKGROUND

Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses.

RESULTS

We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses.

CONCLUSIONS

Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Early life stress affects mortality rate more than social behavior, gene expression or oxidative damage in honey bee workers

Early life stressors can affect aging and life expectancy in positive or negative ways. Individuals can adjust their behavior and molecular physiology based on early life experiences but relatively few studies have connected such mechanisms to demographic patterns in social organisms. Sociality buffers individuals from environmental influences and it is unclear how much early life stress affects later life history. Workers of the honey bee (Apis mellifera L.) were exposed to two stressors, Varroa parasitism and Paraquat exposure, early in life. Consequences were measured at the molecular, behavioral, and demographic level. While treatments did not significantly affect levels of oxidative damage, expression of select genes, and titers of the common deformed wing virus, most of these measures were affected by age. Some of the age effects, such as declining levels of deformed wing virus and oxidative damage, were opposite to our predictions but may be explained by demographic selection. Further analyses suggested some influences of worker behavior on mortality and indicated weak treatment effects on behavior. The latter effects were inconsistent among the two experiments. However, mortality rate was consistently reduced by Varroa mite stress during development. Thus, mortality was more responsive to early life stress than our other response variables. The lack of treatment effects on these measures may be due to the social organization of honey bees that buffers the individual from the impact of stressful developmental conditions.

Effect of Varroa destructor, Wounding and Varroa Homogenate on Gene Expression in Brood and Adult Honey Bees

Honey bee (Apis mellifera) gene expression related to immunity for hymenoptaecin (AmHym) and defensin-1 (AmDef-1), longevity for vitellogenin (AmVit2) and stem cell proliferation for poly U binding factor 68 kDa (AmPuf68) was compared following Varroa destructor parasitism, buffer injection and injection of V. destructor compounds in its homogenate. In adults, V. destructor parasitism decreased expression of all four genes, while buffer injection decreased expression of AmHym, AmPuf68 and AmVit2, and homogenate injection decreased expression of AmPuf68 and AmVit2 but increased expression of AmDef-1 relative to their respective controls. The effect of V. destructor parasitism in adults relative to the controls was not significantly different from buffer injection for AmHym and AmVit2 expression, and it was not significantly different from homogenate injection for AmPuf68 and AmVit2. In brood, V. destructor parasitism, buffer injection and homogenate injection decreased AmVit2 expression, whereas AmHym expression was decreased by V. destructor parasitism but increased by buffer and homogenate injection relative to the controls. The effect of varroa parasitism in brood was not significantly different from buffer or homogenate injection for AmPuf68 and AmVit2. Expression levels of the four genes did not correlate with detectable viral levels in either brood or adults. The results of this study indicate that the relative effects of V. destructor parasitism on honey bee gene expression are also shared with other types of stresses. Therefore, some of the effects of V. destructor on honey bees may be mostly due to wounding and injection of foreign compounds into the hemolymph of the bee during parasitism. Although both brood and adults are naturally parasitized by V. destructor, their gene expression responded differently, probably the result of different mechanisms of host responses during development.

Synthesis of Enantiopure Alicyclic Ethers and Their Activity on the Chemosensory Organ of the Ectoparasite of Honey Bees, Varroa destructor

The preparation of enantiopure conformationally restricted alicyclic ethers and their inhibitory activities on the chemosensory organ of the Varroa destructor, a parasite of honey bees, are reported in this article. We tested the effect of enantiopure ethers of cis-5-(2'-hydroxyethyl)cyclopent-2-en-1-ol on the Varroa chemosensory organ by electrophysiology, for their ability to inhibit the responses to two honey bee-produced odors that are important for the mite to locate its host: nurse bee head space odor and (E)-β-ocimene, a honey bee brood pheromone. Previous work with the racemic compounds showed that they suppress the mite's olfactory response to its bee host, which led to incorrect host choice. Based on a structure-activity relationship, we predicted that the two most active compounds-cis-1-butoxy-5-(2'-methoxyethyl)cyclopent-2-ene, cy{4,1}, and (cis-1-ethoxy-5-(2'ethoxyethyl)cyclopent-2-ene, cy{2,2}-could have opposite active enantiomers. Here we studied the enantiomers of both ethers, whose preparation involved enzymatic resolution of racemic diol cis-5-(2'-hydroxyethyl)cyclopent-2-en-1-ol using Lipase AK with vinyl acetate. The racemic diol was prepared from commercially available 2,5-norbornadiene. We observed that the responses of the chemosensory organ to honey bee head space volatiles were significantly decreased by both enantiomers of cy{4,1} and cy{2,2}, but that responses to (E)-β-ocimene were decreased significantly only by (+)-cy{4,1} (1R,5S) and (-)-cy{2,2} (1S,5R) and not by their respective enantiomers. The importance of this result is that the racemates could be used to inhibit olfactory detection of bee odors by mites, without a loss in activity relative to the more expensive enantiopure compounds.

Pathogenesis of varroosis at the level of the honey bee (Apis mellifera) colony

The parasitic mite Varroa destructor, in interaction with different viruses, is the main cause of honey bee colony mortality in most parts of the world. Here we studied how effects of individual-level parasitization are reflected by the bee colony as a whole. We measured disease progression in an apiary of 24 hives with differing degree of mite infestation, and investigated its relationship to 28 biometrical, physiological and biochemical indicators. In early summer, when the most heavily infested colonies already showed reduced growth, an elevated ratio of brood to bees, as well as a strong presence of phenoloxidase/prophenoloxidase in hive bees were found to be predictors of the time of colony collapse. One month later, the learning performance of worker bees as well as the activity of glucose oxidase measured from head extracts were significantly linked to the timing of colony collapse. Colonies at the brink of collapse were characterized by reduced weight of winter bees and a strong increase in their relative body water content. Our data confirm the importance of the immune system, known from studies of individually-infested bees, for the pathogenesis of varroosis at colony level. However, they also show that single-bee effects cannot always be extrapolated to the colony as a whole. This fact, together with the prominent role of colony-level factors like the ratio between brood and bees for disease progression, stress the importance of the superorganismal dimension of Varroa research.