Detection of selection footprints provides insight into the evolution process and the underlying mechanisms controlling the phenotypic diversity of traits that have been exposed to selection. Selection focused on certain characters, mapping certain genomic regions often shows a loss of genetic diversity with an increased level of homozygosity. Therefore, the runs of homozygosity (ROHs), homozygosity by descent (HBD), and effective population size (N_e) are effective tools for exploring the genetic diversity, understanding the demographic history, foretelling the signature of directional selection, and improving the breeding strategies to use and conserve genetic resources. We characterized the ROH, HBD, N_e, and signature of selection of six Chinese goat populations using single nucleotide polymorphism (SNP) 50K Illumina beadchips. Our results show an inverse relationship between the length and frequency of ROH. A long ROH length, higher level of inbreeding, long HBD segment, and smaller N_e in Guangfeng (GF) goats suggested intensive selection pressure and recent inbreeding in this breed. We identified six reproduction-related genes within the genomic regions with a high ROH frequency, of which two genes overlapped with a putative selection signature. The estimated pair-wise genetic differentiation (F_ST) among the populations is 9.60% and the inter- and intra-population molecular variations are 9.68% and 89.6%, respectively, indicating low to moderate genetic differentiation. Our selection signatures analysis revealed 54 loci harboring 86 putative candidate genes, with a strong signature of selection. Further analysis showed that several candidate genes, including MARF1, SYCP2, TMEM200C, SF1, ADCY1, and BMP5, are involved in goat fecundity. We identified 11 candidate genes by using cross-population extended haplotype homozygosity (XP-EHH) estimates, of which MARF1 and SF1 are under strong positive selection, as they are differentiated in high and low reproduction groups according to the three approaches used. Gene ontology enrichment analysis revealed that different biological pathways could be involved in the variation of fecundity in female goats. This study provides a new insight into the ROHs patterns for maintenance of within breed diversity and suggests a role of positive selection for genetic variation influencing fecundity in Chinese goat.

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. A total of 483 miRNAs (including 97 known, 369 conserved and 17 predicated novel miRNAs), which belong to 183 different miRNA families, were identified in ovaries of Tan sheep and Small Tail Han (STH) sheep. Compared with the three stages of the breeding season, 25 shared significantly differentially expressed (including 19 up- and six down-regulated) miRNAs were identified in ovine anestrus. KEGG Pathway analysis revealed that the target genes for some of the differentially expressed miRNAs were involved in reproductive hormone related pathways (e.g. steroid biosynthesis, androgen and estrogen metabolism and GnRH signaling pathway) as well as follicular/luteal development related pathways. Moreover, the expression of the differentially expressed miRNAs and most of their target genes were negatively correlated in the above pathways. Furthermore, the levels of estrogen, progesterone and LH in ovine anestrus were significantly lower than those in the breeding season. Combining the results of pathway enrichment analysis, expression of target genes and hormone measurement, we suggest that these differentially expressed miRNAs in anestrus might participate in attenuation of ovarian activity by regulating the above pathways. Besides miRNAs, a large and unexpectedly diverse set of piRNAs were also identified. The miRNA profiles of ovine ovaries in anestrus were presented for the first time. The identification and characterization of miRNAs that are differentially expressed between ovine anestrus and the breeding season will help understanding of the role of miRNAs in the regulation of seasonal estrus, and provides candidates for determining miRNAs which could be potentially used to regulate ovine seasonal estrus.

