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

[This corrects the article DOI: 10.3389/fimmu.2021.747848.].

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 ppm 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.

Metatranscriptome Analysis of Sympatric Bee Species Identifies Bee Virus Variants and a New Virus, Andrena-Associated Bee Virus-1.

Bees are important plant pollinators in agricultural and natural ecosystems. High average annual losses of honey bee (Apis mellifera) colonies in some parts of the world, and regional population declines of some mining bee species (Andrena spp.), are attributed to multiple factors including habitat loss, lack of quality forage, insecticide exposure, and pathogens, including viruses. While research has primarily focused on viruses in honey bees, many of these viruses have a broad host range. It is therefore important to apply a community level approach in studying the epidemiology of bee viruses. We utilized high-throughput sequencing to evaluate viral diversity and viral sharing in sympatric, co-foraging bees in the context of habitat type. Variants of four common viruses (i.e., black queen cell virus, deformed wing virus, Lake Sinai virus 2, and Lake Sinai virus NE) were identified in honey bee and mining bee samples, and the high degree of nucleotide identity in the virus consensus sequences obtained from both taxa indicates virus sharing. We discovered a unique bipartite + ssRNA Tombo-like virus, Andrena-associated bee virus-1 (AnBV-1). AnBV-1 infects mining bees, honey bees, and primary honey bee pupal cells maintained in culture. AnBV-1 prevalence and abundance was greater in mining bees than in honey bees. Statistical modeling that examined the roles of ecological factors, including floral diversity and abundance, indicated that AnBV-1 infection prevalence in honey bees was greater in habitats with low floral diversity and abundance, and that interspecific virus transmission is strongly modulated by the floral community in the habitat. These results suggest that land management strategies that aim to enhance floral diversity and abundance may reduce AnBV-1 spread between co-foraging 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.

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

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

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

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

Honey Bee and Bumble Bee Antiviral Defense.

Bees are important plant pollinators in both natural and agricultural ecosystems. Managed and wild bees have experienced high average annual colony losses, population declines, and local extinctions in many geographic regions. Multiple factors, including virus infections, impact bee health and longevity. The majority of bee-infecting viruses are positive-sense single-stranded RNA viruses. Bee-infecting viruses often cause asymptomatic infections but may also cause paralysis, deformity or death. The severity of infection is governed by bee host immune responses and influenced by additional biotic and abiotic factors. Herein, we highlight studies that have contributed to the current understanding of antiviral defense in bees, including the Western honey bee (Apis mellifera), the Eastern honey bee (Apis cerana) and bumble bee species (Bombus spp.). Bee antiviral defense mechanisms include RNA interference (RNAi), endocytosis, melanization, encapsulation, autophagy and conserved immune pathways including Jak/STAT (Janus kinase/signal transducer and activator of transcription), JNK (c-Jun N-terminal kinase), MAPK (mitogen-activated protein kinases) and the NF-κB mediated Toll and Imd (immune deficiency) pathways. Studies in Dipteran insects, including the model organism Drosophila melanogaster and pathogen-transmitting mosquitos, provide the framework for understanding bee antiviral defense. However, there are notable differences such as the more prominent role of a non-sequence specific, dsRNA-triggered, virus limiting response in honey bees and bumble bees. This virus-limiting response in bees is akin to pathways in a range of organisms including other invertebrates (i.e., oysters, shrimp and sand flies), as well as the mammalian interferon response. Current and future research aimed at elucidating bee antiviral defense mechanisms may lead to development of strategies that mitigate bee losses, while expanding our understanding of insect antiviral defense and the potential evolutionary relationship between sociality and immune function.

Honey bee (Apis mellifera) colony health and pathogen composition in migratory beekeeping operations involved in California almond pollination.

Honey bees are important pollinators of agricultural crops. Pathogens and other factors have been implicated in high annual losses of honey bee colonies in North America and some European countries. To further investigate the relationship between multiple factors, including pathogen prevalence and abundance and colony health, we monitored commercially managed migratory honey bee colonies involved in California almond pollination in 2014. At each sampling event, honey bee colony health was assessed, using colony population size as a proxy for health, and the prevalence and abundance of seven honey bee pathogens was evaluated using PCR and quantitative PCR, respectively. In this sample cohort, pathogen prevalence and abundance did not correlate with colony health, but did correlate with the date of sampling. In general, pathogen prevalence (i.e., the number of specific pathogens harbored within a colony) was lower early in the year (January-March) and was greater in the summer, with peak prevalence occurring in June. Pathogen abundance in individual honey bee colonies varied throughout the year and was strongly associated with the sampling date, and was influenced by beekeeping operation, colony health, and mite infestation level. Together, data from this and other observational cohort studies that monitor individual honey bee colonies and precisely account for sampling date (i.e., day of year) will lead to a better understanding of the influence of pathogens on colony mortality and the effects of other factors on these associations.

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

Recent high annual losses of honey bee colonies are associated with many factors, including RNA virus infections. Honey bee antiviral responses include RNA interference and immune pathway activation, but their relative roles in antiviral defense are not well understood. To better characterize the mechanism(s) of honey bee antiviral defense, bees were infected with a model virus in the presence or absence of dsRNA, a virus associated molecular pattern. Regardless of sequence specificity, dsRNA reduced virus abundance. We utilized next generation sequencing to examine transcriptional responses triggered by virus and dsRNA at three time-points post-infection. Hundreds of genes exhibited differential expression in response to co-treatment of dsRNA and virus. Virus-infected bees had greater expression of genes involved in RNAi, Toll, Imd, and JAK-STAT pathways, but the majority of differentially expressed genes are not well characterized. To confirm the virus limiting role of two genes, including the well-characterized gene, dicer, and a probable uncharacterized cyclin dependent kinase in honey bees, we utilized RNAi to reduce their expression in vivo and determined that virus abundance increased, supporting their involvement in antiviral defense. Together, these results further our understanding of honey bee antiviral defense, particularly the role of a non-sequence specific dsRNA-mediated antiviral pathway.

