Dynamic expression and functional analysis of circRNA in granulosa cells during follicular development in chicken

CircRNAs involved in the red light of effect on follicle selection in pigeons

Red light (RL) can enhance egg production in poultry. CircRNAs play a crucial role by serving as transcriptional regulators. However, their role in influencing follicle development in White King pigeons remains unexplored. In this study, 54 paired White King pigeons were chosen and divided into RL and white light (WL) groups, each with 3 subgroups. The egg production of paired pigeons in each replicate was recorded for 45 d, and the characteristics of follicle development were monitored during the laying interval (LI). The granulosa cell layer from follicles of the second-largest follicle (F2) was collected, and high-throughput sequencing was performed to elucidate the molecular mechanism of follicle development in pigeons. The study confirmed that RL enhances egg production in pigeons. Additionally, under RL, the F2 follicle was selected, while under WL, small follicles were kept on the third day (LI3). A total of 5,510 circRNAs were identified across all samples, revealing differentially expressed circRNAs (DECs) in various comparisons: 627 in RF1 vs. WF1, 900 in RF2 vs. WF2, 606 in RF1 vs. RF2, and 937 in WF1 vs. WF2. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that host genes of DECs were enriched in pathways like steroid hormone biosynthesis, oocyte meiosis, GnRH signaling pathway, and apoptosis pathway. Moreover, circRNA_5497, circRNA_2016, and circRNA_3328 were common DECs across 4 groups, sharing miRNA binding sites with follicle selection-associated genes. In conclusion, our findings suggest that RL promotes egg production by stimulating follicle selection during LI, offering insights into the regulatory mechanisms of circRNAs in follicle selection under RL. This knowledge can help enhance the reproductive performance of pigeons.

Diaph1 knockout inhibits mouse primordial germ cell proliferation and affects gonadal development

BACKGROUND

Exploring the molecular mechanisms of primordial germ cell (PGC) migration and the involvement of gonadal somatic cells in gonad development is valuable for comprehending the origins and potential treatments of reproductive-related diseases.

METHODS

Diaphanous related formin 1 (Diaph1, also known as mDia1) was screened by analyzing publicly available datasets (ATAC-seq, DNase-seq, and RNA-seq). Subsequently, the CRISPR-Cas9 technology was used to construct Diaph1 knockout mice to investigate the role of Diaph1 in gonad development.

RESULTS

Based on data from public databases, a differentially expressed gene Diaph1, was identified in the migration of mouse PGC. Additionally, the number of PGCs was significantly reduced in Diaph1 knockout mice compared to wild type mice, and the expression levels of genes related to proliferation (Dicer1, Mcm9), adhesion (E-cadherin, Cdh1), and migration (Cxcr4, Hmgcr, Dazl) were significantly decreased. Diaph1 knockout also inhibited Leydig cell proliferation and induced apoptosis in the testis, as well as granulosa cell apoptosis in the ovary. Moreover, the sperm count in the epididymal region and the count of ovarian follicles were significantly reduced in Diaph1 knockout mice, resulting in decreased fertility, concomitant with lowered levels of serum testosterone and estradiol. Further research found that in Diaph1 knockout mice, the key enzymes involved in testosterone synthesis (CYP11A1, 3β-HSD) were decreased in Leydig cells, and the estradiol-associated factor (FSH receptor, AMH) in granulosa cells were also downregulated.

CONCLUSIONS

Overall, our findings indicate that the knockout of Diaph1 can disrupt the expression of factors that regulate sex hormone production, leading to impaired secretion of sex hormones, ultimately resulting in damage to reproductive function. These results provide a new perspective on the molecular mechanisms underlying PGC migration and gonadal development, and offer valuable insights for further research on the causes, diagnosis, and treatment of related diseases.

Identification and Functional Analysis of circRNAs during Goat Follicular Development

Litter size is a crucial quantitative trait in animals, closely linked to follicular development. Circular RNA (circRNA), a type of single-stranded closed-loop endogenous RNA with stable expression, plays pivotal roles in various biological processes, yet its function in goat follicular development remains unclear. In this study, we collected large (follicle diameter > 3 mm) and small (1 mm < follicle diameter < 3 mm) follicles from black goats in the Chuanzhong region for circRNA sequencing, with the aim of elucidating the functional circRNAs that influence follicle development in goats. Differential analysis revealed that 17 circRNAs were upregulated in large follicles, and 28 circRNAs were upregulated in small follicles. Functional enrichment analysis revealed significant enrichment of pathways related to reproduction, including cellular response to follicle-stimulating hormone stimulus, the PI3K-Akt signaling pathway, the MAPK signaling pathway, and the Notch signaling pathway. Based on the ceRNA mechanism, 45 differentially expressed circRNAs were found to target and bind a total of 418 miRNAs, and an intercalation network including miR-324-3p (circRNA2497, circRNA5650), miR-202-5p (circRNA3333, circRNA5501), and miR-493-3p (circRNA4995, circRNA5508) was constructed. In addition, conservation analysis revealed that 2,239 circRNAs were conserved between goats and humans. Prediction of translation potential revealed that 154 circRNAs may potentially utilize both N6-methyladenosine (m6A) and internal ribosome entry site (IRES) translation mechanisms. Furthermore, the differential expression and circularization cleavage sites of five circRNAs were validated through RT-qPCR and DNA sequencing. Our study constructed a circRNA map in goat follicle development, offering a theoretical foundation for enhancing goat reproductive performance.

Multi-tissue transcriptome profiling linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken

Heterosis has been widely utilized in chickens. The nonadditive inheritance of genes contributes to this biological phenomenon. However, the role of circRNAs played in the heterosis is poorly determined. In this study, we observed divergent heterosis for residual feed intake (RFI) between 2 crossbreds derived from a reciprocal cross between White Leghorns and Beijing You chickens. Then, circRNA landscape for 120 samples covering the hypothalamus, liver, duodenum mucosa and ovary were profiled to elucidate the regulatory mechanisms of heterosis. We detected that a small proportion of circRNAs (7.83-20.35%) were additively and non-additively expressed, in which non-additivity was a major inheritance of circRNAs in the crossbreds. Tissue-specific expression of circRNAs was prevalent across 4 tissues. Weighted gene co-expression network analysis revealed circRNA-mRNA co-expression modules associated with feed intake and RFI in the hypothalamus and liver, and the co-expressed genes were enriched in oxidative phosphorylation pathway. We further identified 8 nonadditive circRNAs highly correlated with 16 nonadditive genes regulating negative heterosis for RFI in the 2 tissues. Circ-ITSN2 was validated in the liver tissue for its significantly positive correlation with PGPEP1L. Moreover, the bioinformatic analysis indicated that candidate circRNAs might be functioned by binding the microRNAs and interacting with the RNA binding proteins. The integration of multi-tissue transcriptome firstly linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken, which provide novel insights into the molecular mechanism underlying heterosis for feed efficiency. The validated circRNAs can act as potential biomarkers for predicting RFI and its heterosis.

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.

Circular RNA circRPS19 promotes chicken granulosa cell proliferation and steroid hormone synthesis by interrupting the miR-218-5p/INHBB axis

Ovarian follicle development is an important physiological activity for females and makes great significance in maintaining female health and reproduction performance. The development of ovarian follicle is mainly affected by the granulosa cells (GCs), whose growth is regulated by a variety of factors. Here, we identified a novel circular RNA (circRNA) derived from the Ribosomal protein S19 (RPS19) gene, named circRPS19, which is differentially expressed during chicken ovarian follicle development. Further explorations identified that circRPS19 promotes GCs proliferation and steroid hormone synthesis. Furthermore, circRPS19 was found to target and regulate miR-218-5p through a competitive manner with endogenous RNA (ceRNA). Functionals investigation revealed that miR-218-5p attenuates GCs proliferation and steroidogenesis, which is opposite to that of circRPS19. In addition, we also confirmed that circRPS19 upregulates the expression of Inhibin beta B subunit (INHBB) by binding with miR-218-5p to facilitate GCs proliferation and steroidogenesis. Overall, this study revealed that circRPS19 regulates GCs development by releasing the repression of miR-218-5p on INHBB, which suggests a novel mechanism in respect to circRNA and miRNA regulation in ovarian follicle development.

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.

