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

Transcriptional control of a metabolic switch regulating cellular methylation reactions is part of a common response to stress in divergent bee species

Recent global declines in bee health have elevated the need for a more complete understanding of the cellular stress mechanisms employed by diverse bee species. We recently uncovered the biomarker lethal (2) essential for life [l(2)efl] genes as part of a shared transcriptional program in response to a number of cell stressors in the western honey bee (Apis mellifera). Here, we describe another shared stress-responsive gene, glycine N-methyltransferase (Gnmt), which is known as a key metabolic switch controlling cellular methylation reactions. We observed Gnmt induction by both abiotic and biotic stressors. We also found increased levels of the GNMT reaction product sarcosine in the midgut after stress, linking metabolic changes with the observed changes in gene regulation. Prior to this study, Gnmt upregulation had not been associated with cellular stress responses in other organisms. To determine whether this novel stress-responsive gene would behave similarly in other bee species, we first characterized the cellular response to endoplasmic reticulum (ER) stress in lab-reared adults of the solitary alfalfa leafcutting bee (Megachile rotundata) and compared this with age-matched honey bees. The novel stress gene Gnmt was induced in addition to a number of canonical gene targets induced in both bee species upon unfolded protein response (UPR) activation, suggesting that stress-induced regulation of cellular methylation reactions is a common feature of bees. Therefore, this study suggests that the honey bee can serve as an important model for bee biology more broadly, although studies on diverse bee species will be required to fully understand global declines in bee populations.

Ontogeny of immunity and natural viral infection in Apis mellifera drones and workers

The most common viral diseases affecting honey bees (Apis mellifera) in Israel include deformed wing viruses (DWV-A and DWV-B) and acute paralysis viruses (ABPV and IAPV). These viruses are transmitted within and between colonies, both horizontally and vertically. All members of the colony contribute to this transmission, on the other hand individual and social immunity, particularly hygienic behaviour, may affect the outcome of the process. In this study, we evaluated the ontogeny of natural infections of DWV-A, DWV-B, ABPV and IAPV, their prevalence and loads, in workers and drones from high (H) and low (L) hygienic colonies. In parallel, we evaluated the expression of two immune genes: peptidoglycan recognition protein S2(PGRP-S2) and hymenoptaecin. The prevalence of DWV-B and IAPV increased with age and was higher in workers than in drones. ABPV was not detected in drones. The expression of both immune genes was significantly affected by age and sex. Drones from H colonies had higher expression of these genes. The increased expression of immune genes with drones' age, particularly in hygienic colonies, suggest additional value of honey bee breeding for hygienic behaviour for sustainable beekeeping.

Alterations in Histone Methylation States Increased Profusion of Lethal(2)-Essential-for-Life-Like (l(2)elf), Trithorax and Polycomb Genes in Apis mellifera under Heat Stress

Histone post-translational modifications (PTMs) represent a key mechanism in the thermal adaptation of the honeybee Apis mellifera. In this study, a chromatin immunoprecipitation assay and qPCR were employed to explore the changes in the methylation states of H3K4m2, H3K4m3, H3K27m2 and H3K27m3 associated with l2efl (ID: 72474, 724405, 724488), histone methyltransferases (HMTs) ((trx) and PR-set7) and Polycomb (Pc) and (Su(z)12) genes in A. m. jemenitica (tolerant subspecies) and A. m. carnica (susceptible subspecies) in response to heat treatment (42 °C for 1 h). The results revealed significant enrichment fold changes in the methylation/demethylation of most H3K4 and H3K27 marks at all targeted genes. These changes increased the profusion of l2efl (ID: 72474, 724405, 724488), histone methyltransferases (HMTs) (trx) and Polycomb (Pc) and Su(z)12 and decreased the profusion of HMT (PR-set7) in both honeybee subspecies. The changes in the methylation enrichment folds of histone methyltransferases (HMTs) ((trx), PR-set) and Polycomb (Pc), Su(z)12 genes demonstrate the well-harmonized coordination of epigenetic gene regulation in response to heat treatment. Compared to the control, the changes in the methylation enrichment folds of H3K4m3 at Polycomb Su(z)12 were about 30× and 100× higher in treated A. m. jemenitica and A.m. carnica, respectively. Similarly, changes in the methylation/demethylation enrichment folds of HMT (trx) and Polycomb (Pc) and Su(z)12 were 2-3× higher in A. m. carnica than in A. m. jemenitica after treatment (42 °C). It is evident that post-translational chromatin modification in both honeybee subspecies can diminish heat stress impact by (I) increasing the transcriptional provision of l2efl associated with survival and (II) increasing the silencing of genes associated with general cellular activities.

Alcohol extract of the gypsy mushroom (Cortinarius caperatus) inhibits the development of Deformed wing virus infection in western honey bee (Apis mellifera)

Deformed wing virus (DWV) transmitted by the parasitic mite Varroa destructor is one of the most significant factors contributing to massive losses of managed colonies of western honey bee (Apis mellifera) subspecies of European origin reported worldwide in recent decades. Despite this fact, no antiviral treatment against honey bee viruses is currently available for practical applications and the level of viral infection can only be controlled indirectly by reducing the number of Varroa mites in honey bee colonies. In this study, we investigated the antiviral potential of the gypsy mushroom (Cortinarius caperatus) to reduce DWV infection in honey bees. Our results indicate that the alcohol extract of C. caperatus prevented the development of DWV infection in cage experiments as well as after direct application to honey bee colonies in a field experiment. The applied doses did not shorten the lifespan of honey bees. The reduced levels of DWV in C. caperatus-treated honey bees in cage experiments were accompanied by significant changes in the gene expression of Tep7, Bap1, and Vago. The C. caperatus treatment was not effective against the trypanosomatid Lotmaria passim. No residues of C.caperatus were found in honey harvested in the spring from colonies supplemented with the mushroom extract for their winter feeding. These findings suggest that C. caperatus alcohol extract could be a potential natural remedy to treat DWV infection in honey bees.

Age-related response to mite parasitization and viral infection in the honey bee suggests a trade-off between growth and immunity

Host age at parasites' exposure is often neglected in studies on host-parasite interactions despite the important implications for epidemiology. Here we compared the impact of the parasitic mite Varroa destructor, and the associated pathogenic virus DWV on different life stages of their host, the western honey bee Apis mellifera. The pre-imaginal stages of the honey bee proved to be more susceptible to mite parasitization and viral infection than adults. The higher viral load in mite-infested bees and DWV genotype do not appear to be the drivers of the observed difference which, instead, seems to be related to the immune-competence of the host. These results support the existence of a trade-off between immunity and growth, making the pupa, which is involved in the highly energy-demanding process of metamorphosis, more susceptible to parasites and pathogens. This may have important implications for the evolution of the parasite's virulence and in turn for honey bee health. Our results highlight the important role of host's age and life stage at exposure in epidemiological modelling. Furthermore, our study could unravel new aspects of the complex honey bee-Varroa relationship to be addressed for a sustainable management of this parasite.

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.

A conserved RNAi molecule Ago2 involved in antiviral immunity of oyster Crassostrea gigas

