Divergent forms of endoplasmic reticulum stress trigger a robust unfolded protein response in honey bees

Preliminary Analysis of Genetic Markers for Functional Ethanol Tolerance in Honey Bees (Apis mellifera) Using a Free-Flying Paradigm

Honey bees are a commonly used species for alcohol research due to their genome being fully sequenced, their behavioral changes following consumption, and their preference for alcohol. The purpose of this article is to provide a preliminary examination of the genetic expression of heat shock protein 70 (HSP70) and big potassium ion channel protein (BKP) in honey bees following the consumption of either 0%, 2.5%, 5%, or 10% ethanol (EtOH) solutions. The foraging behaviors of the bees were observed and recorded through their return and drinking times. There were significant differences in the return and drinking times between some of the groups. The bees in the 10% condition took significantly longer to return compared to the other groups. Additionally, the bees in the 5% group spent significantly more time drinking compared to the bees in the control (0%) group. There were no significant differences in HSP70 or BKP between the different ethanol groups. Cumulatively, these findings suggest that, while bees may exhibit behavioral differences, the differences in gene expression may not be observed at the transcriptional level.

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.

The response of heat shock proteins in honey bees to abiotic and biotic stressors

Honey bees, Apis mellifera, are the most important managed pollinators worldwide. They are highly impacted by various abiotic and biotic stressors, especially temperature extremes, which can lead to cellular damage and death. The induction of heat shock proteins (HSPs) has been recorded in honey bees as a response to various types of stressors. HSPs are classified into different gene families according to their molecular weights. HSPs play an important role in maintaining cellular protein homeostasis due to their contribution as molecular chaperones or co-chaperones. HSPs in honey bees have complex functions with induction even under normal colony conditions. Previous studies have suggested various functions of HSPs to protect cells from damage under exposure to environmental stressors, pollutants, and pathogens. Surprisingly, HSPs have also been found to play roles in larval development and age-related tasks. The expression of HSPs varies depending on tissue type, developmental stage, age, and stress period. This article reviews studies on HSPs (sHSPs, HSP40, HSP60, HSP70, and HSP90) in honey bees and highlights gaps in the available knowledge. This review is crucial for honey bee research, particularly in the face of climate change challenges.

Mixture toxic effects of thiacloprid and cyproconazole on honey bees (Apis mellifera L.)

Pesticide exposure remains one of the main factors in the population decline of insect pollinators. It is urgently necessary to assess the effects of mixtures on pollinator risk assessments because they are often exposed to numerous agrochemicals. In the present study, we explored the mixture toxic effects of thiacloprid (THI) and cyproconazole (CYP) on honey bees (Apis mellifera L.). Our findings revealed that THI possessed higher acute toxicity to A. mellifera (96-h LC₅₀ value of 216.3 mg a.i. L^-1) than CYP (96-h LC₅₀ value of 601.4 mg a.i. L^-1). It's worth noting that the mixture of THI and CYP exerted an acute synergistic effect on honey bees. At the same time, the activities of detoxification enzyme cytochrome P450s (CYP450s) and neuro target enzyme Acetylcholinesterase (AChE), as well as the expressions of seven genes (CRBXase, CYP306A1, CYP6AS14, apidaecin, defensing-2, vtg, and gp-93) associated with detoxification metabolism, immune response, development, and endoplasmic reticulum stress, were significantly altered in the combined treatment compared with the corresponding individual exposures of THI or CYP. These data indicated that a mixture of THI and CYP could disturb the physiological homeostasis of honey bees. Our study provides a theoretical basis for in-depth studies on the impacts of pesticide mixtures on the health of honey bees. Our study also provides important guidance for the rational application of pesticide mixtures to protect pollinators in agricultural production effectively.

Halofuginone triggers a transcriptional program centered on ribosome biogenesis and function in honey bees

