Whole-Genome Resequencing of Ujimqin Sheep Identifies Genes Associated with Vertebral Number

The number of vertebrae is a crucial economic trait that can significantly impact the carcass length and meat production in animals. However, our understanding of the quantitative trait loci (QTLs) and candidate genes associated with the vertebral number in sheep (Ovis aries) remains limited. To identify these candidate genes and QTLs, we collected 73 Ujimqin sheep with increased numbers of vertebrae (T13L7, T14L6, and T14L7) and 23 sheep with normal numbers of vertebrae (T13L6). Through high-throughput genome resequencing, we obtained a total of 24,130,801 effective single-nucleotide polymorphisms (SNPs). By conducting a selective-sweep analysis, we discovered that the most significantly selective region was located on chromosome 7. Within this region, we identified several genes, including VRTN, SYNDIG1L, LTBP2, and ABCD4, known to regulate the spinal development and morphology. Further, a genome-wide association study (GWAS) performed on sheep with increased and normal vertebral numbers confirmed that ABCD4 is a candidate gene for determining the number of vertebrae in sheep. Additionally, the most significant SNP on chromosome 7 was identified as a candidate QTL. Moreover, we detected two missense mutations in the ABCD4 gene; one of these mutations (Chr7: 89393414, C > T) at position 22 leads to the conversion of arginine (Arg) to glutamine (Gln), which is expected to negatively affect the protein's function. Notably, a transcriptome expression profile in mouse embryonic development revealed that ABCD4 is highly expressed during the critical period of vertebral formation (4.5-7.5 days). Our study highlights ABCD4 as a potential major gene influencing the number of vertebrae in Ujimqin sheep, with promising prospects for future genome-assisted breeding improvements in sheep.

Effects of novel variants in BMP15 gene on litter size in Mongolia and Ujimqin sheep breeds

Reproductive traits, such as ovulation rate and litter size, are important factors influencing the sheep industry. The bone morphogenetic protein 15 (BMP15) is a major gene affecting the reproductive traits in sheep, and multiple mutations in BMP15 gene could affect the ovulation rate and litter size in many sheep breeds, showing high breed specificity. However, identification of novel variations and seeking breed-specific markers associated with litter size in other sheep breeds are still important. In this study, we sequenced the BMP15 gene of Mongolia sheep, and 12 novel variants were detected by direct sequencing and whole-genome resequencing. Among them, the g.50985975 G > A polymorphism in intron and synonymous c.755 T > C (Leu252Pro) in exon 2 of BMP15 were significantly associated with the litter sizes of Mongolia ewes (P < 0.01 and P < 0.01, respectively), as well as the g.50988478C > A and g.50987863G > A in the promoter region of BMP15 were significantly associated with the litter sizes of Ujimqin ewes (P < 0.05 and P < 0.01, respectively). The c.755 T > C mutation is predicted to change the tertiary structure of BMP15. Our findings may provide potentially useful genetic markers for increasing litter size in sheep.

Association between novel variants in BMPR1B gene and litter size in Mongolia and Ujimqin sheep breeds

Prolificacy is an important trait of animals, specifically for sheep. The Bone morphogenetic protein receptor 1B (BMPR1B) is a major gene affecting the litter size of many sheep breeds. The well-known FecB mutation (Q249R) was associated fully with the hyper prolific phenotype of Booroola Merino. However, the identification of variation in all exonic regions of BMPR1B was rare. In this study, we sequenced all exonic regions of BMPR1B gene of Mongolia sheep breed, and ten novel variants were detected by direct sequencing. Among them, the litter size of the Mongolia ewes with the CC genotype was significantly higher (0.34 additional lambs, p < .05) than those with the TT genotype of the g.29346567C>T single nucleotide polymorphism (SNP). The litter size of the Mongolia ewes with the TT genotype was significantly higher (0.19 additional lambs, p < .05 and .31 additional lambs, p < .01, respectively) than those with the GT and GG genotypes of the c.1470G>T SNP. The silent c.1470G>T mutation is predicted to increase the stability of the mRNA secondary structure through reducing minimum free energy and is predicted to change the mRNA secondary structure of BMPR1B. Our findings may give potentially useful genetic markers for increasing litter size in sheep.