High-throughput genotyping methods have increased the analytical power to study complex traits but high cost has remained a barrier for large scale use in animal improvement. We have adapted genotyping-by-sequencing (GBS) used in plants for genotyping 47 animals representing 7 taurine and indicine breeds of cattle from the US and Africa. Genomic DNA was digested with different enzymes, ligated to adapters containing one of 48 unique bar codes and sequenced by the Illumina HiSeq 2000. PstI was the best enzyme producing 1.4 million unique reads per animal and initially identifying a total of 63,697 SNPs. After removal of SNPs with call rates of less than 70%, 51,414 SNPs were detected throughout all autosomes with an average distance of 48.1 kb, and 1,143 SNPs on the X chromosome at an average distance of 130.3 kb, as well as 191 on unmapped contigs. If we consider only the SNPs with call rates of 90% and over, we identified 39,751 on autosomes, 850 on the X chromosome and 124 on unmapped contigs. Of these SNPs, 28,843 were not tightly linked to other SNPs. Average marker density per autosome was highly correlated with chromosome size (coefficient of correlation = ？0.798, r2 = 0.637) with higher density in smaller chromosomes. Average SNP call rate was 86.5% for all loci, with 53.0% of the loci having call rates >90% and the average minor allele frequency being 0.212. Average observed heterozygosity ranged from 0.046–0.294 among individuals, and from 0.064–0.197 among breeds, with Brangus showing the highest diversity as expected. GBS technique is novel, flexible, sufficiently high-throughput, and capable of providing acceptable marker density for genomic selection or genome-wide association studies at roughly one third of the cost of currently available genotyping technologies.
References
[1]
Pryce J, Hayes B (2012) A review of how dairy farmers can use and profit from genomic technologies. Animal Production Science 52: 180–184.
[2]
Andersson L (2001) Genetic dissection of phenotypic diversity in farm animals Nature Review Genetics. 2: 130–138.
[3]
Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, et al. (2009) Development and characterization of a high density SNP genotyping assay for cattle. PLOS ONE 4: e5350.
[4]
Snelling WM, Allan MF, Keele JW, Kuehn LA, McDaneld T, et al. (2010) Genome-wide association study of growth in crossbred beef cattle. Journal of Animal Science 88: 837–848.
[5]
Peters SO, Kizilkaya K, Garrick DJ, Fernando RL, Reecy JM, et al.. (2012) Bayesian genome wide association analyses of growth and yearling ultrasound measures of carcass traits in Brangus heifers. Journal of Animal Science.
Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, et al. (2011) Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nature Review Genetics 12: 499–510.
[8]
Bovine HapMap Consortium (2009) Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science 324: 528–532.
[9]
Zhan B, Fadista J, Thomsen B, Hedegaard J, Panitz F, et al. (2011) Global assessment of genomic variation in cattle by genome resequencing and high-throughput genotyping. BMC Genomics 12: 557.
[10]
Perez-Enciso M, Ferretti L (2010) Massive parallel sequencing in animal genetics: wherefroms and wheretos. Animal Genetics 41: 561–569.
[11]
Van Tassell CP, Smith TP, Matukumalli LK, Taylor JF, Schnabel RD, et al. (2008) SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nature Methods 5: 247–252.
[12]
Albrechtsen A, Nielsen FC, Nielsen R (2010) Ascertainment biases in SNP chips affect measures of population divergence. Molecular Biology and Evolution 27: 2534–2547.
[13]
Schaeffer LR (2006) Strategy for applying genome-wide selection in dairy cattle. Journal of Animal Breeding and Genetics 123: 218–223.
[14]
Pryce JE, Daetwyler HD (2011) Designing dairy cattle breeding schemes under genomic selection – a review of international research. Animal Production Science 52(3): 107–114.
[15]
Schefers JM, Weigel KA (2012) Genomic selection in dairy cattle: integration of DNA testing into breeding programs. Animal Frontiers 2: 4–9.
[16]
Pryce JE, Goddard ME, Raadsma HW, Hayes BJ (2010) Deterministic models of breeding scheme designs that incorporate genomic selection. Journal of Dairy Science 93: 5455–5466.
[17]
Van Eenennaam AL, van der Werf JHJ, Goddard ME (2011) The value of using DNA markers for beef bull selection in the seed stock sector. Journal of Animal Science 89: 307–320.
[18]
Boichard D, Chung H, Dassonneville R, David X, Eggen A, et al. (2012) Design of a bovine low-density SNP array optimized for imputation. PLOS ONE 7: e34130.
[19]
Hawken RJ, Zhang YD, Fortes MR, Collis E, Barris WC, et al. (2012) Genome-wide association studies of female reproduction in tropically adapted beef cattle. Journal of Animal Science 90: 1398–1410.
[20]
Parameswaran P, Jalili R, Tao L, Shokralla S, Gharizadeh B, et al. (2007) A pyrosequencing-tailored nucleotide barcode design unveils opportunities for large-scale sample multiplexing. Nucleic Acids Research 35: e130.
[21]
Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, et al. (2008) Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLOS ONE 3: e3376.
[22]
Sansaloni C, Petroli C, Jaccoud D, Carling J, Detering F, et al. (2011) Diversity Arrays Technology (DArT) and next-generation sequencing combined: genome-wide, high throughput, highly informative genotyping for molecular breeding of Eucalyptus.. BMC Proceedings 5(S7): P54.
[23]
van Orsouw NJ, Hogers RCJ, Janssen A, Yalcin F, Snoeijers S, et al. (2007) Complexity Reduction of Polymorphic Sequences (CRoPS？): A novel approach for large-scale polymorphism discovery in complex genomes. PLOS ONE 2: e1172.
[24]
Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, et al. (2011) A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLOS ONE 6: e19379.
[25]
Gore MA, Chia JM, Elshire RJ, Sun Q, Ersoz ES, et al. (2009) A first-generation haplotype map of maize. Science 326: 1115–1117.
[26]
Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754–1760.
[27]
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Molecular Biology and Evolution 28: 2731–2739.
[28]
Seroussi E, Glick G, Shirak A, Yakobson E, Weller JI, et al. (2010) Analysis of copy loss and gain variations in Holstein cattle autosomes using BeadChip SNPs. BMC Genomics 11: 673.
[29]
Cole JB, VanRaden PM, O'Connell JR, Van Tassell CP, Sonstegard TS, et al. (2009) Distribution and location of genetic effects for dairy traits. Journal of Dairy Science 92: 2931–2946.
[30]
Sharma RK, Romay C, Harriman J, Elshire RJ, Glaubitz J, et al.. (2011) The degree of ascertainment bias in SSRs, SNPs and genotyping-bysequencing (GBS). Abstract presented in the 53rd Annual Maize Genetics Conference, St. Charles, Illinois, March 17–20.
[31]
Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, et al. (2006) Global variation in copy number in the human genome. Nature 444: 444–454.
[32]
Hou Y, Liu GE, Bickhart DM, Matukumalli LK, Li C, et al. (2012) Genomic regions showing copy number variations associate with resistance or susceptibility to gastrointestinal nematodes in Angus cattle. Functional and Integrative Genomics 12: 81–92.
[33]
Ku CS, Loy EY, Salim A, Pawitan Y, Chia KS (2010) The discovery of human genetic variations and their use as disease markers: past, present and future. American Journal of Human Genetics 55: 403–415.
[34]
Barris W, Harrison BE, McWilliam S, Bunch RJ, Goddard ME, et al. (2012) Next generation sequencing of African and indicine cattle to identify single nucleotide polymorphisms. Animal Production Science 52: 133–142.
[35]
Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, et al. (2009) The Sorghum bicolor genome and the diversification of grasses. Nature 457: 551–556.
[36]
Bovine Genome Sequencing, Analysis Consortium (2009) The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science 324: 522–528.
[37]
Qanbari S, Pimentel EC, Tetens J, Thaller G, Lichtner P, et al. (2010) The pattern of linkage disequilibrium in German Holstein cattle. Animal Genetics 41(4): 346–56.
[38]
Poland JA, Brown PJ, Sorrells ME, Jannink J-L (2012) Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLOS ONE 7(2): e32253.
