Differential circular RNAs expression in ovary during oviposition in honey bees

Unraveling the modulatory manner and function of circRNAs in the Asian honey bee larval guts

Circular RNAs (circRNAs) are a class of non-coding RNAs (ncRNAs) that can participate in biological processes such as gene expression, growth, and development. However, little has been explored about the function of circRNAs in the development of Apis cerana larval guts. By using our previously gained deep sequencing data from the guts of A. cerana worker larvae at 4-, 5-, and 6-day-old (Ac4, Ac5, and Ac6 groups), the expression pattern and regulatory role of circular RNAs (circRNAs) during the development process was comprehensively investigated, with a focus on differentially expressed circRNAs (DEcircRNAs) relevant to immunity pathways and developmental signaling pathways, followed by validation of the binding relationships among a key competing endogenous RNA (ceRNA) axis. Here, 224 (158) DEcircRNAs were detected in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group. It's suggested that 172 (123) parental genes of DEcircRNAs were involved in 26 (20) GO terms such as developmental process and metabolic process and 138 (136) KEGG pathways like Hippo and Wnt signaling pathways. Additionally, ceRNA network analysis indicated that 21 (11) DEcircRNAs could target seven (three) DEmiRNAs, further targeting 324 (198) DEmRNAs. These DEmRNAs can be annotated to 33 (26) GO terms and 168 (200) KEGG pathways, including 12 (16) cellular and humoral immune pathways (endocytosis, lysosome, Jak-STAT, etc.) and 10 (nine) developmental signaling pathways (Hippo, mTOR, Hedgehog, etc.). Interestingly, DEcircRNAs in these two comparison groups could target the same ace-miR-6001-y, forming complex sub-networks. The results of PCR and Sanger sequencing confirmed the back-splicing sites within four randomly selected DEcircRNAs. RT-qPCR detection of these four DEcircRNAs verified the reliability of the used transcriptome data. The results of dual-luciferase reporter assay verified the binding relationships between novel_circ_001627 and ace-miR-6001-y and between ace-miR-6001-y and apterous-like. Our data demonstrated that DEcircRNAs were likely to modulate the developmental process of the A. cerana worker larval guts via regulation of parental gene transcription and ceRNA network, and novel_circ_001627/ace-miR-6001-y/apterous-like was a potential regulatory axis in the larval gut development. Findings from this work offer a basis and a candidate ceRNA axis for illustrating the circRNA-modulated mechanisms underlying the A. cerana larval guts.

Involvement of circular RNAs in the control of porcine ovarian cell functions: Upregulation by ciR-00596 and downregulation by ciR-00646

The current understanding of the role of circular RNAs (circRNAs) in regulating ovarian functions is inadequate. To assess the impact of ciR-00596 and ciR-00646 on the regulation of basic porcine ovarian granulosa cell functions, we conducted upregulation (utilizing overexpressing vectors) and downregulation (utilizing shRNA vectors) of these circRNAs. The relative expression of both circRNAs, cell viability and proliferation (accumulation of PCNA, cyclin B1, and XTT-positive cells), cytoplasmic (accumulation of bax and caspase-3) and nuclear (DNA fragmentation) apoptosis, and the release of progesterone, testosterone, estradiol, IGF-I, and oxytocin were evaluated. Transfection of cells with the ciR-00596 overexpression vector resulted in increases in cell viability and proliferation and the release of progesterone and IGF-I, while it decreased the cytoplasmic and nuclear apoptosis, testosterone, estradiol, and oxytocin output. CiR-00596 inhibition had the opposite effects. The overexpression of ciR-00646 decreased cell viability and proliferation, and the release of progesterone, IGF-I, and oxytocin, while increasing cytoplasmic and nuclear apoptosis and the output of testosterone and estradiol. Our findings are the first to show the stimulatory action of ciR-00596 and the inhibitory effect of ciR-00646 on ovarian cell functions, including the cell cycle, apoptosis, and secretory activity.

High expression circRALGPS2 in atretic follicle induces chicken granulosa cell apoptosis and autophagy via encoding a new protein

BACKGROUND

The reproductive performance of chickens mainly depends on the development of follicles. Abnormal follicle development can lead to decreased reproductive performance and even ovarian disease among chickens. Chicken is the only non-human animal with a high incidence of spontaneous ovarian cancer. In recent years, the involvement of circRNAs in follicle development and atresia regulation has been confirmed.

