Adult honeybee queens and workers drastically differ in ovary state and ovary size. However, this reproductive bias is only partially understood from the view of a single RNA type. In this study, we predicted 10,271 mRNAs, 7235 lncRNAs, 11,794 circRNAs, and 164 miRNAs in the ovary of honeybee workers through bioinformatics. Combining RNA sequencing data of honeybee virgin queens, 4385 mRNAs, 2390 lncRNAs, 5602 circRNAs, and 75 miRNAs were differentially expressed in workers compared with virgins. Compared with egg-laying queens, 6536 mRNAs, 3130 lncRNAs, 5751 circRNAs, and 81 miRNAs were differentially expressed in workers. Further, functional annotation revealed that neural regulation was closely related to ovary state. Moreover, the potential interactions among circRNAs, miRNAs, lncRNAs, and mRNAs revealed that vitellogenin, ecdysone-induced protein 74, ame_circ_0001176, and ame_circ_0001243 might play critical roles in the competing endogenous RNA network. These findings suggest that the integrative RNA networks have potential effects in ovarian phenotype differences in honeybees.

中文翻译：

---

蜜蜂工蜂和蜂后卵巢中编码和非编码 RNA 的全基因组表征

成年蜂王和工蜂的子房状态和子房大小差异很大。然而，从单一 RNA 类型的角度只能部分理解这种生殖偏差。在这项研究中，我们通过生物信息学预测了蜜蜂工蚁卵巢中的 10,271 个 mRNA、7235 个 lncRNA、11,794 个 circRNA 和 164 个 miRNA。结合蜜蜂处女皇后的 RNA 测序数据，与处女相比，工蜂有 4385 个 mRNA、2390 个 lncRNA、5602 个 circRNA 和 75 个 miRNA。与产卵蜂王相比，6536个mRNAs、3130个lncRNAs、5751个circRNAs和81个miRNAs在工蚁中差异表达。此外，功能注释显示神经调节与卵巢状态密切相关。此外，circRNA、miRNA、lncRNA 和 mRNA 之间的潜在相互作用表明，卵黄蛋白原、蜕皮激素诱导蛋白 74、ame_circ_0001176 和 ame_circ_0001243 可能在竞争的内源性 RNA 网络中发挥关键作用。这些发现表明，整合的 RNA 网络对蜜蜂的卵巢表型差异具有潜在影响。