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

Identification and expression profile analysis of circRNAs associated with goat uterus with different fecundity during estrous cycle

BACKGROUND

The Yunshang Black Goat, a distinguished meat goat breed native to China, is renowned for its superior reproductive capabilities. Despite this, there is considerable phenotypic variability within the breed. During the reproductive cycle, the uterus plays a pivotal role, with its functions evolving in line with the different stages of the cycle. This study focuses on the uterine tissues, including both the endometrium and myometrium, of Yunshang Black Goats with high fecundity (HF) and low fecundity (LF) during the proliferative (FP) and secretory (LP) phases of the estrous cycle. By examining these tissues, we aim to elucidate the underlying molecular and physiological mechanisms of the observed differences in reproductive success.

RESULTS

High-throughput sequencing was conducted, followed by bioinformatics analysis to identify the expression profiles of circRNAs. A total of 7,445 circRNAs were identified through the integration of findings from find_circ and CIRI2 software. Comparative analyses between the FPLF vs. FPHF and LPLF vs. LPHF revealed 149 differentially expressed (DE) circRNAs (94 up-regulated and 55 down-regulated) and 276 DE circRNAs (56 up-regulated and 220 down-regulated), respectively. The enrichment analysis indicated that the primary pathways involved were the Sphingolipid signaling pathway, MAPK signaling pathway, and GnRH signaling pathway, all of which are closely associated with cellular growth and development. Additionally, several key candidate genes were identified, such as FGF2 and MBTPS1. We also predicted a total of 281 miRNA-circRNA binding pairs, encompassing 263 circRNAs and 60 miRNAs, and simultaneously, 14 coding circRNAs were anticipated.

CONCLUSION

Based on the analysis, we have established the expression profiles of circRNAs during the follicular and luteal phases, respectively. Furthermore, using various analytical methods and data from high- and low-yield experimental control groups over different periods, we have identified multiple circRNAs that affect the high reproductive capacity of goats. Through enrichment analysis of the host genes of these circRNAs, we have discovered several key candidate genes. These findings provide fundamental data for the study of the molecular mechanisms underlying the fecundity of goats and pave the way for future genetic improvement strategies.

CircMYH9/miR-133a-3p/CXCR4 axis: a novel regulatory network in sperm fertilization and embryo development

This study aimed to investigate the influence of sperm miRNAs on fertilization rates (FR) in in vitro fertilization (IVF) and to explore potential regulatory mechanisms in sperm-mediated fertilization and embryo development. Through high-throughput sequencing, we identified differentially expressed miRNAs in sperm, with miR-133a-3p significantly upregulated in samples associated with low FR and available embryo rate (AER). Key regulatory circRNAs and mRNAs were further identified via the Starbase database, intersected with differentially expressed RNA, and analyzed through GO, KEGG, and PPI analyses. The circMYH9/miR-133a-3p/CXCR4 axis emerged as a critical regulatory network. In vitro assays using the GC-2 spd mouse spermatogenic cell line revealed that miR-133a-3p inhibited cell growth and proliferation while promoting apoptosis. circMYH9, acting as a competing endogenous RNA (ceRNA) for miR-133a-3p, modulated CXCR4 expression, enhancing GC-2 spd cell growth and inhibiting apoptosis through the miR-133a-3p/CXCR4 axis. In vivo experiments using a mouse model confirmed that circMYH9 overexpression increased IVF success rates and promoted embryo development via this axis. Mechanistically, miR-133a-3p suppresses sperm fertilization and embryo development by targeting the circMYH9/miR-133a-3p/CXCR4 axis. These findings suggest that this regulatory network could serve as a novel biomarker for assessing fertilization potential and embryo quality in clinical settings and as a potential therapeutic target to improve IVF outcomes and address infertility. This study provides valuable insights into the molecular mechanisms governing sperm function and early embryonic development.

Identification of functional circRNAs regulating ovarian follicle development in goats

BARKGROUND

