The objective of this study was to investigate association between GDF9 and BMP15 gene polymorphism and litter size in fat-tailed sheep, a total of 97 mature ewes from four breeds (Afshari=19; Baluchi=18; Makui=30 and Mehraban=30) were genotyped for the BMP15 HinfI and GDF9 HhaI polymorphisms by PCR-RFLP technique. The highest and lowest mutant allele frequencies were found in Makui (0.27) and Afshari (0.10) sheep for the BMP15 gene and in Afshari (0.24) and Mehraban (0.18) sheep for the GDF9 gene, respectively. Litter size was significantly influenced by genotype of the ewe for two genes (P < 0.01). Heterozygous genotypes for both loci showed higher litter size than homozygous genotypes (P < 0.01). None of the individuals carried homozygous genotype for both of the GDF9 and BMP15 variants in these breeds. The individuals carrying the mutant allele for one of the investigated candidate gene still showed fertile phenotype. Thus, existence of homozygosity at one of the BMP15 and GDF9 variant is not probably able to block normal hormonal pathway of reproduction in fat-tailed sheep.
A major gene for bovine ovulation rate has been mapped to a 1.2 Mb region of chromosome 10. Screening of coding regions of positional candidate genes within this region failed to reveal a causative polymorphism, leading to the hypothesis that the phenotype results from differences in candidate gene expression rather than alteration of gene structure. This study tested differences in expression of positional candidate genes in granulosa cells between carriers and noncarriers of the high fecundity allele, as well as characterizing differences in the transcriptomic profile between genotypes. Five carriers and five noncarriers, female descendants of "Trio," a carrier of the high fecundity allele were initially used in an RNA-seq analysis of gene expression. Four of ten samples were contaminated with theca cells, so that six samples were used in the final analysis (three of each genotype). Of 14 973 genes expressed, 143 were differentially expressed (false discovery rate P < 0.05) in carriers versus noncarriers. Among the positional candidate genes, SMAD6 was 6.6-fold overexpressed in the carriers compared to noncarriers (P < 5 × 10-5). This result was replicated in an independent group of 12 females (7 carriers and 5 noncarriers) using quantitative real-time PCR; SMAD6 was 9.3-fold overexpressed in carriers versus noncarriers (P = 1.17 × 10-6). Association of overexpression of SMAD6, an inhibitor of the BMP/SMAD signaling pathway, with high ovulation rate corresponds well with disabling mutations in ligands (BMP15 and GDF9) and a receptor (BMPR1B) of this pathway that cause increased ovulation rate in sheep.