FragGeneScan: predicting genes in short and error-prone reads.
The advances of next-generation sequencing technology have facilitated metagenomics research that attempts to determine directly the whole collection of genetic material within an environmental sample (i.e. the metagenome). Identification of genes directly from short reads has become an important yet challenging problem in annotating metagenomes, since the assembly of metagenomes is often not available. Gene predictors developed for whole genomes (e.g. Glimmer) and recently developed for metagenomic sequences (e.g. MetaGene) show a significant decrease in performance as the sequencing error rates increase, or as reads get shorter. We have developed a novel gene prediction method FragGeneScan, which combines sequencing error models and codon usages in a hidden Markov model to improve the prediction of protein-coding region in short reads. The performance of FragGeneScan was comparable to Glimmer and MetaGene for complete genomes. But for short reads, FragGeneScan consistently outperformed MetaGene (accuracy improved ∼62% for reads of 400 bases with 1% sequencing errors, and ∼18% for short reads of 100 bases that are error free). When applied to metagenomes, FragGeneScan recovered substantially more genes than MetaGene predicted (>90% of the genes identified by homology search), and many novel genes with no homologs in current protein sequence database.
Rho
Mina
M
School of Informatics and Computing, Indiana University, Bloomington, IN 47408, USA.
Tang
Haixu
H
Ye
Yuzhen
Y
D000465:Algorithms; D003062:Codon; D018965:Frameshifting, Ribosomal; D005796:Genes; D054892:Metagenome; D056186:Metagenomics / Q000379:methods*; D015233:Models, Statistical; D017422:Sequence Analysis, DNA; D012984:Software
D003062:Codon
Proc Natl Acad Sci U S A. 1974 Apr;71(4):1342-6
4598299
Curr Opin Microbiol. 2002 Jun;5(3):240-5
12057676
Genome Res. 2002 Apr;12(4):656-64
11932250
Proc Natl Acad Sci U S A. 2009 Apr 7;106(14):5725-30
19307581
BMC Genomics. 2009 Nov 12;10:520
19909532
DNA Res. 2008 Dec;15(6):387-96
18940874
Nucleic Acids Res. 2009 Jul;37(Web Server issue):W101-5
19429689
Genome Biol. 2009;10(8):R83
19682367
Nucleic Acids Res. 1998 Jan 15;26(2):544-8
9421513
Bioinformatics. 2006 Jul 15;22(14):e281-9
16873483
Nucleic Acids Res. 1999 Dec 1;27(23):4636-41
10556321
Nat Rev Genet. 2010 Jan;11(1):31-46
19997069
BMC Bioinformatics. 2008 Sep 19;9:386
18803844
Bioinformatics. 2009 Apr 1;25(7):962-3
19254921
Microbiol Mol Biol Rev. 2004 Dec;68(4):669-85
15590779
Genome Biol. 2002;3(2):REVIEWS0003
11864374
Annu Rev Genet. 2004;38:525-52
15568985
Nucleic Acids Res. 2002 Oct 1;30(19):4264-71
12364605
Proc Natl Acad Sci U S A. 2009 Sep 15;106(37):15527-33
19805210
DNA Res. 2009 Feb;16(1):1-12
19147530
Microbiol Mol Biol Rev. 2008 Dec;72(4):557-78, Table of Contents
19052320
Nature. 2009 Jan 22;457(7228):480-4
19043404
BMC Genomics. 2006 Jul 04;7:171
16820057
Appl Environ Microbiol. 2008 Mar;74(5):1453-63
18192407
BMC Genomics. 2008 Feb 08;9:75
18261238
Nucleic Acids Res. 2010 Jan;38(Database issue):D340-5
19892825
DNA Res. 2007 Aug 31;14(4):169-81
17916580
Nucleic Acids Res. 2010 Jul;38(12):e132
20403810
Nucleic Acids Res. 2006 Jul 1;34(Web Server issue):W588-90
16845076
Nucleic Acids Res. 1998 Feb 15;26(4):1107-15
9461475
Annu Rev Genomics Hum Genet. 2009;10:135-51
19715439
Nucleic Acids Res. 2010 Jan;38(Database issue):D468-72
19906710
Nucleic Acids Res. 2006;34(19):5623-30
17028096
BMC Bioinformatics. 2008 Mar 25;9:160
18366730
BMC Bioinformatics. 2008 Apr 10;9:182
18402669
BMC Bioinformatics. 2008 Apr 28;9:217
18442389
PLoS Comput Biol. 2006 May;2(5):e57
16710451
J Mol Biol. 1990 Oct 5;215(3):403-10
2231712
PLoS One. 2008 Oct 08;3(10):e3373
18841204
Science. 2006 Jun 2;312(5778):1355-9
16741115
Microbiol Rev. 1995 Mar;59(1):143-69
7535888
Annu Rev Microbiol. 2003;57:369-94
14527284