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

Comprehensive circular RNA profiling in various sheep tissues

Despite the scientific relevance of circular RNAs (circRNAs), the study of these RNAs in non-model organisms, especially in sheep, is still in its infancy. On the other hand, while some studies have focused on sheep circRNA identification in a limited number of tissues, there is a lack of comprehensive analysis that profile circRNA expression patterns across the tissues not yet investigated. In this study, 61 public RNA sequencing datasets from 12 different tissues were uniformly analyzed to identify circRNAs, profile their expression and investigate their various characteristics. We reported for the first time a circRNA expression landscape with functional annotation in sheep tissues not yet investigated including hippocampus, BonMarrowMacrophage, left-ventricle, thymus, ileum, reticulum and 23-day-embryo. A stringent computational pipeline was employed and 8919 exon-derived circRNAs with high confidence were identified, including 88 novel circRNAs. Tissue-specificity analysis revealed that 3059 circRNAs were tissue-specific, which were also more specific to the tissues than linear RNAs. The highest number of tissue-specific circRNAs was found in kidney, hippocampus and thymus, respectively. Co-expression analysis revealed that expression of circRNAs may not be affected by their host genes. While most of the host genes produced more than one isoform, only one isoform had dominant expression across the tissues. The host genes of the tissue-specific circRNAs were significantly enriched in biological/pathways terms linked to the important functions of their corresponding tissues, suggesting potential roles of circRNAs in modulating physiological activity of those tissues. Interestingly, functional terms related to the regulation and various signaling pathways were significantly enriched in all tissues, suggesting some common regulatory mechanisms of circRNAs to modulate the physiological functions of tissues. Finding of the present study provide a valuable resource for depicting the complexity of circRNAs expression across tissues of sheep, which can be useful for the field of sheep genomic and veterinary research.

Hybridization promotes growth performance by altering rumen microbiota and metabolites in sheep

Hybridization can substantially improve growth performance. This study used metagenomics and metabolome sequencing to examine whether the rumen microbiota and its metabolites contributed to this phenomenon. We selected 48 approximately 3 month-old male ♂Hu × ♀Hu (HH, n = 16), ♂Poll Dorset × ♀Hu (DH, n = 16), and ♂Southdown × ♀Hu (SH, n = 16) lambs having similar body weight. The sheep were fed individually under the same nutritional and management conditions for 95 days. After completion of the trial, seven sheep close to the average weight per group were slaughtered to collect rumen tissue and content samples to measure rumen epithelial parameters, fermentation patterns, microbiota, and metabolite profiles. The final body weight (FBW), average daily gain (ADG), and dry matter intake (DMI) values in the DH and SH groups were significantly higher and the feed-to-gain ratio (F/G) significantly lower than the value in the HH group; additionally, the papilla height in the DH group was higher than that in the HH group. Acetate, propionate, and total volatile fatty acid (VFA) concentrations in the DH group were higher than those in the HH group, whereas NH3-N concentration decreased in the DH and SH groups. Metagenomic analysis revealed that several Prevotella and Fibrobacter species were significantly more abundant in the DH group, contributing to an increased ability to degrade dietary cellulose and enrich their functions in enzymes involved in carbohydrate breakdown. Bacteroidaceae bacterium was higher in the SH group, indicating a greater ability to digest dietary fiber. Metabolomic analysis revealed that the concentrations of rumen metabolites (mainly lysophosphatidylethanolamines [LPEs]) were higher in the DH group, and microbiome-related metabolite analysis indicated that Treponema bryantii and Fibrobacter succinogenes were positively correlated with the LPEs. Moreover, we found methionine sulfoxide and N-methyl-4-aminobutyric acid were characteristic metabolites in the DH and SH groups, respectively, and are related to oxidative stress, indicating that the environmental adaptability of crossbred sheep needs to be further improved. These findings substantially deepen the general understanding of how hybridization promotes growth performance from the perspective of rumen microbiota, this is vital for the cultivation of new species and the formulation of precision nutrition strategies for sheep.

Transcriptomic Analysis Reveals Differentially Expressed Circular RNAs Associated with Fecundity in the Sheep Hypothalamus with Different FecB Genotypes

