Association of SNPs within PTPN3 gene with wool production and growth traits in a dual-purpose sheep population.

Protein tyrosine phosphatase non-receptor type 3 (PTPN3), a member of the membrane-associated non-receptor protein tyrosine phosphatase (PTP) family, plays significant roles in the cytoplasm and affects the development and growth of skin and hair. A recent study identified the PTPN3 as the potential gene related to sheep wool quality. To detect single-nucleotide polymorphisms (SNPs) of PTPN3 and elucidate its association with wool production and growth traits in fine wool sheep a total of 644 healthy SG (South African mutton merino♂ × Gansu alpine fine-wool sheep♀, SG) and SSG (South African mutton merino♂ × SG♀, SSG) hybrid sheep were selected. Pooled-DNA sequencing and SNPscan methods were used to scan and genotype SNPs within PTPN3. Association analyses between SNPs and wool production and growth traits were implemented. Consequently, the results revealed that PTPN3 has six SNPs (two missense mutations, one synonymous mutation, and three intron mutations), of which four loci (SNP2, SNP3, SNP4, and SNP5) were significantly positively correlated with growth and wool traits (p < 0.05). SNP4 was significantly (p < 0.05) linked with thigh wool length, and SNP6 was significantly (p < 0.05) associated with abdomen wool length. Moreover, one strongly linked SNP block was identified to be correlated with wool production and growth traits (body weight and body size). The significant SNPs founded by this study could serve as useful genetic markers for breeding fine-wool sheep.®.

Developmental changes of nutrient digestion in young lambs are influenced by weaning and associated with intestinal microbiota.

Understanding the effects of weaning on the changes of digestive function could help to design efficient intervention strategies for promoting the development of the lamb during the early stages of life. In this study, 24 lambs were divided into two groups (control group, lambs were not weaned; and weaning group, lambs were weaned at 21 days of age). The growth, nutrient digestion, gastrointestinal enzyme activity, plasma biochemical indicators, and intestinal microbiota at 7-49 days were determined, as well as the impact of early weaning. The nutrient digestion changed rapidly with age, especially at 14-28 days (p < 0.05). Weaning reduced the dry matter (DM), crude protein (CP), and ether extract (EE) intake and digestion, but increased the starch, neutral detergent fiber (NDF), and acid detergent fiber (ADF) intake and digestion (p < 0.05). Weaning did not affect the overall jejunal microbiota (p > 0.05), but affected the relative abundance of certain bacteria taxa (p < 0.05). Lactic acid-producing bacteria, such as Olsenella, Bacillus, Sharpea, and Bifidobacterium are closely related to CP or EE digestion and growth performance (p < 0.05). In summary, we delineated the pattern of nutrient digestion and intestinal microbiota development in young lambs, and the impact of early weaning.

Differences of the anti-oxidative capability, GPX3, and Cu/ZnSOD expression in Hu sheep testis with different size at six-month-old.

This study aimed to investigate the differences in the anti-oxidant capabilities and related gene expressions of six-month-old Hu sheep with different testis sizes. A total of 201 Hu ram lambs were fed up to 6 months in the same environment. Based on their testis weight and sperm count, 18 individuals were selected and divided into large (n = 9) and small (n = 9) groups, with an average testis weight of 158.67 g ± 5.21 g and 44.58 g ± 4.14 g, respectively. The total antioxidant capacity (T-AOC), total superoxide dismutase (T-SOD), and malondialdehyde (MDA) concentration in testis tissue were tested. The localization of antioxidant-related genes, GPX3 and Cu/ZnSOD in testis were detected by immunohistochemistry. The GPX3, Cu/ZnSOD expression, and relative mitochondrial DNA (mtDNA) copy number were detected by quantitative real-time PCR. Compared with the small group, the T-AOC (2.69 ± 0.47 vs. 1.16 ± 0.22 U/mgprot) and T-SOD (22.35 ± 2.59 vs. 9.92 ± 1.62 U/mgprot) in the large group were significantly higher, whereas the MDA (0.72 ± 0.13 vs. 1.34 ± 0.17 nM/mgprot) and relative mtDNA copy number in the large group was significantly lower (p < .05). Immunohistochemistry results indicated that the GPX3 and Cu/ZnSOD were expressed in Leydig cells and seminiferous tubule. The expressions of GPX3 and Cu/ZnSOD mRNA in the large group were significantly higher than those in the small group (p < .05). In conclusion, Cu/ZnSOD and GPX3 widely expressed in the Leydig cells and seminiferous tubule, high expression of Cu/ZnSOD and GPX3 in a large group has a higher potential in addressing oxidative stress and contribute to spermatogenesis.

