Litter size is one of the most important economic traits in sheep. GDF9 and BMPR1B are major genes affecting the litter size of sheep. In this study, the whole coding region of GDF9 was sequenced and all the SNPs (single nucleotide polymorphisms) were determined in Luzhong mutton ewes. The FecB mutation was genotyped using the Sequenom MassARRAY^®SNP assay technology. Then, the association analyses between polymorphic loci of GDF9 gene, FecB, and litter size were performed using a general linear model procedure. The results showed that eight SNPs were detected in GDF9 of Luzhong mutton sheep, including one novel mutation (g.41769606 T > G). The g.41768501A > G, g.41768485 G > A in GDF9 and FecB were significantly associated with litter size in Luzhong mutton ewes. The g.41768485 G > A is a missense mutation in the mature GDF9 protein region and is predicted to affect the tertiary structure of the protein. The results preliminarily demonstrated that GDF9 was a major gene affecting the fecundity of Luzhong mutton sheep and the two loci g.41768501A > G and g.41768485 G > A may be potential genetic markers for improving litter size.

The FecB gene has been discovered as an important gene in sheep for its high relationship with the ovulation rate, but its regulatory mechanism remains unknown. In the present study, liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS) techniques were adopted to detect the metabolic effects of FecB gene in follicular fluid (FF) and ovarian vein serum (OVS) in Small Tail Han (STH) sheep. ANOVA and random forest statistical methods were employed for the identification of important metabolic pathways and biomarkers. Changes in amino acid metabolism, redox environment, and energy metabolism were observed in FF from the three FecB genotype STH ewes. Principal component analysis (PCA) and hierarchical clustering analysis (HCA) showed that metabolic effects of FecB gene are more pronounced in FF than in OVS. Therefore, the difference of the metabolic profile in FF is also affected by the FecB genotypes. In Spearman correlation analysis, key metabolites (e.g., glucose 6-phosphate, glucose 1-phosphate, aspartate, asparagine, glutathione oxidized (GSSG), cysteine-glutathione disulfide, γ-glutamylglutamine, and 2-hydrosybutyrate) in ovine FF samples showed a significant correlation with the ovulation rate. Our findings will help to explain the metabolic mechanism of high prolificacy ewes and benefit fertility identification.

Photoperiod is one of the important factors leading to seasonal reproduction of sheep. However, the molecular mechanisms underlying the photoperiod regulation of seasonal reproduction remain poorly understood. In this study, we compared the expression profiles of mRNAs, lncRNAs, and circRNAs in uterine tissues from Sunite sheep during three different photoperiods, namely, the short photoperiod (SP), short transfer to long photoperiod (SLP), and long photoperiod (LP). The results showed that 298, 403, and 378 differentially expressed (DE) mRNAs, 171, 491, and 499 DE lncRNAs, and 124, 270, and 400 DE circRNAs were identified between SP and LP, between SP and SLP, and between LP and SLP, respectively. Furthermore, functional enrichment analysis showed that the differentially expressed RNAs were mainly involved in the GnRH signaling pathway, thyroid hormone synthesis, and thyroid hormone signaling pathway. In addition, co-expression networks of lncRNA–mRNA were constructed based on the correlation analysis between the differentially expressed RNAs. Our study provides new insights into the expression changes of RNAs in different photoperiods, which might contribute to understanding the molecular mechanisms of seasonal reproduction in sheep.

The litter size of domestic goats and sheep is an economically important trait that shows variation within breeds. Strenuous efforts have been made to understand the genetic mechanisms underlying prolificacy in goats and sheep. However, there has been a paucity of research on the genetic convergence of prolificacy between goats and sheep, which likely arose because of similar natural and artificial selection forces. Here, we performed comparative genomic and transcriptomic analyses to identify the genetic convergence of prolificacy between goats and sheep. By combining genomic and transcriptomic data for the first time, we identified this genetic convergence in (1) positively selected genes (CHST11 and SDCCAG8), (2) differentially expressed genes (SERPINA14, RSAD2, and PPIG at follicular phase, and IGF1, GPRIN3, LIPG, SLC7A11, and CHST15 at luteal phase), and (3) biological pathways (genomic level: osteoclast differentiation, ErbB signaling pathway, and relaxin signaling pathway; transcriptomic level: the regulation of viral genome replication at follicular phase, and protein kinase B signaling and antigen processing and presentation at luteal phase). These results indicated the potential physiological convergence and enhanced our understanding of the overlapping genetic makeup underlying litter size in goats and sheep.

Circular RNA (circRNA), as an emerging class of noncoding RNA, has been found to play key roles in many biological processes. However, its expression profile in the hypothalamus, a powerful organ initiating the reproductive process, has not yet been explored. Therefore, we used RNA sequencing to explore the expression of circRNAs in the hypothalamus of sheep with the FecB ++ genotype. We totally identified 41,863 circRNAs from sheep hypothalamus, in which 333 (162 were upregulated, while 171 were downregulated) were differentially expressed in polytocous sheep in the follicular phase versus monotocous sheep in the follicular phase (PF vs. MF), moreover, 340 circRNAs (163 were upregulated, while 177 were downregulated) were differentially expressed in polytocous sheep in the luteal phase versus monotocous sheep in the luteal sheep (PL vs. ML). We also identified several key circRNAs including oar_circ_0018794, oar_circ_0008291, oar_circ_0015119, oar_circ_0012801, oar_circ_0010234, and oar_circ_0013788 through functional enrichment analysis and oar_circ_0012110 through a competing endogenous RNA network, most of which may participate in reproduction by influencing gonadotropin-releasing hormone (GnRH) activities or affecting key gene expression, indirectly or directly. Our study explored the overall expression profile of circRNAs in sheep hypothalamus, which potentially provides an alternative insight into the mechanism of sheep prolificacy without the effects of FecB mutation.

Ovine litter size (LS) is an important trait showing variability within breeds. It remains largely unknown whether inbreeding depression on LS exists based on genomic homozygous regions, and whether the homozygous regions resulted from inbreeding are significantly associated with LS in sheep. We here reanalyze a set of single nucleotide polymorphism (SNP) chip of six breeds to characterize the patterns of runs of homozygosity (ROH), to evaluate inbreeding levels and inbreeding depressions on LS, and to identify candidate homozygous regions responsible for LS. Consequently, unique ROH patterns were observed among six sheep populations. Inbreeding depression on LS was only found in Hu sheep, where a significant reduction of 0.016, 0.02, and 0.02 per 1% elevated inbreeding F_ROH4–8, F_ROH>8 and the total inbreeding measure was observed, respectively. Nine significantly homozygous regions were found for LS in Hu sheep, where some promising genes for LS possibly via regulation of the development of oocytes (NGF, AKT1, and SYCP1), fertilization (SPAG17, MORC1, TDRD9, ZFYVE21, ADGRB3, and CKB), embryo implantation (PPP1R13B, INF2, and VANGL1) and development (DPPA2, DPPA4, CDCA4, CSDE1, and ADSSL1), and reproductive health (NRG3, BAG5, CKB, and XRCC3) were identified. These results from the present study would provide insights into the genetic management and complementary understandings of LS in sheep.