- Open Access
Computational identification of strain-, species- and genus-specific proteins
BMC Bioinformatics volume 6, Article number: 279 (2005)
The identification of unique proteins at different taxonomic levels has both scientific and practical value. Strain-, species- and genus-specific proteins can provide insight into the criteria that define an organism and its relationship with close relatives. Such proteins can also serve as taxon-specific diagnostic targets.
A pipeline using a combination of computational and manual analyses of BLAST results was developed to identify strain-, species-, and genus-specific proteins and to catalog the closest sequenced relative for each protein in a proteome. Proteins encoded by a given strain are preliminarily considered to be unique if BLAST, using a comprehensive protein database, fails to retrieve (with an e-value better than 0.001) any protein not encoded by the query strain, species or genus (for strain-, species- and genus-specific proteins respectively), or if BLAST, using the best hit as the query (reverse BLAST), does not retrieve the initial query protein. Results are manually inspected for homology if the initial query is retrieved in the reverse BLAST but is not the best hit. Sequences unlikely to retrieve homologs using the default BLOSUM62 matrix (usually short sequences) are re-tested using the PAM30 matrix, thereby increasing the number of retrieved homologs and increasing the stringency of the search for unique proteins. The above protocol was used to examine several food- and water-borne pathogens. We find that the reverse BLAST step filters out about 22% of proteins with homologs that would otherwise be considered unique at the genus and species levels. Analysis of the annotations of unique proteins reveals that many are remnants of prophage proteins, or may be involved in virulence. The data generated from this study can be accessed and further evaluated from the CUPID (C ore and U nique P rotein Id entification) system web site (updated semi-annually) at http://pir.georgetown.edu/cupid.
CUPID provides a set of proteins specific to a genus, species or a strain, and identifies the most closely related organism.
Over 200 pathogenic and non-pathogenic bacteria have been completely sequenced , including multiple strains from several species. The availability of sequence data from related genomes has facilitated comparative genomic analysis, which not only allows the study of major evolutionary processes, but also the determination of proteins conserved across — or unique to — different species [2–4]. Information on the presence or absence of genes is a powerful tool to gain knowledge about the metabolism, pathogenicity, physiology and behavior of different organisms [2–4]. It can also provide the basis for the detection of pathogens in a given sample and to distinguish between pathogenic and non-pathogenic relatives — critical to combating disease and to the emerging biodefense field.
Two recent studies (ORFanage  and Procom ) have focused on species- and clade-specific genes, both noting a substantial number of unique genes encoded in specific organisms. The ORFanage database includes about 32,000 unique ORFs (ORFans) from 84 fully sequenced microbial genomes. An ORFan is a protein that failed to hit any other protein encoded within those 84 genomes with an E-value better than 10e-3 (or 10e-5 for alignments of <80 residues) using BLAST . The Procom database compared proteins from thirty completely sequenced eukaryotic genomes. The proteins were pair-wise compared using WU-BLASTP  with a threshold E-value of 1.
We developed a general protocol to identify proteins unique to different taxa, and applied it to a set of food- and water-borne pathogens. Specifically, we: a) identify proteins that are unique to a particular strain, species, or genus; b) extract the set of proteins common to two or more strains or species; and c) determine the organism most closely related to a particular genus, species or strain. The information generated from this study is available from the CUPID (C ore and U nique P rotein Id entification) system via a flexible and easy-to-use web interface.
Construction and content
Definition of terms
Throughout this paper the following terms apply:
The leaf-most taxonomic node indicated for a given fully-sequenced genome. The node can be at the species, strain, or sub-strain taxonomic level.
The source strain, species or genus of a BLAST query when considering strain-, species-, and genus-specific proteins, respectively. For example, a protein encoded by any Bacillus is considered self at the genus level when the query organism is Bacillus anthracis strain Ames.
The following terms refer to the protein sets generated by this study:
Proteins that are present in all selected organisms. "Selected organisms" will include all the completely-sequenced genomes in the data set that are taxonomically identical to the query organism at the strain, species or genus level.
Proteins that have no related sequences in non-self organisms. Related sequences are mined from all sources (not just from completely-sequenced genomes). Related sequences from the query organism (paralogs) are ignored as self.
Proteins that have related sequences in all selected organisms, but not in non-selected organisms.
The following terms refer to the evidence provided for the indicated uniqueness status on the results page. Unless otherwise indicated, the evidence is based on BLAST results:
The query retrieves a protein from a non-self organism with an e-value better than 0.001.
The query sequence is identical to a sequence from a non-self organism.
The best non-self hit has the original query as its best non-self hit. Reciprocal best hits represent an approximation of orthology .
No external hits
The query fails to hit a protein from a non-self organism using 0.001 as a cutoff.
No reverse hit
The best hit does not retrieve the original query when testing for reciprocal hits in a reverse BLAST. Additional evidence is provided for such cases. Forward gives the E-value for the hit between the query and subject. Maxrev gives the worst E-value reported when the subject is used as query.
The best hit retrieves the original query in a reverse BLAST, but it is not the reciprocal best hit. The E-values for the forward and reverse BLASTS are given.
"Unclear" cases that were resolved after manual checking of the data.
The results shown are for BLAST using the PAM30 matrix. The default matrix is BLOSUM62.
Data source and algorithm
Complete proteomes were retrieved from the NREF database at PIR , a comprehensive (~2.1 million sequences) non-redundant protein database. Identical sequences from a given species are merged in this database. Sequence searches were conducted using the BLAST program. An initial set of organism-specific proteins was generated computationally. Proteins that could not be confirmed as unique by computer analysis were manually checked by several methods , including multiple sequence alignment , hidden Markov models , and PSI-BLAST . Proteins specific to a set of organisms were identified by determining which of the conserved core proteins [13, 14] are unique (using the same procedure as used to determine organism-specific proteins). A flowchart of the process is shown in Figure 1.
Certain sequences failed to retrieve themselves using the BLOSUM62 matrix. It is known that the PAM matrices may increase the information content in the alignment and hence optimize the alignment score . The PAM matrices were evaluated to identify the optimal word size, gap existence cost and gap extension cost values to be used. BLAST results using the different PAM matrices and different parameters were manually checked to identify which sequences should be analyzed using a different matrix (data not shown). We found that sequences that could not retrieve themselves with an expect value of 1e-14 or better in the forward BLAST using the BLOSUM62 matrix (usually short sequences <30 aa) should be queried again using the PAM30 matrix (word size = 2, gap existence cost = 9 and gap extension cost = 1) to improve the ability to retrieve homologs. The inclusion of the reverse BLAST step and optimized search parameters both yield an analysis method that is skewed toward tagging proteins as not unique whenever possible.
The information generated for each organism at each taxonomic level (genus, species, or strain) is stored in a flat-file. Data for each protein is given on a single machine-parsable line, and is tracked using a unique protein identifier to allow integration with other PIR resources. The line contains the protein ID, an indication of whether homologs were found (+ for yes, ! for no), evidence for homology or lack thereof, identity and source organism for any core protein(s) found, and an indication of whether the query protein is part of the core. The results can be queried and browsed through the CUPID web interface. Additional analysis of the data is facilitated by links to PIR protein resources  and other bioinformatics servers.
Unique proteins for common food- and water-borne pathogens will be re-determined every six months using pre-computed BLAST (generated quarterly). De novo re-computation is necessary due to the constant increase in (and possible change to) protein sequences. We will add selected NIAID category pathogens with each semi-annual update cycle. In addition, we provide a mechanism to request additional proteome analyses. Users may select from a list of organisms, accessible from the web interface, that is populated by proteomes with stable sequence annotation in the UniProt Knowledgebase (UniProtKB) . Analysis time will depend on several factors, including proteome size and evolutionary distance from other proteomes. The user will be notified upon completion of the analysis. Results from user-selected proteomes will be made available from the CUPID web site. Such analyses will be stamped with the date of the last computation, but will not be part of the regular update cycle.
Usage and utility
Retrieving strain-, species-, or genus-specific proteins
Unique proteins are retrieved by clicking an organism listed on the CUPID homepage (Figure 2), followed by clicking on the taxonomic level of interest. By default, the initial retrieval returns only core unique proteins (Figure 3). To retrieve just the unique proteins, one should select the taxonomic level of interest and "Unique" from the pull-down menus. For example, selecting "Genus" level "Unique" proteins for E. coli O157:H7 will display all proteins in E. coli O157:H7 that may or may not be present in other Escherichia, but are not present in any other genus. The "Unique" column will either have a plus sign (+), be left blank, or have a question mark (?). A plus sign means that the protein is unique by "no external hits", by no reverse BLAST hit, or by curator judgment. Blank means that the protein is not unique, and a question mark means that the relationship is of unclear status (that is, it cannot be unambiguously determined by computational means alone, and human judgment is required). The core status column indicates if the given protein has homologs in all proteomes considered self, while the core hits column lists the best hit protein from each of the self proteomes. The user can save the retrieved list in tab-delimited format or save the protein sequences in fasta format.
Retrieving the set of proteins specific to two or more organisms
Proteins unique to two or more organisms can be retrieved from the species- or genus-specific proteins results page. This is done, in effect, by redefining the proteomes to be considered for the core. For example, once the results for E. coli O157:H7 are displayed, one can display the proteins unique to E. coli O157:H7 and E. coli O157:H7 EDL933 by checking the box for "Core Hits must *only* be from" and the box for Tax ID 155864. This will retrieve proteins that are conserved in both E. coli O157:H7 and E. coli O157:H7 EDL933, but are not found anywhere else. Alternatively, selecting "Show Genus level All proteins" with the same checkboxes selected as above will display all the proteins conserved between E. coli O157:H7 and E. coli O157:H7 EDL933 strains – not just the unique ones – but will exclude those found in other Escherichia proteomes.
Retrieving all the proteins from a particular organism
All the proteins from an organism can be retrieved from any results page by selecting "Show XX level All proteins" (where XX can be Genus, Species, or Strain).
Determining the closest relative
The genus, species, or strain that is most closely related to the selected organism based on the best BLAST hits of its entire proteome can be identified by clicking on the TaxID for the top non-self hits on the results page. For example, from the result page of genus-specific proteins of E. coli O157:H7, the top non-Escherichia hits are to Shigella flexneri, and from the result page of strain-specific proteins for E. coli O157:H7 the top non-Escherichia coli O157:H7 hits are to E. coli O157:H7 EDL933 (not shown).
Description of selected results
We identified the unique proteins in thirty organisms (predominantly food- or water-borne bacterial pathogens and their close relatives). The results for ten common water and wastewater strains of Escherichia, Salmonella and Helicobacte r are described below. E. coli is the primary indicator organism for biologically contaminated food and water. Several strains of E. coli, such as the laboratory strain K12 , are completely harmless, while E. coli CFT073 is an extra-intestinal uropathogenic bacterium and the O157:H7 strains are enterohemorrhagic pathogens. S. enterica Typhi is a human pathogen that causes enteric typhoid fever. The Ty2 strain of S. enterica was isolated before the advent of antibiotics and has no plasmids, whereas the CT18 strain harbors a multiple-drug-resistance plasmid and a cryptic plasmid. S. typhimurium LT2 is a strain of S. enterica (synonym: S. enterica subsp. enterica serovar Typhimurium strain LT2). H. pylori is associated with peptic ulcers and several types of gastric cancer. H. hepaticus causes chronic hepatitis and is also a recognized carcinogen. A summary of computer-assisted manual analyses of these organisms is shown in Table 1.
The number of unique proteins in each of four E. coli strains (K12, CFT073, O157:H7 and O157:H7 EDL933) did not show any apparent correlation with the pathogenic nature of the organism. For example, E. coli CFT073 (pathogenic) has only three unique proteins compared to 35 in K12 (non-pathogenic) and 117 and 73 in O157:H7 (pathogenic) and O157:H7 EDL933 (pathogenic) strains, respectively. High numbers of unique sequences in some organisms could be due to the unavailability of a closely related non-self genome or, as is likely in this example, unique plasmids that may be present in the query genome.
The number of genus/species-specific proteins for each E. coli is not identical as one might first expect; rather, they range from 12 to 15. This phenomenon can also be seen in other organisms where there are multiple genomes from the same species or genus. Three explanations are possible: i) each strain has different numbers of paralogs; ii) a frameshift produces two annotated open reading frames in one strain where another has one; iii) the comprehensive dataset used contains largely duplicate but non-identical entries.
The number of unique proteins for several E. coli strains drops as one compares strain level to genus level (except for the CFT073 strain), a trend that is also apparent for the Helicobacter strains. However, the opposite is true for the Salmonella. The progression is not necessarily linear or consistent. For example, the number of unique proteins from the CT18 strain of Salmonella enterica first drops from 38 at the strain level to 7 at the species level before again rising to 101 at the genus level.
A majority of the unique proteins detected in this study lacked meaningful functional annotation (i.e., were annotated simply as "hypothetical protein") (Table 2). The available annotation for the remainder indicates a preponderance of proteins with some relationship to pathogenesis or virulence (22%), or derived from phages (25%). In addition, a combination of annotation and subcellular localization prediction using PSORT-B  indicates that 26% of the proteins are external to the cell (cell wall attached or secreted).
It is evident from this study that a major reason that few unique proteins are found in some cases is the presence of sequence data for closely-related organisms and, by extension, the peculiarities of taxonomic designations. Therefore, only a general trend in the number of strain-, species-, or genus-specific proteins can be established. Strains with closely related sequenced strains tend to have relatively few unique proteins at that level while the converse is true for those without close relatives (compare Helicobacter pylori strains with Helicobacter hepaticus). The trend also holds true at the genus level (compare Escherichia or Helicobacter genus-specific proteins with Salmonella). We note that a proteome from a closely-related genus is represented in the protein database for both Escherichia (Shigella) and Helicobacter (Campylobacter). Shigella is so similar to E. coli that there are recommendations to consider them different species within the same genus , while Helicobacter pylori was once Campylobacter pylori.
The number of apparently unique genes encoded in specific organisms can depend on the definition of "unique" and the parameters and underlying database used to identify these proteins. For example, CUPID identifies 110 proteins unique to H. pylori 26695 at the species level, whereas the ORFanage database lists 260 proteins. The ORFanage database considers the protein HP0052 to be unique. BLAST with HP0052 (UniProt ID: O24893_HELPY) retrieves several non-Helicobacter proteins (Moraxella nonliquefaciens - E = 9e-22, Ehrlichia canis - E = 3e-18, Mycoplasma mycoides - E = 2e-15). None of these non-self hits have a complete genome sequence deposited in the database, and therefore were not considered during the construction of the ORFanage. In addition, other lists of unique proteins consider a protein to be unique if it did not produce a BLAST hit with the E value cut-off set to 10-3 [20, 21], and do not further assess cases that have lower similarity. However, the contribution of low-similarity homologs is significant. For example, 950 proteins from the organisms presented in Table 1 that were considered not unique at the genus level were considered so because of a reciprocal best hit. Overall, we found that about 22% of proteins were removed from the list of unique proteins based on the reciprocal best hit criterion. In general, many sequences are considered to be singleton ORFs by other studies but not considered unique in this study because: a) a limited protein set was used for BLAST comparison as opposed to a comprehensive database; b) the definition of self is different (strain differences may not be considered); c) only unidirectional BLAST hits (without reverse BLAST) are considered to find homologs; and d) the cutoff used to define homolog is less conservative. Using a thorough computational method and a comprehensive database leads to a more conservative estimate of the number of unique proteins.
Despite the conservative approach used here, one must be mindful of certain pitfalls in deriving lists of unique proteins. First, a protein might be labeled as unique only because homologs from other organisms were missed upon submission of the sequence, or because of some other conceptual translation problem. In all cases, the short list of potential unique proteins should be further screened computationally at the DNA level using tBLASTn. Running tBLASTn using a protein of interest will make sure that the gene is indeed unique, at least with respect to the current pool of submitted sequences (to aid in this task, we have included a link to NCBI's tBLASTn web page). Second, many of the proteins identified here are remnants of prophage proteins . The leading role of bacteriophages in shaping the E. coli O157:H7 genome is evident by the presence of 24 prophages and prophage-like elements, and several genes that have been laterally transferred . At the strain level, labeling of such proteins as unique may be more a reflection of a gap in whole-genome sequence information than of true specificity. Accordingly, discrimination between individual strains (isolates) may require laboratory comparison methods such as pulsed-field gel electrophoresis or whole-cell fatty acid analysis. In contrast, identifying species- or genus-specific proteins can be done with confidence when multiple representatives have been sequenced. In such cases, conservation within multiple strains of a species (for example) gives confidence in the "reality" of the uniqueness because that status has been conserved over time (core unique proteins).
Precise identification of pathogens is important so that adequate action can be taken to either eliminate or reduce the threat of infection. One use of the CUPID system is to help identify diagnostic targets specific to a particular clade of these pathogens. The unique proteins form a short list of diagnostic targets to be validated in the laboratory. Proteins predicted to be external to the bacterial cell – possibly involved in host interactions and virulence – may be used to develop protein-based detection systems. In addition, it should be possible to use the DNA encoding these proteins as the basis for diagnostics. However, it is important to note that protein-identified DNA probes must be verified as to their uniqueness at the DNA level.
The salient features of CUPID are: a) provides sets of proteins unique to a strain, species, and genus level; b) includes a check for additional homologs based on reciprocal hits; c) uses different parameters for short sequences; d) provides the identity of the nearest non-self neighbor; and e) allows retrieval of unique, core, and core unique proteins at different taxonomic levels.
Availability and requirements
CUPID is freely accessible from the PIR website at http://pir.georgetown.edu/cupid.
C ore and U nique P rotein Id entification system
P rotein I nformation R esource
Bernal A, Ear U, Kyrpides N: Genomes OnLine Database (GOLD): a monitor of genome projects world-wide. Nucleic Acids Res 2001, 29: 126–127. 10.1093/nar/29.1.126
Eisen JA, Fraser CM: Phylogenomics: intersection of evolution and genomics. Science 2003, 300: 1706–1707. 10.1126/science.1086292
Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families. Science 1997, 278: 631–637. 10.1126/science.278.5338.631
Siew N, Fischer D: Analysis of singleton ORFans in fully sequenced microbial genomes. Proteins 2003, 53: 241–251. 10.1002/prot.10423
Siew N, Azaria Y, Fischer D: The ORFanage: an ORFan database. Nucleic Acids Res 2004, 32 Database issue: D281–3. 10.1093/nar/gkh116
Li JB, Zhang M, Dutcher SK, Stormo GD: Procom: a web-based tool to compare multiple eukaryotic proteomes. Bioinformatics 2005, 21: 1693–1694. 10.1093/bioinformatics/bti161
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25: 3389–3402. 10.1093/nar/25.17.3389
Gish W: WU-BLASTP http://blast.wustl.edu. 1996.
Wu CH, Yeh LS, Huang H, Arminski L, Castro-Alvear J, Chen Y, Hu Z, Kourtesis P, Ledley RS, Suzek BE, Vinayaka CR, Zhang J, Barker WC: The Protein Information Resource. Nucleic Acids Res 2003, 31: 345–347. 10.1093/nar/gkg040
Mazumder R, Iyer LM, Vasudevan S, Aravind L: Detection of novel members, structure-function analysis and evolutionary classification of the 2H phosphoesterase superfamily. Nucleic Acids Res 2002, 30: 5229–5243. 10.1093/nar/gkf645
Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol 2000, 302: 205–217. 10.1006/jmbi.2000.4042
Eddy SR: Profile hidden Markov models. Bioinformatics 1998, 14: 755–763. 10.1093/bioinformatics/14.9.755
Zafar N, Mazumder R, Seto D: CoreGenes: a computational tool for identifying and cataloging "core" genes in a set of small genomes. BMC Bioinformatics 2002, 3: 12. 10.1186/1471-2105-3-12
Mazumder R, Kolaskar A, Seto D: GeneOrder: comparing the order of genes in small genomes. Bioinformatics 2001, 17: 162–166. 10.1093/bioinformatics/17.2.162
Altschul SF: A protein alignment scoring system sensitive at all evolutionary distances. J Mol Evol 1993, 36: 290–300. 10.1007/BF00160485
Bairoch A, Apweiler R, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LS: The Universal Protein Resource (UniProt). Nucleic Acids Res 2005, 33 Database Issue: D154–9.
Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF, Gregor J, Davis NW, Kirkpatrick HA, Goeden MA, Rose DJ, Mau B, Shao Y: The complete genome sequence of Escherichia coli K-12. Science 1997, 277: 1453–1474. 10.1126/science.277.5331.1453
Gardy JL, Spencer C, Wang K, Ester M, Tusnady GE, Simon I, Hua S, deFays K, Lambert C, Nakai K, Brinkman FS: PSORT-B: Improving protein subcellular localization prediction for Gram-negative bacteria. Nucleic Acids Res 2003, 31: 3613–3617. 10.1093/nar/gkg602
Pupo GM, Lan R, Reeves PR: Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics. Proc Natl Acad Sci U S A 2000, 97: 10567–10572. 10.1073/pnas.180094797
Janssen PJ, Audit B, Ouzounis CA: Strain-specific genes of Helicobacter pylori: distribution, function and dynamics. Nucleic Acids Res 2001, 29: 4395–4404. 10.1093/nar/29.21.4395
Charlebois RL, Clarke GD, Beiko RG, St Jean A: Characterization of species-specific genes using a flexible, web-based querying system. FEMS Microbiol Lett 2003, 225: 213–220. 10.1016/S0378-1097(03)00512-3
Krylov VN: Role of horizontal gene transfer by bacteriophages in the origin of pathogenic bacteria. Genetika 2003, 39: 595–620.
Hayashi T, Makino K, Ohnishi M, Kurokawa K, Ishii K, Yokoyama K, Han CG, Ohtsubo E, Nakayama K, Murata T, Tanaka M, Tobe T, Iida T, Takami H, Honda T, Sasakawa C, Ogasawara N, Yasunaga T, Kuhara S, Shiba T, Hattori M, Shinagawa H: Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12. DNA Res 2001, 8: 11–22. 10.1093/dnares/8.1.11
Wu CH, Huang H, Nikolskaya A, Hu Z, Barker WC: The iProClass integrated database for protein functional analysis. Comput Biol Chem 2004, 28: 87–96. 10.1016/j.compbiolchem.2003.10.003
Kall L, Krogh A, Sonnhammer EL: A combined transmembrane topology and signal peptide prediction method. J Mol Biol 2004, 338: 1027–1036. 10.1016/j.jmb.2004.03.016
Development of CUPID was supported by grant AWD4220707 from the District of Columbia Water and Sewer Authority. The bioinformatics infrastructure on which CUPID is based was supported in part by grant U01-HG02712 from the National Institutes of Health. We would like to thank Drs. Winona Barker and Hongzhan Huang for reviewing the manuscript and providing useful comments.
RM conceived, designed and coordinated the study, developed a general outline for the algorithm and drafted the manuscript. DN developed the specific algorithm and was responsible for software design and implementation, and participated in the writing of the manuscript. CW and SM participated in the design and evaluation of the study and manuscript writing. RT and DN developed the web interface. All authors read and approved the final manuscript.
Raja Mazumder, Darren A Natale contributed equally to this work.