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

LH/hCG Regulation of Circular RNA in Mural Granulosa Cells during the Periovulatory Period in Mice

Ovarian follicles undergo a series of dynamic changes following the ovulatory surge of luteinizing hormone including cumulus expansion, oocyte maturation, ovulation, and luteinization. Post-transcriptional gene regulatory events are critical for mediating LH follicular responses, and among all RNA isoforms, circular RNA (circRNA) is one of the most abundant forms present in cells, yet they remain the least studied. Functionally, circRNA can act as miRNA sponges, protein sponges/decoys, and regulators of transcription and translation. In the context of ovarian follicular development, the identity and roles of circRNA are relatively unknown. In the present study, high throughput RNA sequencing of granulosa cells immediately prior to and 4-h after the LH/hCG surge identified 42,381 circRNA originating from 7712 genes. A total of 54 circRNA were identified as differentially expressed between 0-h and 4-h time points (Fold Change ± 1.5, FDR ≤ 0.1), among them 42 circRNA were upregulated and 12 circRNA were downregulated. All differentially expressed circRNA between the 0-h and 4-h groups were subjected to circinteractome analysis and identified networks of circRNA-protein and circRNA-miRNA were further subjected to "micro-RNA target filter analysis" in Ingenuity Pathway Analyses, which resulted in the identification of miRNA targeted mRNAs. A comparison of these circRNA target mRNAs with LH-induced mRNAs identified Runx2, Egfr, Areg, Sult1el, Cyp19a1, Cyp11a1, and Hsd17b1 as targets of circKif2, circVcan, circMast4, and circMIIt10. These newly identified LH/hCG-induced circRNA, their target miRNA and protein networks provide new insights into the complex interactions associated with periovulatory follicular development.

Differential expression and functional analysis of circRNA in the ovaries of Yili geese at different egg-laying stages

BACKGROUND

Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. Circular RNA (circRNA) has been proven to play a significant role in the regulation of animal reproduction.

OBJECTIVE

Studying the expression characteristics of circRNAs in goose ovaries at various egg-laying stages can provide a reference for the molecular regulation mechanism of ovary development in geese that is mediated by circRNAs.

METHODS

In this study, the expression profiles of circRNAs were compared in ovary tissues from Yili geese in three different breeding periods, namely the prelaying period (KL), laying period (CL), and ceased period (XL), and differentially expressed circRNAs related to ovarian development in Yili geese were screened. The potential biological functions of differential circRNAs were predicted by bioinformatics, and the differential circRNA-miRNA regulatory network was constructed.

RESULTS

The results showed that a total of 4483 circRNAs were identified in 12 ovarian tissue samples from Yili geese at different laying stages. In the KL vs. CL, XL vs. CL, and XL vs. KL groups, 159, 455, and 383 differentially expressed circRNAs were identified, respectively. The host genes of the differential circRNAs were mostly enriched in the signal transduction, metabolism, and other related pathways, such as those for phototransduction, glycerophospholipid metabolism, aminoacyl-tRNA biosynthesis, and retinol metabolism. Finally, we constructed circRNA-miRNA regulation networks. Nine differential circRNAs were randomly selected for qRT-PCR verification, and the expression trends were consistent with the sequencing results.

CONCLUSIONS

Our results indicated that significant differences in the expression profiles of circRNAs in the ovaries of Yili geese at different egg-laying stages. Meanwhile, through analyzing the differential circRNA-miRNA interaction network, core regulators such as circRNA NW_013186107.1:36835|52,574 and gga-miR-34b-5p were screened. This study provides a reference for the further analysis of the molecular regulatory mechanism of the circRNAs regulating goose ovary development and enriches the theory of genetic regulation during goose ovary development.

CircRNAs in Xiang pig ovaries among diestrus and estrus stages

BACKGROUND

The fecundity of sows is a trait of major economic in pig industry. The molecular regulation of estrus cycles can affect the fecundity of female animals. Compared with the other pig breeds, Xiang pig exhibits the special estrus behaviors. CircRNAs are thought to involve in regulation of multiple biological processes. However, the potential roles of circRNAs in ovary regulation on Xiang pig estrus are largely unknown.

RESULTS

8,937 circRNAs were identified from eight libraries constructed from the ovarian samples of Xiang pig at estrus and diestrus stages by RNA sequencing method. Of which, 1,995 were high confidence circRNAs detected at least two junction reads in each ovary sample and seven circRNAs were validated by RT-PCR method. Furthermore, we identified 290 upregulated and 15 downregulated circRNAs in estrus ovaries. These differentially expressed circRNAs (DECs) derived from 273 host genes. And 207 miRNAs were identified to be targets sponged by 156 DECs with 432 binding sites, containing more than one miRNA binding site in each circRNA. Function enrichment analysis revealed that the host genes and the targets of miRNAs sponged by DECs were enriched in several reproduction-related signaling pathways, such as ovarian steroidogenesis, oocyte maturation, circadian rhythm, estrogen signaling pathway, GnRH signaling pathway, circadian entrainment, and oocyte meiosis. The circRNA-miRNA-mRNA networks revealed that 153 miRNAs interacting with 122 DECs and 86 miRNAs interacting with 84 DECs were involved in ovarian functions and ovarian circadian entrainment and circadian rhythm respectively. The DEC-miRNA-DEG (differentially expressed gene, DEG) networks associated with reproduction-related signaling pathways contained 22 DECs,18 miRNAs and 7 DEGs. 22 DECs were recognized as hub circRNAs during the estrus phase of Xiang pigs.

CONCLUSIONS

The circRNAs that function as miRNA sponges could play a key role in post-transcriptional regulation of gene expression during Xiang pig's estrus cycle.

Integrated analysis of lncRNA, miRNA and mRNA reveals novel insights into the fertility regulation of large white sows

Background

Improving sow fertility is extremely important as it can lead to increased reproductive efficiency and thus profitability for swine producers. There are considerable differences in fertility rates among individual animals, but the underlying molecular mechanisms remain unclear. In this study, by using different types of RNA libraries, we investigated the complete transcriptome of ovarian tissue during the luteal (L) and follicular (F) phases of the estrous cycle in Large White pigs with high (H) and low (L) fecundity, and performed a comprehensive analysis of long noncoding RNAs (lncRNAs), mRNAs and micro RNAs (miRNAs) from 16 samples by combining RNA sequencing (RNA-seq) with bioinformatics.

Results

In total, 24,447 lncRNAs, 27,370 mRNAs, and 216 known miRNAs were identified in ovarian tissues. The genomic features of lncRNAs, such as length distribution and number of exons, were further analyzed. We selected a threshold of P < 0.05 and |log₂ (fold change)| ≥ 1 to obtain the differentially expressed lncRNAs, miRNAs and mRNAs by pairwise comparison (LH vs. LL, FH vs. FL). Bioinformatics analysis of these differentially expressed RNAs revealed multiple significantly enriched pathways (P < 0.05) that were closely involved in the reproductive process, such as ovarian steroidogenesis, lysosome, steroid biosynthesis, and the estrogen and GnRH signaling pathways. Moreover, bioinformatics screening of differentially expressed miRNAs that share common miRNA response elements (MREs) with lncRNAs and their downstream mRNA targets were performed. Finally, we constructed lncRNA–miRNA–mRNA regulation networks. The key genes in these networks were verified by Reverse Transcription Real-time Quantitative PCR (RT-qRCR), which were consistent with the results from RNA-Seq data.

Conclusions

These results provide further insights into the fertility of pigs andcan contribute to further experimental investigation of the functions of these genes.