A Profile of Single Nucleotide Polymorphisms in Fecundity Genes Among Iranian Sheep Breeds by Using Polymerase Chain Reaction–Restriction Fragment Length Polymorphism (PCR-RFLP) Method

Document Type: Research Articles


1 Department of Genomics, Agricultural Biotechnology Research Institute of Iran (ABRII), North Branch, Rasht, Iran

2 Department of Animal Science, Faculty of Agricultural Science, University of Guilan, Rasht, Iran

3 Department of Animal Biotechnology, Agricultural Biotechnology Research Institute of Iran (ABRII), North Branch, Rasht, Iran


Ovulation rate and litter size are traits of fecundity and have a very remarkable economic value that leads to an increase in the production of meat, wool and milk. Therefore, exploring on the candidate and functional genes for these traits is very important. Growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and bone morphogenetic protein receptor-1B (BMPR1B) are the most important of these genes containing single nucleotide polymorphisms (SNPs) markers that are associated with high fertility in sheeps. The aim of this study was to determine the allelic and genotypic frequencies of seven of these SNPs including; G1, G4, G7 and G8 in GDF9 gene, B2 and B4 in BMP15 gene and FecB in BMPR1B gene among 532 sheeps from eleven of Iranian native breeds using polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP) method. For each mutation region, polymerase chain reactions and digestion reactions with specific endonuclease enzymes, were performed. The results showed that G1 mutation presents among ten breeds except Zandi, G4 mutation in all breeds, and the FecB mutation in two breeds; Kalakuei and Zandi. The G7, G8, B2 and B4 mutations were not observed in any of breeds, and all sheeps were homozygote wild type for these four positions.



The domesticated sheep breeds are probably originated from the wild sheep species such as Asian mouflons (Ovis orientalis) and Asiatic urial (Ovis vignei), which are still found in Iran National Parks and protected natural areas (Valizadeh, 2010). Iran has 28 native sheep breeds and in terms of breed diversity in sheep, has occupied the first rank in the world (‎Akbarinejad et al. 2014). According to the latest livestock census, conducted in 2017 by Ministry of Jihad-e-Agriculture of Iran (Agricultural Statistics of Iran, 2017), the Iranian sheep population is 47638600 animals that encompass 63.60% of Iranian animals (sheeps/lambs, goats/kids, cattles/calves, buffalos and camels). The meat production trait has taken into account in most previous Iranian sheep breeding program and fertility traits have played the main role to improve lamb and mutton production. Based on previous investigations, different mutations in the transforming growth factor-β (TGFβ) superfamily genes have been reported to affect ovulation rate and litter size in sheep. These genes include bone morphogenetic protein receptor, type 1B (BMPR-1B), bone morphogenetic protein 15 (BMP15), and growth differentiation factor 9 (GDF9) (‎Jansson, 2014). The Booroola's ancestor is Garole breed sheep that is indigenous to Bengal area in India (‎Abraham and Thomas, 2012). The Garole is a small sheep that have a mean litter size of 1.74 to 2.27 (Davis et al. 2002). The Booroola phenotype, which increased the litter size, was originally selected from special lines of Australian merino sheep. Booroola fecundity gene is the first major gene for prolificacy identified in sheep ‎(Wilson et al. 2001). This gene is located on chromosome 6 of sheep and is effective on the maturation of the granulosa cells, the development of the immature oocyte and its function (Abraham and Thomas, 2012). This gene is known as the activin receptor-like kinase (ALK6) and its mutation called FecB. A change from adenine (in wild type) to guanine (in mutant types) at the nucleotide position 746, results in a change from glutamine amino acid (in wild type) to arginine amino acid (in mutant types) at coding region 249 (Mulsant et al. 2001). This single autosomal locus can improve the ovulation rate and litter size. The FecB mutation, in addition to the Merino Booroola breed, also presents in some other sheep breeds such as; Hu ‎(Chu et al. 2011; Wang et al. 2018‎), Small Tailed Han (Chu et al. 2006; Zhang et al. 2011‎) and Bayanbulak in China (‎Zuo et al. 2013), Kendrapada in India (‎Dash et al. 2017) and etc and effects on high fertility in these breeds. Bone morphogenetic protein 15 (BMP15) (FecX) gene is located on the chromosome X of sheep and has two exons separated by an intron, which is about 5400 nucleotides. The protein encoded by this gene has 1179 nucleotides and its immature peptide has 393 amino acids. The mature peptide has 125 amino acids (Galloway et al. 2000). The gene, which was first found in Romney sheep, is also called GDF9B or FecX. This gene, and also GDF9, encode separate proteins that express in ovarian tissue. This mechanism is very effective in mammalian fertilization (‎Hanrahan et al. 2004). Until now, several mutations of this gene including Inverdale (FecXI), Hanna (FecXH), Belclare (FecXB), Galway (FecXG), Lacaune (FecXL), Rasa aragonesa (FecXR), Grivette (FecXGr) and Olkuska (FecXO) have been identified in different sheep breeds. The first six mutations increase the rate of ovulation in heterozygous individuals (‎Qiuyue et al. 2014). For instance, the estimated effect of FecXG is ranged from 0.77 in Belclare ewes to 1.18 in Cambridge ewes, and for FecXB is 2.38 in Belclare ewes (‎‎Hanrahan et al. 2004). Growth differentiation factor 9 (GDF9) (FecG) is essential for folliculogenesis in sheep ovary and has a highly significant role from the first stage of follicular growth until ovulation. GDF9 gene is located on chromosome 5 of sheep and has about 2500 nucleotides. Two exons are separated by an intron with 1126 nucleotides. The two exons encode an immature peptide with 453 amino acids. The mature peptide has 135 residues (‎‎Hanrahan et al. 2004). FecGH or G8 is one of eight mutations (G1 to G8) in this gene that has been identified in Cambridge and Belclare sheeps and led to an increase in ovulation rates in heterozygotes and infertility in mutant homozygotes ewes. The effect of that is 1.79 in Belclare ewes and 2.35 in Cambridge ewes. Up to now, the FecGH has been reported in breeds of sheep like as Chios and Karagouniki in Greece too (‎Liandris et al. 2012). So far, in addition to the FecGH, some other mutations in this gene such as; FecGV in Ile de France, FecGE in Embrapa and FecGT in Icelandic Thoka sheep breeds have also been identified that demonstrated an association with high ovulation rate and litter size (‎Qiuyue et al. 2014). The purpose of this study was to investigate the presence or absence of seven SNPs in the three mentioned genes and to determine their allelic and genotypic frequencies among these eleven breeds.



