Sequence and Expression Characteristics of Long Noncoding RNAs in Honey Bee Caste Development – Potential Novel Regulators for Transgressive Ovary Size

Division of labor in social insect colonies relies on a strong reproductive bias that favors queens. Although the ecological and evolutionary success attained through caste systems is well sketched out in terms of ultimate causes, the molecular and cellular underpinnings driving the development of caste phenotypes are still far from understood. Recent genomics approaches on honey bee developmental biology revealed a set of genes that are differentially expressed genes in larval ovaries and associated with transgressive ovary size in queens and massive cell death in workers. Amongst these, two contigs called special attention, both being over 200 bp in size and lacking apparent coding potential. Herein, we obtained their full cDNA sequences. These and their secondary structure characteristics placed in evidence that they are bona fide long noncoding RNAs (lncRNA) differentially expressed in larval ovaries, thus named lncov1 and lncov2. Genomically, both map within a previously identified QTL on chromosome 11, associated with transgressive ovary size in honey bee workers. As lncov1 was over-expressed in worker ovaries we focused on this gene. Real-time qPCR analysis on larval worker ovaries evidenced an expression peak coinciding with the onset of autophagic cell death. Cellular localization analysis through fluorescence in situ hybridization revealed perinuclear spots resembling omega speckles known to regulate trafficking of RNA-binding proteins. With only four lncRNAs known so far in honey bees, two expressed in the ovaries, these findings open a novel perspective on regulatory factors acting in the fine tuning of developmental processes underlying phenotypic plasticity related to social life histories.