Genotyping of Novel SNPs in BMPR1B,BMP15, and GDF9 Genes for Association with Prolificacy in Seven Indian Goat Breeds

Goats form the backbone of rural livelihood and financial security systems in India but their population is showing decreasing trend. Improvement of reproductive traits such as prolificacy offers a solution to stabilize the decreasing goat population and to meet the nutritional needs of growing human population. In the present study, six novel SNPs in three candidate genes for prolificacy (BMPR1B, BMP15, and GDF9) were genotyped in seven breeds of Indian goats to evaluate their association with litter size. Tetra primer ARMS-PCR and PCR-RFLP based protocols were developed for genotyping six novel SNPs, namely, T(-242)C in BMPR1B; G735A and C808G in BMP15; and C818T, A959C, and G1189A in GDF9 gene. The effect of breed was highly significant (p ≤ 0.01) on litter size but the effect of genotype was nonsignificant. The effect of parity on litter size was also significant in the prolific Black Bengal breed. The litter size differences observed between breeds are attributed to breed differences. Novel mutations observed at different loci in GDF9, BMP15, and BMPR1B genes do not contribute to the reproductive capability of the investigated breeds. Further studies with more number of breeds and animals exploring association of these novel SNPs with reproductive traits may be fruitful.     

小尾寒羊高繁殖力候选基因BMP15和GDF 9的研究

以控制Belclare和Cambridge绵羊高繁殖力的骨形态发生蛋白 15 (bonemorphogeneticprotein 15 ,BMP15 )基因和生长分化因子 9(growthdifferentiationfactor 9,GDF9)基因为候选基因 ,采用PCR RFLP技术检测BMP15基因和GDF9基因在高繁殖力绵羊品种 (小尾寒羊、湖羊 )以及低繁殖力绵羊品种 (多赛特羊、特克塞尔羊、德国肉用美利奴羊 )中的单核苷酸多态性 ,同时研究这两个基因对小尾寒羊高繁殖力的影响。结果表明 :在 5个绵羊品种中都没有检测到GDF9基因的G8突变 (C→T) ,也没有检测到BMP15基因的B4突变 (G→T)。高繁殖力的小尾寒羊在BMP15基因编码序列第 718位碱基处发生了与Belclare绵羊和Cambridge绵羊相同的B2突变 (C→T) ,而其余 4个绵羊品种则没有发生这种突变。对于BMP15基因的B2突变 ,在小尾寒羊中检测到AA、AB两种基因型 ,A等位基因频率为 0 734,B等位基因频率为 0 2 6 6。小尾寒羊与其余 4个绵羊品种间B2突变基因型分布差异极显著 (P <0 0 0 1)。突变杂合基因型 (AB)小尾寒羊平均产羔数比野生纯合基因型 (AA)多 0 6 2只 (P <0 0 1)。研究结果表明 ,BMP15B2突变对小尾寒羊高繁殖力影响作用十分明显 ,同时排除了GDF9G8突变和BMP15B4突变影响小尾寒羊高繁殖力的可能性     

Patterns of expression of messenger RNAs encoding GDF9, BMP15, TGFBR1, BMPR1B, and BMPR2 during follicular development and characterization of ovarian follicular populations in ewes carrying the Woodlands FecX2W mutation

Woodlands sheep have a putative genetic mutation (FecX2(W)) that increases ovulation rate. At present, the identity of FecX2(W) is unknown. The trait does not appear to be due to the previously described mutations in bone morphogenetic protein 15 (BMP15), growth differentiation factor 9 (GDF9), or bone morphogenetic protein receptor type 1B (BMPR1B) that affect ovulation rate in sheep. Potentially, FecX2(W) could be an unidentified genetic mutation in BMP15 or in the closely related GDF9, which interacts with BMP15 to control ovarian function. Alternatively, FecX2(W) may affect ovulation rate by changing the expression patterns in the molecular pathways activated by genes known to regulate ovulation rate. The objectives of these experiments were to sequence the complete coding region of the BMP15 and GDF9 genes, determine the patterns of expression of mRNAs encoding GDF9, BMP15, TGFBR1, BMPR1B, and BMPR2 during follicular development, and characterize the follicular populations in ewes heterozygous for the Woodlands mutation and their wild-type contemporaries. No differences in the coding sequences of BMP15 or GDF9 genes were identified that were associated with enhanced ovulation rate. The expression patterns of GDF9 and BMPR2 mRNAs were not different between genotypes. However, expression of BMP15 mRNA was less in oocytes of FecX2(W) ewes in large preantral and antral follicles. Expression of ALK5 mRNA was significantly higher in the oocytes of FecX2(W) ewes, whereas expression of BMPR1B was decreased in both oocytes and granulosa cells of FecX2(W) ewes. FecX2(W) ewes also had increased numbers of antral follicles <1 mm in diameter. These follicles were smaller in average diameter, with the oocytes also being of a smaller mean diameter. Given that a mutation in BMP15 or BMPR1B results in increased ovulation rates in sheep, the differences in expression levels of BMP15 and BMPR1B may play a role in the increase in ovulation rate observed in Woodlands ewes with the FecX2(W) mutation.     