Blood samples were collected from 532 individuals belong to 11 Iranian native sheep breeds including; Afshari (N=55), Farahani (N=55), Grey Shiraz (N=55), Turki-Ghashghaei (N=55), Naeini (N=48), Kalakuei (N=47), Moghani (N=46), Shal (N=46), Taleshi (N=45), Makuei (N=41) and Zandi (N=39). The detailed information are given in Table 1. The detailed information of SNPs (Table 2) and the position of amino acid changes on amino acid sequences in genes (Figure 1) are given. The blood samples was collected between 2005 and 2007 according to breed distribution areas over a large geographic area. During sampling, attempts were made to collect animals as unrelative as possible. A modified protocol was applied to isolate genomic DNA from blood samples using GPP kit (Gene Pajoohan Pooya). The extracted DNA was qualified and quantified using NanoDrop (2000) and agarose gel electrophoresis methods. Subsequently, polymerase chain reactions were performed. Some of the primers were deliberately introduced a point mismatch (shown by underlines in Table 3) resulting in PCR products containing restriction site for each endonuclease enzyme (‎Davis et al. 2002; ‎‎‎Hanrahan et al. 2004). Details of each primer and PCR products (Table 3), PCR reactions (Table 4) and the temperature program for SNPs (Table 5) are given. The total volume of each reaction was considered 15 µL and the applied biosystems (AB) thermal cycler (model 9700) was used for amplification.


Table 1 Introduction of eleven breeds in this study



Table 2 Informations of the seven single nucleotide polymorphisms (SNPs)



Figure 1 a: four amino acid changes in GDF9; b: two amino acid changes in BMP15 and c: one amino acid change in BMPR1B



Table 3 details of primers


SNP: single nucleotide polymorphism and PCR:polymerase chain reaction.


Table 4 Polymerase chain reaction reaction components



Table 5 Polymerase chain reaction (PCR) program in thermal cycler


SNP: single nucleotide polymorphism.


Figure 2 indicates the sequences of seven PCR products. The sequences highlighted in red are annealing sites of the primers, the green nucleotide is the position of mutation and the underlined sequences are cutting sites of each enzyme. Restriction digestion process was carried out by using a specific restriction enzymes for each mutation. Information of each restriction enzyme and the size of fragments obtained from the products of the enzymes are given in Tables 6 and 7, respectively. The digested products were electrophoresed in 2% agarose gel and the DNA bands were visualized by ethidium bromide staining technique. Existence or absence of a restriction site is the source of polymorphism. For genotyping, each band represents an allele having a specific size in base pairs (bp). The data were subjected to appropriate analyses using POPGENE software (ver. 1.31).



A total of 532 individuals from the 11 Iranian sheep breeds were genotyped with the PCR-RFLP approach. The electrophoretic profiles of genomic DNA and PCR products of seven loci are shown in Figures 3 to 10. Digestion with restriction enzymes in seven loci on agarose geles (2%) have been shown in Figures 11 to 17. Subsequently, statistical analaysis were carried out. Frequencies of genotypes and allelesfor G1 mutation (G→ A), G4 mutation (G→A), FecB mutation (A→G) and other SNPs are given in Tables 8 to 11, respectively. Due to the fact that polymorphism was observed only in the G1, G4 and FecB loci, the Hardy-Weinberg equilibrium for these three positions was examined based on the chi-square test. The HWE test showed that all three polymorphic loci did not deviate from HWE (P>0.05) while considering each breed individually.


Figure 2 The PCR product sequences

a: 462 bp; b: 161 bp; c: 156 bp and d: 139 bp in GDF9; e: 141 bp and f: 153 bp in BMP15; g: 190 bp in BMPR1B



Table 6 Information of restriction enzymes



Table 7 Digestion reaction components



But when all the breeds were combined into one population, deviation from HWE happened only in FecB mutation (P<0.05). Table 12 showes the p-value of HWE test for each mutation. Genotyping of G1 showed that there was no recessive homozygote individuals in all breeds. The highest and lowest frequencies of heterozygotes (GA) obtained in Taleshi and Zandi breeds with 0.156 and 0.000, respectively. The mean frequency of heterozygotes in all sheeps as one population was 0.068. The highest and lowest frequency of mutant allele (A) was observed in Taleshi and Zandi breeds with 0.078 and 0.000, respectively. Mean of mutant allele frequency per all animals was 0.034. Further details about this mutation to be extracted from different references are given in Table 13. The G1 mutation was identified in Cambridge and Belclare sheeps for the first time but without any association with ovulation rate (‎Hanrahan et al. 2004). With looking at the Table 13, it can be seen that, in other studies such as ours, which aimed only at evaluating allelic and genotypic frequencies, no sample of mutant homozygote is reported. It indicates that the frequency of this genotype is very low in the main populations of these breeds. In other studies with sheep fecundity records (and probably, the samples have also been selected based on their records), the number of mutant homozygous is low.


Figure 3 Some of genomic DNA in this study



Figure 4 G1 locus; a 462 bp band



Figure 5 G4 locus; a 161 bp band



Figure 6 G7 locus; a 156 bp band



The highest frequency of mutant homozygotes has been reported in the Greek Chios breed with 0.087 (‎Liandris et al. 2012). The homozygous (AA) ewes in this breed displayed litter size of 2.25 implying a significant effect of the A allele, as high as 0.33 lambs. The heterozygous ewes showed negative dominance deviation with a statistically significant decrease in litter size of 0.47 lambs.