Dynamic mRNA expression during chicken ovarian follicle development

Ovarian follicle development is a complex and well-orchestrated biological process of great economic significance for poultry production. Specifically, understanding the molecular mechanisms underlying follicular development is essential for high-efficiency follicular development can benefit the entire industry. In addition, domestic egg-laying hens often spontaneously develop ovarian cancer, providing an opportunity to study the genetic, biochemical, and environmental risk factors associated with the development of this cancer. Here, we provide high-quality RNA sequencing data for chicken follicular granulosa cells across 10 developmental stages, which resulted in a total of 204.57 Gb of clean sequencing data (6.82 Gb on average per sample). We also performed gene expression, time-series, and functional enrichment analyses across the 10 developmental stages. Our study revealed that SWF (small while follicle), F1 (F1 hierarchical follicles), and POFs (postovulatory follicles) best represent the transcriptional changes associated with the prehierarchical, preovulatory, and postovulatory stages, respectively. We found that the preovulatory stage F1 showed the greatest divergence in gene expression from the POF stage. Our research lays a foundation for further elucidation of egg-laying performance of chicken and human ovarian disease.

Transcriptome Analysis of Granulosa Cells Reveals Regulatory Mechanisms Related to Chicken Follicle Development

In this study, we aimed to better understand the difference between the functions of the two types of granulosa cells and sought to discover more key genes involved in follicle development and follicle selection. Herein, we separately collected pre-hierarchical follicle granulosa cells (PHGCs) and preovulatory follicle granulosa cells (POGCs) for RNA extraction; the transcriptomes of the two groups were compared via RNA-seq. A total of 5273 differentially expressed genes (DEGs) were identified between the PHGCs and POGCs; 2797 genes were up-regulated and 2476 were down-regulated in the PHGCs compared with the POGCs. A qPCR analysis confirmed that the expression patterns of 16 randomly selected DEGs were highly consistent with the RNA-seq results. In the POGCs, many of the genes with the most significant increase in expression were related to steroid hormone synthesis. In addition, the genes with the most significant decline in expression, including AMH and WT1, were related to the inhibition of steroid hormone synthesis. These results suggest that steroid hormones play a key role in follicle development. Furthermore, a Gene Ontology (GO) analysis revealed that these DEGs were mainly involved in the primary metabolic process, the carbohydrate metabolic process, the cellular process, ribosomes, the cytoplasm, and intracellular processes. A Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the DEGs were mainly enriched in steroid biosynthesis, the cell cycle, ribosomes, the TGF-beta signaling pathway, focal adhesion, and so on. We also observed the morphology of the follicles at different developmental stages, and the results showed that the thickness of the granular layer of the small yellow follicles (SYFs) decreased significantly with further development. In addition, we also found that the thickness of the granulosa layer of hens over 300 days old was significantly lower than that of 200-day-old hens. In short, these data indicate that the tissue morphology and function of granulosa cells change throughout follicle development.

microRNA transcriptome analysis of granulosa cells predicts that the Notch and insulin pathways affect follicular development in chickens

MicroRNAs (miRNAs) have been documented to play critical roles in chicken reproduction. Granulosa cell (GC) development of the follicle is closely related to hierarchical follicle ordering, making it an important factor in determining laying performance. Thus, it is meaningful to mine follicular development-related miRNAs. To identify regulatory miRNAs and the biological mechanisms by which they control follicular development, we conducted small RNA sequencing of GCs isolated from prehierarchical follicles named small yellow follicle (SYFG), the smallest hierarchical follicle (F6G), and the largest hierarchical follicle (F1G). A total of 99, 196, and 110 differentially expressed miRNAs (DEMs) were identified in SYFG.vs.F6G, SYFG.vs.F1G, and F6G.vs.F1G, respectively. Of these, 22 miRNAs, including miR-223, miR-103a, miR-449c-3p, and miR-203a, were ubiquitously identified as DEMs in three stages. Target gene prediction suggested that these miRNAs are associated with the MAPK, TGF-β, and Wnt signaling pathways, which are all associated with follicular development. The Notch and insulin signaling pathways were commonly enriched in all three comparisons. RT-qPCR analysis further indicated that the expression levels of PSEN2, which encodes an essential factor regulating Notch and insulin signaling, was significantly changed in SYFG, F6G, and F1G. The current study provides basic data and offers a new foundation for further exploration of the roles of miRNAs in follicular development in chickens.

Transcriptome identification and characterization of long non-coding RNAs in the ovary of hens at four stages

Long non-coding RNAs (lncRNAs) play important roles in transcriptional and post-transcriptional regulation. LncRNAs, which are defined as non-coding RNAs more than 200 bp in length, are involved in key biological processes, such as cell proliferation and differentiation, epigenetic regulation, and gene transcriptional translation. Recent studies have shown that lncRNAs also play major regulatory roles in the reproduction of mammals. However, knowledge of the roles of lncRNAs in the chicken ovary lacking. In this study, we performed RNA-seq analyses of ovarian tissue from Hy-Line brown laying hens at four physiological stages [15, 20, 30, and 68 weeks of age (W)]. We identified 657 lncRNA transcripts that were differentially expressed during ovarian development, the number of down-regulated lncRNAs was higher than the number of up-regulated lncRNAs during development. We predicted the cis and trans target genes of the DE lncRNAs and constructed a lncRNA-mRNA interaction network, which indicated that the DE genes (DEGs) and the target genes of the DE lncRNAs are mainly involved in signaling pathways associated with ovarian development, including oocyte meiosis, calcium signaling pathways, ECM-receptor interactions, and ribosome and focal adhesion. Overall, we found that twelve lncRNAs were strongly involved in ovarian development: LNC_013443, LNC_001029, LNC_005713, LNC_016762, ENSGALT00000101857, LNC_003913, LNC_013692, LNC_012219, LNC_004140, ENSGALT00000096941, LNC_009356, and ENSGALT00000098716. In summary, our study utilized RNA-seq analysis of hen ovaries to explore key lncRNAs involved in ovarian development and function. Furthermore, the comprehensive analysis identified the target genes of these lncRNAs providing a better understanding of the mechanisms underlying ovarian development in hens and a theoretical basis for further research.

Circular RNAs are associated with the resistance to Newcastle disease virus infection in duck cells

Introduction

Newcastle disease virus (NDV) is prevalent worldwide with an extensive host range. Among birds infected with velogenic NDV strains, chickens experience high pathogenicity and mortality, whereas ducks mostly experience mild symptoms or are asymptomatic. Ducks have a unique, innate immune system hypothesized to induce antiviral responses. Circular RNAs (circRNAs) are among the most abundant and conserved eukaryotic transcripts. These participate in innate immunity and host antiviral response progression.

Methods

In this study, circRNA expression profile differences post-NDV infection in duck embryo fibroblast (DEF) cells were analyzed using circRNA transcriptome sequencing. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to reveal significant enrichment of differentially expressed (DE) circRNAs. The circRNA-miRNA-mRNA interaction networks were used to predict the related functions of circRNAs. Moreover, circ-FBXW7 was selected to determine its effect on NDV infection in DEFs.

Results

NDV infection altered circRNA expression profiles in DEF cells, and 57 significantly differentially expressed circRNAs were identified post-NDV infection. DEF responded to NDV by forming circRNAs to regulate apoptosis-, cell growth-, and protein degradation-related pathways via GO and KEGG enrichment analyses. circRNA-miRNA-mRNA interaction networks demonstrated that DEF cells combat NDV infection by regulating cellular pathways or apoptosis through circRNA-targeted mRNAs and miRNAs. circ-FBXW7 overexpression and knockdown inhibited and promoted viral replication, respectively. DEF cells mainly regulated cell cycle alterations or altered cellular sensing to combat NDV infection.

Conclusion