The proregion of mouse BMP15 regulates the cooperative interactions of BMP15 and GDF9

Bone morphogenetic protein 15 (BMP15) and growth and differentiation factor 9 (GDF9) are secreted by the mammalian oocyte and are essential for ovarian follicular development, ovulation, and fertility. However, the secreted forms of the BMP15 and GDF9 proteins and the nature of cooperative molecular interactions between BMP15 and GDF9 previously reported have not been fully characterized. In this study, we found that recombinant mouse BMP15 and GDF9 are secreted as cleaved mature and proregion proteins, with BMP15 also secreted as uncleaved promature protein. Noncovalent interactions were identified between the mature and proregion proteins of each growth factor. Moreover, GDF9 mature protein was found to coimmunoprecipitate with the BMP15 proregion, suggestive of a heteromeric association between BMP15 and GDF9. Mouse GDF9 was found to exist mostly as a dimer of mature protein, in both the presence and absence of BMP15. In contrast, BMP15 formed mostly multimers of proregion and mature protein when combined with GDF9, providing further evidence for heteromeric interaction. Mouse BMP15 was found to act cooperatively with GDF9 in a rat granulosa cell thymidine incorporation bioassay and to signal through the BMPR2 and ACVR1B/TGFBR1/ACVR1C receptor-mediated pathways. Immunoneutralization experiments using GDF9 mature protein antibody indicated that these cooperative interactions are species specific. Additionally, immunoneutralization with proregion antibodies highlighted the involvement of the BMP15 proregion in BMP15/GDF9 cooperative interactions. Taken together, these findings support a novel hypothesis where the extracellular cooperative interactions of recombinant mouse BMP15 and GDF9 are multimeric, involving the proregion of BMP15, and may well be species specific.     

绵羊GDF9和BMP15基因多态性检测

以绵羊GDF9基因FecG^H突变和BMP15基因FecX^B和FecX^G突变为候选基因,采用PCR—RFLP方法研究其在湖羊、夏洛来、陶赛特、萨福克、中国美利奴肉用多胎品系、中国美利奴羊和罗米丽羊7个品种中的多态性．结果发现,湖羊群体中存在GDF和BMP15基因的FecG^H和FecX^B突变,但发生率极低,分别为0．645％(2／310)和0．968％(3／310)：而其它品种中则没有发现GDF9和BMP15基因的相应突变．这一发现对于建立绵羊基因标记辅助选择方法具有重要意义．     

The effects of immunizing sheep with different BMP15 or GDF9 peptide sequences on ovarian follicular activity and ovulation rate

The aims of these studies were to determine the abilities of antisera against different regions of ovine bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) to inhibit ovarian follicular activity, estrus (mating), and ovulation in sheep. The 9-15-mer peptides were conjugated to keyhole limpet hemocyanin (KLH) and used to generate antibodies against the flexible N-terminal regions of the mature protein as well as against regions in which dimerization of the protein or interaction with a type 1 BMP or a type 2 TGFB or BMP receptor was predicted to occur. Ewes (n = 10 per treatment group) were vaccinated with KLH or the KLH-BMP15 (n = 9 different peptides) or KLH-GDF9 (n = 10) peptides in Freund adjuvant at five consecutive monthly intervals. Overall, antisera generated against peptides that corresponded to amino acid residues 1-15 of the N-terminus of the BMP15 or GDF9 mature protein or GDF9 amino acid residues 21-34 were the most potent at inhibiting ovulation following primary and single booster vaccination. Several other BMP15 (8/9) or GDF9 (6/10) treatment groups, but not KLH alone, also produced significant reductions in the numbers of animals that ovulated, although 2, 3 or 4 booster vaccinations were required. Anovulation was commonly associated with the inhibition of normal ovarian follicular development and anestrus. The in vitro neutralization studies with IgG from the BMP15 or GDF9 immunized ewes showed that the mean inhibition of BMP15 plus GDF9 stimulation of (3)H-thymidine uptake by rat granulosa cells was approximately 70% for animals without corpora lutea (CL), whereas for animals with one to three CL or more than three CL, the inhibition was 24%-33% or 27%-42%, respectively. In summary, these data suggest that reagents that block the biological actions of BMP15 or GDF9 at their N-termini have potential as contraceptives or sterilizing agents.     

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)