RESULTS

In the present study, we used healthy and atretic chicken follicles for circRNA RNC-seq. The results showed differential expression of circRALGPS2. It was then confirmed that circRALGPS2 can translate into a protein, named circRALGPS2-212aa, which has IRES activity. Next, we found that circRALGPS2-212aa promotes apoptosis and autophagy in chicken granulosa cells by forming a complex with PARP1 and HMGB1.

CONCLUSIONS

Our results revealed that circRALGPS2 can regulate chicken granulosa cell apoptosis and autophagy through the circRALGPS2-212aa/PARP1/HMGB1 axis.

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.

CircKDM5B sponges miR-128 to regulate porcine blastocyst development by modulating trophectoderm barrier function

Circular RNAs (circRNAs), which exert critical functions in the regulation of transcriptional and post-transcriptional gene expression, are found in mammalian cells but their functions in mammalian preimplantation embryo development remain poorly understood. Here, we showed that circKDM5B mediated miRNA-128 (miR-128) to regulate porcine early embryo development. We screened circRNAs potentially expressed in porcine embryos through an integrated analysis of sequencing data from mouse and human embryos, as well as porcine oocytes. An authentic circRNA originating from histone demethylase KDM5B (referred to as circKDM5B) was abundantly expressed in porcine embryos. Functional studies revealed that circKDM5B knockdown not only significantly reduced blastocyst formation but also decreased the number of total cells and trophectoderm (TE) cells. Moreover, the knockdown of circKDM5B resulted in the disturbance of tight junction assembly and impaired paracellular sealing within the TE epithelium. Mechanistically, miR-128 inhibitor injection could rescue the early development of circKDM5B knockdown embryos. Taken together, the findings revealed that circKDM5B functions as a miR-128 sponge, thereby facilitating early embryonic development in pigs through the modulation of gene expression linked to tight junction assembly.

The Potential Biological Roles of Circular RNAs in the Immune Systems of Insects to Pathogen Invasion

Circular RNAs (circRNAs) are a newly discovered class of endogenously expressed non-coding RNAs (ncRNAs). They are highly stable, covalently closed molecules that frequently exhibit tissue-specific expression in eukaryotes. A small number of circRNAs are abundant and have been remarkably conserved throughout evolution. Numerous circRNAs are known to play important biological roles by acting as microRNAs (miRNAs) or protein inhibitors ('sponges'), by regulating the function of proteins, or by being translated themselves. CircRNAs have distinct cellular functions due to structural and production differences from mRNAs. Recent advances highlight the importance of characterizing circRNAs and their targets in a variety of insect species in order to fully understand how they contribute to the immune responses of these insects. Here, we focus on the recent advances in our understanding of the biogenesis of circRNAs, regulation of their abundance, and biological roles, such as serving as templates for translation and in the regulation of signaling pathways. We also discuss the emerging roles of circRNAs in regulating immune responses to various microbial pathogens. Furthermore, we describe the functions of circRNAs encoded by microbial pathogens that play in their hosts.

Circular RNA ame_circ_2015 Function as microRNA Sponges in Regulating Egg-Laying of Honeybees (Apis mellifera)

Honeybees (Apis mellifera) are critical to maintaining ecological balance and are important pollinators. The oviposition behavior in honeybees is important and complex. Circular RNAs (circRNAs) are found to form circRNA-miRNA crosstalk and play important roles in reproduction processes. Here, dual luciferase reporter was used to confirm the crosstalk between ame_circ_2015 and ame_miR-14-3p. Functional experiments in vitro and in vivo were performed to investigate the biological functions of ame_circ_2015 in egg-laying of queens. The results showed that ame_circ_2015 directly target ame_miR-14-3p, and the expression of ame_circ_2015 was negatively correlated with ame_miR-14-3p expression. Overexpression results showed that ame_circ_2015 promoted the number of eggs laid and knockdown of ame_circ_2015 suppressed the number of eggs laid. It demonstrates that up-regulated ame_circ_2015 promotes the number of eggs laid by sponging ame_miR-14-3p. The study will provide information towards a better understanding of circRNA-miRNA crosstalk in egg-laying in honeybees.

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.

Identification of circular RNAs and corresponding regulatory networks reveals potential roles in the brains of honey bee workers exposed to dinotefuran