Transcriptome profiling reveals differential gene expression in the rumen of Hu lambs at different developmental stages.

The development of the rumen is a critical physiological challenge in newborn ruminants. However, the molecular mechanism underlying different stages of rumen development in sheep remains poorly understood. Here, RNA sequencing and bioinformatics analysis were performed to compare the transcription profiles of rumen development at 7, 28 and 56 days of birth (D7, D28 and D56). We identified 1246, 2257 and 627 differentially expressed genes (DEGs) between D7 and D28, between D7 and D56, between D28 and D56, respectively. Also, 70 DGEs were co-expressed at these three time points. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated most DEGs mainly related to transporter activity, channel activity and metabolism pathways. Noteworthy, the expression levels of most genes (CA4, CA9, CA12 and CA14) in nitrogen metabolic pathways were negatively correlated with the papilla length and width, but the papilla length and width were positively correlated with the expression of genes (PLA2G3, SLC26A9, SLC34A3) in ion transport pathway, suggesting that these genes may be involved in nitrogen metabolic and ion transport pathway and thus affect rumen development. These results provide new insight into the changes in RNA expression at different time points of rumen development in Hu sheep.

Transcriptomic analysis reveals that non-forage or forage fiber source promotes rumen development through different metabolic processes in lambs.

Dietary fiber supplementation can stimulate rumen development in lambs during the pre-weaning period. However, it is unclear whether different sources of fiber have varying effects on rumen development. This study aimed to investigate the molecular mechanism of rumen morphological and functional development based on non-forage or forage as a starter dietary fiber source. Twenty-four male Hu lambs with similar body weights (BW, 3.67 ± 0.08 kg) were selected and divided into two groups that received diets supplemented with either alfalfa hay (AH) or soybean hull (SH). At the age of 70 days, six lambs were slaughtered from each treatment group for rumen fermentation and morphological analyses. Three samples of the rumen tissue from the ventral sac were collected for transcriptomic analysis. The results identified 633 differentially expressed genes (DEGs), of which 210 were upregulated and 423 were downregulated in the SH group compared with those in the AH group. The upregulated DEGs were most enriched in the immune function and proteolysis pathways, whereas the downregulated DEGs were mainly involved in cell proliferation, apoptosis, and differentiation pathways. These findings indicated that non-forage as a starter dietary fiber source improved immune function and enhanced nitrogen utilization, whereas forage facilitated rumen morphological development.

Variants of ADPGK gene and its effect on the male reproductive organ parameters and sperm count in Hu sheep.

ADP-dependent glucokinase (ADPGK) plays an important role instead of hexokinase in regulating energy metabolism via the Embden-Meyerhof-Parnas Pathway. And energy provided via glycolysis promotes testis development and spermatogenesis. In this study, 466 Hu sheep were screened for mutations in the ADPGK gene to examine the association of the ADPGK gene polymorphisms with the testis traits and spermatogenesis. The NC_056060.1: g.31295 C > T SNP was found in the 3'-UTR region, resulting in two genotypes CC and TC type with genotypic frequencies of 0.66 and 0.34, respectively. This mutation was significantly associated with testis weight, testis long circumference, testis short girth, epididymis weight, and sperm concentration (p < 0.05). Moreover, TC genotype individuals had an increased tendency in the expression of the ADPGK gene and had significant reproductive performance advantages compared with CC genotype individuals in the study. And compared with the small testes (<50 g), the ADPGK gene expression of big testes (>160 g) increased significantly. This indicates an association between the ADPGK gene and reproductive organ parameters and sperm count in selected Hu sheep breed, and this SNP may serve as an effective DNA molecular marker for marker-assisted selection in Hu sheep breeding programs.

Effects of starter feeding and early weaning on developmental expressions of IGF-I gene in liver and IGF-IR gene in rumen of lambs from birth to eighty-four days of age.

