Search

Tibetan sheep have developed unique adaptations for survival in the Qinghai–Tibet Plateau environment. However, the functional genes and molecular mechanisms that regulate hypoxia adaptation have not been fully characterized. In this study, based on the whole-genome resequencing data for Tibetan sheep at different altitudes, the population differentiation index (FST) and nucleotide diversity ratio (θπ ratio) were evaluated in populations of 20 Oula sheep (3501 m altitude, OL), 20 Zashijia sheep (4369 m altitude, ZSJ), and 20 Awang sheep (4643 m altitude, AW) to reveal candidate loci related to high-altitude hypoxia. We screened 728 and 524 candidate genes in the AW vs. OL and ZSJ vs. OL groups, respectively, of which 134 genes were jointly screened. Candidate genes were mainly enriched in the Ras, melanoma, melanogenesis, VEGF, and PPAR signaling pathways. HIF1AN, PDGFA, PDGFD, ANXA2, SOCS2, NOXA1, WNT7B, MMP14, GNG2, ATF6, PGAM2, PPP3R1, GSTCD, and PPARA may play important roles in the high-altitude adaptation of Tibetan sheep. In conclusion, this study provides valuable insights into the genes and molecular mechanisms that underlie high-altitude hypoxia adaptation in Tibetan sheep.

Numerous genes govern male reproduction, modulating testicular development and spermatogenesis. Our study leveraged RNA-Seq to explore candidate genes and pivotal pathways influencing fecundity in an F1 hybrid of Southdown × Hu sheep testes across four developmental milestones: M0 (0 months old, newborn), M3 (3 months old, sexually immature), M6 (6 months old, sexually mature), and Y1 (1 years old, adult). Histological examination using hematoxylins and eosin staining revealed that the cross-sectional area of the spermatid tubules and the number of supportive cells increased in the other groups, as compared to the M0 group. The cross-sectional area of the vasculature and the number of supporting cells were found to be significantly increased in all other groups in comparison to the M0 group. We conducted GO and KEGG analyses of the differentially expressed genes (DEGs) in the three comparison groups and identified key pathways, including cAMP, MAPK, ECM–receptor interactions, PI3K-Akt, and FOXO signaling, which are closely related to testicular development and spermatogenesis. Notably, alternative splicing (AS) events were markedly elevated in M6 and Y1 stages. Key genes like GATA4, GATA6, SMAD4, SOX9, YAP1, ITGB1 and MAPK1 emerged as significantly enriched in these pathways, potentially orchestrating the transition from immature to mature testes in sheep. These findings offer valuable insights into male reproductive potential and can inform strategies for optimizing animal breeding.

SNPs associated with important traits of fine-wool sheep that were previously obtained through genome-wide association analysis screening were verified and analyzed. A total of 20 SNPs related to birth weight, bundle strength, cleaning rate, and fiber diameter were screened using whole-genome resequencing, and the SNPshot assay was used to detect and analyze polymorphisms. This study found that, among the 20 SNPs associated with important traits in Alpine Merino sheep, 8 were monomorphic and 12 were polymorphic, of which 6 showed moderate polymorphisms and 6 showed low polymorphisms. The heterozygosity of the 12 polymorphic loci ranged from 0.10 to 0.49, the effective number of alleles ranged from 1.11 to 1.98, and the polymorphic information content ranged from 0.09 to 0.37. The chi-square test showed that only RHPN2:g.42678119T>G and RALYL:g.90030866A>G were in Hardy–Weinberg disequilibrium (p < 0.05); the other loci were in equilibrium (p > 0.05). These SNPs associated with important traits in Alpine Merino sheep provide a theoretical basis for genomic selection and molecular design breeding.

The SLIT3 gene encodes a secreted protein, and the ZNF280B gene is a member of the transcription factor family. Both genes have multiple biological functions. This study was conducted to investigate the association between SLIT3 and ZNF280B gene polymorphisms and wool fiber diameter and to determine potential molecular marker sites for breeding sheep with fine wool. We used Kompetitive Allele-Specific PCR to type the single nucleotide polymorphism (SNP) loci in the SLIT3 and ZNF280B genes within 1081 Alpine Merino sheep and associated these SNPs with wool fiber diameter. The results revealed one SNP in SLIT3 and ZNF280B, which were each related to sheep fiber diameter. The wool fiber diameters of sheep with the CC genotype in SLIT3 g.478807C>G and AA genotype in ZNF280B g.677G>A were the smallest and differed significantly from the diameters of other genotypes (p < 0.05). These results suggest potential molecular marker sites for fine-wool sheep breeding.

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.