Argonaute (Ago) is the core component of RNA-induced silencing complex to play a crucial role in the antiviral immunity, which always cooperates with Dicer in RNA interference (RNAi) to silence the target genes. In the present study, an Ago homologue (CgAgo2) was identified in the Pacific oyster Crassostrea gigas. There were four classical functional domains in the predicted CgAgo2 protein, including an N-terminal domain, a PAZ domain, a Mid domain, and a PIWI domain. The deduced amino acid sequence of CgAgo2 shared 63.52%-84.27% identity with other Agos. Transcriptome analysis showed that CgAgo2 was highly expressed in embryonic period and gradually decreased from blastula to gastrula. The transcripts of CgAgo2 were detectable in all the examined tissues of adult oysters, with the highest expression in haemocytes (36.61-fold of that in adductor muscle, p < 0.001). The expression level of CgAgo2 mRNA in haemocytes increased significantly at 12 h after poly (I:C) and dsRNA stimulation, which were 2.71-fold (p < 0.05) and 58.00-fold (p < 0.001) of that in the control group respectively. Immunocytochemistry assay revealed that CgAgo2 proteins were mainly distributed in the cytoplasm and nucleus of haemocytes. The interaction between the recombinant CgAgo2 protein (rCgAgo2) and cleavage protein rCgDicer was observed in vitro by BLI and pull-down assays. These results indicated that CgAgo2 participated in the antiviral immunity of oyster by functioning as a component of RNA-induced silencing complex in RNAi.

Large-Scale Application of Double-Stranded RNA Shows Potential for Reduction of Sacbrood Virus Disease in Apis cerana Apiaries

Sacbrood virus (SBV) infection has emerged as a remarkable threat to Apis cerana colonies in South Korea, necessitating prompt control measures. In this study, RNA interference (RNAi) targeting the VP3 gene was developed to assess its safety and efficacy in protecting and treating SBV in vitro and in infected colonies in South Korean apiaries. The efficacy of VP3 double-stranded RNA (dsRNA) was demonstrated in laboratory-based experiments, wherein infected larvae treated with VP3 dsRNA exhibited a 32.7% increase in survival rate compared to untreated larvae. Data from a large-scale field trial indicate the efficacy of dsRNA treatment since none of the treated colonies had symptomatic SBV infections, whereas disease was observed in 43% (3/7) of the control colonies. In the 102 colonies exhibiting symptoms of SBV disease, RNAi treatment provided partial protection with weekly treatment, prolonging the survival period of colonies to 8 months compared to 2 months in colonies treated at 2- and 4-week intervals. Therefore, this study demonstrated that RNAi is a valuable tool for preventing SBV disease outbreaks in healthy and low-level SBV-infected colonies.

Can immune gene silencing via dsRNA feeding promote pathogenic viruses to control the globally invasive Argentine ant?

Pest control methods that can target pest species with limited environmental impacts are a conservation and economic priority. Species-specific pest control using RNA interference is a challenging but promising avenue in developing the next generation of pest management. We investigate the feasibility of manipulating a biological invader's immune system using double-stranded RNA (dsRNA) in order to increase susceptibility to naturally occurring pathogens. We used the invasive Argentine ant as a model, targeting the immunity-associated genes Spaetzle and Dicer-1 with dsRNA. We show that feeding with Spaetzle dsRNA can result in partial target gene silencing for up to 28 days in the laboratory and 5 days in the field. Dicer-1 dsRNA only resulted in partial gene knockdown after 2 days in the laboratory. Double-stranded RNA treatments were associated with significant gene expression disruptions across immune pathways in the laboratory and to a lower extent in the field. In total, 12 viruses and four bacteria were found in these ant populations. Some changes in viral loads in dsRNA-treated groups were observed. For example, Linepithema humile Polycipivirus 2 (LhuPCV2) loads increased after 2 days of treatment with Spaetzle and Dicer-1 dsRNA treatments in the laboratory. After treatment with the dsRNA in the field, after 5 days the virus Linepithema humile toti-like virus 1 (LhuTLV1) was significantly more abundant. However, immune pathway disruption did not result in a consistent increase in microbial infections, nor did it alter ant abundance in the field. Some viruses even declined in abundance after dsRNA treatment. Our study explored the feasibility of lowering a pest's immunity as a control tool. We demonstrate that it is possible to alter immune gene expression of pest species and pathogen loads, although in our specific system the affected pathogens did not appear to influence pest abundance. We provide suggestions on future directions for dsRNA-mediated immune disruption in pest species, including potential avenues to improve dsRNA delivery as well as the importance of pest and pathogen biology. Double-stranded RNA targeting immune function might be especially useful for pest control in systems in which viruses or other microorganisms are prevalent and have the potential to be pathogenic.

Differential Gene Expression of Innate Immune Response Genes Consequent to Solenopsis invicta Virus-3 Infection

The red imported fire ant Solenopsis invicta Buren (fire ant hereafter) is a global pest that inflicts billions of dollars in damages to the United States economy and poses a major threat on a global scale. Concerns with the broad-spectrum application of insecticides have facilitated the hunt for natural enemy-mediated controls. One of these, the virus Solenopsis invicta virus-3 (SINV-3 hereafter) is exceptionally virulent in laboratory settings. However, despite high mortality rates in the laboratory and documented widespread SINV-3 prevalence in the southern United States, the fire ant remains a major pest. To explore this paradox, we document the immune response elicited by the fire ant when infected with SINV-3. We sequence the fire ant transcriptome prior to and following infection with SINV-3, and identify and discuss in detail genes in immune response pathways differentially expressed following infection with SINV-3. This information provides insights into genes and pathways involved in the SINV-3 infection response in the fire ant and offers avenues to pursue, to suppress key immune response genes and force the fire ant to succumb to SINV-3 infection in the field.

Addressing the challenges of symbiont-mediated RNAi in aphids

Because aphids are global agricultural pests and models for bacterial endosymbiosis, there is a need for reliable methods to study and control their gene function. However, current methods available for aphid gene knockout and knockdown of gene expression are often unreliable and time consuming. Techniques like CRISPR-Cas genome editing can take several months to achieve a single gene knockout because they rely on aphids going through a cycle of sexual reproduction, and aphids often lack strong, consistent levels of knockdown when fed or injected with molecules that induce an RNA interference (RNAi) response. In the hopes of addressing these challenges, we attempted to adapt a new method called symbiont-mediated RNAi (smRNAi) for use in aphids. smRNAi involves engineering a bacterial symbiont of the insect to continuously supply double-stranded RNA (dsRNA) inside the insect body. This approach has been successful in thrips, kissing bugs, and honeybees. We engineered the laboratory Escherichia coli strain HT115 and the native aphid symbiont Serratia symbiotica CWBI-2.3^T to produce dsRNA inside the gut of the pea aphid (Acyrthosiphon pisum) targeting salivary effector protein (C002) or ecdysone receptor genes. For C002 assays, we also tested co-knockdown with an aphid nuclease (Nuc1) to reduce RNA degradation. However, we found that smRNAi was not a reliable method for aphid gene knockdown under our conditions. We were unable to consistently achieve the expected phenotypic changes with either target. However, we did see indications that elements of the RNAi pathway were modestly upregulated, and expression of some targeted genes appeared to be somewhat reduced in some trials. We conclude with a discussion of the possible avenues through which smRNAi, and aphid RNAi in general, could be improved in the future.

