Polymorphism of the GnRHR gene and its association with litter size in Boer goats
Wucai YangI; Keqiong TangI; Chunyan ZhangI; Deqing XuII; Qunying WenIII; Liguo YangI,#
IKey lab of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education. Huazhong Agricultural University, Wuhan, 430070, China
IIBoer Goat Breeding Station, Yidu, 443300, China
IIIShiyan Municipal Bureau of Animal Husbandry of Hubei Province, Shiyan, 442000, China
]]>
ABSTRACT
The polymorphism in the gonadotropin-releasing hormone receptor (GnRHR) gene was detected by PCR single-strand conformational (PCR-SSCP) and DNA sequencing methods in 224 Boer goats. Two single nucleotide polymorphisms (SNPs), G891T and G757A, were found in the exon I of the goat GnRHR gene and their associations with litter size in Boer goats were evaluated. Association analysis showed that G891T and G757A had significant effects on litter size. In locus G891T, individuals with GT genotype had significantly larger litter sizes compared to GG genotype in the first and third parity. In locus G757A, ewes with GG genotype had significantly higher litter sizes than those of GA in the fourth parity. The present study primarily revealed an association between allele T at 891 locus of the GnRHR gene and high litter size in the first and third parity in Boer goats. These results provide evidence that the GnRHR gene may be significantly correlated with a high prolificacy in goats.
Keywords: GnRHR, litter size, boer goat, SSCP
The Boer goat was developed in South Africa in the early 1900s (Casey & Van Niekerk, 1988). The first importation to China was during 1995, and it is now a popular breed for meat production. The breed is well adapted and established in every region of China. Numerous studies have confirmed that farmers prefer the Boer to other goat breeds for its rapid growth characteristics (Vatankhah & Talebi, 2008; Zhang et al., 2008). It has been accepted, in general, that the economic efficiency of a flock of goats is dependent on total productivity (Khanum et al., 2007), and this productivity is more dependent on fertility and prolificacy of the dams than any other component (Boujenane, 2002). It will, therefore, be meaningful to use a candidate gene approach to identify genes that influence litter size in Boer goats. The hypothalamic gonadotropin-releasing hormone receptor (GnRHR) is a key regulator of the reproductive system, which triggers the synthesis and release of the luteinizing hormones (LH) and follicle-stimulating hormone (FSH) in the pituitary gland (Millar, 2005; Rispoli & Nett, 2005). Previous studies demonstrated that PKC and various downstream signal, transduction cascades were activated following the binding of GnRH to GnRHR (Kraus et al., 2001; Millar & Pawson, 2004). These pathways then differentially regulate the synthesis and release of the gonadotropic hormones; LH and FSH, that in turn regulate the production of gametes and gonadal hormones (Naor, 2009). In the absence of sufficient FSH or LH, follicles fail to develop beyond the early antral stage, and therefore ovulation does not occur (Barnes et al., 2002; Clarkson & Herbiso, 2009). Due to the central role of GnRHR in regulating gonadotropin synthesis and release, the GnRHR gene appeared to be a good candidate for mutations associated with reproductive performance (Bédécarrats et al., 2003). The aim of the present study was to detect new SNP in exon 1 of the GnRHR gene in Boer goats by PCR-SSCP, and to investigate its association with litter size.
Genomic DNA samples were obtained from 224 Boer goats, with records in the Yidu Boer Goat Breeding Station, under common rearing conditions from 2003 to 2007. Genomic DNA was extracted from white blood cells by using the standard phenol-chloroform extraction protocol. Records of litter size for Boer goats from the first to fifth parity, were collected for statistical analysis (Table 1).
According to the DNA sequence of the sheep GnRHR gene (GenBank accession number: L42937.1), a pair of primers (sense: 5-CTCGTTCGCTTTAGCACCC-3; anti-sense: 5- CTGTGGTCCAGCAAA GATG-3) was designed to amplify exon 1 of the goat GnRHR gene. PCR reactions were performed in a 20 µL reaction mixture containing 10 µmol of each primer, 200 µM of dNTP (deoxyribo nucleotide triphosphate), 2.5 µL of 10 × reaction buffer which contained 1.5 mM of MgCl2, 0.5 unit of Taq-DNA polymerase (Promega, Madison, WI) and 50 ng of genomic DNA as template. After denaturation at 94 ºC for 5 min, 32 amplification cycles were performed comprising denaturation at 94 ºC for 45 s, annealing at 62.0 ºC for 45 s, and extension at 72 ºC for 45 s, followed by a further 10 min extension at 72 ºC. PCR products were electrophoretically detected on 1.2% agarose gel (5 V/cm) by staining with ethidium bromide.
Aliquots of 4 µL of the PCR products were mixed with 10 µL of a denaturing solution (95% formamide, 25 mM EDTA, 0.025% xylene-cyanole and 0.025% bromophenol blue), heated for 10 min at 99 ºC and chilled on ice. Denatured DNA was subjected to 10% polyacrylamide gel electrophoresis (PAGE) (39:1 acrilamyde:bisacrylamide) in 1 × tris-borate EDTA (TBE) buffer at constant voltage (120 V) for 10 - 16 h. The gel was stained with 0.1% silver nitrate and the PCR products, with different SSCP patterns, were sequenced in an ABI PRIZM 377 DNA sequencer to identify the SNP. The PCR products were digested with MspI and EarI (TaKaRa, Tokyo, Japan) according to the sequencing results. The digestion mixture contained 4 µL PCR products, 1 × digestion buffer, 4.0 U of each enzyme and digested at 37 ºC overnight. Fragments were separated on 2% agarose gels and visualized with ethidium bromide staining.
]]> The associations of GnRHR genotypes with litter size were calculated by the General Linear Model of SAS 8.1 (SAS Inst. Inc., Cary, NC). The variables, which were measured in one to five parities, were considered as separate traits (Noguera et al., 2002). The fixed effect involved genotype, parity, age of dam at kidding, kidding year, kidding season, and their interaction. Since the goats were breeding at random, the fixed effects of sire and dam within sire were not taken into account (Wu et al., 2009). Their effects on litter size were analyzed using the GLM procedure and compared by using Duncan's Multiple Range Test (SAS 8.1). Only factors (genotype, age, parity and kidding year) that affected the records significantly (P <0.05) were fitted in the final statistical model shown as the following:Yjbcs = µ+ Gj + Ab + Pc + Ks + ejbc
Where Yjbcs is the phenotypic value of traits; µ is the mean; Gj, Ab, PC and Ks are the fixed effects of the genotype, age, parity and kidding year, respectively; ejbcs is the random residual error.
Polymorphisms of the goat GnRHR gene were detected by PCR-SSCP (Figure 1) and DNA sequencing. In order to better describe the different genotypes of G﹥T and G﹥A mutations of the goat GnRHR gene, we named the alleles 891-G, 891-T, 757-G and 757-A, respectively. In locus G891T, the mutation G﹥T (GAAGAG to GAATAG) destroys an EarI endonuclease restriction site (GAAGAG). Therefore, the EarI digestion of amplified products shows two fragments (425 and 324 bp) for 891-G allele and one fragment (749 bp) for 891-T allele (Figure 2). In locus G757A, the mutation G﹥A (CCGG to CCAG) destroys a MspI endonuclease restriction site (CCGG). Therefore, MspI digestion of amplified products shows two fragments (449 and 300 bp) for 757-G allele and one fragment (749 bp) for 757-A allele (Figure 2). Frequencies of 891-G, 891-T, 757-G and 757-A, alleles were 0.951, 0.049, 0.925 and 0.075 in the analyzed population (Table 2). In the present study, the homozygous TT and AA genotypes were not detected in the GnRHR gene, similar to a previous finding in bovines (Yang et al., 2011). We speculated that these genotypes may have a negative effect on reproduction, so that ewes with this genotype might have been eliminated through the process of artificial selection and breeding programme, due to their extremely poor performance in reproduction. It is also possible that TT and AA genotypes were not found in the study population due to its low allele frequency. Therefore, our results should be further validated in larger goat populations.
Litter size was significantly influenced by age, parity and kidding year (P <0.05). The parity of the dam had a significant effect on litter size (P <0.01), least-square means of litter size was the smallest at the first parity, and increased with each parity from 1 to 4, and then decreased in the fifth parity (Table 3). The results of the association analysis between GnRHR genotypes and litter size are given in Table 4. The statistical analysis revealed that those two SNPs (G891T and G757A) only had a significant effect on litter size in certain parities in Boer goats. In locus G891T, individuals with GT genotype significantly outperformed individuals with the GG genotype (P <0.05) with regard to litter sizes of 0.27 and 0.22 kids, respectively, in the first and third parity (Table 4). In locus G757A, ewes with GG genotype had larger (P <0.05) litter sizes by 0.23 kids in the fourth parity than ewes with the GA genotype (Table 4). This may contribute to the limited sample size of the Boer goat in the present study and the results should be further validated in larger goat populations. Although those mutations do not change the amino acid sequences, it probably affects the density of GnRHR on the gonadotropes by influencing the expression of the GnRHR gene and the stability of the GnRHR transcription, and could then result in a decrease or increase of gonadotrophins concentrations (Schubert et al., 2000; Leanos-Miranda et al., 2002). It has been observed that the synonymous SNP can affect in vivo protein folding and consequently its function, as well as gene expression and phenotype (Kimchi-Sarfaty et al., 2007; Ren et al., 2010; Tang et al., 2011).
In summary, our results identified two new polymorphisms (G891T and G757A) in the goat GnRHR gene, and showed an association with allele T at 891 locus between the GnRHR gene and the high litter size in the first and third parity in Boer goats. Due to the lack of functional data, the conclusion was not final/irrefutable and further studies are necessary for the validation of these polymorphisms.
Acknowledgments
This work was supported by programme 863 (No.2008AA101010) and National Key Technology R&D Programme (No. 2006BAD04A02-11).
]]> References
Barnes, R.B., Namnoum, A.B., Rosenfield, R.L. & Layman, L.C., 2002. The role of LH and FSH in ovarian androgen secretion and ovarian follicular development: clinical studies in a patient with isolated FSH deficiency and multicystic ovaries. Hum. Reprod. 17, 88-91. [ Links ]
Bédécarrats, G.Y. & Kaiser, U.B., 2007. Mutations in the human gonadotropin-releasing hormone receptor: insights into receptor biology and function. Semin. Reprod. Med. 25, 368-378. [ Links ]
Boujenane, I., 2002. Development of the DS synthetic breed of sheep in Morocco: ewe reproduction and lamb preweaning growth and survival. Small Rumin. Res. 45, 61-66. [ Links ]
Casey, N.H. & Van Niekerk, W.A., 1988. The boer goat. I. Origin, adaptability, performance testing, reproduction and milk production. Small Rumin. Res. 1, 291-302. [ Links ]
Clarkson, J. & Herbison, A.E., 2009. Oestrogen, kisspeptin, GPR54 and the pre-ovulatory luteinising hormone surge. J. Neuroendocrinol. 21, 305-311. [ Links ]
Khanum, S.A., Hussain, M. & Kausar, R., 2007. Assessment of reproductive parameters in female Dwarf goat (Capra hircus) on the basis of progesterone profiles. Anim. Reprod. Sci. 102, 267-275. [ Links ]
Kimchi-Sarfaty, C., Oh, J.M., Kim, I.W., Sauna, Z.E., Calcagno, A.M., Ambudkar, S.V. & Gottesman, M.M., 2007. A "silent" polymorphism in the MDR1 gene changes substrate specificity. Science 315, 525-528. [ Links ]
Kraus, S., Naor, Z.S. & Eger, R., 2001. Intracellular signaling pathways mediated by the gonadotropin-releasing hormone (GnRH) receptor. Arch. Med. Res. 32, 499-509. [ Links ]
Leanos-Miranda, A., Janovick, J.A. & Conn, P.M., 2002. Receptor-misrouting: an unexpectedly prevalent and rescuable etiology in gonadotropin-releasing hormone receptor-mediated hypogonadotropic hypogonadism. J. Clin. Endocrinol. Metab. 87, 4825-4828. [ Links ]
Millar, R.P. & Pawson, A.J., 2004. Outside in and inside out signalling: the new concept that selectivity of ligand binding at the Gonadotropin-Releasing Hormone receptor is modulated by the intracellular environment. Endocrinology 145, 3590-3593. [ Links ]
Millar, R.P., 2005. GnRHs and GnRH receptors. Anim. Reprod. Sci. 88, 5-28. [ Links ]
Naor, Z., 2009. Signaling by G-protein-coupled receptor (GPCR): Studies on the GnRH receptor. Fron. Neuroendocrinology 30, 10-29. [ Links ]
Noguera, J.L., Varona, L., Babot, D. & Estany, J., 2002. Multivariate analysis of litter size for multiple parities with production traits in pigs. I. Bayesian variance component estimation. J. Anim. Sci. 80, 2540-2547. [ Links ]
Ren, G., Chen, H., Zhang, L.Z., Lan, X.Y., Wei, T.B., Li, M.J., Jing, Y.J., Lei, C.Z. & Wang, J.Q., 2010. A coding SNP of LHX4 gene is associated with body weight and body length in bovine. Mol. Biol. Rep. 37, 417-422. [ Links ]
Rispoli, L.A. & Nett, T.M., 2005. Pituitary gonadotropin-releasing hormone (GnRH) receptor: structure, distribution and regulation of expression. Anim. Reprod. Sci. 88, 57-74. [ Links ]
Schubert, U., Anton, L.C., Gibbs, J., Norbury, C.C., Yewdell, J.W. & Bennink, J.R., 2000. Rapid degradation of a large fraction of newly synthesized proteins by proteasomes. Nature 404, 770-774. [ Links ]
Tang, K.Q., Li, S.J., Yang, W.C., Yu, J.N., Han, L., Li, X. & Yang, L.G., 2010. An MspI polymorphism in the inhibin alpha gene and its associations with superovulation traits in Chinese Holstein cows. Mol. Biol. Rep. 38, 17-21. [ Links ]
Vatankhah, M. & Talebi, M.A., 2008. Genetic parameters of body weight and fat-tail measurements in lambs. Small Rumin. Res. 75, 1-6. [ Links ]
Wu, W.S., Hua, G.H., Yang, L.G., Wen, Q.Y., Zhang, C.Y., Khairy, M.Z. & Chen, S.L., 2009. Association analysis of the INHA gene with litter size in Boer goats. Small Rumin. Res. 82, 139-143. [ Links ]
Yang, W.C., Li, S.J., Xie, Y.H., Tang, K.Q., Hua, G.H., Zhang, C.Y. & Yang, L.G., 2011. Two novel SNPs of the type I gonadotropin releasing hormone receptor gene and their associations with superovulation traits in Chinese Holstein cows. Livest. Sci. 136, 164-168. [ Links ]
Zhang, C.Y., Yang, L.G. & Shen, Z., 2008. Variance components and genetic parameters for weight and size at birth in the Boer goat. Livest. Sci. 115, 73-79. [ Links ]
Copyright resides with the authors in terms of the Creative Commons Attribution 2.5 South African Licence. See: http://creativecommons.org/licenses/by/2.5/za
Condition of use: The user may copy, distribute, transmit and adapt the work, but must recognise the authors and the South African Journal of Animal Science.
# Correspondence author: yangliguo2006@yahoo.com.cn
