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

1H NMR Profiling of Honey Bee Brains across Varying Ages and Seasons

Honey bees (Apis mellifera) provide a useful model for studying aging because of the differences in longevity between the relatively short-lived summer and long-lived winter bees, as well as bees lacking signs of cognitive senescence as they age. Bee brains were dissected from newly emerged, 14-day-, and 28-day- old bees in mid- and late summer, as well as brood nest bees in fall, winter, and spring, before, during, and after overwintering, respectively. Brains were examined with nuclear magnetic resonance (NMR) spectroscopy to analyze their metabolome. Nine variable importance in projection (VIP) variables were identified, primarily amino acids and choline derivatives. Differences in metabolite concentrations were found with different ages of summer bees, mostly between newly emerged and 14-day- old bees, such as a decrease in phenylalanine and an increase in β-alanine, but there were also changes in older adults, such as o-phosphocholine that declined in 28-day- old bees. Differences in brood nest bees were observed, including a decline in tryptophan and an increase in β-alanine. These may provide distinct metabolomic signatures with age and season. Such research holds promise for a better understanding of the complex interplays between bee physiology, development, and aging, which has implications for improving bee health and management.

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.

Learning performance and GABAergic pathway link to deformed wing virus in the mushroom bodies of naturally infected honey bees

Viral infections can be detrimental to the foraging ability of the western honey bee, Apis mellifera. The deformed wing virus (DWV) is the most common honey bee virus and has been proposed as a possible cause of learning and memory impairment. However, evidence for this phenomenon so far has come from artificially infected bees, while less is known about the implications of natural infections with the virus. Using the proboscis extension reflex (PER), we uncovered no significant association between a simple associative learning task and natural DWV load. However, when assessed through a reversal associative learning assay, bees with higher DWV load performed better in the reversal learning phase. DWV is able to replicate in the honey bee mushroom bodies, where the GABAergic signalling pathway has an antagonistic effect on associative learning but is crucial for reversal learning. Hence, we assessed the pattern of expression of several GABA-related genes in bees with different learning responses. Intriguingly, mushroom body expression of selected genes was positively correlated with DWV load, but only for bees with good reversal learning performance. We hypothesise that DWV might improve olfactory learning performance by enhancing the GABAergic inhibition of responses to unrewarded stimuli, which is consistent with the behavioural patterns that we observed. However, at higher disease burdens, which might be induced by an artificial infection or by a severe, natural Varroa infestation, this DWV-associated increase in GABA signalling could impair associative learning as previously reported by other studies.

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.

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.

Signatures of Adaptation, Constraints, and Potential Redundancy in the Canonical Immune Genes of a Key Pollinator

All organisms require an immune system to recognize, differentiate, and defend against pathogens. From an evolutionary perspective, immune systems evolve under strong selective pressures exerted by fast-evolving pathogens. However, the functional diversity of the immune system means that different immune components and their associated genes may evolve under varying forms of selection. Insect pollinators, which provide essential ecosystem services, are an important system in which to understand how selection has shaped immune gene evolution as their populations are experiencing declines with pathogens highlighted as a potential contributing factor. To improve our understanding of the genetic variation found in the immune genes of an essential pollinator, we performed whole-genome resequencing of wild-caught Bombus terrestris males. We first assessed nucleotide diversity and extended haplotype homozygosity for canonical immune genes finding the strongest signatures of positive selection acting on genes involved in pathogen recognition and antiviral defense, possibly driven by growing pathogen spread in wild populations. We also identified immune genes evolving under strong purifying selection, highlighting potential constraints on the bumblebee immune system. Lastly, we highlight the potential loss of function alleles present in the immune genes of wild-caught haploid males, suggesting that such genes are potentially less essential for development and survival and represent redundancy in the gene repertoire of the bumblebee immune system. Collectively, our analysis provides novel insights into the recent evolutionary history of the immune system of a key pollinator, highlighting targets of selection, constraints to adaptation, and potential redundancy.

Unraveling the Role of Antimicrobial Peptides in Insects