Dietary Supplement of Grape Wastes Enhances Honeybee Immune System and Reduces Deformed Wing Virus (DWV) Load

Ingredients rich in phenolic compounds and antioxidants of winemaking wastes, which play an important role in the prevention of various diseases and the control of viruses, are being explored. Currently, there is a concern about honeybee population loss, with deformed wing virus (DWV) being the most common virus infecting apiaries and one of the main causes of honeybee decline. Hence, the effect of grape pomace powder (GPP) as a dietary supplement to enhance the immune system of honeybees affected by DWV was evaluated. The characteristics of the ingredient GPP, obtained by spray-drying, revealed a high anthocyanin content (1102.45 mg 100 g^-1), and it was applied at doses of 0.5, 1, 2.5 and 5% as a dietary supplement for bees infected by DWV. The results showed that the GPP treatments strengthened the immune response of honeybees against DWV. Moreover, the expression of the Relish gene was significantly higher in bees fed with GPP compared to the infected control. This study, which is framed in the search of food waste valorization for environmental sustainability, proves the feasibility of using grape wastes as dietary supplements for pollinators, and provides knowledge of the influence of polyphenols on the expression profiles of immune-related genes in honeybees.

Combined transcriptome and metabolite profiling analyses provide insights into the chronic toxicity of carbaryl and acetamiprid to Apis mellifera larvae

Despite many studies have revealed that developing honey bee (Apis mellifera) larvae are posting a high risk on exposure to insecticides, the toxicology information on bee larvae remain limited. The present study demonstrated the first assessment of the effects of no observed adverse effect concentration (NOAEC) of carbaryl (CR) and acetamiprid (ACE) on transcriptome and metabolome in honeybee larvae reared in vitro. Chronic exposure to carbaryl caused transcriptional disorders associated with oxidative stress. In addition, a series of metabolic homeostasis were disrupted by carbaryl stress, such amino acid metabolism, purine and pyrimidine metabolism and flavone and flavonol biosynthesis. The activities of enzymic biomarkers including GST, P450, CAT, AChE and SOD were not influenced by ACE stress, while the CR exposure slightly decreased the activity of CAT and SOD. Our results clearly show that ACE and CR display different potential to modulate transcriptome and metabolome associated with their different toxicity against bee larvae.

Fertility costs of cryptic viral infections in a model social insect

Declining insect populations emphasize the importance of understanding the drivers underlying reductions in insect fitness. Here, we investigated viruses as a threat to social insect reproduction, using honey bees as a model species. We report that in two independent surveys (N = 93 and N = 54, respectively) of honey bee (Apis mellifera) queens taken from a total of ten beekeeping operations across British Columbia, high levels of natural viral infection are associated with decreased ovary mass. Failed (poor quality) queens displayed higher levels of viral infection, reduced sperm viability, smaller ovaries, and altered ovary protein composition compared to healthy queens. We experimentally infected queens with Israeli acute paralysis virus (IAPV) and found that the ovary masses of IAPV-injected queens were significantly smaller than control queens, demonstrating a causal relationship between viral infection and ovary size. Queens injected with IAPV also had significantly lower expression of vitellogenin, the main source of nutrition deposited into developing oocytes, and higher levels of heat-shock proteins, which are part of the honey bee’s antiviral response. This work together shows that viral infections occurring naturally in the field are compromising queen reproductive success.

Paromomycin Reduces Vairimorpha (Nosema) ceranae Infection in Honey Bees but Perturbs Microbiome Levels and Midgut Cell Function

Paromomycin is a naturally occurring aminoglycoside antibiotic that has effects on both prokaryotic and eukaryotic microbes. However, previous reports have indicated that it has little effect on microsporidia, including Vairimorpha (Nosema) ceranae, in cell culture models. V. ceranae is one of a number of microsporidia species that cause disease in honey bees and substantial efforts to find new treatment strategies for bees that are infected with these pathogens are ongoing. When testing compounds for potential activity against V. ceranae in whole organisms, we found that paromomycin reduces the infection intensity of this parasite. Critically, the necessary doses of paromomycin have high activity against the bacteria of the honey bee microbiome and cause evident stress in bees. Microsporidia have been shown to lack an essential binding site on the ribosome that is known to allow for maximal inhibition by paromomycin. Thus, it is possible that paromomycin impacts parasite levels through non-cell autonomous effects on microsporidia infection levels via effects on the microbiome or midgut cellular function. As paromomycin treatment could cause widespread honey bee health issues in agricultural settings, it does not represent an appropriate anti-microsporidia agent for use in the field.

Meta-Analysis of Immune Induced Gene Expression Changes in Diverse Drosophila melanogaster Innate Immune Responses

Simple Summary

Organisms can be infected by a wide range of pathogens, including bacteria, viruses, and parasites. Following infection, the host mounts an immune response to attempt to eliminate the pathogen. These responses are often specific to the type of pathogen and mediated by the expression of specialized genes. We have characterized the expression changes induced in host Drosophila fruit flies following infection by multiple types of pathogens, and identified a small number of genes that show expression changes in each infection. This includes genes that are known to be involved in pathogen resistance, and others that have not been previously studied as immune response genes. These findings provide new insight into transcriptional changes that accompany Drosophila immunity. They may suggest possible roles for the differentially expressed genes in innate immune responses to diverse classes of pathogens, and serve to identify candidate genes for further empirical study of these processes.

Abstract

Organisms are commonly infected by a diverse array of pathogens and mount functionally distinct responses to each of these varied immune challenges. Host immune responses are characterized by the induction of gene expression, however, the extent to which expression changes are shared among responses to distinct pathogens is largely unknown. To examine this, we performed meta-analysis of gene expression data collected from Drosophila melanogaster following infection with a wide array of pathogens. We identified 62 genes that are significantly induced by infection. While many of these infection-induced genes encode known immune response factors, we also identified 21 genes that have not been previously associated with host immunity. Examination of the upstream flanking sequences of the infection-induced genes lead to the identification of two conserved enhancer sites. These sites correspond to conserved binding sites for GATA and nuclear factor κB (NFκB) family transcription factors and are associated with higher levels of transcript induction. We further identified 31 genes with predicted functions in metabolism and organismal development that are significantly downregulated following infection by diverse pathogens. Our study identifies conserved gene expression changes in Drosophila melanogaster following infection with varied pathogens, and transcription factor families that may regulate this immune induction.

Aryl hydrocarbon receptor and Krüppel like factor 10 mediate a transcriptional axis modulating immune homeostasis in mosquitoes

