- Open Access
Cross-species cluster co-conservation: a new method for generating protein interaction networks
Genome Biology volume 8, Article number: R185 (2007)
Co-conservation (phylogenetic profiles) is a well-established method for predicting functional relationships between proteins. Several publicly available databases use this method and additional clustering strategies to develop networks of protein interactions (cluster co-conservation (CCC)). CCC has previously been limited to interactions within a single target species. We have extended CCC to develop protein interaction networks based on co-conservation between protein pairs across multiple species, cross-species cluster co-conservation.
The exponential increase in sequence information has widened the gap between the number of predicted and experimentally characterized proteins. At present, about 400 microbial genomes are fully sequenced. The prediction of protein function from sequence is a critical issue in genome annotation efforts. Currently, the best established method for function prediction is based on sequence similarity to proteins of known function. Unfortunately, homoogy-based prediction is of limited use due to the large number of homologous protein families with no known function for any member. An alternative method for predicting protein function is the phylogenetic profiles approach, also known as the co-conservation (CC) method first introduced by Pellegrini et al. . Co-conservation predicts interactions between pairs of proteins by determining whether both proteins are consistently present or absent across diverse genomes [2–8]. CC methods have been shown to be more powerful than sequence similarity alone at predicting protein function.
Even though all CC methods rely on the premise that functionally related proteins are gained or lost together over the course of evolution, several different strategies for performing CC studies have been reported. For example, Date et al.  used real BLASTP best hit E-values normalized across 11 bins instead of binary classification for conservation, while Zheng and coworkers  constructed phylogenetic profiles using presence/absence of neighboring gene pairs. Alternatively, Pagel et al.  constructed phylogenetic profiles between domains, instead of genes, and then created domain interaction maps. Barker et al.  applied maximum likelihood statistical modeling for predicting functional gene linkages based on phylogenetic profiling. Their method detected independent instances of protein pair correlated gain or loss on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny .
Currently, several web-based databases that compile predictions of protein-protein interactions are available, for example, PLEX , String , Prolinks , and Predictome . These databases use various methods, including CC, to organize groups of proteins within individual species into clusters (cluster co-conservation (CCC)) that represent predicted protein interaction networks. Here, we have investigated the degree to which these within-species clusters are conserved across different species, using an automated method for comparing phylogenetic profiling based CCC across multiple species (CS-CCC; Figure 1). CS-CCC is essentially a meta-analysis of CCC that automates the identification of interactions that are uniquely present or absent across different species, which cannot be easily accomplished using existing methods. We have shown that this method increased groupings among proteins that function in distinct but coordinate processes and decreased groupings among proteins with unknown functions. This suggests that CS-CCC, in comparison to CCC, allows one to extend the network to better understand pathways involving proteins with multiple functions. Our intention for CS-CCC was that the identity of proteins present or absent in co-conserved clusters when evaluated across multiple species would facilitate the assignment of protein function, enable the development of novel and testable biological hypotheses, and provide experimentalists with the scientific justification required to test these hypotheses. We show these features through a number of different examples involving complex biological phenomena (that is, flagellum, chemotaxis, and biofilm proteins).
Cross-species clustered co-conservation
CS-CCC is based on the use of CC methods simultaneously across several species. As such, the reliability of the CS-CCC method is directly linked to the reliability of existing CC methods, which has been extensively documented [2–8]. Specifically, since CC methods produce protein-protein interactions involving proteins with previously uncharacterized functions, CC methods perform better than sequence similarity methods alone at predicting protein function. Here, we performed the same comparison to assess the performance of CS-CCC (up to six species) when compared to CCC alone (one species) (Figure 2a). The reliability of predicted protein interaction pairs was evaluated by using a combination of Clusters of Orthologous Groups (COG) functional categories, and The Institute for Genomic Research (TIGR) role categories (Additional data file 1). As the number of species included in our CS-CCC analysis increased, the number of predicted interactions involving proteins with unclassified functions decreased (yellow bars). Interestingly, at the lowest confidence level, the number of predicted interactions involving proteins from different functional categories increased with the number of included species. At the highest confidence level, grouping between proteins from the same functional category increased. For example, 56% of Escherichia coli K12 protein pairs (confidence level of 0.6) consisted of proteins within the same COG functional group, 19% of protein pairs were in different functional categories, and 25% had at least one unclassified member due to limited experimental data. As the number of species is expanded, these percentages range from 54-62%, 30-45%, and 0-10%, respectively. At the highest confidence level (0.8), the inclusion of 6 species resulted in almost 80% of the predicted interactions involving proteins from the same functional category. These results suggest that expanding the number of species included in the analysis, as provided for by CS-CCC, not only predicts interactions that are not predicted at different confidence levels used in CCC analysis, but also that the nature of such predicted interactions is fundamentally different. One explanation for such observations is that CS-CCC has improved capabilities for extending the protein interaction network to include the various functions required in complex biological processes (that is, regulatory relationships, nutrient transport/catabolism links, and so on). As an example of this possibility, in the CS-CCC analysis using all 6 bacterial species at confidence level 0.8 (the green bar on the far right on Figure 2a), there were 6 co-conserved protein pairs involving 9 total proteins that were not in the same COG functional category. When the larger network that these pairs fall into was extracted (Figure 2b), it became apparent that each of the proteins in question function within the context of two larger, coherent networks involving related processes. For example, rpoA and rpsD encode proteins of differing functions, yet their interaction is well conserved across multiple species within a 12-gene network of related functions. The remaining seven proteins of varying functions were also well conserved across multiple species in a larger network. These data suggest that the addition of multiple species to the analysis adds confidence to predicted interactions among proteins from different functional categories (that is, a meta-analysis). This point is exemplified via the color-coded, species specific arcs in Figure 2b, where it is clear that addition of multiple species both adds new interactions (that is, unique sub-networks) and reinforces the interactions predicted for comparison species.
CS-CCC identifies interactions that could not be identified by CCC
Our analysis of CCC across six bacterial species indicated that CS-CCC revealed unique and useful information not provided by CCC alone. As one example, CS-CCC uniquely revealed that amino-acid biosynthesis and flagellar networks are connected via FliY (Figure 3c), a component of the flagella motor-switch complex that is predicted to transport amino acids . Both E. coli and Pseudomonas aeruginosa ArgT networks revealed connections with the FliY protein (Figure 3a,b), but such networks did not include the extensive set of additional flagellar protein interactions predicted in the Bacillus subtilis network. Such information can be used to not only develop more precise hypotheses about protein function but also to provide the justification required to test such hypotheses. A second example of information uniquely revealed by CS-CCC suggests how the process of chemotaxis has evolved across species. A CS-CCC comparison of chemotaxis in E. coli K12 and Salmonella revealed that Salmonella lacks Tap, which transports maltose, but has Tcp, which transports citrate. In contrast, E. coli has Tap but lacks Tcp. CCC analysis alone does not capture this difference in chemotaxis responsiveness. As a final example, extending this CS-CCC analysis of chemotaxis proteins to include P. aeruginosa indicated new links among type IV pili and biofilm formation proteins [13, 14], suggesting that the process of chemotaxis has evolved different functional relationships in different species. These three examples provide a simple demonstration of the ability of CS-CCC to predict unique and biologically informative interactions when compared to CCC alone. The next several sections elaborate upon the specific types of interactions that CS-CCC is uniquely suited at identifying.
CS-CCC reveals how proteins that function in distinct but coordinated processes may have evolved
Chemotaxis proteins are co-conserved across the examined bacteria (Figure 4). Three classes of proteins are essential for chemotaxis: transmembrane receptors, cytoplasmic signaling components, and enzymes for adaptive methylation. The transmembrane receptors are two-component signal transduction complexes called methyl-accepting chemotaxis proteins (MCPs). E. coli MCPs are Tsr, Tar, Trg, Tap, and Aer, and each recognizes specific sugars, amino acids or dipeptides (Figure 4a,c). Even though different bacteria have different MCPs, they are highly co-conserved among Gram-negative and positive bacteria. For example, Salmonella lacks Tap, which recognizes maltose, but has Tcp, a citrate sensor , which is co-conserved with the other Salmonella MCPs (Figure 4b,c). The cytoplasmic signaling components transmit signal between the MCP receptors and the flagellar apparatus. These proteins are CheA, CheW, CheY and CheZ, and they are not co-conserved among the bacteria. CheZ is not co-conserved because it has no homology across many bacteria . CheY is likely not co-conserved because it functions with CheZ. CheA and CheW are sometimes co-conserved and sometimes not, which may suggest that they function independently in different bacteria. The enzymes for adaptive methylation, CheB and CheR, modulate signaling of the cytoplasmic proteins, and both of these proteins are highly co-conserved among all six bacteria. Thus, chemotaxis analysis illustrates two important points. First, the CS-CCC method reveals species differences in protein interaction, including co-conserved pairs that are unique to a given species or that are common across select species (Figure 4c). For instance, the sequences of CheA and CheW are conserved but the proteins are not co-conserved, suggesting that their interactions and functions may differ among bacterial species. Second, the CS-CCC method yields information that functional assays do not. For instance, different MCPs recognize different ligands and yet are co-conserved because they function in the same pathway.
Figure 4 shows a cluster containing proteins that function in distinct but inter-dependent processes. For instance, in P. aerginosa, flagella, chemotaxis machinery, and type IV pili are important for bacterial biofilm formation [13, 14] and are co-conserved. Type IV pili mediate twitching motility, which is important for subsequent spreading of the bacteria over the surface and the formation of microcolonies within a developing biofilm . Twitching motility proteins PilJ and PilK are co-conserved within this cluster and are highly interconnected with flagella and chemotaxis proteins. Flagellar motility appears to be required for approaching surfaces, and 17 flagellar proteins are co-conserved (Figure 4c). Chemotaxis is required for the bacteria to swim towards nutrients associated with a surface. P. aerginosa has two chemotaxis signaling systems, and proteins representing both are in the biofilm cluster (CheR1, CheR2, CheA, CheW, PA0173, PA0178; PctA, PctB, PctC). These data suggest that chemotaxis, flagella, and pili proteins may be co-conserved because they all contribute to biofilm formation. Moreover, the inclusion of P. aerginosa in the CS-CCC analysis brought pili proteins into the biofilm cluster, suggesting that in some bacteria, all of these processes co-evolved. Thus, CS-CCC can identify co-conserved networks of proteins that function in biochemically distinct pathways but that contribute to complex biological phenomenon.
RpoN connects RpoN-regulated proteins with flagella and with type III secretion system proteins
In some of the bacteria studied, RpoN (also known as σ54 or SigL) clustered with RpoN-regulated proteins and flagella proteins are clustered with type III secretion system proteins (Figure 4c). Flagellar proteins are cluster co-conserved with specific components of type III secretion systems (T3SS), which are important for virulence in Salmonella enterica serotype Typhimurium LT2, E. coli O157, Shigella flexneri and P. aerginosa  (Table 1). The T3SS of Shigella is not chromosomally encoded and so was not included in our analysis. The three subunits of the T3SS and flagella that are co-conserved are integral inner membrane proteins of the flagellar or T3SS export apparatus that forms the channel through which proteins are secreted . S. typhimurium LT2 and E. coli O157 both encode two T3SSes, and the corresponding proteins from each are within this cluster. In E. coli K12, S. typhimurium LT2, and B. subtilis, RpoN connects the RpoN-regulated and the flagellar/T3SS clusters. This is consistent with experimental data that flagellar genes (flhA and flhB) are activated by RpoN . Thus, RpoN likely connects two distinct clusters because it regulates proteins in both clusters. This demonstrates that because CS-CCC examines multiple genomes simultaneously, it has the power to show that proteins unique to particular organisms may function with proteins common to multiple organisms, enabling the placement of unstudied proteins within a broader biological context.
CS-CCC can be used to assign function to unstudied proteins
Genes that function in biofilm formation
Figure 5a shows two large clusters of proteins built around YegE or YfiN in E. coli K12 and P. aeruginosa. These clusters are co-conserved with variable numbers of proteins among all of our Gram-negative bacteria. Even though most of these proteins have unknown function, many have GGDEF (Gly-Gly-Asp-Glu-Phe) or EAL (Glu-Ala-Leu) domains, which have been implicated in expression of biofilm phenotypes . Interestingly, each protein of known function within this cluster in PAO1 (WspR, MorA, and FimX) has also been implicated in biofilm phenotypes. WspR is a response regulator that activates pili adhesion genes required for biofilm formation . MorA is a membrane-localized negative regulator of the timing of flagellar formation and plays a role in the establishment of biofilms . FimX is required for a type of twitching motility critical to biofilm formation . FimX is a signal sensing protein with phosphotransfer activity and a GGDEF domain. GGDEF encodes a dinucleotide cyclase that generates cyclic di-GMP and is present in all proteins known to be involved in the regulation of cellulose synthesis. Cyclic di-GMP is a novel bacterial second messenger that directs the transition from sessility to motility . Cyclic di-GMP is degraded by proteins with EAL domains, which are cyclic dinuclotide phosphodiesterases . Proteins containing the GGDEF and EAL domain can regulate biofilm formation and/or cell aggregation by controlling the levels of cyclic di-GMP . Interestingly, most of the proteins in these large clusters have GGDEF or EAL domains. Of the 44 known P. aeruginosa proteins with GGDEF or EAL domains , 34 are in this cluster; 19 have GGDEF and 15 have EAL domains. E. coli K12 has a similar cluster of GGDEF and EAL domains (Figure 5a). The 25 proteins within this cluster are highly interconnected. Of the 38 E. coli K12 known GGDEF or EAL domain containing proteins , 24 are co-conserved within this cluster. EvgS is a sensor protein for a two component regulatory system  that is also within this cluster. Evgs is involved in quorum sensing and may be important in biofilm establishment or maintenance. Over-expression of evgS causes abnormal biofilm architecture  and previous studies also noted that quorum sensing is involved in biofilm formation . Our experimental data show that four of the GGDEF domain containing proteins in the network of Figure 5a that previously had no known function do indeed mediate biofilm formation . Similar biofilm clusters were identified by the CS-CCC method in all of the Gram-negative bacteria we examined. Thus, by clustering together unstudied proteins, whether or not they have sequence homology, CS-CCC suggests that these proteins may function in a common phenomenon.
Small clusters can contain proteins that function in the same processes
Examination of small protein clusters revealed that most pairs or triplets contain proteins that function in the same processes. To further test this observation, we experimentally examined the triplet containing YcgB, YeaH, and YeaG, which cluster together across different bacteria (Figure 5b). Because independent data indicate that yeaH, but not yeaG, contributes to antimicrobial peptide resistance in S. typhimurium , we determined whether strains lacking ycgB have a similar phenotype. Strains lacking ycgB were indeed sensitive to antimicrobial peptides (unpublished data). Thus, CS-CCC analyses revealed previously unknown protein interactions that provided sufficient justification to test a specific biological hypothesis suggested by these interactions.
When proteins are not identified as co-conserved using CS-CCC
In this study, we have shown that CS-CCC of proteins provides important information. Both the presence and the absence of clustered co-conservation for any given protein are informative. There are at least two reasons why proteins that function together are not co-conserved in a species: first, a protein is found only in certain organisms or a protein function is performed by different proteins in different organisms; and second, a result is a false negative.
A protein is found only in certain organisms: T3SS effectors
Effector proteins are secreted by T3SS machinery and function to alter host cell physiology . A bacterial species can have many effectors but they generally do share apparent sequence homology, either within or between bacteria . We examined 49 known SPI2 and SPI1 effectors in S. typhimurium LT2 and 40 known effectors in P. aeruginosa and found that none of these proteins are co-conserved. In contrast, some of the known translocon T3SS proteins, which form the secretion apparatus, are highly co-conserved (Figure 4c). Thus, while CS-CCC offers insights into the function of proteins that are co-conserved, our results show that some of the non co-conserved proteins, such as effectors, are organism specific.
A result is a false negative: flagella and RpoN
Our analysis of false negatives reveals that the CS-CCC method produces some false negatives. For instance, there is no co-conservation between RpoN and flagella in E. coli 0157, S. flexneri and P. aeruginosa (Figure 4c). However, it has been experimentally shown in P. aeruginosa that many flagellar genes, such as flhA and flhB, are regulated by RpoN . In addition, an RpoN consensus sequence is located in the intergenic region between flhB and flhA . These data suggest that the absence of co-clustering of RpoN with flagellar proteins in P. aeruginosa is a false negative result. Thus, when proteins are not co-conserved, it cannot be concluded that they are functionally unrelated. This result further underlines the value of developing and comparing interaction networks from multiple genomes when attempting to infer function.
There are also some situations in which a result is both a false negative and the protein in question is found only in certain organisms. The bacterial flagellum is a complex molecular system with multiple components required for functional motility. It extends from the cytoplasm to the cell exterior. Not only are flagella organelles of locomotion, but they also play important roles in attachment and biofilm formation. There are common themes in flagellar protein control and assembly, but there also appears to be variation among organisms. Some of the flagellar proteins are not co-conserved in any of the bacteria of our study, such as, three ring proteins (FlgH, FlgI, and FliF), and some of the axle-like proteins FliE, FlgB, FlgF, FlgL, and FliD. FliE has been shown to physically interact with FlgB . The stator motor proteins MotA and MotB are also not co-conserved. Thus, CS-CCC analysis of the flagellar cluster yields both false negative results and is also a consequence of species-specific proteins. This also illustrates that determining why proteins are not co-conserved can be difficult, without additional information.
Large volumes of data make computational methods feasible, exciting, and preferable to gene-by-gene homology searches. We have shown that use of CS-CCC expands protein interaction networks to include proteins with distinct functions that are involved in coherent biological processes, offers insight into the function of uncharacterized proteins, reveals unique information about each genome examined, and gives insight into the process of evolution.
Protein co-conservation can be a result of many factors, including vertical inheritance or functional selection. Thus, we have examined patterns of CCC within and across several bacteria using CS-CCC. Our analysis showed that this computational approach provides us with more information than the traditional homology approaches or CCC. Homology approaches to protein function are based on similarity to other proteins with known functions and are limited by the fact that many proteins have unknown functions. While homology-based methods can be effective for predicting the functions of remote homologs, these methods perform poorly as the evolutionary distance between homologous proteins increases. Even a sophisticated homology-based method fails to successfully assign functions to most of the proteins for a particular organism. CCC, on the other hand, is not strictly based on homology but is limited by its ability to analyze only a single species at a time. In contrast, CS-CCC examines each cluster across multiple species and reveals interactions that both homology-based methods and CCC fail to identify. Use of CS-CCC allows researchers to extend the protein interaction network to better understand pathways involving multiple proteins with multiple functions. Therefore, the CS-CCC method is a significant advance and will be useful for researches in many different fields of biology.
Prediction by CS-CCC provided us with global views of six complete bacterial genomes. Identification by CS-CCC of proteins that cluster together enabled more accurate predictions of the biological roles that proteins with previously unstudied functions may play. For instance, proteins that function in distinct but coordinated processes can be co-conserved across species even though not all processes occur in all bacteria (Figure 4c). In addition, in large, highly interconnected clusters in which most of the proteins have unknown functions, it is likely that they all function together in a common phenomenon. The GGDEF/EAL cluster is an example of this, as many of the previously unknown proteins in this cluster play roles in biofilm formation (Figure 5a). Even small protein clusters identified by CS-CCC are likely to consist of proteins that function in the same process, as shown by COG/TIGR analysis and experimentally (Figure 5b). These analyses provide evidence that the CS-CCC method is a reliable predictor of functional relationships.
For any given method, there are advantages and disadvantages. The number of false positives and false negatives is a key measurement of accuracy. In our case, the number of false negatives is not possible to estimate without performing many additional laboratory experiments. However, our evaluation of CS-CCC showed that the number of false positives was low. Since this method was evaluated based on our selected bacteria, there may be some bias toward overestimation of accuracy when applied to other organisms, and this remains to be tested. In addition, we have shown that our results can be sensitive to the number of bacteria included in our analysis. Finally, there may be some aspects of the bacteria we chose that are not representative of other bacteria, further reducing the generality of these results. Thus, while the report here represents a compelling demonstration of the value of performing CCC across multiple species, future efforts should be focused on developing better understanding of which and how many organisms to include in CS-CCC studies.
Materials and methods
Bacteria used to create CS-CCC graphs
We chose to focus on the Gamma subgroup of proteobacteria because members of this subgroup are among the best characterized, including whole genome sequences and curated datasets of protein functions and interactions. The genomes of five closely related Gamma Gram-negative and one low G+C bacteria (B. subtilis) were used to evaluate the CCC method. Substantial experimental data exist for all six bacteria. The gammaproteobacteria included E. coli (K12 and O157-O157:H7 EDL933), S. flexneri (2a str. 2457T), S. typhimurium (LT2), and P. aeruginosa (PAO1). E. coli (K12) is the most intensively studied Gram-negative bacteria and is the closest studied relative of P. aeruginosa, and S. typhimurimum LT2. E. coli (O157-O157:H7 EDL933) is a clinical isolate from raw hamburger meat implicated in hemorrhagic colitis outbreak, and S. typhimurium LT2 causes enteritis in humans. P. aeruginosa is an opportunistic pathogen and is the major cause of morbidity and mortality in patients with cystic fibrosis; P. aeruginosa PAO1 was isolated from a wound . P. aeruginosa is a versatile Gram-negative bacterium that also thrives in soil, marshes and coastal marine habitats, and on plant tissues . E. coli K12 diverged 4.5 million years ago (MYA) from O157, an estimated 100 MYA from Salmonella, 200 MYA from Pseudomonas, and 1,200 MYA from Bacillus. Thus, we examined a combination of pathogenic and non-pathogenic organisms that range from closely to distantly related.
Construction of CS-CCC graphs
We began construction of CS-CCC graphs (Figure 1) using predictions of pairwise protein-protein interactions based on phylogenetic profiles (CC methods; Figure 1a). Currently, several databases that compile predictions are available, including Prolinks , String , and Predictome . We used the Prolinks Database 2.0, which contains a total of 168 microbial genomes, including 10 eukaryotes, 16 Archaea, and 142 Bacteria . Even though ProLinks provides predicted interactions based on a number of different methods (that is, Rosetta stone, gene neighbors, and so on), we have used only interactions prediction by the phylogenetic profiling method in this study. We chose not to use the STRING database as a source of predictions because it conflates co-conservation with orthology information from the COG database ; we used COG functional category and TIGR functional role category data to evaluate purely co-conservation inferences. Predictome  was not used because it does not provide statistical measures to evaluate the accuracy of each prediction. For each pair assignment (CC), we required a confidence scheme using phylogenetic profiling of at least 60% according to the Prolinks scoring scheme . An E-value of less than 10-10 was used as the threshold for BLASTP in Prolinks to define a homolog of a query protein to be present in a secondary genome. For each bacterial genome analyzed, the number of assigned pairs is shown in Table 2. For each bacterial species, we mapped accession IDs from Prolinks predicted protein pairs to NCBI  and then to EcoCyc  for E. coli K12, P. aeruginosa  for P. aeruginosa and B. subtilis  for B. subtilis. We matched corresponding proteins between species by protein name or synonym. We then constructed CCC graphs using the pairwise links for each species (Figures 1b and 6) using a binary adjacency matrix where 1 indicates the corresponding pair was co-conserved, and 0 otherwise. Networks were represented by graphs in which each node represents a protein and each edge represents an interaction that links two proteins. Network graphs were visualized using Cytoscape , an open-source, platform-independent environment software. The lengths of the lines connecting proteins hold no meaning and vary to facilitate viewing of the network. Each network is color-coded based on protein function categories, as described in the corresponding figure legends. The assignment of putative functions was based on EcoCyc, Pseudomonas.com, NCBI and SubtiList, as given in the links above. For separation of connected components of the network and building clusters of proteins, we used breadth-first search (BFS) graph algorithms.
Finally, for comparison of each cluster across different species (CS-CCC), we used BFS to build a network (source network) for a set of target proteins from the source genome. We then built networks for each additional organism that contained proteins with the same name as at least one of the proteins from the source networks. This process identifies proteins and protein interactions that are consistently identified across multiple species (colored gray in Figure 1c) or that are unique to individual species (colored red in Figure 1c). This same method can be used to further parse such networks to identify combined, common and unique networks present for specific proteins across a collection of organisms (Figure 1c). In this way, CS-CCC builds on information generated by CCC (Figure 1b) to provide more accurate and genome-specific protein function assignment. We used protein name to map links across conserved species (thus, links are not explicitly based on orthology) [37–39]. Like all methods, the use of protein names has both advantages and disadvantages. Here, protein name was chosen in order to validate that CS-CCC provides new and biologically informative data not accessible by CCC alone. For this purpose, we chose to validate this method using named proteins where functional information was available. While this is appropriate for method validation, the disadvantage is that there are problems with annotation due in part to a lack of standardization, which would limit the number of proteins for which this analysis can be reliably performed. In light of this limitation, we considered using reciprocal homology as an alternative to protein name. We found that this introduces unacceptable levels of cross-talk, much of which is likely noise. Addressing this limitation is an important area for continued effort.
Data are available upon request.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a figure that shows the reliability of predicted protein interaction pairs using TIGR role categories at three different confidence levels.
Clusters of Orthologous Groups
cross-species clustered co-conservation
methyl-accepting chemotaxis protein
million years ago
The Institute for Genomic Research
type III secretion systems.
Pellegrini M, Marcotte EM, Thompson MJ, Eisenberg D, Yeates TO: Assigning protein functions by comparative genome analysis: protein phylogenetic profiles. Proc Natl Acad Sci USA. 1999, 96: 4285-4288. 10.1073/pnas.96.8.4285.
Marcotte EM, Pellegrini M, Ng HL, Rice DW, Yeates TO, Eisenberg D: Detecting protein function and protein-protein interactions from genome sequences. Science. 1999, 285: 751-753. 10.1126/science.285.5428.751.
Strong M, Mallick P, Pellegrini M, Thompson MJ, Eisenberg D: Inference of protein function and protein linkages in Mycobacterium tuberculosis based on prokaryotic genome organization: a combined computational approach. Genome Biol. 2003, 4: R59-10.1186/gb-2003-4-9-r59.
Huynen M, Snel B, Lathe W, Bork P: Predicting protein function by genomic context: quantitative evaluation and qualitative inferences. Genome Res. 2000, 10: 1204-1210. 10.1101/gr.10.8.1204.
Mellor JC, Yanai I, Clodfelter KH, Mintseris J, DeLisi C: Predictome: a database of putative functional links between proteins. Nucleic Acids Res. 2002, 30: 306-309. 10.1093/nar/30.1.306.
Bowers PM, Pellegrini M, Thompson MJ, Fierro J, Yeates TO, Eisenberg D: Prolinks: a database of protein functional linkages derived from coevolution. Genome Biol. 2004, 5: R35-10.1186/gb-2004-5-5-r35.
Date SV, Marcotte EM: Protein function prediction using the Protein Link EXplorer (PLEX). Bioinformatics. 2005, 21: 2558-2559. 10.1093/bioinformatics/bti313.
von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B: STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003, 31: 258-261. 10.1093/nar/gkg034.
Zheng Y, Roberts RJ, Kasif S: Genomic functional annotation using co-evolution profiles of gene clusters. Genome Biol. 2002, 3: RESEARCH0060-10.1186/gb-2002-3-11-research0060.
Pagel P, Wong P, Frishman D: A domain interaction map based on phylogenetic profiling. J Mol Biol. 2004, 344: 1331-1346. 10.1016/j.jmb.2004.10.019.
Barker D, Pagel M: Predicting functional gene links from phylogenetic-statistical analyses of whole genomes. PLoS Comput Biol. 2005, 1: e3-10.1371/journal.pcbi.0010003.
Mytelka DS, Chamberlin MJ: Escherichia coli fliAZY operon. J Bacteriol. 1996, 178: 24-34.
Korber DR, Lawrence JR, Caldwell DE: Effect of motility on surface colonization and reproductive success of Pseudomonas fluorescens in dual-dilution continuous culture and batch culture systems. Appl Environ Microbiol. 1994, 60: 1421-1429.
Pratt LA, Kolter R: Genetic analysis of Escherichia coli biofilm formation: roles of flagella, motility, chemotaxis and type I pili. Mol Microbiol. 1998, 30: 285-293. 10.1046/j.1365-2958.1998.01061.x.
Manson MD, Armitage JP, Hoch JA, Macnab RM: Bacterial locomotion and signal transduction. J Bacteriol. 1998, 180: 1009-1022.
Galan JE, Collmer A: Type III secretion machines: bacterial devices for protein delivery into host cells. Science. 1999, 284: 1322-1328. 10.1126/science.284.5418.1322.
Mecsas JJ, Strauss EJ: Molecular mechanisms of bacterial virulence: type III secretion and pathogenicity islands. Emerg Infect Dis. 1996, 2: 270-288.
Fleiszig SM, Arora SK, Van R, Ramphal R: FlhA, a component of the flagellum assembly apparatus of Pseudomonas aeruginosa, plays a role in internalization by corneal epithelial cells. Infect Immun. 2001, 69: 4931-4937. 10.1128/IAI.69.8.4931-4937.2001.
Simm R, Morr M, Kader A, Nimtz M, Romling U: GGDEF and EAL domains inversely regulate cyclic di-GMP levels and transition from sessility to motility. Mol Microbiol. 2004, 53: 1123-1134. 10.1111/j.1365-2958.2004.04206.x.
Spiers AJ, Bohannon J, Gehrig SM, Rainey PB: Biofilm formation at the air-liquid interface by the Pseudomonas fluorescens SBW25 wrinkly spreader requires an acetylated form of cellulose. Mol Microbiol. 2003, 50: 15-27. 10.1046/j.1365-2958.2003.03670.x.
Choy WK, Zhou L, Syn CK, Zhang LH, Swarup S: MorA defines a new class of regulators affecting flagellar development and biofilm formation in diverse Pseudomonas species. J Bacteriol. 2004, 186: 7221-7228. 10.1128/JB.186.21.7221-7228.2004.
Huang B, Whitchurch CB, Mattick JS: FimX, a multidomain protein connecting environmental signals to twitching motility in Pseudomonas aeruginosa. J Bacteriol. 2003, 185: 7068-7076. 10.1128/JB.185.24.7068-7076.2003.
Karp PD, Riley M, Saier M, Paulsen IT, Collado-Vides J, Paley SM, Pellegrini-Toole A, Bonavides C, Gama-Castro S: The EcoCyc Database. Nucleic Acids Res. 2002, 30: 56-58. 10.1093/nar/30.1.56.
Georgellis D, Kwon O, De Wulf P, Lin EC: Signal decay through a reverse phosphorelay in the Arc two-component signal transduction system. J Biol Chem. 1998, 273: 32864-32869. 10.1074/jbc.273.49.32864.
Tenorio E, Saeki T, Fujita K, Kitakawa M, Baba T, Mori H, Isono K: Systematic characterization of Escherichia coli genes/ORFs affecting biofilm formation. FEMS Microbiol Lett. 2003, 225: 107-114. 10.1016/S0378-1097(03)00507-X.
Ren D, Bedzyk LA, Ye RW, Thomas SM, Wood TK: Differential gene expression shows natural brominated furanones interfere with the autoinducer-2 bacterial signaling system of Escherichia coli. Biotechnol Bioeng. 2004, 88: 630-642. 10.1002/bit.20259.
Lynch MD, Warnecke T, Gill RT: SCALEs: multiscale analysis of library enrichment. Nat Methods. 2007, 4: 87-93. 10.1038/nmeth946.
Erickson KD, Detweiler CS: The Rcs phosphorelay system is specific to enteric pathogens/commensals and activates ydeI, a gene important for persistent Salmonella infection of mice. Mol Microbiol. 2006, 62: 883-894. 10.1111/j.1365-2958.2006.05420.x.
Gophna U, Ron EZ, Graur D: Bacterial type III secretion systems are ancient and evolved by multiple horizontal-transfer events. Gene. 2003, 312: 151-163. 10.1016/S0378-1119(03)00612-7.
Hueck CJ: Type III protein secretion systems in bacterial pathogens of animals and plants. Microbiol Mol Biol Rev. 1998, 62: 379-433.
Saijo-Hamano Y, Uchida N, Namba K, Oosawa K: In vitro characterization of FlgB, FlgC, FlgF, FlgG, and FliE, flagellar basal body proteins of Salmonella. J Mol Biol. 2004, 339: 423-435. 10.1016/j.jmb.2004.03.070.
Stover CK, Pham XQ, Erwin AL, Mizoguchi SD, Warrener P, Hickey MJ, Brinkman FS, Hufnagle WO, Kowalik DJ, Lagrou M, et al: Complete genome sequence of Pseudomonas aeruginosa PA01, an opportunistic pathogen. Nature. 2000, 406: 959-964. 10.1038/35023079.
NCBI Genbank Protein Annotation. [http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi]
Winsor GL, Lo R, Sui SJ, Ung KS, Huang S, Cheng D, Ching WK, Hancock RE, Brinkman FS: Pseudomonas aeruginosa Genome Database and PseudoCAP: facilitating community-based, continually updated, genome annotation. Nucleic Acids Res. 2005, D338-343. 33 Database
Institut Pasteur. [http://genolist.pasteur.fr/SubtiList/]
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.
Snel B, van Noort V, Huynen MA: Gene co-regulation is highly conserved in the evolution of eukaryotes and prokaryotes. Nucleic Acids Res. 2004, 32: 4725-4731. 10.1093/nar/gkh815.
Penkett CJ, Morris JA, Wood V, Bahler J: YOGY: a web-based, integrated database to retrieve protein orthologs and associated Gene Ontology terms. Nucleic Acids Res. 2006, W330-334. 10.1093/nar/gkl311. 34 Web Server
Fitch WM: Homology a personal view on some of the problems. Trends Genet. 2000, 16: 227-231. 10.1016/S0168-9525(00)02005-9.
We thank Norman Pace for excellent discussions, Daniel Barker and Sonia M Leach for reading the manuscript and helpful comments. We also thank Kevin B Cohen for helpful comments. This study was supported by NSF grant BES0228584, and NIH grants K25_AI064338, R01-AI-072492A, and R01-LM-008111.
AK implemented the methods and analyzed the data. CSD interpreted the results. The manuscript was written by AK, CSD and edited by RTG and LH. KDE performed experiments. RTG oversaw all biological aspects of the work and LH supervised the computational aspect.
Anis Karimpour-Fard, Corrella S Detweiler contributed equally to this work.
Electronic supplementary material
Additional data file 1: Comparison of TIGR functional categories of predicted pairs at three different confidence levels. The first method (1) used only E. coli K12. Each subsequent method added an additional (underlined) bacterial strain. 1, E. coli K12; 2, E. coli K12 and E. coli O157; 3, E. coli K12, E. coli O157 and S. flexneri; 4, E. coli K12, E. coli O157, S. flexneri, and S. typhimurium LT2; 5, E. coli K12, E. coli O157, S. flexneri, S. typhimurium LT2, and P. aeruginosa; 6, E. coli K12, E. coli O157, S. flexneri, S. typhimurium LT2, P. aeruginosa, and B. subtilis. Same functional category (blue); different functional category (green); at least one protein is unclassified (yellow). (PDF 17 KB)
Authors’ original submitted files for images
About this article
Cite this article
Karimpour-Fard, A., Detweiler, C.S., Erickson, K.D. et al. Cross-species cluster co-conservation: a new method for generating protein interaction networks. Genome Biol 8, R185 (2007) doi:10.1186/gb-2007-8-9-r185
- Additional Data File
- Protein Interaction Network
- Protein Pair
- Phylogenetic Profile
- Chemotaxis Protein