Honey bees are important and highly efficient pollinators of agricultural crops and have been negatively affected by insecticides in recent years. Circular RNA (circRNA) plays an important role in the regulation of multiple biological and pathological processes; however, its role in the honey bee brain after exposure to dinotefuran is not well understood. Here, the expression profiles and potential modulation networks of circRNAs in the brains of workers (Apis mellifera) were comprehensively investigated using RNA sequencing and bioinformatics. In total, 33, 144, and 211 differentially expressed (DE) circRNAs were identified on the 1st, 5th and 10th days after exposure to dinotefuran, respectively. Enrichment analyses revealed that the host genes of DE circRNAs were enriched in the Hippo signaling pathway-fly, Wnt signaling pathway, and neuroactive ligand-receptor interaction. circ_0002266, circ_0005080, circ_0010239 and circ_0005415 were found to have translational potential due to the presence of an internal ribosome entry site (IRES). An integrated analysis of the DE circRNA-miRNA-mRNA networks suggest that circ_0008898 and circ_0001829 may participate in the immune response to dinotefuran exposure by acting as miRNA sponges. Our results provide invaluable basic data on A. mellifera brain circRNA patterns and a molecular basis for further study of the biological function of circRNAs in the development and immune response of honey bees.

Identification and characterization of circRNAs in the liver of blunt snout bream (Megalobrama amblycephala) infected with Aeromonas hydrophila

Circular RNAs (circRNAs), a class of non-coding RNAs, play an important role in regulating various biological processes. In the present study, circRNAs from the Megalobrama amblycephala liver were identified at five different time points post Aeromonas hydrophila using RNA-seq technology. A total of 250 circRNAs were identified, of which 106 were differentially expressed (DE) in ten pairwise comparisons. GO and KEGG analyses showed that the parental genes of DE circRNAs were enriched in phagocytosis, complement and coagulation cascades, and Fc gamma R-mediated phagocytosis pathways. According to ceRNA hypothesis, the interaction network of circRNAs, miRNAs and mRNAs was constructed. Moreover, WGCNA was conducted, and five specific modules significantly related to bacterial infection were identified. All the above results reveal the important role of circRNAs in immune response, which enriches the information of circRNAs in teleost, and helps to understand the immune response mechanism of M. amblycephala to A. hydrophila.

Comprehensive CircRNA Profiling and Selection of Key CircRNAs Reveal the Potential Regulatory Roles of CircRNAs throughout Ovarian Development and Maturation in Cynoglossus semilaevis

Simple Summary

CircRNAs: as molecules involved in gene regulation, have become a new research hotspot in the non-coding RNA field. CircRNAs show tissue- or developmental stage-specific patterns of expression and can influence the expression levels of their parental genes. Recent studies have documented the potential biological roles of circRNAs in the growth, development, reproduction and health of humans and animals. Tongue sole (Cynoglossus semilaevis) is a marine flatfish that is an economically important farmed species in China. The commercial aquaculture of tongue sole has developed in the last few years because wild resources have decreased. Reproduction is regulated by brain-pituitary-gonad-liver axis which limits the development of artificial tongue sole culture. However, the roles of circRNAs in the ovarian development and maturation of tongue sole has never been reported. The identification of the potential functions of circRNAs provides a foundation for understanding the genetic mechanisms that regulate oocyte growth and maturation, which will allow the efficiency of tongue sole reproduction to be improved. Moreover, our findings extend the knowledge about a new type of endogenous RNA involved in regulating the ovarian development and maturation of tongue sole.

Abstract

CircRNAs are novel endogenous non-coding small RNAs involved in the regulation of multiple biological processes. However, little is known regarding circRNAs in ovarian development and maturation of fish. Our study, for the first time, provides the genome-wide overview of the types and relative abundances of circRNAs in tongue sole tissues during three ovarian developmental stages. We detected 6790 circRNAs in the brain, 5712 in the pituitary gland, 4937 in the ovary and 4160 in the liver. Some circRNAs exhibit tissue-specific expression, and qRT-PCR largely confirmed 6 differentially expressed (DE) circRNAs. Gene Ontology and KEGG pathway analyses of DE mRNAs were performed. Some DE circRNA parental genes were closely associated with biological processes in key signalling pathways and may play essential roles in ovarian development and maturation. We found that the selected circRNAs were involved in 10 pathways. RNase R digestion experiment and Sanger sequencing verified that the circRNA had a ring structure and was RNase R resistant. qRT-PCR results largely confirmed differential circRNA expression patterns from the RNA-seq data. These findings indicate that circRNAs are widespread in terms of present in production-related tissues of tongue sole with potentially important regulatory roles in ovarian development and maturation.

