Superior kidding rate is an important economic trait in production of meat goat, and ovulation rate is the precondition of kidding rate. MicroRNAs (miRNAs) play critical roles in almost all ovarian biological processes, including folliculogenesis, follicle development, follicle atresia, luteal development and regression. To find out the different ovarian activity and follicle recruitment with miRNA-mediated posttranscriptional regulation, the small RNAs expressed pattern in the ovarian tissues of multiple and uniparous Anhui White goats during follicular phase was analyzed using Solexa sequencing data. 1008 miRNAs co-expressed, 309 and 433 miRNAs specifically expressed in the ovaries of multiple and uniparous goats during follicular phase were identified. The 10 most highly expressed miRNAs in the multiple library were also the highest expressed in the uniparous library, and there were no significantly different between each other. The highest specific expressed miRNA in the multiple library was miR-29c, and the one in the uniparous library was miR-6406. 35 novel miRNAs were predicted in total. GO annotation and KEGG Pathway analyses were implemented on target genes of all miRNA in two libraries. RT-PCR was applied to detect the expression level of 5 randomly selected miRNAs in multiple and uniparous hircine ovaries, and the results were consistent with the Solexa sequencing data. In the present study, the different expression of miRNAs in the ovaries of multiple and uniparous goats during follicular phase were characterized and investigated using deep sequencing technology. The result will help to further understand the role of miRNAs in kidding rate regulation and also may help to identify miRNAs which could be potentially used to increase hircine ovulation rate and kidding rate in the future.

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.

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.

Maintenance of appropriate levels of fecundity is critical for efficient sheep production. Opportunities to increase sheep litter size include identifying single gene mutations with major effects on ovulation rate and litter size. Whole-genome sequencing (WGS) data of 89 Chinese domestic sheep from nine different geographical locations and ten Australian sheep were analyzed to detect new polymorphisms affecting litter size. Comparative genomic analysis of sheep with contrasting litter size detected a novel set of candidate genes. Two SNPs, g.71874104G>A and g.71833755T>C, were genotyped in 760 Small Tail Han sheep and analyzed for association with litter size. The two SNPs were significantly associated with litter size, being in strong linkage disequilibrium in the region 71.80–71.87 Mb. This haplotype block contains one gene that may affect litter size, Histone Cell Cycle Regulator (HIRA). HIRA mRNA levels in sheep with different lambing ability were significantly higher in ovaries of Small Tail Han sheep (high fecundity) than in Sunite sheep (low fecundity). Moreover, the expression levels of HIRA in eight tissues of uniparous Small Tail Han sheep were significantly higher than in multiparous Small Tail Han sheep (p < 0.05). HIRA SNPs significantly affect litter size in sheep and are useful as genetic markers for litter size.

Bamei mutton sheep is a Chinese domestic sheep breed developed by crossing German Mutton Merino sheep and indigenous Mongolian sheep for meat production. Here, we focused on detecting candidate genes associated with the increasing of the litter size in this breeds under recent artificial selection to improve the efficiency of mutton production. We selected five high- and five low-fecundity Bamei mutton sheep for whole-genome resequencing to identify candidate genes for sheep prolificacy. We used the FST and XP-EHH statistical approach to detect the selective sweeps between these two groups. Combining the two selective sweep methods, the reproduction-related genes JUN, ITPR3, PLCB2, HERC5, and KDM4B were detected. JUN, ITPR3, and PLCB2 play vital roles in GnRH (gonadotropin-releasing hormone), oxytocin, and estrogen signaling pathway. Moreover, KDM4B, which had the highest FST value, exhibits demethylase activity. It can affect reproduction by binding the promoters of estrogen-regulated genes, such as FOXA1 (forkhead box A1) and ESR1 (estrogen receptor 1). Notably, one nonsynonymous mutation (p.S936A) specific to the high-prolificacy group was identified at the TUDOR domain of KDM4B. These observations provide a new opportunity to research the genetic variation influencing fecundity traits within a population evolving under artificial selection. The identified genomic regions that are responsible for litter size can in turn be used for further selection.

Alpha 1-antitrypsin (α1-AT) may affect the susceptibility of mastitis in dairy cattle for its possible role in the protection of lactoferrin from proteolytic degradation in the mammary. Milk somatic cell score (SCS) is a logarithmic transformation of the milk somatic cell count widely used as an index to evaluate mastitis. To study the relationships of α1-AT gene and SCS in Chinese Holstein cows, methods of PCR-SSCP, DNA sequencing, PCR-RFLP, and CRS-PCR technologies were used to detect single nucleotide polymorphisms sites in α1-AT gene. Two polymorphic sites at G5503A and G5746C of α1-AT gene were found. AA (0.3633), AB (0.4644) and BB (0.1723) genotypes were detected at G5503A site, CC (0.3483), CD (0.4906) and DD (0.1611) genotypes were found at G5746C in Chinese Holstein cows. Least squares mean of SCS for individuals with BB genotype was significantly lower than that with AA and AB genotype (p < 0.01), and that with AB genotype was significantly lower than that with AA (p < 0.05). There was no significant difference among individuals with CC, CD and DD genotypes (p > 0.05). Least squares mean of SCS for individuals with BBDD genotype combination were significantly lower than those with AACC and AACD (p < 0.05). Statistical analysis indicated that B allele and BBDD genotype combination of α1-AT can improve mastitis resistance in dairy cattle.