Identification of circular RNA in the Bombyx mori silk gland

Dynamics and regulatory role of circRNAs in Asian honey bee larvae following fungal infection

Non-coding RNA (ncRNA) plays a vital part in the regulation of immune responses, growth, and development in plants and animals. Here, the identification, characteristic analysis, and molecular verification of circRNAs in Apis cerana cerana worker larval guts were conducted, followed by in-depth investigation of the expression pattern of larval circRNAs during Ascosphaera apis infection and exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 3178 circRNAs in the larval guts of A. c. cerana were identified, with a length distribution ranging from 15 to 96,007 nt. Additionally, 155, 95, and 86 DEcircRNAs were identified in the in the 4-, 5-, and 6-day-old larval guts following A. apis infection. These DEcircRNAs were predicted to target 29, 25, and 18 parental genes relevant to 12, 20, and 17 GO terms as well as 144, 114, and 61 KEGG pathways, including 5 cellular and 4 humoral immune pathways. Complex competing endogenous RNA (ceRNA) regulatory networks were detected as being formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The target DEmRNAs were engaged in 36, 47, and 47 GO terms as well as 331, 332, and 331 pathways, including 6 cellular and 6 humoral immune pathways. Further, 19 DEcircRNAs, 5 DEmiRNAs, and 3 mRNAs were included in the sub-networks relative to 3 antioxidant enzymes. Finally, back-splicing sites within 15 circRNAs and the difference in the 15 DEcircRNAs' expression between uninoculated and A. apis-inoculated larval guts were confirmed based on molecular methods. These findings not only enrich our understanding of bee host-fungal pathogen interactions but also lay a foundation for illuminating the mechanism underlying the DEcircRNA-mediated immune defense of A. c. cerana larvae against A. apis invasion. KEY POINTS: • The expression pattern of circRNAs was altered in the A. cerana worker larval guts following A. apis infection. • Back-splicing sites within 15 A. cerana circRNAs were verified using molecular approaches. DEcircRNAs potentially modulated immune responses and antioxidant enzymes in A. apis-challenged host guts.

MicroRNA750-3p Targets Processing of Precursor 7 to Suppress Rice Black-Streaked Dwarf Virus Propagation in Vector Laodelphax striatellus

MicroRNAs (miRNAs) are non-coding RNAs, which, as members of the RNA interference pathway, play a pivotal role in antiviral infection. Almost 80% of plant viruses are transmitted by insect vectors; however, little is known about the interaction of the miRNAs of insect vectors with plant viruses. Here, we took rice black-streaked dwarf virus (RBSDV), a devastating virus to rice production in eastern Asia, and the small brown planthopper, (SBPH, Laodelphax striatellus) as a model to investigate the role of microRNA750-3p (miR750-3p) in regulating viral transmission. Our results showed that Ls-miR750-3p was downregulated in RBSDV-infected SBPH and predominately expressed in the midgut of SBPH. Injection with miR750-3p agomir significantly reduced viral accumulation, and the injection with the miR750-3p inhibitor, antagomir-750-3p, dramatically promoted the viral accumulation in SBPH, as detected using RT-qPCR and Western blotting. The processing of precursor 7 (POP7), a subunit of RNase P and RNase MRP, was screened, identified, and verified using a dual luciferase reporter assay as one target of miR750-3p. Knockdown of POP7 notably increased RBSDV viral propagation in SBPH and then increased the viral transmission rate by SBPH. Taken together, our data indicate that miR750-3p targets POP7 to suppress RBSDV infection in its insect vector. These results enriched the role of POP7 in modulating virus infection in host insects and shared new insight into the function of miRNAs in plant virus and insect vector interaction.

Novel Insights into the circRNA-Modulated Developmental Mechanism of Western Honey Bee Larval Guts