Differential Expression and Functional Analysis of CircRNA in the Ovaries of Low and High Fecundity Hanper Sheep

Simple Summary

Litter size is an important trait affecting reproductive capacity and breeding economics in meat sheep. Consequently, revealing its molecular mechanism helps us understand multiple lambs from the genetic perspective. In this study, we provide a genome-wide expression profile of circular RNAs (circRNAs) expression in Hanper sheep, which is a new breed of meat sheep raised by cross and self-group breeding for 15 years. In this study, ovarian circular RNAs and miRNAs associated with high and low fertility Hanper sheep are identified during the follicular and luteal phases of the estrous cycle, and their potential biological functions are predicted through Gene Ontology (GO), KEGG, GSEA, STEM, WGCNA analysis.

Abstract

Litter size is a considerable quality that determines the production efficiency of mutton sheep. Therefore, revealing the molecular regulation of high and low fertility may aid the breeding process to develop new varieties of mutton sheep. CircRNAs are the important factors regulating follicular development, but their mechanism role in the regulation of litter size in Hanper sheep is not clear. In the present study, ovarian tissues from the follicular (F) or luteal phase (L) of Hanper sheep that were either consecutive monotocous (M) or polytocous were collected. Then, we performed transcriptome sequencing to screen for differentially expressed circRNAs (DE-circRNAs) and elucidate their function. In total, 4256 circRNA derived from 2184 host genes were identified in which 183 (146 were upregulated, while 37 were downregulated) were differentially expressed in monotocous sheep in the follicular phase versus polytocous sheep in the follicular phase (MF vs. PF). Moreover, 34 circRNAs (14 were upregulated, while 20 were downregulated) were differentially expressed in monotocous sheep in the luteal phase versus polytocous sheep in the luteal sheep (ML vs. PL). This was achieved through DE-circRNAs function enrichment annotation analysis by GESA, GO, and KEGG, which function through the EGF-EGFR-RAS-JNK, TGF-β and thyroid hormone signaling pathway to affect the litter size of Hanper sheep in MF vs. PF and ML vs. PL. STEM results showed that MAPK signaling pathways play a key role in MF vs. PF and ML vs. PL. Through WGCNA analysis, AKT3 was a core gene in MF vs. PF and ML vs. PL. Moreover, competitive endogenous RNA (ceRNA) network analysis revealed the target binding sites for miRNA such as oar-miR-27a, oar-miR-16b, oar-miR-200a/b/c, oar-miR-181a, oar-miR-10a/b, and oar-miR-432 in the identified DE-cirRNAs.

Circular DDX10 is associated with ovarian function and assisted reproductive technology outcomes through modulating the proliferation and steroidogenesis of granulosa cells

circRNAs are present in human ovarian tissue, but how they regulate ovarian function remains unknown. In the current study, we investigated the levels of circRNAs in granulosa cells (GCs) derived from human follicular fluid, explored their correlation with female ovarian reserve function and clinical outcomes of assisted reproduction technique (ART), and investigated their effects on the biological functions of GC cell lines (COV434) in vitro. We identified that the levels of circDDX10 in GCs decreased gradually with aging (P < 0.01) and was positively correlated with AMH (r = 0.45, P < 0.01) and AFC (r = 0.32, P < 0.01), but not with FSH and estradiol (P > 0.05). Additionally, circDDX10 was related to the number of oocytes obtained, and good quality embryo rates. Silencing circDDX10 in GCs could markedly up-regulate the expression of apoptosis-related factors, reduce cell proliferation activity, inhibit the expression of steroid hormone synthesis-related factors, and prohibit the synthesis of estradiol. On the contrary, over-expression of circDDX10 had the opposite effect. circDDX10 is expected to become a novel biomarker for predicting the outcomes of ART, and may participate in the regulation of ovarian function by affecting the proliferation and apoptosis of GCs and steroid hormone synthesis.

CircRNA profiling reveals circ880 functions as miR-375-3p sponge in medaka gonads