Belclare and Cambridge are prolific sheep breeds, the origins of which involved selecting ewes with exceptionally high litter size records from commercial flocks. The variation in ovulation rate in both breeds is consistent with segregation of a gene (or genes) with a large effect on this trait. Sterile ewes, due to a failure of normal ovarian follicle development, occur in both breeds. New naturally occurring mutations in genes for the oocyte-derived growth factors growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP15) are described. These mutations are associated with increased ovulation rate in heterozygous carriers and sterility in homozygous carriers in both breeds. This is the first time that a mutation in the gene for GDF9 has been found that causes increased ovulation rate and infertility in a manner similar to inactivating mutations in BMP15, and shows that GDF9 is essential for normal folliculogenesis in sheep. Furthermore, it is shown, for the first time in any species, that individuals with mutations in both GDF9 and BMP15 have a greater ovulation rate than sheep with either of the mutations separately.     

Maximum‐likelihood approaches reveal signatures of positive selection in BMP15 and GDF9 genes modulating ovarian function in mammalian female fertility

Bone morphogenetic proteins (BMPs) and the growth factors (GDFs) play an important role in ovarian folliculogenesis and essential regulator of processes of numerous granulosa cells. BMP15 gene variations linked to various ovarian phenotypic consequences subject to the species, from infertility to improved prolificacy in sheep, primary ovarian insufficiency in women or associated with minor subfertility in mouse. To study the evolving role of BMP15 and GDF9, a phylogenetic analysis was performed. To find out the candidate gene associated with prolificacy in mammals, the nucleotide sequence of BMP15 and GDF9 genes was recognized under positive selection in various mammalian species. Maximum‐likelihood approaches used on BMP15 and GDF9 genes exhibited a robust divergence and a prompted evolution as compared to other TGFβ family members. Furthermore, among 32 mammalian species, we identified positive selection signals in the hominidae clade resulting to 132D, 147E, 163Y, 191W, and 236P codon sites of BMP15 and 162F, 188K, 206R, 240A, 244L, 246H, 248S, 251D, 253L, 254F and other codon sites of GDF9. The positively selected amino acid sites such as Alanine, Lucien, Arginine, and lysine are important for signaling. In conclusion, this study evidences that GDF9 and BMP15 genes have rapid evolution than other TGFß family members and was subjected to positive selection in the mammalian clade. Selected sites under the positive selection are of remarkable significance for the particular functioning of the protein and consequently for female fertility.     

Regulation of GDF9 and CDKN1B expression in Tibetan sheep testes during different stages of maturity

Normal spermatogenesis is heavily dependent on the balance of germ cell proliferation, differentiation and apoptosis. Growth differentiation factor 9 (GDF9) and cyclin-dependent kinase inhibitor 1 B (CDKN1B) are strongly associated with cell cycle transition from G0/G1 to S and G2/M phase and hence regulating the growth and development of testicular germ cells and somatic cells. The current study was aimed at seeking out scientific evidence to determine if GDF9 and CDKN1B gene expression functions in the development of Tibetan sheep testes. To this end, developmental testes were derived from three-month-old (pre-puberty), one-year-old (sexual maturity), and three-year-old (adult) Tibetan sheep and then the expression and localization patterns of GDF9 and CDKN1B in these testes were evaluated using quantitative real-time PCR (qRT-PCR), Western blot and immunofluorescence. qRT-PCR and Western blot results showed that GDF9 and CDKN1B were detected in the testes throughout the different developmental stages. The abundance of GDF9 mRNA and protein in the testes of one- and three-year-old Tibetan sheep were higher than that in the testes of three-month-old Tibetan sheep; the mRNA and protein abundance of the CDKN1B gene in three-month-old Tibetan sheep testes were higher than that in the testes of the one-and three-year-old sheep. Moreover, immunofluorescence results suggested that the GDF9 protein was expressed in spermatogonia and Leydig cells, and that the CDKN1B protein was localized mainly in Leydig cells with some in the seminiferous epithelium throughout developmental stages. This indicated a novel role of the GDF9 and CDKN1B genes in Leydig cell development over and above their known roles in germ cell development. These findings have significant implications for our understanding of the molecular mechanisms of GDF9 and CDKN1B genes in Tibetan sheep spermatogenesis.     

Polymorphism study of growth differentiation factor 9B (GDF9B) gene and its association with reproductive traits in sheep