Immune responses require delicate controls to maintain homeostasis while executing effective defense. Aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor. The Krüppel-like factor 10 (KLF10) is a C2H2 zinc-finger containing transcription factor. The functions of mosquito AhR and KLF10 have not been characterized. Here we show that AhR and KLF10 constitute a transcriptional axis to modulate immune responses in mosquito Anopheles gambiae. The manipulation of AhR activities via agonists or antagonists repressed or enhanced the mosquito antibacterial immunity, respectively. KLF10 was recognized as one of the AhR target genes in the context. Phenotypically, silencing KLF10 reversed the immune suppression caused by the AhR agonist. The transcriptome comparison revealed that silencing AhR and KLF10 plus challenge altered the expression of 2245 genes in the same way. The results suggest that KLF10 is downstream of AhR in a transcriptional network responsible for immunomodulation. This AhR–KLF10 axis regulates a set of genes involved in metabolism and circadian rhythms in the context. The axis was required to suppress the adverse effect caused by the overactivation of the immune pathway IMD via the inhibitor gene Caspar silencing without a bacterial challenge. These results demonstrate that the AhR–KLF10 axis mediates an immunoregulatory transcriptional network as a negative loop to maintain immune homeostasis.

Critical View on the Importance of Host Defense Strategies on Virus Distribution of Bee Viruses: What Can We Learn from SARS-CoV-2 Variants?

Bees, both wild and domesticated ones, are hosts to a plethora of viruses, with most of them infecting a wide range of bee species and genera. Although viral discovery and research on bee viruses date back over 50 years, the last decade is marked by a surge of new studies, new virus discoveries, and reports on viral transmission in and between bee species. This steep increase in research on bee viruses was mainly initiated by the global reports on honeybee colony losses and the worldwide wild bee decline, where viruses are regarded as one of the main drivers. While the knowledge gained on bee viruses has significantly progressed in a short amount of time, we believe that integration of host defense strategies and their effect on viral dynamics in the multi-host viral landscape are important aspects that are currently still missing. With the large epidemiological dataset generated over the last two years on the SARS-CoV-2 pandemic, the role of these defense mechanisms in shaping viral dynamics has become eminent. Integration of these dynamics in a multi-host system would not only greatly aid the understanding of viral dynamics as a driver of wild bee decline, but we believe bee pollinators and their viruses provide an ideal system to study the multi-host viruses and their epidemiology.

Context-Dependent Viral Transgenerational Immune Priming in Honey Bees (Hymenoptera: Apidae)

Transgenerational immune priming is the process of increased resistance to infection in offspring due to parental pathogen exposure. Honey bees (Apis mellifera L. (Hymenoptera: Apidae)) are hosts to multiple pathogens, and this complex immune function could help protect against overwhelming infection. Honey bees have demonstrated transgenerational immune priming for the bacterial pathogen Paenibacillus larvae; however, evidence for viral transgenerational immune priming is lacking across insects in general. Here we test for the presence of transgenerational immune priming in honey bees with Deformed wing virus (DWV) by injecting pupae from DWV-exposed queens and measuring virus titer and immune gene expression. Our data suggest that there is evidence for viral transgenerational immune priming in honey bees, but it is highly context-dependent based on route of maternal exposure and potentially host genetics or epigenetic factors.

Establishing RNAi for basic research and pest control and identification of the most efficient target genes for pest control: a brief guide

RNA interference (RNAi) has emerged as a powerful tool for knocking-down gene function in diverse taxa including arthropods for both basic biological research and application in pest control. The conservation of the RNAi mechanism in eukaryotes suggested that it should-in principle-be applicable to most arthropods. However, practical hurdles have been limiting the application in many taxa. For instance, species differ considerably with respect to efficiency of dsRNA uptake from the hemolymph or the gut. Here, we review some of the most frequently encountered technical obstacles when establishing RNAi and suggest a robust procedure for establishing this technique in insect species with special reference to pests. Finally, we present an approach to identify the most effective target genes for the potential control of agricultural and public health pests by RNAi.

Chemical Stimulants and Stressors Impact the Outcome of Virus Infection and Immune Gene Expression in Honey Bees (Apis mellifera)

Western honey bees (Apis mellifera) are ecologically, agriculturally, and economically important plant pollinators. High average annual losses of honey bee colonies in the US have been partially attributed to agrochemical exposure and virus infections. To examine the potential negative synergistic impacts of agrochemical exposure and virus infection, as well as the potential promise of phytochemicals to ameliorate the impact of pathogenic infections on honey bees, we infected bees with a panel of viruses (i.e., Flock House virus, deformed wing virus, or Sindbis virus) and exposed to one of three chemical compounds. Specifically, honey bees were fed sucrose syrup containing: (1) thyme oil, a phytochemical and putative immune stimulant, (2) fumagillin, a beekeeper applied fungicide, or (3) clothianidin, a grower-applied insecticide. We determined that virus abundance was lower in honey bees fed 0.16 ppb thyme oil augmented sucrose syrup, compared to bees fed sucrose syrup alone. Parallel analysis of honey bee gene expression revealed that honey bees fed thyme oil augmented sucrose syrup had higher expression of key RNAi genes (argonaute-2 and dicer-like), antimicrobial peptide expressing genes (abaecin and hymenoptaecin), and vitellogenin, a putative honey bee health and age indicator, compared to bees fed only sucrose syrup. Virus abundance was higher in bees fed fumagillin (25 ppm or 75 ppm) or 1 ppb clothianidin containing sucrose syrup relative to levels in bees fed only sucrose syrup. Whereas, honey bees fed 10 ppb clothianidin had lower virus levels, likely because consuming a near lethal dose of insecticide made them poor hosts for virus infection. The negative impact of fumagillin and clothianidin on honey bee health was indicated by the lower expression of argonaute-2, dicer-like, abaecin, and hymenoptaecin, and vitellogenin. Together, these results indicate that chemical stimulants and stressors impact the outcome of virus infection and immune gene expression in honey bees.

Honey bee sHSP are responsive to diverse proteostatic stresses and potentially promising biomarkers of honey bee stress

The pollination services provided by the honey bee are critical in both natural and agricultural ecosystems. Honey bee colonies in the United States have suffered from an increased rate of die-off in recent years, stemming from a complex set of interacting stresses that remain poorly described. Defining specific common cellular processes and cellular stress responses impacted by multiple stressors represent a key step in understanding these synergies. Proteotoxic stresses negatively impact protein synthesis, folding, and degradation. Diverse proteotoxic stresses induce expression of genes encoding small heat shock proteins (sHSP) of the expanded lethal (2) essential for life (l(2)efl) gene family. In addition to upregulation by the Integrated Stress Response (ISR), the Heat Shock Response (HSR), and the Oxidative Stress Response (OSR), our data provide first evidence that sHSP genes are upregulated by the Unfolded Protein Response (UPR). As these genes appear to be part of a core stress response that could serve as a useful biomarker for cellular stress in honey bees, we designed and tested an RT-LAMP assay to detect increased l(2)efl gene expression in response to heat-stress. While this assay provides a powerful proof of principle, further work will be necessary to link changes in sHSP gene expression to colony-level outcomes, to adapt our preliminary assay into a Point of Care Testing (POCT) assay appropriate for use as a diagnostic tool for use in the field, and to couple assay results to management recommendations.

Honey Bee (Apis mellifera) Immunity