Circular RNAs (circRNAs) have been regarded as regulators in the biological processes of various species. However, there is no report about circRNAs in the gonads of model fish medaka (Oryzias latipes). In this study, 1157 and 1570 circRNAs were obtained in the ovary and testis by RNA-sequencing. The characteristics of circRNAs were explored in sequence length, exon composition, and chromosome position. 24 circRNAs were significantly up or down-regulated in the testis compared to the ovary, 9 of which were verified by qRT-PCR. Interestingly, circ452 was highly expressed in the testis while circ880 expression exhibited sexual dimorphism. In situ hybridization (ISH) revealed that circ452 and circ880 were expressed in meiotic germ cells, and circ880 was also abundant in spermatogonia. In addition, dual-luciferase reporter assay manifested that circ880 and Oldnd can combine with miR-375-3p. Overall, these results provide emerging circRNA libraries and open new avenues for future investigation of circRNAs in medaka.

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.

Investigation of circular RNAs in an ectoparasitic mite Varroa destructor (Acarina: Varroidae) of the honey bee

Circular RNAs (circRNAs) are a large class of non-protein-coding transcripts that are involved in a diverse spectrum of regulatory mechanisms across a broad range of biological processes. To date, however, few studies on circRNAs have investigated their role in the biology of invertebrate parasites. The ectoparasitic mite Varroa destructor is perceived as the principal biotic threat towards global honey bee health. This parasite cannot be sustainably controlled partially due to the lack of knowledge about its basic molecular biology. In this paper, we unveil the circRNA profile of V. destructor for the first time and report the sources, distribution, and features of the identified circRNAs. Exonic, intronic, exon-intron, and intergenic circRNAs were discovered and exon-intron circRNAs were the most abundant within the largest spliced length. Three hundred and eighty-six (8.3%) circRNAs were predicted to possess translational potential. Eleven circRNAs, derived from six parental genes, exhibited strong bonds with miRNAs as sponges, suggesting an efficient post-transcriptional regulation. GO term and KEGG pathway enrichment analyses of the parental genes of the identified circRNAs showed that these non-coding RNAs were mainly engaged in protein processing, signal transduction, and various metabolism processes. To our knowledge, this is the first catalog of a circRNA profile of parasitiformes species, which reveals the prevalence of circRNAs in the parasite and provides biological insights for future genetic studies on this ubiquitous parasitic mite.

Research Progress on Oviposition-Related Genes in Insects

Oviposition-related genes have remained a consistent focus of insect molecular biology. Previous research has gradually clarified our mechanistic understanding of oviposition-related genes, including those related to oviposition-gland-related genes, oogenesis-related genes, oviposition-site-selection-related genes, and genes related to ovulation and hatching. Moreover, some of this research has revealed how the expression of single oviposition-related genes affects the expression of related genes, and more importantly, how individual node genes function to link the expression of upstream and downstream genes. However, the research to date is not sufficient to completely explain the overall interactions among the genes of the insect oviposition system. Through a literature review of a large number of studies, this review provides references for future research on oviposition-related genes in insects and the use of RNAi or CRISPR/Cas9 technology to verify the functions of oviposition-related genes and to prevent and control harmful insects.

Identification of circular RNAs in the ovarian follicles of Meishan and Duroc sows during the follicular phase

Circular RNAs (circRNAs) are a newly discovered class of endogenous non-coding RNAs that play an important role in growth and development by regulating gene expression and participating in a variety of biological processes. However, the role of circRNAs in porcine follicles remains unclear. Therefore, this study examined middle-sized ovarian follicles obtained from Meishan and Duroc sows at day 4 of the follicular phase. High-throughput RNA sequencing (RNA-seq) was utilized to construct circRNAs, and differential expression was identified. The findings were validated using reverse transcription PCR (RT-PCR) and DNA sequencing, GO and KEGG analyses were performed, and potential miRNA targets were identified. The RNA-seq identified a total of 15,866 circRNAs, with 244 differentially expressed in the Meishan relative to the Duroc (111 up-regulated and 133 down-regulated). The RT-PCR finding confirmed the RNA-seq results, and quantitative real-time PCR (qPCR) analysis examining a subset of the circRNAs showed that they are resistant to RNase R digestion. Bioinformatics analysis (GO and KEGG) showed that the host genes associated with the differentially expressed circRNAs are involved in reproduction and follicular development signaling pathways. Furthermore, many of the circRNAs were found to interact with miRNAs that are associated with follicular development. This study presents a new perspective for studying circRNAs and provides a valuable resource for further examination into the potential roles of circRNAs in porcine follicular development.

