Pituitary pars tuberalis (PT) plays an important role as the transmission center in the seasonal reproduction of animals. It helps convert external photoperiod signals into intrinsic seasonal reproduction signals. In sheep PT, specific expression patterns of several genes (including short photoperiod-induced gene CHGA and long photoperiod genes EYA3 and TSHβ) under different photoperiods are crucial characteristics during this signal transduction. Recent studies have revealed the role of epigenetics in regulating the expression of seasonal reproductive key genes. Therefore, we explored whether microRNAs and LncRNAs regulated the expressions of the above key genes. Firstly, the expression of miR-25 and CHGA showed a significant negative correlation in sheep PT. Results of the dual luciferase reporter assay and miR-25 overexpression indicated that miR-25 could inhibit the expression of CHGA by specifically binding to its 3'UTR region in pituitary cells. Then, expression negative correlation and dual luciferase reporter analyses were used to screen and identify the candidate LncRNA (Lnc107153) targeted by miR-25. Finally, the results of fluorescence in situ hybridization and Lnc107153 overexpression suggested that Lnc107153 and miR-25 were involved in the epigenetic regulation of CHGA expression. However, the expressions of EYA3 and TSHβ were not regulated by miRNAs. These results will provide new insights into the epigenetic regulatory network of key genes in sheep seasonal reproduction.

Most sheep are seasonal estrus, and they breed in autumn when the days get shorter. Seasonal estrus is an important factor that affects the productivity and fertility of sheep. The key point to solve this problem is to explore the regulation mechanism of estrus in sheep. Therefore, in this study, transcriptomic sequencing technology was used to identify differentially expressed mRNAs in the hypothalamus, pituitary and ovary of Small Tail Han sheep (year-round estrus) and tan sheep (seasonal estrus) among luteal, proestrus and estrus stages. There were 256,923,304,156 mRNAs being identified in the hypothalamus, pituitary and ovary, respectively. Functional analysis showed that the photosensor, leucine and isoleucine biosynthesis pathways were enriched significantly. It is speculated that photoperiod may initiate estrus by stimulating the corresponding pathways in hypothalamus. ODC1, PRLH, CRYBB2, SMAD5, OPN1SW, TPH1 are believed to be key genes involved in the estrogen process. In conclusion, this study expanded the database of indigenous sheep breeds, and also provided new candidate genes for future genetic and molecular studies on the seasonal estrus trait in sheep.

Seasonal estrus is a critical limiting factor of animal fecundity, and it involves changes in both ovarian biology and hormone secretion in different seasons. Previous studies indicate that two classes of small RNAs (miRNAs and piRNAs) play important regulatory roles in ovarian biology. To understand the roles of small RNA-mediated post-transcriptional regulation in ovine seasonal estrus, the variation in expression patterns of ovarian small RNAs during anestrus and the breeding season were analyzed using Solexa sequencing technology. In addition, reproductive hormone levels were determined during ovine anestrus and the breeding season.

Most ewes in China are seasonally polyestrous with normal ovulatory cycles, which is controlled by photoperiod (length of the daily light phase). These ewes are estrous in the short-day season and anestrus in the long-day season and cannot mate during anestrus. Thus seasonal breeding limits both diversification and intensification of production. If sheep can estrus all round year, it can be mated twice per year, which can greatly improve the economic benefits. To change seasonal estrus at the genetic level and cultivating new sheep breeds, it is important to understand the molecular mechanisms of seasonal breeding trait in sheep. The BCL2-associated agonist of cell death (BAD) gene being a regulator of cellular apoptosis was identified by our previous RNA-Seq, which is associated with follicular development in mammalian ovaries. However, the mechanism how BAD can regulate estrus in sheep was poorly understood. In this study, we characterized ovine BAD, including full-length mRNA cloning and protein sequence prediction, as well as BAD expression profile in Small-tailed Han (STH) sheep. The highest expression levels of BAD were observed in sheep hypothalamus, lung, and pituitary, while the lowest expression was in liver. Functional analysis of BAD was performed in primary granulosa cells of sheep. The concentration of P4 was significantly increased after RNAi interference of BAD, while P4 level was shown to be opposite after BAD overexpression in vitro. It has been found that BAD can reduce progesterone levels by promoting ovarian GC apoptosis, which might be involved in regulating the estrus cycle in sheep.

