title: Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success.
abstract: BACKGROUND
  Single nucleotide polymorphism (SNP) discovery is an important goal of many studies. However, the number of 'putative' SNPs discovered from a sequence resource may not provide a reliable indication of the number that will successfully validate with a given genotyping technology. For this it may be necessary to account for factors such as the method used for SNP discovery and the type of sequence data from which it originates, suitability of the SNP flanking sequences for probe design, and genomic context. To explore the relative importance of these and other factors, we used Illumina sequencing to augment an existing Roche 454 transcriptome assembly for the Antarctic fur seal (Arctocephalus gazella). We then mapped the raw Illumina reads to the new hybrid transcriptome using BWA and BOWTIE2 before calling SNPs with GATK. The resulting markers were pooled with two existing sets of SNPs called from the original 454 assembly using NEWBLER and SWAP454. Finally, we explored the extent to which SNPs discovered using these four methods overlapped and predicted the corresponding validation outcomes for both Illumina Infinium iSelect HD and Affymetrix Axiom arrays.
  
  
  RESULTS
  Collating markers across all discovery methods resulted in a global list of 34,718 SNPs. However, concordance between the methods was surprisingly poor, with only 51.0 % of SNPs being discovered by more than one method and 13.5 % being called from both the 454 and Illumina datasets. Using a predictive modeling approach, we could also show that SNPs called from the Illumina data were on average more likely to successfully validate, as were SNPs called by more than one method. Above and beyond this pattern, predicted validation outcomes were also consistently better for Affymetrix Axiom arrays.
  
  
  CONCLUSIONS
  Our results suggest that focusing on SNPs called by more than one method could potentially improve validation outcomes. They also highlight possible differences between alternative genotyping technologies that could be explored in future studies of non-model organisms.
authors: Department of Animal Behaviour, University of Bielefeld, Postfach 100131, 33501, Bielefeld, Germany. emily.humble@uni-bielefeld.de.|British Antarctic Survey, High Cross, Madingley Road, Cambridge, CB3 OET, UK. emily.humble@uni-bielefeld.de.
authors: Emily
authors: http://orcid.org/0000-0001-6841-1269
authors: E
authors: Humble
authors: British Antarctic Survey, High Cross, Madingley Road, Cambridge, CB3 OET, UK.
authors: Michael A S
authors: MA
authors: Thorne
authors: British Antarctic Survey, High Cross, Madingley Road, Cambridge, CB3 OET, UK.
authors: Jaume
authors: J
authors: Forcada
authors: Department of Animal Behaviour, University of Bielefeld, Postfach 100131, 33501, Bielefeld, Germany.
authors: Joseph I
authors: JI
authors: Hoffman
mesh_terms: D005787:Gene Frequency / Q000235:genetics; D060005:Genotyping Techniques / Q000379:methods*; D006801:Humans; D058977:Molecular Sequence Annotation; D020641:Polymorphism, Single Nucleotide / Q000235:genetics*; D011336:Probability; D015203:Reproducibility of Results; D017422:Sequence Analysis, DNA / Q000379:methods*; D059467:Transcriptome / Q000235:genetics*
keywords: Antarctic fur seal; Arctocephalus gazella; Illumina HiSeq sequencing; Marine mammal; Roche 454 sequencing; Single nucleotide polymorphism; Transcriptome; Validation success