Uncovering the immune responses of Apis mellifera ligustica larval gut to Ascosphaera apis infection utilizing transcriptome sequencing

- This is the first study on transcriptome analysis of immune responses of the bee larval gut to Ascosphaera apis using RNA sequencing.
- Cellular immunity of the larval gut was extremely inhibited by Ascosphaera apis, which may lay some foundation for the outbreak of chalkbrood.
- Activation of the humoral immunity in the western honeybee larval gut was induced by Ascosphaera apis.
- No component of the Jak-STAT signaling pathway was indentified in the challenged larval gut of Apis mellifera ligustica.

Honeybees are susceptible to a variety of diseases, including chalkbrood, which is capable of causing huge losses of both the number of bees and colony productivity. This research is designed to characterize the transcriptome profiles of Ascosphaera apis-treated and un-treated larval guts of Apis mellifera ligustica in an attempt to unravel the molecular mechanism underlying the immune responses of western honeybee larval guts to mycosis. In this study, 24, 296 and 2157 genes were observed to be differentially expressed in A. apis-treated Apis mellifera (4-, 5- and 6-day-old) compared with un-treated larval guts. Moreover, the expression patterns of differentially expressed genes (DEGs) were examined via trend analysis, and subsequently, gene ontology analysis and KEGG pathway enrichment analysis were conducted for DEGs involved in up- and down-regulated profiles. Immunity-related pathways were selected for further analysis, and our results demonstrated that a total of 13 and 50 DEGs were annotated in the humoral immune-related and cellular immune-related pathways, respectively. Additionally, we observed that many DEGs up-regulated in treated guts were part of cellular immune pathways, such as the lysosome, ubiquitin mediated proteolysis, and insect hormone biosynthesis pathways and were induced by A. apis invasion. However, more down-regulated DEGs were restrained. Surprisingly, a majority of DEGs within the Toll-like receptor signaling pathway, and the MAPK signaling pathway were up-regulated in treated guts, while all but two genes involved in the NF-κB signaling pathway were down-regulated, which suggested that most genes involved in humoral immune-related pathways were activated in response to the invasive fungal pathogen. This study's findings provide valuable information regarding the investigation of the molecular mechanism of immunity defenses of A. m. ligustica larval guts to infection with A. apis. Furthermore, these studies lay the groundwork for future researches on key genes controlling the susceptibility of A. m. ligustica larvae to chalkbrood.