Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success

Humble, Emily; Thorne, Michael A.S. ORCID: https://orcid.org/0000-0001-7759-612X; Forcada, Jaume ORCID: https://orcid.org/0000-0002-2115-0150; Hoffman, Joseph I.. 2016 Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success. BMC Research Notes, 9 (1), 418. 12, pp. https://doi.org/10.1186/s13104-016-2209-x

Before downloading, please read NORA policies.

Preview

Text
Humble.pdf - Published Version
Available under License Creative Commons Attribution 4.0.
Download (1MB) | Preview

Official URL: http://dx.doi.org/10.1186/s13104-016-2209-x

Abstract/Summary

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.

Item Type:	Publication - Article
Digital Object Identifier (DOI):	https://doi.org/10.1186/s13104-016-2209-x
Programmes:	BAS Programmes > BAS Programmes 2015 > Ecosystems
ISSN:	1756-0500
Additional Keywords:	transcriptome, Roche 454 sequencing, Illumina HiSeq sequencing, single nucleotide polymorphism, validation success, marine mammal, Antarctic fur seal, Arctocephalus gazella
Date made live:	07 Sep 2016 10:19 +0 (UTC)
URI:	https://nora.nerc.ac.uk/id/eprint/514399

Item Type:

Publication - Article

Digital Object Identifier (DOI):

https://doi.org/10.1186/s13104-016-2209-x

Programmes:

BAS Programmes > BAS Programmes 2015 > Ecosystems

ISSN:

1756-0500

Additional Keywords:

transcriptome, Roche 454 sequencing, Illumina HiSeq sequencing, single nucleotide polymorphism, validation success, marine mammal, Antarctic fur seal, Arctocephalus gazella

Date made live:

07 Sep 2016 10:19 +0 (UTC)

URI:

https://nora.nerc.ac.uk/id/eprint/514399