Circular RNAs (circRNAs) play important regulatory roles in a variety of biological processes in mammals. Multiple birth-traits in goats are affected by several factors, but the expression and function of circRNAs in follicular development of goats are not clear. In this study, we aimed to investigate the possible regulatory mechanisms of circRNA and collected five groups of large follicles (Follicle diameter > 6 mm) and small follicles (1 mm < Follicle diameter < 3 mm) from Leizhou goats in estrus for RNA sequencing.

RESULTS

RNA sequencing showed that 152 circRNAs were differentially expressed in small and large follicles. Among them, 101 circRNAs were up-regulated in large follicles and 51 circRNAs were up-regulated in small follicles. GO and KEGG enrichment analyses showed that parental genes of the differential circRNAs were significantly enriched in important pathways, such as ovarian steroidogenesis, GnRH signaling pathway, animal autophagy and oxytocin signalling pathway. BioSignal analysis revealed that 152 differentially expressed circRNAs could target 91 differential miRNAs including miR-101 family (chi-miR-101-3p, chi-miR-101-5p), miR-202 family (chi-miR-202-5p, chi-miR-202-3p),60 circRNAs with translation potential. Based on the predicted sequencing results, the ceRNA networks chicirc_008762/chi-miR-338-3p/ARHGAP18 and chicirc_040444/chi-miR-338-3p/STAR were constructed in this study. Importantly, the new gene circCFAP20DC was first discovered in goats. The EDU assay and flow cytometry results indicated that circCFAP20DC enhanced the proliferation of follicular granulosa cells(GCs). Real-time quantitative PCR and western blotting assays showed that circCFAP20DC activated the Retinoblastoma(RB) pathway and promoted the progression of granulosa cells from G1 to S phase.

CONCLUSION

Differential circRNAs in goat size follicles may have important biological functions for follicular development. The novel gene circCFAP20DC activates the RB pathway, promoting the progression of GCs from G1 to S phase. This, in turn, enhances the proliferation of follicular GCs in goats.

RNA-Seq Analysis Revealed circRNAs Associated with Resveratrol-Induced Apoptosis of Porcine Ovarian Granulosa Cells

Circular RNAs (circRNAs) are a class of circular non-coding RNAs that play essential roles in the intricate and dynamic networks governing cell growth, development, and apoptosis. Resveratrol (RSV), a non-flavonoid polyphenol, is known to participate in follicular development and ovulation. In our previous research, we established a model using porcine ovarian granulosa cells (POGCs) treated with resveratrol, which confirmed its regulatory effects on long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) within these cells. However, the influence of resveratrol on circRNA expression has not been thoroughly investigated. To explore how resveratrol affects circRNA levels in POGCs, we designed an experiment with three groups: a control group (CON, n = 3, 0 μM RSV), a low-dose RSV group (LOW, n = 3, 50 μM RSV), and a high-dose RSV group (HIGH, n = 3, 100 μM RSV) for circRNA sequencing. We identified a total of 10,045 candidate circRNAs from POGCs treated with different concentrations of resveratrol (0, 50, and 100 μM). Differential expression analysis indicated that 96 circRNAs were significantly altered in the LOW vs. CON group, while 109 circRNAs showed significant changes in the HIGH vs. CON group. These circRNAs were notably enriched in biological processes associated with cell metabolism, apoptosis, and oxidative stress. Functional enrichment analysis of the host genes revealed their involvement in critical signaling pathways, including mTOR, AMPK, and apoptosis pathways. Additionally, we identified potential miRNA sponge candidates among the differentially expressed circRNAs, particularly novel_circ_0012954 and novel_circ_0004762, which exhibited strong connectivity within miRNA-target networks. Our findings provide valuable insights into the regulatory mechanisms of circRNAs in the context of resveratrol-induced apoptosis in POGCs, highlighting their potential as innovative therapeutic targets in reproductive biology.

Transcriptome Analysis Reveals Differentially Expressed circRNAs Associated with Fecundity in Small-Tail Han Sheep Thyroid with Different FecB Genotypes