These results demonstrate that DEF cells exert antiviral responses by forming circRNAs, providing novel insights into waterfowl antiviral responses.

Circular RNA as a Novel Regulator and Promising Biomarker in Polycystic Ovary Syndrome

Polycystic ovary syndrome (PCOS) is a prevalent metabolic and reproductive disorder that causes low fertility in females. Despite its detrimental effects on women's health, care for PCOS has been impeded by its undefined pathogenesis. Thus, there is an urgent need to explore novel biomarkers and therapeutic targets for the diagnosis and treatment of PCOS. Circular RNAs (circRNAs) are a class of noncoding RNAs with covalently closed cyclic structures, present in high abundance, and show development-stage specific expression patterns. Recent studies have demonstrated that circRNAs participate in PCOS progression by modulating various biological functions, including cell proliferation, apoptosis, and steroidogenesis. In addition, circRNAs are widely present in the follicular fluid of women with PCOS, indicating their potential as diagnostic biomarkers and therapeutic targets for PCOS. This review provides the current knowledge of circRNAs in PCOS, including their regulatory functions and molecular mechanisms, and explores their potential as diagnostic biomarkers and therapeutic targets.

Effects of quercetin on granulosa cells from prehierarchical follicles by modulating MAPK signaling pathway in chicken

Quercetin (Que), widely found in a huge variety of plants, plays important roles in ovarian function. However, to data, there have been no reports about Que regulating granulosa cells (GCs) in prehierarchical follicles in chicken. Herein, GCs from follicles diameter from 4 to 8 mm in chicken were treated by Que in vitro culture to investigate how Que exerts its effect on follicular development. GCs treated by Que in concentrations of 10, 100, and 1,000 ng/mL were tested for cell proliferation and progesterone secretion. Eight cDNA libraries were constructed from GCs (4 samples per group) to explore transcriptome expression changes. The role of the MAPK/ERK signaling pathway was validated in this process. Treatment with 100 and 1,000 ng/mL levels of Que significantly promoted cell proliferation and progesterone secretion (P < 0.05). RNA-seq analysis data showed that 402 and 263 differentially expressed genes (DEGs) were up- and down-regulated, respectively. Functional enrichment analysis that the pathways related to follicular development included biosynthesis of amino acids, MAPK signaling pathway, and calcium signaling pathway. Notably, the function exerted in GCs of the different levels of Que was associated with the suppression of the MAPK pathway. In conclusion, our results proved that low levels of Que could promote MAPK signaling pathway, but high levels of Que inhibit MAPK signaling pathway in GCs from the prehierarchical follicles, promote cell proliferation and progesterone secretion, and benefit follicle selection.

Transcriptome RNA Sequencing Reveals That Circular RNAs Are Abundantly Expressed in Embryonic Breast Muscle of Duck

Circular RNAs are widespread in various species and have important roles in myogenesis. However, the circular RNAs involved in breast muscle development in ducks have not yet been studied. Here, to identify circular RNAs during duck skeletal muscle development, three pectorales from Shan Ma ducks at E13 and E19, which represent undifferentiated and differentiated myoblasts, respectively, were collected and subjected to RNA sequencing. A total of 16,622 circular RNAs were identified, of which approximately 80% were exonic circular RNAs and 260 were markedly differentially expressed between E19 and E13. The parental genes of the differentially expressed circular RNAs were significantly enriched in muscle-related biological processes. Moreover, we found that the overexpression of circGAS2-2 promoted cell cycle progression and increased the proliferation viability of duck primary myoblasts; conversely, knockdown of circGAS2-2 retarded the cell cycle and reduced the proliferation viability of myoblasts. Taken together, our results demonstrate that circular RNAs are widespread and variously expressed during the development of duck skeletal muscle and that circGAS2-2 is involved in the regulation of myogenesis.

Interacting Networks of the Hypothalamic-Pituitary-Ovarian Axis Regulate Layer Hens Performance

Egg production is a vital biological and economic trait for poultry breeding. The 'hypothalamic-pituitary-ovarian (HPO) axis' determines the egg production, which affects the layer hens industry income. At the organism level, the HPO axis is influenced by the factors related to metabolic and nutritional status, environment, and genetics, whereas at the cellular and molecular levels, the HPO axis is influenced by the factors related to endocrine and metabolic regulation, cytokines, key genes, signaling pathways, post-transcriptional processing, and epigenetic modifications. MiRNAs and lncRNAs play a critical role in follicle selection and development, atresia, and ovulation in layer hens; in particular, miRNA is known to affect the development and atresia of follicles by regulating apoptosis and autophagy of granulosa cells. The current review elaborates on the regulation of the HPO axis and its role in the laying performance of hens at the organism, cellular, and molecular levels. In addition, this review provides an overview of the interactive network regulation mechanism of the HPO axis in layer hens, as well as comprehensive knowledge for successfully utilizing their genetic resources.

Regulatory effects of circular RNA on hypoxia adaptation in chicken embryos

The Tibetan chicken, a native breed of the Tibetan plateau, is adapted to the high-altitude and hypoxic environment of the plateau. As endogenous molecules, circular RNAs (circRNAs) have been shown to play an important role in the adaptation to hypoxic environments and regulation of angiogenesis. In this study, highland Tibetan and lowland Chahua chicken eggs were incubated in a hypoxic environment and the chorionic allantoic membrane was collected for Ribo-Zero RNA sequencing. A total of 1,414 circRNAs, mostly derived from exons, were identified. Of these, 93 differentially expressed circRNAs were detected between Tibetan and Chahua chickens. Combined with the differentially expressed miRNAs and mRNAs identified in our previous study, we identified four circRNAs (circBRD1, circPRDM2, circPTPRS, and circDENND4C). These circRNAs may act as competing endogenous RNA (ceRNA) to upregulate APOA1 expression by absorbing novel_miR_589, thereby regulating angiogenesis and affecting hypoxia adaptation in chicken embryos. The regulatory circRNAs/novel_miR_589/APOA1 axis provides valuable evidence for a better understanding of the specific functions and molecular mechanisms of circRNAs in plateau hypoxia adaptation in Tibetan chickens.

Unveiling the circRNA-Mediated Immune Responses of Western Honey Bee Larvae to Ascosphaera apis Invasion

Western honey bee (Apis mellifera), a eusocial insect with a superior economic and ecological value, is widely used in the beekeeping industry throughout the world. As a new class of non-coding RNAs (ncRNAs), circular RNAs (circRNAs) participate in the modulation of considerable biological processes, such as the immune response via diverse manners. Here, the identification, characteristic investigation, and molecular verification of circRNAs in the Apis mellifera ligustica larval guts were conducted, and the expression pattern of larval circRNAs during the Ascosphaera apis infection was analyzed, followed by the exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 2083 circRNAs in the larval guts of A. m. ligustcia were identified, with a length distribution ranging from 106 nt to 92,798 nt. Among these, exonic circRNAs were the most abundant type and LG1 was the most distributed chromosome. Additionally, 25, 14, and 30 up-regulated circRNAs as well as 26, 25, and 62 down-regulated ones were identified in the A. apis-inoculated 4-, 5-, and 6-day-old larval guts in comparison with the corresponding un-inoculated larval guts. These DEcircRNAs were predicted to target 35, 70, and 129 source genes, which were relative to 12, 23, and 20 GO terms as well as 11, 10, and 27 KEGG pathways, including 5 cellular and humoral immune pathways containing apoptosis, autophagy, endocytosis, MAPK, Toll, and Imd signaling pathways. Furthermore, complex competing endogenous RNA (ceRNA) regulatory networks were detected to be formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The Target DEmRNAs were engaged in 24, 20, and 25 functional terms as well as 62, 80, and 159 pathways, including several vital immune defense-associated pathways, namely the lysosome, endocytosis, phagosome, autophagy, apoptosis, MAPK, Jak-STAT, Toll, and Imd signaling pathways. Finally, back-splicing sites within 15 circRNAs and the difference in the 9 DEcircRNAs' expression between un-inoculated and A. apis-inoculated larval guts were confirmed utilizing 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. m. ligustica larvae against A. apis invasion.

