Discovery of single nucleotide polymorphisms in bone morphogenetic protein (BMP) genes of goats by double digest restriction-site associated DNA sequencing
Shalu Elizabeth Simon![https://orcid.org/0000-0002-7791-8667](/media/client/orcid_16x16.png)
A Department of Animal Genetics, College of Veterinary and Animal Sciences, Mannuthy, Thrissur 680651, India.
B Department of Veterinary Physiology, College of Veterinary and Animal Sciences, Mannuthy, Thrissur 680651, India.
C Corresponding author. Email: shalusimon4@gmail.com
Animal Production Science 61(7) 630-636 https://doi.org/10.1071/AN20013
Submitted: 8 January 2020 Accepted: 28 January 2021 Published: 4 March 2021
Abstract
Context: Two native goat breeds from Kerala, Malabari and Attappady Black, differ significantly in prolificacy (i.e. no. of kids born/kidding). Prolificacy is an important economic trait and the subject of genetic research showing that bone morphogenetic protein (BMP) genes have a significant effect. Double digest restriction-site associated DNA sequencing (ddRADseq) is a highly efficient and low cost technology for high density discovery of single nucleotide polymorphisms (SNPs), which could serve as predictive markers for animal breeding programs.
Aims: The study was aimed at discovering SNPs in BMP genes that affect prolificacy, using ddRADseq followed by validation of selected SNP.
Methods: Blood DNA samples of 10 highly prolific Malabari and 10 less prolific Attappady Black goats were pooled by group and subjected to ddRADseq. SNPs observed in BMP genes were investigated and compared between groups. A validation study was done for the c.614–32789C>T variant in 100 Malabari and 50 Attappady Black goats by using PCR-RFLP.
Key results: In total, 6333 variants were identified by ddRADseq. Three variants were identified in BMP genes, which included two intronic variants c.614–32789C>T and c.490+6793T>C, in genes BMP6 and BMP5 and a downstream gene variant near the BMPR1B gene. According to ddRADseq data, variants in BMP5 and BMP6 differed in allelic distribution between Malabari and Attappady Black goats. For c.490+6793T>C in BMP5, the CC genotype was predominant in the highly prolific Malabari whereas TC was present in the Attappady Black group. The variant c.614–32789C>T in BMP6 was genotyped as TC in Malabari and CC in Attappady Black goats by ddRADseq. This variant was predicted to have an effect on splicing, according to the tool SplicePort. On the basis of bioinformatics analysis and the role of BMP6 gene in follicular dynamics, the variant in BMP6 was selected for further validation studies. All three genotypes were identified by PCR-RFLP; the C allele was the rare allele in the population with an allele frequency of 0.36. Presence of both alleles C and T and the three genotypes CC, TC and TT in this larger population substantiated the robustness of ddRADseq technique.
Conclusions: The technique discovered high confidence SNPs, which could be used for further validation and association studies to develop markers for selection of animals and for genetic improvement of this complex trait.
Implications: Techniques such as ddRADseq provide a large number of SNPs, and investigation of those polymorphisms found across the genome will help to identify new loci affecting traits of interest. This, in turn, will aid in exploring genetically complex traits in a faster and cheaper manner.
Keywords: goats, ddRADseq, BMPs, PCR-RFLP, prolificacy, SNPs.
References
Abdoli R, Zamani P, Mirhoseini SZ, Ghavi Hossein‐Zadeh N, Nadri S (2016) A review on prolificacy genes in sheep. Reproduction in Domestic Animals 51, 631–637.| A review on prolificacy genes in sheep.Crossref | GoogleScholarGoogle Scholar | 27491513PubMed |
Ahlawat S, Sharma R, Maitra A, Roy M, Tantia MS (2014) Designing, optimization and validation of tetra-primer ARMS PCR protocol for genotyping mutations in caprine Fec genes Meta Gene 2, 439–449.
| Designing, optimization and validation of tetra-primer ARMS PCR protocol for genotyping mutations in caprine Fec genesCrossref | GoogleScholarGoogle Scholar | 25606428PubMed |
Ahlawat S, Sharma R, Maitra A, Raja KN, Verma NK, Tantia MS (2015) Prolificacy in Indian goat breeds is independent of Fec B mutation. The Indian Journal of Animal Sciences 85, 64–67.
Al-Musawi SL, Gladwell RT, Knight PG (2007) Bone morphogenetic protein-6 enhances gonadotrophin-dependent progesterone and inhibin secretion and expression of mRNA transcripts encoding gonadotrophin receptors and inhibin/activin subunits in chicken granulosa cells. Reproduction (Cambridge, England) 134, 293–306.
| Bone morphogenetic protein-6 enhances gonadotrophin-dependent progesterone and inhibin secretion and expression of mRNA transcripts encoding gonadotrophin receptors and inhibin/activin subunits in chicken granulosa cells.Crossref | GoogleScholarGoogle Scholar |
Altshuler D, Pollara VJ, Cowles CR, Van Etten WJ, Baldwin J, Linton L, Lander ES (2000) An SNP map of the human genome generated by reduced representation shotgun sequencing. Nature 407, 513–516.
| An SNP map of the human genome generated by reduced representation shotgun sequencing.Crossref | GoogleScholarGoogle Scholar | 11029002PubMed |
Bashir BP, Venkatachalapathy RT (2017) Study on supply chains of goats in northern Kerala. Advances in Animal and Veterinary Sciences 5, 395–399.
| Study on supply chains of goats in northern Kerala.Crossref | GoogleScholarGoogle Scholar |
Bindu KA, Raghavan KC (2010) Haemoglobin polymorphism in Malabari goats. The Indian Veterinary Journal 77, 966–968.
Bragdon B, Moseychuk O, Saldanha S, King D, Julian J, Nohe A (2011) Bone morphogenetic proteins: a critical review. Cellular Signalling 23, 609–620.
| Bone morphogenetic proteins: a critical review.Crossref | GoogleScholarGoogle Scholar | 20959140PubMed |
Chen DI, Zhao M, Mundy GR (2004) Bone morphogenetic proteins. Growth Factors 22, 233–241.
| Bone morphogenetic proteins.Crossref | GoogleScholarGoogle Scholar | 15621726PubMed |
Davis GH, Montgomery GW, Allison AJ, Kelly RW, Bray AR (1982) Segregation of a major gene influencing fecundity in progeny of Booroola sheep. New Zealand Journal of Agricultural Research 25, 525–529.
| Segregation of a major gene influencing fecundity in progeny of Booroola sheep.Crossref | GoogleScholarGoogle Scholar |
Frota IM, Leito CC, Costa JJ, VandenHurk R, Saraiva MV, Figueiredo JR, Silva JR (2013) Levels of BMP-6 mRNA in goat ovarian follicles and in vitro effects of BMP-6 on secondary follicle development. Zygote (Cambridge, England) 21, 270–278.
| Levels of BMP-6 mRNA in goat ovarian follicles and in vitro effects of BMP-6 on secondary follicle development.Crossref | GoogleScholarGoogle Scholar |
Garcia EV, Valdecantos PA, Barrera D, Roldan-Olarte M, Miceli DC (2014) Bone morphogenetic proteins in the bovine oviduct: differential expression of BMP-5 in the isthmus during the estrous cycle. Theriogenology 81, 1032–1041.
| Bone morphogenetic proteins in the bovine oviduct: differential expression of BMP-5 in the isthmus during the estrous cycle.Crossref | GoogleScholarGoogle Scholar | 24582268PubMed |
Glister C, Kemp CF, Knight PG (2004) Bone morphogenetic protein (BMP) ligands and receptors in bovine ovarian follicle cells: actions of BMP-4,-6 and-7 on granulosa cells and differential modulation of Smad-1 phosphorylation by follistatin. Reproduction 127, 239–254.
| Bone morphogenetic protein (BMP) ligands and receptors in bovine ovarian follicle cells: actions of BMP-4,-6 and-7 on granulosa cells and differential modulation of Smad-1 phosphorylation by follistatin.Crossref | GoogleScholarGoogle Scholar | 15056790PubMed |
Godara A (2010) Polymorphism in BMP-15 gene in Marwari goat. PhD Thesis, Rajasthan University of Veterinary and Animal Sciences, Bikaner, India.
Haldar A, Pal P, Datta M, Paul R, Pal SK, Majumdar D, Biswas CK, Pan S (2014) Prolificacy and is relationship with age, body weight, parity, previous litter size and body linear type traits in meat-type goats. Asian-Australasian Journal of Animal Sciences 27, 628–634.
| Prolificacy and is relationship with age, body weight, parity, previous litter size and body linear type traits in meat-type goats.Crossref | GoogleScholarGoogle Scholar | 25049997PubMed |
Hua GH, Chen SL, Ai JT, Yang LG (2008) None of polymorphism of ovine fecundity major genes FecB and FecX was tested in goat. Animal Reproduction Science 108, 279–286.
| None of polymorphism of ovine fecundity major genes FecB and FecX was tested in goat.Crossref | GoogleScholarGoogle Scholar | 17964743PubMed |
Kai W, Nomura K, Fujiwara A, Nakamura Y, Yasuike M, Ojima N, Masaoka T, Ozaki A, Kazeto Y, Gen K, Nagao J (2014) A ddRAD-based genetic map and its integration with the genome duplication. BMC Genomics 15, 233–239.
| A ddRAD-based genetic map and its integration with the genome duplication.Crossref | GoogleScholarGoogle Scholar | 24669946PubMed |
Kess T, Gross J, Harper F, Boulding EG (2016) Low-cost ddRAD method of SNP discovery and genotyping applied to the periwinkle Littorina saxatilis. The Journal of Molluscan Studies 82, 104–109.
Kolodziejczyk SM, Hall BK (1996) Signal transduction and TGF-β superfamily receptors. The International Journal of Biochemistry & Cell Biology 74, 299–314.
| Signal transduction and TGF-β superfamily receptors.Crossref | GoogleScholarGoogle Scholar |
Marker PC, Seung K, Bland AE, Russell LB, Kingsley DM (1997) Spectrum of Bmp5 mutations from germline mutagenesis experiments in mice. Genetics 145, 435–443.
| Spectrum of Bmp5 mutations from germline mutagenesis experiments in mice.Crossref | GoogleScholarGoogle Scholar | 9071596PubMed |
McNatty KP, O’Keefe LE, Henderson KM, Heath DA, Lun S (1986) 125I-labelled hCG binding characteristics in theca interna and other tissues from Romney ewes and from Booroola×Romney ewes with and without a major gene influencing their ovulation rate. Journal of Reproduction and Fertility 77, 477–488.
| 125I-labelled hCG binding characteristics in theca interna and other tissues from Romney ewes and from Booroola×Romney ewes with and without a major gene influencing their ovulation rate.Crossref | GoogleScholarGoogle Scholar | 3090244PubMed |
Miller MR, Dunham JP, Amores A, Cresko WA, Johnson EA (2007) Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Research 17, 240–248.
| Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers.Crossref | GoogleScholarGoogle Scholar | 17189378PubMed |
Miyoshi T, Otsuka F, Inagaki K, Otani H, Takeda M, Suzuki J, Goto J, Ogura T, Makino H (2007) Differential regulation of steroidogenesis by bone morphogenetic proteins in granulosa cells: involvement of extracellularly regulated kinase signaling and oocyte actions in follicle-stimulating hormone-induced estrogen production. Endocrinology 148, 337–345.
| Differential regulation of steroidogenesis by bone morphogenetic proteins in granulosa cells: involvement of extracellularly regulated kinase signaling and oocyte actions in follicle-stimulating hormone-induced estrogen production.Crossref | GoogleScholarGoogle Scholar | 17008391PubMed |
Monget P, Fabre S, Mulsant P, Lecerf F, Elsen JM, Mazerbourg S, Pisselet C, Monniaux D (2002) Regulation of ovarian folliculogenesis by IGF and BMP system in domestic animals. Domestic Animal Endocrinology 23, 139–154.
| Regulation of ovarian folliculogenesis by IGF and BMP system in domestic animals.Crossref | GoogleScholarGoogle Scholar | 12142233PubMed |
Monteiro A, Costa JM, Lima MJ (2017) Goat system productions: Advantages and disadvantages to the animal, environment and farmer. Goat Science 351–366.
Moore RK, Otsuka F, Shimasaki S (2003) Molecular basis of bone morphogenetic protein-15 signaling in granulosa cells. The Journal of Biological Chemistry 278, 304–310.
| Molecular basis of bone morphogenetic protein-15 signaling in granulosa cells.Crossref | GoogleScholarGoogle Scholar | 12419820PubMed |
Mora-Marquez F, Garcia Olivweres V, Emerson BC, Lopez de Heredia U (2017) ddradseqtools: a softwwere package for in silico simulation and testing of double digest RAD seq experiments. Molecular Ecology Research 17, 230–246.
| ddradseqtools: a softwwere package for in silico simulation and testing of double digest RAD seq experiments.Crossref | GoogleScholarGoogle Scholar |
Naicy T, Venkatachalapathy RT, Aravindakshan TV, Radhika G, Raghavan KC, Mini M, Shyama K (2016) Nerve Growth Factor gene ovarian expression, polymorphism identification, and association with litter size in goats. Theriogenology 86, 2172–2178.
| Nerve Growth Factor gene ovarian expression, polymorphism identification, and association with litter size in goats.Crossref | GoogleScholarGoogle Scholar | 27544869PubMed |
Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstre HE (2012) Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS ONE 7, e37135
| Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species.Crossref | GoogleScholarGoogle Scholar | 22675423PubMed |
Pfendler KC, Yoon J, Taborn GU, Kuehn M, Rm Iannaccone PM (2000) Nodal and bone morphogenetic protein 5 interact in murine mesoderm formation and implantation. Genesis (New York, N.Y.) 28, 1–14.
| Nodal and bone morphogenetic protein 5 interact in murine mesoderm formation and implantation.Crossref | GoogleScholarGoogle Scholar |
Polley S, De S, Batabyal S, Kaushik R, Yadav P, Arora JS, Chattopadhyay S, Pan S, Brahma B, Datta TK, Goswami SL (2009) Polymorphism of fecundity genes(BMPR1B, BMP15 and GDF9) in the Indian prolific Black Bengal goat. Small Ruminant Research 85, 122–129.
| Polymorphism of fecundity genes(BMPR1B, BMP15 and GDF9) in the Indian prolific Black Bengal goat.Crossref | GoogleScholarGoogle Scholar |
Radhika G (2015) Diversity analysis among goat genetic groups of Kerala. PhD Thesis, Kerala Veterinary and Animal Sciences University, Pookode, India.
Shi J, Yoshino O, Osuga Y, Koga K, Hirota Y, Hirata T, Yano T, Nishii O, Taketani Y (2009) Bone morphogenetic protein-6 stimulates gene expression of follicle-stimulating hormone receptor, inhibin/activin β subunits, and anti-Mullerian hormone in human granulosa cells. Fertility and Sterility 92, 1794–1798.
| Bone morphogenetic protein-6 stimulates gene expression of follicle-stimulating hormone receptor, inhibin/activin β subunits, and anti-Mullerian hormone in human granulosa cells.Crossref | GoogleScholarGoogle Scholar | 19539911PubMed |
Shimasaki S, Moore RK, Erickson GF, Otsuka F (2003) The role of bone morphogenetic proteins in ovarian function. Reproduction 61, 323–337.
Stephen M, Raja TV, Sosamma I (2005) Survey and characterization of Attappady Black goats of Kerala, India Animal Genetic Resources 37, 43–52.
| Survey and characterization of Attappady Black goats of Kerala, IndiaCrossref | GoogleScholarGoogle Scholar |
Wang X, Su L, Pan X, Yao J, Li Z, Wang X, Xu B (2015) Effect of BMP-6 on development and maturation of mouse preantral follicles in vitro. Biotechnolology and Biotechnological Equipments 29, 336–344.
| Effect of BMP-6 on development and maturation of mouse preantral follicles in vitro.Crossref | GoogleScholarGoogle Scholar |
Yu H, You X, Li J, Liu H, Meng Z, Xiao L, Zhang H, Lin HR, Zhang Y, Shi Q (2016) Genome-wide mapping of growth-related quantitative trait loci in orange-spotted grouper (Epinepheluscoioides) using double digest restriction-site associated DNAsequencing (ddRADseq). International Journal of Molecular Science 17, 501–509.
| Genome-wide mapping of growth-related quantitative trait loci in orange-spotted grouper (Epinepheluscoioides) using double digest restriction-site associated DNAsequencing (ddRADseq).Crossref | GoogleScholarGoogle Scholar |