Circular RNAs (circRNAs) are a specific type of noncoding RNA, and some have defined roles in cellular and biological processes. However, little is known about the role of circRNAs in follicular development in sheep with FecB (fecundity Booroola) mutations. Here, the expression profiles of circRNAs were investigated using RNA sequencing (RNA-seq) in the follicular phase (F) and the luteal phase (L) of FecB mutant homozygous (BB) and wild-type (WW) Small Tail Han sheep. A total of 38,979 circRNAs were identified, and 314, 343, 336, and 296 of them were differentially expressed (DE) between BB_F and BB_L, WW_F and WW_L, BB_F and WW_F, and BB_L and WW_L, respectively. The length, type, and chromosome distribution of the circRNAs and the expression characteristic between the circRNAs and their host genes in the sheep hypothalamus were ascertained. Enrichment analysis showed that the host genes of DE circRNAs in the follicular and luteal phases were annotated to MAPK, gap junctions, progesterone-mediated oocyte maturation, oocyte meiosis, and other hormone-related signaling pathways, and the different FecB genotypes were annotated to the gap junctions, circadian entrainment, MAPK, and other hormone-related signaling pathways. The competing endogenous RNA network prediction revealed that the 129 target miRNAs might be bound to 336 DE circRNAs. oar_circ_0000523 and oar_circ_0028984, which were specifically expressed during the follicular phase in the BB genotype sheep, probably acted as miRNA sponges involved in the regulation of LH synthesis and secretion. This study reveals the expression profiles and characterization of circRNAs at two phases of follicular development considering different FecB genotypes, thereby providing an improved understanding of the roles of circRNAs in the sheep hypothalamus and their involvement in follicular development and ovulation.

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.

Transcriptome analysis reveals pituitary lncRNA, circRNA and mRNA affecting fertility in high- and low-yielding goats

The pituitary gland serves as the central endocrine regulator of growth, reproduction, and metabolism and plays a crucial role in the reproductive process of female animals. Transcriptome analysis was conducted using pituitary gland samples from Leizhou goats with varying levels of fecundity to investigate the effects of long noncoding RNA (lncRNA), circular RNA (circRNA), and mRNA regulation on pituitary hormone secretion and its association with goat fecundity. The analysis aimed to identify lncRNAs, circRNAs, and mRNAs that influence the fertility of Leizhou goats. GO and KEGG enrichment analyses were performed on differentially expressed lncRNAs, circRNAs, and mRNAs and revealed considerable enrichment in pathways, such as regulation of hormone secretion, germ cell development, and gonadotropin-releasing hormone secretion. The pituitary lncRNAs (ENSCHIT00000010293, ENSCHIT00000010304, ENSCHIT00000010306, ENSCHIT00000010290, ENSCHIT00000010298, ENSCHIT00000006769, ENSCHIT00000006767, ENSCHIT00000006921, and ENSCHIT00000001330) and circRNAs (chicirc_029285, chicirc_026618, chicirc_129655, chicirc_018248, chicirc_122554, chicirc_087101, and chicirc_078945) identified as differentially expressed regulated hormone secretion in the pituitary through their respective host genes. Additionally, differential mRNAs (GABBR2, SYCP1, HNF4A, CBLN1, and CDKN1A) influenced goat fecundity by affecting hormone secretion in the pituitary gland. These findings contribute to the understanding of the molecular mechanisms underlying pituitary regulation of fecundity in Leizhou goats.

Comparative Transcriptomics Identify Key Pituitary Circular RNAs That Participate in Sheep (Ovis aries) Reproduction

CircRNAs have been found to play key roles in many biological processes and have diverse biological functions. There have been studies on circRNAs in sheep pituitary, and some important circRNAs have been found. But there are still few studies on circRNAs in sheep pituitary with different fecundity. In this study, we obtained the circRNAs expression profiles in the pituitary of FecB ++ genotype Small Tail Han sheep with different fecundity and estrous phases. A total of 34,878 circRNAs were identified in 12 pituitary samples, 300 differentially expressed circRNAs (DE circRNAs) (down: 104; up: 196) were identified in polytocous sheep in the follicular phase (PF) and monotocous sheep in the follicular phase (MF) (PF vs. MF), and 347 DE circRNAs (down: 162; up: 185) were identified in polytocous sheep in the luteal phase (PL) and monotocous sheep in the luteal phase (ML) (PL vs. ML). Cortisol synthesis and secretion pathway (follicular phase) and estrogen signaling pathway (luteal phase) were obtained by functional enrichment analysis of circRNAs source genes. Competing endogenous RNA (ceRNA) network analysis of key DE circRNAs revealed that oar-circ-0022776 (source gene ITPR2, follicular phase) targeted oar-miR-432, oar-circ-0009003 (source gene ITPR1, luteal phase) and oar-circ-0003113 (source gene PLCB1, luteal phase) targeted oar-miR-370-3p. We also explored the coding ability of DE circRNAs. In conclusion, our study shows that changes in the pituitary circRNAs may be related to the response of the pituitary to steroid hormones and regulate the reproductive process of sheep by affecting the pituitary function. Results of this study provide some new information for understanding the functions of circRNAs and the fecundity of FecB ++ genotype sheep.

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.