Expression spectrum of new circFBN1 in various tissues and follicles of Taihang chicken

In this study, the new circFBN1 derived from the FBN1 gene was identified. The expression spectrum of circFBN1 in various tissues and follicles of Taihang chickens was verified and analyzed. The effects on proliferation of GCs were checked. The results revealed that circFBN1 is indeed present and was differentially expressed in tissues and follicles and significantly promoted cell proliferation. In conclusion, our results suggested that circFBN1 could affect chicken follicular development by regulating the proliferation of GCs. These results will enable us recognize the molecular mechanisms in animal reproductive regulation.

Differential expression of circRNAs of testes with high and low sperm motility in Yili geese

The aim of this study was to explore the potential biological function of circular RNAs (circRNAs) in the sperm motility traits of Xinjiang Yili geese, and to provide a reference for analyzing the mechanism of regulation of Yili geese sperm motility. The 10 selected Xinjiang Yili Geese with high or low sperm motility (five for each group) were 3 years old, in good health, and were kept in the same feeding conditions. Yili geese were slaughtered for the collection of testicular tissue and high-throughput sequencing technology was used to screen differentially expressed circRNAs for bioinformatics analysis. Combined with the previously screened miRNAs related to the sperm motility of Yili geese, the circRNAs miRNAs regulatory network was constructed. The results showed that a total of 26,311 circRNAs were obtained from testicular tissues with high and low sperm motility, and 173 DECs were screened between the two groups (p < 0.05, |log2Foldchange|>0), of which 82 were up-regulated and 91 were down-regulated. Functional analysis of the source genes of these DECs showed that the source genes were mainly involved in biological processes. KEGG enrichment analysis showed that the source genes of DECs were mainly enriched in autophagy-animal, ubiquinone and other terpenoid-quinone biosynthesis, progesterone-mediated oocyte maturation, regulation of the actin cytoskeleton and other pathways. Furthermore, the visual regulatory network of differential circRNA-miRNA-mRNA was constructed, including 20 circRNAs, 18 miRNAs and 177 mRNAs, and nine core regulatory circRNAs were screened, including novell_circ_0045314, novel_circ_0019994 and novel_circ_0020422, etc., targeting ppy-mir-16, hsa-mir-221–3p, gga-mir-499–5p, etc. The results suggest that circRNAs may interact with miRNAs to further regulate mRNA to regulate sperm motility in Yili geese, so as to provide a reference for analyzing the molecular mechanism of sperm motility regulation.

Expression profile of circular RNAs in continuous light-induced ovarian dysfunction

PURPOSE

This study aims to elucidate the underlying relationship between the expression profiles of circular RNAs (circRNAs) and the ovarian dysfunction induced by continuous light.

METHODS

High-throughput sequencing was used to profile the transcriptome of differentially expressed circRNAs (DEcircRNAs) in rat ovary under continuous light exposure (12 h:12 h light/light cycle, L/L group) and a control cycle (12 h:12 h light/dark cycle, L/D group). Gene ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and circRNAs-microRNAs-messenger RNAs networks were performed to predict the role of DEcircRNAs in biological processes and pathways. A quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) assay was conducted to verify the high-throughput sequencing results and the expression level of circadian rhythm genes.

RESULTS

In total, 305 circRNAs were differentially expressed between the L/L and L/D groups. Among these, 211 circRNAs were up regulated, while 94 were down regulated. Eight candidate circRNAs from 305 DEcircRNAs were verified by qRT-PCR. Further bioinformatics analysis revealed that interactions between DEcircRNAs and a set of microRNAs involved in ovarian dysfunction-related pathways, such as regulation of androgen receptors, gonadotrophin releasing hormone signaling pathway, endocrine resistance, etc. Subsequently, we identified rno_circ:chr2:86868285-86964272 and rno_circ:chr1:62330221-62360073 may participate in the pathophysiology of ovarian dysfunction by constructing circRNAs-microRNAs-messenger RNAs networks. Meanwhile, constant light reduced the expression of circadian rhythm genes CLOCK, BAML1, PER1, and PER2 compared with that of controls. Caspase3 and Bax were up regulated in the L/L group compared with the L/D group, while Bcl-2 was down regulated.

CONCLUSIONS

In summary, the results reveal that the expression profiles and potential functions of DEcircRNAs in rat ovaries may play important roles in continuous light-induced ovarian dysfunction. These findings provide novel clues and molecular targets for studying the mechanisms and clinical therapy of ovarian dysfunction.

Circular RNA expression profiles and CircSnd1-miR-135b/c-foxl2 axis analysis in gonadal differentiation of protogynous hermaphroditic ricefield eel Monopterus albus

BACKGROUND

The expression and biological functions of circular RNAs (circRNAs) in reproductive organs have been extensively reported. However, it is still unclear whether circRNAs are involved in sex change. To this end, RNA sequencing (RNA-seq) was performed in gonads at 5 sexual stages (ovary, early intersexual stage gonad, middle intersexual stage gonad, late intersexual stage gonad, and testis) of ricefield eel, and the expression profiles and potential functions of circRNAs were studied.

RESULTS

Seven hundred twenty-one circRNAs were identified, and the expression levels of 10 circRNAs were verified by quantitative real-time PCR (qRT-PCR) and found to be in accordance with the RNA-seq data, suggesting that the RNA-seq data were reliable. Then, the sequence length, category, sequence composition and the relationship between the parent genes of the circRNAs were explored. A total of 147 circRNAs were differentially expressed in the sex change process, and GO and KEGG analyses revealed that some differentially expressed (such as novel_circ_0000659, novel_circ_0004005 and novel_circ_0005865) circRNAs were closely involved in sex change. Furthermore, expression pattern analysis demonstrated that both circSnd1 and foxl2 were downregulated in the process of sex change, which was contrary to mal-miR-135b. Finally, dual-luciferase reporter assay and RNA immunoprecipitation showed that circSnd1 and foxl2 can combine with mal-miR-135b and mal-miR-135c. These data revealed that circSnd1 regulates foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c.

CONCLUSION

Our results are the first to demonstrate that circRNAs have potential effects on sex change in ricefield eel; and circSnd1 could regulate foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c. These data will be useful for enhancing our understanding of sequential hermaphroditism and sex change in ricefield eel or other teleosts.

Regulation of Non-Coding RNA in the Growth and Development of Skeletal Muscle in Domestic Chickens

Chicken is the most widely consumed meat product worldwide and is a high-quality source of protein for humans. The skeletal muscle, which accounts for the majority of chicken products and contains the most valuable components, is tightly correlated to meat product yield and quality. In domestic chickens, skeletal muscle growth is regulated by a complex network of molecules that includes some non-coding RNAs (ncRNAs). As a regulator of muscle growth and development, ncRNAs play a significant function in the development of skeletal muscle in domestic chickens. Recent advances in sequencing technology have contributed to the identification and characterization of more ncRNAs (mainly microRNAs (miRNAs), long non-coding RNAs (LncRNAs), and circular RNAs (CircRNAs)) involved in the development of domestic chicken skeletal muscle, where they are widely involved in proliferation, differentiation, fusion, and apoptosis of myoblasts and satellite cells, and the specification of muscle fiber type. In this review, we summarize the ncRNAs involved in the skeletal muscle growth and development of domestic chickens and discuss the potential limitations and challenges. It will provide a theoretical foundation for future comprehensive studies on ncRNA participation in the regulation of skeletal muscle growth and development in domestic chickens.

Differentially Expressed Circular RNA Profile Signatures Identified in Prolificacy Trait of Yunshang Black Goat Ovary at Estrus Cycle

