Unveiling the circRNA-Mediated Immune Responses of Western Honey Bee Larvae to Ascosphaera apis Invasion

Abstract

Western honey bee (Apis mellifera), a eusocial insect with a superior economic and ecological value, is widely used in the beekeeping industry throughout the world. As a new class of non-coding RNAs (ncRNAs), circular RNAs (circRNAs) participate in the modulation of considerable biological processes, such as the immune response via diverse manners. Here, the identification, characteristic investigation, and molecular verification of circRNAs in the Apis mellifera ligustica larval guts were conducted, and the expression pattern of larval circRNAs during the Ascosphaera apis infection was analyzed, followed by the exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 2083 circRNAs in the larval guts of A. m. ligustcia were identified, with a length distribution ranging from 106 nt to 92,798 nt. Among these, exonic circRNAs were the most abundant type and LG1 was the most distributed chromosome. Additionally, 25, 14, and 30 up-regulated circRNAs as well as 26, 25, and 62 down-regulated ones were identified in the A. apis-inoculated 4-, 5-, and 6-day-old larval guts in comparison with the corresponding un-inoculated larval guts. These DEcircRNAs were predicted to target 35, 70, and 129 source genes, which were relative to 12, 23, and 20 GO terms as well as 11, 10, and 27 KEGG pathways, including 5 cellular and humoral immune pathways containing apoptosis, autophagy, endocytosis, MAPK, Toll, and Imd signaling pathways. Furthermore, complex competing endogenous RNA (ceRNA) regulatory networks were detected to be formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The Target DEmRNAs were engaged in 24, 20, and 25 functional terms as well as 62, 80, and 159 pathways, including several vital immune defense-associated pathways, namely the lysosome, endocytosis, phagosome, autophagy, apoptosis, MAPK, Jak-STAT, Toll, and Imd signaling pathways. Finally, back-splicing sites within 15 circRNAs and the difference in the 9 DEcircRNAs’ expression between un-inoculated and A. apis-inoculated larval guts were confirmed utilizing molecular methods. These findings not only enrich our understanding of bee host-fungal pathogen interactions, but also lay a foundation for illuminating the mechanism underlying the DEcircRNA-mediated immune defense of A. m. ligustica larvae against A. apis invasion.