Figure 7 G8 locus; a 139 bp band



Figure 8 B2 locus; a 141 bp band



Figure 9 B4 locus; a 153 bp band



Figure 10 FecB locus; a 190 bp band



The G1 has a major effect on litter size in Baluchi sheep breed of Iran. The heterozygous and wild homozygous ewes had 0.35 and 0.21 more lambs than the homozygous mutant ewes, respectively (‎Moradband et al. 2011). Another report says this mutation leads to increased litter size in heterozygotes (1.56) compared with the wild homozygotes (1.25) in Chilota sheeps in Chile (‎Paz et al. 2015). However, some other studies did not report any association between this mutation and high fertility in sheep. Genotyping of G4 revealed that there was no recessive homozygote in all studied breeds.


Figure 11 Polymorphism in G1 locus



Figure 12 Polymorphism in G4 locus



Figure 13 Monomorphism in G7 locus



Figure 14 Monomorphism in G8 locus



The highest and lowest frequency of heterozygotes (GA) obtained in Taleshi and Shal breeds with 0.156 and 0.043, respectively. The mean frequency of heterozygotes in all animals was 0.115.


Figure 15 Monomorphism in B2 locus



Figure 16 Monomorphism in B4 locus



Figure 17 Polymorphism in FecB locus



The highest and lowest frequencies of the mutant allele (A) were observed in Taleshi and Shal breeds with 0.078 and 0.022 respectively. Mean frequency of the mutant allele in all animals was 0.057. This mutation has not been reported in Garole sheep of India (‎Polley et al. 2010) but has been identified in some other breeds around the world. More details are given in Table 14. For the first time, it had been reported that 13 of 29 Belclare sheep, and 1 of 26 Cambridge sheep, had at least one copy of the G4 mutation with no association to fertility (‎‎Hanrahan et al. 2004). In Greek Chios sheep breed, the highest litter size (2.17) was observed in heterozygotes and a significant positive dominance deviation as high as 0.61 lambs was obtained (Liandris et al. 2012).


Table 8 Frequencies of genotypes and alleles for G1 mutation



Table 9 Frequencies of genotypes and alleles for G4 mutation



Table 10 Frequencies of genotypes and alleles for FecB mutation



Likewise, another report from India, indicates that the average litter size was 1.63, 2.00 and 1.91 for wild homozygous, heterozygous and mutant homozygous of Kendrapada ewes, respectively (‎Dash et al. 2017). Genotyping of FecB revealed that this mutation only finds in Kalakuei and Zandi breeds. The frequencies of wild homozygote (AA), heterozygote (AG) and mutant homozygote (GG) obtained 0.809, 0.170 and 0.021 respectively, for Kalakuei breed and 0.949, 0.051 and 0.000 for Zandi breed. The frequencies of the wild and mutant alleles were 0.894 and 0.106 respectively, in Kalakuei and 0.974 and 0.026 in Zandi. Mean frequency of the mutant allele in all samples was calculated as 0.011. More details about presence and absence of this mutation in many sheep breeds has been brought in Table 15. As seen in the Table 15, all reports of the occurrence of this mutation in sheeps show a significant relationship between that and higher fertility traits, except for the Bonapala breed of India.


Table 11 Frequencies of genotypes and alleles for G7, G8, B2 and B4 mutations



Table 12 Hardy-Weinberg equilibrium in three polymorphic sites1


* (P<0.05).

NS: non significant.


In this sheep, although the litter size in mutant homozygotes (1.71) was more than heterozygotes (1.64) and both were more than wild homozygotes (1.50), but the statistical analysis didn’t reveal a significant association (‎Roy et al. 2011). FecB mutation showed an additive effect on litter size in Hu sheep in China. Mutant homozygote and heterozygote ewes significantly have 0.713 and 0.419 more lambs than non carrier homozygotes, respectively (‎Wang et al. 2018). In an investigation, this mutation showed an additive effect on ovulation rate and partially dominant effect on litter size in Small Tailed Han sheep of China. In this study, mean of litter size in homozygote mutant and heterozygote ewes was 2.65 and 2.16 respectively, which were significantly more than 1.14 of wild homozygotes (‎Chu et al. 2011). Another report from China clarified a significant effect on litter size in Small Tailed Han, Poll Dorset and their F1 generation, in total (Jia et al. 2005). The ewes with BB and B+ genotypes had about 1.04 and 0.74 lambs more than wild non-carriers ewes. Also, the mean of litter size in BB ewes was significantly more than ++ ewes in both breeds. FecB mutation also showed a phenotypic effect on litter size in Indian Kendrapada sheep breed. The average litter size of non-carriers, heterozygotes and homozygote carriers were 1.61, 1.80 and 2.06 respectively (‎Dash et al. 2017). Another report suggests a highly significant difference in the average of litter size in different genotypes of Indonesian fat Tailed sheep (‎Maskur et al. 2016). Mutant homozygous ewes had the highest average (1.685) compared with heterozygotes (1.455) and wild types (1.145). In Iranian Kalehkoohi sheep, genotypes BB (1.902) and B+ (1.725) had 0.52 and 0.35 lambs, more than ++ type (1.379), respectively and that was a statistically significant difference (P<0.01). But no significant difference between BB and B+ genotypes were observed (‎Mahdavi et al. 2014). Reffering to the Table 15, it can be seen that all breeds that carry the FecB allele, are distributed in different parts of Asia (China, Indonesia, India and Iran). Most likely, the Merino Booroola breed in Australia also received this gene from the Indian Garole sheep in the 18th century (‎‎Davis et al. 2002). The effect of FecB is additive for ovulation rate and can increase by about 1.6 corpora lutea per cycle for each copy in Booroola sheep (‎Wilson et al. 2001).


Table 13 The polymorphisms of G1 reported in some studies1


1 This table is arranged based on the frequency of the mutant allele, from the highest to the lowest.

SAWF: significant association with fecundity.