Prolificacy is a crucial characteristic of livestock, particularly for species such as sheep that have many births. The objectives of this study were as follows: (1) to investigate the genetic diversity of the 13 new and 7 known variants in the BMPRIB, GDF9, BMP15, LEPR, and B4GALNT2 genes in Ujimqin (UM), the F1 population of Dorper × Ujimqin crossbred (DPU), the F1 population of Suffolk × Ujimqin crossbred (SFKU), Sonid sheep (SN), Tan sheep (Tan), Hu sheep (Hu), and Small-tailed Han sheep (STH) sheep breeds/populations; (2) to perform an association analysis of the above 20 variants with litter size in 325 UM, 304 DPU, and 66 SFKU sheep populations; (3) to compare the frequencies of the litter-size-related alleles of these 20 variants among 8 sheep breeds/populations (the above seven sheep breeds + Mongolia sheep breed). With the use of the Sequenom MassARRAY^®SNP assay technology, these 20 mutations were genotyped. The association analysis results showed that the c.746A>G (FecB) mutation in BMPR1B was significantly associated with the litter size of UM and DPU, the c.994A>G (FecG^A) in GDF9 was significantly associated with the litter size of SFKU, and the c.31_33CTTinsdel (B1) in BMP15 was significantly associated with the litter size of UM. Our findings might provide valuable genetic markers for expanding sheep litter sizes.

Tan sheep are greatly preferred by consumers in China because of their nutritional value and unique flavor. However, the meat quality of Tan sheep has decreased due to the change in feeding systems from grazing to indoor. Studies investigating the mechanisms for the decrease in meat quality are limited. A total of 28 Tan sheep were randomly allocated to two treatments, receiving a concentrated diet, or pasture. Flavor precursors and volatile compounds were analyzed with foodomics. E-nose and E-tongue analyses suggested that the aroma and taste profiles differed between the feeding systems. The grazing lambs had higher levels of linoleic acid and n-3 polyunsaturated fatty acids (p < 0.01). Metabolomics indicated that 25 hydrophilic metabolites active in glycolipid and amino acid metabolisms were changed by the feeding system. Among the 57 volatile compounds identified, methional, γ-butyrolactone, benzaldehyde, and ethyl acetate were at concentrations significantly higher in the grazing lambs than the indoor-fed lambs (p < 0.01). These results reveal key changes in flavor precursors and flavor profiles affected by the feeding system, which may provide an initial view of the reason for consumer preference for the grazing Tan sheep.

The effects of the different feeding systems, graze feeding (GSF), time-limited graze feeding (GF), and stall-feeding (SF)) on the fatty acid content, metabolites, and genes expression of the longissimus dorsi (LD) in Tan lambs were investigated in the present study. Thirty-nine 4-month-old male Tan lambs with similar body weight (24.91 ± 1.74 kg) were selected and divided into the three feeding systems (n = 13) randomly. Lambs were slaughtered after 83 days of the feeding trails, and LD muscle samples were collected for further analysis. The results indicated that different feeding systems have no significant effect on short-chain fatty acids in Tan lambs (p > 0.05). However, the total saturated fatty acids (∑SFA) and monounsaturated fatty acids (∑MUFA) in the GSF and GF groups were lower than those in the SF group (p < 0.001). The total polyunsaturated fatty acids (∑PUFA) in the GSF group were higher than those in the GF and SF groups (p < 0.001). Moreover, in the comparison of both GF vs GSF groups and SF vs GSF groups, metabolomic analysis showed that metabolites such as cis-(6,9,12)-linolenic acid, arachidic acid, acetylcarnitine, and L-carnitine with lower concentration were significantly enriched in the biosynthesis of unsaturated fatty acid pathway (p < 0.05), but metabolites such as phosphorylcholine, glycerophosphocholine, cytidine 5’-diphosphocholine, and glycerol-3-phosphate with higher concentrations were enriched in the glycerophospholipid metabolism pathway. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis of the results indicated that in the comparison of the GSF group with the SF group, differentially expressed genes (DEGs) such as LIPC, ERFE, FABP3, PLA2R1, LDLR, and SLC10A6, were enriched in the steroid biosynthesis and cholesterol metabolism pathways. In addition, differential metabolites and genes showed a significant correlation with the content of ∑SFA, ∑MUFA, and ∑PUFA in lamb meat (p < 0.05). These findings demonstrated that the feeding system was an important factor in regulating fatty acid content by affecting lipid-metabolism-related metabolites and gene expression in muscle, and graze-feeding system provided lamb meat with higher ∑PUFA content than time-limited-grazing and stall-feeding systems.