Abstract:ABSTRACT Nosema ceranae is a widespread fungal parasite for adult honeybees, severely damaging bee health and sustainable development of apiculture. Circular RNAs (circRNAs) are a class of newly discovered noncoding RNAs (ncRNAs) that regulate a number of biological processes such as immune defense and development. In this current work, based on previously obtained whole transcriptome data, 8 199 and 8 711 circRNAs were predicted from the midguts of Apis mellifera ligustica workers at 7 days (AmT1) and 10 days (AmT2) post inoculation (dpi) with N. ceranae using bioinformatics. Additionally, in combination with transcriptome data from uninfected midguts (AmCK1 and AmCK2) (Xiong et al., 2018), 4 464 circRNAs were found to be shared by the aforementioned four groups, whereas the numbers of specifically transcribed circRNAs in each group were 1 389, 1 696, 1 019 and 1 871, respectively. Furthermore, 10 226 circRNAs were homologous to Apis cerana cerana circRNAs, while 20 circRNAs had homology with Homo sapiens circRNAs; in addition, 16 circRNAs were highly conserved in these three species. Differential expression analysis showed that 168 (306) differentially expressed circRNAs (DEcircRNAs) were identified in AmCK1 vs AmT1 (AmCK2 vs AmT2) comparison group, including 61 (143) upregulated circRNAs and 107 (163) downregulated circRNAs. Moreover, RT-qPCR results showed that the expression trend of eight DEcircRNAs was consistent with that of the transcriptome dataset. Based on GO database annotation, we observed that source genes of DEcircRNAs in AmCK1 vs AmT1 (AmCK2 vs AmT2) were engaged in 27 (35) functional terms, including two (two) cell renewal-associated terms, seven (seven) cell structure-associated terms, and one (one) immunity-associated terms. Additionally, DEcircRNA source genes in AmCK1 vs AmT1 were involved in two cell renewal-related pathways, Hippo and Wnt signaling pathways, and three carbohydrate metabolism-related pathways, galactose metabolism, starch and sucrose metabolism, fructose and mannose metabolism, only one energy metabolism-related pathway (oxidative phosphorylation pathway), three cellular immune-related pathways, endocytosis, phagosome, and lysosome, and a humoral immune-related pathway (FoxO signaling pathway). In AmCK2 vs AmT2 comparison group, more source genes of DEcircRNAs were associated with the abovementioned pathways relative to cell renewal, carbohydrate metabolism, and cellular and humoral immune pathways. In addition, 122 (234) DEcircRNAs in the host midgut at 7 dpi (10 dpi) with N. ceranae targeted 82 (106) miRNAs. Furthermore, 75 (103) miRNAs targeted by 86 (178) DEcircRNAs in AmCK1 vs AmT1 (AmCK2 vs AmT2) further bound to 215 (305) mRNAs. These targets could be annotated as an array of functional terms and pathways related to cellular renewal, cellular structure, carbohydrate and energy metabolism, and cellular and humoral immunity. In a word, we for the first time explored immune responses mediated by DEcircRNAs in the midguts of A. m. ligustica workers to N. ceranae infection. Our data provide a foundation for clarifying the molecular mechanism underlying immune response of western honeybee to N. ceranae invasion, but also a new insight into further understanding the host-pathogen interaction during bee microsporidiosis.

Differences between uncapping and removal behaviors in Apis cerana from the perspective of long non-coding RNAs

Recapping Behavior in Apis Cerana : Does It Contribute to Resistance Against Varroa Spp.?

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis Mellifera Ligustica Response to Nosema Ceranae Infection

Comparison of Brain Gene Expression Profiles Associated with Auto-Grooming Behavior between Apis cerana and Apis mellifera Infested by Varroa destructor

Systematic Characterization and Regulatory Role of lncRNAs in Asian Honey Bees Responding to Microsporidian Infestation

Genomic and transcriptomic analysis of the Asian honeybee Apis cerana provides novel insights into honeybee biology

Diverse Regulatory Manners and Potential Roles of lncRNAs in the Developmental Process of Asian Honey Bee (Apis cerana) Larval Guts

Regulatory Roles of Long Non-Coding RNAs Relevant to Antioxidant Enzymes and Immune Responses of Apis cerana Larvae Following Ascosphaera apis Invasion

Deciphering the Mechanism Underlying Circrna-Mediated Immune Responses of Western Honeybees to Nosema Ceranae Infection

Go East for Better Honey Bee Health: Apis cerana Is Faster at Hygienic Behavior than A. mellifera

Dynamics and regulatory role of circRNAs in Asian honey bee larvae following fungal infection

Systematic Identification of Circular RNAs and Corresponding Regulatory Networks Unveil Their Potential Roles in the Midguts of Eastern Honeybee Workers

Study on microRNAs and their regulatory roles in ontogenesis,castes differentiation and labor division in honeybee(Apis)

Extent and complexity of RNA processing in honey bee queen and worker caste development

Possible interactions between gut microbiome and division of labor in honey bees

The Neuroproteomic Basis of Enhanced Perception and Processing of Brood Signals That Trigger Increased Reproductive Investment in Honeybee (Apis mellifera) Workers

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

Regulatory roles of long non-coding RNAs in short-term heat stress in adult worker bees

Comparative identification of microRNAs in Apis cerana cerana workers'midguts responding to Nosema ceranae invasion

Extent and complexity of RNA processing in the development of honey bee queen and worker castes revealed by Nanopore direct RNA sequencing

Morphological and functional characterization of honey bee, Apis mellifera , hemocyte cell communities