Parallel epigenomic and transcriptomic responses to viral infection in honey bees (Apis mellifera)

The role of host plants in driving pathogen susceptibility in insects through chemicals, immune responses and microbiota

In this comprehensive exploration, we delve into the pivotal role of host plants in shaping the intricate interactions between herbivorous insects and their pathogens. Recent decades have seen a surge in studies that demonstrate that host plants are crucial drivers of the interactions between insects and pathogens, providing novel insights into the direct and indirect interactions that shape tri-trophic interactions. These studies have built on a wide range of pathogens, from viruses to bacteria, and from protozoans to fungi. We summarise these studies, and discuss the mechanisms of plant-mediated insect resistance to infection, ranging from the toxicity of plant chemicals to pathogens to enhancement of anti-pathogen immune responses, and modulation of the insect's microbiome. Although we provide evidence for the roles of all these mechanisms, we also point out that the majority of existing studies are phenomenological, describing patterns without addressing the underlying mechanisms. To further our understanding of these tri-trophic interactions, we therefore urge researchers to design their studies to enable them specifically to distinguish the mechanisms by which plants affect insect susceptibility to pathogens.

Molecular antiviral responses, immune priming and inheritance in insects

Viral diseases transmitted by insects to plants cause severe agricultural damage and arboviruses transmitted to humans cause severe disease outbreaks. The interaction between viruses and the insect defences is complex and has evolved into acting-counteracting molecular interplays. Viruses depict complex molecular mechanisms to ensure invasion, replication and exit the insect host cell, to invade other cells. On the other hand, insect cells use molecular strategies to recognize, halt replication and eliminate the invaders. In turn, virus counteract with evasive strategies. The main antiviral defence mechanism RNA interference (RNAi) recognizes and degrades viral RNA, thereby inhibiting viral replication. These in conjunction with other canonical immune pathways, Toll, IMD, JAK/STAT and Akt-ERK developed mainly to combat bacteria, fungi and protozoa, along with mechanisms to eliminate infected cells like apoptosis and phagocytosis comprise a multifactorial system. Insects exposed to an attenuated or sublethal viral infection could respond with faster and enhanced immune responses to the same pathogen (priming), which is like immunological memory in vertebrates. Several mechanisms have been proposed to explain priming, including endoreplication, epigenetic gene modifications by DNA methylation and histone acetylation. Priming could be inherited by the offspring (transgenerational immune priming, TGIP). However, the precise molecular mechanisms underlying TGIP remain to be elucidated. This article reviews the molecular mechanisms employed by insects to combat viral infections, discusses the current information and the outstanding research questions in the area.

Gut-binding peptides as potential tools to reduce virus binding to honey bee gut surface proteins