Circular RNAs (circRNAs) are a class of novel non-coding RNAs (ncRNAs) that play essential roles in the development and growth of vertebrates through multiple manners. However, the mechanism by which circRNAs modulate the honey bee gut development is currently poorly understood. Utilizing the transcriptome data we obtained earlier, the highly expressed circRNAs in the Apis mellifera worker 4-, 5-, and 6-day-old larval guts were analyzed, which was followed by an in-depth investigation of the expression pattern of circRNAs during the process of larval guts development and the potential regulatory roles of differentially expressed circRNAs (DEcircRNAs). In total, 1728 expressed circRNAs were detected in the A. mellifera larval guts. Among the most highly expressed 10 circRNAs, seven (novel_circ_000069, novel_circ_000027, novel_circ_000438, etc.) were shared by the 4-, 5-, and 6-day-old larval guts. In addition, 21 (46) up-regulated and 22 (27) down-regulated circRNAs were, respectively, screened in the Am4 vs. Am5 (Am5 vs. Am6) comparison groups. Additionally, nine DEcircRNAs, such as novel_circ_000340, novel_circ_000758 and novel_circ_001116, were shared by these two comparison groups. These DEcircRNAs were predicted to be transcribed from 14 and 29 parental genes; these were respectively annotated to 15 and 22 GO terms such as biological regulation and catalytic activity as well as 16 and 21 KEGG pathways such as dorsoventral axis formation and apoptosis. Moreover, a complicated competing endogenous RNA (ceRNA) network was observed; novel_circ_000838 in the Am4 vs. Am5 comparison group potentially targeted ame-miR-6000a-3p, further targeting 518 mRNAs engaged in several developmental signaling pathways (e.g., TGF-beta, hedgehog, and wnt signaling pathway) and immune pathways (e.g., phagosome, lysosome, and MAPK signaling pathway). The results demonstrated that the novel_circ_000838-ame-miR-6000a-3p axis may plays a critical regulatory part in the larval gut development and immunity. Furthermore, back-splicing sites of six randomly selected DEcircRNAs were amplified and verified by PCR; an RT-qPCR assay of these six DEcircRNAs confirmed the reliability of the used high-throughput sequencing data. Our findings provide a novel insight into the honey bee gut development and pave a way for illustration of the circRNA-modulated developmental mechanisms underlying the A. mellifera worker larval guts.

Genome-wide analysis of long noncoding RNAs and their association in regulating the metamorphosis of the Sarcophaga peregrina (Diptera: Sarcophagidae)

BACKGROUND

The flesh fly, Sarcophaga peregrina (Diptera: Sarcophagidae), is an important hygiene pest, that causes myiasis in humans and other mammals, typically livestock, and as a vector for various parasitic agents, including bacteria, viruses, and parasites. The role of long non-coding RNAs (lncRNAs) in regulating gene expression during metamorphosis of the flesh fly has not been well established.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we performed genome-wide identification and characterization of lncRNAs from the early pupal stage (1-days pupae), mid-term pupal stage (5-days pupae), and late pupal stage (9-days pupae) of S. peregrina by RNA-seq, and a total of 6921 lncRNAs transcripts were identified. RT-qPCR and enrichment analyses revealed the differentially expressed lncRNAs (DE lncRNAs) that might be associated with insect metamorphosis development. Furthermore, functional analysis revealed that the DE lncRNA (SP_lnc5000) could potentially be involved in regulating the metamorphosis of S. peregrina. RNA interference of SP_lnc5000 caused reduced expression of metamorphosis-related genes in 20-hydroxyecdysone (20E) signaling (Br-c, Ftz-F1), cuticle tanning pathway (TH, DOPA), and chitin related pathway (Cht5). Injection of dsSP_lnc5000 in 3rd instar larvae of S. peregrina resulted in deformed pupae, stagnation of pupal-adult metamorphosis, and a decrease in development time of pupal, pupariation rates and eclosion rates. Hematoxylin-eosin staining (H&E), scanning electron microscope (SEM) observation and cuticle hydrocarbons (CHCs) analysis indicated that SP_lnc5000 had crucial roles in the metamorphosis developmental by modulating pupal cuticular development.

CONCLUSIONS/SIGNIFICANCE

We established that the lncRNA SP_lnc5000 potentially regulates the metamorphosis of S. peregrina by putatively affecting the structure and composition of the pupal cuticle. This study enhances our understanding of lncRNAs as regulators of metamorphosis in S. peregrina, and provide valuable insights into the identification of potential targets for vector control and the development of effective strategies for controlling the spread of myiasis and parasitic diseases.

