Functional associations of proteins in entire genomes by means of exhaustive detection of gene fusions
© Enright and Ouzounis, licensee BioMed Central Ltd 2001
Received: 19 March 2001
Accepted: 2 July 2001
Published: 28 August 2001
It has recently been shown that the detection of gene fusion events across genomes can be used for predicting functional associations of proteins, including physical interaction or complex formation. To obtain such predictions we have made an exhaustive search for gene fusion events within 24 available completely sequenced genomes.
Each genome was used as a query against the remaining 23 complete genomes to detect gene fusion events. Using an improved, fully automatic protocol, a total of 7,224 single-domain proteins that are components of gene fusions in other genomes were detected, many of which were identified for the first time. The total number of predicted pairwise functional associations is 39,730 for all genomes. Component pairs were identified by virtue of their similarity to 2,365 multidomain composite proteins. We also show for the first time that gene fusion is a complex evolutionary process with a number of contributory factors, including paralogy, genome size and phylogenetic distance. On average, 9% of genes in a given genome appear to code for single-domain, component proteins predicted to be functionally associated. These proteins are detected by an additional 4% of genes that code for fused, composite proteins.
These results provide an exhaustive set of functionally associated genes and also delineate the power of fusion analysis for the prediction of protein interactions.
Recent progress in genome analysis has shown that it is possible to predict protein interactions or, more generally, functional associations of proteins using genome sequences alone [1,2,3]. These powerful methods rely on the observation that pairs of genes encoding proteins of known function (usually interacting or forming a complex) tend to be found in other species as a fused gene encoding a single multifunctional protein . This type of event is known as gene fusion and is a well-known process in molecular evolution . Many of these gene fusion events appear to be selectively advantageous by decreasing the regulational load in the cell for a particular process [1,3,5]. Therefore, the detection of gene fusions in one genome (defined as 'composite' proteins) allows the prediction of functional associations between homologous genes that remain separate in another genome (defined as 'component' proteins).
Although gene fusion events appear to be relatively rare, the accurate detection of a gene fusion event in one genome allows interactions to be predicted between many proteins in other genomes. It is this kind of one-to-many relationship that makes this method unique for discovering possible interactions or functional associations between proteins, even for those of unknown function. Unlike previous methods that rely on gene proximity to predict functional coupling , this robust method can also detect distal genes within a genome that may be involved in the same process. Furthermore, we have previously demonstrated  the high precision of our algorithm, which with an additional constraint of minimum alignment overlap has now increased to over 86% (see Materials and methods). This family of sequence-based methods is analogous with and complementary to the experimental approaches for the detection of protein interaction .
Genomes used in the present analysis
Organism name (strain)
Number of ORFs
Aeropyrum pernix (K1)
Aquifex aeolicus (VF5)
Archaeoglobus fulgidus (DSM4304)
Bacillus subtilis (168)
Borrelia burgdorferi (B31) + plasmids
Chlamydia pneumoniae (CWL029)
Chlamydia trachomatis (serovar D)
Escherichia coli (K12-MG1655)
Haemophilus influenzae (KW20)
Helicobacter pylori (26695)
Helicobacter pylori (J99)
Methanococcus jannaschii (DSM 2661)
Methanobacterium thermoautotrophicum (delta)
Mycoplasma genitalium (G-37)
Mycoplasma pneumoniae (M129)
Mycobacterium tuberculosis (H37Rv)
Pyrococcus horikoshii (shinkaj OT3)
Rickettsia prowazekii (Madrid E)
Saccharomyces cerevisiae (S288C)
Synechocystis sp. (PCC 6803)
Thermotoga maritima (MSB8)
Treponema pallidum (Nichols)
Our algorithm was applied individually for each of the 24 genomes, against the remaining 23 genomes which are used as references (see also Materials and methods). Paralogy in the query genome makes it difficult to determine precisely the actual number of possible associations. As we have previously pointed out, paralogy in the query genome increases uncertainty, while paralogy in the reference genome increases the fidelity of the predictions . It is for this reason that detected component and composite proteins from all genomes are subsequently clustered according to sequence similarity . Each cluster should therefore indicate a distinct family of component or composite proteins. The analysis of the distribution of these gene fusion classes among genomes allows us to investigate the dynamics and distribution of this evolutionary process and to assess the extent of the predictive power of the approach.
The detection of gene fusion events yielded 132,812 component and 66,406 composite proteins in an all-against-all genome comparison, but these values represent multiple occurrences of the same proteins across species. Of these, there are 7,224 component and 2,365 composite unique proteins across the 24 species (a 18- and 28-fold reduction respectively). The multiple detection of these cases within or across genomes signifies that the majority of components and composites are observed more than once and therefore represent genuine cases (as opposed to sequencing artifacts, which are usually isolated cases).
Currently, the only species for which predictions can be extensively validated is the yeast Saccharomyces cerevisiae, given the ongoing work on transcript profiling  and two-hybrid technology . For yeast, there are 440 distinct component cases (predicted by all other genomes as reference, excluding some highly paralogous Drosophila melanogaster homologs) involved in 706 predicted interactions, most of which are detected by their homology to composite proteins from Caenorhabditis elegans and D. melanogaster. Two examples of predicted protein pairs that are known to interact are CPA1 (YOR303w) with CPA2 (YJR109c)  and MET3 (YJR010w) with MET14 (YKL001c) , both derived from C. elegans homologs.
We have attempted to test the validity of our predictions by comparing the set of components to a list of potentially interacting gene products, using results from a large-scale two-hybrid experiment . However, there is only one case shared between the 1,004 proteins involved in 957 putative interactions detected by the two-hybrid system and the complete set of 706 pairs in this analysis: YIL033C (SRA1) and YKL166C (TPK3) matching the C. elegans protein C09G4.2 and D. melanogaster protein CT10911. This very low count of common pairs may be expected by the sampling biases of the two rather independent methodologies, given that each approach can only detect a very small subset of the total number of actual interacting pairs in yeast. Interestingly, based on a simple conditional probability calculation, an estimate for the total number of detectable interactions in the yeast cell may be of the order of 675,000.
We have analyzed the S. cerevisiae predictions and detected many interesting cases, which appear to be hitherto undetected functional associations between yeast proteins. Two of these are discussed in some detail here. First, MXR1 (peptide methionine sulfoxide reductase, involved in anti-oxidative processes)  and YCL033C (function unknown) are predicted to be functionally associated by virtue of gene fusion in three other species - Helicobacter pylori (both strains), Haemophilus influenzae and Treponema pallidum. This observation is supported by experimental results . MXR1 is 39% identical to the amino terminus of the H. pylori composite proteins and YCL033C is 38% identical to the carboxyl terminus of these proteins. It appears that YCL033C is a selenoprotein, also homologous to the human SelX protein, which may be involved in scavenging reactive oxygen species . These two proteins may be associated to protect the yeast cell from oxidative damage.
Second, another interesting observation involves yeast proteins MSS4 (phosphatidylinositol 4-phosphate kinase), which is involved in a signaling pathway responsible for the cell-cycle-dependent organization of actin cytoskeleton , and CCT3 (cytoplasmic chaperonin subunit gamma) which is involved in microtubule and actin assembly . A central domain of CCT3 is 25% identical to a large domain of C. elegans protein VF11C1L.1 and the carboxy-terminal domain of MSS4 is 29% identical to its carboxyl terminus. Thus, these two proteins are predicted to cooperate in cell-cycle-dependent cytoskeleton organization and assembly.
This has been a key observation that dictated the all-against-all genome comparison in this analysis. In other words, when species A is used as a query against species B, the resulting set of component and composite proteins is different from that with the reverse comparison, when species B is used as a query against species A. The three principal factors in gene fusion during evolution appear to be paralogy, genome size and phylogenetic distance. For instance, larger genomes have more composite, possibly paralogous, proteins. At the same time, closely related species evidently show similar patterns of gene fusion. The results below address each of these factors in turn and examine their relative contribution to the gene fusion process and their effects on the prediction of functional association of proteins.
For every genome, both sets of component and composite proteins were subsequently clustered , to detect the degree of paralogy for these proteins (Figure 3, green bars; see also sequence clustering in Materials and methods). There are 2,534 component families, with an average of 105 families per genome (Figure 3a, green bars) and 1,323 composite families, with an average of 55 families per genome (Figure 3b, green bars). Comparing these numbers with the number of unique cases, it is evident that there is a paralogy level of two- to threefold per genome for the composite and component proteins, respectively. As mentioned above, this effect contributes to the confidence of the predictions, depending on whether paralogy is observed in the query or the reference genome.
Another characteristic of this process is the redundancy of both sets of component and composite cases: the number of instances of these may be high but they are widely present across species, falling into well-defined protein families. When all components and composites are clustered as a single set (as opposed to within species, above), sequence clustering results in 1,287 single component families and 621 single composite families (as represented in the current analysis for the 24 species). Comparing these numbers with the number of families per species, it is apparent that there is a further twofold reduction for both sets. This result indicates that gene fusion is widespread in evolution but forms a finite set. Different species may contain a common core of composite families, but also provide new families that are used to predict functional association. For instance, D. melanogaster provides far more composite families (more than 200) compared to C. elegans (fewer than 100) (Figure 3b, green bars). Genomes with unique composite families, such as D. melanogaster, contribute strongly to the majority of predicted interactions. It may also be that only certain classes of proteins are involved in gene fusion and that there is an upper limit for the predictive power of this approach obtainable from (currently available) 621 families.
The exhaustive detection of gene fusion events in entire genome sequences allows the prediction of functionally associated components based merely on genome structure. The all-against-all species comparison is a necessary step because we have repeatedly observed fused, composite proteins in taxonomically lower organisms. The landscape of gene fusions appears to be a complex one, affected by paralogy, genome size and phylogenetic distance.
Although gene fusion is widely present across various phylogenetic groups, it is a process that may involve only certain types of proteins. Yet, this approach for the prediction of functional associations of proteins results in robust predictions for physical interactions, pathway involvement, complex formation and other types of functional associations of protein molecules.
With the present analysis, we delineate the available universe of fusion events and detect a set of 621 composite protein families from which predictions may be obtained. This approach results in 39,730 pairs of functionally associated proteins across 24 species, with high precision and coverage. This novel set of predictions is made available to the scientific community for the first time, and we believe that many of these cases can be subsequently verified by experimental methods.
Materials and methods
All 24 genomes were filtered using the CAST compositional bias filtering algorithm , then compared against themselves and each of the other 23 genomes using the BLASTp  sequence similarity searching algorithm with a cut-off E-value of 1 × 10-10. The DifFuse algorithm  was then applied automatically to each genome in turn as a query against the other 23 (reference) genomes. Using other protein databases as reference yields fewer composite cases (for example, the well-known case of the TopA/TopB pair appears multiple times in this analysis), showcasing the extreme bias of annotated databases, such as SwissProt (data not shown). Performing the same computation using the non-redundant sequence database (nrdb) is prohibitively expensive in terms of computation time for an analysis of this size. The detected gene fusion results for each of the 552 comparisons were further automatically filtered for significant overlap of the BLAST alignments of the component proteins. In this case, component proteins that overlap by more than 10% of their total length when aligned together with the composite protein. This step avoids the detection of 'promiscuous domains'  and gene prediction errors, which result in false-positive fusion detection cases. The detected component and composite proteins are far fewer in number than for the two previous reports of E. coli  and S. cerevisiae , due to the much stricter criteria employed in the present analysis and the multi-step protocol we have developed. This analysis was fully automatic and carried out in parallel over a period of four weeks on 11 SUN UltraSPARC CPUs running Solaris 7.
Expression profile analysis
Gene expression ratios for all experiments were transformed into log-odds values so that induction and repression measurements are directly comparable (positive and negative values, respectively). The log-odd values were then normalized across all timepoints for each experiment, using Z-score values for each timepoint. The Z-score values for all time points of each experiment thus allow cross-comparison of gene expression across separate experiments .
Our predicted functional associations for S. cerevisiae with available expression data represent 536 component pairs in total. For each pair of proteins, a Pearson correlation coefficient was calculated between two corresponding experiments and averaged over all experiments. To estimate noise in these data, a control set of 536 randomly selected S. cerevisiae proteins was taken and treated as above (Figure 2).
The distribution of averaged Pearson correlation coefficients for the predicted functional associations was compared against the distribution of coefficients for the control set using a t-test for mean values (where the null hypothesis is that the two means are equal). The test results in a t-value of 3.6 (critical t-value is 1.64), which is highly significant (P-value is 0.000173), indicating that there is a higher average correlation of expression profiles for the predicted functional associations against the background.
All proteins involved in gene fusion events as either component or composite proteins were identified automatically from the results of the fusion analysis. From these data we can obtain raw counts of the number of gene fusion events detected and the number of proteins involved in these events as either composite or component proteins. These figures are skewed however, due to the presence of homology in both the query and reference sets. Proteins involved in gene fusion events as either component or composite genes are then assembled into two lists. These lists are then used to generate two sequence databases, the first one containing all component sequences from the 24 genomes and the second containing all composite sequences.
These sequence databases of component and composite proteins are then compared against themselves using the BLASTp (version 2.0) sequence similarity searching algorithm  (cut-off E-value 1 × 10-10), then clustered according to their similarity using the RAGE algorithm . The RAGE algorithm lists all composite and component proteins in clusters according to similarity and domain structure. Homologous proteins with similar domain structure were clustered together. Each cluster in this case indicates a distinct class of fusion event and cluster members indicate which proteins involved in this type of event from different genomes. These clusters are used to calculate the number of unique fusions detected within and across genomes. This is done by examining how many distinct types of fusion are present in any given genome.
All composite proteins were clustered into 621 families and a distance measure δ was derived according to the sharing of clusters between the 24 species examined. This pairwise distance measure is calculated as δ = (1 - SA,B/TA,B) × 100, where SA,B is the number of shared composite clusters and TA,B is the average of the composite cluster counts from the two species. This measure is reminiscent of a recent genome-wide "ortholog" analysis . This measure was used to calculate a nearest-neighbor dendrogram for the 24 species. Bootstrap values were generated using a 'delete-half jack-knife procedure.
All results of the present analysis are available from the Computational Genomics Group website .
We thank John Aach (Harvard Medical School), Despina Alexandraki (University of Crete and IMBB, Heraklion) and members of the Computational Genomics Group at the EBI for discussions. This work was fully supported by the European Molecular Biology Laboratory (EMBL). C.O. acknowledges further support from the European Commission DGXII (Science, Research and Development), the Medical Research Council (UK) and IBM Research. Patent application filed on behalf of EMBL.
- Enright AJ, Iliopoulos I, Kyrpides NC, Ouzounis CA: Protein interaction maps for complete genomes based on gene fusion events. Nature. 1999, 402: 86-90. 10.1038/47056.PubMedView ArticleGoogle Scholar
- Marcotte EM, Pellegrini M, Thompson MJ, Yeates TO, Eisenberg D: A combined algorithm for genome-wide prediction of protein function. Nature. 1999, 402: 83-86. 10.1038/47048.PubMedView ArticleGoogle Scholar
- Marcotte EM, Pellegrini M, Ng H-L, Rice DW, Yeates TO, Eisenberg D: Detecting protein function and protein-protein interactions from genome sequences. Science. 1999, 285: 751-753. 10.1006/bbrc.2001.5221.PubMedView ArticleGoogle Scholar
- Sali A: Functional links between proteins. Nature. 1999, 402: 23-26. 10.1016/S0027-5107(97)00278-9.PubMedView ArticleGoogle Scholar
- Doolittle RF: Do you dig my groove?. Nat Genet. 1999, 23: 6-8. 10.1038/12597.PubMedView ArticleGoogle Scholar
- Overbeek R, Fonstein M, D'Souza M, Pusch GD, Maltsev N: The use of gene clusters to infer functional coupling. Proc Natl Acad Sci USA. 1999, 96: 2896-2901. 10.1073/pnas.96.6.2896.PubMedPubMed CentralView ArticleGoogle Scholar
- Ito T, Tashiro K, Muta S, Ozawa R, Chiba T, Nishizawa M, Yamamoto K, Kuhara S, Sakaki Y: Toward a protein-protein interaction map of the budding yeast: a comprehensive system to examine two-hybrid interactions in all possible combinations between the yeast proteins. Proc Natl Acad Sci USA. 2000, 97: 1143-1147. 10.1073/pnas.97.3.1143.PubMedPubMed CentralView ArticleGoogle Scholar
- Enright AJ, Ouzounis CA: GeneRAGE: a robust algorithm for sequence clustering and domain detection. Bioinformatics. 2000, 16: 451-457. 10.1093/bioinformatics/16.5.451.PubMedView ArticleGoogle Scholar
- DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278: 680-686. 10.1126/science.278.5338.680.PubMedView ArticleGoogle Scholar
- Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, Lockshon D, Narayan V, Srinivasan M, Pochart P, et al: A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. Nature. 2000, 403: 623-627. 10.1038/35001009.PubMedView ArticleGoogle Scholar
- Lim AL, Powers-Lee SG: Requirement for the carboxyl-terminal domain of Saccharomyces cerevisiae carbamoyl-phosphate synthetase. J Biol Chem. 1996, 271: 11400-11409. 10.1074/jbc.271.19.11400.PubMedView ArticleGoogle Scholar
- Blaiseau PL, Isnard AD, Surdin-Kerjan Y, Thomas D: Met31p and Met32p, two related zinc finger proteins, are involved in transcriptional regulation of yeast sulfur amino acid metabolism. Mol Cell Biol. 1997, 17: 3640-3648.PubMedPubMed CentralView ArticleGoogle Scholar
- Cho RJ, Campbell MJ, Winzeler EA, Steinmetz L, Conway A, Wodicka L, Wolfsberg TG, Gabrielian AE, Landsman D, Lockhart DJ, Davis RW: A genome-wide transcriptional analysis of the mitotic cell cycle. Mol Cell. 1998, 2: 65-73.PubMedView ArticleGoogle Scholar
- Chu S, DeRisi J, Eisen M, Mulholland J, Botstein D, Brown PO, Herskowitz I: The transcriptional program of sporulation in budding yeast [Erratum: Science 282:1421]. Science. 1998, 282: 699-705. 10.1006/jmbi.1998.2134.PubMedView ArticleGoogle Scholar
- Moskovitz J, Berlett BS, Poston JM, Stadtman ER: The yeast peptide-methionine sulfoxide reductase functions as an antioxidant in vivo. Proc Natl Acad Sci USA. 1997, 94: 9585-9589. 10.1073/pnas.94.18.9585.PubMedPubMed CentralView ArticleGoogle Scholar
- Lescure A, Gautheret D, Carbon P, Krol A: Novel selenoproteins identified in silico and in vivo by using a conserved RNA structural motif. J Biol Chem. 1999, 274: 38147-38154. 10.1074/jbc.274.53.38147.PubMedView ArticleGoogle Scholar
- Helliwell SB, Howald I, Barbet N, Hall MN: TOR2 is part of two related signaling pathways coordinating cell growth in Saccharomyces cerevisiae. Genetics. 1998, 148: 99-112.PubMedPubMed CentralGoogle Scholar
- Stoldt V, Rademacher F, Kehren V, Ernst JF, Pearce DA, Sherman F: The Cct eukaryotic chaperonin subunits of Saccharomyces cerevisiae and other yeasts. Yeast. 1996, 12: 523-529. 10.1002/(SICI)1097-0061(199605)12:6<523::AID-YEA962>3.3.CO;2-3.PubMedView ArticleGoogle Scholar
- Kyrpides NC: Genomes OnLine Database (GOLD 1.0): a monitor of complete and ongoing genome projects worldwide. Bioinformatics. 1999, 15: 773-774. 10.1093/bioinformatics/15.9.773.PubMedView ArticleGoogle Scholar
- Promponas VJ, Enright AJ, Tsoka S, Kreil DP, Leroy C, Hamodrakas S, Sander C, Ouzounis CA: CAST: an iterative algorithm for the complexity analysis of sequence tracts. Complexity analysis of sequence tracts. Bioinformatics. 2000, 16: 915-922. 10.1093/bioinformatics/16.10.915.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Snel B, Bork P, Huynen MA: Genome phylogeny based on gene content. Nature Genet. 1999, 21: 108-110. 10.1038/5052.PubMedView ArticleGoogle Scholar
- Computational Genomics Group. [http://www.ebi.ac.uk/research/cgg/diffuse]