circITGB1 Regulates Adipocyte Proliferation and Differentiation via the miR-23a/ARRB1 Pathway

Adipose tissues represent an important energy storage organ in animals and are the largest endocrine organ. It plays an important regulatory role in the pathogenesis of insulin resistance, cardiovascular disease, and metabolic syndrome. Adipose development is a complex biological process involving multiple key genes, signaling pathways, and non-coding RNAs, including microRNAs and circular RNAs. In this study, we characterized circITGB1 and named its host gene ITGB1, which is differentially expressed in sheep of different months based on sequencing data. We collated and analyzed the sequencing data to select miRNA-23a with strong binding to ARRB1. We found that miRNA-23a regulates the development and differentiation of sheep adipocytes by targeting ARRB1. As a competing endogenous RNA, circITGB1 overexpression effectively alleviated the inhibitory effect of miR-23a on ARRB1. Conclusively, we provide evidence that circITGB1 regulates the proliferation and differentiation of sheep adipocytes via the miR-23a/ARRB1 pathway. This study provides a scientific basis for further studies on adipose tissue development at the circRNA level.

Comprehensive Analysis of Differentially Expressed CircRNAs in the Ovaries of Low- and High-Fertility Sheep

CircRNAs are essential in regulating follicle growth and development and the female reproductive system at multiple levels. However, the molecular mechanism by which circRNAs regulate reproduction in sheep is unclear and requires further exploration. In this study, RNA sequencing was performed to reveal the circRNA expression profiles in the ovaries of Cele black sheep and Hetian sheep during estrus. Analysis of the number of circRNAs in their host genes revealed that 5031 genes could produce 20,835 circRNAs. Among the differentially expressed circRNAs (DEcircRNA), 75 were upregulated, and 105 were downregulated. Functional enrichment analysis showed that the host genes of DEcircRNA were involved in several pathways, including the MAPK and Hippo signaling pathway. In addition, we constructed a subnetwork of competitive endogenous RNA (ceRNA) containing 4 mRNAs, 4 microRNAs (miRNAs), and 10 circRNAs, potentially related to follicle development. Functional circRNAs (e.g., novel_circ_0003851, novel_circ_0015526, novel_circ_0008117) were found to act as ceRNAs for follicle growth and development-related mRNAs (CUEDC1, KPNB1, ZFPM2) by sponging functional miRNAs (miR-29a, miR-29b, miR-17-5p). Finally, through an RNA pull-down assay, oar-miR-125b was selected and confirmed as the target miRNA of novel-circ-0041512. We analyzed the overall expression of circRNAs in sheep ovaries. Further, we explored the potential mechanisms underlying the circRNA functions, providing a theoretical basis for the genetic progress of reproductive traits in sheep.

Identification of hub genes associated with follicle development in multiple births sheep by WGCNA

Sheep exhibit a distinct estrous cycle that includes four different phases: proestrus, estrus, late estrus, and luteal phase. As the estrous cycle repeats, follicular development regularly alternates. We thus investigated ovarian transcriptome data from each of the four phases using weighted gene co-expression network analysis (WGCNA) to identify modules, pathways, and genes essential to follicle growth and development. We clustered mRNA and long non-coding RNA (lncRNA) into different modules by WGCNA, and calculated correlation coefficients between genes and Stages of the estrous cycle. Co-expression of the black module (cor = 0.81, P<0.001) and the yellow module (cor = 0.61, P<0.04) was found to be critical for follicle growth and development. A total of 2066 genes comprising the black and yellow modules was used for functional enrichment. The results reveal that these genes are mainly enriched in Cell cycle, PI3K-Akt signaling pathway, Oocyte meiosis, Apoptosis, and other important signaling pathways. We also identified seven hub genes (BUB1B, MAD2L1, ASPM, HSD3B1, WDHD1, CENPA, and MXI1) that may play a role in follicle development. Our study may provide several important new markers allowing in depth exploration of the genetic basis for multiparous reproduction in sheep.

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.