Multi-omics study and ncRNA regulation of anti-BmNPV in silkworms, Bombyx mori: an update

Bombyx mori silkworm is an important economic insect which has a significant contribution to the improvement of the economy. Bombyx mori nucleopolyhedrovirus (BmNPV) is a vitally significant purulent virus that impedes the sustainable and stable development of the silkworm industry, resulting in substantial economic losses. In recent years, with the development of biotechnology, transcriptomics, proteomics, metabolomics, and the related techniques have been used to select BmNPV-resistant genes, proteins, and metabolites. The regulatory networks between viruses and hosts have been gradually clarified with the discovery of ncRNAs, such as miRNA, lncRNA, and circRNA in cells. Thus, this paper aims to highlight the results of current multi-omics and ncRNA studies on BmNPV resistance in the silkworm, providing some references for resistant strategies in the silkworm to BmNPV.

Effect of fluvalinate on the expression profile of circular RNA in brain tissue of Apis mellifera ligustica workers

Fluvalinate is widely used in apiculture as an acaricide for removing Varroa mites, but there have been growing concerns about the negative effects of fluvalinate on honeybees in recent years. Previous research revealed changes in the miRNA and mRNA expression profiles of Apis mellifera ligustica brain tissues during fluvalinate exposure, as well as key genes and pathways. The role of circRNAs in this process, however, is unknown. The goal of this study was to discover the fluvalinate-induced changes in circular RNA (circRNA) expression profiles of brain tissue of A. mellifera ligustica workers. A total of 10,780 circRNAs were detected in A. mellifera ligustica brain tissue, of which eight were differentially expressed between at least two of the four time periods before and after fluvalinate administration, and six circRNAs were experimentally verified to be structurally correct, and their expression patterns were consistent with transcriptome sequencing results. Furthermore, ceRNA analysis revealed that five differentially expressed circRNAs (DECs) (novel_circ_012139, novel_circ_011690, novel_circ_002628, novel_circ_004765, and novel_circ_010008) were primarily involved in apoptosis-related functions by competitive binding with miRNAs. This study discovered changes in the circRNA expression profile of A. mellifera ligustica brain tissue caused by fluvalinate exposure, and it provides a useful reference for the biological function study of circRNAs in A. mellifera ligustica.

Aal-circRNA-407 regulates ovarian development of Aedes albopictus, a major arbovirus vector, via the miR-9a-5p/Foxl axis

Aedes albopictus shows a rapid global expansion and dramatic vectorial capacity for various arboviruses, thus posing a severe threat to global health. Although many noncoding RNAs have been confirmed to play functional roles in various biological processes in Ae. albopictus, the roles of circRNA remain a mystery. In the present study, we first performed high-throughput circRNA sequencing in Ae. albopictus. Then, we identified a cysteine desulfurase (CsdA) superfamily gene-originated circRNA, named aal-circRNA-407, which was the third most abundant circRNA in adult females and displayed a fat body highly expressed manifestation and blood feeding-dependent onset. SiRNA-mediated knockdown of circRNA-407 resulted in a decrease in the number of developing follicles and a reduction in follicle size post blood meal. Furthermore, we demonstrated that circRNA-407 can act as a sponge of aal-miR-9a-5p to promote the expression of its target gene Foxl and eventually regulate ovarian development. Our study is the first to report a functional circRNA in mosquitoes, expanding our current understanding of important biological roles in mosquitoes and providing an alternative genetic strategy for mosquito control.

CircRNA-regulated immune responses of asian honey bee workers to microsporidian infection

