Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Wool fineness affects the quality of wool, and some studies have identified about forty candidate genes that affect sheep wool fineness, but these genes often reveal only a certain proportion of the variation in wool thickness. We further explore additional genes associated with the fineness of sheep wool. Whole-genome resequencing of eight sheep breeds was performed to reveal selection signals associated with wool fineness, including four coarse wool and four fine/semi-fine wool sheep breeds. Multiple methods to reveal selection signals (Fst and θπ Ratio and XP-EHH) were applied for sheep wool fineness traits. In total, 269 and 319 genes were annotated in the fine wool (F vs. C) group and the coarse wool (C vs. F) group, such as LGR4, PIK3CA, and SEMA3C and NFIB, OPHN1, and THADA. In F vs. C, 269 genes were enriched in 15 significant GO Terms (p < 0.05) and 38 significant KEGG Pathways (p < 0.05), such as protein localization to plasma membrane (GO: 0072659) and Inositol phosphate metabolism (oas 00562). In C vs. F, 319 genes were enriched in 21 GO Terms (p < 0.05) and 16 KEGG Pathways (p < 0.05), such as negative regulation of focal adhesion assembly (GO: 0051895) and Axon guidance (oas 04360). Our study has uncovered genomic information pertaining to significant traits in sheep and has identified valuable candidate genes. This will pave the way for subsequent investigations into related traits.

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, F_st & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

Hu sheep and Tibetan sheep are two commonly raised local sheep breeds in China, and they have different morphological characteristics, such as tail type and adaptability to extreme environments. A fat tail in sheep is the main adipose depot in sheep, whereas the liver is an important organ for fat metabolism, with the uptake, esterification, oxidation, and secretion of fatty acids (FAs). Meanwhile, adaptations to high-altitude and arid environments also affect liver metabolism. Therefore, in this study, RNA-sequencing (RNA-seq) technology was used to characterize the difference in liver fat metabolism between Hu sheep and Tibetan sheep. We identified 1179 differentially expressed genes (DEGs) (Q-value < 0.05) between the two sheep breeds, including 25 fat-metabolism-related genes. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, 16 pathways were significantly enriched (Q-value < 0.05), such as the proteasome, glutamatergic synapse, and oxidative phosphorylation pathways. In particular, one of these pathways was enriched to be associated with fat metabolism, namely the thermogenesis pathway, to which fat-metabolism-related genes such as ACSL1, ACSL4, ACSL5, CPT1A, CPT1C, SLC25A20, and FGF21 were enriched. Then, the expression levels of ACSL1, CPT1A, and FGF21 were verified in mRNA and protein levels via qRT-PCR and Western blot analysis between the two sheep breeds. The results showed that the mRNA and protein expression levels of these three genes were higher in the livers of Tibetan sheep than those of Hu sheep. The above genes are mainly related to FAs oxidation, involved in regulating the oxidation of liver FAs. So, this study suggested that Tibetan sheep liver has a greater FAs oxidation level than Hu sheep liver. In addition, the significant enrichment of fat-metabolism-related genes in the thermogenesis pathway appears to be related to plateau-adaptive thermogenesis in Tibetan sheep, which may indicate that liver- and fat-metabolism-related genes have an impact on adaptive thermogenesis.

Reducing fat deposition in sheep (Ovis aries) tails is one of the most important ways to combat rising costs and control consumer preference. Our previous studies have shown that oar-miR-432 is differentially expressed in the tail adipose tissue of Hu (a fat-tailed sheep breed) and Tibetan (a thin-tailed sheep breed) sheep and is a key factor in the negative regulation of fat deposition through BMP2 in ovine preadipocytes. This study investigated the effect of oar-miR-432 and its target genes in ovine preadipocytes. A dual luciferase assay revealed that DDI1 is a direct target gene of oar-miR-432. We transfected an oar-miR-432 mimic and inhibitor into preadipocytes to analyze the expression of target genes. Overexpression of oar-miR-432 inhibits DDI1 expression, whereas inhibition showed the opposite results. Compared with thin-tailed sheep, DDI1 was highly expressed in the fat-tailed sheep at the mRNA and protein levels. Furthermore, we transfected the overexpression and knockdown target genes into preadipocytes to analyze their influence after inducing differentiation. Knockdown of DDI1 induced ovine preadipocyte differentiation into adipocytes but suppressed oar-miR-432 expression. Conversely, the overexpression of DDI1 significantly inhibited differentiation but promoted oar-miR-432 expression. DDI1 overexpression also decreased the content of triglycerides. Additionally, DDI1 is a nested gene in intron 1 of PDGFD. When DDI1 was overexpressed, the PDGFD expression also increased, whereas DDI1 knockdown showed the opposite results. This is the first study to reveal the biological mechanisms by which oar-miR-432 inhibits preadipocytes through DDI1 and provides insight into the molecular regulatory mechanisms of DDI1 in ovine preadipocytes. These results have important applications in animal breeding and obesity-related human diseases.