At the individual level, honey bees (Apis mellifera) rely on innate immunity, which operates through cellular and humoral mechanisms, to defend themselves against infectious agents and parasites. At the colony level, honey bees have developed collective defense mechanisms against pathogens and pests, such as hygienic and grooming behaviors. An understanding of the immune responses of honey bees is critical to implement strategies to reduce mortality and increase colony productivity. The major components and mechanisms of individual and social immunity of honey bees are discussed in this review.

The Honey Bee Gene Bee Antiviral Protein-1 Is a Taxonomically Restricted Antiviral Immune Gene

Insects have evolved a wide range of strategies to combat invading pathogens, including viruses. Genes that encode proteins involved in immune responses often evolve under positive selection due to their co-evolution with pathogens. Insect antiviral defense includes the RNA interference (RNAi) mechanism, which is triggered by recognition of non-self, virally produced, double-stranded RNAs. Indeed, insect RNAi genes (e.g., dicer and argonaute-2) are under high selective pressure. Honey bees (Apis mellifera) are eusocial insects that respond to viral infections via both sequence specific RNAi and a non-sequence specific dsRNA triggered pathway, which is less well-characterized. A transcriptome-level study of virus-infected and/or dsRNA-treated honey bees revealed increased expression of a novel antiviral gene, GenBank: MF116383, and in vivo experiments confirmed its antiviral function. Due to in silico annotation and sequence similarity, MF116383 was originally annotated as a probable cyclin-dependent serine/threonine-protein kinase. In this study, we confirmed that MF116383 limits virus infection, and carried out further bioinformatic and phylogenetic analyses to better characterize this important gene-which we renamed bee antiviral protein-1 (bap1). Phylogenetic analysis revealed that bap1 is taxonomically restricted to Hymenoptera and Blatella germanica (the German cockroach) and that the majority of bap1 amino acids are evolving under neutral selection. This is in-line with the results from structural prediction tools that indicate Bap1 is a highly disordered protein, which likely has relaxed structural constraints. Assessment of honey bee gene expression using a weighted gene correlation network analysis revealed that bap1 expression was highly correlated with several immune genes-most notably argonaute-2. The coexpression of bap1 and argonaute-2 was confirmed in an independent dataset that accounted for the effect of virus abundance. Together, these data demonstrate that bap1 is a taxonomically restricted, rapidly evolving antiviral immune gene. Future work will determine the role of bap1 in limiting replication of other viruses and examine the signal cascade responsible for regulating the expression of bap1 and other honey bee antiviral defense genes, including coexpressed ago-2, and determine whether the virus limiting function of bap1 acts in parallel or in tandem with RNAi.

Thermal stress induces tissue damage and a broad shift in regenerative signaling pathways in the honey bee digestive tract

Honey bee colonies in the USA have suffered from increased die-off in the last few years with a complex set of interacting stresses playing a key role. With changing climate, an increase in the frequency of severe weather events, such as heat waves, is anticipated. Understanding how these changes may contribute to stress in honey bees is crucial. Individual honey bees appear to have a high capacity to endure thermal stress. One reason for this high-level endurance is likely their robust heat shock response (HSR), which contributes to thermotolerance at the cellular level. However, less is known about other mechanisms of thermotolerance, especially those operating at the tissue level. To elucidate other determinants of resilience in this species, we used thermal stress coupled with RNAseq and identified broad transcriptional remodeling of a number of key signaling pathways in the honey bee, including those pathways known to be involved in digestive tract regeneration in the fruit fly such as the Hippo and JAK/STAT pathways. We also observed cell death and shedding of epithelial cells, which likely leads to induction of this regenerative transcriptional program. We found that thermal stress affects many of these pathways in other tissues, suggesting a shared program of damage response. This study provides important foundational characterization of the tissue damage response program in this key pollinating species. In addition, our data suggest that a robust regeneration program may also be a critical contributor to thermotolerance at the tissue level, a possibility which warrants further exploration in this and other species.

Investigating Virus-Host Interactions in Cultured Primary Honey Bee Cells

Honey bee (Apis mellifera) health is impacted by viral infections at the colony, individual bee, and cellular levels. To investigate honey bee antiviral defense mechanisms at the cellular level we further developed the use of cultured primary cells, derived from either larvae or pupae, and demonstrated that these cells could be infected with a panel of viruses, including common honey bee infecting viruses (i.e., sacbrood virus (SBV) and deformed wing virus (DWV)) and an insect model virus, Flock House virus (FHV). Virus abundances were quantified over the course of infection. The production of infectious virions in cultured honey bee pupal cells was demonstrated by determining that naïve cells became infected after the transfer of deformed wing virus or Flock House virus from infected cell cultures. Initial characterization of the honey bee antiviral immune responses at the cellular level indicated that there were virus-specific responses, which included increased expression of bee antiviral protein-1 (GenBank: MF116383) in SBV-infected pupal cells and increased expression of argonaute-2 and dicer-like in FHV-infected hemocytes and pupal cells. Additional studies are required to further elucidate virus-specific honey bee antiviral defense mechanisms. The continued use of cultured primary honey bee cells for studies that involve multiple viruses will address this knowledge gap.

Conservation genomics reveals pesticide and pathogen exposure in the declining bumble bee Bombus terricola

In recent years, many pollinators have experienced large population declines, which threaten food security and the stability of natural ecosystems. Bumble bees are particularly important because their ability to “buzz” pollinate and tolerate cooler conditions make them critical pollinators for certain plants and regions. Here, we apply a conservation genomics approach to study the vulnerable Bombus terricola. We sequenced RNA from 30 worker abdomens, 18 of which were collected from agricultural sites and 12 of which were collected from nonagricultural sites. We found transcriptional signatures associated with exposure to insecticides, with gene expression patterns suggesting that bumble bees were exposed to neonicotinoids and/or fipronil—two compounds known to negatively impact bees. We also found transcriptional signatures associated with pathogen infections. In addition to the transcriptomic analysis, we carried out a metatranscriptomic analysis and detected five pathogens in the abdomens of workers, three of which are common in managed honey bee and bumble bee colonies. Our conservation genomics study provides functional support for the role of pesticides and pathogen spillover in the decline of B. terricola. We demonstrate that conservation genomics is an invaluable tool which allows researchers to quantify the effects of multiple stressors that impact pollinator populations in the wild.

Non-Target Effects of dsRNA Molecules in Hemipteran Insects

Insect pest control by RNA interference (RNAi)-mediated gene expression knockdown can be undermined by many factors, including small sequence differences between double-stranded RNA (dsRNA) and the target gene. It can also be compromised by effects that are independent of the dsRNA sequence on non-target organisms (known as sequence-non-specific effects). This study investigated the species-specificity of RNAi in plant sap-feeding hemipteran pests. We first demonstrated sequence-non-specific suppression of aphid feeding by dsRNA at dietary concentrations ≥0.5 µg µL⁻¹. Then we quantified the expression of NUC (nuclease) genes in insects administered homologous dsRNA (with perfect sequence identity to the target species) or heterologous dsRNA (generated against a related gene of non-identical sequence in a different insect species). For the aphids Acyrthosiphon pisum and Myzus persicae, significantly reduced NUC expression was obtained with the homologous but not heterologous dsRNA at 0.2 µg µL⁻¹, despite high dsNUC sequence identity. Follow-up experiments demonstrated significantly reduced expression of NUC genes in the whitefly Bemisia tabaci and mealybug Planococcus maritimus administered homologous dsNUCs, but not heterologous aphid dsNUCs. Our demonstration of inefficient expression knockdown by heterologous dsRNA in these insects suggests that maximal dsRNA sequence identity is required for RNAi targeting of related pest species, and that heterologous dsRNAs at appropriate concentrations may not be a major risk to non-target sap-feeding hemipterans.