We previously found that pharmacological inhibition of prolyl-tRNA synthetase by halofuginone has potent activity against Nosema ceranae, an important pathogen of honey bees. However, we also observed that prolyl-tRNA synthetase inhibition is toxic to bees, suggesting further work is necessary to make this a feasible therapeutic strategy. As expected, we found that pharmacological inhibition of prolyl-tRNA synthetase activity resulted in robust induction of select canonical ATF4 target genes in honey bees. However, our understanding of this and other cellular stress responses in general in honey bees is incomplete. Thus, we used RNAseq to identify novel changes in gene expression after halofuginone treatment and observed induction of genes involved in ribosome biogenesis, translation, tRNA synthesis, and ribosome-associated quality control (RQC). These results suggest that halofuginone, potentially acting through the Integrated Stress Response (ISR), promotes a transcriptional response to ribosome functional impairment in honey bees rather than the response designed to oppose amino acid limitation, which has been observed in other organisms after ISR induction. In support of this idea, we found that cycloheximide (CHX) administration also induced all tested target genes, indicating that this gene expression program could be induced by ribosome stalling in addition to tRNA synthetase inhibition. Only a subset of halofuginone-induced genes was upregulated by Unfolded Protein Response (UPR) induction, suggesting that mode of activation and cross-talk with other cellular signaling pathways significantly influence ISR function and cellular response to its activation. Future work will focus on understanding how the apparently divergent transcriptional output of the ISR in honey bees impacts the health and disease of this important pollinator species.

Honey bee sHSP are responsive to diverse proteostatic stresses and potentially promising biomarkers of honey bee stress

The pollination services provided by the honey bee are critical in both natural and agricultural ecosystems. Honey bee colonies in the United States have suffered from an increased rate of die-off in recent years, stemming from a complex set of interacting stresses that remain poorly described. Defining specific common cellular processes and cellular stress responses impacted by multiple stressors represent a key step in understanding these synergies. Proteotoxic stresses negatively impact protein synthesis, folding, and degradation. Diverse proteotoxic stresses induce expression of genes encoding small heat shock proteins (sHSP) of the expanded lethal (2) essential for life (l(2)efl) gene family. In addition to upregulation by the Integrated Stress Response (ISR), the Heat Shock Response (HSR), and the Oxidative Stress Response (OSR), our data provide first evidence that sHSP genes are upregulated by the Unfolded Protein Response (UPR). As these genes appear to be part of a core stress response that could serve as a useful biomarker for cellular stress in honey bees, we designed and tested an RT-LAMP assay to detect increased l(2)efl gene expression in response to heat-stress. While this assay provides a powerful proof of principle, further work will be necessary to link changes in sHSP gene expression to colony-level outcomes, to adapt our preliminary assay into a Point of Care Testing (POCT) assay appropriate for use as a diagnostic tool for use in the field, and to couple assay results to management recommendations.

Proteasome Inhibition Is an Effective Treatment Strategy for Microsporidia Infection in Honey Bees

The microsporidia Nosema ceranae is an obligate intracellular parasite that causes honey bee mortality and contributes to colony collapse. Fumagillin is presently the only pharmacological control for N. ceranae infections in honey bees. Resistance is already emerging, and alternative controls are critically needed. Nosema spp. exhibit increased sensitivity to heat shock, a common proteotoxic stress. Thus, we hypothesized that targeting the Nosema proteasome, the major protease removing misfolded proteins, might be effective against N. ceranae infections in honey bees. Nosema genome analysis and molecular modeling revealed an unexpectedly compact proteasome apparently lacking multiple canonical subunits, but with highly conserved proteolytic active sites expected to be receptive to FDA-approved proteasome inhibitors. Indeed, N. ceranae were strikingly sensitive to pharmacological disruption of proteasome function at doses that were well tolerated by honey bees. Thus, proteasome inhibition is a novel candidate treatment strategy for microsporidia infection in honey bees.

Diet-derived transmission of MicroRNAs from host plant into honey bee Midgut

Background

MicroRNA (miRNA) is a class of small noncoding RNAs, which targets on thousands of mRNA and thus plays important roles in many biological processes. It has been reported that miRNA has cross-species regulation functions between parasitoid-host, or plant-animal, etc. For example, several plant miRNAs enter into the honey bees and regulate gene expression. However, whether cross-species regulation function of miRNAs is a universal mechanism remains a debate question.

Results

