Revealing the Modulatory Role of Microsporidian circRNAs in the Infection of Honey Bee Workers

Vairimorpha ceranae (formerly Nosema ceranae) is an obligate intracellular parasite that poses a major threat to the health of the honey bee. Circular RNAs (circRNAs) have been recognized as key regulators in gene expression and pathogen–host interactions. However, their expression patterns and regulatory roles in V. ceranae infection remain largely unexplored. In this study, we performed circRNA profiling in V. ceranae spores (NcCK) and the midguts of Apis mellifera ligustica workers at 7 d post inoculation (dpi) and 10 dpi (Nc7T and Nc10T) based on transcriptome sequencing, followed by in-depth investigation of the regulatory roles of differentially expressed circRNAs (DEcircRNAs). In total, 243 circRNAs were identified in V. ceranae, with lengths predominantly ranging from 201 to 400 nucleotides. Comparative analysis screened 70 and 192 DEcircRNAs in the NcCK vs. Nc7T and NcCK vs. Nc10T comparison groups, respectively, with a significant majority being downregulated. The parental genes of these DEcircRNAs were significantly enriched in fundamental cellular processes and critical pathways such as protein processing in the endoplasmic reticulum and ribosome biogenesis. Additionally, we constructed a competing endogenous RNA network, suggesting that DEcircRNAs could potentially interact with DEmiRNAs to modulate mRNAs associated with fungal proliferation-relevant signaling pathways like MAPK, PI3K–Akt, and cAMP. Moreover, numerous DEcircRNAs were predicted to contain internal ribosome entry site elements, indicative of their potential for protein coding. The back-splicing junctions and expression trends of selected DEcircRNAs were successfully validated by RT-PCR and qRT-PCR. Our data not only offer a valuable resource for future functional studies but also provide a basis for elucidating the circRNA-mediated mechanisms underlying microsporidian pathogenesis and host–pathogen interactions.