Systematic identification of circular RNAs and corresponding regulatory networks unveil their potential roles in the midguts of eastern honeybee workers.

Currently, knowledge of circular RNAs (circRNAs) in insects including honeybee is extremely limited. Here, differential expression profiles and regulatory networks of circRNAs in the midguts of Apis cerana cerana workers were comprehensively investigated using transcriptome sequencing and bioinformatics. In total, 9589 circRNAs (201-800 nt in length) were identified from 8-day-old and 11-day-old workers' midguts (Ac1 and Ac2); among them, 5916 (61.70%) A. cerana cerana circRNAs showed conservation with our previously indentified circRNAs in Apis mellifera ligucstica workers' midguts (Xiong et al., Acta Entomologica Sinica 61:1363-1375, 2018). Five circRNAs were confirmed by RT-PCR and Sanger sequencing. Interestingly, novel_circ_003723, novel_circ_002714, novel_circ_002451, and novel_circ_001980 were highly expressed in both Ac1 and Ac2. In addition, the source genes of circRNAs were involved in 34 GO terms including organelle and cellular process and 141 pathways such as endocytosis and Wnt signaling pathway. Moreover, 55 DEcircRNAs including 34 upregulated and 21 downregulated circRNAs were identified in Ac2 compared with Ac1. circRNA-miRNA regulatory networks indicated that 1060 circRNAs can target 74 miRNAs; additionally, the DEcircRNA-miRNA-mRNA networks suggested that 13 downregulated circRNAs can bind to eight miRNAs and 29 miRNA-targeted mRNAs, while 16 upregulated circRNAs can link to 9 miRNAs and 29 miRNA-targeted mRNAs. These results indicated that DEcircRNAs as ceRNAs may play a comprehensive role in the growth, development, and metabolism of the worker's midgut via regulating source genes and interacting with miRNAs. Notably, eight DEcircRNAs targeting miR-6001-y were likely to be key participants in the midgut development. Our findings not only offer a valuable resource for further studies on A. cerana cerana circRNA and novel insights into understanding the molecular mechanisms underlying the midgut development of eastern honeybee but also provide putative circRNA candidates for functional research in the near future and novel biomarkers for identification of eastern honeybee species including A. cerana cerana and honeybee diseases such as chalkbrood and microsporidiosis.

Genome-Wide Identification of Circular RNAs in Fungal Parasite Nosema ceranae.

Circular RNAs (circRNAs) are newly discovered endogenous non-coding RNAs (ncRNAs) that play key roles in microRNA function and transcriptional regulation. Though a large number of circRNAs had been identified in animals and plants, however, little is known regarding circRNAs in Nosema ceranae, a widespread fungal parasite of honeybee. In this study, using deep sequencing technology and bioinformatic analysis, we predicted 204 circRNAs from N. ceranae spore samples, including 174 exonic circRNAs and 30 intergenic circRNAs. In addition, the expression of seven N. ceranae circRNAs was confirmed by RT-PCR assay. Furthermore, regulation networks of circRNAs were constructed, and 15 circRNAs were found to act as sponges of the corresponding three miRNAs. GO categorization and pathway enrichment analysis suggested that the circRNAs are likely to play significant roles in N. ceranae spore. This is the first report of circRNAs generated by a microsporidia species. Our results provide novel insights into understanding the basic biology of N. ceranae.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection.

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.