Sheep litter size heredity basis using genome-wide selective analysis

Sheep are important herbivorous domestic animal globally, and the Chinese indigenous sheep breed has a multitude of economically significant variations due to the diverse geographical and ecological conditions. In particular, certain native breeds exhibit a visible high litter size phenotype due to the selection pressure of natural and artificial for thousands of years, offering an ideal animal model for investigating sheep's fecundity. In this study, selective signal analysis was performed on public whole-genome sequencing data from 60 sheep across eight breeds to identify candidate genes related to litter size. Results revealed that a total of 34,065,017 single-nucleotide polymorphisms (SNPs) were identified from all sheep, and 65 candidate genes (CDGs) were pinpointed from the top 1% of interacted windows and SNPs between the pairwise fixation index (F_ST, >0.149543) and cross-population extended haplotype homozygosity (XP-EHH, >0.701551). A total of 41 CDGs (e.g. VRTN, EYA2 and MCPH1) were annotated to 576 GO terms, of which seven terms were directly linked to follicular and embryonic development (e.g. TBXT, BMPR1B, and BMP2). In addition, 73 KEGG pathways were enriched by 21 CDGs (e.g. ENTPD5, ABCD4 and RXFP2), mainly related to Hippo (TCF4, BMPR1B and BMP2), TGF-β (BMPR1B and BMP2), PI3K-Akt (ITGB4, IL4R and PPP2R5A) and Jak-STAT signalling pathways (IL20RA and IL4R). Notably, a series of CDGs was under strong selection in sheep with high litter size traits. These findings result could improve the comprehension of the genetic underpinnings of sheep litter size. Furthermore, it provides valuable CDGS for future molecular breeding.