RNA Interference-Mediated Knockdown of Genes Encoding Spore Wall Proteins Confers Protection against Nosema ceranae Infection in the European Honey Bee, Apis mellifera

Nosema ceranae (Opisthosporidia: Microsporidia) is an emergent intracellular parasite of the European honey bee (Apis mellifera) and causes serious Nosema disease which has been associated with worldwide honey bee colony losses. The only registered treatment for Nosema disease is fumagillin-b, and this has raised concerns about resistance and off-target effects. Fumagillin-B is banned from use in honey bee colonies in many countries, particularly in Europe. As a result, there is an urgent need for new and effective therapeutic options to treat Nosema disease in honey bees. An RNA interference (RNAi)-based approach can be a potent strategy for controlling diseases in honey bees. We explored the therapeutic potential of silencing the sequences of two N. ceranae encoded spore wall protein (SWP) genes by means of the RNAi-based methodology. Our study revealed that the oral ingestion of dsRNAs corresponding to SWP8 and SWP12 used separately or in combination could lead to a significant reduction in spore load, improve immunity, and extend the lifespan of N. ceranae-infected bees. The results from the work completed here enhance our understanding of honey bee host responses to microsporidia infection and highlight that RNAi-based therapeutics are a promising treatment for honey bee diseases.

Comparative response of Spodoptera litura challenged per os with Serratia marcescens strains differing in virulence

Host plays an important role in influencing virulence of a pathogen and efficacy of a biopesticide. The present study was aimed to characterize the possible factors present in Spodoptera litura that influenced pathogenecity of orally ingested S. marcescens strains, differing in their virulence. Fifth instar larvae of S. litura responded differently as challenged by two Serratia marcescens strains, SEN (virulent strain, LC₅₀ 7.02 10³ cfu/ml) and ICC-4 (non-virulent strain, LC₅₀ 1.19 10¹² cfu/ml). Considerable increase in activity of lytic enzymes protease and phospholipase was recorded in the gut and hemolymph of larvae fed on diet supplemented with S. marcescens strain ICC-4 as compared to the larvae treated with S. marcescens strain SEN. However, a significant up-regulation of antioxidative enzymes SOD (in foregut and midgut), CAT (in the midgut) and GST (in the foregut and hemolymph) was recorded in larvae fed on diet treated with the virulent S. marcescens strain SEN in comparison to larvae fed on diet treated with the non-virulent S. marcescens strain ICC-4. Activity of defense related enzymes lysozyme and phenoloxidase activity were also higher in the hemolymph of larvae fed with diet treated with S. marcescens strain SEN as compared to hemolymph of S. marcescens strain ICC-4 treated larvae. More number of over-expressed proteins was observed in the gut and hemolymph of S. marcescens strains ICC-4 and SEN treated larvae, respectively. Identification of the selected differentially expressed proteins indicated induction of proteins involved in insect innate immune response (Immunoglobulin I-set domain, Apolipophorin III, leucine rich repeat and Titin) in S. marcescens strain SEN treated larvae. Over-expression of two proteins, actin related protein and mt DNA helicase, were noted in S. marcescens treated larvae with very high levels observed in the non-virulent strain. Up-regulation of homeobox protein was noted only in S. marcescens strain ICC-4 challenged larvae. This study indicated that ingestion of non-virulent S. marcescens strain ICC-4 induced strong immune response in insect gut while there was weak response to the virulent S. marcescens strain SEN which probably resulted in difference in their virulence.

The Herbal Supplements NOZEMAT HERB® and NOZEMAT HERB PLUS®: An Alternative Therapy for N. ceranae Infection and Its Effects on Honey Bee Strength and Production Traits

Honey bees (Apis mellifera L.) are the most effective pollinators for different crops and wild flowering plants, thus maintaining numerous ecosystems in the world. However, honey bee colonies often suffer from stress or even death due to various pests and diseases. Among the latter, nosemosis is considered to be one of the most common diseases, causing serious damage to beekeeping every year. Here, we present, for the first time, the effects from the application of the herbal supplements NOZEMAT HERB^® (NH) and NOZEMAT HERB PLUS^® (NHP) for treating N. ceranae infection and positively influencing the general development of honey bee colonies. To achieve this, in autumn 2019, 45 colonies were selected based on the presence of N. ceranae infections. The treatment was carried out for 11 months (August 2019–June 2020). All colonies were sampled pre- and post-treatment for the presence of N. ceranae by means of light microscopy and PCR analysis. The honey bee colonies’ performance and health were evaluated pre- and post-treatment. The obtained results have shown that both supplements have exhibited statistically significant biological activity against N. ceranae in infected apiaries. Considerable enhancement in the strength of honey bee colonies and the amount of sealed workers was observed just one month after the application of NH and NHP. Although the mechanisms of action of NH and NHP against N. ceranae infection are yet to be completely elucidated, our results suggest a new holistic approach as an alternative therapy to control nosemosis and to improve honey bee colonies’ performance and health.

Transcriptomic Responses of the Honey Bee Brain to Infection with Deformed Wing Virus

Managed colonies of European honey bees (Apis mellifera) are under threat from Varroa destructor mite infestation and infection with viruses vectored by mites. In particular, deformed wing virus (DWV) is a common viral pathogen infecting honey bees worldwide that has been shown to induce behavioral changes including precocious foraging and reduced associative learning. We investigated how DWV infection of bees affects the transcriptomic response of the brain. The transcriptomes of individual brains were analyzed using RNA-Seq after experimental infection of newly emerged adult bees with DWV. Two analytical methods were used to identify differentially expressed genes from the ~15,000 genes in the Apis mellifera genome. The 269 genes that had increased expression in DWV infected brains included genes involved in innate immunity such as antimicrobial peptides (AMPs), Ago2, and Dicer. Single bee brain NMR metabolomics methodology was developed for this work and indicates that proline is strongly elevated in DWV infected brains, consistent with the increased presence of the AMPs abaecin and apidaecin. The 1361 genes with reduced expression levels includes genes involved in cellular communication including G-protein coupled, tyrosine kinase, and ion-channel regulated signaling pathways. The number and function of the downregulated genes suggest that DWV has a major impact on neuron signaling that could explain DWV related behavioral changes.

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

Simple Summary

Honey bees and other pollinators are threatened by numerous stressors, including virus infections. Currently, no effective treatments are available, stressing the importance of natural defenses. These defenses may be enhanced through selective breeding. This study sought to evaluate the potential for breeding, while also testing a few potential mechanisms of natural immune responses and assessing how widespread viruses are in commercial honey bee queens in the U.S. We identified significant differences in survival of virus infection among and within U.S honey bee stocks, indicating that selective breeding may be able to decrease the virus susceptibility of honey bees. Survival differences may be related to differences in the natural immune system of honey bees and could relate to how much virus stress bees have experienced in the past.