We have evaluated transmission of miRNAs from sunflower and sedr plants into the midgut of honey bee using RNA-Seq analyses complemented with confirmation by RT-qPCR. The results showed that at least 11 plant miRNAs were found in the midgut of honey bee feeding by sunflower and sedr pollen. Among which, nine miRNAs, including miR-30d, miR-143, miR-148a, miR-21, let-7 g, miR-26a, miR-126, miR-27a, and miR-203, were shared between the sunflower- and sedr-fed honey bees, suggesting they might have essential roles in plant-insect interactions. Moreover, existence of these co-shared miRNAs presents a strong evidence to support the successful transmission of miRNAs into the midgut of the insect. In total, 121 honeybee mRNAs were predicted to be the target of these 11 plant-derived miRNAs. Interestingly, a sedr-derived miRNA, miR-206, targets on 53 honeybee genes. Kyoto Encyclopedia of Genes and Genome (KEGG) analyses showed that these target genes are significantly involved in hippo signaling pathway-fly, Wnt signaling pathway, and N-Glycan biosynthesis.

Conclusions

In summary, these results provide evidence of cross-species regulation function of miRNA between honeybee and flowering host plants, extending our understanding of the molecular interactions between plants and animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07916-4.

The IRE1 pathway regulates honey bee Unfolded Protein Response gene expression

Our molecular understanding of honey bee cellular stress responses is incomplete. Previously, we sought to identify and began functional characterization of the components of the Unfolded Protein Response (UPR) in honey bees. We observed that UPR stimulation resulted in induction of target genes upon IRE1 pathway activation, as assessed by splicing of Xbp1 mRNA. However, we were not able to determine the relative role of the various UPR pathways in gene activation. Our understanding of honey bee signal transduction and transcriptional regulation has been hampered by a lack of tools. After using RNA-seq to expand the known UPR targets in the honey bee, we used the Drosophila melanogaster S2 cell line and honey bee trans and cis elements to investigate the role of the IRE1 pathway in the transcriptional activation of one of these targets, the honey bee Hsc70-3 gene. Using a luciferase reporter, we show that honey bee Hsc70 promoter activity is inducible by UPR activation. In addition, we show that this activation is IRE1-dependent and relies on specific cis regulatory elements. Experiments using exogenous honey bee or fruit fly XBP1S proteins demonstrate that both factors can activate the Hsc70-3 promoter and further support a role for the IRE1 pathway in control of Hsc70-3 expression in the honey bee. By providing foundational knowledge about the UPR in the honey bee and demonstrating the usefulness of a heterologous cell line for molecular characterization of honey bee pathways, this work stands to improve our understanding of this critical species.

Molecular characterization and expression pattern of inositol-requiring enzyme 1 (IRE1) in blunt snout bream (Megalobrama amblycephala): its role of IRE1 involved in inflammatory response induced by lipopolysaccharide

This study aimed to characterize the full-length cDNA of IRE1 from fish Megalobrama amblycephala and investigate its role in the pro-inflammatory response. A full-length cDNA coding IRE1 was cloned from blunt snout bream by RT-PCR and RACE approaches. The cDNA obtained covered 3665 bp with an open reading frame of 3096 bp encoding 1031 amino acids. Sequence alignment and phylogenetic analysis revealed a high degree of conservation (74-92%) among various species, retaining one signal peptide, one luminal domain, one serine/threonine kinase domain, one RNase domain, one activation loop, two N-linked glycosylation sites, and several phosphorylation sites. The highest IRE1 expression was observed in the trunk kidney followed by the brain and spleen, whereas relatively low expression levels were detected in the liver, intestine, adipose, skin, and heart. After lipopolysaccharide (LPS) challenge, the expressions of glucose-regulated protein 78 (GRP78), inositol-requiring enzyme 1 (IRE1), spliced X-box binding protein 1 (XBP1s), C/EBP homologous protein (CHOP), nuclear factor kappa B (NF-κB), tumor necrosis factor alpha (TNFα), and interleukin-6 (IL-6) all increased remarkably in the spleen and brain at different sampling time points, while LPS also upregulated all the genes tested in the intestine except C/EBP homologous protein. Overall, the results indicated that the IRE1 gene of Megalobrama amblycephala shared a high similarity compared with other vertebrates including several bony fish species. Its expression in three tissues was induced remarkably by the LPS challenge, which indicated that IRE1 played a vital role in LPS-induced inflammation on fish.

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.