CircRNAs acting as miRNA sponges play important roles in the growth process of animal individuals. The prolificacy trait of goats is involved in many pathways, however, the variation of circRNA expression profiles in the different phases of the estrus cycle at high and low fecundity groups is still unknown. Here, we analyzed the circRNA profiles of ovarian tissues among high and low fecundity groups in the follicular phase (HF vs LF), high and low fecundity groups in the luteal phase (HL vs LL), and high and low fecundity in the whole estrus cycle (HF vs HL and LF vs LL) using RNA-seq. A total of 283 (114 upregulated and 169 downregulated), 559 (299 upregulated and 260 downregulated), 449 (254 upregulated and 195 downregulated), and 314 (210 upregulated and 104 downregulated) differentially expressed (DE) circRNAs were screened in HF vs LF, HF vs HL, HL vs LL, and LF vs LL groups, respectively. Enrichment analysis suggested that the targeting of DE circRNAs was mainly enriched in oocyte meiosis, the GnRH signaling pathway, and estrogen signaling pathway. After integrating our previous study of miRNA-seq, there were 56 miRNAs that could target to 192 DE circRNAs, including the miR-133 family (including miR-133a-3p and miR-133b), miR-129-3p, and miR-21, which also had important influence on the prolificacy trait of goats. Then, 18 circRNAs with coding potential were obtained by four software predictions, and 9 circRNAs were validated by RT-qPCR. Together, circRNAs play a key role in the prolificacy trait and the transformation of the follicular phase to the luteal phase in the estrus cycle of goats.

Transcriptome analyses of potential regulators of pre- and post-ovulatory follicles in the pigeon (Columba livia)

To identify the dominant genes controlling follicular maturation, ovulation and regression for pigeon, we used RNA-seq to explore the gene expression profiles of pre- and post-ovulatory follicles of pigeon. We obtained total of 4.73million (96% of the raw data) high-quality clean reads, which could be aligned with 20282 genes. Gene expression profile analysis identified 1461 differentially expressed genes (DEGs) between the pre- (P4) and post-ovulatory follicles (P5). Of these, 843 genes were upregulated, and 618 genes were down-regulated. Furthermore, many DEGs were significantly enriched in some pathways closely related to follicle maturation, ovulation and regression, such as ECM-receptor interaction, vascular smooth muscle contraction, progesterone-mediated oocyte maturation, phagosome. Importantly, the DGEs in ECM-receptor interaction pathway included COL1A1 , COL1A2 , COL4A1 , COL4A2 , ITGA11 , ITGB3 and SDC3 , in the progesterone-mediated oocyte maturation pathway involved CDK1 , CDC25A , CCNB3 , CDC20 and Plk1 , and in the vascular smooth muscle contraction covered CALD1 , KCNMA1 , KCNMB1 , CACNA1 , ACTA2 , MYH10 , MYL3 , MYL6 , MYL9 , closely related to promoting follicular maturation and ovulation in pre-ovulatory follicles. Moreover, it seems that the lysosomal cathepsin family has a decisive role in the regression of early stage of post-ovulatory follicle. Taken together, these data enrich the research of molecular mechanisms of pigeon follicular activities at the transcriptional level and provide novel insight of breeding-related physiology for birds.

CircRNA expression in chicken granulosa cells illuminated with red light

Red light (RL) can improve egg production in Jinghai Yellow hens. Circular RNAs (circRNAs) are novel, non-coding RNAs, but the molecular mechanism underlying circRNA function during follicular development in hens under monochromatic light has not been established. Herein, we compared expression profiles of granulosa cells (GCs) from small yellow follicles (SYFs) from hens under RL and white light (WL). A total of 2,468 circRNAs were identified, of which 22 were differentially expressed (DE) in the RL and WL groups. DE circRNA host genes were enriched in ovarian steroidogenesis, and MAPK and PI3K-Akt signaling pathways. Furthermore, DE circRNA_0320 and circRNA_0185 interacted with miR-143-3p, which targets the follicle-stimulating hormone receptor and is essential for GC differentiation and follicle development. These findings will facilitate further analysis of the molecular mechanism leading to GC development in hens raised under monochromatic light, which could lead to increased egg production.

Expression Profile Analysis to Identify Circular RNA Expression Signatures in the Prolificacy Trait of Yunshang Black Goat Pituitary in the Estrus Cycle

The pituitary gland is an important organ. It is a complex area of the brain involved in endocrine function and reproductive regulation. However, the function of the pituitary in goat reproduction is still unclear. Herein, RNA sequencing was used to explore the expression patterns of circle RNAs (circRNAs) in the pituitary of Yunshang black goats during the various estrus phases. Then the host genes of the circRNAs were predicted, and a competing endogenous RNA (ceRNA) network was constructed. The results showed a total of 6,705 circRNAs in the pituitary of Yunshang black goats, among which 388 differentially expressed (DE) circRNAs (214 were upregulated, while 174 were downregulated) were identified between high- and low-yield Yunshang black goats in the follicular phase (HF vs. LF); moreover, 361 DE circRNAs (136 were upregulated, while 225 were downregulated) were identified between high- and low-yield Yunshang black goats in the luteal phase (HL vs. LL). There were 65 DE circRNAs targeting 40 miRNAs in the HF vs. LF comparison and 46 DE circRNAs targeting 31 miRNAs in the HL vs. LL comparison. We identified chi_circ_0030920, chi_circ_0043017, chi_circ_0008353, chi_circ_0041580, and chi_circ_0016478 as the key circRNAs through functional enrichment analysis. The ceRNA network analysis showed that chi_circ_0031209 and chi_circ_0019448 might play an important role in reproduction by influencing the expression of prolactin receptor (PRLR) in high- and low-yielding goats during the luteal phase, whereas chi_circ_0014542 regulates the expression of WNT5A during the follicular phase. Our study provided the overall expression profiles of circRNAs in the goat pituitary during the estrus phase, which provides new insight into the mechanism of high-yield goats, which can be helpful to guide goat breeding.

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

Background

Circular RNAs (circRNAs) are a group of noncoding RNAs that participate in gene expression regulation in various pathways. The essential roles of circRNAs have been revealed in many species. However, knowledge of circRNAs in fungi is still not comprehensive.

Results

Trichophyton rubrum (T. rubrum) is considered a model organism of human pathogenic filamentous fungi and dermatophytes. In this study, we performed a genome-wide investigation of circRNAs in T. rubrum based on high-throughput sequencing and ultimately identified 4254 circRNAs. Most of these circRNAs were specific to the conidial or mycelial stage, revealing a developmental stage-specific expression pattern. In addition, 940 circRNAs were significantly differentially expressed between the conidial and mycelial stages. PCR experiments conducted on seven randomly selected differentially expressed (DE-) circRNAs confirmed the circularized structures and relative expression levels of these circRNAs. Based on their genome locations, most circRNAs originated from intergenic regions, unlike those in plants and animals. Furthermore, we constructed circRNA-miRNA-mRNA regulatory networks that included 661 DE-circRNAs targeting 140 miRNAs and further regulating 2753 mRNAs. The relative expression levels of two randomly selected circRNA-miRNA-mRNA axes were investigated by qRT-PCR, and the competing endogenous RNA (ceRNA) network theory was validated. Functional enrichment analysis of the target genes suggested that they were significantly involved in posttranscriptional processes and protein synthesis as well as some small-molecule metabolism processes. CircRNAs are relatively more conserved in closely related dermatophytes but rarely conserved in distantly related species. Tru_circ07138_001 is a highly conserved circRNA that was conserved in all ten dermatophytes analyzed in our study and three distantly related species. Its host gene TERG_07138 was also highly conserved in two of these distantly related species Gallus gallus and Caenorhabditis elegans. The specific role of this circRNA deserves further exploration.

Conclusions

Our study is the first to provide a global profile of circRNAs in T. rubrum as well as dermatophytes. These results could serve as valuable resources for research on circRNA regulatory mechanisms in fungi and reveal new insights for further investigation of the physical characteristics of these significant human fungal pathogens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08184-y.

Genome-Wide Gene Expression Profiles Reveal Distinct Molecular Characteristics of the Goose Granulosa Cells