Abstract

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.

The Interplay Between Viruses and RNAi Pathways in Insects

As an overarching immune mechanism, RNA interference (RNAi) displays pathogen specificity and memory via different pathways. The small interfering RNA (siRNA) pathway is the primary antiviral defense mechanism against RNA viruses of insects and plays a lesser role in defense against DNA viruses. Reflecting the pivotal role of the siRNA pathway in virus selection, different virus families have independently evolved unique strategies to counter this host response, including protein-mediated, decoy RNA-based, and microRNA-based strategies. In this review, we outline the interplay between insect viruses and the different pathways of the RNAi antiviral response; describe practical application of these interactions for improved expression systems and for pest and disease management; and highlight research avenues for advancement of the field.

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.

Egg transcriptome profile responds to maternal virus infection in honey bees, Apis mellifera

Trans-generational disease effects include vertical pathogen transmission but also immune priming to enhance offspring immunity. Accordingly, the survival consequences of maternal virus infection can vary and its molecular consequences during early development are poorly understood. The honey bee queen is long-lived and represents the central hub for vertical virus transmission as the sole reproductive individual in her colony. Even though virus symptoms in queens are mild, viral infection may have severe consequences for the offspring. Thus, transcriptome patterns during early developmental are predicted to respond to maternal virus infection. To test this hypothesis, gene expression patterns were compared among pooled honey bee eggs laid by queens that were either infected with Deformed wing virus (DWV¹), Sacbrood virus (SBV²), both viruses (DWV and SBV), or no virus. Whole transcriptome analyses revealed significant expression differences of a few genes, some of which have hitherto no known function. Despite the paucity of single gene effects, functional enrichment analyses revealed numerous biological processes in the embryos to be affected by virus infection. Effects on several regulatory pathways were consistent with maternal responses to virus infection and correlated with responses to DWV and SBV in honey bee larvae and pupae. Overall, effects on egg transcriptome patterns were specific to each virus and the results of dual-infection samples suggested synergistic effects of DWV and SBV. We interpret our results as consequences of maternal infections. Thus, this first study to document and characterize virus-associated changes in the transcriptome of honey bee eggs represents an important contribution to understanding trans-generational virus effects, although more in-depth studies are needed to understand the detailed mechanisms of how viruses affect honey bee embryos.

Longitudinal monitoring of honey bee colonies reveals dynamic nature of virus abundance and indicates a negative impact of Lake Sinai virus 2 on colony health

Honey bees (Apis mellifera) are important pollinators of plants, including those that produce nut, fruit, and vegetable crops. Therefore, high annual losses of managed honey bee colonies in the United States and many other countries threaten global agriculture. Honey bee colony deaths have been associated with multiple abiotic and biotic factors, including pathogens, but the impact of virus infections on honey bee colony population size and survival are not well understood. To further investigate seasonal patterns of pathogen presence and abundance and the impact of viruses on honey bee colony health, commercially managed colonies involved in the 2016 California almond pollination event were monitored for one year. At each sample date, colony health and pathogen burden were assessed. Data from this 50-colony cohort study illustrate the dynamic nature of honey bee colony health and the temporal patterns of virus infection. Black queen cell virus, deformed wing virus, sacbrood virus, and the Lake Sinai viruses were the most readily detected viruses in honey bee samples obtained throughout the year. Analyses of virus prevalence and abundance revealed pathogen-specific trends including the overall increase in deformed wing virus abundance from summer to fall, while the levels of Lake Sinai virus 2 (LSV2) decreased over the same time period. Though virus prevalence and abundance varied in individual colonies, analyses of the overall trends reveal correlation with sample date. Total virus abundance increased from November 2015 (post-honey harvest) to the end of the almond pollination event in March 2016, which coincides with spring increase in colony population size. Peak total virus abundance occurred in late fall (August and October 2016), which correlated with the time period when the majority of colonies died. Honey bee colonies with larger populations harbored less LSV2 than weaker colonies with smaller populations, suggesting an inverse relationship between colony health and LSV2 abundance. Together, data from this and other longitudinal studies at the colony level are forming a better understanding of the impact of viruses on honey bee colony losses.

Injection of seminal fluid into the hemocoel of honey bee queens (Apis mellifera) can stimulate post-mating changes

Honey bee queens undergo dramatic behavioral (e.g., reduced sexual receptivity), physiological (e.g., ovary activation, ovulation, and modulation of pheromone production) and transcriptional changes after they complete mating. To elucidate how queen post-mating changes are influenced by seminal fluid, the non-spermatozoa-containing component of semen, we injected queens with semen or seminal fluid alone. We assessed queen sexual receptivity (as measured by likelihood to take mating flights), ovary activation, worker retinue response (which is influenced by queen pheromone production), and transcriptional changes in queen abdominal fat body and brain tissues. Injection with either seminal fluid or semen resulted in decreased sexual receptivity, increased attractiveness of queens to workers, and altered expression of several genes that are also regulated by natural mating in queens. The post-mating and transcriptional changes of queens receiving seminal fluid were not significantly different from queens injected with semen, suggesting that components in seminal fluid, such as seminal fluid proteins, are largely responsible for stimulating post-mating changes in queens.

Identification of a novel host protein interacting with the structural protein VP2 of deformed wing virus by yeast two-hybrid screening

Deformed wing virus (DWV) interacting with Varroa destructor is a possible cause of honeybee colony mortality. VP2 is the structural protein of DWV but its function remains unknown. To clarify the function of VP2 and screen for novel binding proteins that interact with VP2, we carried out a membrane protein yeast two-hybrid screening using VP2 as bait. Subsequently, the interaction between VP2 and the host interacting protein [heat shock protein 10 (Hsp10)] was further verified using glutathione S-transferase pull-down assay in vitro and co-immunoprecipitation assay in cells. Furthermore, fluorescence confocal microscopy revealed that VP2 and Hsp10 were mainly co-localized in the cytoplasm. Using real-time polymerase chain reaction, we found that Hsp10 expression in DWV-infected worker honey bees were downregulated compared with that in healthy honey bees. Additionally, we showed that overexpression of VP2 protein could reduce the expression of Hsp10. These results suggest that Hsp10 plays a vital role in host immunity and antiviral effects.

Immunopathology and immune homeostasis during viral infection in insects

Organisms clear infections by mounting an immune response that is normally turned off once the pathogens have been cleared. However, sometimes this immune response is not properly or timely arrested, resulting in the host damaging itself. This immune dysregulation may be referred to as immunopathology. While our knowledge of immune and metabolic pathways in insects, particularly in response to viral infections, is growing, little is known about the mechanisms that regulate this immune response and hence little is known about immunopathology in this important and diverse group of organisms. In this chapter we focus both on documenting the molecular mechanisms described involved in restoring immune homeostasis in insects after viral infections and on identifying potential mechanisms for future investigation. We argue that learning about the immunopathological consequences of an improperly regulated immune response in insects will benefit both insect and human health.

