The fungus Ascosphaera apis, an obligate fungal pathogen of honey bee brood, causes chalkbrood disease in honey bee larvae worldwide. Biological characteristics of the fungal pathogen and the molecular interactions between A. apis and honey bees have been studied extensively. However, little is known about the effects of A. apis infection on antioxidant enzyme activities and metabolic profiles of the gut of honey bee larvae. In this study, sandwich enzyme-linked immunosorbent assay and LC-MS based untargeted metabolomic analysis were employed to determine the changes in the specific activities of antioxidant enzymes and the metabolomic profiles in gut tissues of A. apis-infected larvae (10⁵ A. apis spores per larva) and controls. Results showed that specific activities of superoxide dismutase, catalase and glutathione S-transferase were significantly higher in the guts of the control larvae than in the guts of the A. apis-infected larvae. The metabolomic data revealed that levels of 28 and 52 metabolites were significantly higher and lower, respectively, in the guts of A. apis-infected larvae than in the guts of control larvae. The 5-oxo-ETE level in the infected larvae was two times higher than that in the control larvae. Elevated 5-oxo-ETE levels may act as a potential metabolic biomarker for chalkbrood disease diagnosis, suggesting that A. apis infection induced obvious oxidative stress in the honey bee larvae. The levels of metabolites such as taurine, docosahexaenoic acid, and L-carnitine involved in combating oxidative stress were significantly decreased in the gut of A. apis-infected larvae. Overall, our results suggest that A. apis infection may compromise the ability of infected larvae to cope with oxidative stress, providing new insight into changing patterns of physiological responses to A. apis infection in honey bee larvae by concurrent use of conventional biochemical assays and untargeted metabolomics.

RNA viruses that contain single-stranded RNA genomes of positive sense make up the largest group of pathogens infecting honey bees. Sacbrood virus (SBV) is one of the most widely distributed honey bee viruses and infects the larvae of honey bees, resulting in failure to pupate and death. Among all of the viruses infecting honey bees, SBV has the greatest number of complete genomes isolated from both European honey bees Apis mellifera and Asian honey bees A. cerana worldwide. To enhance our understanding of the evolution and pathogenicity of SBV, in this study, we present the first report of whole genome sequences of two U.S. strains of SBV. The complete genome sequences of the two U.S. SBV strains were deposited in GenBank under accession numbers: MG545286.1 and MG545287.1. Both SBV strains show the typical genomic features of the Iflaviridae family. The phylogenetic analysis of the single polyprotein coding region of the U.S. strains, and other GenBank SBV submissions revealed that SBV strains split into two distinct lineages, possibly reflecting host affiliation. The phylogenetic analysis based on the 5′UTR revealed a monophyletic clade with the deep parts of the tree occupied by SBV strains from both A. cerane and A. mellifera, and the tips of branches of the tree occupied by SBV strains from A. mellifera. The study of the cold stress on the pathogenesis of the SBV infection showed that cold stress could have profound effects on sacbrood disease severity manifested by increased mortality of infected larvae. This result suggests that the high prevalence of sacbrood disease in early spring may be due to the fluctuating temperatures during the season. This study will contribute to a better understanding of the evolution and pathogenesis of SBV infection in honey bees, and have important epidemiological relevance.

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

The honey bee (Apis mellifera) plays vital ecological roles in the pollination of crops and the maintenance of ecological balance, and adult honey bees may be exposed to exogenous chemicals including heavy metals during their foraging activities. Cadmium (Cd) is regarded as a nonessential toxic metal and is readily accumulated in plants; honey bees can therefore acquire Cd through the collection of contaminated nectar. In the present study, honey bees were chronically exposed to Cd to investigate the effects of sublethal cadmium doses on the olfactory learning and brain gene expression profiles of honey bees. The results showed that Cd-treated bees exhibited significantly impaired olfactory learning performances in comparison with control bees. Moreover, the head weight was significantly lower in Cd-treated bees than in control bees after chronic exposure to Cd. Gene expression profiles between the Cd treatment and the control revealed that 79 genes were significantly differentially expressed. Genes encoding chemoreceptors and olfactory proteins were downregulated, whereas genes involved in response to oxidative stress were upregulated in Cd-treated bees. The results suggest that Cd exposure exerts oxidative stress in the brain of honey bees, and the dysregulated expression of genes encoding chemoreceptors, olfactory proteins, and cytochrome P450 enzymes is probably associated with impaired olfactory learning in honey bees.

Sacbrood virus (SBV) is a common honey bee virus disease. SBV variants and strains identified in Asian honey bees, Apis cerana, have created confusion in identifications. Although the regional names indicated the expansions of the virus in new regions, pathogenesis, and genomes of these variants are not distinct enough to be a separate virus species. However, current SBV qPCR methods may not detect newly identified A. cerana SBV variants (Ac SBV) according to the genome sequences. Since these Ac SBV can naturally infect A. mellifera and possibly other hymenopterans, ignorance of Ac SBV variants in detection methods is simply unwise. In this report, we updated the qPCR method based on Blanchard’s design that used conserved regions of VP1 to design a TaqMan method with an MGB (minor groove binder) probe. We tested the method in bees and hornets, including A. mellifera, A. cerana, and Vespa velutina. The updated primers and the probe can match published SBV and Ac SBV genomes in databases, and this updated method has reasonable sensitivity and flexibility to be applied as a detection and quantification method before the discovery of variants with more mutated VP1 gene.

MicroRNAs (miRNAs) play a vital role in the nerve regulation of honey bees (Apis mellifera). This study aims to investigate the differences in expression of miRNAs in a honey bee’s brain for olfactory learning tasks and to explore their potential role in a honey bee’s olfactory learning and memory. In this study, 12 day old honey bees with strong and weak olfactory performances were utilized to investigate the influence of miRNAs on olfactory learning behavior. The honey bee brains were dissected, and a small RNA-seq technique was used for high-throughput sequencing. The data analysis of the miRNA sequences revealed that 14 differentially expressed miRNAs (DEmiRNAs) between the two groups, strong (S) and weak (W), for olfactory performance in honey bees were identified, which included seven up-regulated and seven down-regulated. The qPCR verification results of the 14 miRNAs showed that four miRNAs (miR-184-3p, miR-276-3p, miR-87-3p, and miR-124-3p) were significantly associated with olfactory learning and memory. The target genes of these DEmiRNAs were subjected to the GO database annotation and KEGG pathway enrichment analyses. The functional annotation and pathway analysis showed that the neuroactive ligand-receptor interaction pathway, oxidative phosphorylation, biosynthesis of amino acids, pentose phosphate pathway, carbon metabolism, and terpenoid backbone biosynthesis may be a great important pathway related to olfactory learning and memory in honey bees. Our findings together further explained the relationship between olfactory performance and the brain function of honey bees at the molecular level and provides a basis for further study on miRNAs related to olfactory learning and memory in honey bees.