RESUMEN
Circular RNAs (circRNAs) are a class of novel non-coding RNAs (ncRNAs) that play essential roles in the development and growth of vertebrates through multiple manners. However, the mechanism by which circRNAs modulate the honey bee gut development is currently poorly understood. Utilizing the transcriptome data we obtained earlier, the highly expressed circRNAs in the Apis mellifera worker 4-, 5-, and 6-day-old larval guts were analyzed, which was followed by an in-depth investigation of the expression pattern of circRNAs during the process of larval guts development and the potential regulatory roles of differentially expressed circRNAs (DEcircRNAs). In total, 1728 expressed circRNAs were detected in the A. mellifera larval guts. Among the most highly expressed 10 circRNAs, seven (novel_circ_000069, novel_circ_000027, novel_circ_000438, etc.) were shared by the 4-, 5-, and 6-day-old larval guts. In addition, 21 (46) up-regulated and 22 (27) down-regulated circRNAs were, respectively, screened in the Am4 vs. Am5 (Am5 vs. Am6) comparison groups. Additionally, nine DEcircRNAs, such as novel_circ_000340, novel_circ_000758 and novel_circ_001116, were shared by these two comparison groups. These DEcircRNAs were predicted to be transcribed from 14 and 29 parental genes; these were respectively annotated to 15 and 22 GO terms such as biological regulation and catalytic activity as well as 16 and 21 KEGG pathways such as dorsoventral axis formation and apoptosis. Moreover, a complicated competing endogenous RNA (ceRNA) network was observed; novel_circ_000838 in the Am4 vs. Am5 comparison group potentially targeted ame-miR-6000a-3p, further targeting 518 mRNAs engaged in several developmental signaling pathways (e.g., TGF-beta, hedgehog, and wnt signaling pathway) and immune pathways (e.g., phagosome, lysosome, and MAPK signaling pathway). The results demonstrated that the novel_circ_000838-ame-miR-6000a-3p axis may plays a critical regulatory part in the larval gut development and immunity. Furthermore, back-splicing sites of six randomly selected DEcircRNAs were amplified and verified by PCR; an RT-qPCR assay of these six DEcircRNAs confirmed the reliability of the used high-throughput sequencing data. Our findings provide a novel insight into the honey bee gut development and pave a way for illustration of the circRNA-modulated developmental mechanisms underlying the A. mellifera worker larval guts.
RESUMEN
MiRNAs, as a kind of key regulators in gene expression, play vital roles in numerous life activities from cellular proliferation and differentiation to development and immunity. However, little is known about the regulatory manner of miRNAs in the development of Asian honey bee (Apis cerana) guts. Here, on basis of our previously gained high-quality transcriptome data, transcriptome-wide identification of miRNAs in the larval guts of Apis cerana cerana was conducted, followed by investigation of the miRNAs' differential expression profile during the gut development. In addition to the regulatory network, the potential function of differentially expressed miRNAs (DEmiRNAs) was further analyzed. In total, 330, 351, and 321 miRNAs were identified in the 4-, 5-, and 6-day-old larval guts, respectively; among these, 257 miRNAs were shared, while 38, 51, and 36 ones were specifically expressed. Sequences of six miRNAs were confirmed by stem-loop RT-PCR and Sanger sequencing. Additionally, in the "Ac4 vs. Ac5" comparison group, there were seven up-regulated and eight down-regulated miRNAs; these DEmiRNAs could target 5041 mRNAs, involving a series of GO terms and KEGG pathways associated with growth and development, such as cellular process, cell part, Wnt, and Hippo. Comparatively, four up-regulated and six down-regulated miRNAs detected in the "Ac5 vs. Ac6" comparison group and the targets were associated with diverse development-related terms and pathways, including cell, organelle, Notch and Wnt. Intriguingly, it was noticed that miR-6001-y presented a continuous up-regulation trend across the developmental process of larval guts, implying that miR-6001-y may be a potential essential modulator in the development process of larval guts. Further investigation indicated that 43 targets in the "Ac4 vs. Ac5" comparison group and 31 targets in the "Ac5 vs. Ac6" comparison group were engaged in several crucial development-associated signaling pathways such as Wnt, Hippo, and Notch. Ultimately, the expression trends of five randomly selected DEmiRNAs were verified using RT-qPCR. These results demonstrated that dynamic expression and structural alteration of miRNAs were accompanied by the development of A. c. cerana larval guts, and DEmiRNAs were likely to participate in the modulation of growth as well as development of larval guts by affecting several critical pathways via regulation of the expression of target genes. Our data offer a basis for elucidating the developmental mechanism underlying Asian honey bee larval guts.
RESUMEN
Asian honey bee Apis cerana is the original host for Nosema ceranae, a unicellular fungal parasite that causes bee nosemosis throughout the world. Currently, interaction between A. cerana and N. ceranae is largely unknown. Our group previously prepared A. c. cerana workers' midguts at 7 days post inoculation (dpi) and 10 dpi with N. ceranae spores as well as corresponding un-inoculated workers' midguts, followed by cDNA library construction and a combination of RNAs-seq and small RNA-seq. Meanwhile, we previously prepared clean spores of N. ceranae, which were then subjected to cDNA library construction and deep sequencing. Here, based on the gained high-quality transcriptome datasets, N. ceranae differentially expressed mRNAs (DEmiRNAs) targeted by host DEmiRNAs, and A. c. cerana DEmRNAs targeted by microsporidian DEmiRNAs were deeply investigated, with a focus on targets involved in N. ceranae glycolysis/glyconeogenesis as well as virulence factors, and A. c. cerana energy metabolism and immune response. In A. c. cerana worker's midguts at 7 (10) dpi (days post inoculation), eight (seven) up-regulated and six (two) down-regulated miRNAs were observed to target 97 (44) down-regulated and 60 (15) up-regulated N. ceranae mRNAs, respectively. Additionally, two up-regulated miRNAs (miR-60-y and miR-676-y) in host midgut at 7 dpi could target genes engaged in N. ceranae spore wall protein and glycolysis/gluconeogenesis, indicating potential host miRNA-mediated regulation of microsporidian virulence factor and energy metabolism. Meanwhile, in N. ceranae at 7 (10) dpi, 121 (110) up-regulated and 112 (104) down-regulated miRNAs were found to, respectively, target 343 (247) down-regulated and 138 (110) down-regulated mRNAs in A. c. cerana workers' midguts. These targets in host were relevant to several crucial cellular and humoral immune pathways, such as phagasome, endocytosis, lysosomes, regulation of autophagy, and Jak-STAT signaling pathway, indicative of the involvement of N. ceranae DEmiRNAs in regulating these cellular and humoral immune pathways. In addition, N. ceranae miR-21-x was up-regulated at 7 dpi and had a target relative to oxidative phosphorylation, suggesting that miR-21-x may be used as a weapon to modulate this pivotal energy metabolism pathway. Furthermore, potential targeting relationships between two pairs of host DEmiRNAs-microsporidian DEmRNAs and two pairs of microsporidian DEmiRNAs-host DEmRNAs were validated using RT-qPCR. Our findings not only lay a foundation for exploring the molecular mechanism underlying cross-kingdom regulation between A. c. cerana workers and N. ceranae, but also offer valuable insights into Asian honey bee-microsporidian interaction.