Automated protein subfamily identification and classification.
Function prediction by homology is widely used to provide preliminary functional annotations for genes for which experimental evidence of function is unavailable or limited. This approach has been shown to be prone to systematic error, including percolation of annotation errors through sequence databases. Phylogenomic analysis avoids these errors in function prediction but has been difficult to automate for high-throughput application. To address this limitation, we present a computationally efficient pipeline for phylogenomic classification of proteins. This pipeline uses the SCI-PHY (Subfamily Classification in Phylogenomics) algorithm for automatic subfamily identification, followed by subfamily hidden Markov model (HMM) construction. A simple and computationally efficient scoring scheme using family and subfamily HMMs enables classification of novel sequences to protein families and subfamilies. Sequences representing entirely novel subfamilies are differentiated from those that can be classified to subfamilies in the input training set using logistic regression. Subfamily HMM parameters are estimated using an information-sharing protocol, enabling subfamilies containing even a single sequence to benefit from conservation patterns defining the family as a whole or in related subfamilies. SCI-PHY subfamilies correspond closely to functional subtypes defined by experts and to conserved clades found by phylogenetic analysis. Extensive comparisons of subfamily and family HMM performances show that subfamily HMMs dramatically improve the separation between homologous and non-homologous proteins in sequence database searches. Subfamily HMMs also provide extremely high specificity of classification and can be used to predict entirely novel subtypes. The SCI-PHY Web server at http://phylogenomics.berkeley.edu/SCI-PHY/ allows users to upload a multiple sequence alignment for subfamily identification and subfamily HMM construction. Biologists wishing to provide their own subfamily definitions can do so. Source code is available on the Web page. The Berkeley Phylogenomics Group PhyloFacts resource contains pre-calculated subfamily predictions and subfamily HMMs for more than 40,000 protein families and domains at http://phylogenomics.berkeley.edu/phylofacts/.
Brown
Duncan P
DP
Department of Bioengineering, University of California, Berkeley, California, United States of America.
Krishnamurthy
Nandini
N
Sjölander
Kimmen
K
D000465:Algorithms; D000595:Amino Acid Sequence; D001185:Artificial Intelligence; D008390:Markov Chains; D008969:Molecular Sequence Data; D010363:Pattern Recognition, Automated / Q000379:methods*; D011506:Proteins / Q000737:chemistry* / Q000145:classification*; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D016415:Sequence Alignment / Q000379:methods*; D020539:Sequence Analysis, Protein / Q000379:methods*
D011506:Proteins
Proc Natl Acad Sci U S A. 1994 Dec 6;91(25):12091-5
7991589
BMC Bioinformatics. 2005 Aug 05;6:198
16083500
Proteins. 2003;53 Suppl 6:491-6
14579338
Bioinformatics. 1998;14(10):846-56
9927713
Proteins. 2000 Oct 1;41(1):98-107
10944397
Bioinformatics. 1999 May;15(5):391-412
10366660
Methods Enzymol. 1996;266:460-80
8743700
Proteins. 2005;61 Suppl 7:201-213
16187363
Math Biosci. 2005 Feb;193(2):223-34
15748731
J Mol Biol. 1998 Dec 11;284(4):1201-10
9837738
BMC Bioinformatics. 2002 May 16;3:14
12028595
J Mol Biol. 1995 Apr 7;247(4):536-40
7723011
Bioinformatics. 2006 Jul 1;22(13):1658-9
16731699
Nucleic Acids Res. 2005 Jan 1;33(Database issue):D192-6
15608175
J Mol Biol. 1994 Nov 4;243(4):574-8
7966282
Nucleic Acids Res. 2003 Jul 1;31(13):3829-32
12824430
J Mol Biol. 2001 Dec 14;314(5):1041-52
11743721
Science. 2004 Apr 2;304(5667):66-74
15001713
Nucleic Acids Res. 2000 Jan 1;28(1):254-6
10592239
Nucleic Acids Res. 2004 Jan 1;32(Database issue):D138-41
14681378
J Mol Evol. 2001 Jun;52(6):540-2
11443357
Nucleic Acids Res. 1995 Jul 25;23(14):2715-23
7651832
Nucleic Acids Res. 2001 Jan 1;29(1):346-9
11125133
Bioinformatics. 2005 Aug 1;21(15):3241-7
15905279
Pac Symp Biocomput. 2005;:322-33
15759638
PLoS Comput Biol. 2005 Oct;1(5):e45
16217548
Bioinformatics. 2003 Dec 12;19(18):2484-5
14668239
J Mol Biol. 2000 Oct 13;303(1):61-76
11021970
Trends Genet. 2001 Aug;17(8):429-31
11485799
Nucleic Acids Res. 2003 Jul 1;31(13):3799-803
12824422
Arch Immunol Ther Exp (Warsz). 2002;50(3):159-67
12098931
Bioinformatics. 2004 Jan 22;20(2):170-9
14734307
Nucleic Acids Res. 2004 Jul 1;32(Web Server issue):W313-7
15215401
Bioinformatics. 2002 Jul;18(7):908-21
12117788
Biochemistry. 2006 Feb 28;45(8):2545-55
16489747
Annu Rev Genet. 2004;38:525-52
15568985
Trends Genet. 1999 Apr;15(4):132-3
10203816
Nucleic Acids Res. 2004 Sep 30;32(17):5231-8
15459293
Bioinformatics. 2002 Jan;18(1):77-82
11836214
Nucleic Acids Res. 2004 Jan 1;32(Database issue):D189-92
14681391
Nucleic Acids Res. 2007 Jul;35(Web Server issue):W27-32
17488835
In Silico Biol. 1998;1(1):55-67
11471243
Proc Natl Acad Sci U S A. 1987 Jul;84(13):4355-8
3474607
Nucleic Acids Res. 2003 Jan 1;31(1):294-7
12520006
Proteins. 2003 Nov 15;53(3):683-92
14579359
Genome Biol. 2006;7(9):R83
16973001
J Mol Biol. 1994 Feb 4;235(5):1501-31
8107089
Biochem Biophys Res Commun. 2001 Apr 6;282(3):832-8
11401539
Nucleic Acids Res. 2004 Jan 1;32(Database issue):D115-9
14681372
PLoS Comput Biol. 2006 Jun 30;2(6):e77
16846248
Genome Res. 2003 Sep;13(9):2178-89
12952885
Methods Biochem Anal. 2003;44:501-7
12647401
Comput Appl Biosci. 1996 Aug;12(4):327-45
8902360
Proc Int Conf Intell Syst Mol Biol. 1998;6:165-74
9783222
Curr Opin Chem Biol. 1999 Oct;3(5):548-56
10508675
Science. 2001 Feb 16;291(5507):1304-51
11181995
Bioinformatics. 1998;14(9):755-63
9918945
BMC Bioinformatics. 2004 Nov 18;5:178
15550167
Bioinformatics. 2002 Jan;18(1):147-59
11836223
Genome Res. 1998 Mar;8(3):163-7
9521918
Bioinformatics. 2002 Jan;18(1):92-9
11836216
Nucleic Acids Res. 1998 Sep 15;26(18):4291-300
9722651
Mol Biol Evol. 2001 Aug;18(8):1435-41
11470834
Nucleic Acids Res. 2004 Mar 19;32(5):1792-7
15034147
Brief Bioinform. 2006 Sep;7(3):225-42
16772267