Nosema ceranae is a widespread fungal parasite for honey bees, causing bee nosemosis. Based on deep sequencing and bioinformatics, identification of circular RNAs (circRNAs) in Apis cerana workers’ midguts and circRNA-regulated immune response of host to N. ceranae invasion were conducted in this current work, followed by molecular verification of back-splicing sites and expression trends of circRNAs. Here, 10185 and 7405 circRNAs were identified in the midguts of workers at 7 days (AcT1) and 10 days (AcT2) post inoculation days post-inoculation with N. ceranae. PCR amplification result verified the back-splicing sites within three specific circRNAs (novel_circ_005123, novel_circ_007177, and novel_circ_015140) expressed in N. ceranae-inoculated midgut. In combination with transcriptome data from corresponding un-inoculated midguts (AcCK1 and AcCK2), 2266 circRNAs were found to be shared by the aforementioned four groups, whereas the numbers of specific ones were 2618, 1917, 5691, and 3723 respectively. Further, 83 52) differentially expressed circRNAs (DEcircRNAs) were identified in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group. Source genes of DEcircRNAs in workers’ midgut at seven dpi were involved in two cellular immune-related pathways such as endocytosis and ubiquitin mediated proteolysis. Additionally, competing endogenous RNA (ceRNA) network analysis showed that 23 13) DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group could target 18 14) miRNAs and further link to 1111 (1093) mRNAs. These target mRNAs were annotated to six cellular immunity pathways including endocytosis, lysosome, phagosome, ubiquitin mediated proteolysis, metabolism of xenobiotics by cytochrome P450, and insect hormone biosynthesis. Moreover, 284 164) internal ribosome entry site and 54 26) ORFs were identified from DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group; additionally, ORFs in DEcircRNAs in midgut at seven dpi with N. ceranae were associated with several cellular immune pathways including endocytosis and ubiquitin-mediated proteolysis. Ultimately, RT-qPCR results showed that the expression trends of six DEcircRNAs were consistent with those in transcriptome data. These results demonstrated that N. ceranae altered the expression pattern of circRNAs in A. c. cerana workers’ midguts, and DEcircRNAs were likely to regulate host cellular and humoral immune response to microsporidian infection. Our findings lay a foundation for clarifying the mechanism underlying host immune response to N. ceranae infection and provide a new insight into interaction between Asian honey bee and microsporidian.

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

Vairimorpha ceranae is a widespread fungal parasite of adult honey bees that leads to a serious disease called nosemosis. Circular RNAs (circRNAs) are newly discovered non-coding RNAs (ncRNAs) that regulate biological processes such as immune defense and development. Here, 8199 and 8711 circRNAs were predicted from the midguts of Apis mellifera ligustica workers at 7 d (Am7T) and 10 d (Am10T) after inoculation (dpi) with V. ceranae spores. In combination with transcriptome data from corresponding uninoculated midguts (Am7CK and Am10CK), 4464 circRNAs were found to be shared by these four groups. Additionally, 16 circRNAs were highly conserved among A. m. ligustica, Apis cerana cerana, and Homo sapiens. In the Am7CK vs. Am7T (Am10CK vs. Am10T) comparison group, 168 (306) differentially expressed circRNAs (DEcircRNAs) were identified. RT-qPCR results showed that the expression trend of eight DEcircRNAs was consistent with that in the transcriptome datasets. The source genes of DEcircRNAs in Am7CK vs. Am7T (Am10CK vs. Am10T) were engaged in 27 (35) GO functional terms, including 1 (1) immunity-associated terms. Moreover, the aforementioned source genes were involved in three cellular immune-related pathways. Moreover, 86 (178) DEcircRNAs in workers' midguts at 7 (10) dpi could interact with 75 (103) miRNAs, further targeting 215 (305) mRNAs. These targets were associated with cellular renewal, cellular structure, carbohydrate and energy metabolism, and cellular and humoral immunity. Findings in the present study unraveled the mechanism underlying circRNA-mediated immune responses of western honey bee workers to V. ceranae invasion, but also provided new insights into host-microsporidian interaction during nosemosis.

Circular RNA, circ1-3p, is Involved in Cyflumetofen Resistance by Acting as a Competitive RNA against miR-1-3p in Tetranychus cinnabarinus