This experiment was conducted to study the effects of dietary energy level on the growth performance and meat quality of weaned Alpine Merino lambs. The study ran for a total of 104 days (20-day pretrial, 84-day trial). From three groups of test lambs, we randomly selected ten lambs per group to determine slaughter performance, meat quality characteristics, and organ indexes. The slaughter performances of the lambs improved as the dietary energy level increased. The live weight before the slaughter of the lambs was significantly higher in the high group than in the low and medium groups. The carcass weight was significantly higher in the high group than in the low group. Dietary energy level had little effect on the organ weight of lambs. Meat quality differed among the three dietary energy levels. The muscle yellowness and redness scores decreased significantly as the energy levels increased. The C18:0, C21:0, C20:1, C18:2n6c, and C20:2 contents in the muscle were significantly higher in the high group than in the medium and low groups. The C18:3n6 content in the muscle was significantly higher in the low group than in the medium group. The C20:5n3 content in the longissimus dorsi was significantly higher in the high group than in the medium and low groups. The monounsaturated and unsaturated fatty acid contents in the muscle were significantly higher in the high group than in the low group. A dietary energy level of 10.5 MJ/kg is suitable for fattening weaned male Alpine Merino lambs.

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.

Inner Mongolia and Liaoning cashmere goats in China are well-known for their cashmere quality and yield. Thus, they are great models for identifying genomic regions associated with cashmere traits. Herein, 53 Inner Mongolia cashmere goats, Liaoning cashmere goats and Huanghuai goats were genotyped, and 53,347 single-nucleotide polymorphisms (SNPs) were produced using the Illumina Caprine 50K SNP chip. Additionally, we identified some positively selected SNPs by analyzing Fst and XP-EHH. The top 5% of SNPs had selection signatures. After gene annotation, 222 and 173 candidate genes were identified in Inner Mongolia and Liaoning cashmere goats, respectively. Several genes were related to hair follicle development, such as TRPS1, WDR74, LRRC14, SPTLC3, IGF1R, PADI2, FOXP1, WNT10A and CSN3. Gene enrichment analysis of these cashmere trait-associated genes related 67 enriched signaling pathways that mainly participate in hair follicle development and stem cell pluripotency regulation. Furthermore, we identified 20 overlapping genes that were selected in both cashmere goat breeds. Among these overlapping genes, WNT10A and CSN3, which are associated with hair follicle development, are potentially involved in cashmere production. These findings may improve molecular breeding of cashmere goats in the future.

Tibetan goat is an ancient breed, which inhabits the adverse conditions of the plateaus in China. To investigate the role of selection in shaping its genomes, we genotyped Tibetan goats (Nagqu Prefecture, above 4500 m) and three lowland populations (Xinjiang goats, Taihang goats and Huanghuai goats). The result of PCA, neighbor-joining (N-J) tree and model-based clustering showed that the genetic structure between the Tibetan goat and the three lowland populations has significant difference. As demonstrated by the di statistic, we found that some genes were related to the high-altitude adaptation of Tibetan goats. Functional analysis revealed that these genes were enriched in the VEGF (vascular endothelial growth factor) signaling pathway and melanoma, suggesting that nine genes (FGF2, EGFR, AKT1, PTEN, MITF, ENPEP, SIRT6, KDR, and CDC42) might have important roles in the high-altitude adaptation of Nagqu Tibetan goats. We also found that the LEPR gene was under the strongest selection (di value = 16.70), and it could induce upregulation of the hypoxic ventilatory response. In addition, five genes (LEPR, LDB1, EGFR, NOX4 and FGF2) with high di values were analyzed using q-PCR. Among them, we found that LEPR, LDB1 and FGF2 exhibited higher expression in the lungs of the Tibetan goats; LEPR, EGFR and LDB1 exhibited higher expression in the hearts of the Huanghuai goat. Our results suggest that LEPR, LDB1, EGFR and FGF2 genes may be related to the high-altitude adaptation of the goats. These findings improve our understanding of the selection of the high-altitude adaptability of the Nagqu Tibetan goats and provide new theoretical knowledge for the conservation and utilization of germplasm resources.

