A repertoire of single nucleotide polymorphisms (SNPs) of major fecundity BMPR1B gene among 75 sheep breeds worldwide

The BMPR1B gene is a major determinant of sheep reproductive capacity. Previous studies revealed that Q249R (FecB) is a profound variant of BMPR1B that influences the ovulation rate and litter size in sheep. However, unlike Q249R locus, the full spectrum of single nucleotide polymorphisms (SNPs) within BMPR1B has not been extensively studied. A systematic screen of SNPs in BMPR1B would facilitate the discovery of novel variants that are associated with litter size. This study aimed to investigate SNPs in the BMPR1B gene via whole genome sequence (WGS) data from 2409 individuals of 75 sheep breeds worldwide. Herein, a total of 9688 variants were screened, among which 15 were coding variants and 8 were novel changes. Specifically, we presented the most comprehensive frequency distribution map of the well-known FecB mutation to date. Besides, among the abovementioned SNPs, one synonymous mutation (g.30050773C > T) was found to be likely under selection and is potentially associated with fecundity in Duolang sheep. Thus, our study greatly expands the variation repertoire of the ovine BMPR1B gene and provides a valuable resource for exploring causative mutations and genetic markers associated with litter size.