Sequence and structural properties of circular RNAs in the brain of nurse and forager honeybees (Apis mellifera)

Abstract

BackgroundThe honeybee (Apis mellifera) represents a model organism for social insects displaying behavioral plasticity. This is reflected by an age-dependent task allocation. The most protruding tasks are performed by young nurse bees and older forager bees that take care of the brood inside the hive and collect food from outside the hive, respectively. The molecular mechanism leading to the transition from nurse bees to foragers is currently under intense research. Circular RNAs, however, were not considered in this context so far. As of today, this group of non-coding RNAs was only known to exist in two other insects, Drosophila melanogaster and Bombyx mori. Here we complement the state of circular RNA research with the first characterization in a social insect.ResultsWe identified numerous circular RNAs in the brain of A. mellifera nurse bees and forager bees using RNA-Seq with exonuclease enrichment. Presence and circularity were verified for the most abundant representatives. Back-splicing in honeybee occurs further towards the end of transcripts and in transcripts with a high number of exons. The occurrence of circularized exons is correlated with length and CpG-content of their flanking introns. The latter coincides with increased DNA-methylation in the respective loci. For two prominent circular RNAs the abundance in worker bee brains was quantified in TaqMan assays. In line with previous findings of circular RNAs in Drosophila, circAmrsmep2 accumulates with increasing age of the insect. In contrast, the levels of circAmrad appear age-independent and correlate with the bee’s task. Its parental gene is related to amnesia-resistant memory.ConclusionsWe provide the first characterization of circRNAs in a social insect. Many of the RNAs identified here show homologies to circular RNAs found in Drosophila and Bombyx, indicating that circular RNAs are a common feature among insects. We find that exon circularization is correlated to DNA-methylation at the flanking introns. The levels of circAmrad suggest a task-dependent abundance that is decoupled from age. Moreover, a GO term analysis shows an enrichment of task-related functions. We conclude that circular RNAs could be relevant for task allocation in honeybee and should be investigated further in this context.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-5402-6) contains supplementary material, which is available to authorized users.

Single-cell transcriptomic analysis of honeybee brains identifies vitellogenin as caste differentiation-related factor

“…Previous studies investigating neuronal cells in the bee brain suggested that there are remarkable differences in the transcriptomes ( Behura and Whitfield, 2010 ; Liu et al., 2011 , 2012 ; Tholken et al., 2019 ), proteomes ( Hernandez et al., 2012 ), and methylomes ( Foret et al., 2012 ; Herb et al., 2012 ; Lyko et al., 2010 ) between workers and queens. To better understand the difference between queens and workers, we established high-quality single-cell transcriptomic atlases of female honeybee brains.…”

Genome-wide identification and functional analysis of circRNAs in Trichophyton rubrum conidial and mycelial stages

“…These results illustrate the complexity of gene expression regulation. In addition, transcription factors (TFs) have been found to participate in this regulatory process [ 51 ]. In our previous study, the miRNA Tru-miR5658 was suggested to bind to three TFs, including the C6 transcription factor, bZIP transcription factor and transcription factor TFIIIB component [ 52 ].…”

Deciphering the CircRNA-Regulated Response of Western Honey Bee (Apis mellifera) Workers to Microsporidian Invasion

“…In honey bees, only three circRNA-associated studies have been reported [ 28 , 29 , 39 ]. In a previous work, we conducted deep sequencing of the midguts of A. m. ligustica workers at 7 dpi and 10 dpi with V. ceranae and corresponding uninoculated midguts using strand-specific cDNA library-based RNA-seq and identified 10,833 novel circRNAs based on transcriptome data from control groups; it was found that these circRNAs were 15–1000 nt in length and the main type of circularization was annotated exon circRNA [ 40 ].…”