Pathogen prevalence and abundance in honey bee colonies involved in almond pollination.

Honey bees are important pollinators of agricultural crops. Since 2006, US beekeepers have experienced high annual honey bee colony losses, which may be attributed to multiple abiotic and biotic factors, including pathogens. However, the relative importance of these factors has not been fully elucidated. To identify the most prevalent pathogens and investigate the relationship between colony strength and health, we assessed pathogen occurrence, prevalence, and abundance in Western US honey bee colonies involved in almond pollination. The most prevalent pathogens were Black queen cell virus (BQCV), Lake Sinai virus 2 (LSV2), Sacbrood virus (SBV), Nosema ceranae, and trypanosomatids. Our results indicated that pathogen prevalence and abundance were associated with both sampling date and beekeeping operation, that prevalence was highest in honey bee samples obtained immediately after almond pollination, and that weak colonies had a greater mean pathogen prevalence than strong colonies.

Antiviral Defense Mechanisms in Honey Bees.

Honey bees are significant pollinators of agricultural crops and other important plant species. High annual losses of honey bee colonies in North America and in some parts of Europe have profound ecological and economic implications. Colony losses have been attributed to multiple factors including RNA viruses, thus understanding bee antiviral defense mechanisms may result in the development of strategies that mitigate colony losses. Honey bee antiviral defense mechanisms include RNA-interference, pathogen-associated molecular pattern (PAMP) triggered signal transduction cascades, and reactive oxygen species generation. However, the relative importance of these and other pathways is largely uncharacterized. Herein we review the current understanding of honey bee antiviral defense mechanisms and suggest important avenues for future investigation.

Honey Bee Infecting Lake Sinai Viruses.

Honey bees are critical pollinators of important agricultural crops. Recently, high annual losses of honey bee colonies have prompted further investigation of honey bee infecting viruses. To better characterize the recently discovered and very prevalent Lake Sinai virus (LSV) group, we sequenced currently circulating LSVs, performed phylogenetic analysis, and obtained images of LSV2. Sequence analysis resulted in extension of the LSV1 and LSV2 genomes, the first detection of LSV4 in the US, and the discovery of LSV6 and LSV7. We detected LSV1 and LSV2 in the Varroa destructor mite, and determined that a large proportion of LSV2 is found in the honey bee gut, suggesting that vector-mediated, food-associated, and/or fecal-oral routes may be important for LSV dissemination. Pathogen-specific quantitative PCR data, obtained from samples collected during a small-scale monitoring project, revealed that LSV2, LSV1, Black queen cell virus (BQCV), and Nosema ceranae were more abundant in weak colonies than strong colonies within this sample cohort. Together, these results enhance our current understanding of LSVs and illustrate the importance of future studies aimed at investigating the role of LSVs and other pathogens on honey bee health at both the individual and colony levels.

VPg of murine norovirus binds translation initiation factors in infected cells.

BACKGROUND: Norovirus genomic and subgenomic RNAs are covalently linked at the 5' nucleotide to a 15 kD protein called VPg. VPg of two human norovirus strains binds translation initiation factor eIF3 and other eIFs in vitro, suggesting VPg functions in initiation of protein synthesis on viral RNA. Human norovirus strains are not cultivable, and thus experimental evidence of interactions between VPg and eIFs in infected cells has been lacking. We used the cultivable murine norovirus MNV-1 as a model to study interactions between VPg and eIFs in infected cells. RESULTS: As shown previously for human norovirus VPg, MNV-1 VPg bound eIF3, eIF4GI, eIF4E, and S6 ribosomal protein in cell extracts by GST pull-down assay. Importantly, MNV-1 VPg co-precipitated eIF4GI and eIF4E from infected macrophages, providing evidence that VPg interacts with components of the translation machinery in norovirus infected cells. CONCLUSION: The interactions between MNV-1 VPg and eIFs completely mimic those reported for the human norovirus VPg, illustrating the utility of MNV-1 as a relevant molecular model to study mechanisms of human norovirus replication.

The genome-linked protein VPg of the Norwalk virus binds eIF3, suggesting its role in translation initiation complex recruitment.

The positive-strand RNA genomes of caliciviruses are not capped, but are instead covalently linked at their 5' ends to a viral protein called VPg. The lack of a cap structure typical of eukaryotic mRNA and absence of an internal ribosomal entry site suggest that VPg may function in translation initiation on calicivirus RNA. This hypothesis was tested by analyzing binding of Norwalk virus VPg to translation initiation factors. The eIF3d subunit of eIF3 was identified as a binding partner of VPg by yeast two-hybrid analysis. VPg bound to purified mammalian eIF3 and to eIF3 in mammalian cell lysates. To test the effects of the VPg- eIF3 interaction on translation, VPg was added to cell-free translation reactions programmed with either capped reporter RNA, an RNA containing an EMCV internal ribosomal entry site (IRES) or an RNA with a cricket paralysis virus IRES. VPg inhibited translation of all reporter RNAs in a dose-dependent manner. Together, the data suggest that VPg may play a role in initiating translation on calicivirus RNA through unique protein-protein interactions with the translation machinery.