The early weaning and starter feeding have significant effects on lamb growth and digestive tract development. However, it is not clear whether the expression of IGF-1/IGF-1R genes could be affected by feeding starter and weaning. Therefore, a total of 102 Hu male neonatal lambs were randomly divided into 3 groups: fed with starter at 42d + weaned at 56d (group A), fed with starter at 7d + weaned at 28d (group B) and fed with starter at 7d + weaned at 56d (group C), to explore the effects of starter feeding and weaning age on developmental expressions of IGF-I gene in liver and IGF-IR gene in rumen of Hu sheep. The results showed that IGF-I and IGF-IR genes were expressed extensively in various tissues of lambs, the expression of IGF-I was significantly higher in liver (p < 0.01), while the expression of IGF-IR was higher in rumen among gastrointestinal tissues. The mRNA level of IGF-I of group C was strikingly higher than that of group A at 28, 70, and 84d (p < 0.01) in liver, respectively. Rumen mRNA level of IGF-IR of group C was prominently higher than that of group A at 70d and 84d (p < 0.05), respectively. However, IGF-IR expression of group A was higher than that of group C at 14 and 42d (p < 0.05). The mRNA level of IGF-I of group B was significantly higher than that of group C at 42, 56, 70, and 84d (p < 0.01) in liver, while IGF-IR expression of group B was higher than that of group C at 70d (p < 0.01). In conclusion, early weaning and starter feeding affected liver IGF-I mRNA and rumen IGF-IR mRNA expression.

Effect of feed efficiency on growth performance, body composition, and fat deposition in growing Hu lambs.

The present study aimed to investigate the relationship between growth performance, body composition, and fat deposition factors, and feed efficiency in growing lambs. We measured average daily feed intake (ADFI) and body weight (BW) from 653 Hu sheep that were fed a pellet diet. The residual feed intake (RFI) not significantly genetic and phenotypic correlated with the metabolic body weight (MBW) and average daily gain (ADG), but it was significantly genetic and phenotypic correlated with ADFI and the feed conversion ratio (FCR) (p < 0.01). However, the FCR was significantly associated with growth traits (p < 0.01). With the same ADG, body fat deposition was greater in animals with low feed efficiency compared with high feed efficiency. Therefore, excessive fat deposition can affect the feed efficiency of the body, and organ weight and gut-weight have a greater impact on the feed efficiency of lambs. The reticulum stomach and jejunum of lambs with a low RFI were smaller compared with that in the high RFI, indicating that lambs with a low RFI have less intake and a higher absorption rate. Small organs, such as the liver, of lambs with high FE might be associated with low energy expenditure and slow metabolism. This study provides a new perspective to study the biological processes responsible for feed efficiency variation in lambs.

Polymorphism and expression of the HMGA1 gene and association with tail fat deposition in Hu sheep.

Hu sheep is an excellent short fat-tailed breed in China. Fat deposition in Hu sheep tail affects carcass quality and consumes a lot of energy, leading to an increase in feed cost. The objective of this study was to analyze the effects of HMGA1 polymorphism on tail fat weight in Hu sheep. Partial coding and non-coding sequences of HMGA1 were amplified with PCR and single nucleotide polymorphisms (SNP) of HMGA1 in 1163 Hu sheep were detected using DNA sequencing and KASPar technology. RT-qPCR analysis was performed to test HMGA1 expression in different tissues. The results showed that the expression of HMGA1 was higher in the duodenum, liver, spleen, kidney, and lung than in the heart, muscle, rumen, tail fat, and lymph. A mutation, g.5312 C > T, was detected in HMGA1; g.5312 C > T was significantly associated with tail fat weight, relative weight of tail fat (body weight), and relative weight of tail fat (carcass) (p < 0.05). The tail fat weight of the TT genotype was remarkably higher than that of the CC and TC genotypes. Therefore, HMGA1 can be used as a genetic marker for marker-assisted selection of tail fat weight in Hu sheep.

Polymorphisms in the ovine GP5 gene associated with blood physiological indices.

This study aimed to analyze the effects of polymorphisms in GP5 on blood physiological indices of 1065 sheep. The coefficients of variation of the red blood cell count (RBC), hemoglobin concentration (HGB), mean platelet volume (MPV), mean erythrocyte hemoglobin content (MCH), and red blood cell distribution-coefficient of variation (RDW-CV) were greater than 10%, and there was a very significant correlation between the main indices such as RBC, white blood cell, and platelet count (PLT) and most other indices (p < 0.01). qRT-PCR showed that GP5 was expressed in the heart, liver, spleen, lung, kidney, rumen, duodenum, muscle, tail fat, and lymph tissue, with significantly higher expression in the lymph. Subsequently, we detected single nucleotide polymorphisms (SNPs) in GP5 from group, which identified synonymous mutation g.657 T > C in the first exon of GP5. Association analysis showed significant correlations between the SNP and the physiological traits (p < 0.05), in which the RBC, neutrophilic granulocyte (NEUT) and RDW-CV values in sheep with the TC genotype and TT genotype were markedly lower than those in the CC genotype (p < 0.05). Thus, GP5 polymorphisms could be candidate biomarkers to screen blood physiological indices.

