Association analysis between variants in KISS1 gene and litter size in goats
© An et al.; licensee BioMed Central Ltd. 2013
Received: 5 December 2012
Accepted: 25 July 2013
Published: 2 August 2013
Kisspeptins are the peptide products of KISS1 gene, which operate via the G - protein-coupled receptor GPR54. These peptides have emerged as essential upstream regulators of neurons secreting gonadotropin-releasing hormone (GnRH), the major hypothalamic node for the stimulatory control of the hypothalamic–pituitary– gonadal (HPG) axis. The present study detected the polymorphisms of caprine KISS1 gene in three goat breeds and investigated the associations between these genetic markers and litter size.
Three goat breeds (n = 680) were used to detect single nucleotide polymorphisms (SNPs) in the coding regions with their intron–exon boundaries and the proximal flanking regions of KISS1 gene by DNA sequencing and PCR–RFLP. Eleven novel SNPs (g.384G>A, g.1147T>C, g.1417G>A, g.1428_1429delG, g.2124C>T, g.2270C>T, g.2489T>C, g.2510G>A, g.2540C>T, g.3864_3865delCA and g.3885_3886insACCCC) were identified. It was shown that Xinong Saanen and Guanzhong goat breeds were in Hardy-Weinberg disequilibrium at g.384G>A locus (P < 0.05). Both g.2510G>A and g.2540C>T loci were closely linked in Xinong Saanen (SN), Guanzhong (GZ) and Boer (BG) goat breeds (r 2 > 0.33). The g.384G>A, g.2489T>C, g.2510G>A and g.2540C>T SNPs were associated with litter size (P<0.05). Individuals with AATTAATT combinative genotype of SN breed (SC) and TTAATT combinative genotype of BG breed (BC) had higher litter size than those with other combinative genotypes in average parity. The results extend the spectrum of genetic variation of the caprine KISS1 gene, which might contribute to goat genetic resources and breeding.
This study explored the genetic polymorphism of KISS1 gene, and indicated that four SNPs may play an important role in litter size. Their genetic mechanism of reproduction in goat breeds should be further investigated. The female goats with SC1 (AATTAATT) and BC7 (TTAATT) had higher litter size than those with other combinative genotypes in average parity and could be used for the development of new breeds of prolific goats. Further research on a large number of animals is required to confirm the link with increased prolificacy in goats.
KeywordsCombinative genotype SNP PCR-RFLP Candidate gene
Kisspeptins are the peptide products of KISS1 gene, which operate via the G - protein-coupled receptor GPR54 (also known as KISS1R). These peptides have emerged as essential upstream regulators of neurons secreting gonadotropin-releasing hormone (GnRH), the major hypothalamic node for the stimulatory control of the hypothalamic–pituitary– gonadal (HPG) axis . They are potent elicitors of gonadotropin secretion in various species and physiological settings. Moreover, KISS1 neurons in the hypothalamus participate in crucial features of reproductive maturation and function, such as brain-level sex differentiation, puberty onset and the neuroendocrine regulation of gonadotropin secretion and ovulation . Irwig et al. (2004) and Navarro et al. (2004) have provided evidences in rats that kisspeptin-expressing neurons are targets for regulation by sex steroids [3, 4], furthermore, these neurons are directly regulated by the negative and positive feedback actions of sex steroids in distinct regions of the forebrain .
Mutations of KISS1R are associated with hypogonadotrophic hypogonadism in humans [6, 7], a phenotype which is also observed in mice carrying inactivating mutations of KISS1 or KISS1R genes . In addition, to their prominent expression at hypothalamic levels, fragmentary evidences suggest that KISS1 and/or KISS1R mRNAs or proteins are also present in several peripheral reproductive tissues including the ovary [9, 10], oviduct  and testes . In humans, Pinto et al. (2012) reported kisspeptin modulated sperm progressive motility causing a biphasic (stimulatory and inhibitory) response and also induced transient sperm hyperactivation . One novel nonsynonymous single nucleotide polymorphism (G54650055T) substituting one amino acid in kisspeptin (P110T) was found to be statistically related to central precocious puberty (P<0.025) in humans . In sheep, KISS1 mRNA-expressing cells are found in the arcuate nucleus (ARC) and dorsallateral preoptic area and both appear to mediate the positive feedback effect of estradiol to generate the preovulatory GnRH/LH surge . The LH surge has been associated with an increase in the LH response to kisspeptin in humans and sheep [16, 17], indicating the surge may be generated by increased kisspeptin output and sensitivity. These findings indicate that KISS1 gene is an excellent candidate gene for reproductive traits in human and livestock.
Based on above considerations, here we detected the polymorphisms of caprine KISS1 gene in three goat breeds and investigated the associations between these genetic markers and litter size. This study provides some useful information on goat genetic resources and breeding.
SNPs identification and genotypes
Association analysis of SNPs with litter size
Least square means and standard errors of the litter size of SN and GZ breeds at g .384 G > A locus
1stparity litter size
2ndparity litter size
3rdparity litter size
4thparity litter size
Average litter size
Least square means and standard errors of the litter size of three breeds at g .2489 T > C , g .2510 G > A and g .2540 C > T loci
1stparity litter size
2ndparity litter size
3rdparity litter size
4thparity litter size
Average litter size
Effects of combinative genotypes on litter size
Association analysis of combinative genotypes of g.384G>A, g.2489T>C, g.2510G>A and g.2540C>T loci with litter size was done in SN and GZ breeds. In SN breed (Additional file 1: Table S4), individuals with SC3 (AATCGACT) had the greatest litter size in comparison with other combinative genotypes in the first parity. Individuals with SC1 (AATTAATT) and SC3 (AATCGACT) had greater litter size than those with SC4 (AATTGACT) and SC10 (GGTTAATT) in the second parity (P < 0.05). Individuals with SC1 (AATTAATT) had greater litter size than those with SC2 (AATTGGCC) and SC4 (AATTGACT) in the fourth parity (P < 0.05). Individuals with SC3 (AATCGACT) had greater litter size than those with SC4 (AATTGACT), SC5 (GATCGGCC) and SC13 (GATTGGCC) in average parity (P < 0.05). In GZ breed (Additional file 1: Table S5), Individuals with GC2 (GATCGGCC) had greater litter size than those with GC12 (GGTCGACT) and GC14 (GGTCGATT) in the first parity (P < 0.05). Individuals with GC14 (GGTCGATT) had the lowest litter size in comparison with other combinative genotypes in the third parity. Association analysis of combinative genotypes of g.2489T>C, g.2510G>A and g.2540C>T loci with litter size was done in BG breeds. In BG breed (Additional file 1: Table S6), individuals with BC7(TTAATT) had greater litter size than those with BC6(TTGACT) in the third, fourth and average parity (P < 0.05). Individuals with BC6 (TTGACT) had the lowest litter size in comparison with other combinative genotypes in average parity.
According to the classification of PIC (low polymorphism if PIC value < 0.25, moderate polymorphism if 0.25 < PIC value < 0.50, and high polymorphism if PIC > 0.50), SN and GZ breeds at g.384G>A locus had moderate genetic diversity, and SN, GZ and BG breeds had also moderate genetic diversity at g.2510G>A and g.2540C>T loci. The g.384G>A locus was in Hardy–Weinberg disequilibrium in SN and GZ breeds (P<0.05), which showed the genotypic frequencies were affected by selection, mutation or migration.
Identification of the candidate genes that are responsible for variation in continuous traits or quantitative traits has been a challenge in modern genetics. So far, there have been some studies of KISS1 gene as a candidate gene on reproductive traits in animals, which revealed that KISS1 gene plays an important role in animal reproduction [18–21]. Both c.374C>T and c.422C>G mutations were identified in human KISS1 gene, and the c.374C>T variant was associated with higher kisspeptin resistance to degradation in comparison with the wild type, suggesting a role for this mutation in the precocious puberty phenotype . Huijbregts et al. (2012) detected three SNPs (c.638insT, c.641C>G and c.645G>CA) in the 3′UTR of human KISS1 gene, and the c.645G>CA mutation was associated with central precocious puberty . In Jining Grey goat KISS1 gene, there are two mutations (G3433A and C3688A) in exon 3, three mutations (G296C, G454T and T505A) in intron 1 and an 18 bp deletion/insertion (1960–1977) in intron 2 and no mutations in exon 2 . Feng et al. (2009) detected the polymorphism in exon 2 of goat KISS1 gene and did not find polymorphism . In the current study, ten polymorphisms were detected in KISS1 gene of three goat breeds (g.1147T>C, g.1417G>A, g.1428_1429delG, g.2124C>T, g.2270C>T, g.2489T>C, g.2510G>A, g.2540C>T, g.3864_3865delCA and g.3885_3886insACCCC). The g.384G>A mutation was detected in SN and GZ breeds. Cao et al. (2010) indicated an association between allele C of the 296 locus and allele deletion of the 1960–1977 locus in KISS1 gene and high litter size in Jining Grey goats . Hou et al. (2011) identified T2643C and 8bp base deletions (2677AGTTCCCC) in the intron 2 of goat KISS1 gene and the T2643C had significant effects on litter size (P < 0.05) .
The reproductive traits are complex quantitative traits involving multiple genes, loci and interactions, so it is important to analyze the combined effect of multiple genes or loci on reproductive traits. In the study, association between multiple locus and litter size from the first to the fourth parity was analyzed. Mean litter size of goat tended to increase in later parity. Individuals with SC1 (AATTAATT) had higher litter size than those with SC4 (AATTGACT) and SC10 (GGTTAATT) in the second parity of SN breed. In addition, individuals with SC1 (AATTAATT) had higher litter size than those with SC10 (GGTTAATT) in average parity of SN breed. The litter size at second kidding is often a valuable index to determine whether a goat is prolific . Therefore, SC1 (AATTAATT) can be used in marker-assisted selection to select the individuals with higher litter size. Accumulating evidence further showed that central or peripheral administration of kisspeptin stimulates GnRH-dependent luteinizing hormone (LH) and follicle-stimulating hormone (FSH) secretion in various mammalian species from rodents to humans [27–29], suggesting that kisspeptin plays an essential role in governing reproductive functions throughout species. The biochemical and physiological functions, together with the results obtained in our study, indicate that KISS1 gene could be as a molecular breeding marker in goats.
This study explored the genetic polymorphism of KISS1 gene, and indicated that four SNPs may play an important role in litter size. Their genetic mechanism of reproduction in goat breeds should be further investigated. The female goats with SC1(AATTAATT) and BC7(TTAATT) had higher litter size than those with other combination genotypes in average parity and could be used for the development of new breeds of prolific goats. Further research on a large number of animals is required to confirm the link with increased prolificacy in goats.
Animals and genomic DNA isolation
Blood samples were obtained from 680 goats belonging to three breeds: Xinong Saanen goat (SN; n=257), Guanzhong goat (GZ, n=231) and Boer goat (BG; n=192). They were reared in Qianyang, Zhouzhi and Liuyou county of Shaanxi province, respectively. All diets were based on alfalfa, corn silage, and a combination of concentrates including corn, soya meal, and bone meal. Health, fertility and production records were maintained by the dairymen and veterinarians. The litter size from the first to fourth parity was obtained from production records. Five milliliters blood per goat were collected aseptically from the jugular vein and kept in a tube containing anticoagulant ACD (citric acid:sodium citrate:dextrose – 10: 27: 38). All samples were delivered back to the laboratory in an ice box. The genomic DNA was extracted from white blood cells using standard phenol-chloroform extraction protocol. All experiments were performed in accordance with the National Institute of Health Guide for the Care and Use of Laboratory Animals.
SNPs investigation and genotyping
According to caprine KISS1 gene (GenBank accession no. GU142847), Seven pairs of primers were designed to amplify caprine KISS1 gene. Their optimal annealing temperatures are showed in Additional file 1: Table S7. Herein we screened them for identifying SNPs of this gene by DNA pooling sequencing assay . Five microliters DNA (100ng/μl) per sample were collected to create a DNA pool for each goat breed. PCR products were sent to Beijing Genomics Institute (Beijing, China) to sequence in both directions. Discovery of SNPs was conducted using Chromas 2.31 and DNAstar 7.0 software.
The SNP in 5'UTR (89-409bp) of KISS1 gene was genotyped using primer–introduced restriction analysis–polymerase chain reaction (PIRA-PCR) . Other SNPs of KISS1 gene were genotyped with polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP). The 25 μL volume contained 50 ng genomic DNA, 12.5 μL 2 × reaction mix (including 500 μM dNTP each; 20 mM Tris–HCl; pH 9; 100 mM KCl; 3 mM MgCl2 ), 0.5 μM of each primer, and 0.5 units of Taq DNA polymerase. The cycling protocol was 5 min at 95°C, 35 cycles of denaturing at 94°C for 30 s, annealing at X°C (Additional file 1: Table S7) for 30 s, extending at 72°C for 35 s, with a final extension at 72°C for 10 min. PCR products (5μl) of different primer pairs were mixed with 0.7 μl 10 × Buffer, 2.5 U restriction enzyme (NEB, Ipswich, Britain) and 3.8 μl sterilized ddH2O, and then incubated for 1.5 h at 37°C. The restriction enzymes were showed in Additional file 1: Table S1. Digestion products were subjected to 3.5% horizontal agarose gel electrophoresis or 12% polyacrylamide gel electrophoresis (PAGE). The agarose and polyacrylamide gels were stained with ethidium bromide and 0.1% silver nitrate, respectively, and then the genotypes were observed.
The allelic frequencies, heterozygosity (He) and polymorphism information content (PIC) were calculated using Popgene (version 1.31). The linkage disequilibrium was performed by SHEsis software . Association analysis of combinative genotypes of different loci in KISS1 gene with litter size was done in three goat breeds. SCn ,GCn and BCn represented different combinative genotypes of SN, GZ and BG breeds, respectively. Statistical analysis was performed using univariate analysis in the general linear model procedure of SPSS 16 statistical software. Multiple comparisons of the means were performed using the least significant difference method or Dunnett’s T3. The model applied was: Yikm = μ + Gi + Sk+ Eikm, where Yikm is the trait measured on each of the ikmth animal, μ is the overall population mean, Gi is the fixed effect associated with ith genotype or combinative genotype, Sk is the fixed effect associated with the kth sire, and Eikm is the random error. Effects associated with farm, birth year and season of birth are not matched in the linear model, as the preliminary statistical analyses indicated that these effects did not have a significant influence on variability of traits in the analyzed populations.
This study was supported by the National Support Program of China (2011BAD28B05-3) and Science and Technology Innovation Project of Shaanxi Province (2011KTCL02-09).
- Navarro VM, Tena-Sempere M: Neuroendocrine control by kisspeptins: role in metabolic regulation of fertility. Nat Rev Endocrinol. 2012, 8: 40-53.View ArticleGoogle Scholar
- Caraty A, Franceschini I, Hoffman GE: Kisspeptin and the preovulatory gonadotrophin-releasing hormone/luteinising hormone surge in the ewe: basic aspects and potential applications in the control of ovulation. J Neuroendocrinol. 2010, 22: 710-715.PubMedGoogle Scholar
- Irwig MS, Fraleyb GS, Smith JT, Acohido BV, Popa SM, Cunningham MJ, Gottsch ML, Clifton DK, Steiner RA: Kisspeptin activation of gonadotropin releasing hormone neurons and regulation of KiSS-1 mRNA in the male rat. Neuroendocrinology. 2004, 80: 264-272. 10.1159/000083140.View ArticlePubMedGoogle Scholar
- Navarro VM, Castellano JM, Fernandez-Fernandez R, Barreiro ML, Roa J, Sanchez-Criado JE, Aguilar E, Dieguez C, Pinilla L, Tena-Sempere M: Developmental and hormonally regulated messenger ribonucleic acid expression of KiSS-1 and its putative receptor, GPR54, in rat hypothalamus and potent luteinizing hormone-releasing activity of KiSS-1 peptide. Endocrinology. 2004, 145: 4565-4574. 10.1210/en.2004-0413.View ArticlePubMedGoogle Scholar
- Gottsch M, Clifton D, Steiner R: Kisspeptin-GPR54 signaling in the neuroendocrine reproductive axis. Biol Reprod. 2006, 254–255: 91-96.Google Scholar
- de Roux N, Genin E, Carel JC, Matsuda F, Chaussain JL, Milgrom E: Hypogonadotropic hypogonadism due to loss of function of the KiSS1-derived peptide receptor GPR54. P Natl Acad Sci USA. 2003, 100: 10972-10976. 10.1073/pnas.1834399100.View ArticleGoogle Scholar
- Seminara SB, Messager S, Chatzidaki EE, Thresher RR, Acierno JS, Shagoury JK, Bo-Abbas Y, Kuohung W, Schwinof KM, Hendrick AG, et al: The GPR54 gene as a regulator of puberty. New Engl J Med. 2003, 349: 1614-U1618. 10.1056/NEJMoa035322.View ArticlePubMedGoogle Scholar
- Tena-Sempere M: Roles of kisspeptins in the control of hypothalamic-gonadotropic function: focus on sexual differentiation and puberty onset. Endocr Dev. 2010, 17: 52-62.View ArticlePubMedGoogle Scholar
- Castellano JM, Gaytan M, Roa J, Vigo E, Navarro VM, Bellido C, Dieguez C, Aguilar E, Sanchez-Criado JE, Pellicer A, et al: Expression of KiSS-1 in rat ovary: Putative local regulator ovulation?. Endocrinology. 2006, 147: 4852-4862. 10.1210/en.2006-0117.View ArticlePubMedGoogle Scholar
- Gaytan F, Gaytan M, Castellano JM, Romero M, Roa J, Aparicio B, Garrido N, Sanchez-Criado JE, Millar RP, Pellicer A, et al: KiSS-1 in the mammalian ovary: distribution of kisspeptin in human and marmoset and alterations in KiSS-1 mRNA levels in a rat model of ovulatory dysfunction. Am J Physiol-Endoc M. 2009, 296: E520-E531.Google Scholar
- Gaytan M, Castellano JM, Roa J, Sanchez-Criado JE, Tena-Sempere M, Gaytan F: Expression of KiSS-1 in rat oviduct: possible involvement in prevention of ectopic implantation?. Cell Tissue Res. 2007, 329: 571-579. 10.1007/s00441-007-0426-2.View ArticlePubMedGoogle Scholar
- Ohtaki T, Shintani Y, Honda S, Matsumoto H, Hori A, Kanehashi K, Terao Y, Kumano S, Takatsu Y, Masuda Y, et al: Metastasis suppressor gene KiSS-1 encodes peptide ligand of a G-protein-coupled receptor. Nature. 2001, 411: 613-617. 10.1038/35079135.View ArticlePubMedGoogle Scholar
- Pinto FM, Cejudo-Roman A, Ravina CG, Fernandez-Sanchez M, Martin-Lozano D, Illanes M, Tena-Sempere M, Candenas ML: Characterization of the kisspeptin system in human spermatozoa. Int J Androl. 2012, 35: 63-73. 10.1111/j.1365-2605.2011.01177.x.View ArticlePubMedGoogle Scholar
- Luan XH, Zhou YX, Wang W, Yu H, Li P, Gan XH, Wei DZ, Xiao JH: Association study of the polymorphisms in the KISS1 gene with central precocious puberty in Chinese girls. Eur J Endocrinol. 2007, 157: 113-118. 10.1530/EJE-07-0061.View ArticlePubMedGoogle Scholar
- Smith JT, Li Q, Yap KS, Shahab M, Roseweir AK, Millar RP, Clarke IJ: Kisspeptin is essential for the full preovulatory LH surge and stimulates GnRH release from the isolated ovine median eminence. Endocrinology. 2011, 152: 1001-1012. 10.1210/en.2010-1225.View ArticlePubMedGoogle Scholar
- Dhillo WS, Chaudhri OB, Thompson EL, Murphy KG, Patterson M, Ramachandran R, Nijher GK, Amber V, Kokkinos A, Donaldson M, et al: Kisspeptin-54 stimulates gonadotropin release most potently during the preovulatory phase of the menstrual cycle in women. J Clin Endocr Metab. 2007, 92: 3958-3966. 10.1210/jc.2007-1116.View ArticlePubMedGoogle Scholar
- Smith JT, Saleh SNH, Clarke IJ: Seasonal and cyclical change in the luteinizing hormone response to kisspeptin in the ewe. Neuroendocrinology. 2009, 90: 283-291. 10.1159/000227806.View ArticlePubMedGoogle Scholar
- Ardlie KG, Kruglyak L, Seielstad M: Patterns of linkage disequilibrium in the human genome. Nat Rev Genet. 2002, 3: 299-309. 10.1038/nrg777.View ArticlePubMedGoogle Scholar
- Shi YY, He L: SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005, 15: 97-98. 10.1038/sj.cr.7290272.View ArticlePubMedGoogle Scholar
- Huijbregts L, Roze C, Bonafe G, Houang M, Le Bouc Y, Carel JC, Leger J, Alberti P, de Roux N: DNA polymorphisms of the KiSS1 3 ' Untranslated region interfere with the folding of a G-rich sequence into G-quadruplex. Mol Cell Endocrinol. 2012, 351: 239-248. 10.1016/j.mce.2011.12.014.View ArticlePubMedGoogle Scholar
- Tomikawa J, Homma T, Tajima S, Shibata T, Inamoto Y, Takase K, Inoue N, Ohkura S, Uenoyama Y, Maeda K, et al: Molecular characterization and estrogen regulation of hypothalamic KISS1 gene in the pig. Biol Reprod. 2010, 82: 313-319. 10.1095/biolreprod.109.079863.View ArticlePubMedGoogle Scholar
- Silveira LG, Noel SD, Silveira-Neto AP, Abreu AP, Brito VN, Santos MG, Bianco SDC, Kuohung W, Xu S, Gryngarten M, et al: Mutations of the KISS1 gene in disorders of puberty. J Clin Endocr Metab. 2010, 95: 2276-2280. 10.1210/jc.2009-2421.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao GL, Chu MX, Fang L, Di R, Feng T, Li N: Analysis on DNA sequence of KiSS-1 gene and its association with litter size in goats. Mol Biol Rep. 2010, 37: 3921-3929. 10.1007/s11033-010-0049-7.View ArticlePubMedGoogle Scholar
- Feng T, Zhao YZ, Chu MX, Zhang YJ, Fang L, Di R, Cao GL, Li N: Association between sexual precocity and alleles of KISS-1 and GPR54 genes in goats. Anim Biotechnol. 2009, 20: 172-176. 10.1080/10495390903004493.View ArticlePubMedGoogle Scholar
- Hou JX, An XP, Wang JG, Song YX, Cui YH, Wang YF, Chen QJ, Cao BY: New genetic polymorphisms of KiSS-1 gene and their association with litter size in goats. Small Ruminant Res. 2011, 96: 106-110. 10.1016/j.smallrumres.2010.11.013.View ArticleGoogle Scholar
- Yuqin W, Li YX, Zhang NN, Wang ZB, Bai JY: Polymorphism of exon 2 of BMP15 gene and its relationship with litter size of two Chinese goats. Asian Austral J Anim. 2011, 24: 905-911. 10.5713/ajas.2011.10432.View ArticleGoogle Scholar
- Gottsch ML, Cunningham MJ, Smith JT, Popa SM, Acohido BV, Crowley WF, Seminara S, Clifton DK, Steiner RA: A role for kisspeptins in the regulation of gonadotropin secretion in the mouse. Endocrinology. 2004, 145: 4073-4077. 10.1210/en.2004-0431.View ArticlePubMedGoogle Scholar
- Navarro VM, Castellano JM, Fernandez-Fernandez R, Tovar S, Roa J, Mayen A, Nogueiras R, Vazquez MJ, Barreiro ML, Magni P, et al: Characterization of the potent luteinizing hormone-releasing activity of KiSS-1 peptide, the natural ligand of GPR54. Endocrinology. 2005, 146: 156-163.View ArticlePubMedGoogle Scholar
- Shahab M, Mastronardi C, Seminara SB, Crowley WF, Ojeda SR, Plant TM: Increased hypothalamic GPR54 signaling: a potential mechanism for initiation of puberty in primates. Proc Natl Acad Sci USA. 2005, 102: 2129-2134. 10.1073/pnas.0409822102.PubMed CentralView ArticlePubMedGoogle Scholar
- Bansal A, van den Boom D, Kammerer S, Honisch C, Adam G, Cantor CR, Kleyn P, Braun A: Association testing by DNA pooling: an effective initial screen. Proc Natl Acad Sci USA. 2002, 99: 16871-16874. 10.1073/pnas.262671399.PubMed CentralView ArticlePubMedGoogle Scholar
- Ke X, Collins A, Ye S: PIRA PCR designer for restriction analysis of single nucleotide polymorphisms. Bioinformatics. 2001, 17: 838-839. 10.1093/bioinformatics/17.9.838.View ArticlePubMedGoogle Scholar