Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses, having switched hosts from the eastern to western honey bees (Apis cerana to Apis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of deformed wing viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here, we use a phylogeographic analysis based on whole-genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in East Asia and spread in the mid-20th century. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside East Asia and appears absent from the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector's host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalization poses urgent threats to biodiversity and food security.

Globalization has provided opportunities for parasites/pathogens to cross geographic boundaries and expand to new hosts. Recent studies showed that Nosema ceranae, originally considered a microsporidian parasite of Eastern honey bees, Apis cerana, is a disease agent of nosemosis in European honey bees, Apis mellifera, along with the resident species, Nosema apis. Further studies indicated that disease caused by N. ceranae in European honey bees is far more prevalent than that caused by N. apis. In order to gain more insight into the epidemiology of Nosema parasitism in honey bees, we conducted studies to investigate infection of Nosema in its original host, Eastern honey bees, using conventional PCR and duplex real time quantitative PCR methods. Our results showed that A. cerana was infected not only with N. ceranae as previously reported [Fries, I., Feng, F., Silva, A.D., Slemenda, S.B., Pieniazek, N.J., 1996. Nosema ceranae n. sp. (Microspora, Nosematidae), morphological and molecular characterization of a microsporidian parasite of the Asian honey bee Apis cerana (Hymenoptera, Apidae). Eur. J. Protistol. 32, 356-365], but also with N. apis. Both microsporidia produced single and mixed infections. Overall and at each location alone, the prevalence of N. ceranae was higher than that of N. apis. In all cases of mixed infections, the number of N. ceranae gene copies (corresponding to the parasite load) significantly out numbered those of N. apis. Phylogenetic analysis based on a variable region of small subunit ribosomal RNA (SSUrRNA) showed four distinct clades of N. apis and five clades of N. ceranae and that geographical distance does not appear to influence the genetic diversity of Nosema populations. The results from this study demonstrated that duplex real-time qPCR assay developed in this study is a valuable tool for quantitative measurement of Nosema and can be used to monitor the progression of microsprodian infections of honey bees in a timely and cost efficient manner.

Recent steep declines in honey bee health have severely impacted the beekeeping industry, presenting new risks for agricultural commodities that depend on insect pollination. Honey bee declines could reflect increased pressures from parasites and pathogens. The incidence of the microsporidian pathogen Nosema ceranae has increased significantly in the past decade. Here we present a draft assembly (7.86 MB) of the N. ceranae genome derived from pyrosequence data, including initial gene models and genomic comparisons with other members of this highly derived fungal lineage. N. ceranae has a strongly AT-biased genome (74% A+T) and a diversity of repetitive elements, complicating the assembly. Of 2,614 predicted protein-coding sequences, we conservatively estimate that 1,366 have homologs in the microsporidian Encephalitozoon cuniculi, the most closely related published genome sequence. We identify genes conserved among microsporidia that lack clear homology outside this group, which are of special interest as potential virulence factors in this group of obligate parasites. A substantial fraction of the diminutive N. ceranae proteome consists of novel and transposable-element proteins. For a majority of well-supported gene models, a conserved sense-strand motif can be found within 15 bases upstream of the start codon; a previously uncharacterized version of this motif is also present in E. cuniculi. These comparisons provide insight into the architecture, regulation, and evolution of microsporidian genomes, and will drive investigations into honey bee-Nosema interactions.

To clarify the mechanisms of Nosema ceranae parasitism, we deep-sequenced both honey bee host and parasite mRNAs throughout a complete 6-day infection cycle. By time-series analysis, 1122 parasite genes were significantly differently expressed during the reproduction cycle, clustering into 4 expression patterns. We found reactive mitochondrial oxygen species modulator 1 of the host to be significantly down regulated during the entire infection period. Our data support the hypothesis that apoptosis of honey bee cells was suppressed during infection. We further analyzed genome-wide genetic diversity of this parasite by comparing samples collected from the same site in 2007 and 2013. The number of SNP positions per gene and the proportion of non-synonymous substitutions per gene were significantly reduced over this time period, suggesting purifying selection on the parasite genome and supporting the hypothesis that a subset of N. ceranae strains might be dominating infection.

Transmission routes impact pathogen virulence and genetics, therefore comprehensive knowledge of these routes and their contribution to pathogen circulation is essential for understanding host-pathogen interactions and designing control strategies. Deformed wing virus (DWV), a principal viral pathogen of honey bees associated with increased honey bee mortality and colony losses, became highly virulent with the spread of its vector, the ectoparasitic mite Varroa destructor. Reproduction of Varroa mites occurs in capped brood cells and mite-infested pupae from these cells usually have high levels of DWV. The removal of mite-infested pupae by worker bees, Varroa Sensitive Hygiene (VSH), leads to cannibalization of pupae with high DWV loads, thereby offering an alternative route for virus transmission. We used genetically tagged DWV to investigate virus transmission to and between worker bees following pupal cannibalisation under experimental conditions. We demonstrated that cannibalization of DWV-infected pupae resulted in high levels of this virus in worker bees and that the acquired virus was then transmitted between bees via trophallaxis, allowing circulation of Varroa-vectored DWV variants without the mites. Despite the known benefits of hygienic behaviour, it is possible that higher levels of VSH activity may result in increased transmission of DWV via cannibalism and trophallaxis.

Microsporidia comprise a phylum of single cell, intracellular parasites and represent the earliest diverging branch in the fungal kingdom. The microsporidian parasite Nosema ceranae primarily infects honey bee gut epithelial cells, leading to impaired memory, suppressed host immune responses and colony collapse under certain circumstances. As the genome of N. ceranae is challenging to assembly due to very high genetic diversity and repetitive region, the genome was re-sequenced using long reads. We present a robust 8.8 Mbp genome assembly of 2,280 protein coding genes, including a high number of genes involved in transporting nutrients and energy, as well as drug resistance when compared with sister species Nosema apis. We also describe the loss of the critical protein Dicer in approximately half of the microsporidian species, giving new insights into the availability of RNA interference pathway in this group. Our results provided new insights into the pathogenesis of N. ceranae and a blueprint for treatment strategies that target this parasite without harming honey bees. The unique infectious apparatus polar filament and transportation pathway members can help to identify treatments to control this parasite.

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.

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.

Annual losses of honey bee colonies remain high and pesticide exposure is one possible cause. Dangerous combinations of pesticides, plant-produced compounds and antibiotics added to hives may cause or contribute to losses, but it is very difficult to test the many combinations of those compounds that bees encounter. We propose a mechanism-based strategy for simplifying the assessment of combinations of compounds, focusing here on compounds that interact with xenobiotic handling ABC transporters. We evaluate the use of ivermectin as a model substrate for these transporters. Compounds that increase sensitivity of bees to ivermectin may be inhibiting key transporters. We show that several compounds commonly encountered by honey bees (fumagillin, Pristine, quercetin) significantly increased honey bee mortality due to ivermectin and significantly reduced the LC50 of ivermectin suggesting that they may interfere with transporter function. These inhibitors also significantly increased honey bees sensitivity to the neonicotinoid insecticide acetamiprid. This mechanism-based strategy may dramatically reduce the number of tests needed to assess the possibility of adverse combinations among pesticides. We also demonstrate an in vivo transporter assay that provides physical evidence of transporter inhibition by tracking the dynamics of a fluorescent substrate of these transporters (Rhodamine B) in bee tissues. Significantly more Rhodamine B remains in the head and hemolymph of bees pretreated with higher concentrations of the transporter inhibitor verapamil. Mechanism-based strategies for simplifying the assessment of adverse chemical interactions such as described here could improve our ability to identify those combinations that pose significantly greater risk to bees and perhaps improve the risk assessment protocols for honey bees and similar sensitive species.

Varroa mites are widely considered the biggest honey bee health problem worldwide. Until recently, Varroa jacobsoni has been found to live and reproduce only in Asian honey bee (Apis cerana) colonies, while V. destructor successfully reproduces in both A. cerana and A. mellifera colonies. However, we have identified an island population of V. jacobsoni that is highly destructive to A. mellifera, the primary species used for pollination and honey production. The ability of these populations of mites to cross the host species boundary potentially represents an enormous threat to apiculture, and is presumably due to genetic variation that exists among populations of V. jacobsoni that influences gene expression and reproductive status. In this work, we investigate differences in gene expression between populations of V. jacobsoni reproducing on A. cerana and those either reproducing or not capable of reproducing on A. mellifera, in order to gain insight into differences that allow V. jacobsoni to overcome its normal species tropism.

The small hive beetle (SHB), Aethina tumida, is a major pest of managed honey bee (Apis mellifera) colonies in the United States and Australia, and an emergent threat in Europe. While strong honey bee colonies generally keep SHB populations in check, weak or stressed colonies can succumb to infestations. This parasite has spread from a sub-Saharan Africa to three continents, leading to immense management and regulatory costs. We performed a transcriptomic analysis involving deep sequencing of multiple life stages and both sexes of this species. The assembled transcriptome appears to be nearly complete, as judged by conserved insect orthologs and the ability to find plausible homologs for 11,952 proteins described from the genome of the red flour beetle. Expressed genes include each of the major metabolic, developmental and sensory groups, along with genes for proteins involved with immune defenses and insecticide resistance. We also present a total of 23,085 high-quality SNP's for the assembled contigs. We highlight potential differences between this beetle and its honey bee hosts, and suggest mechanisms of future research into the biology and control of this species. SNP resources will allow functional genetic analyses and analyses of dispersal for this invasive pest. All resources are posted as Supplemental Tables at https://data.nal.usda.gov/dataset/data-transcriptomic-and-functional-resources-small-hive-beetle-aethina-tumida-worldwide, and at NCBI under Bioproject PRJNA256171.

As scientists continue to pursue various 'omics-based research, there is a need for high quality data for the most fundamental 'omics of all: genomics. The bacterium Paenibacillus larvae is the causative agent of the honey bee disease American foulbrood. If untreated, it can lead to the demise of an entire hive; the highly social nature of bees also leads to easy disease spread, between both individuals and colonies. Biologists have studied this organism since the early 1900s, and a century later, the molecular mechanism of infection remains elusive. Transcriptomics and proteomics, because of their ability to analyze multiple genes and proteins in a high-throughput manner, may be very helpful to its study. However, the power of these methodologies is severely limited without a complete genome; we undertake to address that deficiency here.

The small hive beetle (SHB), Aethina tumida, has emerged as a worldwide threat to honey bees in the past two decades. These beetles harvest nest resources, feed on larval bees, and ultimately spoil nest resources with gelatinous slime together with the fungal symbiont Kodamaea ohmeri.

The ongoing decline of honey bee health worldwide is a serious economic and ecological concern. One major contributor to the decline are pathogens, including several honey bee viruses. However, information is limited on the biology of bee viruses and molecular interactions with their hosts. An experimental protocol to test these systems was developed, using injections of Israeli Acute Paralysis Virus (IAPV) into honey bee pupae reared ex-situ under laboratory conditions. The infected pupae developed pronounced but variable patterns of disease. Symptoms varied from complete cessation of development with no visual evidence of disease to rapid darkening of a part or the entire body. Considerable differences in IAPV titer dynamics were observed, suggesting significant variation in resistance to IAPV among and possibly within honey bee colonies. Thus, selective breeding for virus resistance should be possible. Gene expression analyses of three separate experiments suggest IAPV disruption of transcriptional homeostasis of several fundamental cellular functions, including an up-regulation of the ribosomal biogenesis pathway. These results provide first insights into the mechanisms of IAPV pathogenicity. They mirror a transcriptional survey of honey bees afflicted with Colony Collapse Disorder and thus support the hypothesis that viruses play a critical role in declining honey bee health.

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.

Honey bee larvae are highly susceptible to the bacterial pathogen Paenibacillus larvae only during the first instar of larval development. Transcript levels were measured for genes encoding two antimicrobial peptides, abaecin and defensin, as well as for two candidates in the immune response cascade (PGRP-LD and masquerade) in control larvae and larvae exposed to the pathogen. Transcripts for all four are present throughout development. This suggests that other physiological or dietary factors may better explain the age-based change in vulnerability to this pathogen. One of these genes, abaecin, shows significant up-regulation 24 h following oral inoculation with P. larvae, precisely when the bacterium surmounts the midgut epithelia of bees. Expression of both antimicrobial peptides varied by 1000-fold across different nestmate bees, indicating an allelic component to their expression. The implications of these results for current hypotheses related to disease tolerance in social insects are discussed, along with implications for breeding bees resistant to this important disease.

Israeli acute paralysis virus (IAPV) is a widespread RNA virus of honey bees that has been linked with colony losses. Here we describe the transmission, prevalence, and genetic traits of this virus, along with host transcriptional responses to infections. Further, we present RNAi-based strategies for limiting an important mechanism used by IAPV to subvert host defenses. Our study shows that IAPV is established as a persistent infection in honey bee populations, likely enabled by both horizontal and vertical transmission pathways. The phenotypic differences in pathology among different strains of IAPV found globally may be due to high levels of standing genetic variation. Microarray profiles of host responses to IAPV infection revealed that mitochondrial function is the most significantly affected biological process, suggesting that viral infection causes significant disturbance in energy-related host processes. The expression of genes involved in immune pathways in adult bees indicates that IAPV infection triggers active immune responses. The evidence that silencing an IAPV-encoded putative suppressor of RNAi reduces IAPV replication suggests a functional assignment for a particular genomic region of IAPV and closely related viruses from the Family Dicistroviridae, and indicates a novel therapeutic strategy for limiting multiple honey bee viruses simultaneously and reducing colony losses due to viral diseases. We believe that the knowledge and insights gained from this study will provide a new platform for continuing studies of the IAPV-host interactions and have positive implications for disease management that will lead to mitigation of escalating honey bee colony losses worldwide.

Honey bees, Apis mellifera, face many parasites and pathogens and consequently rely on a diverse set of individual and group-level defenses to prevent disease. One route by which honey bees and other insects might combat disease is through the shielding effects of their microbial symbionts. Bees carry a diverse assemblage of bacteria, very few of which appear to be pathogenic. Here we explore the inhibitory effects of these resident bacteria against the primary bacterial pathogen of honey bees, Paenibacillus larvae.

The ectoparasitic mite, Varroa destructor, is the most destructive parasite of managed honeybee colonies worldwide. Since V. destructor transfers pathogens to honeybees, it may be adaptive for bees to respond to mite infestation by upregulating their immune responses. Mites, however, may overcome the host's immune responses by suppressing them, which could facilitate the mite's ability to feed on hemolymph. A humoral immune response of bees parasitized by V. destructor may be detected by studying the expression levels of antibacterial peptides, such as abaecin and defensin, known to be immune-responsive. Expression levels for these two antibacterial peptides changed non-linearly with respect to the number of mites parasitizing honeybee pupae. Bees exposed to low or moderate number of mites had fewer immune-related transcripts than pupae that were never parasitized or pupae with high mite loads. Although many of the pupae tested indicated the presence of bacteria, no correlation with mite numbers or immune-response levels existed. All bees tested negative for acute paralysis and Kashmir bee viruses known to be vectored by V. destructor.

Nosemosis C, a Nosema disease caused by microsporidia parasite Nosema ceranae, is a significant disease burden of the European honey bee Apis mellifera which is one of the most economically important insect pollinators. Nevertheless, there is no effective treatment currently available for Nosema disease and the disease mechanisms underlying the pathological effects of N. ceranae infection in honey bees are poorly understood. Iron is an essential nutrient for growth and survival of hosts and pathogens alike. The iron tug-of-war between host and pathogen is a central battlefield at the host-pathogen interface which determines the outcome of an infection, however, has not been explored in honey bees. To fill the gap, we conducted a study to investigate the impact of N. ceranae infection on iron homeostasis in honey bees. The expression of transferrin, an iron binding and transporting protein that is one of the key players of iron homeostasis, in response to N. ceranae infection was analysed. Furthermore, the functional roles of transferrin in iron homeostasis and honey bee host immunity were characterized using an RNA interference (RNAi)-based method. The results showed that N. ceranae infection causes iron deficiency and upregulation of the A. mellifera transferrin (AmTsf) mRNA in honey bees, implying that higher expression of AmTsf allows N. ceranae to scavenge more iron from the host for its proliferation and survival. The suppressed expression levels of AmTsf via RNAi could lead to reduced N. ceranae transcription activity, alleviated iron loss, enhanced immunity, and improved survival of the infected bees. The intriguing multifunctionality of transferrin illustrated in this study is a significant contribution to the existing body of literature concerning iron homeostasis in insects. The uncovered functional role of transferrin on iron homeostasis, pathogen growth and honey bee's ability to mount immune responses may hold the key for the development of novel strategies to treat or prevent diseases in honey bees.