Expression and polymorphisms of CD8B gene and its associations with body weight and size traits in sheep.

The growth traits are economically important traits in sheep. Improving growth rates will increase the profitability of producers. The aim of this study was to identify alleles of CD8B (encoding T-cell surface glycoprotein CD8 beta chain) that are aberrantly expressed in different tissues and to assess the effects and associations of its different genotypes on weight and size traits in sheep. Using quantitative real-time reverse transcription PCR arrays, expression profiling of CD8B was performed in various organs and tissues. CD8B was ubiquitously expressed, with very high expression in the lung, spleen, lymph, duodenum, and liver. One intronic mutation (chr3:62,718,030 (Oar_rambouillet_v1.0, same below) G > A) was identified using pooled DNA sequencing. Subsequently, the variants (AA, AG, and GG) were genotyped using the KASPar® PCR single nucleotide polymorphism (SNP) genotyping system. The results of association analysis with body weight and body size traits in 1304 sheep showed that increases in multiple phenotypic traits correlated with the AA genotype (body weight, p < 0.05; body length, p < 0.05). Thus, SNP chr3:62,718,030 G > A is a promising molecular marker for marker-assisted selection in sheep breeding.

Polymorphism of sheep PRKAA2 gene and its association with growth traits.

Small ruminants farming plays an important role in the livelihood of a large part of the population. Herein we aimed to analyze the effects of single nucleotide polymorphisms in PRKAA2 gene on the growth-related traits of Hu sheep and Dorper sheep. The body weight and body type of 1254 sheep were measured at 80, 100, 120, 140, 160 and 180d, and 37620 phenotypic data were collected. RT-qPCR analysis was performed to test PRKAA2 gene expressed in different tissues of sheep, with the highest expression level in spleen, followed by kidney. In the present study, the PRKAA2 gene sequencing revealed one polymorphism located on Chr1 (Oar_rambouillet_v1.0), termed as chr1:32832382 G > A, and were significantly associated with growth traits of sheep (p < 0.05). The body weight, body length, chest circumference, and cannon circumference of individuals with AA genotype were significantly higher than those with the GG and GA genotypes (p < 0.05). Our findings reveal that PRKAA2 gene could be used as a marker-assisted selection to improve the growth-related traits of sheep.

Estimation of genetic correlations of two key feed efficiency traits with production traits in male Hu sheep.

Increased body growth and feed efficiency are important for breeding due to expensive feed costs. This study estimated the genetic parameters for two feed efficiency traits [the feed conversion ratio (FCR) and the residual feed intake (RFI)] and other important economic traits in male Hu sheep using 1642 male Hu lambs. The heritability of the RFI was estimated at 0.31 (±0.10)-0.54 (±0.1) and the heritability of the FCR was 0.05 (±0.07)-0.77 (±0.12). There was a negative genetic correlation between rib-eye muscle area traits and RFI. A positive genetic correlation was found between fat deposition traits and feed efficiency traits. The digestive tract weight decreased with the direction of RFI: however, the FCR showed the opposite: A low FCR was accompanied by a high digestive tract weight. The genetic correlation between feed efficiency (RFI and FCR) and rumen weight was negative. Heritability estimates for feed efficiency generate fluctuations with different experimental intervals. The weight of other traits should be reasonably allocated to avoid losses in production and feed efficiency in the selection of feed efficiency traits. The present study added to our understanding of the genetic parameters of sheep under the condition of house feeding.

Identification and characterization of circular RNAs in association with the feed efficiency in Hu lambs.