Antimicrobial peptides (AMPs) are short, mainly positively charged, amphipathic molecules. AMPs are important effectors of the immune response in insects with a broad spectrum of antibacterial, antifungal, and antiparasitic activity. In addition to these well-known roles, AMPs exhibit many other, often unobvious, functions in the host. They support insects in the elimination of viral infections. AMPs participate in the regulation of brain-controlled processes, e.g., sleep and non-associative learning. By influencing neuronal health, communication, and activity, they can affect the functioning of the insect nervous system. Expansion of the AMP repertoire and loss of their specificity is connected with the aging process and lifespan of insects. Moreover, AMPs take part in maintaining gut homeostasis, regulating the number of endosymbionts as well as reducing the number of foreign microbiota. In turn, the presence of AMPs in insect venom prevents the spread of infection in social insects, where the prey may be a source of pathogens.

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.

The adverse effects of synthetic acaricide tau-fluvalinate (tech.) on winter adult honey bees

Currently several pyrethroids (e.g., flumethrin and tau-fluvalinate) are used in apiculture worldwide as acaricides/miticides. The long half-lives of pyrethroids in synthetic acaricides applied to hive matrices, may adversely affect the health of bee colony. The potentially adverse effects of synthetic acaricide/miticide tau-fluvalinate (tech.) on winter honeybees were assessed in this study (OECD 245 2017). No dose-dependent mortality in in vitro reared winter honeybees was observed after chronic oral 10-day exposure to syrup (50% w/v) spiked with a maximum concentration of 750 μg a.i./kg diet and its 1/10 concentration. The No Observed Effect Concentration is ≥ 750 μg a.i./kg diet. Tau-fluvalinate testing for the sublethal effects on bee immune system showed up-regulated gene expression encoding abaecin, lysozyme, and defensin in both tested groups, however the expression of hymenoptaecin gene was reduced. Moreover, tau-fluvalinate significantly induced levels of DNA damage in exposed bees, which can result in adverse genotoxic effect.

An Insight Into the microRNA Profile of the Ectoparasitic Mite Varroa destructor (Acari: Varroidae), the Primary Vector of Honey Bee Deformed Wing Virus

The remarkably adaptive mite Varroa destructor is the most important honey bee ectoparasite. Varroa mites are competent vectors of deformed wing virus (DWV), and the Varroa-virus complex is a major determinant of annual honey bee colony mortality and collapse. MicroRNAs (miRNAs) are 22-24 nucleotide non-coding RNAs produced by all plants and animals and some viruses that influence biological processes through post-transcriptional regulation of gene expression. Knowledge of miRNAs and their function in mite biology remains limited. Here we constructed small RNA libraries from male and female V. destructor using Illumina's small RNA-Seq platform. A total of 101,913,208 and 91,904,732 small RNA reads (>18 nucleotides) from male and female mites were analyzed using the miRDeep2 algorithm. A conservative approach predicted 306 miRNAs, 18 of which were upregulated and 13 downregulated in female V. destructor compared with males. Quantitative real-time PCR validated the expression of selected differentially-expressed female Varroa miRNAs. This dataset provides a list of potential miRNA targets involved in regulating vital Varroa biological processes and paves the way for developing strategies to target Varroa and their viruses.

Methodology for Single Bee and Bee Brain 1H-NMR Metabolomics

The feasibility of metabolomic 1H NMR spectroscopy is demonstrated for its potential to help unravel the complex factors that are impacting honeybee health and behavior. Targeted and non-targeted 1H NMR metabolic profiles of liquid and tissue samples of organisms could provide information on the pathology of infections and on environmentally induced stresses. This work reports on establishing extraction methods for NMR metabolic characterization of Apis mellifera, the European honeybee, describes the currently assignable aqueous metabolome, and gives examples of diverse samples (brain, head, body, whole bee) and biologically meaningful metabolic variation (drone, forager, day old, deformed wing virus). Both high-field (600 MHz) and low-field (80 MHz) methods are applicable, and 1H NMR can observe a useful subset of the metabolome of single bees using accessible NMR instrumentation (600 MHz, inverse room temperature probe) in order to avoid pooling several bees. Metabolite levels and changes can be measured by NMR in the bee brain, where dysregulation of metabolic processes has been implicated in colony collapse. For a targeted study, the ability to recover 10-hydroxy-2-decenoic acid in mandibular glands is shown, as well as markers of interest in the bee brain such as GABA (4-aminobutyrate), proline, and arginine. The findings here support the growing use of 1H NMR more broadly in bees, native pollinators, and insects.

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.