Granulosa cells (GCs) are decisive players in follicular development. In this study, the follicle tissues and GCs were isolated from the goose during the peak-laying period to perform hematoxylin-eosin staining and RNA-seq, respectively. Moreover, the dynamic mRNA and lncRNA expression profiles and mRNA-lncRNA network analysis were integrated to identify the important genes and lncRNAs. The morphological analysis showed that the size of the GCs did not significantly change, but the thickness of the granulosa layer cells differed significantly across the developmental stages. Subsequently, 14,286 mRNAs, 3,956 lncRNAs, and 1,329 TUCPs (transcripts with unknown coding potential) were detected in the GCs. We identified 37 common DEGs in the pre-hierarchical and hierarchical follicle stages, respectively, which might be critical for follicle development. Moreover, 3,089 significant time-course DEGs (Differentially expressed genes) and 13 core genes in 4 clusters were screened during goose GCs development. Finally, the network lncRNA G8399 with CADH5 and KLF2, and lncRNA G8399 with LARP6 and EOMES were found to be important for follicular development in GCs. Thus, the results would provide a rich resource for elucidating the reproductive biology of geese and accelerate the improvement of the egg-laying performance of geese.

Follicular development in livestock: Influencing factors and underlying mechanisms

Livestock farming development has become increasingly important in recent years. It not only provides us with meat nutrition and pet feeding but also increases the economic value by providing numerous employment opportunities, which improves our life quality. The livestock farming development depends on successful animal reproduction. As a vital process in animal reproduction, folliculogenesis and its influencing factors as well as their underlying mechanisms need to be understood thoroughly. This review is aimed at summarizing the factors such as cellular processes, gene regulation, noncoding RNAs and other endocrine or paracrine regulatory factors that affect follicular development, and their underlying mechanisms of action in livestock in order to provide novel insights for future studies. The above factors were found as significant determinants influencing the follicular development in livestock through various signaling pathways.

Transcriptome analysis reveals that long noncoding RNAs contribute to developmental differences between medium-sized ovarian follicles of Meishan and Duroc sows

Ovulation rate is an extremely important factor affecting litter size in sows. It differs greatly among pig breeds with different genetic backgrounds. Long non-coding RNAs (lncRNAs) can regulate follicle development, granulosa cell growth, and hormone secretion, which in turn can affect sow litter size. In this study, we identified 3554 lncRNAs and 25,491 mRNAs in M2 follicles of Meishan and Duroc sows. The lncRNA sequence and open reading frame lengths were shorter than mRNAs, and lncRNAs had fewer exons, were less abundant, and more conserved than protein-coding RNAs. Furthermore, 201 lncRNAs were differentially expressed (DE) between breeds, and quantitative trait loci analysis of DE lncRNAs were performed. A total of 127 DE lncRNAs were identified in 119 reproduction trait-related loci. In addition, the potential target genes of lncRNAs in cis or trans configurations were predicted. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that some potential target genes were involved in follicular development and hormone secretion-related biological processes or pathways, such as progesterone biosynthetic process, estrogen metabolic process, ovarian steroidogenesis, and PI3K-Akt signaling pathway. Furthermore, we also screened 19 differentially expressed lncRNAs in the PI3K-Akt signaling pathway as candidates. This study provides new insights into the roles of lncRNAs in follicular growth and development in pigs.

Dynamic Expression and Regulatory Network of Circular RNA for Abdominal Preadipocytes Differentiation in Chicken (Gallus gallus)

Circular RNA (circRNA), as a novel endogenous biomolecule, has been emergingly demonstrated to play crucial roles in mammalian lipid metabolism and obesity. However, little is known about their genome-wide identification, expression profile, and function in chicken adipogenesis. In present study, the adipogenic differentiation of chicken abdominal preadipocyte was successfully induced, and the regulatory functional circRNAs in chicken adipogenesis were identified from abdominal adipocytes at different differentiation stages using Ribo-Zero RNA-seq. A total of 1,068 circRNA candidates were identified and mostly derived from exons. Of these, 111 differentially expressed circRNAs (DE-circRNAs) were detected, characterized by stage-specific expression, and enriched in several lipid-related pathways, such as Hippo signaling pathway, mTOR signaling pathway. Through weighted gene co-expression network analyses (WGCNA) and K-means clustering analyses, two DE-circRNAs, Z:35565770|35568133 and Z:54674624|54755962, were identified as candidate regulatory circRNAs in chicken adipogenic differentiation. Z:35565770|35568133 might compete splicing with its parental gene, ABHD17B, owing to its strictly negative co-expression. We also constructed competing endogenous RNA (ceRNA) network based on DE-circRNA, DE-miRNA, DE-mRNAs, revealing that Z:54674624|54755962 might function as a ceRNA to regulate chicken adipogenic differentiation through the gga-miR-1635-AHR2/IRF1/MGAT3/ABCA1/AADAC and/or the novel_miR_232-STAT5A axis. Translation activity analysis showed that Z:35565770|35568133 and Z:54674624|54755962 have no protein-coding potential. These findings provide valuable evidence for a better understanding of the specific functions and molecular mechanisms of circRNAs underlying avian adipogenesis.

A genome-wide circular RNA transcriptome in rat

Circular RNAs (circRNAs) are a novel class of noncoding RNAs that back-splice from 5ʹ donor site and 3ʹ acceptor sites to form a circular structure. A number of circRNAs have been discovered in model organisms including human, mouse, Drosophila, among other organisms. There are a few candidate-based studies on circRNAs in rat, a well-studied model organism as well. A number of pipelines have been published to identify the back splice junctions for the discovery of circRNAs but studies comparing these tools have suggested that a combination of tools would be a better approach to identify high-confidence circRNAs. The availability of a recent dataset of transcriptomes encompassing 11 tissues, 4 developmental stages, and 2 genders motivated us to explore the landscape of circRNAs in the organism in this context. In order to understand the difference among different pipelines, we employed five different combinations of tools to identify circular RNAs from the dataset. We compared the results of the different combination of tools/pipelines with respect to alignment, total number of circRNAs identified and read-coverage. In addition, we identified tissue-specific, development-stage specific and gender-specific circRNAs and further independently validated 16 circRNA junctions out of 24 selected candidates in 5 tissue samples and estimated the quantitative expression of five circRNA candidates using real-time polymerase chain reaction and our analysis suggests three candidates as tissue-enriched. This study is one of the most comprehensive studies which provides a map of circRNAs transcriptome as well as to understand the difference among different computational pipelines in rat.

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.

circ_C20orf11 enhances DDP resistance by inhibiting miR-527/YWHAZ through the promotion of extracellular vesicle-mediated macrophage M2 polarization in ovarian cancer

Ovarian cancer is a fatal gynecologic tumor, and conventional treatment is mainly limited by chemoresistance. The mechanism contributing to chemoresistance in ovarian cancer has yet to be established. This study aimed to investigate the specific role of circ_C20orf11 in regulating chemoresistance to cisplatin (DDP)in ovarian cancer. We first established two DDP-resistant ovarian cancer cell lines. Then, we identified the effect of circ_C20orf11 on specific cellular characteristics (proliferation, apoptosis, DDP resistance) via a series of experiments. The binding sites between circ_C20orf11 and miR-527 and between miR-527 and YWHAZ were predicted using a bioinformatics tool and confirmed with a dual-luciferase reporter assay. Furthermore, extracellular vesicles (EVs) derived from DDP-resistant cell lines were identified, and the effect of EVs on macrophage polarization was examined. circ_C20orf11 was upregulated in ovarian cancer. Increased circ_C20orf11 expression enhanced DDP resistance and cell proliferation and reduced cell apoptosis in DDP-resistant cell lines after DDP treatment by sponging miR-527 and promoting YWHAZ expression. In addition, we found that DDP-resistant cell-derived EVs can induce macrophage M2 polarization, whereas silencing of circ_C20orf11 inhibited EV-induced macrophage M2 polarization. Consistent with these results, silencing of circ_C20orf11 enhanced sensitivity to DDP in vivo. Importantly, we proved that circ_C20orf11 expression was upregulated in EVs extracted from the serum of DDP-resistant patients. Our study demonstrated that silencing circ_C20orf11 sensitizes ovarian cancer to DDP by promoting miR-527/YWHAZ signaling and EV-mediated macrophage M2 polarization.