BACKGROUND: Circular RNA (circRNA), as a new members of noncoding RNA family, have vital functions in many biological processes by as microRNA sponges or competing endogenous RNAs (ceRNAs). However, little has been reported about the genetic mechanism of circRNAs regulation of feed efficiency in sheep. RESULTS: This study aimed to explore the expression of circRNAs in the liver of Hu sheep with High-RFI (High residual feed intake) and Low-RFI (Low residual feed intake) using transcriptome sequencing. A total of 20,729 circRNAs were identified in two groups, in which 219 circRNAs were found as significantly differentially expressed. Several circRNAs were validated by using RT-PCR, sanger sequencing and RT-qPCR methods. These results demonstrated that the RNA-seq result and expression level of circRNAs identified are reliable. Subsequently, GO and KEGG enrichment analysis of the parental genes of the differentially expressed (DE) circRNAs were mainly involved in immunity response and metabolic process. Finally, the ceRNA regulatory networks analysis showed that the target binding sites for miRNA such as novel_41, novel_115, novel_171 and oar-miR-485-3p in the identified DE cirRNAs. Importantly, two metabolic (SHISA3 and PLEKHH2) and four (RTP4, CD274, OAS1, and RFC3) immune-related target mRNAs were identified from 4 miRNAs. Association analysis showed that the polymorphism (RTP4 c.399 A > G) in the target gene RTP4 were significantly associated with RFI (P < 0.05). CONCLUSIONS: Analysis of sequencing data showed some candidate ceRNAs that may play key roles in the feed efficiency in sheep by regulating animal immune and metabolic. These results provide the basis data for further study of the biological functions of circRNAs in regulating sheep feed efficiency.

A mutation modulating DDX3Y gene expression cosegregates with the major Y-chromosomal haplogroups and with testis size in Hu sheep.

Variations in the Y-chromosome are usually correlated with male-specific traits. However, this condition has been described only sporadically, even in human genetics. The present study was conducted to clone the full-length gene sequence of ovine DEAD-box helicase 3, Y-linked (DDX3Y), and investigate the effect of the expression and variation of DDX3Y on the reproductive traits of Hu sheep. Consequently, we identified the full coding sequence and genomic sequence of ovine DDX3Y. Quantitative PCR (qPCR) analysis showed that ovine DDX3Y was highly expressed in testis, and the expression level increased during testicular development. Furthermore, individuals with larger testis at 6 months expressed significantly more DDX3Y mRNA in the testis than individuals with smaller testis. Notably, a novel SNP (g. 12657 C>A) in the 3' untranslated region was identified in Hu sheep and Tan sheep according to the investigation of the full DDX3Y genomic sequence of 1069 individuals from nine sheep breeds. Association analysis revealed that the SNP was significantly related to testis size in Hu sheep. Meanwhile, Hu rams with the derived C allele showed significantly higher expression levels of DDX3Y in testis than those with the ancestral A allele. In addition, data mining in a previous study showed that the C allele cosegregated with the globally major Y-chromosomal haplogroups y-HA and y-HC, and the A allele is found in all rams with haplogroups y-HB1, y-HB2 and y-HD. This study suggests that the association of the Y-chromosomal haplogroups with testis size in Hu sheep can be extrapolated to the sheep population worldwide.

α-Linolenic acid induced TM4 Sertoli cells proliferation and enhanced total antioxidant capacity.

The present study was undertaken to elucidate the direct ALA effects on mice TM4 Sertoli cells proliferation in vitro. Our results showed that TM4 cells viability was significantly stimulated by ALA (p < 0.05). The 50 µM ALA increased the concentration of total antioxidant capacity, induced the mitochondrial membrane hyperpolarized, and markedly decreased the number of apoptosis cells (p < 0.05). ALA also up-regulated G2/Mitotic-specific cyclin-B1 gene and apoptosis suppressive gene Bcl2 expression (p < 0.05). In conclusion, those results indicated that ALA could increase TM4 Sertoli cells antioxidant capacity, induced the mitochondrial membrane hyperpolarized, inhibited cells apoptosis and stimulated TM4 Sertoli cells proliferation in vitro.

Expression profile of FASN gene and association of its polymorphisms with intramuscular fat content in Hu sheep.