Litter size is an economically important trait in sheep, and it is a complex trait controlled by multiple genes in multiple organs. Among them, the regulation of lamb number trait by the thyroid gland is a very important part. However, the molecular mechanisms of the thyroid gland in sheep reproduction remain unclear. Here, RNA-seq was used to detect transcriptome expression patterns in the thyroid gland between follicular phase (FP) and luteal phase (LP) in FecB BB (MM) and FecB ++ (ww) STH sheep, respectively, and to identify differentially expressed circRNAs (DECs) associated with reproduction. Bioinformatic analysis of the source genes of these DECs revealed that they can be enriched in multiple signaling pathways involved in the reproductive process of animals. We found that the source genes of these DECs, such as GNAQ, VEGFC, MAPK1, STAT1, and HSD17B7, may play important roles in the reproductive process of animals. To better understand the function of these DECs, we constructed circRNA-miRNA co-expression networks. Dual luciferase reporter assays suggested that a ceRNA regulatory mechanism between circ_0003259-oar-miR-133-TXLNA and circ_0012128-oar-miR-370-3p-FGFR1 may hold. All of these DEC expression profiles in the thyroid gland provide a novel resource for elucidating the regulatory mechanisms underlying STH sheep prolificacy.

Thyroid Transcriptomics Revealed the Reproductive Regulation of miRNA in the Follicular and Luteal Phases in Small-Tail Han Sheep with Different FecB Genotypes

MicroRNA (miRNA) is a type of endogenous short-stranded ncRNA that influences many biological processes such as animal growth, development and metabolism. The thyroid gland is an important endocrine gland in sheep, and an increasing number of studies have shown that the thyroid gland plays an important role in animal reproduction, but the molecular mechanisms of the thyroid gland in sheep reproduction are poorly understood. In this study, RNA-seq was used to detect transcriptome expression patterns in the thyroid gland between the follicular phase (FP) and luteal phase (LP) in FecB BB (MM) and FecB ++ (ww) small-tail Han (STH) sheep, respectively, and to identify differentially expressed miRNAs (DEMs) associated with reproduction. Bioinformatic analysis of the target genes of these DEMs revealed that they can be enriched in multiple GO terms associated with the reproductive process in animals and in the KEGG signaling pathway. The miRNA-mRNA coexpression network revealed that oar-miR-133 and oar-miR-370-3p may play an important role in sheep reproduction. The results of the dual-luciferase reporter assay suggest a possible targeting relationship between novel-51 and TARBP2. These results provided a novel resource for elucidating regulatory mechanisms underlying STH sheep prolificacy.

Comparative proteomics of ovaries elucidated the potential targets related to ovine prolificacy

Small Tail Han (STH) sheep, a unique Chinese breed, is recognized for its early maturity, year-round estrus, and prolificacy. However, the molecular mechanism of its high prolificacy has not been fully elucidated. The Proteomics approach is feasible and effective to reveal the proteins involved in the complex physiological processes of any organism. Given this, we performed the protein expression profiling of ovarian tissues during the luteal phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PL vs. ML), and the follicular phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PF vs. MF), respectively. Parallel Reaction Monitoring (PRM) was conducted to validate the differentially abundant proteins (DAPs). The tandem mass tag (TMT) quantitative proteomic results showed that a total of 5,237 proteins were identified, of which 49 and 44 showed differential abundance in the PL vs. ML and PF vs. MF groups, respectively. Enrichments analyses indicated that the DAPs including TIA1 cytotoxic granule-associated RNA-binding protein-like 1 (TIAL1), nicotinamide phosphoribosyltransferase (NAMPT), and cellular retinoic acid-binding protein 1 (CRABP1) were enriched at the luteal phase, while TIAL1, inhibin beta-a-subunit (A2ICA4), and W5PG55 were enriched at the follicular phase, potentially mediating reproductive processes in polytocous ewes. Furthermore, six DAPs were verified using PRM, confirming the accuracy of the TMT data acquired in this study. Together, our work expanded the database of indigenous sheep breeds and provided new ovarian candidate molecular targets, which will help in the study of the genetic mechanisms of ovine prolificacy.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

The Role of Circular RNAs in the Physiology and Pathology of the Mammalian Ovary