BMPRIB gene would be inactivated partially, leading to an advanced differentiation of granulosa cells and also an advanced maturation of ovulatory follicles in homozygous mutant ewes of Australian Booroola merino (‎Mulsant et al. 2001). Genotyping of G7 revealed that this SNP didn’t occur in any of these eleven breeds and all sheeps have dominant homozygote genotypes. Our report about the absence of this mutation is consistent with several reports of other sheep breeds and of course, is not compatible with a few reports of other breeds. More details about this mutation in several sheep breeds are given in Table 16. For the first time, It had been reported that at least one copy of G7 mutation presents in 2 of 19 Belclare and 7 of 24 Cambridge sheeps with no linkage to high fertility ‎(‎‎Hanrahan et al. 2004). A similar study reported that there was no statistical significant association between mutation and ovulation rate in Finnish Landrace sheep, although a significant association was observed in the Belclare sheep. Heterozygote ewes had higher ovulation than wild homozygotes (0.17) and ovulation rates in heterozygotes were significantly lower than mean of homozygotes (‎Mullen and Hanrahan et al. 2014).


Table 14 The polymorphisms of G4 reported in some studies1


1 This table is arranged based on the frequency of the mutant allele, from the highest to the lowest.

SAWF: significant association with fecundity.


Based on the estimated breeding values for the Norwegian White sheeps, daughters of homozygous mutant rams produce additional lambs (0.46 to 0.57) compared to daughters of wild homozygous rams. Also, daughters of heterozygous rams produce more lambs (0.20 to 0.25) than daughters of homozygous wild rams (Vage et al. 2013). Also mentioned that the average litter size for ewes at first year of age is about 2.061 and at second year is about 2.671, showing that the average litter size is higher at the second year and the genotype showed a significant effect on litter size at both years (‎Siddiqua, 2013). In Olkuska sheep, only the presence of one copy of mutant allele, resulted in an increase (0.55 lamb) in litter size. Heterozygous ewes demonstrated high annual lamb production of 346%, with the same for dominant homozygous ewes at 236% . It is clear from Table 16, that those breeds which have this mutation are natives of the northern Europe or raised in that region. Genotyping of G8 (FecGH) revealed that this mutation didn’t occur in any of the studied breeds and all sheeps were dominant homozygotes. Our report about the absence of this mutation is not compatible with some reports of other sheep breeds but is consistent with many reports of other sheep breeds. More details are given in Table 17. It has been identified that this mutation is essential for folliculogenesis and results in high fertility in heterozygotes and infertility in mutant homozygotes in Cambridge and Belclare sheeps. Estimation of the effect of the mutation on the ovulation rate showed that Belclare ewes had an increase of about 1.79 and Cambridge ewes had about 2.35. Based on progeny test data in Belclare sheeps, it has been shown that G8 can increase the ovulation rate of about 0.83 in ewe lambs and 1.75 in adult ewes, which both were significant. Overall, the mutation has led to an increase in ovulation rate about 1.39 in both breeds ‎(‎‎Hanrahan et al. 2004). Also, another study on Chios breed in Greece reported a statistically positive dominance deviation (0.44 lambs) for this mutation in heterozygous ewes but with no evidence of infertility in homozygous mutant ewes (‎Liandris et al. 2012). It has been suggested that the Lleyn sheep breed is the most likely source of this SNP (and also FecXG) in Belclare and Cambridge sheeps (‎Mullen et al. 2013). Genotyping of B2 (FecXG) revealed that this mutation did not occur in any of eleven breeds and all animals were dominant homozygotes. More details are given in Table 18. In Cambridge and Belclare breeds, the homozygous ewes for this mutation and also, heterozygous ewes for this mutation and B4 simultaneously, were sterile. The mutation significantly increases ovulation rate about 0.77 in Belclare and 1.18 in Cambridge ewes. Based on progeny test data in Belclare sheep, it has been shown that the mutation increases the ovulation rate of about 0.62 in ewe lambs and 0.72 in adult ewes, which both were significant. Overall, this mutation leads to an increase in ovulation rate of about 0.70 in both breeds. The effects of this mutation with B4 was about 0.28 in the Belclare breed ‎(‎‎Hanrahan et al. 2004). A study on Small Tailed Han sheep showed that the polymorphism of this mutation was not significantly associated with litter size but tended to be associated with that trait. In Hu sheep, the heterozygous ewes had 0.52 more lambs than the wild genotype ones. In total population, the heterozygous Hu and Small Tailed Han ewes showed 0.30 more lambs than those with wild genotype (‎Wang et al. 2015). Also, based on a study, the heterozygous Small Tailed Han ewes had 0.55 more lambs than wild homozygotes in China (‎Chu et al. 2006). However, no significant allele effects were obtained in this mutation polymorphisms in Chios and Karagouniki breeds of Greece )‎Liandris et al. 2012(. It has been suggested that the Lleyn sheep breed is the most likely source of this mutation (and FecGH) in Belclare and Cambridge sheeps (Mullen et al. 2013). Genotyping of B4 (FecXB) revealed that this SNP was absent in all of the eleven breeds and all sheeps were dominant homozygotes. More details of this mutation in some studies are given in Table 19. This SNP has been identified in Belclare sheepbut not in Cambridge sheep.


Table 15 The polymorphisms of FecB reported in some studies1


Continued Table 15


1 This table is arranged based on the frequency of the mutant allele, from the highest to the lowest.

SAWF: significant association with fecundity.


Table 16 The polymorphisms of G7 reported in some studies1


1 This table is arranged based on the frequency of the mutant allele, from the highest to the lowest.

SAWF: significant association with fecundity.


The mutation significantly increases ovulation rate and its effect was estimated about 2.38. Homozygous mutant ewes and also heterozygous ewes for this mutation and B2 simultaneously, were sterile. Based on progeny test data in Belclare sheep, this SNP increases the ovulation rate of about 0.76 in ewe lambs and 1.11 in adult ewes, which both were significant. Overall, the mutation leads to an increase in ovulation rate of about 0.97 in both breeds. The effect of this mutation with B2 was about 0.28 in the Belclare breed ‎(‎‎Hanrahan et al. 2004). Based on some finding, the FecXB mutation that found among the hyper-prolific ewes in Ireland, has come from a high fertility line that had developed with the use of prolific ewes that choosed from commercial flocks in that country and then used for the genesis of the Belclare sheep breed (‎Mullen et al. 2013). In Greek Chios and Karagouniki breeds, no significant allele effect weas obtained in the polymorphisms analysis of this mutation (‎Liandris et al. 2012). Survey of linkage between G1 and G4 indicates that all sheeps that were heterozygous by G1, were also heterozygous by G4. On the contrary, sheeps that were heterozygous for G4 did not necessarily had the G1. Hence, the frequency of the G4 was more than the G1. It seems that there is a local haplotype for these SNPs in most breeds. As it is shown in Table 20, Taleshi, Shal and Naeini breeds have the highest linkage and the Zandi breed shows the lowest linkage in the nucleotide position of two mutations. Considering all animals as a population, three local haplotypes were identified. A similar report indicates a high degree of linkage disequilibrium between these SNPs in Mehraban sheep breed that has been confirmed by the high correlation coefficient values (r2=0.91) (‎Talebi et al. 2018). Based on the observed polymorphism in the G1, G4 and FecB mutations of GDF9 and BMPR1B genes, all animals used in the present study were classified in six genotypic groups (Table 21).


Table 17 The polymorphisms of G8 (FecGH) reported in some studies1


1 This table is arranged based on the frequency of the mutant allele, from the highest to the lowest.

SAWF: significant association with fecundity.


Table 18 The polymorphisms of B2 (FecXG) reported in some studies1


1 This table is arranged based on the frequency of the mutant allele, from the highest to the lowest.

SAWF: significant association with fecundity.


Table 19 The polymorphisms of B4 (FecXB) reported in some studies1


1 This table is arranged based on the frequency of the mutant allele, from the highest to the lowest.

SAWF: significant association with fecundity.


Table 20 G1 and G4 mutations and three observed local haplotypes



Table 21 Genotypic classification of Iranian sheep breeds based on three point mutations at GDF9 and BMPR1B genes




Fertility is largely influenced by environmental conditions. Having larger populations, minimizing environmental factors, accurate, multiple and seasonal recordings and repeatability estimation are factors that can better determine the association or disassociation between mutants and the reproductive potential of sheeps. In the present study, the presence or absence of seven reported SNPs known as G1, G4, G7 and G8 in GDF9 gene, B2 and B4 in BMP15 gene and FecB in BMPR1B gene were examined among eleven Iranian native sheep breeds that are different in litter size. The results showed that G7, G8, B2, and B4 mutations were not in any of these breeds, and all of the sheeps were monomorphic with dominant homozygote genotype. FecB in the two breeds of Kalakuei and Zandi, G1 in all breeds except the Zandi and G4 mutation in all eleven breeds were identified. The absence of reported major mutations (such as G8, B2, B4, and FecB), suggests that minor mutations of fecundity genes and their cumulative effects may have pivotal role in monitoring of fertility traits among Iranian sheeps. It is suggested that to achieve a deep insight on how to control the fertility characteristics of the native sheep in Iran, exploring the major genes controlling fecundity trait in Iranian indigenous sheep breeds is performed through whole genome sequencing and SNP calling approaches.



The authors are grateful to the local farmers of the target provinces to provide the herd for sampling.

Abraham A. and Thomas N. (2012). Role of fecundity genes in prolificacy of small ruminants. J. Ind. Vet. Assoc. Kerala. 10(3), 34-37.

Agricultural Statistics of Iran. (2017). Ministry of Jihad-e-Agriculture, Deputy Director of Planning and Economics. Center for Information and Communication Technology. Tehran, Iran.

Ahmadi A., Afraz F., Talebi R., Farahavar A. and Vahidi S.M.F. (2016). Investigation of GDF9 and BMP15 polymorphisms in Mehraban sheep to find the missenses as impact on protein. Iranian J. Appl. Anim. Sci. 6(4), 863-872.

Akbarinejad V., Kazempooor R., Shojaei M., Rezaee M. and Oji R. (2014). Atlas of Iranian Sheep Breeds. Noorbakhsh Press upon the Request of Biotechnology Development Council. Tehran, Iran.

Akbarpour M., Houshmand M., Ghorashi A. and Hayatgheibi H. (2008). Screening for FecGH mutation of growth differentiation factor 9 gene in Iranian Ghezel sheep population. Royan Inst. Int. J. Fertil. Steril. 2(3), 139-144.

Al-Barzinji Y.M. and Othman G.U. (2013). Genetic polymorphism in FecB gene in Iraqi sheep breeds using RFLP-PCR technique. IOSR J. Agric. Vet. Sci. 2(4), 46-48.

Ali A.S., Ibrahim M.T., Mohammed M.M., Elobied A.A. and Lühken G. (2016). Growth differentiation factor 9 gene variants in Sudanese desert sheep ecotypes. South African J. Anim. Sci. 46(4), 373-379.

Bahrami Y., Bahrami S., Mohammadi H.R., Chekani-Azar V. and Mousavizadeh S.A. (2014). The polymorphism of GDF-9 gene in Hisari sheep. Biol. Forum. 6(2), 46-52.

Chu M., Jia L., Zhang Y., Jin M., Chen H., Fang L., Di R., Cao G., Feng T., Tang Q., Ma Y. and Li K. (2011). Polymorphisms of coding region of BMPR-IB gene and their relationship with litter size in sheep. Mol. Biol. Rep38(6), 4071-4076.

Chu M.X., Liu Z.H., Jiao C.L., He Y.Q., Fang L., Ye S.C., Chen G.H. and Wang J.Y. (2006). Mutations in BMPR-IB and BMP-15 genes are associated with litter size in Small Tailed Han sheep (Ovis aries). J. Anim. Scence. 85, 598-603.

Dash S., Maity A., Bisoi P.C., Palai T.K., Polley S., Mukherjee A. and De S. (2017). Coexistance of polymorphism in fecundity genes BMPR1B and GDF9 of Indian Kendrapada sheep. Exp. Anim. Med. Res7(1), 33-38.

Davis G.H., Farquhar P.A., O'Connell A.R., Everett-Hincks J.M., Wishart P.J., Galloway S.M. and Dodds K.G. (2006). A putative autosomal gene increasing ovulation rate in Romney sheep. Anim. Reprod. Sci. 92(1), 65-73.

Davis G.H., Galloway S.M., Ross K.I., Gergan S.M., Ward J., Nimbkar B.V., Ghalsasi P.M., Nimbkar C., Gray G.D., Inounu I., Tiesnamurti B., Hanrahan J.P. and Wilson T. (2002). DNA tests in prolific sheep from eight countries provide new evidence on origin of the Booroola (FecB) mutation. Biol. Reprod. 66, 1869-1874.

Debnath J. and Singh R.V. (2014). Genetic polymorphism of boorola FecB gene and its association with litter size in Balangir, Shahabadi and Bonpala sheep breeds. Indian J. Anim. Res48(4), 307-314.

Dincel D., Ardicli S., Samli H. and Balci F. (2018). Genotype frequency of FecXB (Belclare) mutation of BMP15 gene in Chios (Sakiz) sheep. Uludag University J. Fac. Vet. Med37(2), 1-5.

Eghbalsaied S., Amini H., Shahmoradi S. and Farahi M. (2014). Simultaneous presence of G1 and G4 mutations in growth differentiation factor 9 gene of Iranian sheep. Iranian J. Appl. Anim. Sci. 4(4), 781-785.

El-Hanafy A.A. and El-Saadani M.A. (2009). Fingerprinting of FecB gene in five Egyptian sheep breeds. Biotechnol. Anim. Husb25(3), 205-212.

Escobar-Chaparro R.A., Guillén G., Espejo-Galicia L.U., Meza-Villalvazo V.M., Pena-Castro J.M. and Abad-Zavaleta J. (2017). qPCR and HRM-based diagnosis of SNPs on growth differentiation factor 9 (GDF9), a gene associated with sheep (Ovis aries) prolificacy. 3 Biotech7(3), 204-211.

Esen Gursel F., Akis I., Durak H., Mengi A. and Oztabak K. (2011). Determination of BMP-15, BMPR-1B and GDF-9 gene mutations of the indigenous sheep breeds in Turkey. Kafkas Univ. Vet. Fak. Derg. 17(5), 725-729.

Galloway S.M., McNatty K.P., Cambridge L.M., Laitinen M.P., Juengel J.L., Jokiranta T.S., McLaren R.J., Luiro K., Dodds K.G., Montgomery G.W., Beattie A.E., Davis G.H. and Ritvos O. (2000). Mutations in an oocyte-derived growth factor gene (BMP15) cause increased ovulation rate and infertility in a dosage-sensitive manner. Nat. Genet. 25(3), 279-283.

Ghaderi A., Beigi nasiri M.T., Mirzadeh K.H., Fayazi J. and Sadr A.S. (2010). Identification of the GDF9 mutation in two sheep breeds by using polymerase chain reaction- restriction fragment length polymorphism (PCR-RFLP) technique. African J. Biotechnol. 9(47), 8020-8022.

Hanrahan J.P., Gregan S.M., Mulsant P., Mullen M., Davis G.H., Powell R. and Galloway S.M. (2004). Mutations in the genes for oocyte-derived growth factors GDF9 and BMP15 are associated with both increased ovulation rate and sterility in Cambridge and Belclare sheep (Ovis aries). Biol. Reprod. 70, 900-909.

Holanda G.M.L., Oliveira J.C., Silva D.M.F., Rocha S.S.N., Pandolfi V., Adriao M. and Wischral A. (2017). Survey of mutations in prolificacy genes in Santa Ines and Morada Nova sheep. Arq. Bras. Med. Vet. Zootec. 69(4), 1047-1053.

Jamshidi R., Kasirian M.M. and Rahimi G.A. (2013). Application of PCR-RFLP technique to determine Booroola gene polymorphism in the Sangsari sheep breed of Iran. Turkish J. Vet. Anim. Sci37(2), 129-133.

Jansson T. (2014). Genes involved in ovulation rate and litter size in sheep. BScThesis. Examensarbete, Swedish University of Agricultural Science, Uppsala, Sweden.

Jawasreh K.I.Z., Awawdeh F.T. and Al-Qaisy A. (2014). Association between GDF9, FecB and prolactin gene polymorphisms and prolificacy of Awassi sheep. Pp. 25-32 in 10th World Congr. Genet. Appl. Livest. Prod. Vancouver, Canada.

Jia C., Li N., Zhao X., Zhu X. and Jia Z. (2005). Association of single nucleotide polymorphisms in exon 6 region of BMPRIB gene with litter size traits in sheep. Asian Australian J. Anim. Sci18(10), 1375-1378.

Karsli T., Sahin E., Karsli B.A., Alkan S. and Balcioglu M.S. (2012). An investigation of mutations (FecXG, FecXI, FecXH, FecXB) on BMP-15 gene in some local sheep breeds raised in Turkey. Akdeniz Üniv. Ziraat Fak. 25(1), 29-33.

Khaldari M. (2014). Sheep and Goat Husbandry. Publishing Organization, Tehran, Iran.

Khanahmadi A., Gharehveysi S., Khatamnejad R. and Arab J. (2014). Polymorphic variants of G1 and B4 from GDF9 and BMP15 genes of Dalagh sheep. Res. Anim. Prod. 5(10), 148-155.

Khatamnejhad R. and Khanahmadi A. (2012). Genetic polymorphism in GDF9 and FecB genes in Dalagh sheep breed of Iran. J. Anim. Vet. Adv. 11(6), 766-768.

Khodabakhshzadeh R., Mohammadabadi M.R., Esmailizadeh A.K., Shahrebabak H.M., Bordbar F. and Namin S.A. (2016). Identification of point mutations in exon 2 of GDF9 gene in Kermani sheep. Polish J. Vet. Sci19(2), 281-289.

Kumar R., Taraphder S., Sahoo A., Senpati P., Paul R., Roy M., Sarkar U., Raja D., Samanta I. and Baidya S. (2016). Polymorphism and nucleotide sequencing of BMP-15 gene in an organised herd of Garole sheep (Ovis aries). Indian J. Anim. Hlth55(1), 77-86.

Kumar S., Mishra A.K., Kolte A.P., Dash S.K. and Karim S.A. (2008). Screening for Booroola (FecB) and Galway (FecXG) mutations in Indian sheep. Small Rumin. Res80(1), 57-61.

Liandris E., Kominakis A., Andreadou M., Kapeoldassi K., Chadio S., Tsiligianni T., Gazouli M. and Ikonomopoulos I. (2012). Associations between single nucleotide polymorphisms of GDF9 and BMP15 genes and litter size in two dairy sheep breeds of Greece. Small Rumin. Res. 107, 16-21.

Lopez-Ramirez R.B., Magana-Sevilla H.F., Zamora-Bustillos R., Ramon-Ugalde J.P. and Gonzalez-Mendoza D. (2014). Analysis of the 3′ end regions of the GDF9 and BMPR1B genes in Blackbelly sheep from Yucatán, Mexico. Cien. Inv. Agr. 41(1), 123-128.

Mahdavi M., Nanekarani S. and Hosseini S.D. (2014). Mutation in BMPR-IB gene is associated with litter size in Iranian Kalehkoohi sheep. Anim. Reprod. Sci147(3), 93-98.

Mahrous K.F., Hassanane M.S., Abdel Mordy M., Shafey H.I. and Rushdi H.E. (2015). Polymorphism of some genes associated with meat-related traits in Egyptian sheep breeds. Iranian J. Appl. Anim. Sci5(3), 655-663.

Maskur M., Tapaul R. and Kasip L. (2016). Genetic polymorphism of bone morphogenetic protein receptor 1B (BMPR-1B) gene and its association with litter size in Indonesian fat-Tailed sheep. African J. Biotechnol. 15(25), 1315-1319.

Moradband F., Rahimi G. and Gholizadeh M. (2011). Association of polymorphisms in fecundity genes of GDF9, BMP15 and BMP15-1B with litter size in Iranian Baluchi sheep. Asian-Australasian J. Anim. Sci. 24(9), 1179-1183.

Mullen M.P. and Hanrahan J.P. (2014). Direct evidence on the contribution of a missense mutation in GDF9 to variation in ovulation rate of Finnsheep. PloS One. 9(4), e95251.

Mullen M.P., Hanrahan J.P., Howard D.J. and Powell R. (2013). Investigation of prolific sheep from UK and Ireland for evidence on origin of the mutations in BMP15 (FecXG, FecXB) and GDF9 (FecGH) in Belclare and Cambridge sheep. PloS One8(1), e53172.

Mulsant P., Lecerf F., Fabre S., Schibler L., Monget P., Lanneluc I., Pisselet C., Riquet J., Monniaux D., Callebaut I., Cribiu E., Thimonier J., Teyssier J., Bodin L., Cognié Y., Chitour N. and Elsen J.M. (2001). Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes. Proc. Natl. Acad. Sci98(9), 5104-5109.

Nagdy H., Mahmoud K.G.M., Kandiel M.M., Helmy N.A., Ibrahim S.S., Nawito M.F. and Othman O.E. (2018). PCR-RFLP of bone morphogenetic protein 15 (BMP15/FecX) gene as a candidate for prolificacy in sheep. Int. J. Vet. Sci. Med. 6, 68-72.

Nanekarani S., Goodarzi M., Khederzadeh S., Torabi S. and Landy N. (2016). Detection of polymorphism in Booroola gene and growth differentiation factor 9 in Lori sheep breed. Trop. J. Pharm. Res15(8), 1605-1611.

Oraon T., Singh D.K., Ghosh M., Kullu S.S., Kumar R. and Singh L.B. (2016). Allelic and genotypic frequencies in polymorphic Booroola fecundity gene and their association with multiple birth and postnatal growth in Chhotanagpuri sheep. Vet. World9(11), 1294-1299.

Pan Z., Zhang B., Hu W., Di R., Lin Q., Wang X., Yin D., Wang P. and Chu M. (2015). The genetic diversity of both FecB gene and microsatellite GC101 is associated with reproduction selection in sheep. Turkish J. Vet. Anim. Sci39(3), 254-260.

Paz E., Quinones J., Bravo S., Montaldo H.H. and Sepulveda N. (2015). Genotyping of BMPR1B, BMP15 and GDF9 genes in Chilean sheep breeds and association with prolificacy. Anim. Genet. 46(1), 98-109.

Paz E., Quinones J., Bravo S., Rodero E., Gonzalez A. and Sepulveda N. (2014). Identification of G1 and G8 polymorphisms of GDF9 gene in Araucano creole sheep. Arch. Med. Vet. 46, 327-331.

Pirzadi Y., Nanakarani S. and Godarzi M. (2014). Genetic polymorphism of GDF9 gene in Lori-Romanov crossbreeds sheep by PCR-RFLP. Pp. 45-51 in 1st Int. Conf. New Ideas. Agric. Islamic Azad University Khorasgan Branch, Isfahan, Iran.

Polley S., De S., Brahma B., Mukherjee A., Vinesh P.V., Batabyal S., Arora J.S., Pan S., Samanta A.K., Datta T.K. and Goswami S.L. (2010). Polymorphism of BMPR1B, BMP15 and GDF9 fecundity genes in prolific Garole sheep. Trop. Anim. Health Prod. 42, 985-993.

Qiuyue L., Zhangyuan P., Xiangyu W., Wenping H., Ran D., Yaxing Y. and Mingxing C. (2014). Progress on major genes for high fecundity in ewes. Front. Agric. Sci. Engin. 1(4), 282-290.

Roy J., Polley S., De S., Mukherjee A., Batabyal S., Pan S., Brahma B., Datta T.K. and Goswami S.L. (2011). Polymorphism of fecundity genes (FecB, FecX and FecG) in the Indian Bonpala sheep. Anim. Biotechnol. 22(3), 151-162.

Saadatnouri M. and Siahmansour S. (2011). Principles of Sheep Husbandry. Ashrafi Publications, Tehran, Iran.

Shafieiyan Z., Mohammadi G., Jolodarzadeh A. Amiri S. (2013). No mutations of FecB and FecGH in Iranian Lory sheep. Vet. Res. Forum. 4(4), 265-268.

Shi H., Bai J., Niu Z., Muniresha Fen L. and Jia B. (2010). Study on candidate gene for fecundity traits in xingjiang cele black sheep. African J. Biotechnol. 9(49), 8498-8505.

Siddiqua S.A. (2013). Phenotypic effects of a fertility mutation in Norwgain White sheep. MS Thesis. Norwegian University of Life Science, Ås, Norway.

Somarny W.W., Erin A.R., Suhaimi A.H.M.S., Nurulhuda M.O. and Hifzan R.M. (2013). A study of major prolificacy genes in Malin and Dorper sheep in Malaysia. J. Trop. Agric. Food Sci. 41(2), 265-272.

Souza C.J.H., McNeilly A.S., Benavides M.V., Melo E.O. and Moraes J.C.F. (2014). Mutation in the protease cleavage site of GDF 9 increases ovulation rate and litter size in heterozygous ewes and causes infertility in homozygous ewes. Anim. Genet45(5), 732-739.

Sudhakar A. (2009). Molecular polymorphism analysis of BMPR-IB, IGFBP-3 and POULF1 genes in Nilagiri and Mecheri sheep. Ph D. Thesis. Tamil Nadu Veterinary and Animal Sciences Univ., Chennai, India.

Talebi R., Ahmadi A., Afraz F., Sarry J., Woloszyn F. and Fabre S. (2018). Detection of single nucleotide polymorphisms at prolificacy major genes in the Mehraban sheep and association with litter size. Ann. Anim. Sci. 18(5), 685-698.

Tavakolian J. (2000). An introduction to Genetic Resources of Native Farm Animals in Iran. Animal Science Genetic Research Institute Press, Tehran, Iran.

Vacca G.M., Dhaouadi A., Rekik M., Carcangiu V., Pazzola M. and Dettori M.L. (2010). Prolificacy genotypes at BMPR 1B, BMP15 and GDF9 genes in North African sheep breeds. Small Rumin. Res. 88(1), 67-71.

Vage D.I., Husdal M., Matthew P.K., Klemetsdal G. and Boman I.A. (2013). A missense mutation in growth differentiation factor 9 (GDF9) is strongly associated with litter size in sheep. BMC Genet. 14, 1-10.

Valizadeh R. (2010). Iranian sheep and goat industry at a glance. Pp. 1-6 in Proc. Stress Manag. Small Rumin. Prod. Product Proc. Jaipur, India.

Vera J., Bravo S., Paz E. and Sepulveda Y.N. (2015). Detection of polymorphisms major gene GDF9 in corriedale flocks in Aysen region, Chile. Chilean J. Agric. Anim. Sci., ex Agro-Cien. 31(1), 15-20.

Vivekanand J.P. (2010). Screening of local sheep breeds of gujarat for presence of fecundity gene mutations by PCR-RFLP. Ph D. Thesis. Anand Agricultural Univ., Gujarat, India.

Wang W., La Y., Zhou X., Zhang X., Li F. and Liu B. (2018). The genetic polymorphisms of TGFβ superfamily genes are associated with litter size in a Chinese indigenous sheep breed (Hu sheep). Anim. Reprod. Sci189, 19-29.

Wang W., Liu S., Li F., Pan X., Li C., Zhang X., Ma Y., La Y., Xi R. and Li T. (2015). Polymorphisms of the ovine BMPR-IB, BMP-15 and FSHR and their associations with litter size in two Chinese indigenous sheep breeds. Int. J. Mol. Sci16(5), 11385-11397.

Wilson T., Wu X.Y., Juengel J.L., Ross I.K., Lumsden J.M., Lord, E.A., Dodds K.G., Walling G.A., McEwan J.C., O’Connell A.R., McNatty K.P. and Montgomery G.W. (2001). Highly prolific Booroola sheep have a mutation in the intracellular kinase domain of bone morphogenetic protein IB receptor (ALK-6) that is expressed in both oocytes and granulosa cells. Biol. Reprod. 64(4), 1225-1235.

Zamani P., Abdoli R., Deljou A. and Rezvan H. (2015). Polymorphism and bioinformatics analysis of growth differentiation factor 9 gene in Lori sheep. Ann. Anim. Sci. 15(2), 337-348.

Zhang C.S., Geng L.Y., Du L.X., Liu Z.Z., Fu Z.X., Feng M.S. and Gong Y.F. (2011). Polymorphic Study of FecXG, FecGH and FecB mutations in four domestic sheep breeds in the lower yellow river valley of China. J. Anim. Vet. Adv. 10(17), 2198-2201.

Zuo B.Y., Qian H.G., Liu J.S., Li Y.H., Wuyun B., Aiken J.,Ying S.J., Wang Z.Y., Noor N., Wang X., Guo Z.Q. and Wang F. (2012). Detection of the 10 mutations of BMPRIB, BMP15 and GDF9 gene in German Mutton Merino sheep. J. Nanjing Agric. Univ. 35(3), 114-120.

Zuo B.Y., Qian H.G., Wang Z.Y., Wang X., Noor N., Bayier A., Ying S.J., Hu X.L., Gong C.G., Guo Z.Q. and Wang F. (2013). A study on BMPR-IB genes of Bayanbulak sheep. Asian-Australasian J. Anim. Sci. 26(1), 36-42.