GDF9B protein plays a critical role in growth and differentiation of early ovarian follicles. In Inverdale and Hanna sheep, mutations in exon-2 of GDF9B gene have been recorded to show increased ovulation rate in heterozygous condition whereas homozygotes are infertile. Present screen study was carried out to explore the presence of these reported mutations in Corriedale and Local Kashmir Valley sheep with high rate of twinning. Exon-2 of GDF9B gene was amplified and the polymorphism was explored by SSCP technique. In the process three different bandings were observed. Later on these patterns corresponded with three different allelic forms on nucleotide sequencing. Phylogenetic analysis revealed that the nucleotide sequences of alleles observed in the present study and that of a published sequence of sheep were having the same point of origin. The results were also compared with goats, large ruminants and humans. The allelic frequencies of allele A and B were 0.64 and 0.36, respectively in Corriedale sheep whereas the allelic frequencies of all the three alleles in Kashmir Valley sheep were 0.60, 0.34 and 0.06. SNP “C” of the designated genotype AC was observed to pronounce a significant effect on litter size with average litter size going up by 0.63 as compared with the nearest genotype AB wherein the litter size was 1.29 ± 0.05. The average litter size between AA and AB genotypes did not vary significantly.     

Quantitative analysis of bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) gene expression in calf and adult bovine ovaries

Background  

It has been reported that calf oocytes are less developmentally competent than oocytes obtained from adult cows. Bone morphogenetic protein 15 (BMP15) and growth and differentiation factor 9 (GDF9) play critical roles in folliculogenesis, follicular development and ovulation in mammalian ovaries. In the present study, we attempted to compare the expression patterns of BMP15 and GDF9 in the cells of calf and cow ovaries to determine a relationship between the level of these genes and the low developmental competence of calf oocytes.     

Allelic frequency and genotypes of prion protein at codon 136 and 171 in Iranian Ghezel sheep breeds

PrP genotypes at codons 136 and 171 in 120 Iranian Ghezel sheep breeds were studied using allele-specific PCR amplification and compared with the well-known sheep breeds in North America, the United States and Europe. The frequency of V allele and VV genotype at codon 136 of Ghezel sheep breed was significantly lower than AA and AV. At codon 171, the frequency of allele H was significantly lower than Q and R. Despite the similarities of PrP genotypes at codons 136 and 171 between Iranian Ghezel sheep breeds and some of the studied breeds, significant differences were found with others. Planning of effective breeding control and successful eradication of susceptible genotypes in Iranian Ghezel sheep breeds will not be possible unless the susceptibility of various genotypes in Ghezel sheep breeds to natural or experimental scrapie has been elucidated.Key words: scrapie, Ghezel sheep breed, PrP genotyping, allele specific amplification, codon 136, codon 171Scrapie was first described in England in 1732,¹ and it is an infectious neurodegenerative fatal disease of sheep and goats belonging to the group of transmissible subacute spongiform encephalopathies (TSEs), along with bovine spongiform encephalopathy (BSE), chronic wasting disease and Creutzfeldt-Jakob disease.^2,3 The term prion, proteinaceous infectious particles, coined by Stanley B. Prusiner, was introduced, and he presents the idea that the causal agent is a protein.⁴ Prion proteins are discovered in two forms, the wild-type form (PrP^c) and the mutant form (PrP^Sc).⁵ Although scrapie is an infectious disease, the susceptibility of sheep is influenced by genotypes of the prion protein (PrP) gene.^2,6 Researchers have found that the PrP allelic variant alanine/arginine/arginine (ARR) at codons 136, 154 and 171 is associated with resistance to scrapie in several breeds.^7–14 Most of the sheep populations in the Near East and North African Region (84% of the total population of 255 million) are raised in Iran, Turkey, Pakistan, Sudan, Algeria, Morocco, Afghanistan, Syria and Somalia.¹⁵ In 2003, the Iranian sheep population was estimated at 54,000,000 head. The Ghezel sheep breed, which also is known as Kizil-Karaman, Mor-Karaman, Dugli, Erzurum, Chacra, Chagra, Chakra, Gesel, Gezel, Kazil, Khezel, Khizel, Kizil, Qezel, Qizil and Turkish Brown, originated in northwestern Iran and northeastern Turkey. By considering sheep breeds as one of the main sources of meat, dairy products and related products, a global screening attempt is started in different areas. In compliance with European Union Decision 2003/100/EC, each member state has introduced a breeding program to select for resistance to TSEs in sheep populations to increase the frequency of the ARR allele. A similar breeding program is established in United States and Canada. The Near East and North African Region still needs additional programs to help the global plan of eradication of scrapie-susceptible genotypes. The current study was the first to assess the geographical and molecular variation of codons 136 and 171 polymorphism between Iranian Ghezel sheep breed and well-known sheep breeds.Polymorphism at codon 136 is associated with susceptibility to scrapie in both experimental and natural models.^10,11,13,16 17 and Austrian Carynthian sheep.¹⁸ Swiss White Alpine showed higher frequency of allele V at position 136 than Swiss Oxford Down, Swiss Black-Brown Mountain and Valais Blacknose.¹⁹ Comparison of polymorphism at codon 136 in the current study with some of other breeds (20 some flock of Hampshire sheep²¹ with current study, but the frequency of it is higher than that of some other breeds.

Table 1

Comparison of PrP allelic and genotype frequencies at codon 136 in different breeds