Genome-wide characterization of coding and non-coding RNAs in the ovary of honeybee workers and queens

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.