Photoperiods induced the circRNA differential expression in the thyroid gland of OVX+E2 ewes

Circular RNAs (circRNAs) are non-coding RNAs newly identified and play important roles in RNA regulation. However, little is known regarding photoperiods induced circRNAs in the thyroid gland. In this study, we performed a comprehensive analysis of circRNA profiles in the thyroid gland of OVX+E₂ ewes at different photoperiods by whole transcriptome sequencing. A total of 37,470 novel circRNAs were detected in different photoperiods (42 days of short photoperiod treatment, SP42; 42 days of long photoperiod treatment, LP42; SP42 transfer to LP42, SPLP42), with a total of 817 circRNAs for SP42-LP42 (down: 132; up: 114), LP42-SPLP42 (down: 136; up: 112) and SP42-SPLP42 (down: 182; up: 141) having differentially expressed. Functional enrichment annotation analysis of DE-circRNAs for GO and KEGG by R package, features that influence photoperiod response in Sunite ewes through the Inositol phosphate metabolism, cGMP-PKG signaling pathway, Calcium signaling pathway, MAPK signaling pathway, and Oocyte meiosis. In addition, competitive endogenous RNA (ceRNA) network analysis revealed target binding sites for identified miRNAs in DE-cirRNAs such as oar-miR-10b, oar-miR-200c, oar-miR-21, oar-miR-370-3p, oar-miR-377-3p, oar-miR-181a, oar-miR-432, and oar-miR-495-3p. These results of this study will provide some new information for understanding circRNA function as well as the changes in the sheep thyroid gland under different photoperiods.

Analysis of Expression Profiles of CircRNA and MiRNA in Oviduct during the Follicular and Luteal Phases of Sheep with Two Fecundity (FecB Gene) Genotypes

CircRNA and miRNA, as classes of non-coding RNA, have been found to play pivotal roles in sheep reproduction. There are many reports of circRNA and miRNA in the ovary and uterus, but few in the oviduct. In this study, RNA-Seq was performed to analyze the expression profile of circRNA and miRNA in the oviduct during the follicular phase and luteal phase of sheep with FecB^BB and FecB⁺⁺ genotypes. The results showed that a total of 3223 circRNAs and 148 miRNAs were identified. A total of 15 DE circRNAs and 40 DE miRNAs were found in the comparison between the follicular phase and luteal phase, and 1 DE circRNA and 18 DE miRNAs were found in the comparison between the FecB^BB genotype and FecB⁺⁺ genotype. GO and KEGG analyses showed that the host genes of DE circRNAs were mainly enriched in the Rap1 signaling pathway, PI3K-Akt signaling pathway and neuroactive ligand-receptor interactions. Novel_circ_0004065, novel_circ_0005109, novel_circ_0012086, novel_circ_0014274 and novel_circ_0001794 were found to be possibly involved in the oviductal reproduction process. GO and KEGG analyses showed that the target genes of DE miRNAs were mainly enriched in insulin secretion, the cAMP signaling pathway, the cGMP-PKG signaling pathway, the Rap1 signaling pathway and the TGF-β signaling pathway, and the target genes LPAR1, LPAR2, FGF18, TACR3, BMP6, SMAD4, INHBB, SKP1 and TGFBR2 were found to be associated with the reproductive process. Miranda software was used to identify 27 miRNAs that may bind to 13 DE circRNAs, including miR-22-3p (target to novel_circ_0004065), miR-127, miR-136 (target to novel_circ_0000417), miR-27a (target to novel_circ_0014274) and oar-miR-181a (target to novel_circ_ 0017815). The results of this study will help to elucidate the regulatory mechanisms of circRNAs and miRNAs in sheep reproduction. Our study, although not establishing direct causal relationships of the circRNA and miRNA changes, enriches the sheep circRNA and miRNA database and provides a basis for further studies on sheep reproduction.