currently , the uniprotkb protein sequence database contains > 17 million protein sequences ( 1 ) .
this wealth of data is helping us to understand biology at an ever increasing rate .
a large fraction of these sequences can be grouped into a few thousand common protein families .
proteins within these families often share common functions that can allow information experimentally gleaned on one protein to be transferred to uncharacterized ones .
this process of transitive annotation is essential to make sense of the rapidly growing amount of sequence data .
there are concerns about transitive annotation not being robust and thus leading to numerous annotation errors ( 2 ) .
although this phenomenon does occur it seems clear that high - quality manual curation of the protein sequence databases , the careful use of databases of protein families for annotation and feedback from users of protein databases have largely kept the gross errors in check .
for example , incorrect protein function assignments from large - scale genome projects in general have not been transferred to hundreds or thousands of other proteins as feared . on the other hand ,
subtler misannotations such as assigning an incorrect but related enzymatic activity to a protein ( for example phosphorylating the wrong substrate ) occur . due to the lack of experimental work on most proteins
a recent estimate for six large enzyme superfamilies studied suggested a range of 563% of incorrect annotations ( 3 ) .
a further source of error in the sequence databases is the prediction of spurious genes ( 4 ) .
automatic gene prediction methods in prokaryotes are increasingly accurate , glimmer3 , for example , both improves start site prediction relative to glimmer2 and reduces the high false - positive rate for high gc genomes .
however , given the large number of proteins being deposited in the sequence databases , it is still likely that many thousands of the included sequences are either wholly spurious or improperly extended , past their true start sites . as the capacity to manually curate gene predictions diminishes ,
it has been noted that certain alternate reading frames seem more likely to give rise to long spurious open reading frames ( orfs ) ( 6 ) .
normark et al . ( 7 ) found that frames + 3 and 1 were most likely to give rise to long spurious orfs .
although alternative overlapping reading frames are used in viral genomes , there are relatively few confirmed cases found in prokaryotes or eukaryotes . during the construction of the pfam database of protein families ( 8) ,
we have occasionally been alerted to the presence of families that were entirely composed of spuriously predicted orfs .
once one gene has been spuriously predicted and put in the sequence database , there is a danger that future genome projects will annotate new protein - coding genes by similarity to the first spurious orf .
this can lead to entire families of spurious orfs . in the worst - case scenario ,
( 9 ) where a spuriously predicted gene on the opposite strand of ribosomal rna had been given the incorrect function of cell wall hydrolase ( pf10695 ) .
it may seem surprising that spuriously predicted orfs would appear to have conservation like bona fide proteins .
however , at the protein level the alignment of spurious orfs can look like a normal protein alignment . in figure 1
, we show the multiple sequence alignment for former pfam family pf10695 , showing a protein - like conservation pattern .
this conservation is actually due to the selective forces conserving the opposite strand rrna sequence and structure .
once these errors are propagated to pfam and other databases , then there is a danger that the error will be widely transferred and hence difficult to correct .
the first example ( figure 2a ) shows a pair of proteins with correctly identified homology domains but with an uncharacteristically long tail - to - tail overlap .
the second is an example of a hidden markov model ( hmm)-based domain definition identifying a region in a spurious gene call that overlaps a true gene ( figure 2b ) . 
 figure 1.seed alignment for the antifam family derived from pf10695 .
amino acids are colored by average similarity according to the blosum62 amino acid substitution matrix from most similar ( light blue ) to less similar ( gray ) .
 s and  e in the first row stand for sequence start and sequence end , respectively .
