Integration of lncRNA–miRNA–mRNA reveals novel insights into oviposition regulation in honey bees

Abstract

BackgroundThe honey bee (Apis mellifera) is a highly diverse species commonly used for honey production and pollination services. The oviposition of the honey bee queen affects the development and overall performance of the colony. To investigate the ovary activation and oviposition processes on a molecular level, a genome-wide analysis of lncRNAs, miRNAs and mRNA expression in the ovaries of the queens was performed to screen for differentially expressed coding and noncoding RNAs. Further analysis identified relevant candidate genes or RNAs.ResultsThe analysis of the RNA profiles in different oviposition phase of the queens revealed that 740 lncRNAs, 81 miRNAs and 5,481 mRNAs were differently expressed during the ovary activation; 88 lncRNAs, 13 miRNAs and 338 mRNAs were differently expressed during the oviposition inhibition process; and finally, 100 lncRNAs, four miRNAs and 497 mRNAs were differently expressed during the oviposition recovery process. In addition, functional annotation of differentially expressed RNAs revealed several pathways that are closely related to oviposition, including hippo, MAPK, notch, Wnt, mTOR, TGF-beta and FoxO signaling pathways. Furthermore, in the QTL region for ovary size, 73 differentially expressed genes and 14 differentially expressed lncRNAs were located, which are considered as candidate genes affecting ovary size and oviposition. Moreover, a core set of genes served as bridges among different miRNAs were identified through the integrated analysis of lncRNA-miRNA-mRNA network.ConclusionThe observed dramatic expression changes of coding and noncoding RNAs suggest that they may play a critical role in honey bee queens’ oviposition. The identified candidate genes for oviposition activation and regulation could serve as a resource for further studies of genetic markers of oviposition in honey bees.

Population Genomics Provide Insights into the Evolution and Adaptation of the Eastern Honey Bee (Apis cerana)

“…2016 ), the TGF-beta signaling pathway, the Wnt signaling pathway, the Hedgehog signaling pathway, the Jak-STAT signaling pathway, and the Notch signaling pathway that were shown to be involved in ovary activity ( Duncan et al. 2016 ; Chen et al. 2017 ).…”

Identification of Immune Response to Sacbrood Virus Infection in Apis cerana Under Natural Condition

“…Small RNAs have been confirmed to be related to a diversity of biological processes including cell proliferation and apoptosis through the post-transcriptional regulation of gene expression (Charkhpour et al, 2011;Zhang et al, 2014;Yang et al, 2017). RNA profile analyses suggested that these sRNAs were involved in biological processes associated with development and immunity of honeybees (Chen et al, 2017). The key initiator of the RNAi pathway is dsRNA, which produced by proliferation process are processed into vsiRNA duplexes by dicer-2 (Xu and Zhou, 2017).…”

Analysis of dynamic and widespread lncRNA and miRNA expression in fetal sheep skeletal muscle

“…We estimated the PCR efficiency of each gene using standard curve calculation of the four cDNA serial dilution points. Cycle threshold (Ct) values were transformed to quantities using the comparative Ct method described by Chen et al (2017). We carried out lncRNA and mRNA data normalization using the GAPDH reference gene, and miRNA data normalization using the U6 reference gene.…”