Immune Response of Eastern Honeybee Worker to Nosema ceranae Infection Revealed by Transcriptomic Investigation

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.

CircRNA-regulated immune responses of asian honey bee workers to microsporidian infection

Nosema ceranae is a widespread fungal parasite for honey bees, causing bee nosemosis. Based on deep sequencing and bioinformatics, identification of circular RNAs (circRNAs) in Apis cerana workers’ midguts and circRNA-regulated immune response of host to N. ceranae invasion were conducted in this current work, followed by molecular verification of back-splicing sites and expression trends of circRNAs. Here, 10185 and 7405 circRNAs were identified in the midguts of workers at 7 days (AcT1) and 10 days (AcT2) post inoculation days post-inoculation with N. ceranae. PCR amplification result verified the back-splicing sites within three specific circRNAs (novel_circ_005123, novel_circ_007177, and novel_circ_015140) expressed in N. ceranae-inoculated midgut. In combination with transcriptome data from corresponding un-inoculated midguts (AcCK1 and AcCK2), 2266 circRNAs were found to be shared by the aforementioned four groups, whereas the numbers of specific ones were 2618, 1917, 5691, and 3723 respectively. Further, 83 52) differentially expressed circRNAs (DEcircRNAs) were identified in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group. Source genes of DEcircRNAs in workers’ midgut at seven dpi were involved in two cellular immune-related pathways such as endocytosis and ubiquitin mediated proteolysis. Additionally, competing endogenous RNA (ceRNA) network analysis showed that 23 13) DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group could target 18 14) miRNAs and further link to 1111 (1093) mRNAs. These target mRNAs were annotated to six cellular immunity pathways including endocytosis, lysosome, phagosome, ubiquitin mediated proteolysis, metabolism of xenobiotics by cytochrome P450, and insect hormone biosynthesis. Moreover, 284 164) internal ribosome entry site and 54 26) ORFs were identified from DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group; additionally, ORFs in DEcircRNAs in midgut at seven dpi with N. ceranae were associated with several cellular immune pathways including endocytosis and ubiquitin-mediated proteolysis. Ultimately, RT-qPCR results showed that the expression trends of six DEcircRNAs were consistent with those in transcriptome data. These results demonstrated that N. ceranae altered the expression pattern of circRNAs in A. c. cerana workers’ midguts, and DEcircRNAs were likely to regulate host cellular and humoral immune response to microsporidian infection. Our findings lay a foundation for clarifying the mechanism underlying host immune response to N. ceranae infection and provide a new insight into interaction between Asian honey bee and microsporidian.

In-depth investigation of microRNA-mediated cross-kingdom regulation between Asian honey bee and microsporidian