Dynamic expression and functional analysis of circular RNA in the gonads of Chinese soft-shelled turtles (Pelodiscus sinensis)

Circular RNA (circRNA) is a noncoding RNA that can regulate a variety of biological processes. CircRNAs can regulate gene expression posttranscriptionally by acting as microRNA sponges. Many turtle species are remarkable organisms due to their reproductive processes. However, information on circRNA in the gonads of turtles is limited. In this study, 6, 121 circRNAs were identified in the testes and ovaries of Chinese soft-shelled turtles (Pelodiscus sinensis) using the Illumina platform, and 710 circRNAs were significantly differentially expressed (DE). The DE circRNAs included 541 upregulated and 169 downregulated circRNAs in the testes. GO and KEGG pathway analysis indicated that the DE circRNAs were enriched in several signaling pathways, including GnRH, Wnt, FoxO, Progesterone mediated oocyte maturation, and mTOR signaling pathways. Five DE circRNAs were randomly selected, and their relative expression levels in ovaries and testes were detected by quantitative real-time PCR. All of these circRNAs were differentially expressed. In addition, 9, 883 interactions between circRNAs and miRNAs were predicted in the turtles. Target genes of the miRNAs include a range of genes regulating gonadal development. Seven ceRNA networks (DE circRNAs-DE miRNAs-DE mRNAs), including 7 DE circRNAs, 11 DE miRNAs and 20 DE mRNAs, were constructed. The networks included Cdc6, the miR-1 family, the miR-203 family, and the miR-302 family. The expression profile of gonadal circRNAs might help to elucidate the roles of nonprotein coding RNAs in turtle gonadal development.

Dynamic Expression Profiles of Circular RNAs during Brown to White Adipose Tissue Transformation in Goats (Capra hircus)

Simple Summary

In our study, we launched RNA-seq in order to investigate the potential functions of circRNA during brown adipose tissue (BAT) to white adipose tissue (WAT) transformation. As a result, 6610 circRNAs and 61 differentially expressed circRNAs (DEcircRNAs) were identified. Moreover, 65 miRNAs were detected that could potentially interact with DEcircRNAs. The present study provides a detailed circRNA expression landscape and evidence for circRNA functions in the transformation from BAT to WAT.

Abstract

Adipose tissues are mainly divided into brown adipose tissue (BAT) and white adipose tissue (WAT). WAT mainly functions to buffer excess calories, whereas BAT plays a role in the non-shivering thermogenesis to maintain body temperature and energy balance. Moreover, circRNAs play important roles in various biological processes. However, knowledge of the expression profile and function of circRNAs from BAT to WAT remains largely unknown. In this study, a total of 6610 unique circRNAs were identified in the perirenal adipose tissues of 1-day, 30-days, and 1-year goats. Functional annotation revealed that host genes of circRNAs were involved in some BAT-related pathways, such as the thyroid hormone signaling pathway, MAPK signaling pathway, and VEGF signaling pathway. Furthermore, a total of 61 DEcircRNAs were detected across three stages. Additionally, five selected circRNAs were validated by RNase R assay, qPCR, and Sanger sequencing. Finally, the circRNA–miRNA network was constructed between the DEcircRNAs and their miRNA binding sites.

Circular RNAs: Novel potential regulators in embryogenesis, female infertility, and pregnancy-related diseases

Circular RNAs (circRNAs) are endogenous noncoding RNAs with unique cyclic structures. Although they were previously considered as nonfunctional transcription byproducts, numerous studies have demonstrated that circRNAs regulate gene transcription and expression via different mechanisms. Reproductive health influences the quality of life and affects offspring propagation in women. CircRNAs have been found to modify pregnancy-related diseases, gynecologic cancers, polycystic ovary syndrome, aging, gamete, and embryo development. It's promising for circRNAs to be the novel diagnostic and therapeutic targets for multiple reproductive diseases. With the widespread application of assisted reproduction technology (ART), it has been revealed that circRNA identification contributes to estimating the quality of gametes and embryos, reflecting the success rate of ART. CRISPR-Cas9 gene editing technology has enabled the discovery of new roles of circRNAs. So far, the roles of circRNAs in the reproductive system remain poorly defined. In this review, we describe the classification and functions of circRNAs in embryogenesis and the female reproductive system diseases, revealing potential roles of circRNAs physiologically and pathologically. In so-doing, we provide ideas for developing circRNA-based therapeutic treatment and clinical application of various female reproductive system diseases.

Characteristics and expression profiles of circRNAs during abdominal adipose tissue development in Chinese Gushi chickens

Circular RNAs (circRNAs) play important roles in adipogenesis. However, studies on circRNA expression profiles associated with the development of abdominal adipose tissue are lacking in chickens. In this study, 12 cDNA libraries were constructed from the abdominal adipose tissue of Chinese domestic Gushi chickens at 6, 14, 22, and 30 weeks. A total of 1,766 circRNAs were identified by Illumina HiSeq 2500 sequencing. These circRNAs were primarily distributed on chr1 through chr10 and sex chromosomes, and 84.95% of the circRNAs were from gene exons. Bioinformatic analysis showed that each circRNA has 35 miRNA binding sites on average, and 62.71% have internal ribosome entry site (IRES) elements. Meanwhile, these circRNAs were primarily concentrated in TPM < 0.1 and TPM > 60, and their numbers accounted for 18.90% and 80.51%, respectively, exhibiting specific expression patterns in chicken abdominal adipose tissue. In addition, 275 differentially expressed (DE) circRNAs were identified by comparison analysis. Functional enrichment analysis showed that the parental genes of DE circRNAs were primarily involved in biological processes and pathways related to lipid metabolism, such as regulation of fat cell differentiation, fatty acid homeostasis, and triglyceride homeostasis, as well as fatty acid biosynthesis, fatty acid metabolism, and glycerolipid metabolism. Furthermore, ceRNA regulatory networks related to abdominal adipose development were constructed. The results of this study indicated that circRNAs can regulate lipid metabolism, adipocyte proliferation and differentiation, and cell junctions during abdominal adipose tissue development in chickens through complex ceRNA networks between circRNAs, miRNAs, genes, and pathways. The results of this study may help to expand the number of known circRNAs in abdominal adipose tissue and provide a valuable resource for further research on the function of circRNAs in chicken abdominal adipose tissue.

LncRNA lnc_13814 promotes the cells apoptosis in granulosa cells of duck by acting as apla-miR-145-4 sponge

Follicle development is a vital factor which determines the reproductive performance of poultry. Long noncoding RNAs (lncRNAs) have been reported to maintain animal reproductive function and play key roles in ovarian development and hormone secretion. But the regulatory mechanism of lncRNAs in duck follicle development has seldom been reported. In this study, to better explore the molecular mechanism of follicle development in ducks, the follicular lncRNA was sequenced and analyzed. A total of 9,551 lncRNAs were predicted in the duck follicles. Four hundred and forty-five lncRNAs were differentially expressed between the white follicles and yellow follicles. The results of our studies showed that lnc_13814 promoted cell apoptosis in duck GCs. Furthermore, the bioinformatics analysis results demonstrated that lnc_13814 was involved in a lncRNA-miRNA-mRNA coexpression network and it was observed to sponge two follicle-related miRNAs by a luciferase activity assay. Moreover, we found that overexpression of lnc_13814 significantly increased DNA damage inducible transcript 3 (DDIT3) expression and downregulated GCs apoptosis. Finally, we found that lnc_13814 directly binds to and inhibits apla-mir-145-4; then, lnc_13814 increases the expression of DDIT3 and up-regulates GCs apoptosis. Taken together, our findings demonstrate that lncRNAs have potential effects on duck ovarian follicles and lncRNAs may represent a new approach to understand follicular development.

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.

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.

The mediation of pigeon egg production by regulating the steroid hormone biosynthesis of pigeon ovarian granulosa cells