The content of intramuscular fat (IMF) is one of the most important factors that has a large impact on meat quality, and it is an effective way to improve IMF according to marker-assisted selection (MAS). Fatty-acid synthase (FASN) is a key gene in meat lipid deposition and fatty acid composition. Thus, this study was conducted to investigate the expression profile of FASN in mRNA and protein levels using real-time quantitative PCR (RT-qPCR) and western-blot methods. In addition, single nucleotide polymorphisms (SNPs) within FASN in 921 Hu rams with IMF content records were investigated using DNA-pooling sequencing and improved multiple ligase detection reaction (iMLDR) methods. Consequently, the highest mRNA expression level of FASN was observed in the perinephric fat, and the lowest in the liver among the 11 tissues analyzed, while no significant difference was found in mRNA and protein expression levels in longissimus dorsi among individuals with different IMF contents. A total of 10 putative SNPs were identified within FASN, and 9 of them can be genotyped by iMLDR method. Notably, two SNPs were significantly associated with IMF content, including NC_040262.1: g.5157 A > G in intron 5 (p = 0.046) and NC_040262.1: g.9413 T > C in intron 16 (p = 0.041), which supply molecular markers for improving meat quality in sheep breeding.

The functional development of the rumen is influenced by weaning and associated with ruminal microbiota in lambs.

Rumen development is critical for the development of early lambs. This work aims to evaluate the effects of abrupt weaning at day 21 on rumen fermentation, histomorphological traits and the ruminal microbiota compared with continuous suckling. Twelve pairs of artificially reared full-sib neonatal male Hu lambs were allocated to two groups, one of which was weaned at day 21 (EW group) and the other which was not weaned (CON group). At day 26 and day 49, six lambs from each group were randomly selected and sacrificed to collect ruminal contents and rumen tissue samples. Results showed that weaning influenced the fermentation parameters in the rumen, and altered the microbial community composition on day 49 (p < 0.05). Several genera were associated with rumen fermentation parameters (p < 0.05). Volatile fatty acid (VFA) concentration is the key parameter impacting microbiota composition. Weaning influenced the expression of genes associated with VFA metabolism and regulation of cell proliferation (p < 0.05). In conclusion, weaning significantly influenced the morphological and functional development of the rumen, and bacterial community composition. The microbial community composition was strongly associated with rumen weight and fermentation profiles, but not with morphological development.

The association of polymorphisms in the ovine PPARGC1B and ZEB2 genes with body weight in Hu sheep.

The aims of this study were to analyze the effects of PPARGC1B and ZEB2 polymorphisms on the body weight of Hu sheep. DNA sequencing and KASPar technologies were used to detect single nucleotide polymorphisms (SNPs) within the PPARGC1B and ZEB2 genes of Hu sheep (n = 207). Two synonymous mutations, g.300 G > A and g.645 C > T, were detected in PPARGC1B and ZEB2, respectively. The body weights of sheep were recorded at 80, 100, 120, 140, 160 and 180-days, and associations between these polymorphisms and body weight changes were analyzed. Association analysis demonstrated that the polymorphisms in PPARGC1B and ZEB2 significantly associated with body weight (p < 0.05). At the g.300 G > A locus, individuals with the GG genotype had significantly higher body weight than those with the AA genotype, and at the g.645 C > T locus, individuals with the TT genotype had significantly higher body weight than those with the TC genotype. Individuals with both polymorphisms exhibited significantly different growth (p < 0.05). These data suggest that polymorphisms in the PPARGC1B and ZEB2 genes can be used as candidate molecular markers for the breeding of desirable growth traits in Hu sheep.

Whole-genome resequencing identified candidate genes associated with the number of ribs in Hu sheep.

The number of ribs is an important economic trait in the sheep industry when the sheep are raised for mutton. However, in sheep, the genetic mechanisms regulating rib number are poorly understood. In the present study, we aimed to identify important candidate genes that affect the increase in rib number in sheep. Whole-genome resequencing of 36 Hu sheep with an increased number of ribs (R14) and 36 sheep with normal (R13) rib numbers was carried out. Analysis using three methods (fixation index (FST), Fisher's exact test, and Chi-squared test) showed that 219 single nucleotide polymorphism sites overlapped among the results of the three methods, which represented 206 genes. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that the genes were mainly associated with regulation of developmental process, inorganic anion transport, cellular biosynthetic process, tight junction, the oxytocin signaling pathway, and arrhythmogenic right ventricular cardiomyopathy. Four mutations were selected according to the significantly selected genomic regions and important pathways for genotyping and association analysis. The result demonstrated that three synonymous mutations correlated significantly with the rib number. Importantly, we revealed that the CPOX (encoding coproporphyrinogen oxidase), KCNH1 (encoding potassium voltage-gated channel subfamily H member 1), and CPQ (encoding carboxypeptidase Q) genes have a combined effect on rib number in Hu sheep. Our results identified candidate molecular markers for rib number in sheep breeding.