To date, molecular genetics and population studies in Tibetan sheep (Ovis aries) have been limited, and little is known about the phylogenetic evolution and phylogeography of Tibetan sheep populations. The aim of the present research was to explore phylogeography and phylogenetic evolution of Tibetan sheep populations, on the basis of mitochondrial DNA (mtDNA) gene MT-CYB (1140 base pairs). Our dataset consisted of 641 MT-CYB sequences from the same amount of animals belonging to 15 populations of Tibetan sheep living in the Qinghai–Tibetan Plateau, China. Haplotype and nucleotide diversities were 0.748 ± 0.010 and 0.003 ± 0.001, respectively. The analysis of phylogeography revealed the presence of two formerly described haplogroups in 15 populations of Tibetan sheep, however only one haplogroup was present in Awang sheep. Moreover, 641 Tibetan sheep were distributed into a minimum of two clusters by clustering analysis. The 15 Tibetan sheep populations and 19 reference populations of 878 individuals were separated into six main groups based on their substitutions per site, from which we constructed a phylogenetic tree. Minor differences in branching order of various taxa between trees acquired from either gene were observed. This study provides insights on the origins and phylogenetic evolution of populations residing in the Qinghai–Tibetan Plateau, which will aid information of future conservation programs aimed at conserving this valuable genetic resource.

Bone morphogenetic protein 2 (BMP2) is strongly selected in both fat-tailed and thin-tailed sheep and may be a candidate gene for sheep tail type selection. However, the mechanism of action of BMP2 in sheep tail fat deposition remains unclear. This study investigated genetic variation and haplotype combinations of the BMP2 gene in sheep with different tail types, aiming to reveal the molecular mechanism of BMP2 in sheep tail fat deposition. We detected a total of three single nucleotide polymorphisms (SNPs) (g.48401619 T > A, g.48401272 C > A, and g.48401136 C > T) among 533 sheep. The alleles and genotype frequencies of these SNPs were in Hardy–Weinberg equilibrium and showed significant correlations with tail length. Linkage disequilibrium existed between the g.48401272 C > A and g.48401136 C > T sites, where CACT was the predominant genotype. At the cellular level, the expression levels of peroxisome proliferator-activated receptor gamma (PPARγ) and lipoprotein lipase (LPL) were upregulated after BMP2 overexpression; there were significantly higher levels of PPARγ than controls at 0 d and 1 d, and of LPL than controls at 1 d and 7 d. These results indicate that the BMP2 gene may participate in sheep tail fat deposition and could be used for molecular-marker-assisted selection of sheep tail type.

Body weight is an important economic trait for sheep and it is vital for their successful production and breeding. Therefore, identifying the genomic regions and biological pathways that contribute to understanding variability in body weight traits is significant for selection purposes. In this study, the genome-wide associations of birth, weaning, yearling, and adult weights of 460 fine-wool sheep were determined using resequencing technology. The results showed that 113 single nucleotide polymorphisms (SNPs) reached the genome-wide significance levels for the four body weight traits and 30 genes were annotated effectively, including AADACL3, VGF, NPC1, and SERPINA12. The genes annotated by these SNPs significantly enriched 78 gene ontology terms and 25 signaling pathways, and were found to mainly participate in skeletal muscle development and lipid metabolism. These genes can be used as candidate genes for body weight in sheep, and provide useful information for the production and genomic selection of Chinese fine-wool sheep.