The aim of this study was to determine the molecular mechanism of miR-205b targeting 11β-hydroxysteroid dehydrogenase type 1 (HSD11B1) on the apoptosis and proliferation of granulosa cells (GC) of pigeons. Our previous studies suggested that HSD11B1 was the target gene of miR-205b and played a key role in steroid hormone biosynthesis and GC development. The adenovirus-miR-205b recombinant virus and adenovirus-cli-miR-205b-sh recombinant virus were generated, verified, and their characteristics determined. The recombinant viruses were used to infect the GC of pigeons, with real time quantitative PCR used to examine the expressions of HSD11B1 and related genes. The HSD11B1 antibody was obtained and verified, and Western blotting was used to detect the protein level of HSD11B1. The Cell Counting Kit-8 assay kit was used to detect cell viability, and the Annexin V-FITC/PI kit was used for the apoptosis assays. The expression of HSD11B1 was significantly lower in the overexpression (OE) than in OE negative control (OE-NC) treatments and significantly higher in short hairpin (SH) than in SH negative control (SH-NC) treatments. The expression levels of cytochrome P4503A5 was significantly higher in SH and lower in OE treatments, and the rhythms of cytochrome P450 aromatase mRNA levels were similar. The mRNA level of cytochrome P450scc in OE was lower than in OE-NC treatments and higher in SH than in SH-NC treatments. The protein expressions of HSD11B1 were decreased in the GC of OE, whereas increased in the SH group. The Cell Counting Kit-8 assay revealed that overexpression of miR-205b significantly suppressed proliferation of the GC of pigeons, whereas interference of miR-205b significantly induced the proliferation of the GC. The overexpression and the interference of miR-205b did not have a significant effect on cell cycle. The overexpression of miR-205b significantly increased the number of apoptotic cells, whereas the interference of miR-205b decreased the number of apoptotic cells. These findings indicated that miR-205b mediated pigeon egg production by regulating the steroid hormone biosynthesis of the pigeon ovarian GC by targeting HSD11B1, which may be useful in increasing pigeon egg production.

lncRNA H19 acts as a ceRNA to regulate the expression of CTGF by targeting miR-19b in polycystic ovary syndrome

The etiology of polycystic ovary syndrome (PCOS) is complex and the pathogenesis is not fully understood. Some studies have shown that dysregulation of ovarian granulosa cells may be related to abnormal follicles and excessive androgen in women with PCOS. Our team has also confirmed the high expression status of H19 in PCOS patients in the early stage. However, the relationship between H19 and miR-19b in the development of PCOS is still unknown. Therefore, we used bioinformatics to predict the binding sites of human H19 and miR-19b, and of miR-19b and CTGF genes. After the silencing and overexpression of H19, real-time polymerase chain reaction (PCR) was used to detect the expressions of H19, miR-19b, and CTGF. Western blotting was used to detect CTGF protein. Proliferation of KGN cells after H19 silencing was detected by CCK8. Flow cytometry was used to detect the apoptosis of KGN cells after H19 silencing. After the overexpression of H19, it was found that the expression of miR-19b gene decreased and the expression of CTGF increased, whereas silencing of H19 did the opposite. In addition, H19 could promote cell proliferation and decrease cell apoptosis. Finally, luciferase reporter assays showed that the 3′-end sequences of lncRNA H19 and CTGF contained the binding site of miR-19b. In conclusion, our study indicated that lncRNA H19 acted as a ceRNA to bind to miR-19b via a “sponge” to regulate the effect of CTGF on KGN cells, which may play a vital role in PCOS.

Integrative analysis of circRNAs, miRNAs, and mRNAs profiles to reveal ceRNAs networks in chicken intramuscular and abdominal adipogenesis

BACKGROUND

Tissue-specific fat deposition is regulated by a series of complex regulatory mechanisms. Reports indicate that epigenetic regulators, such as circular RNAs (circRNAs), are crucial in diseases progression, animal development, metabolism, and adipogenesis. In this study, to assess the functional roles of circRNAs in adipogenesis and tissue-specific fat deposition, we comprehensively analyzed the Ribo-Zero RNA-Seq and miRNAs data during chicken intramuscular and abdominal adipogenic differentiation.

RESULTS

circRNAs and miRNAs profiles during chicken adipogenic differentiation were found in adipocytes derived from various adipose tissues. It was also discovered that high levels of downregulated miRNAs potentially promote adipogenesis by activating their target genes which are associated with fatty acid metabolism and adipogenic differentiation. Through analysis of the correlation between the expression levels of circRNAs and adipogenic genes, as well as the dynamic expression patterns of circRNAs during adipogenic differentiation, several candidate circRNAs were identified. Moreover, competing endogenous RNA (ceRNAs) networks were constructed during chicken intramuscular and abdominal adipogenesis by combining miRNAs with mRNAs data. Several candidate circRNAs potentially influence adipogenesis by regulating miRNAs via PPAR and fatty acid metabolism-related pathways were identified, such as circLCLAT1, circFNDC3AL, circCLEC19A and circARMH1.

CONCLUSION

In conclusion, our findings reveal that circRNAs and the circRNA-miRNAs-mRNAs-ceRNAs network may play important roles in chicken adipocytes differentiation and tissue-specific fat deposition.

Ovarian transcriptomic analysis and follicular development of Leizhou black duck

This study investigated the factors that caused the differences in egg production during the development of ovarian follicles in Leizhou black ducks. Leizhou black ducks population was divided into 2 groups as high-yield group (HG) and low-yield group (LG). The number of eggs (NE), age at first egg (AFE), weight at first egg, and egg weight (EW) of both groups were recorded, and differences were analyzed using the t test. The logistic model was used to simulate the egg production curves to analyze the production rules. The ovarian follicles of both duck groups were collected to count the number of different grades sized follicles, weigh the ovaries, and observe follicular sections to analyze the developmental differences. Ovarian transcriptomic sequencing was performed to investigate differentially expressed genes and signal pathways in both duck groups. The results revealed a significant difference (P < 0.01) in the NE laid, AFE, and EW between both groups. Comparatively, HG had significantly more (P < 0.01) large yellow follicles (LYF) than LG. The density of medullary layer cells of the follicle section was greater in HG than LG ducks. Transcriptome sequencing revealed a total of 1,027 differentially expressed genes between the HG and LG ducks of which 495 genes were upregulated, and 532 genes were downregulated. Fifty genes were related to reproduction and reproductive processes. Kyoto Encyclopedia of Genes and Genomes–enriched signaling pathways revealed 274 signal pathways enriched in these differentially expressed genes of which the steroid biosynthesis pathway was significantly enriched. Analysis (Q < 0.05) showed that HSD3β → gonadotropin-releasing hormone (GnRH) and estrogen receptor (ESR) → LHβ/ERK1/2 were enriched in the steroid biosynthesis signal pathway. Follicle-stimulating hormone signal pathway mediated by HSD3β → GnRH and ESR → LHβ/ERK1/2 may be involved in ovarian follicle development to regulate LYF reserve process and affect its ovulation cycle, which in turn influence the egg production of Leizhou black ducks.

Transcriptome Analysis of circRNA and mRNA in Theca Cells during Follicular Development in Chickens

Development of ovarian follicles requires interactions between granulosa cells, theca cells, and oocytes. Multiple transcription levels are involved but information about the role of noncoding RNAs, especially circular RNAs (circRNAs), is lacking. Here, we used RNA sequencing to profile circRNAs and mRNAs in theca cells from three types of follicle: small yellow follicles (SYF), the smallest hierarchical follicles (F6), and the largest hierarchical follicles (F1). Using bioinformatics analysis, we identified a total of 14,502 circRNAs in all theca cells, with 5622 widely distributed in all stages of development. Differential expression analysis suggested that some genes display differential isoforms during follicular development. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed enrichment of both differentially expressed circRNAs and mRNAs in pathways associated with reproduction, including the TGF-β signaling pathway, oocyte meiosis, and vascular smooth muscle contraction. Our study provides the first visual information about circRNAs and mRNAs in theca cells during follicle development in chickens and adds to the growing body of knowledge about theca cells.