As a newly recognized type of noncoding RNA, circular RNA can mediate a variety of physiological changes in mammals by regulating the post-transcriptional expression level of genes. However, the function of circRNA in the evolution of pesticide resistance in arthropods is still unknown. In this study, 2546 circRNAs were identified in Tetranychus cinnabarinus by transcriptome sequencing. The differentially expressed gene analysis indicated that 44 circRNAs were overexpressed in a cyflumetofen-resistant strain, of which a circRNA (named circ1-3p) was found to contain the response elements of miR-1-3p, an miRNA that is involved in cyflumetofen resistance by targeting TcGSTm04. The circular structure of circ1-3p was further determined using a divergent primer. The results of different molecular assays in vitro and in vivo showed that circ1-3p can compete with TcGSTm04 in miR-1-3p binding. The colocalization of circ1-3p and miR-1-3p was found using fluorescence in situ hybridization, suggesting that circ1-3p can directly sponge miR-1-3p in T. cinnabarinus. In addition, silencing the expression of circ1-3p resulted in the upregulation of miR-1-3p and the downregulation of TcGSTm04 as well as a significant increase in the sensitivity of T. cinnabarinus to cyflumetofen. All these pieces of evidence indicates that overexpressed circ1-3p promotes the expression of TcGSTm04 through sponging miR-1-3p, thereby involving in the evolution of cyflumetofen resistance in T. cinnabarinus.

Integrative analysis of circRNA/miRNA/mRNA regulatory network reveals the potential immune function of circRNAs in the Bombyx mori fat body

Bombyx mori nucleopolyhedrosis virus (BmNPV) is one of the greatest threats to sustainable development of the sericulture industry. Circular RNA (circRNA), a type of non-coding RNA, has been shown to play important roles in gene expression regulation, immune response, and diseases. The fat body is a tissue with both metabolic and immune functions. To explore the potential immune function of circRNAs, we analyzed differentially expressed (DE)circRNAs, microRNAs(miRNAs), and mRNAs in the B. mori fat body in response to BmNPV infection using high-throughput RNA sequencing. A total of 77 DEcircRNAs, 32 DEmiRNAs, and 730 DEmRNAs that are associated with BmNPV infection were identified. We constructed a DEcircRNA/DEmiRNA/DEmRNA and DEcircRNA/DEmiRNA/BmNPV gene regulatory network and validated the differential expression of circ_0001432 and its corresponding miRNA (miR-2774c and miR-3406-5p) and mRNA (778467 and 101745232) in the network. Tissue-specific expression of circ_0001432 and its expression at different time points were also examined. KEGG pathway analysis of DEmRNAs, target genes of DEmiRNAs, and host genes of DEcircRNAs in the network showed that these genes were enriched in several metabolic pathways and signaling pathways, which could play important roles in insect immune responses. Our results suggest that circRNA could be involved in immune responses of the B. mori fat body and help in understanding the molecular mechanisms underlying silkworm-pathogen interactions.

Systematic Identification and Functional Analysis of Circular RNAs During Rice Black-Streaked Dwarf Virus Infection in the Laodelphax striatellus (Fallén) Midgut

Circular RNAs (circRNAs) are endogenous RNAs that have critical regulatory roles in numerous biological processes. However, it remains largely unknown whether circRNAs are induced in response to plant virus infection in the insect vector of the virus as well as whether the circRNAs regulate virus infection. Rice black-streaked dwarf virus (RBSDV) is transmitted by Laodelphax striatellus (Fallén) in a persistent propagative manner and causes severe losses in East Asian countries. To explore the expression and function of circRNAs in the regulation of virus infection, we determined the circRNA expression profile in RBSDV-free or RBSDV-infected L. striatellus midgut tissues by RNA-Seq. A total of 2,523 circRNAs were identified, of which thirteen circRNAs were differentially expressed after RBSDV infection. The functions of these differentially circRNAs were predicted by GO and KEGG pathway analyses. The expression changes of five differentially expressed circRNAs and eight parental genes were validated by RT-qPCR. The circRNAs-microRNAs (miRNAs) interaction networks were analyzed and two miRNAs, which were predicted to bind circRNAs, were differentially expressed after virus infection. CircRNA2030 was up-regulated after RBSDV infection in L. striatellus midgut. Knockdown of circRNA2030 by RNA interference inhibited the expression of its predicted parental gene phospholipid-transporting ATPase (PTA) and enhanced RBSDV infection in L. striatellus. However, none of the six miRNAs predicting to bind circRNA2030 was up-regulated after circRNA2030 knockdown. The results suggested that circRNA2030 might affect RBSDV infection via regulating PTA. Our results reveal the expression profile of circRNAs in L. striatellus midgut and provide new insight into the roles of circRNAs in virus-insect vector interaction.