The aim of this study was to examine the correlation between the platelet-derived growth factor-D (PDGF-D) gene and sheep tail type character and explore the potential underlying mechanism. A total of 533 sheep were included in this study. Polymorphic sites were examined by Pool-seq, and individual genotype identification and correlation analysis between tail type data were conducted using the matrix-assisted laser desorption/ionization time-of-flight mass spectrometer (MALDI-TOF-MS) method. JASPART website was used to predict transcription factor binding sites in the promoter region with and without PDGF-D gene mutation. The effect of PDGF-D on adipogenic differentiation of sheep preadipocytes was investigated. Two single nucleotide polymorphism sites were identified: g.4122606 C > G site was significantly correlated with tail length, and g.3852134 C > T site was significantly correlated with tail width. g.3852134 C > T was located in the promoter region. Six transcription factor binding sites were eliminated after promoter mutation, and three new transcription factor binding sites appeared. Expression levels of peroxisome proliferator-activated receptor gamma (PPARγ) and lipoproteinlipase (LPL) were significantly up-regulated upon PDGF-D overexpression. Oil red O staining showed increased small and large oil drops in the PDGF-D overexpression group. Together these results indicate the PDGF-D gene is an important gene controlling sheep tail shape and regulating sheep tail fat deposition to a certain degree.

With the intensified and large-scale development of sheep husbandry and global warming, sheep heat stress has become an increasingly important issue. However, little is known about the molecular mechanisms related to sheep responses to heat stress. In this study, transcriptomic analysis of liver tissues of sheep in the presence and absence of heat stress was conducted, with the goal of identifying genes and pathways related to regulation when under such stress. After a comparison with the sheep reference genome, 440,226,436 clean reads were obtained from eight libraries. A p-value ≤ 0.05 and fold change ≥ 2 were taken as thresholds for categorizing differentially expressed genes, of which 1137 were identified. The accuracy and reliability of the RNA-Seq results were confirmed by qRT-PCR. The identified differentially expressed genes were significantly associated with 419 GO terms and 51 KEGG pathways, which suggested their participation in biological processes such as response to stress, immunoreaction, and fat metabolism. This study’s results provide a comprehensive overview of sheep heat stress-induced transcriptional expression patterns, laying a foundation for further analysis of the molecular mechanisms of sheep heat stress.

Clinical mastitis is still an intractable problem for sheep breeding. The natural immunologic mechanisms of the mammary gland against infections are not yet understood. For a better understanding of the disease-associated proteins during clinical mastitis in meat sheep, we performed two-dimensional electrophoresis (2-DE)-based comparative proteomic analyses of mammary tissues, including from healthy mammary tissues (HMTs) and from mammary tissues with clinical mastitis (CMMTs). The 2-DE results showed that a total of 10 up-regulated and 16 down-regulated proteins were identified in CMMTs when compared to HMTs. Of these, Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses revealed that most proteins were associated with immune responses or metabolisms. The results of qRT-PCR and Western blot for randomly selected four differentially expressed proteins (DEPs) including superoxide dismutase [Mn] (SOD2), annexin A2 (ANAX2), keratin 10 (KRT10) and endoplasmic reticulum resident protein 29 (ERP29) showed that their expression trends were consistent with 2-DE results except ANXA2 mRNA levels. This is an initial report describing the 2-DE-based proteomics study of the meat sheep mammary gland with clinical mastitis caused by natural infection, which provides additional insight into the immune and metabolic mechanisms during sheep mastitis.

Long noncoding RNAs (lncRNAs) are increasingly being recognized as key regulators in many cellular processes. However, few reports of them in livestock have been published. Here, we describe the identification and characterization of lncRNAs in ovine skeletal muscle. Eight libraries were constructed from the gastrocnemius muscle of fetal (days 85 and 120), newborn and adult Texel and Ujumqin sheep. The 2002 identified transcripts shared some characteristics, such as their number of exons, length and distribution. We also identified some coding genes near these lncRNA transcripts, which are particularly associated with transcriptional regulation- and development-related processes, suggesting that the lncRNAs are associated with muscle development. In addition, in pairwise comparisons between the libraries of the same stage in different breeds, a total of 967 transcripts were differentially expressed but just 15 differentially expressed lncRNAs were common to all stages. Among them, we found that TCONS_00013201 exhibited higher expression in Ujumqin samples, while TCONS_00006187 and TCONS_00083104 were higher in Texel samples. Moreover, TCONS_00044801, TCONS_00008482 and TCONS_00102859 were almost completely absent from Ujumqin samples. Our results suggest that differences in the expression of these lncRNAs may be associated with the muscular differences observed between Texel and Ujumqin sheep breeds.