Small-tailed Han sheep, with different FecB genotypes, manifest distinct ovulation rates and fecundities, which are due to differences in reproductive hormones secreted by the hypothalamic-pituitary-ovarian axis. Nevertheless, the function of the hypothalamus against a FecB mutant background on increasing ovulation rate is rarely reported. Therefore, we determined the expression profiles of hypothalamus tissue collected from six wild-type (WW) and six FecB mutant homozygous (BB) ewes at the follicular and luteal phases by whole-transcriptome sequencing. We identified 53 differentially expressed mRNAs (DEGs) and 40 differentially expressed long non-coding RNAs (DELs) between the two estrus states. Functional annotation analysis revealed that one of the DEGs, PRL, was particularly enriched in the hypothalamic function, hormone-related, and reproductive pathways. The lncRNA-target gene interaction networks and KEGG analysis in combination suggest that the lncRNAs LINC-676 and WNT3-AS cis-acting on DRD2 and WNT9B in different phases may induce gonadotropin-releasing hormone (GnRH) secretion. Furthermore, there were differences of regulatory elements and WNT gene family members involved in the follicular-luteal transition in the reproductive process between wild-type (WNT7A) and FecB mutant sheep (WNT9B). We combined the DEG and DEL data sets screened from different estrus states and genotypes. The overlap of these two sets was identified to select the mRNAs and lncRNAs that have major effects on ovulation. Among the overlapping molecules, seven DEGs and four DELs were involved in the follicular-luteal transition regulated by FecB mutation. Functional annotation analysis showed that two DEGs (FKBP5 and KITLG) were enriched in melanogenesis, oxytocin, and GnRH secretion. LINC-219386 and IGF2-AS were highly expressed in the BB ewes compared with WW ewes, modulating their target genes (DMXL2 and IGF2) to produce more GnRH during follicular development, which explains why mutated ewes produced more mature follicles. These results from expression profiling of the hypothalamus with the FecB mutation at different estrus states provide new insights into how the hypothalamus regulates ovulation under the effect of the FecB mutation.

The expression characteristics of the hypothalamic-pituitary-gonadal (HPG) axis-related candidate genes, DIO2, EYA3, KISS1 and GPR54, were analyzed in year-round estrous rams (small-tail Han sheep, STH) and seasonally estrous rams (Sunite sheep, SNT) using qPCR. The results were as follows: DIO2 was mainly expressed in pituitary, and KISS1 was specifically expressed in hypothalamus in the two groups. However, EYA3 and GPR54 were widely expressed in the cerebrum, cerebellum, hypothalamus, pituitary, testis, epididymis, vas deferens and adrenal gland tissues in both breeds, with significant differences in the cerebellum, hypothalamus, pituitary, testis and vas deferens tissues. We speculated that DIO2 and KISS1 may have positive roles in different regions in ram year-round estrus. Moreover, the expression patterns of EYA3 and GPR54 suggested that they may regulate the estrous mode of ram via testis and vas deferens. This is the first study to systematically analyze the expression patterns of HPG axis-related genes in rams.

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.

The circadian rhythm is a biological rhythm that is closely related to the rhythmic expression of a series of clock genes. Results from several studies have indicated that clock genes are associated with the estrous cycle in female animals. Until now, the relationship between estrus cycle transition and clock gene expression in reproductive-axis-related tissues has remained unknown in Small-tailed Han (STH) sheep. This study was conducted to analyze the expression patterns of six canonical clock genes (Clock, BMAL1, Per1, Per2, Cry1, and Cry2) in the follicle phase and luteal phase of STH sheep. We found that all six genes were expressed in the brain, cerebellum, hypothalamus, pituitary, ovary, uterus, and oviduct in follicle and luteal phases. The results indicated that Clock expression was significantly higher in the cerebellum, hypothalamus, and uterus of the luteal phase than that of the follicle phase, whereas BMAL1 expression was significantly higher in the hypothalamus of the luteal phase than that of the follicle phase. Per1 expression was significantly higher in the brain, cerebellum, hypothalamus, and pituitary of the luteal phase than that of the follicle phase, and Per2 expression was significantly higher in the hypothalamus, pituitary, and uterus of the luteal phase than that of the follicle phase. Cry1 expression was significantly higher in the brain, cerebellum, and hypothalamus of the luteal phase than that of the follicle phase, whereas Cry2 expression was significantly higher in the pituitary of the luteal phase than that of the follicle phase. The clock gene expression in all tissues was different between follicle and luteal phases, but all clock gene mRNA levels were found to exhibit higher expression among seven tissues in the luteal phase. Our results suggest that estrous cycles may be associated with clock gene expression in the STH sheep. This is the first study to systematically analyze the expression patterns of clock genes of different estrous cycle in ewes, which could form a basis for further studies to develop the relationship between clock genes and the estrous cycle.