The average twin lambing rate of Bayanbulak sheep is 2% to 3%. However, a flock of sheep with a close genetic relationship and an average of 2 to 3 lambs per birth has been found recently. To determine the major genes controlling the prolificacy of the flock in the present study, the flock was designated A while 100 normal Bayanbulak sheep were randomly selected to comprise the control flock B. Ligase detection reaction method was applied to detect and analyze the 10 mutational loci of the 3 candidate prolificacy genes including bone morphogenetic protein type I receptors, bone morphogenetic protein 15, and growth differentiation factor 9. The 10 mutational loci are as follows: FecB locus of the BMPR-IB gene; $FecX^I$, $FecX^B$, $FecX^L$, $FecX^H$, $FecX^G$, and $FecX^R$ of the BMP15 gene; and G1, G8, and FecTT of the GDF9 gene. Two mutations including BMPR-IB/FecB and GDF9/G1 were found in Bayanbulak sheep. Independence test results of the two flocks demonstrate that the FecB locus has a significant effect on the lambing number of Bayanbulak sheep. However, the mutation frequency of the G1 locus in GDF9 is very low. Independence test results demonstrate that the GDF9 locus does not have a significant impact on the lambing performance of Bayanbulak sheep. Among the 10 detected loci, BMPR-IB/FecB is the major gene that influences the high lambing rate of Bayanbulak sheep.