( a ) shows two proteins from the corynebacterium efficiens genome that encode components of a restriction system .
( b ) two highly overlapping predicted proteins from the rhodopirellula baltica genome coded on opposite strands of dna . the q7uy10 protein contains two pfam duf1596 domains .
green boxes represent regions matched by pfam families , the red shaded areas represent transmembrane domains predicted by phobius ( 10 ) and the blue shaded areas represent regions of low complexity ( 11 ) .
amino acids are colored by average similarity according to the blosum62 amino acid substitution matrix from most similar ( light blue ) to less similar ( gray ) . 
s and  e in the first row stand for sequence start and sequence end , respectively .
( a ) shows two proteins from the corynebacterium efficiens genome that encode components of a restriction system .
( b ) two highly overlapping predicted proteins from the rhodopirellula baltica genome coded on opposite strands of dna .
green boxes represent regions matched by pfam families , the red shaded areas represent transmembrane domains predicted by phobius ( 10 ) and the blue shaded areas represent regions of low complexity ( 11 ) .
antifam matches to predicted proteins in some cases will suggest that modifications to the extent of the coding region are needed rather than complete deletion of the protein from the sequence database .
most prfb genes , encoding the bacterial translation release factor 2 , have a + 1 programmed frameshift early in the coding region ( 12 ) .
the region downstream of the frameshift site is easily identified by gene finders , but unreconstructed extension 5 to the frameshift results in translation of the wrong reading frame .
some of the multiple sequence alignments used are taken from the pfam database seed alignments for families identified as spurious orfs .
these alignments are kept as they appeared in the final release of pfam before they were withdrawn ( table 1 ) .
several additional custom families have been created to identify other commonly recurring spurious orfs ( table 2 ) .
the profile hmms have been constructed using the hmmer3 package with default parameters ( 13 ) .
the profile hmm library can be searched against any set of protein sequences using the 
due to the speed of the hmmer3 package , searching a sequence database such as uniprotkb will take a few minutes and searching a complete proteome will take seconds . 

table 1.antifam entries derived from pfam familiespfam accession number ( identifier)last pfam release presentreason for deleting from pfamno . of matches in uniprotno . of matches in metagenomics data
evidence for homology to known protein on opposite strand30pf07616 ( duf1578)15.0proteins may not be expressed .
evidence for homology to known protein on opposite strand66pf07630 ( duf1591)15.0proteins may not be expressed .
evidence for homology to known protein on opposite strand60pf07633 ( duf1594)15.0proteins may not be expressed .
evidence for homology to known protein on opposite strand50pf11370 ( duf3170)25.0probable shadow orf of clp protease167pf11194 ( duf2825)25.0probable crispr repeat regions15918pf11664 ( duf3264)25.0probable crispr repeat regions2113pf10695 ( cw - hydrolase)25.0antisense to rrna ( 9)2251,654pf10919 ( duf2699)26.0shadow orf of pf00665 ( integrase core domain 1)2511pf07641 ( duf1596)26.0dubious genome annotation .
family comprises only three sequences from rhodopirellula baltica , two overlapping30the final two columns show the number of matches of each antifam entry to uniprotkb and to a metagenomic data set.the metagenomic set of sequences is the same as that used by pfam ( 14).crispr , clustered regularly interspaced short palindromic repeats . 

table 2.antifam entries derived from custom multiple sequence alignmentidentifiertype of spurious familyno . of matches in uniprotno . of matches in metagenomics data setspurious_orf_10translated bacterial trna , trna01196795spurious_orf_11translated bacterial trna , trna0289170spurious_orf_12translated bacterial trna , trna03143408spurious_orf_13translated bacterial trna , trna0477671spurious_orf_14translated bacterial trna , trna05156191spurious_orf_15translated bacterial trna , trna063163spurious_orf_16translated bacterial trna , trna074017spurious_orf_17translated bacterial trna , trna08510spurious_orf_18translated bacterial trna , trna09439spurious_orf_19translated bacterial trna , trna10712spurious_orf_20translated bacterial trna , trna114328spurious_orf_21prfb frameshift245spurious_orf_22from a lncrna , linc00174261the metagenomic set of sequences is the same as that used by pfam ( 14 ) .
antifam entries derived from pfam families the final two columns show the number of matches of each antifam entry to uniprotkb and to a metagenomic data set .
the metagenomic set of sequences is the same as that used by pfam ( 14 ) .
antifam entries derived from custom multiple sequence alignment the metagenomic set of sequences is the same as that used by pfam ( 14 ) .
antifam is primarily a tool that is aimed at bioinformaticians to be used as part of genome annotation projects .
the antifam alignments and profile hmms can be downloaded from the following url : ftp://ftp.sanger.ac.uk/pub/databases/pfam/antifam/ of the 1310 proteins identified in uniprotkb as probably being spurious the large majority were from trembl , the unreviewed part of uniprotkb .
they had been automatically created from the records in the european nucleotide archive , genbank or dna data bank of japan ( ddbj ) .
one spurious protein found in the reviewed swiss - prot section of uniprotkb was y114_chlmu ( q9pli5 ) that is an uncharacterized protein from chlamydia muridarum .
this belonged to the previously mentioned spurious cw - hydrolase family and was removed in uniprot release 2011_10 .
an additional 13 spurious proteins in the reviewed portion of uniprotkb are also identified , of which 8 are due to non - coding rna translations : 
 o67358.1 aquifex aeolicus trigger factor contains frameshift extension;p19773.1 mycobacterium tuberculosis protein matching duf2699;p47080.1 yeast protein yjl007c product of a dubious gene prediction;p92540.1 arabidopsis protein;q04100.1 yeast protein ydr445c product of dubious gene prediction and partly overlaps ydr444w;q52m62.3 human product of a dubious coding sequence ( cds ) prediction .