circEgg regulates histone H3K9me3 by sponging bmo-miR-3391-5p and encoding circEgg-P122 protein in the silkworm, Bombyx mori

A large number of circular RNAs (circRNAs) have been found in different organisms; however, their function in the regulation of histone modification remains unknown. In this study, we found that the circRNA circEgg, cyclized by the 9th-13th exon of Bombyx mori histone-lysine N-methyltransferase eggless (BmEgg) gene, mainly distributes in the cytoplasm, its expression levels changed with silkworm developmental stages, and the linear transcript level of the BmEgg gene was decreased when circEgg was overexpressed. Moreover, circEgg was found to repress histone H3 lysine 9 methylation (H3K9me3), promote histone H3 lysine 9 acetylation (H3K9ac), and positively regulate histone deacetylase (HDAC) Rpd3 (BmHDAC Rpd3) gene expression by sponging the microRNA bmo-miR-3391-5p. Furthermore, circEgg encodes a circEgg-P122 protein which appears to inhibit H3K9me3. These results suggest that circEgg regulates histone modification by sponging bmo-miR-3391-5p and encoding circEgg-P122 protein. To our knowledge, this is the first report showing that a circRNA produced by BmEgg plays an important role in histone epigenetic modification.

Post-transcriptional regulation of insect metamorphosis and oogenesis

Metamorphic transformation from larvae to adults along with the high fecundity is key to insect success. Insect metamorphosis and reproduction are governed by two critical endocrines, juvenile hormone (JH), and 20-hydroxyecdysone (20E). Recent studies have established a crucial role of microRNA (miRNA) in insect metamorphosis and oogenesis. While miRNAs target genes involved in JH and 20E-signaling pathways, these two hormones reciprocally regulate miRNA expression, forming regulatory loops of miRNA with JH and 20E-signaling cascades. Insect metamorphosis and oogenesis rely on the coordination of hormones, cognate genes, and miRNAs for precise regulation. In addition, the alternative splicing of genes in JH and 20E-signaling pathways has distinct functions in insect metamorphosis and oogenesis. We, therefore, focus in this review on recent advances in post-transcriptional regulation, with the emphasis on the regulatory role of miRNA and alternative splicing, in insect metamorphosis and oogenesis. We will highlight important new findings of miRNA interactions with hormonal signaling and alternative splicing of JH receptor heterodimer gene Taiman.

Systematic identification of circular RNAs and corresponding regulatory networks unveil their potential roles in the midguts of eastern honeybee workers

Currently, knowledge of circular RNAs (circRNAs) in insects including honeybee is extremely limited. Here, differential expression profiles and regulatory networks of circRNAs in the midguts of Apis cerana cerana workers were comprehensively investigated using transcriptome sequencing and bioinformatics. In total, 9589 circRNAs (201-800 nt in length) were identified from 8-day-old and 11-day-old workers' midguts (Ac1 and Ac2); among them, 5916 (61.70%) A. cerana cerana circRNAs showed conservation with our previously indentified circRNAs in Apis mellifera ligucstica workers' midguts (Xiong et al., Acta Entomologica Sinica 61:1363-1375, 2018). Five circRNAs were confirmed by RT-PCR and Sanger sequencing. Interestingly, novel_circ_003723, novel_circ_002714, novel_circ_002451, and novel_circ_001980 were highly expressed in both Ac1 and Ac2. In addition, the source genes of circRNAs were involved in 34 GO terms including organelle and cellular process and 141 pathways such as endocytosis and Wnt signaling pathway. Moreover, 55 DEcircRNAs including 34 upregulated and 21 downregulated circRNAs were identified in Ac2 compared with Ac1. circRNA-miRNA regulatory networks indicated that 1060 circRNAs can target 74 miRNAs; additionally, the DEcircRNA-miRNA-mRNA networks suggested that 13 downregulated circRNAs can bind to eight miRNAs and 29 miRNA-targeted mRNAs, while 16 upregulated circRNAs can link to 9 miRNAs and 29 miRNA-targeted mRNAs. These results indicated that DEcircRNAs as ceRNAs may play a comprehensive role in the growth, development, and metabolism of the worker's midgut via regulating source genes and interacting with miRNAs. Notably, eight DEcircRNAs targeting miR-6001-y were likely to be key participants in the midgut development. Our findings not only offer a valuable resource for further studies on A. cerana cerana circRNA and novel insights into understanding the molecular mechanisms underlying the midgut development of eastern honeybee but also provide putative circRNA candidates for functional research in the near future and novel biomarkers for identification of eastern honeybee species including A. cerana cerana and honeybee diseases such as chalkbrood and microsporidiosis.