Ovarian Circular RNAs Associated with High and Low Fertility in Large White Sows during the Follicular and Luteal Phases of the Estrous Cycle

In this study, the ovarian tissues of Large White pigs were mined for novel circular RNAs (circRNAs), following which, their molecular characteristics and potential mechanisms for fertility regulation were examined. RNA sequencing was used for transcriptome analysis of ovarian follicles and corpora lutea in Large White sows with high (H) and low (L) fertility during the follicular (F) and luteal (L) phases of the estrous cycle. In total, 21,386 circRNA derived from 4535 host genes were identified. Differentially expressed circRNAs were detected in the LH vs. LL (1079) and in the FH vs. FL (1077) comparisons, and their host genes were enriched in steroid biosynthesis and forkhead box O (FOXO), thyroid hormone, cell cycle, and tumor growth factor (TGF)-beta signaling pathways. Protein–protein interaction networks were constructed on the basis of the host genes that were significantly enriched in pathways related to reproductive processes, with AKT3 and PP2CB serving as the hub genes in the networks of the LH vs. LL and FH vs. FL comparisons, respectively. The microRNA (miRNA) binding sites of the differentially expressed circRNAs were predicted, and 128 (LH vs. LL) and 113 (FH vs. FL) circRNA–miRNA pairs were identified. Finally, circRNA–miRNA negative regulatory networks were established on the basis of the gene expression profiles and bioinformatic analyses. In the current study, differentially expressed circRNAs were observed in ovarian tissues between the H and L fertility groups in both F and L phases of the estrous cycle, which suggested roles in pig fertility regulation. These findings provide new clues for elucidating fertility differences in pigs.

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.

Identification and characterization of human ovary-derived circular RNAs and their potential roles in ovarian aging

Circular RNAs (circRNAs) have recently been shown to exert effects on multiple pathological processes by acting as miRNA sponges. However, the roles of circRNAs in ovarian senescence are largely unknown. The objective of this study was to identify the circRNAs involved in ovarian aging and predict their potential biological functions. We first performed RNA-sequencing to generate ovarian circRNA expression profiles from young (n = 3) and aging (n = 3) groups. In total, 48,220 circRNAs were identified, of which 194 circRNAs were significantly up-regulated and 207 circRNAs were down-regulated during aging (fold change > 2, P < 0.05). Bioinformatics analysis demonstrated that the metabolic process, regulated secretory pathway, oxidation-reduction process, steroid hormone biosynthesis, and insulin secretion pathways, which may be associated with ovarian aging, were significantly enriched (P < 0.05). The biological characteristics of ovary-derived circRNA, such as back-splicing, RNase R resistance, stability, and alternative splicing, were further validated. Bioinformatics predicted that most of the circRNAs harboured miRNA binding sites, of which circDDX10-miR-1301-3p/miR-4660-SIRT3 axis may be involved in the regulation of ovarian function. Our study indicates that circRNAs are aberrantly expressed in the aging ovary and may play potential roles in the development of ovarian senescence.

Circular RNA profiling in the oocyte and cumulus cells reveals that circARMC4 is essential for porcine oocyte maturation

Thousands of circular RNAs (circRNAs) have been recently discovered in cumulus cells and oocytes from several species. However, the expression and function of circRNA during porcine oocyte meiotic maturation have been never examined. Here, we separately identified 7,067 and 637 circRNAs in both cumulus cells and oocytes via deep sequencing and bioinformatic analysis. Further analysis revealed that a faction of circRNAs is differentially expressed (DE) in a developmental stage-specific manner. The host genes of DE circRNAs are markedly enriched to multiple signaling pathways associated with cumulus cell function and oocyte maturation. Additionally, most DE circRNAs harbor several miRNA targets, suggesting that these DE circRNAs potentially act as miRNA sponge. Importantly, we found that maternal circARMC4 knockdown by siRNA microinjection caused a severely impaired chromosome alignment, and significantly inhibited first polar body extrusion and early embryo development. Taken together, these results demonstrate for the first time that circRNAs are abundantly and dynamically expressed in a developmental stage-specific manner in cumulus cells and oocytes, and maternally expressed circARMC4 is essential for porcine oocyte meiotic maturation and early embryo development.