Acute thiamethoxam toxicity in honeybees is not enhanced by common fungicide and herbicide and lacks stress-induced changes in mRNA splicing

Securing food supply for a growing population is a major challenge and heavily relies on the use of agrochemicals to maximize crop yield. It is increasingly recognized, that some neonicotinoid insecticides have a negative impact on non-target organisms, including important pollinators such as the European honeybee Apis mellifera. Toxicity of neonicotinoids may be enhanced through simultaneous exposure with additional pesticides, which could help explain, in part, the global decline of honeybee colonies. Here we examined whether exposure effects of the neonicotinoid thiamethoxam on bee viability are enhanced by the commonly used fungicide carbendazim and the herbicide glyphosate. We also analysed alternative splicing changes upon pesticide exposure in the honeybee. In particular, we examined transcripts of three genes: (i) the stress sensor gene X box binding protein-1 (Xbp1), (ii) the Down Syndrome Cell Adhesion Molecule (Dscam) gene and iii) the embryonic lethal/abnormal visual system (elav) gene, which are important for neuronal function. Our results showed that acute thiamethoxam exposure is not enhanced by carbendazim, nor glyphosate. Toxicity of the compounds did not trigger stress-induced, alternative splicing in the analysed mRNAs, thereby leaving dormant a cellular response pathway to these man-made environmental perturbations.

Robust Transcriptional Response to Heat Shock Impacting Diverse Cellular Processes despite Lack of Heat Shock Factor in Microsporidia

The majority of fungal species prefer the 12° to 30°C range, and relatively few species tolerate temperatures higher than 35°C. Our understanding of the mechanisms underpinning the ability of some species to grow at higher temperatures is incomplete. Nosema ceranae is an obligate intracellular fungal parasite that infects honey bees and can cause individual mortality and contribute to colony collapse. Despite a reduced genome, this species is strikingly thermotolerant, growing optimally at the colony temperature of 35°C. In characterizing the heat shock response (HSR) in N. ceranae, we found that this and other microsporidian species have lost the transcriptional regulator HSF and possess a reduced set of putative core HSF1-dependent HSR target genes. Despite these losses, N. ceranae demonstrates robust upregulation of the remaining HSR target genes after heat shock. In addition, thermal stress leads to alterations in genes involved in various metabolic pathways, ribosome biogenesis and translation, and DNA repair. These results provide important insight into the stress responses of microsporidia. Such a new understanding will allow new comparisons with other pathogenic fungi and potentially enable the discovery of novel treatment strategies for microsporidian infections affecting food production and human health.IMPORTANCE We do not fully understand why some fungal species are able to grow at temperatures approaching mammalian body temperature. Nosema ceranae, a microsporidium, is a type of fungal parasite that infects honey bees and grows optimally at the colony temperature of 35°C despite possessing cellular machinery for responding to heat stress that is notably simpler than that of other fungi. We find that N. ceranae demonstrates a robust and broad response to heat shock. These results provide important insight into the stress responses of this type of fungus, allow new comparisons with other pathogenic fungi, and potentially enable the discovery of novel treatment strategies for this type of fungus.

Inositol-requiring enzyme 1 involved in regulating hemocyte apoptosis upon heat stress in Patinopecten yessoensis