Genome-Wide Identification of Circular RNAs Revealed the Dominant Intergenic Region Circularization Model in Apostichopus japonicus

Circular RNAs (circRNAs) were recently recognized to act as competing endogenous RNAs and play roles in gene expression regulation. Previous studies in humans and silkworms have shown that circRNAs take part in immune regulation. Here, we conducted coelomocyte circRNA sequencing to explore its immune functions in healthy and skin ulceration syndrome (SUS)-diseased sea cucumbers. A total of 3,592 circRNAs were identified in libraries with diversified circularization patterns compared with animal models. The common intron-pairing-driven circularization models are not popular in sea cucumber genome, which was replaced with intergenic region circularization. The accuracy of these identified circRNAs was further validated by Sanger sequencing and RNase R-treated assays. Expression profile analysis indicated that 117 circRNAs were upregulated and 144 circRNAs were downregulated in SUS-diseased condition, of which 71.6% were intergenic-type circRNAs. The interaction network of differentially expressed circRNAs and microRNAs (miRNAs) was constructed and showed that miR-2008 and miR-31, detected with significantly differential expression in SUS-affected samples in a previous study, were predicted to be regulated by 10 and 11 differentially expressed circRNAs with more than 10 binding sites, respectively. Moreover, seven circRNAs were further validated by quantitative real-time PCR, whose variation trends were consistent with circRNA sequencing. All our results supported that intergenic-type circRNAs might have a dominant function in Apostichopus japonicas immune response by acting as miRNA regulators.

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

Background

The 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.

Results

We 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.

Conclusions

We 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 material

The online version of this article (10.1186/s12864-018-5402-6) contains supplementary material, which is available to authorized users.

Systematic investigation of circular RNAs in Ascosphaera apis, a fungal pathogen of honeybee larvae

Ascosphaera apis is a widespread fungal pathogen of honeybee larvae, which causes heavy losses in apiculture. To date, knowledge about non-coding RNA (ncRNA) including circular RNA (circRNA) in A. apis is lacking. In this study, A. apis mycelia and spores were sequenced using RNA-seq technology. A total of 551 circRNAs were predicted on the basis of bioinformatic analyses, and most of the circRNAs were 200-600 bp in length, which were different from animal and plant circRNAs. In addition, the expression of six circRNAs in A. apis were confirmed using divergent and convergent PCR. Moreover, circRNA-microRNA regulation networks in A. apis were constructed, and further investigation showed that A. apis circRNAs could regulate gene expression by competitively binding miRNAs. GO and KEGG pathway enrichment analyses of the miRNAs target genes of circRNAs demonstrated that these A. apis circRNAs are likely to play key roles in metabolism, environmental response and gene expression.

Circular transcriptome sequencing of the middle silk gland and posterior silk gland in the Bombyx mori

Circular transcriptome sequencing of the middle silk gland (MSG) and posterior silk gland (PSG) in the Bombyx mori are presented. The middle silk gland and posterior silk gland were collected from the third day of fifth-instar B.mori larvae. The circular RNAs enriched by using RNase R to degrade the linear RNA molecules, and circular RNA sequencing (circRNA-seq) was performed using an Illumina Hiseq. 2500 sequencing platform. Samples are described in the SRA portal (SRP100385) and FASTQ files have been deposited in Sequence Read Archive (accession numbers: SRX2577343 and SRX2577342). The interpretation of these data is presented in the following research article: "Identification of circular RNA in the Bombyx mori silk gland" [1] (Gan et al., 2017).