Asian honey bee Apis cerana is the original host for Nosema ceranae, a unicellular fungal parasite that causes bee nosemosis throughout the world. Currently, interaction between A. cerana and N. ceranae is largely unknown. Our group previously prepared A. c. cerana workers’ midguts at 7 days post inoculation (dpi) and 10 dpi with N. ceranae spores as well as corresponding un-inoculated workers’ midguts, followed by cDNA library construction and a combination of RNAs-seq and small RNA-seq. Meanwhile, we previously prepared clean spores of N. ceranae, which were then subjected to cDNA library construction and deep sequencing. Here, based on the gained high-quality transcriptome datasets, N. ceranae differentially expressed mRNAs (DEmiRNAs) targeted by host DEmiRNAs, and A. c. cerana DEmRNAs targeted by microsporidian DEmiRNAs were deeply investigated, with a focus on targets involved in N. ceranae glycolysis/glyconeogenesis as well as virulence factors, and A. c. cerana energy metabolism and immune response. In A. c. cerana worker’s midguts at 7 (10) dpi (days post inoculation), eight (seven) up-regulated and six (two) down-regulated miRNAs were observed to target 97 (44) down-regulated and 60 (15) up-regulated N. ceranae mRNAs, respectively. Additionally, two up-regulated miRNAs (miR-60-y and miR-676-y) in host midgut at 7 dpi could target genes engaged in N. ceranae spore wall protein and glycolysis/gluconeogenesis, indicating potential host miRNA-mediated regulation of microsporidian virulence factor and energy metabolism. Meanwhile, in N. ceranae at 7 (10) dpi, 121 (110) up-regulated and 112 (104) down-regulated miRNAs were found to, respectively, target 343 (247) down-regulated and 138 (110) down-regulated mRNAs in A. c. cerana workers’ midguts. These targets in host were relevant to several crucial cellular and humoral immune pathways, such as phagasome, endocytosis, lysosomes, regulation of autophagy, and Jak–STAT signaling pathway, indicative of the involvement of N. ceranae DEmiRNAs in regulating these cellular and humoral immune pathways. In addition, N. ceranae miR-21-x was up-regulated at 7 dpi and had a target relative to oxidative phosphorylation, suggesting that miR-21-x may be used as a weapon to modulate this pivotal energy metabolism pathway. Furthermore, potential targeting relationships between two pairs of host DEmiRNAs-microsporidian DEmRNAs and two pairs of microsporidian DEmiRNAs-host DEmRNAs were validated using RT-qPCR. Our findings not only lay a foundation for exploring the molecular mechanism underlying cross-kingdom regulation between A. c. cerana workers and N. ceranae, but also offer valuable insights into Asian honey bee-microsporidian interaction.

Molecular Detection and Differentiation of Arthropod, Fungal, Protozoan, Bacterial and Viral Pathogens of Honeybees

The honeybee Apis mellifera is highly appreciated worldwide because of its products, but also as it is a pollinator of crops and wild plants. The beehive is vulnerable to infections due to arthropods, fungi, protozoa, bacteria and/or viruses that manage to by-pass the individual and social immune mechanisms of bees. Due to the close proximity of bees in the beehive and their foraging habits, infections easily spread within and between beehives. Moreover, international trade of bees has caused the global spread of infections, several of which result in significant losses for apiculture. Only in a few cases can infections be diagnosed with the naked eye, by direct observation of the pathogen in the case of some arthropods, or by pathogen-associated distinctive traits. Development of molecular methods based on the amplification and analysis of one or more genes or genomic segments has brought significant progress to the study of bee pathogens, allowing for: (i) the precise and sensitive identification of the infectious agent; (ii) the analysis of co-infections; (iii) the description of novel species; (iv) associations between geno- and pheno-types and (v) population structure studies. Sequencing of bee pathogen genomes has allowed for the identification of new molecular targets and the development of specific genotypification strategies.

Comparative Transcriptome Investigation of Nosema ceranae Infecting Eastern Honey Bee Workers

Simple Summary

At present, interaction between Nosema ceranae and Apis cerana is poorly understood, though A. cerana is the original host for N. ceranae. Here, comparative investigation was conducted using transcriptome data from N. ceranae infecting Apis cerana cerana workers at seven days post inoculation (dpi) and 10 dpi (NcT1 and NcT2 groups) as well as N. ceranae spores (NcCK group). There were 1411, 604, and 38 DEGs identified in NcCK vs. NcT1, NcCK vs. NcT2, and NcT1 vs. NcT2 comparison groups. Additionally, 10 upregulated genes and nine downregulated ones were shared by above-mentioned comparison groups. GO classification and KEGG pathway analysis suggested that these DEGs were engaged in a number of key functional terms and pathways such as cell part and glycolysis. Further analysis indicated that most of virulence factor-encoding genes were upregulated, while a few were downregulated during the fungal infection. Findings in this current work provide a basis for clarifying the molecular mechanism udnerlying N. ceranae infection and host-microsporidian interaction during bee nosemosis.

Abstract