The inositol-requiring enzyme 1 (IRE1), one of the primary endoplasmic reticulum (ER) transmembrane receptor proteins, is involved in regulating unfolded protein response (UPR) signaling pathway and plays an import role in maintaining cell homeostasis. In the present study, an IRE1 homologue was identified from Patinopecten yessoensis (designated as PyIRE1). The cDNA of PyIRE1 was of 3314 bp with a 2646 bp open reading frame (ORF) of IRE1 encoding a polypeptide of 881 amino acids. There was a signal peptide, four pyrrolo-quinoline quinine (PPQ) domains, a transmembrane helix region, a Serine/Threonine protein kinases domain (S_TKc) and a protein kinases or N-glycanases containing protein domain (PUG) in the deduced amino acid sequence of PyIRE1. The PyIRE1 mRNA was constitutively expressed in all the tested tissues, with the highest expression level in gills. PyIRE1 protein was mainly located in the ER of P. yessoensis hemocytes. The expression profiles of PyIRE1, glucose-regulated protein 94 (designated as PyGRP94) and glucose-regulated protein 78 (designated as PyGRP78) were determined by SYBR Green qRT-PCR after heat shock treatment. The mRNA expression levels of all these three genes were significantly up-regulated and reached their peak values at 2 h (3.97-fold, p < 0.05), 8 h (19.67-fold, p < 0.05) and 4 h (27.37-fold, p < 0.05) in hemocytes, 2 h (3.55-fold, p < 0.05), 12 h (8.58-fold, p < 0.05) and 8 h (35.31-fold, p < 0.05) in gills after heat shock treatment, respectively. After the injection with PyIRE1 dsRNA, the mRNA expression of pro-apoptotic B-cell lymphoma-2 (Bcl-2) family member PyBax and the activity of caspase-3 significantly decreased in comparison with the control group (p < 0.05) after heat shock treatment. These results collectively suggested that PyIRE1, as an ER stress sensor, was potentially involved in the response upon heat stress by regulating the expression of PyBax and apoptosis of hemocytes in P. yessoensis.

The heat shock response and humoral immune response are mutually antagonistic in honey bees

The honey bee is of paramount importance to humans in both agricultural and ecological settings. Honey bee colonies have suffered from increased attrition in recent years, stemming from complex interacting stresses. Defining common cellular stress responses elicited by these stressors represents a key step in understanding potential synergies. The proteostasis network is a highly conserved network of cellular stress responses involved in maintaining the homeostasis of protein production and function. Here, we have characterized the Heat Shock Response (HSR), one branch of this network, and found that its core components are conserved. In addition, exposing bees to elevated temperatures normally encountered by honey bees during typical activities results in robust HSR induction with increased expression of specific heat shock proteins that was variable across tissues. Surprisingly, we found that heat shock represses multiple immune genes in the abdomen and additionally showed that wounding the cuticle of the abdomen results in decreased expression of multiple HSR genes in proximal and distal tissues. This mutually antagonistic relationship between the HSR and immune activation is unique among invertebrates studied to date and may promote understanding of potential synergistic effects of disparate stresses in this critical pollinator and social insects more broadly.

Exposure of honey bees (Apis mellifera) to different classes of insecticides exhibit distinct molecular effect patterns at concentrations that mimic environmental contamination

Pesticides are implicated in the decline of honey bee populations. Many insecticides are neurotoxic and act by different modes of actions. Although a link between insecticide exposure and changed behaviour has been made, molecular effects underlying these effects are poorly understood. Here we elucidated molecular effects at environmental realistic concentrations of two organophosphates, chlorpyrifos and malathion, the pyrethroid cypermethrin, and the ryanodine receptor activator, chlorantraniliprole. We assessed transcriptional alterations of selected genes at three exposure times (24 h, 48 h, 72 h) in caged honey bees exposed to different concentrations of these compounds. Our targeted gene expression concept focused on several transcripts, including nicotinic acetylcholine receptor α 1 and α 2 (nAChRα1, nAChRα2) subunits, the multifunctional gene vitellogenin, immune system related genes of three immune system pathways, genes belonging to the detoxification system and ER stress genes. Our data indicate a dynamic pattern of expressional changes at different exposure times. All four insecticides induced strong alterations in the expression of immune system related genes suggesting negative implications for honey bee health, as well as cytochrome P450 enzyme transcripts suggesting an interference with metabolism. Exposure to neurotoxic chlorpyrifos, malathion and cypermethrin resulted in up-regulation of nAChRα1 and nAChRα2. Moreover, alterations in the expression of vitellogenin occurred, which suggests implications on foraging activity. Chlorantraniliprole induced ER stress which may be related to toxicity. The comparison of all transcriptional changes indicated that the expression pattern is rather compound-specific and related to its mode of action, but clusters of common transcriptional changes between different compounds occurred. As transcriptional alterations occurred at environmental concentrations our data provide a molecular basis for observed adverse effects of these insecticides to bees.