o67358.1 aquifex aeolicus trigger factor contains frameshift extension ; p19773.1 mycobacterium tuberculosis protein matching duf2699 ; p47080.1 yeast protein yjl007c product of a dubious gene prediction ; p92540.1 arabidopsis protein ; q04100.1 yeast protein ydr445c product of dubious gene prediction and partly overlaps ydr444w ; q52m62.3 human product of a dubious coding sequence ( cds ) prediction .
probable non - coding rna ; q9cjr2.1 pasteurella multocida trna - derived match ; q9cmd0.1 p. multocida trna - derived match ; q9cmx0.1 p. multocida trna - derived match ; and q9cmz6.1 p.
in addition to the families reported by pfam users , we tried to identify if further spurious families existed .
the large majority of proteins in the trembl portion of uniprotkb come from translations found in entries in the european nucleotide archive , genbank or ddbj .
thus , we scanned trembl entries to identify uniprotkb entries that overlapped with each other in the nucleotide entry .
we confined our scan to the prokaryotic entries because the nature of overlaps is relatively simple compared to the complex patterns of interlacing and nesting found in eukaryotic gene structures .
this list of proteins was then used to identify further pfam families that contained numerous overlapping genes .
this measure means that large well - known families that are likely to have many overlaps by chance are not at the top of the list .
the first release of antifam contains only a modest number of families . however , we see a number of ways to increase this in the future .
the first of these is to increase the number of non - coding rna - based families .
we currently have only one ribosomal rna - based family and we can add many further families .
we can identify proteins related to ribosomal rnas initially using tblastn , which compares a protein to a nucleotide sequence considering all six reading frames .
in addition , we could also consider comparing a large database of rna sequences to the protein sequence databases to identify further potentially spurious proteins . to date , we have only been able to investigate the pfam families with the highest fraction of overlapping proteins .
but in the coming months , we will investigate this list more thoroughly to identify if any further pfam families should be deleted and added to antifam .
the first release of antifam contains 23 families derived from pfam as well as a small number of non - coding rnas that were erroneously translated into protein sequences .
we expect that this number will grow in the future and we have several ideas to help us to achieve this .
this should increase the power of antifam to reduce the number of spurious orfs finding their way into the sequence databases .
we hope that antifam will become an indispensible tool for quality control in metagenomic and genomic studies .
we are particularly keen for biocurators and experimental biologists to remain vigilant and alert us to new cases of spurious orfs so that we can add them to this resource .
wellcome trust ( grant number wt077044/z/05/z ) ; national human genome research institute ( grant number r01 hg004881 ) . funding for open access charge : wellcome trust ( grant number wt077044/z/05/z ) .
 as the deluge of genomic dna sequence grows the fraction of protein sequences that have been manually curated falls . in turn , as the number of laboratories with the ability to sequence genomes in a high - throughput manner grows , the informatics capability of those labs to accurately identify and annotate all genes within a genome may often be lacking . 
 these issues have led to fears about transitive annotation errors making sequence databases less reliable . during the lifetime of the pfam protein families 
 database a number of protein families have been built , which were later identified as composed solely of spurious open reading frames ( orfs ) either on the opposite strand or in a different , overlapping reading frame with respect to the true protein - coding or non - coding rna gene . 
 these families were deleted and are no longer available in pfam . 
 however , we realized that these may perform a useful function to identify new spurious orfs . 
 we have collected these families together in antifam along with additional custom - made families of spurious orfs . 
 this resource currently contains 23 families that identified 1310 spurious proteins in uniprotkb and a further 4119 spurious proteins in a collection of metagenomic sequences . 
 uniprot has adopted antifam as a part of the uniprotkb quality control process and will investigate these spurious proteins for exclusion . 
Introduction
Description of the resource
Identification of problematic Pfam families
Future plans
Conclusions
Supplementary Data
Funding