title: Application of massive parallel sequencing to whole genome SNP discovery in the porcine genome.
abstract: BACKGROUND
  Although the Illumina 1 G Genome Analyzer generates billions of base pairs of sequence data, challenges arise in sequence selection due to the varying sequence quality. Therefore, in the framework of the International Porcine SNP Chip Consortium, this pilot study aimed to evaluate the impact of the quality level of the sequenced bases on mapping quality and identification of true SNPs on a large scale.
  
  
  RESULTS
  DNA pooled from five animals from a commercial boar line was digested with DraI; 150-250-bp fragments were isolated and end-sequenced using the Illumina 1 G Genome Analyzer, yielding 70,348,064 sequences 36-bp long. Rules were developed to select sequences, which were then aligned to unique positions in a reference genome. Sequences were selected based on quality, and three thresholds of sequence quality (SQ) were compared. The highest threshold of SQ allowed identification of a larger number of SNPs (17,489), distributed widely across the pig genome. In total, 3,142 SNPs were validated with a success rate of 96%. The correlation between estimated minor allele frequency (MAF) and genotyped MAF was moderate, and SNPs were highly polymorphic in other pig breeds. Lowering the SQ threshold and maintaining the same criteria for SNP identification resulted in the discovery of fewer SNPs (16,768), of which 259 were not identified using higher SQ levels. Validation of SNPs found exclusively in the lower SQ threshold had a success rate of 94% and a low correlation between estimated MAF and genotyped MAF. Base change analysis suggested that the rate of transitions in the pig genome is likely to be similar to that observed in humans. Chromosome X showed reduced nucleotide diversity relative to autosomes, as observed for other species.
  
  
  CONCLUSIONS
  Large numbers of SNPs can be identified reliably by creating strict rules for sequence selection, which simultaneously decreases sequence ambiguity. Selection of sequences using a higher SQ threshold leads to more reliable identification of SNPs. Lower SQ thresholds can be used to guarantee sufficient sequence coverage, resulting in high success rate but less reliable MAF estimation. Nucleotide diversity varies between porcine chromosomes, with the X chromosome showing less variation as observed in other species.
authors.affiliation: Animal Breeding and Genomics Centre, Wageningen University, Wageningen 6700 AH, The Netherlands. andreia.fonseca@wur.nl
authors.forename: Andreia J
authors.initials: AJ
authors.lastname: Amaral
authors.forename: Hendrik-Jan
authors.initials: HJ
authors.lastname: Megens
authors.forename: Hindrik H D
authors.initials: HH
authors.lastname: Kerstens
authors.forename: Henri C M
authors.initials: HC
authors.lastname: Heuven
authors.forename: Bert
authors.initials: B
authors.lastname: Dibbits
authors.forename: Richard P M A
authors.initials: RP
authors.lastname: Crooijmans
authors.forename: Johan T
authors.initials: JT
authors.lastname: den Dunnen
authors.forename: Martien A M
authors.initials: MA
authors.lastname: Groenen
mesh_terms: D000465:Algorithms; D000818:Animals; D002874:Chromosome Mapping / Q000379:methods; D033481:Chromosomes, Mammalian / Q000235:genetics; D016678:Genome; D015698:Genomic Library; D005838:Genotype; D008297:Male; D010865:Pilot Projects; D020641:Polymorphism, Single Nucleotide; D016415:Sequence Alignment; D017422:Sequence Analysis, DNA / Q000379:methods*; D034421:Sus scrofa / Q000235:genetics*
references.citation: Genet Sel Evol. 2008 Jan-Feb;40(1):103-28
references.pmid: 18096118
references.citation: Science. 1998 May 15;280(5366):1077-82
references.pmid: 9582121
references.citation: Genetics. 2002 Jan;160(1):247-56
references.pmid: 11805060
references.citation: Genetics. 2007 Apr;175(4):1911-21
references.pmid: 17287527
references.citation: Curr Opin Genet Dev. 2006 Dec;16(6):545-52
references.pmid: 17055251
references.citation: Nat Rev Genet. 2002 Nov;3(11):862-71
references.pmid: 12415316
references.citation: Genetics. 2008 May;179(1):569-79
references.pmid: 18493072
references.citation: Genetics. 2004 Dec;168(4):2261-9
references.pmid: 15611190
references.citation: Science. 2001 Feb 16;291(5507):1304-51
references.pmid: 11181995
references.citation: Gene. 2002 Oct 30;300(1-2):31-42
references.pmid: 12468083
references.citation: Nat Methods. 2008 Mar;5(3):247-52
references.pmid: 18297082
references.citation: Nat Methods. 2008 Feb;5(2):183-8
references.pmid: 18204455
references.citation: Theor Popul Biol. 1975 Apr;7(2):256-76
references.pmid: 1145509
references.citation: Proc Natl Acad Sci U S A. 2005 Mar 15;102(11):4063-7
references.pmid: 15728731
references.citation: Genetics. 2001 Oct;159(2):907-11
references.pmid: 11693127
references.citation: Nucleic Acids Res. 2005 Nov 07;33(19):e171
references.pmid: 16275781
references.citation: BMC Genomics. 2009 Jan 06;10:4
references.pmid: 19126189
references.citation: Anim Biotechnol. 2005;16(1):41-54
references.pmid: 28881173
references.citation: Genome Res. 2008 Nov;18(11):1851-8
references.pmid: 18714091
references.citation: Nature. 2000 Sep 28;407(6803):513-6
references.pmid: 11029002
references.citation: BMC Genet. 2008 Dec 04;9:81
references.pmid: 19055830
references.citation: Nat Rev Genet. 2004 May;5(5):335-44
references.pmid: 15143316
references.citation: Genome Res. 2005 Dec;15(12):1767-76
references.pmid: 16339375
references.citation: Comp Funct Genomics. 2005;6(4):251-5
references.pmid: 18629187
references.citation: Genome Res. 1998 Mar;8(3):175-85
references.pmid: 9521921
references.citation: Clin Chim Acta. 2006 Jan;363(1-2):83-94
references.pmid: 16165119