Colonies of the western honey bee, Apis mellifera, are severely impacted by a wide range of stressors, with Varroa mites and associated viruses being among the most serious threats to honey bee health. Viral load plays an important role in colony demise, with the iflavirus Deformed wing virus (DWV) and the dicistrovirus Israeli acute paralysis virus (IAPV) being of particular concern. By feeding adult honey bees on a phage display library to identify gut-binding peptides (R. Mishra, Y. Guo, P. Kumar, P. E. Cantón, C. S. Tavares, R. Banerjee, S. Kuwar, and B. C. Bonning, Curr Res Insect Sci, 1:100012, https://doi.org/10.1016/j.cris.2021.100012), we identified Bee midgut-Binding Peptide (BBP2.1), which shares 75% and 85% identity with regions on the DWV capsid protein and IAPV ORFx protein, respectively. These viral protein domains are likely to be instrumental in virus interaction with the honey bee gut. Pull-down assays with honey bee gut brush border membrane vesicles were used to confirm peptide-mCherry binding to the gut for BBP2.1 and the two similar virus-derived sequences, peptides BBP2.1DWV and BBP2.1IAPV. In vitro competition assays showed that all three peptides compete with both IAPV and DWV virions for binding to honey bee gut-derived brush border membrane vesicles, suggesting that the three peptides and the two viruses bind to the same proteins. Ingestion of BBP2.1 reduced the movement of DWV, but not IAPV from the honey bee gut into the body and did not rescue IAPV-associated mortality. These results are discussed in relation to the biological function of IAPV ORFx and the potential utility of virus-blocking peptides for suppression of virus infection to reduce virus load and virus-associated honey bee mortality.IMPORTANCEEach year, approximately 40% of managed honey bee hives in the United States are lost due to a variety of environmental stressors. Although increases in virus infection are among the most important factors resulting in colony loss, there are currently no effective tools for the management of virus infection in honey bees. In this study, we identified a peptide that binds to the gut of the honey bee and competes with two of the most important honey bee viruses, Israeli acute paralysis virus of bees (IAPV) and Deformed wing virus (DWV), for binding to gut proteins. In vivo competition between this peptide and DWV demonstrates the potential utility of gut-binding peptides for the protection of honey bees from virus infection for reduced virus-associated honey bee mortality.

An RPA-CRISPR/Cas12a based platform for rapid, sensitive, and visual detection of Apis mellifera filamentous virus

Apis mellifera filamentous virus (AmFV) is an emerging DNA virus significantly affecting honey bee health. AmFV infections weaken bee resistance to other pathogens, and can cause tissue lysis and death. Early, accurate detection of AmFV is crucial for timely intervention and preventing large-scale outbreaks. Current AmFV detection relies largely on polymerase chain reaction (PCR)-based methods. To enable rapid field detection of AmFV, we developed a rapid and ultrasensitive detection platform using recombinase polymerase amplification (RPA) combined with clustered regularly interspaced short palindromic repeats (CRISPR) / CRISPR-associated nuclease 12a (Cas12a) technology. A CRISPR RNA (crRNA1) specifically targeting the AmFV Bro gene was designed, ensuring no cross-reactivity with other insect DNA viruses or uninfected honey bees. After optimization of the reaction time, the platform generated results within 35 min: 20 min for the RPA reaction and 15 min for CRISPR-mediated cleavage. Two visualization approaches, fluorescence-based and lateral flow dipstick, were used to display the detection results. The detection sensitivity of both approaches was as few as 10 copies of the AmFV genome. Validation with field-collected honey bee samples demonstrated consistency with conventional PCR, revealing widespread latent AmFV infections in the field. Taken together, we successfully developed an RPA-CRISPR/Cas12 platform for rapid, specific, and sensitive detection of AmFV in Apis mellifera and Apis cerana. This platform holds promise as a simple, accurate, and cost-effective tool for point-of-care AmFV diagnosis in the field.

Genetic Diversity in Candidate Single-Nucleotide Polymorphisms Associated with Resistance in Honeybees in the Czech Republic Using the Novel SNaPshot Genotyping Panel

Background/Objectives: The increasing pressure from pathogens and parasites on Apis mellifera populations is resulting in significant colony losses. It is desirable to identify resistance-associated single-nucleotide polymorphisms (SNPs) and their variability for the purpose of breeding resilient honeybee lines. This study examined the genetic diversity of 13 SNPs previously studied for associations with various resistance-providing traits, including six linked to Varroa-specific hygiene, five linked to suppressed mite reproduction, one linked to immune response, and one linked to chalkbrood resistance. Methods: Genotyping was performed using a novel SNaPshot genotyping panel designed for this study. The sample pool consisted of 308 honeybee samples in total, covering all 77 administrative districts of the Czech Republic. Results: All examined loci were polymorphic. The frequency of positive alleles in our population is medium to low, depending on the specific SNP. An analysis of genotype frequencies revealed that most loci exhibited the Hardy-Weinberg equilibrium. A comparison of the allele and genotype frequencies of the same locus between samples from hives and samples from flowers revealed no significant differences. The genetic diversity, as indicated by the heterozygosity values, ranged from 0.05 to 0.50. The fixation index (F) was, on average, close to zero, indicating minimal influence of inbreeding or non-random mating on the genetic structure of the analyzed samples. Conclusions: The obtained results provide further insights into the genetic variation of SNPs associated with the immune response and resistance to pathogens in honeybee populations in the Czech Republic. This research provides a valuable foundation for future studies of honeybee diversity and breeding.

Bees as environmental and toxicological bioindicators in the light of pesticide non-targeted exposure

Pesticides have a significant impact on the environment, harming valuable non-target organisms like bees. Honeybees, in particular, are ideal bioindicators of pesticide exposure due to extensive research on how pesticides affect their behavior, immunity, development, biomolecules, and detoxification. However, wild pollinators are less studied in terms of pesticide exposure, and their inclusion is essential for a comprehensive risk assessment. Additionally, food chain organisms, such as the Asian hornet, could serve as indicators of pesticide bioaccumulation. Addressing gaps in honeybee toxicology, understanding the limitations, and exploring the role of wild pollinators and insects as complementary indicators, along with advancements in risk assessment methodologies, could enhance predictive models. These models would help anticipate environmental pesticide impacts while reducing the need for costly, time-consuming research.

Epigenetic Regulation in Insect-Microbe Interactions

Insects have evolved diverse interactions with a variety of microbes, such as pathogenic fungi, bacteria, and viruses. The immune responses of insect hosts, along with the dynamic infection process of microbes in response to the changing host environment and defenses, require rapid and fine-tuned regulation of gene expression programs. Epigenetic mechanisms, including DNA methylation, histone modifications, and noncoding RNA regulation, play important roles in regulating the expression of genes involved in insect immunity and microbial pathogenicity. This review highlights recent discoveries and insights into epigenetic regulatory mechanisms that modulate insect-microbe interactions. A deeper understanding of these regulatory mechanisms underlying insect-microbe interactions holds promise for the development of novel strategies for biological control of insect pests and mitigation of vector-borne diseases.

Gut microbiota contribute to variations in honey bee foraging intensity

Gut microbiomes are increasingly recognized for mediating diverse biological aspects of their hosts, including complex behavioral phenotypes. Although many studies have reported that experimental disruptions to the gut microbial community result in atypical host behavior, studies that address how gut microbes contribute to adaptive behavioral trait variation are rare. Eusocial insects represent a powerful model to test this, because of their simple gut microbiota and complex division of labor characterized by colony-level variation in behavioral phenotypes. Although previous studies report correlational differences in gut microbial community associated with division of labor, here, we provide evidence that gut microbes play a causal role in defining differences in foraging behavior between European honey bees (Apis mellifera). We found that gut microbial community structure differed between hive-based nurse bees and bees that leave the hive to forage for floral resources. These differences were associated with variation in the abundance of individual microbes, including Bifidobacterium asteroides, Bombilactobacillus mellis, and Lactobacillus melliventris. Manipulations of colony demography and individual foraging experience suggested that differences in gut microbial community composition were associated with task experience. Moreover, single-microbe inoculations with B. asteroides, B. mellis, and L. melliventris caused effects on foraging intensity. These results demonstrate that gut microbes contribute to division of labor in a social insect, and support a role of gut microbes in modulating host behavioral trait variation.

Quantitative trait loci mapping for survival of virus infection and virus levels in honey bees

Israeli acute paralysis virus (IAPV) is a highly virulent, Varroa-vectored virus that is of global concern for honey bee health. Little is known about the genetic basis of honey bees to withstand infection with IAPV or other viruses. We set up and analyzed a backcross between preselected honey bee colonies of low and high IAPV susceptibility to identify quantitative trait loci (QTL) associated with IAPV susceptibility. Experimentally inoculated adult worker bees were surveyed for survival and selectively sampled for QTL analysis based on SNPs identified by whole-genome resequencing and composite interval mapping. Additionally, natural titers of other viruses were quantified in the abdomen of these workers via qPCR and also used for QTL mapping. In addition to the full dataset, we analyzed distinct subpopulations of susceptible and non-susceptible workers separately. These subpopulations are distinguished by a single, suggestive QTL on chromosome 6, but we identified numerous other QTL for different abdominal virus titers, particularly in the subpopulation that was not susceptible to IAPV. The pronounced QTL differences between the susceptible and non-susceptible subpopulations indicate either an interaction between IAPV infection and the bees' interaction with other viruses or heterogeneity among workers of a single cohort that manifests itself as IAPV susceptibility and results in distinct subgroups that differ in their interaction with other viruses. Furthermore, our results indicate that low susceptibility of honey bees to viruses can be caused by both, virus tolerance and virus resistance. QTL were partially overlapping among different viruses, indicating a mixture of shared and specific processes that control viruses. Some functional candidate genes are located in the QTL intervals, but their genomic co-localization with numerous genes of unknown function delegates any definite characterization of the underlying molecular mechanisms to future studies.

Jute (Corchorus olitorius L.) Nanocrystalline Cellulose Inhibits Insect Virus via Gut Microbiota and Metabolism

Natural plant nanocrystalline cellulose (NCC), exhibiting a number of exceptional performance characteristics, is widely used in food fields. However, little is known about the relationship between NCC and the antiviral effect in animals. Here, we tested the function of NCC in antiviral methods utilizing honey bees as the model organism employing Israeli acute paralysis virus (IAPV), a typical RNA virus of honey bees. In both the lab and the field, we fed the IAPV-infected bees various doses of jute NCC (JNCC) under carefully controlled conditions. We found that JNCC can reduce IAPV proliferation and improve gut health. The metagenome profiling suggested that IAPV infection significantly decreased the abundance of gut core bacteria, while JNCC therapy considerably increased the abundance of the gut core bacteria Snodgrassella alvi and Lactobacillus Firm-4. Subsequent metabolome analysis further revealed that JNCC promoted the biosynthesis of fatty acids and unsaturated fatty acids, accelerated the purine metabolism, and then increased the expression of antimicrobial peptides (AMPs) and the genes involved in the Wnt and apoptosis signaling pathways against IAPV infection. Our results highlighted that JNCC could be considered as a prospective candidate agent against a viral infection.

Influence of social lifestyles on host-microbe symbioses in the bees

Microbiomes are increasingly recognised as critical for the health of an organism. In eusocial insect societies, frequent social interactions allow for high-fidelity transmission of microbes across generations, leading to closer host-microbe coevolution. The microbial communities of bees with other social lifestyles are less studied, and few comparisons have been made between taxa that vary in social structure. To address this gap, we leveraged a cloud-computing resource and publicly available transcriptomic data to conduct a survey of microbial diversity in bee samples from a variety of social lifestyles and taxa. We consistently recover the core microbes of well-studied corbiculate bees, supporting this method's ability to accurately characterise microbial communities. We find that the bacterial communities of bees are influenced by host location, phylogeny and social lifestyle, although no clear effect was found for fungal or viral microbial communities. Bee genera with more complex societies tend to harbour more diverse microbes, with Wolbachia detected more commonly in solitary tribes. We present a description of the microbiota of Euglossine bees and find that they do not share the "corbiculate core" microbiome. Notably, we find that bacteria with known anti-pathogenic properties are present across social bee genera, suggesting that symbioses that enhance host immunity are important with higher sociality. Our approach provides an inexpensive means of exploring microbiomes of a given taxa and identifying avenues for further research. These findings contribute to our understanding of the relationships between bees and their associated microbial communities, highlighting the importance of considering microbiome dynamics in investigations of bee health.

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.

Consequences of microsporidian prior exposure for virus infection outcomes and bumble bee host health

Host-parasite interactions do not occur in a vacuum, but in connected multi-parasite networks that can result in co-exposures and coinfections of individual hosts. These can affect host health and disease ecology, including disease outbreaks. However, many host-parasite studies examine pairwise interactions, meaning we still lack a general understanding of the influence of co-exposures and coinfections. Using the bumble bee Bombus impatiens, we study the effects of larval exposure to a microsporidian Nosema bombi, implicated in bumble bee declines, and adult exposure to Israeli Acute Paralysis Virus (IAPV), an emerging infectious disease from honey bee parasite spillover. We hypothesize that infection outcomes will be modified by co-exposure or coinfection. Nosema bombi is a potentially severe, larval-infecting parasite, and we predict that prior exposure will result in decreased host resistance to adult IAPV infection. We predict double parasite exposure will also reduce host tolerance of infection, as measured by host survival. Although our larval Nosema exposure mostly did not result in viable infections, it partially reduced resistance to adult IAPV infection. Nosema exposure also negatively affected survival, potentially due to a cost of immunity in resisting the exposure. There was a significant negative effect of IAPV exposure on survivorship, but prior Nosema exposure did not alter this survival outcome, suggesting increased tolerance given the higher IAPV infections in the bees previously exposed to Nosema. These results again demonstrate that infection outcomes can be non-independent when multiple parasites are present, even when exposure to one parasite does not result in a substantial infection.

Transcriptomic Responses Underlying the High Virulence of Black Queen Cell Virus and Sacbrood Virus following a Change in Their Mode of Transmission in Honey Bees (Apis mellifera)

BACKGROUND

Over the last two decades, honey bees (Apis mellifera) have suffered high rates of colony losses that have been attributed to a variety of factors, chief among which are viral pathogens, such as deformed wing virus (DWV), whose virulence has increased because of vector-based transmission by the invasive, ectoparasitic varroa mite (Varroa destructor). A shift in the experimental mode of transmission of the black queen cell virus (BQCV) and sacbrood virus (SBV) from fecal/food-oral (direct horizontal) to vector-mediated (indirect horizontal) transmission also results in high virulence and viral titers in pupal and adult honey bees. Agricultural pesticides represent another factor that acts independently or in interaction with pathogens, and they are also thought to cause colony loss. Understanding the molecular mechanisms underlying the higher virulence following a vector-based mode of transmission provides deeper insight into honey bee colony losses, as does determining whether or not host-pathogen interactions are modulated by exposure to pesticides.

METHODS

Through an experimental design with controlled laboratory, we investigated the effects of the modes of transmission of BQCV and SBV (feeding vs. vector-mediated via injection) alone or in combination with chronic exposure to sublethal and field-realistic concentrations of flupyradifurone (FPF), a novel agricultural insecticide, on honey bee survival and transcription responses by using high-throughput RNA sequencing (RNA-seq) analysis.

RESULTS

Co-exposure to viruses via feeding (VF) or injection (VI) and FPF insecticide had no statistically significant interactive effect on their survival compared to, respectively, VF or VI treatments alone. Transcriptomic analysis revealed a distinct difference in the gene expression profiles of bees inoculated with viruses via injection (VI) and exposed to FPF insecticide (VI+FPF). The number of differentially expressed genes (DEGs) at log2 (fold-change) > 2.0 in VI bees (136 genes) or/and VI+FPF insecticide (282 genes) was very high compared to that of VF bees (8 genes) or the VF+FPF insecticide treatment (15 genes). Of these DEGs, the expression in VI and VI+FPF bees of some immune-related genes, such as those for antimicrobial peptides, Ago2, and Dicer, was induced. In short, several genes encoding odorant binding proteins, chemosensory proteins, odor receptors, honey bee venom peptides, and vitellogenin were downregulated in VI and VI+FPF bees.

CONCLUSIONS

Given the importance of these suppressed genes in honey bees' innate immunity, eicosanoid biosynthesis, and olfactory associative function, their inhibition because of the change in the mode of infection with BQCV and SBV to vector-mediated transmission (injection into haemocoel) could explain the high virulence observed in these viruses when they were experimentally injected into hosts. These changes may help explain why other viruses, such as DWV, represent such a threat to colony survival when transmitted by varroa mites.

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.

Expression of Molecular Markers of Resilience against Varroa destructor and Bee Viruses in Ethiopian Honey Bees (Apis mellifera simensis) Focussing on Olfactory Sensing and the RNA Interference Machinery

Varroa destructor mites and the viruses it vectors are two major factors leading to high losses of honey bees (Apis mellifera) colonies worldwide. However, honey bees in some African countries show resilience to varroa infestation and/or virus infections, although little is known about the mechanisms underlying this resilience. In this study, we investigated the expression profiles of some key molecular markers involved in olfactory sensing and RNA interference, as these processes may contribute to the bees' resilience to varroa infestation and virus infection, respectively. We found significantly higher gene expression of the odorant binding protein, OBP14, in the antennae of Ethiopian bees compared to Belgian bees. This result suggests the potential of OBP14 as a molecular marker of resilience to mite infestation. Scanning electron microscopy showed no significant differences in the antennal sensilla occurrence and distribution, suggesting that resilience arises from molecular processes rather than morphological adaptations. In addition, seven RNAi genes were upregulated in the Ethiopian honey bees and three of them-Dicer-Drosha, Argonaute 2, and TRBP2-were positively correlated with the viral load. We can conclude that the antiviral immune response was triggered when bees were experiencing severe viral infection and that this might contribute to the bees' resilience to viruses.

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.

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.

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.

Antibiotics-induced changes in intestinal bacteria result in the sensitivity of honey bee to virus

Antibiotics are omnipresent in the environment due to their widespread use, and they have wide-ranging negative impacts on organisms. Virus resistance differs substantially between domesticated Apis mellifera and wild Apis cerana, although both are commonly raised in China. Here, we investigated whether antibiotics can increase the sensitivity of honey bees to viral infection using the Israeli acute paralysis virus (IAPV) and tetracycline as representative virus and antibiotic. Although IAPV multiplied to lower levels in A. cerana than A. mellifera, resulting in decreased mortality (P < 0.01), there was no significant difference in immune responses to viral infection between the two species. Adult worker bees (A. cerana and A. mellifera) were treated with or without tetracycline to demonstrate the prominent role of gut microbiota against viral infection, and found Lactobacillus played a vital antiviral role in A. cerana. In A. cerana but not A. mellifera, tetracycline treatment reduced clearly bee survival and increased susceptibility to IAPV infection (P < 0.01). Our findings revealed that long-term antibiotic treatment in A. mellifera had altered the native gut microbiome and promoted the sensitivity to viral infection. We highlight the effects of antibiotics exposure on resistance to microbial and viral infection.

Current Knowledge on Bee Innate Immunity Based on Genomics and Transcriptomics

As important pollinators, bees play a critical role in maintaining the balance of the ecosystem and improving the yield and quality of crops. However, in recent years, the bee population has significantly declined due to various pathogens and environmental stressors including viruses, bacteria, parasites, and increased pesticide application. The above threats trigger or suppress the innate immunity of bees, their only immune defense system, which is essential to maintaining individual health and that of the colony. In addition, bees can be divided into solitary and eusocial bees based on their life traits, and eusocial bees possess special social immunities, such as grooming behavior, which cooperate with innate immunity to maintain the health of the colony. The omics approach gives us an opportunity to recognize the distinctive innate immunity of bees. In this regard, we summarize innate bee immunity from a genomic and transcriptomic perspective. The genetic characteristics of innate immunity were revealed by the multiple genomes of bees with different kinds of sociality, including honeybees, bumblebees, wasps, leaf-cutter bees, and so on. Further substantial transcriptomic data of different tissues from diverse bees directly present the activation or suppression of immune genes under the infestation of pathogens or toxicity of pesticides.

Acaricide flumethrin-induced sublethal risks in honeybees are associated with gut symbiotic bacterium Gilliamella apicola through microbe-host metabolic interactions

Flumethrin is one of the few acaricides that permit the control of Varroa disease or varroosis in bee colonies. However, flumethrin accumulates in hive products. We previously discovered that sublethal doses of flumethrin induce significant physiological stress in honeybees (Apis mellifera L.), however its potential impacts on the honeybee gut microenvironment remains unknown. To fill this gap, honeybees were exposed to a field-relevant concentration of flumethrin (10 μg/L) for 14 d and its potential impacts on gut system were evaluated. The results indicated that flumethrin triggered immune responses in the gut but had limited effects on survival and gut microbial composition. However, survival stress drastically increased in bees exposed to antibiotics, suggesting that the gut microbiota is closely related to flumethrin-induced dysbiosis in the bee gut. Based on a non-targeted metabolomics approach, flumethrin at 10 μg/L considerably altered the composition of intestinal metabolites, and we discovered that this metabolic stress was closely linked with a reduction of gut core bacterial endosymbiont Gilliamella spp. through a combination of microbiological and metabolomics investigations. Finally, an in vitro study showed that while flumethrin does not directly inhibit the growth of Gilliamella apicola isolates, it does have a significant impact on the glycerophospholipid metabolism in bacteria cells, which was also observed in host bees. These findings indicated that even though flumethrin administered at environmental relevant concentrations does not significantly induce death in honeybees, it still alters the metabolism balance between honeybees and the gut symbiotic bacterium, G. apicola. The considerable negative impact of flumethrin on the honeybee gut microenvironment emphasizes the importance of properly monitoring acaricide to avoid potential environmental concerns, and further studies are needed to illustrate the mode of action of bee health-gut microbiota-exogenous pesticides.

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.

Comparative studies of genomic and epigenetic factors influencing transcriptional variation in two insect species

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.

Colony-Level Viral Load Influences Collective Foraging in Honey Bees

Nutrition is an important component of social insect colony health especially in the face of stressors such as parasitism and viral infections. Honey bees are known to preferentially select nectar and pollen based on macronutrient and phytochemical contents and in response to pathogen loads. However, given that honey bees live in colonies, collective foraging decisions may be impacted directly by forager infection status but also by colony health. This field experiment was conducted to determine if honey bee viral infections are correlated with pollen and nectar foraging and if these associations are impacted more by colony or forager infection. By comparing regressions with and without forager and colony variables and through structural equation models, we were able to determine the relative contributions of colony and forager virus loads on forager decisions. We found that foragers had higher numbers and levels of BQCV and CBPV but lower levels of DWV viruses than their respective colonies. Overall, individuals appeared to forage based a combination of their own and colony health but with greater weight given to colony metrics. Colony parasitism by Varroa mites, positively correlated with both forager and colony DWV-B levels, was negatively associated with nectar weight. Further, colony DWV-B levels were negatively associated with individually foraged pollen protein: lipid ratios but positively correlated with nectar weight and sugar content. This study shows that both colony and forager health can simultaneously mediate individual foraging decisions and that the importance of viral infections and parasite levels varies with foraging metrics. Overall, this work highlights the continued need to explore the interactions of disease, nutrition, and genetics in social interactions and structures.

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.

Current understanding and perspectives on the potential mechanisms of immune priming in beetles

Beetles are the most diverse group of insects in Insecta which can be found in almost every habitat and environment on Earth. The possessing of the rapid and effective immune defenses is one of the important factors for their success. It is generally recognized that beetles only rely on the non-specific innate immune defense, without immunological memory, to fight against pathogens. However, there was cumulative evidence for the innate immune memory in invertebrates, including beetles, over the last decades, implying that insect innate immunity is more complex and has more features than previously thought. In beetles, it has been well documented that the specific or nonspecific enhanced immunocompetence can persist throughout development within generations and can even be transferred to the descendents in the next generation. Although insect immune priming might be shaped by epigenetic modifications and transferring effectors, mRNA and microbial signals, the solid experimental evidence to support the causal relationship between any of them and immune priming is still scarce. The combined usage of 'omics' approaches and CRISPR/Cas9 in the appropriate insect models with well-known genetic background, Tribolium castaneum and Tenebrio molitor, will help us to decipher the molecular mechanisms by which immune priming occurs in beetles in depth.

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.

Tools and applications for integrative analysis of DNA methylation in social insects

DNA methylation is a common epigenetic signalling tool and an important biological process which is widely studied in a large array of species. The presence, level and function of DNA methylation vary greatly across species. In some insects, DNA methylation systems are minimal, and overall methylation rates tend to be low in all studied insect species. Low methylation levels probed by whole-genome bisulphite sequencing require great care with respect to data quality control and interpretation. Here, we introduce BWASP/R, a complete workflow that allows efficient, scalable and entirely reproducible analyses of raw DNA methylation sequencing data. Consistent application of quality control filters and analysis parameters provides fair comparisons among different studies and an integrated view of all experiments on one species. We describe the capabilities of the BWASP/R workflow by re-analysing several publicly available social insect WGBS data sets, comprising 70 samples and cumulatively 147 replicates from four different species. We show that the CpG methylome comprises only about 1.5% of CpG sites in the honeybee genome and that the cumulative data are consistent with genetic signatures of site accessibility and physiological control of methylation levels.

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.

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.

Multi-tiered analyses of honey bees that resist or succumb to parasitic mites and viruses

Background

Varroa destructor mites, and the numerous viruses they vector to their honey bee hosts, are among the most serious threats to honey bee populations, causing mortality and morbidity to both the individual honey bee and colony, the negative effects of which convey to the pollination services provided by honey bees worldwide. Here we use a combination of targeted assays and deep RNA sequencing to determine host and microbial changes in resistant and susceptible honey bee lineages. We focus on three study sets. The first involves field sampling of sympatric western bees, some derived from resistant stock and some from stock susceptible to mites. The second experiment contrasts three colonies more deeply, two from susceptible stock from the southeastern U.S. and one from mite-resistant bee stock from Eastern Texas. Finally, to decouple the effects of mites from those of the viruses they vector, we experimentally expose honey bees to DWV in the laboratory, measuring viral growth and host responses.

Results

We find strong differences between resistant and susceptible bees in terms of both viral loads and bee gene expression. Interestingly, lineages of bees with naturally low levels of the mite-vectored Deformed wing virus, also carried lower levels of viruses not vectored by mites. By mapping gene expression results against current ontologies and other studies, we describe the impacts of mite parasitism, as well as viruses on bee health against two genetic backgrounds. We identify numerous genes and processes seen in other studies of stress and disease in honey bee colonies, alongside novel genes and new patterns of expression.

Conclusions

We provide evidence that honey bees surviving in the face of parasitic mites do so through their abilities to resist the presence of devastating viruses vectored by these mites. In all cases, the most divergence between stocks was seen when bees were exposed to live mites or viruses, suggesting that gene activation, rather than constitutive expression, is key for these interactions. By revealing responses to viral infection and mite parasitism in different lineages, our data identify candidate proteins for the evolution of mite tolerance and virus resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08032-z.

Honey bee viruses are highly prevalent but at low intensities in wild pollinators of cucurbit agroecosystems

Managed and wild bee populations are in decline around the globe due to several biotic and abiotic stressors. Pathogenic viruses associated with the Western honey bee (Apis mellifera) have been identified as key contributors to losses of managed honey bee colonies, and are known to be transmitted to wild bee populations through shared floral resources. However, little is known about the prevalence and intensity of these viruses in wild bee populations, or how bee visitation to flowers impacts viral transmission in agroecosystems. This study surveyed honey bee, bumble bee (Bombus impatiens) and wild squash bee (Eucera (Peponapis) pruinosa) populations in Cucurbita agroecosystems across Pennsylvania (USA) for the prevalence and intensity of five honey bee viruses: acute bee paralysis virus (ABPV), deformed wing virus (DWV), Israeli acute paralysis virus (IAPV), Kashmir bee virus (KBV), and slow bee paralysis virus (SBPV). We investigated the potential role of bee visitation rate to flowers on DWV intensity among species in the pollinator community, with the expectation that increased bee visitation to flowers would increase the opportunity for transmission events between host species. We found that honey bee viruses are highly prevalent but in lower titers in wild E. pruinosa and B. impatiens than in A. mellifera populations throughout Pennsylvania (USA). DWV was detected in 88% of B. impatiens, 48% of E. pruinosa, and 95% of A. mellifera. IAPV was detected in 5% of B. impatiens and 4% of E. pruinosa, compared to 9% in A. mellifera. KBV was detected in 1% of B. impatiens and 5% of E. pruinosa, compared to 32% in A. mellifera. Our results indicate that DWV titers are not correlated with bee visitation in Cucurbita fields. The potential fitness impacts of these low viral titers detected in E. pruinosa remain to be investigated.

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.

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.

Lineage and Parent-of-Origin Effects in DNA Methylation of Honey Bees (Apis mellifera) Revealed by Reciprocal Crosses and Whole-Genome Bisulfite Sequencing

Parent-of-origin methylation arises when the methylation patterns of a particular allele are dependent on the parent it was inherited from. Previous work in honey bees has shown evidence of parent-of-origin-specific expression, yet the mechanisms regulating such pattern remain unknown in honey bees. In mammals and plants, DNA methylation is known to regulate parent-of-origin effects such as genomic imprinting. Here, we utilize genotyping of reciprocal European and Africanized honey bee crosses to study genome-wide allele-specific methylation patterns in sterile and reproductive individuals. Our data confirm the presence of allele-specific methylation in honey bees in lineage-specific contexts but also importantly, though to a lesser degree, parent-of-origin contexts. We show that the majority of allele-specific methylation occurs due to lineage rather than parent-of-origin factors, regardless of the reproductive state. Interestingly, genes affected by allele-specific DNA methylation often exhibit both lineage and parent-of-origin effects, indicating that they are particularly labile in terms of DNA methylation patterns. Additionally, we re-analyzed our previous study on parent-of-origin-specific expression in honey bees and found little association with parent-of-origin-specific methylation. These results indicate strong genetic background effects on allelic DNA methylation and suggest that although parent-of-origin effects are manifested in both DNA methylation and gene expression, they are not directly associated with each other.

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.

Reproductive activation in honeybee (Apis mellifera) workers protects against abiotic and biotic stress

Social insect reproductives exhibit exceptional longevity instead of the classic trade-off between somatic maintenance and reproduction. Even normally sterile workers experience a significant increase in life expectancy when they assume a reproductive role. The mechanisms that enable the positive relation between the antagonistic demands of reproduction and somatic maintenance are unclear. To isolate the effect of reproductive activation, honeybee workers were induced to activate their ovaries. These reproductively activated workers were compared to controls for survival and gene expression patterns after exposure to Israeli Acute Paralysis Virus or the oxidative stressor paraquat. Reproductive activation increased survival, indicating better immunity and oxidative stress resistance. After qPCR analysis confirmed our experimental treatments at the physiological level, whole transcriptome analysis revealed that paraquat treatment significantly changed the expression of 1277 genes in the control workers but only two genes in reproductively activated workers, indicating that reproductive activation preemptively protects against oxidative stress. Significant overlap between genes that were upregulated by reproductive activation and in response to paraquat included prominent members of signalling pathways and anti-oxidants known to affect ageing. Thus, while our results confirm a central role of vitellogenin, they also point to other mechanisms to explain the molecular basis of the lack of a cost of reproduction and the exceptional longevity of social insect reproductives. Thus, socially induced reproductive activation preemptively protects honeybee workers against stressors, explaining their longevity. This article is part of the theme issue 'Ageing and sociality: why, when and how does sociality change ageing patterns?'

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.

Tissue-specific transcription patterns support the kinship theory of intragenomic conflict in honey bees (Apis mellifera)

Kin selection may act differently on genes inherited from parents (matrigenes and patrigenes), resulting in intragenomic conflict. This conflict can be observed as differential expression of matrigenes and patrigenes, or parent-specific gene expression (PSGE). In honey bees (Apis mellifera), intragenomic conflict is hypothesized to occur in multiple social contexts. Previously, we found that patrigene-biased expression in reproductive tissues was associated with increased reproductive potential in worker honey bees, consistent with the prediction that patrigenes are selected to promote selfish behaviour in this context. Here, we examined brain gene expression patterns to determine if PSGE is also found in other tissues. As before, the number of transcripts showing patrigene expression bias was significantly greater in the brains of reproductive vs. sterile workers, while the number of matrigene-biased transcripts was not significantly different. Twelve transcripts out of the 374 showing PSGE in either tissue showed PSGE in both brain and reproductive tissues; this overlap was significantly greater than expected by chance. However, the majority of transcripts show PSGE only in one tissue, suggesting the epigenetic mechanisms mediating PSGE exhibit plasticity between tissues. There was no significant overlap between transcripts that showed PSGE and transcripts that were significantly differentially expressed. Weighted gene correlation network analysis identified modules which were significantly enriched in both types of transcripts, suggesting that these genes may influence each other through gene networks. Our results provide further support for the kin selection theory of intragenomic conflict, and provide valuable insights into the mechanisms which may mediate this process.

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.

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.

The novel insecticides flupyradifurone and sulfoxaflor do not act synergistically with viral pathogens in reducing honey bee (Apis mellifera) survival but sulfoxaflor modulates host immunocompetence

The decline of insect pollinators threatens global food security. A major potential cause of decline is considered to be the interaction between environmental stressors, particularly between exposure to pesticides and pathogens. To explore pesticide-pathogen interactions in an important pollinator insect, the honey bee, we used two new nicotinic acetylcholine receptor agonist insecticides (nACHRs), flupyradifurone (FPF) and sulfoxaflor (SULF), at sublethal and field-realistic doses in a fully crossed experimental design with three common viral honey bee pathogens, Black queen cell virus (BQCV) and Deformed wing virus (DWV) genotypes A and B. Through laboratory experiments in which treatments were administered singly or in combination to individual insects, we recorded harmful effects of FPF and pathogens on honey bee survival and immune gene expression. Though we found no evidence of synergistic interactions among stressors on either honey bee survival or viral load, the combined treatment SULF and DWV-B led to a synergistic upregulation of dicer-like gene expression. We conclude that common viral pathogens pose a major threat to honey bees, while co-exposure to these novel nACHR insecticides does not significantly exacerbate viral impacts on host survival in the laboratory.

A salivary chitinase of Varroa destructor influences host immunity and mite's survival

Varroa destructor is an ectoparasite of honey bees and an active disease vector, which represents one of the most severe threats for the beekeeping industry. This parasitic mite feeds on the host’s body fluids through a wound in the cuticle, which allows food uptake by the mother mite and its progeny, offering a potential route of entrance for infecting microorganisms. Mite feeding is associated with saliva injection, whose role is still largely unknown. Here we try to fill this gap by identifying putative host regulation factors present in the saliva of V. destructor and performing a functional analysis for one of them, a chitinase (Vd-CHIsal) phylogenetically related to chitinases present in parasitic and predatory arthropods, which shows a specific and very high level of expression in the mite’s salivary glands. Vd-CHIsal is essential for effective mite feeding and survival, since it is apparently involved both in maintaining the feeding wound open and in preventing host infection by opportunistic pathogens. Our results show the important role in the modulation of mite-honey bee interactions exerted by a host regulation factor shared by different evolutionary lineages of parasitic arthropods. We predict that the functional characterization of Varroa sialome will provide new background knowledge on parasitism evolution in arthropods and the opportunity to develop new bioinspired strategies for mite control based on the disruption of their complex interactions with a living food source.

Varroa destructor is a parasitic mite of honey bees and a major driver of honey bee colony losses. The feeding mites inject a salivary blend of poorly known molecules, which regulate host physiology. Here, we have identified in silico putative host regulation factors occurring in Varroa saliva and characterized the functional role of a highly expressed chitinase, which is conserved across different evolutionary lineages of parasitic arthropods. This enzyme influences host immune response and mite’s survival. An in-depth functional analysis of Varroa saliva will shed light on parasitism evolution in arthropods and will pave the way towards the development of new bioinspired strategies for mite control.

bsAS, an antisense long non-coding RNA, essential for correct wing development through regulation of blistered/DSRF isoform usage

Natural Antisense Transcripts (NATs) are long non-coding RNAs (lncRNAs) that overlap coding genes in the opposite strand. NATs roles have been related to gene regulation through different mechanisms, including post-transcriptional RNA processing. With the aim to identify NATs with potential regulatory function during fly development, we generated RNA-Seq data in Drosophila developing tissues and found bsAS, one of the most highly expressed lncRNAs in the fly wing. bsAS is antisense to bs/DSRF, a gene involved in wing development and neural processes. bsAS plays a crucial role in the tissue specific regulation of the expression of the bs/DSRF isoforms. This regulation is essential for the correct determination of cell fate during Drosophila development, as bsAS knockouts show highly aberrant phenotypes. Regulation of bs isoform usage by bsAS is mediated by specific physical interactions between the promoters of these two genes, which suggests a regulatory mechanism involving the collision of RNA polymerases transcribing in opposite directions. Evolutionary analysis suggests that bsAS NAT emerged simultaneously to the long-short isoform structure of bs, preceding the emergence of wings in insects.

Identification of Immune Response to Sacbrood Virus Infection in Apis cerana Under Natural Condition

Chinese sacbrood virus (CSBV) is a serious threat to eastern honeybees (Apis cerana), especially larvae. However, the pathological mechanism of this deadly disease remains unclear. Here, we employed mRNA and small RNA (sRNA) transcriptome approach to investigate the microRNAs (miRNAs) and small interfering RNAs (siRNAs) expression changes of A. cerana larvae infected with CSBV under natural condition. We found that serine proteases involved in immune response were down-regulated, while the expression of siRNAs targeted to serine proteases were up-regulated. In addition, CSBV infection also affected the expression of larvae cuticle proteins such as larval cuticle proteins A1A and A3A, resulting in increased susceptibility to CSBV infection. Together, our results provide insights into sRNAs that they are likely to be involved in regulating honeybee immune response.

The gut microbiome defines social group membership in honey bee colonies

In the honey bee, genetically related colony members innately develop colony-specific cuticular hydrocarbon profiles, which serve as pheromonal nestmate recognition cues. Yet, despite high intracolony relatedness, the innate development of colony-specific chemical signatures by individual colony members is largely determined by the colony environment, rather than solely relying on genetic variants shared by nestmates. Therefore, it is puzzling how a nongenic factor could drive the innate development of a quantitative trait that is shared by members of the same colony. Here, we provide one solution to this conundrum by showing that nestmate recognition cues in honey bees are defined, at least in part, by shared characteristics of the gut microbiome across individual colony members. These results illustrate the importance of host-microbiome interactions as a source of variation in animal behavioral traits.

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.