Apis cerana is the original host for Nosema ceranae, a widespread fungal parasite resulting in honey bee nosemosis, which leads to severe losses to the apiculture industry throughout the world. However, knowledge of N. ceranae infecting eastern honey bees is extremely limited. Currently, the mechanism underlying N. ceranae infection is still largely unknown. Based on our previously gained high-quality transcriptome datasets derived from N. ceranae spores (NcCK group), N. ceranae infecting Apis cerana cerana workers at seven days post inoculation (dpi) and 10 dpi (NcT1 and NcT2 groups), comparative transcriptomic investigation was conducted in this work, with a focus on virulence factor-associated differentially expressed genes (DEGs). Microscopic observation showed that the midguts of A. c. cerana workers were effectively infected after inoculation with clean spores of N. ceranae. In total, 1411, 604, and 38 DEGs were identified from NcCK vs. NcT1, NcCK vs. NcT2, and NcT1 vs. NcT2 comparison groups. Venn analysis showed that 10 upregulated genes and nine downregulated ones were shared by the aforementioned comparison groups. The GO category indicated that these DEGs were involved in a series of functional terms relevant to biological process, cellular component, and molecular function such as metabolic process, cell part, and catalytic activity. Additionally, KEGG pathway analysis suggested that the DEGs were engaged in an array of pathways of great importance such as metabolic pathway, glycolysis, and the biosynthesis of secondary metabolites. Furthermore, expression clustering analysis demonstrated that the majority of genes encoding virulence factors such as ricin B lectins and polar tube proteins displayed apparent upregulation, whereas a few virulence factor-associated genes such as hexokinase gene and 6-phosphofructokinase gene presented downregulation during the fungal infection. Finally, the expression trend of 14 DEGs was confirmed by RT-qPCR, validating the reliability of our transcriptome datasets. These results together demonstrated that an overall alteration of the transcriptome of N. ceranae occurred during the infection of A. c. cerana workers, and most of the virulence factor-related genes were induced to activation to promote the fungal invasion. Our findings not only lay a foundation for clarifying the molecular mechanism underlying N. ceranae infection of eastern honey bee workers and microsporidian–host interaction.

Reproductive Potential Impacts Body Maintenance Parameters and Global DNA Methylation in Honeybee Workers (Apis mellifera L.)

Simple Summary

The queens and sterile workers arise from genetically identical eggs but as imagoes, they differ in their life span, DNA methylation, and their functions. In the absence of the queen, the larvae develop into rebels, i.e., workers with increased reproductive potential. We assumed that since rebels are similar to the queen in many anatomical and behavioral features, they live longer and have lower levels of global DNA methylation, even when infected, e.g., by Nosema spp. Rebels always lived longer in comparison in normal workers and unexpectedly extended longevity of normal workers when they were together, similarly as the presence of a queen did. Rebels became infected more easily but tolerated the infection better. They also had lower level of global DNA methylation than normal workers. These features expand possibilities of the use of honeybees as a model for studies on senescence, nosemosis, eusocial evolution, and epigenetics.

Abstract

The widely accepted hypothesis in life history evolution about the trade-off between fecundity and longevity is not confirmed by long-living and highly fecund queens in eusocial insects. The fact that the queens and facultatively sterile workers usually arise from genetically identical eggs but differ in DNA methylation makes them a good model for studies on senescence, eusocial evolution, and epigenetics. Therefore, honeybees seem to be especially useful here because of long living rebel-workers (RW) with high reproductive potential recently described. Longevity, ovariole number, nosema tolerance, and global DNA methylation have been assayed in normal workers (NW) versus RW in hives and cages. RW always lived longer than NW and unexpectedly extended longevity of NW when they were together, similarly as the presence of a queen did. RW lived longer despite the fact that they had higher Nosema spore load; surprisingly they became infected more easily but tolerated the infection better. Global DNA methylation increased with age, being lower in RW than in NW. Therefore, RW are queen-like considering global DNA methylation and the link between fecundity, longevity, and body maintenance. Presented features of RW expands possibilities of the use of honeybees as a model for studies on senescence, nosemosis, eusocial evolution, and epigenetics.