Assessing the Risks of Topically Applied dsRNA-Based Products to Non-target Arthropods

RNA interference (RNAi) is a powerful technology that offers new opportunities for pest control through silencing of genes that are essential for the survival of arthropod pests. The approach relies on sequence-specificity of applied double-stranded (ds) RNA that can be designed to have a very narrow spectrum of both the target gene product (RNA) as well as the target organism, and thus allowing highly targeted pest control. Successful RNAi has been reported from a number of arthropod species belonging to various orders. Pest control may be achieved by applying dsRNA as foliar sprays. One of the main concerns related to the use of dsRNA is adverse environmental effects particularly on valued non-target species. Arthropods form an important part of the biodiversity in agricultural landscapes and contribute important ecosystem services. Consequently, environmental risk assessment (ERA) for potential impacts that plant protection products may have on valued non-target arthropods is legally required prior to their placement on the market. We describe how problem formulation can be used to set the context and to develop plausible pathways on how the application of dsRNA-based products could harm valued non-target arthropod species, such as those contributing to biological pest control. The current knowledge regarding the exposure to and the hazard posed by dsRNA in spray products for non-target arthropods is reviewed and suggestions are provided on how to select the most suitable test species and to conduct laboratory-based toxicity studies that provide robust, reliable and interpretable results to support the ERA.

Pesticide-Virus Interactions in Honey Bees: Challenges and Opportunities for Understanding Drivers of Bee Declines

Honey bees are key agricultural pollinators, but beekeepers continually suffer high annual colony losses owing to a number of environmental stressors, including inadequate nutrition, pressures from parasites and pathogens, and exposure to a wide variety of pesticides. In this review, we examine how two such stressors, pesticides and viruses, may interact in additive or synergistic ways to affect honey bee health. Despite what appears to be a straightforward comparison, there is a dearth of studies examining this issue likely owing to the complexity of such interactions. Such complexities include the wide array of pesticide chemical classes with different modes of actions, the coupling of many bee viruses with ectoparasitic Varroa mites, and the intricate social structure of honey bee colonies. Together, these issues pose a challenge to researchers examining the effects pesticide-virus interactions at both the individual and colony level.

The Heat Shock Response in the Western Honey Bee (Apis mellifera) is Antiviral

Honey bees (Apismellifera) are an agriculturally important pollinator species that live in easily managed social groups (i.e., colonies). Unfortunately, annual losses of honey bee colonies in many parts of the world have reached unsustainable levels. Multiple abiotic and biotic stressors, including viruses, are associated with individual honey bee and colony mortality. Honey bees have evolved several antiviral defense mechanisms including conserved immune pathways (e.g., Toll, Imd, JAK/STAT) and dsRNA-triggered responses including RNA interference and a non-sequence specific dsRNA-mediated response. In addition, transcriptome analyses of virus-infected honey bees implicate an antiviral role of stress response pathways, including the heat shock response. Herein, we demonstrate that the heat shock response is antiviral in honey bees. Specifically, heat-shocked honey bees (i.e., 42 °C for 4 h) had reduced levels of the model virus, Sindbis-GFP, compared with bees maintained at a constant temperature. Virus-infection and/or heat shock resulted in differential expression of six heat shock protein encoding genes and three immune genes, many of which are positively correlated. The heat shock protein encoding and immune gene transcriptional responses observed in virus-infected bees were not completely recapitulated by administration of double stranded RNA (dsRNA), a virus-associated molecular pattern, indicating that additional virus-host interactions are involved in triggering antiviral stress response pathways.

Meta-analysis of honey bee neurogenomic response links Deformed wing virus type A to precocious behavioral maturation

Crop pollination by the western honey bee Apis mellifera is vital to agriculture but threatened by alarmingly high levels of colony mortality, especially in Europe and North America. Colony loss is due, in part, to the high viral loads of Deformed wing virus (DWV), transmitted by the ectoparasitic mite Varroa destructor, especially throughout the overwintering period of a honey bee colony. Covert DWV infection is commonplace and has been causally linked to precocious foraging, which itself has been linked to colony loss. Taking advantage of four brain transcriptome studies that unexpectedly revealed evidence of covert DWV-A infection, we set out to explore whether this effect is due to DWV-A mimicking naturally occurring changes in brain gene expression that are associated with behavioral maturation. Consistent with this hypothesis, we found that brain gene expression profiles of DWV-A infected bees resembled those of foragers, even in individuals that were much younger than typical foragers. In addition, brain transcriptional regulatory network analysis revealed a positive association between DWV-A infection and transcription factors previously associated with honey bee foraging behavior. Surprisingly, single-cell RNA-Sequencing implicated glia, not neurons, in this effect; there are relatively few glial cells in the insect brain and they are rarely associated with behavioral plasticity. Covert DWV-A infection also has been linked to impaired learning, which together with precocious foraging can lead to increased occurrence of infected bees from one colony mistakenly entering another colony, especially under crowded modern apiary conditions. These findings provide new insights into the mechanisms by which DWV-A affects honey bee health and colony survival.

Kinome Analysis of Honeybee (Apis mellifera L.) Dark-Eyed Pupae Identifies Biomarkers and Mechanisms of Tolerance to Varroa Mite Infestation

The mite Varroa destructor is a serious threat to honeybee populations. Selective breeding for Varroa mite tolerance could be accelerated by biomarkers within individual bees that could be applied to evaluate a colony phenotype. Previously, we demonstrated differences in kinase-mediated signaling between bees from colonies of extreme phenotypes of mite susceptibility. We expand these findings by defining a panel of 19 phosphorylation events that differ significantly between individual pupae from multiple colonies with distinct Varroa mite tolerant phenotypes. The predictive capacity of these biomarkers was evaluated by analyzing uninfested pupae from eight colonies representing a spectrum of mite tolerance. The pool of biomarkers effectively discriminated individual pupae on the basis of colony susceptibility to mite infestation. Kinome analysis of uninfested pupae from mite tolerant colonies highlighted an increased innate immune response capacity. The implication that differences in innate immunity contribute to mite susceptibility is supported by the observation that induction of innate immune signaling responses to infestation is compromised in pupae of the susceptible colonies. Collectively, biomarkers within individual pupae that are predictive of the susceptibility of colonies to mite infestation could provide a molecular tool for selective breeding of tolerant colonies.

Engineered symbionts activate honey bee immunity and limit pathogens

Honey bees are essential pollinators threatened by colony losses linked to the spread of parasites and pathogens. Here, we report a new approach for manipulating bee gene expression and protecting bee health. We engineered a symbiotic bee gut bacterium, Snodgrassella alvi, to induce eukaryotic RNA interference (RNAi) immune responses. We show that engineered S. alvi can stably recolonize bees and produce double-stranded RNA to activate RNAi and repress host gene expression, thereby altering bee physiology, behavior, and growth. We used this approach to improve bee survival after a viral challenge, and we show that engineered S. alvi can kill parasitic Varroa mites by triggering the mite RNAi response. This symbiont-mediated RNAi approach is a tool for studying bee functional genomics and potentially for safeguarding bee health.