Comparative transcriptome investigation of Nosema ceranae infecting eastern honeybee workers

Apis cerana is the original host for Nosema ceranae, a widespread fungal parasite resulting in bee nosemosis, which leads to severe losses for apiculture industry throughout the world. However, knowledge of N. ceranae infecting eastern honeybees is extremely limited. Currently, the mechanism underlying N. ceranae infection is still largely unknown. Based on our previously gained high-quality transcriptome datasets, comparative transcriptomic investigation was conducted in this work, with a focus on virulence factor-associated differentially expressed genes (DEGs). Microscopic observation showed that A. c. cerana workers midguts were effectively infected after inoculation with clean spores of N. ceranae. Totally, 1411, 604, and 38 DEGs were identified from NcCK vs. NcT1, NcCK vs. NcT2 and NcT1 vs. NcT2 comparison groups. Venn analysis showed that ten up-regulated genes and nine down-regulated ones were shared by aforementioned comparison groups. GO category indicated 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 biosynthesis of secondary metabolites. Further, expression clustering analysis demonstrated that majority of genes encoding virulence factors such as ricin B lectins and polar tube proteins displayed apparent up-regulation, whereas a few virulence factor-associated genes such as hexokinase gene and 6-phosphofructokinase gene presented down-regulation 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. ceranae workers, and most of 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 honeybee workers, but also shed light on developing novel targets for microsporidiosis control.