No abstract

piRNAs play pivotal roles in maintaining genome stability, regulating gene expression, and modulating development and immunity. However, there are few piRNA-associated studies on honey-bees, and the regulatory role of piRNAs in the development of bee guts is largely unknown. Here, the differential expression pattern of piRNAs during the developmental process of the European honey-bee (Apis mellifera) larval guts was analyzed, followed by investigation of the regulatory network and the potential function of differentially expressed piRNAs (DEpiRNAs) in regulating gut development. A total of 843 piRNAs were identified in the larval guts of A. mellifera; among these, 764 piRNAs were shared by 4- (Am4 group), 5- (Am5 group), and 6-day-old (Am6 group) larval guts, while 11, 67, and one, respectively, were unique. The first base of piRNAs in each group had a cytosine (C) bias. Additionally, 61 up-regulated and 17 down-regulated piRNAs were identified in the “Am4 vs. Am5” comparison group, further targeting 9, 983 genes, which were involved in 50 GO terms and 142 pathways, while two up-regulated and five down-regulated piRNAs were detected in the “Am5 vs. Am6” comparison group, further targeting 1, 936 genes, which were engaged in 41 functional terms and 101 pathways. piR-ame-742536 and piR-ame-856650 in the “Am4 vs. Am5” comparison group as well as piR-ame-592661 and piR-ame-31653 in the “Am5 vs. Am6” comparison group were found to link to the highest number of targets. Further analysis indicated that targets of DEpiRNAs in these two comparison groups putatively regulate seven development-associated signaling pathways, seven immune-associated pathways, and three energy metabolism pathways. Moreover, the expression trends of five randomly selected DEpiRNAs were verified based on stem-loop RT-PCR and RT-qPCR. These results were suggestive of the overall alteration of piRNAs during the larval developmental process and demonstrated that DEpiRNAs potentially modulate development-, immune-, and energy metabolism-associated pathways by regulating the expression of corresponding genes via target binding, further affecting the development of A. mellifera larval guts. Our data offer a novel insight into the development of bee larval guts and lay a basis for clarifying the underlying mechanisms.

Vairimorpha ceranae is a widespread fungal parasite of adult honey bees that leads to a serious disease called nosemosis. Circular RNAs (circRNAs) are newly discovered non-coding RNAs (ncRNAs) that regulate biological processes such as immune defense and development. Here, 8199 and 8711 circRNAs were predicted from the midguts of Apis mellifera ligustica workers at 7 d (Am7T) and 10 d (Am10T) after inoculation (dpi) with V. ceranae spores. In combination with transcriptome data from corresponding uninoculated midguts (Am7CK and Am10CK), 4464 circRNAs were found to be shared by these four groups. Additionally, 16 circRNAs were highly conserved among A. m. ligustica, Apis cerana cerana, and Homo sapiens. In the Am7CK vs. Am7T (Am10CK vs. Am10T) comparison group, 168 (306) differentially expressed circRNAs (DEcircRNAs) were identified. RT-qPCR results showed that the expression trend of eight DEcircRNAs was consistent with that in the transcriptome datasets. The source genes of DEcircRNAs in Am7CK vs. Am7T (Am10CK vs. Am10T) were engaged in 27 (35) GO functional terms, including 1 (1) immunity-associated terms. Moreover, the aforementioned source genes were involved in three cellular immune-related pathways. Moreover, 86 (178) DEcircRNAs in workers’ midguts at 7 (10) dpi could interact with 75 (103) miRNAs, further targeting 215 (305) mRNAs. These targets were associated with cellular renewal, cellular structure, carbohydrate and energy metabolism, and cellular and humoral immunity. Findings in the present study unraveled the mechanism underlying circRNA-mediated immune responses of western honey bee workers to V. ceranae invasion, but also provided new insights into host–microsporidian interaction during nosemosis.

Nosema ceranae is a widespread fungal parasite for honey bees, causing bee nosemosis. Based on deep sequencing and bioinformatics, identification of circular RNAs (circRNAs) in Apis cerana workers’ midguts and circRNA-regulated immune response of host to N. ceranae invasion were conducted in this current work, followed by molecular verification of back-splicing sites and expression trends of circRNAs. Here, 10185 and 7405 circRNAs were identified in the midguts of workers at 7 days (AcT1) and 10 days (AcT2) post inoculation days post-inoculation with N. ceranae. PCR amplification result verified the back-splicing sites within three specific circRNAs (novel_circ_005123, novel_circ_007177, and novel_circ_015140) expressed in N. ceranae-inoculated midgut. In combination with transcriptome data from corresponding un-inoculated midguts (AcCK1 and AcCK2), 2266 circRNAs were found to be shared by the aforementioned four groups, whereas the numbers of specific ones were 2618, 1917, 5691, and 3723 respectively. Further, 83 52) differentially expressed circRNAs (DEcircRNAs) were identified in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group. Source genes of DEcircRNAs in workers’ midgut at seven dpi were involved in two cellular immune-related pathways such as endocytosis and ubiquitin mediated proteolysis. Additionally, competing endogenous RNA (ceRNA) network analysis showed that 23 13) DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group could target 18 14) miRNAs and further link to 1111 (1093) mRNAs. These target mRNAs were annotated to six cellular immunity pathways including endocytosis, lysosome, phagosome, ubiquitin mediated proteolysis, metabolism of xenobiotics by cytochrome P450, and insect hormone biosynthesis. Moreover, 284 164) internal ribosome entry site and 54 26) ORFs were identified from DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group; additionally, ORFs in DEcircRNAs in midgut at seven dpi with N. ceranae were associated with several cellular immune pathways including endocytosis and ubiquitin-mediated proteolysis. Ultimately, RT-qPCR results showed that the expression trends of six DEcircRNAs were consistent with those in transcriptome data. These results demonstrated that N. ceranae altered the expression pattern of circRNAs in A. c. cerana workers’ midguts, and DEcircRNAs were likely to regulate host cellular and humoral immune response to microsporidian infection. Our findings lay a foundation for clarifying the mechanism underlying host immune response to N. ceranae infection and provide a new insight into interaction between Asian honey bee and microsporidian.

No abstract