Circular RNAs (circRNAs) are an abundant class of endogenous non-coding RNAs (ncRNAs) generated from exonic, intronic, or untranslated regions of protein-coding genes or intergenic regions. The diverse, stable, and specific expression patterns of circRNAs and their possible functions through cis/trans regulation and protein-coding mechanisms make circRNA a research hotspot in various biological and pathological processes. It also shows practical value as biomarkers, diagnostic indicators, and therapeutic targets. This review summarized the characteristics, classification, biogenesis and elimination, detection and confirmation, and functions of circRNAs. We focused on research advances circRNAs in the mammalian ovary under conditions including ovarian cancer, polycystic ovarian syndrome (PCOS), and maternal aging, as well as during reproductive status, including ovarian follicle development and atresia. The roles of circRNAs in high reproductive traits in domestic animals were also summarized. Finally, we outlined some obstructive factors and prospects to work with circRNA, aiming to provide insights into the functional research interests of circRNAs in the reproduction and gynecology areas.

Characterization of circular RNA profiles of oviduct reveal the potential mechanism in prolificacy trait of goat in the estrus cycle

The mammalian oviduct is functionally highly diverse during the estrus cycle. It provides a suitable milieu for oocyte maturation, sperm capacitation, fertilization, early embryo development and transportation. While there have been many studies of molecular mechanisms on the kidding number of goats, a systematic analysis by which the underlying circular RNAs (circRNAs) changes in the oviduct related to prolificacy traits is lacking. Herein, we present a comprehensive circRNA atlas of the oviduct among high- and low-fecundity goats in the follicular phase (FH vs. FL), luteal phase (LH vs. LL), and estrus cycle (FH vs. LH; FL vs. LL) to unravel their potential regulatory mechanisms in improving kidding number. We generated RNA sequencing data, and identified 4,078 circRNAs from twenty sampled Yunshang black goats. Many of these circRNAs are exon-derived and differentially expressed between each comparison group. Subsequently, eight differentially expressed (DE) circRNAs were validated by RT‒qPCR, which was consistent with the RNA-seq data. GO and KEGG enrichment analyses suggested that numerous host genes of DE circRNAs were involved in the hormone secretion, gamete production, fertilization, and embryo development processes. The competing endogenous RNA (ceRNA) interaction network analysis revealed that 2,673 circRNA–miRNA–mRNA axes (including 15 DE circRNAs, 14 miRNAs, and 1,699 mRNAs) were formed, and several target genes derived from the ceRNA network were associated with oviduct functions and reproduction, including SMAD1, BMPR1B, IGF1, REV1, and BMP2K. Furthermore, miR-15a-5p, miR-181b-5p, miR-23b-5p, miR-204-3p, and miR-145-5p might play important roles in reproduction. Finally, a novel circRNA, circIQCG, was identified as potentially involved in embryo development. Overall, our study provides a resource of circRNAs to understand the oviductal function and its connection to prolificacy trait of goats in the differentiation estrus cycle.

Genetic Diversity and Selection Signatures in Jianchang Black Goats Revealed by Whole-Genome Sequencing Data

Understanding the genetic composition of indigenous goats is essential to promote the scientific conservation and sustainable utilization of these breeds. The Jianchang Black (JC) goat, a Chinese native breed, is solid black and exhibits crude feed tolerance, but is characterized by a low growth rate and small body size. Based on the whole-genome sequencing data for 30 JC, 41 Jintang Black (JT), and 40 Yunshang Black (YS) goats, and 21 Bezoar ibexes, here, we investigated the genetic composition of JC goats by conducting analyses of the population structure, runs of homozygosity (ROH), genomic inbreeding, and selection signature. Our results revealed that JT and YS showed a close genetic relationship with a non-negligible amount of gene flows but were genetically distant from JC, apart from Bezoars. An average of 2039 ROHs were present in the autosomal genome per individual. The ROH-based inbreeding estimates in JC goats generally showed moderate values ranging from 0.134 to 0.264, mainly due to rapid declines in the effective population size during recent generations. The annotated genes (e.g., IL2, IL7, and KIT) overlapping with ROH islands were significantly enriched in immune-related biological processes. Further, we found 61 genes (e.g., STIM1, MYO9A, and KHDRBS2) under positive selection in JC goats via three complementary approaches, which may underly genetic adaptations to local environmental conditions. Our findings provided references for the conservation and sustainable utilization of JC goats.