Annotation-based genome-wide SNP discovery in the large and complex Aegilops tauschii genome using next-generation sequencing without a reference genome sequence.
BACKGROUND
Many plants have large and complex genomes with an abundance of repeated sequences. Many plants are also polyploid. Both of these attributes typify the genome architecture in the tribe Triticeae, whose members include economically important wheat, rye and barley. Large genome sizes, an abundance of repeated sequences, and polyploidy present challenges to genome-wide SNP discovery using next-generation sequencing (NGS) of total genomic DNA by making alignment and clustering of short reads generated by the NGS platforms difficult, particularly in the absence of a reference genome sequence.


RESULTS
An annotation-based, genome-wide SNP discovery pipeline is reported using NGS data for large and complex genomes without a reference genome sequence. Roche 454 shotgun reads with low genome coverage of one genotype are annotated in order to distinguish single-copy sequences and repeat junctions from repetitive sequences and sequences shared by paralogous genes. Multiple genome equivalents of shotgun reads of another genotype generated with SOLiD or Solexa are then mapped to the annotated Roche 454 reads to identify putative SNPs. A pipeline program package, AGSNP, was developed and used for genome-wide SNP discovery in Aegilops tauschii-the diploid source of the wheat D genome, and with a genome size of 4.02 Gb, of which 90% is repetitive sequences. Genomic DNA of Ae. tauschii accession AL8/78 was sequenced with the Roche 454 NGS platform. Genomic DNA and cDNA of Ae. tauschii accession AS75 was sequenced primarily with SOLiD, although some Solexa and Roche 454 genomic sequences were also generated. A total of 195,631 putative SNPs were discovered in gene sequences, 155,580 putative SNPs were discovered in uncharacterized single-copy regions, and another 145,907 putative SNPs were discovered in repeat junctions. These SNPs were dispersed across the entire Ae. tauschii genome. To assess the false positive SNP discovery rate, DNA containing putative SNPs was amplified by PCR from AL8/78 and AS75 and resequenced with the ABI 3730 xl. In a sample of 302 randomly selected putative SNPs, 84.0% in gene regions, 88.0% in repeat junctions, and 81.3% in uncharacterized regions were validated.


CONCLUSIONS
An annotation-based genome-wide SNP discovery pipeline for NGS platforms was developed. The pipeline is suitable for SNP discovery in genomic libraries of complex genomes and does not require a reference genome sequence. The pipeline is applicable to all current NGS platforms, provided that at least one such platform generates relatively long reads. The pipeline package, AGSNP, and the discovered 497,118 Ae. tauschii SNPs can be accessed at (http://avena.pw.usda.gov/wheatD/agsnp.shtml).
Department of Plant Sciences, University of California, Davis, CA 95616, USA. jdvorak@ucdavis.edu
Frank M
FM
You
Naxin
N
Huo
Karin R
KR
Deal
Yong Q
YQ
Gu
Ming-Cheng
MC
Luo
Patrick E
PE
McGuire
Jan
J
Dvorak
Olin D
OD
Anderson
D018745:Genome, Plant / Q000235:genetics*; D058977:Molecular Sequence Annotation / Q000379:methods*; D006109:Poaceae / Q000235:genetics*; D020641:Polymorphism, Single Nucleotide / Q000235:genetics*; D017422:Sequence Analysis, DNA / Q000379:methods*
ISME J. 2009 Nov;3(11):1314-7
19587772
Methods Mol Biol. 2007;406:137-59
18287691
Science. 2009 Nov 20;326(5956):1115-7
19965431
Bioinformatics. 2010 Mar 1;26(5):589-95
20080505
J Hered. 1946 Mar;37:81 107
20985728
PLoS Comput Biol. 2009 May;5(5):e1000386
19461883
Genome Biol. 2007;8(7):R143
17659080
Genome Biol. 2009;10(3):R32
19327155
Nucleic Acids Res. 2006 Jan 1;34(Database issue):D741-4
16381971
Plant J. 2006 Nov;48(3):463-74
17010109
BMC Bioinformatics. 2010 Mar 15;11:130
20230626
Nature. 2000 Sep 28;407(6803):513-6
11029002
Genome. 2006 Sep;49(9):1131-9
17110993
Plant Biotechnol J. 2010 Feb;8(2):196-210
20078842
Nat Methods. 2008 Mar;5(3):247-52
18297082
BMC Genomics. 2008 Jun 30;9:312
18590545
Bioinformatics. 2009 Jul 15;25(14):1754-60
19451168
Proc Natl Acad Sci U S A. 2009 Sep 15;106(37):15780-5
19717446
BMC Bioinformatics. 2008 May 29;9:253
18510760
BMC Bioinformatics. 2010 Apr 13;11:187
20388221
Nucleic Acids Res. 2004 Jan 1;32(Database issue):D360-3
14681434
BMC Genomics. 2010 Aug 16;11:475
20712881
Theor Appl Genet. 1992 Jul;84(3-4):419-29
24203203
Nucleic Acids Res. 2010 Jul;38(Web Server issue):W313-20
20497996
Nat Genet. 1999 Dec;23(4):452-6
10581034
Genome. 2003 Jun;46(3):431-7
12834059
Plant J. 2007 Sep;51(5):910-8
17662031
Plant J. 2004 Nov;40(4):500-11
15500466
Plant Biotechnol J. 2009 May;7(4):334-46
19207216
Genome Res. 2008 May;18(5):821-9
18349386
Genome. 2009 Jun;52(6):576-87
19483776
Nucleic Acids Res. 2008 Jan;36(Database issue):D1028-33
18089549
BMC Genomics. 2010 Dec 14;11:702
21156062
Bioinformatics. 2009 Aug 15;25(16):2078-9
19505943
Cold Spring Harb Protoc. 2010 Apr;2010(4):pdb.prot5409
20360363
Genetics. 2007 May;176(1):685-96
17339218
BMC Genomics. 2010 Jan 15;11:38
20078886