Extraction of phylogenetic network modules from the metabolic network
© Yamada et al; licensee BioMed Central Ltd. 2006
Received: 22 May 2005
Accepted: 13 March 2006
Published: 13 March 2006
In bio-systems, genes, proteins and compounds are related to each other, thus forming complex networks. Although each organism has its individual network, some organisms contain common sub-networks based on function. Given a certain sub-network, the distribution of organisms common to it represents the diversity of its function.
We extracted such "common" sub-networks, defined as "phylogenetic network modules," using phylogenetic profiles and cluster analysis. The enzymes in the same "phylogenetic network module" have similar phylogenetic profiles and related functions. These modules are shown to be phylogenetic building blocks. Furthermore, the network of the modules illustrated hierarchical feature as well as the network of enzymes involved in the metabolism.
We conclude that phylogenetic network modules are evolutionary conserved functional units in the metabolic network. We claim that our concept of phylogenetic modules provides a more accurate understanding of the evolution of biological networks.
Bio-systems are constructed by various relationships between genes, compounds and proteins. These elements are intertwined, thus forming a complex network. Such genetic or physical associations contribute to the complicated cellular functions.
In order to characterize whole networks or to find unknown interactions between elements (gene/compound/protein), many research groups have applied the integration of qualitatively different interactions. For example, Rison et al attempted to find unknown functions of genes using genome location, sequence similarity, and relative position in the metabolic network , and Yamanishi et al similarly attempted to apply a kernel method . Such integrative methods infer common features and new categories of genes. Classically, such new categories were determined manually using expert knowledge. The KEGG metabolic pathway  is one such example. On the other hand, recent studies have attempted to define new categories automatically using comprehensive data sets such as microarray data and metabolic networks . Such analyses can classify elements by globally assessing many relationships. Many research groups advocate various functional categories of genes and proteins, such as transcriptional clusters [5, 6], protein complexes , metabolic network [3, 8], and others [9, 10].
From the perspective of graph theory, it is suggested that networks composed of genes and proteins should have modularity. Ravasz et al suggested a hierarchical network model incorporating modularity in the metabolic network . Following this work, studies were attempted using the modularity in various networks and topological features . However, these topological investigations were not based on phylogeny; Snel et al pointed out problems in which they found conflicts between evolution and function .
In order to resolve these conflicts, the main purpose of our research is the extraction of the phylogenetic primary units of the metabolic network for understanding the evolutionary process. To our knowledge, we were the first to attempt to extract evolution-based enzyme modules. In our previous work, we attempted to extract phylogenetic units using prokaryotic genomes . Here, we expanded the analysis by adding eukaryotic data, and by considering the topological properties of the network of phylogenetic units.
The metabolic network was considered not only to be composed of several particular categories, but also to be a large network consisting of chemical compounds connected by enzymes. For example, KEGG/PATHWAY is a collection of manually drawn pathway maps. However, we can regard all of these pathway maps as a "metabolic network." Our analysis is based on the KEGG/PATHWAY database, which integrates the metabolic network from many organisms. Therefore, in this paper, we use the term "metabolic network" to refer to this integrated metabolic network. Therefore, some parts of the metabolic network are found only in particular organisms, and the number of organisms varies depending on the sub-networks. The phylogenetic profile is very useful to deal with such diversity of enzymes. The phylogenetic profile of an enzyme is the string that encodes the presence or absence of the enzyme in the fully sequenced genome .
Pathway distance and the similarity of phylogenetic profiles
Distribution of the number of enzymes included in a module
The white points in Figure 2-B indicate the frequency of the "phylogenetic network modules" according to their sizes (the number of enzymes). The horizontal axis represents the number of enzymes in a module, and the vertical axis represents its frequency. In total, 1179 modules were extracted, but over 900 modules contain only one enzyme. There are a few large modules, and many small modules. Similar results were obtained using different distance measures such as correlation coefficient and Hamming distance (data not shown).
Furthermore, in order to prove that the enzymes form module structures in the metabolic network, we compared the frequency with the average of 100 randomly relabeled networks, shown as the black points in Figure 2-B. The difference of slope indicates that larger modules tend to appear more frequently in actual networks; that is, they hardly appear in the average of relabeled networks. This implies that enzymes with similar phylogenetic profiles tend to aggregate in the network, and that phylogenetic module structures are indeed present in the metabolic network.
Comparing modules with the categories in the KEGG/PATHWAY database
KEGG/PATHWAY is a collection of manually drawn diagrams called the KEGG reference pathway diagrams (maps), each corresponding to a known network of functional significance. These diagrams are drawn manually representing particular functions of biological processes. We mapped our phylogenetic network modules onto these diagrams.
Module 1 is contained in many organisms, spreading to several amino acid biosynthesis pathways such as histidine metabolism (lower part of Figure 3). Module 2 is specialized to prokaryotes, so it links to the prokaryote-specific Peptidoglycan biosynthesis pathway. Module 3 is also contained in relatively many organisms, and module 4 is specialized for eukaryotes. Almost all of the diagrams in KEGG/PATHWAY include multiple phylogenetic network modules, and some modules spread to many diagrams. In fact, the lysine biosynthesis pathway contains several modules, and module 1 and 2 spread to other diagrams. Thus, using our phylogenetic network modules, we could detect phylogenetic relationships between known functional categories (across pathways).
Global network of phylogenetic network modules
Three features are conspicuous in this figure. The first is a large central module. All of the modules seem to gather around the largest red module as a network core. The big module contains 51 enzymes that spread to a large part of nucleotide metabolism, and a part of amino acid metabolism. Furthermore, this is contained in almost all organisms so that it is not surprising that the module constructs a core of the metabolic network. This module connects other modules with each other and sustains the whole network. Additionally, relatively large modules are also observed in parts of the metabolic network: some glycan related pathway diagrams for eukaryotes and a part of peptidoglycan biosynthesis for bacteria.
As a second feature, modules also aggregate. In Figure 4, some of the modules tend to assemble together into clusters. The dotted circles in the figure indicate the locations of what we call "super-modules," which are modules that are relatively aggregated. The relationship between modules in this figure is based on the connectivity of existing enzymes in the metabolic network, so the "super-modules" represent modules of related functions even though some spread to multiple pathway diagrams (only representative diagrams are annotated in Figure 4). For example, the super module structure just above the network core (the biggest red module) in Figure 4 is composed of many diagrams of Amino acid metabolism, and it contains many links to other modules. Basically, the enzymes in this super-module are highly conserved and are closely related to those in the network core. Thus, we claim that enzymes in these modules use many metabolites synthesized in the core structure and that they also supply many metabolites to the other modules. The third feature is that "linker" modules are scattered around the network. As a topological feature, they have relatively low clustering coefficients compared to the other modules with the same degree (see Methods section). Therefore they do not belong to any particular module cluster, and instead, link module clusters to each other. We found 30 linker modules by the criterion defined in the Methods section. Biologically, linker modules tend to be intermediates for the input and output compounds of functional modules. For example, the linker module indicated by the right arrow in Figure 4 connects Amino Acid Metabolism, Lipid Metabolism, as well as a few others. This module is composed of three enzymes (3-hydoroxyacyl-CoA dehydrogenase, acyl-CoA dehydorogenase and enoyl-CoA hydrase) that catalyze reactions between acetoacetyl-CoA and crotonyl-CoA. This reaction chain produces acetyl-CoA in Amino Acid Metabolism and consumes it in Lipid Metabolism, thus linking these two pathways together. As another example, the module indicated by the left arrow in Figure 4 contains three enzymes, glucose-6-phosphate isomerase, phosphoglucomutase and glucose-6-phosphate 1-dehydrogenase. Although these enzymes play roles in Glycolysis, they connect compounds to other Carbohydrate Metabolism pathways such as Aminosugars metabolism and Fructose and mannose metabolism.
Functionally related phylogenetic module
Known functional modules do not completely coincide with evolutionary modules. Snel et al. investigated whether known functional modules are also evolutionary modules and suggested that all the members of the same functional module do not have co-evolutionary tendencies . This means that evolutionary modules are not in complete agreement with functional modules.
Basically, we agree with Snel's opinion because of the difference of enzyme distributions in the phylogenetic network modules due to the addition of enzyme connectivity in the metabolic network. Before adding information of enzyme connectivity, enzyme clusters were constructed using only the similarity between phylogenetic profiles (Fig. 1B). After that, enzymes were re-clustered within each cluster using connectivity in the metabolic network (Fig. 1C,D,E). We call these sub-clusters phylogenetic network modules. Obviously, the number of enzymes in a phylogenetic network module was less than those in the original clusters. It is clear that enzyme connectivity in the metabolic network subdivides the enzyme clusters. This explains precisely the conflict between functional modules and phylogenetic modules. By definition, the phylogenetic network modules are the evolutionarily conserved and functionally related enzyme modules in the metabolic pathway. As a result, we claim that the phylogenetic network module is the basic functional unit in the metabolic pathway.
Modularity and hierarchy in the evolutionary process of the metabolic network
Historically, there have been many analyses of network evolution , which has led to advanced theories on network evolution. Two major assumptions are generally thought to be the main contenders. One is a retrograde model , and the other is a patchwork model . In the retrograde model, network evolves "backwards" from a key metabolite. This model expands the network due to the acquirement of new enzymes, which synthesize a molecule used up in the environment from other molecules. On the other hand, in the patchwork model, network evolution is based on the concept that enzymes exhibit broad substrate specificity and catalyze multiple reactions. Such enzymes with broad specificity form reaction chains to a key metabolite. The important thing is that those two models are not mutually exclusive, and they are reviewed in .
In any case, these evolutionary models are based on the relationship between enzymes and substrates (metabolites). In this paper, phylogenetic network module corresponds to one function which is a group of chemical reactions catalyzing metabolites into others. Correspondingly, we consider the enzyme module as an extended enzyme function in these models. We claim that the concept of representing multiple nodes as a single node, the enzyme module in this case, is important for understanding the evolutionary process of metabolic network. A similar concept for the network integrating other types of relations was proposed in . Our methodology is based on this concept, and our results support the utility of it.
Future direction and perspective on network modularity
We hierarchically clustered enzymes to extract phylogenetic network modules. According to this method, enzymes are allocated to a particular module. However, there are some cases where enzymes could belong to multiple functional modules. For example, there are many enzymes catalyzing reactions related with pyruvate or acetyle-CoA. These enzymes have numerous relationships to other enzymes in the metabolic network. In such cases, the phylogenetic relationships between them are complicated, and enzyme allocation to a particular module is difficult. While these enzymes play an important role, it is difficult to divide into a particular module. Given this problem, it may be important to somehow allow enzyme redundancy in the modules, or remove these enzymes from the network. In a module network, linker modules correspond to this case. These modules have a characteristic feature of connecting many modules (module groups), so it will be difficult to determine particular module groups for linker modules.
Although our work focused on enzymes, the components of the metabolic network consist of both enzymes and chemical compounds. There have been a few attempts to investigate the relationships between chemical compounds and phylogeny. Hattori et al defined the similarity between compounds in the metabolic network, and generated clusters according to this similarity measure . They also attempted to compare these clusters with operon like structures. However, since the operon data they used was very limited, it was insufficient to attempt phylogenetic analysis. Our phylogenetic network modules are suitable for the analysis, and comprehensive analysis of chemical structures and their evolution is our next research focus.
We extracted "phylogenetic network modules" from the metabolic network. We claim that these modules are the evolutionary building blocks as well as the basic functional units of the metabolic network. Furthermore, we showed that the module network has the hierarchical character, which is also conserved in the enzyme network of metabolism.
Barabasi et al illustrated the hierarchical structure of the metabolic network using topological property. Our result suggests that the relationship between enzymes based on the similarity of phylogenetic profiles is one of the factors forming the hierarchical structure. Other biological systems, such as protein-protein interaction networks, have been reported to have the similar topological property. Thus, the phylogenetic relationship may be a foundation of network evolution including other biological systems.
Construction of the phylogenetic profiles and the network of enzymes
A phylogenetic profile is a bit string that encodes the absence (0) or presence (1) of an enzyme in fully sequenced genomes. We utilized KEGG Orthology (KO) for constructing the phylogenetic profiles. KO is a database of ortholog groups that are defined manually according to the similarity of amino acid sequences, as well as bidirectional best hit information in pairwise genome comparison, and annotated functions of genes in KEGG pathways . Therefore, some proteins with low sequence similarity may be put into the same ortholog group (the same node in the KEGG pathway). In this study, we used KO entries for enzyme genes annotated in the KEGG/PATHWAY database. Phylogenetic profiles were calculated using the KO entries constructed from 174 fully sequenced genomes. To reduce the effect of bias in the organism distribution, these 174 genomes were merged into 59 taxa (11 eukaryotes, 36 bacteria and 12 archaea) according to the NCBI taxonomy . That is, logical OR was applied on the phylogenetic profiles of the organisms of the same taxa. As a result, we obtained 1672 phylogenetic profiles consisting of 59 bits.
The KEGG/PATHWAY database stores metabolic and regulatory pathway information with their functional classification as XML files . The pathway information (the network of enzymes) is described as a collection of binary relationships of enzymes. The binary relationship is defined when the two enzymes are adjacently located on the pathway diagram. We extracted the binary relationships of the enzymes from the KEGG/PATHWAY database and treated them as a network or a graph of enzymes. All the data is accessible from the KEGG website .
Similarity measure between phylogenetic profiles
Similarity measures between phylogenetic profiles are required in cluster analysis. We adopted the Jaccard coefficient (JC)  and the correlation coefficient (CC) as the similarity measures of phylogenetic profiles. The Jaccard coefficient between profile A and B is defined as A ∩ B/A ∪ B. A ∩ B is the number of organisms which include enzyme A and B, and A ∪ B is the number of organisms which include either enzyme A or B. The correlation coefficient used was Pearson's product moment correlation coefficient using a phylogenetic profile as a vector.
Pathway distance and the corresponding average of Jaccard coefficient
Pathway distance is defined as the smallest number of steps between two enzymes in the metabolic network. We calculated the pathway distance of all against all enzymes in KEGG/PATHWAY and further calculated the JC average for each pathway distance. We also generated 100 relabeled enzyme networks as a control. In this paper, we define a relabeled network as one where the phylogenetic profiles are randomly re-assigned to enzymes in the network with the topology unchanged. The JC average of each pathway distance was calculated for each relabeled network. The average was then taken of all the JC averages of each relabeled network.
Extraction of the "phylogenetic network module"
An overview of the extraction of the "phylogenetic network module" is illustrated in Figure 1. We performed a complete linkage clustering of the enzymes (KO entries) in our dataset (Fig. 1B) based on their phylogenetic profile similarities. To determine the threshold of this clustering, we utilized the distribution of all-against-all similarities between the phylogenetic profiles of enzymes. From the distribution, we estimated the top 1, 2.5 and 5 percentile points of the similarities as the significant point (each JC corresponds to 0.76, 0.64 and 0.51 respectively, see additional files 2-3). After that, by linking the enzymes according to their connection in the metabolic pathway (Fig. 1C), some "phylogenetic network modules" were extracted in each cluster (Fig. 1D). That is, each group of linked enzymes within each cluster was extracted as a single module. There was no difference in the distribution shape of the module sizes for different significant points. All figures in this paper are based on the 5 percentile point. The same operations were applied to the 100 relabeled networks, and their distribution of the average number of enzymes per module was obtained.
Network of modules
Finally, we constructed a network of modules, where each node represents a module, and each edge represents the chemical compounds between the enzymes contained in the corresponding modules (Fig. 1E). To infer the structure of the module network, we investigated the relationship between each node's degree and its clustering coefficient. The degree of node i is the number of connected nodes. The clustering coefficient of node i is defined as C i = 2n i /k i (k i -l), where n i denotes the number of links connecting the k i neighbors of node i to each other.
Based on the clustering coefficient and the degree, linker modules were selected. These linker modules tend to have high degree and relatively low clustering coefficients compared with other modules of the same degree. We defined the condition of the linker as having degree > 35, and the deviation of the clustering coefficient of the same degree being < -0.1.
Authors are grateful to Dr. K. F. Aoki Kinoshita and Dr. C. Wheelock for helpful comments on an earlier draft of the manuscript. This work was supported by grants from the Ministry of Education, Culture, Sports, Science and Technology of Japan, the Japan Society for Promotion of Science and the Japan Science and Technology Agency. Computational resource was provided by the Bioinformatics Center, Institute for Chemical Research, Kyoto University.
- Rison SC, Teichmann SA, Thornton JM: Homology, pathway distance and chromosomal localization of the small molecule metabolism enzymes in Escherichia coli. J Mol Biol 2002, 318(3):911–932. 10.1016/S0022-2836(02)00140-7View ArticlePubMedGoogle Scholar
- Yamanishi Y, Vert JP, Kanehisa M: Protein network inference from multiple genomic data: a supervised approach. Bioinformatics 2004, 20(Suppl 1):I363–1370. 10.1093/bioinformatics/bth910View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res 2004, (32 Database):D277–280. 10.1093/nar/gkh063
- DeRisi J, Penland L, Brown PO, Bittner ML, Meltzer PS, Ray M, Chen Y, Su YA, Trent JM: Use of a cDNA microarray to analyse gene expression patterns in human cancer. Nat Genet 1996, 14(4):457–460. 10.1038/ng1296-457View ArticlePubMedGoogle Scholar
- Ihmels J, Friedlander G, Bergmann S, Sarig O, Ziv Y, Barkai N: Revealing modular organization in the yeast transcriptional network. Nat Genet 2002, 31(4):370–377.PubMedGoogle Scholar
- Salgado H, Santos-Zavaleta A, Gama-Castro S, Millan-Zarate D, Diaz-Peredo E, Sanchez-Solano F, Perez-Rueda E, Bonavides-Martinez C, Collado-Vides J: RegulonDB (version 3.2): transcriptional regulation and operon oganization in Escherichia coli K-12. Nucleic Acids Res 2001, 29(1):72–74. 10.1093/nar/29.1.72PubMed CentralView ArticlePubMedGoogle Scholar
- Gavin AC, Bosche M, Krause R, Grandi P, Marzioch M, Bauer A, Schultz J, Rick JM, Michon AM, Cruciat CM, Remor M, Hofert C, Schelder M, Brajenovic M, Ruffner H, Merino A, Klein K, Hudak M, Dickson D, Rudi T, Gnau V, Bauch A, Bastuck S, Huhse B, Leutwein C, Heurtier MA, Copley RR, Edelmann A, Querfurth E, Rybin V, Drewes G, Raida M, Bouwmeester T, Bork P, Seraphin B, Kuster B, Neubauer G, Superti-Furga G: Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature 2002, 415(6868):141–147. 10.1038/415141aView ArticlePubMedGoogle Scholar
- 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(1):56–58. 10.1093/nar/30.1.56PubMed CentralView ArticlePubMedGoogle Scholar
- Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de sur la Cruz N, Tonellato P, Jaiswal P, Seigfried T, White R: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res 2004, (32 Database):D258–261.
- von Mering C, Zdobnov EM, Tsoka S, Ciccarelli FD, Pereira-Leal JB, Ouzounis CA, Bork P: Genome evolution reveals biochemical networks and functional modules. Proc Natl Acad Sci USA 2003, 100(26):15428–15433. 10.1073/pnas.2136809100PubMed CentralView ArticlePubMedGoogle Scholar
- Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL: Hierarchical organization of modularity in metabolic networks. Science 2002, 297(5586):1551–1555. 10.1126/science.1073374View ArticlePubMedGoogle Scholar
- Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks. Science 2002, 298(5594):824–827. 10.1126/science.298.5594.824View ArticlePubMedGoogle Scholar
- Snel B, Huynen MA: Quantifying modularity in the evolution of biomolecular systems. Genome Res 2004, 14(3):391–397. 10.1101/gr.1969504PubMed CentralView ArticlePubMedGoogle Scholar
- Yamada Y, Goto S, Kanehisa M: Extraction of phylogenetic network modules from prokaryote metabolic pathways. Genome informatics 2004, 15(1):249–258.PubMedGoogle Scholar
- 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(8):4285–4288. 10.1073/pnas.96.8.4285PubMed CentralView ArticlePubMedGoogle Scholar
- Lazcano A, Miller SL: On the origin of metabolic pathways. J Mol Evol 1999, 49(4):424–431. 10.1007/PL00006565View ArticlePubMedGoogle Scholar
- Horowitz HN: On the evolution of biochemical synthesis. Proc Natl Acad Sci U S A 1945, 31: 153–157.PubMed CentralView ArticlePubMedGoogle Scholar
- Jensen RA: Enzyme recruitment in evolution of new function. Annu Rev Microbiol 1976, 30: 409–425. 10.1146/annurev.mi.30.100176.002205View ArticlePubMedGoogle Scholar
- Rison SC, Thornton JM: Pathway evolution, structurally speaking. Cur Opin Struct Biol 2002, 12(3):374–382. 10.1016/S0959-440X(02)00331-7View ArticleGoogle Scholar
- Lee I, Date SV, Adai AT, Marcotte EM: A probabilistic functional network of yeast genes. Science 2004, 306(5701):1555–1558. 10.1126/science.1099511View ArticlePubMedGoogle Scholar
- Hattori M, Okuno Y, Goto S, Kanehisa M: Development of a chemical structure comparison method for integrated analysis of chemical and genomic information in the metabolic pathways. J Am Chem Soc 2003, 125(39):11853–11865. 10.1021/ja036030uView ArticlePubMedGoogle Scholar
- Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Church DM, DiCuccio M, Edgar R, Federhen S, Helmberg W, Kenton DL, Khovayko O, Lipman DJ, Madden TL, Maglott DR, Ostell J, Pontius JU, Pruitt KD, Schuler GD, Schriml LM, Sequeira E, Sherry ST, Sirotkin K, Starchenko G, Suzek TO, Tatusov R, Tatusova TA, Wagner L, Yaschenko E: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res 2005, (33 Database):D39–45.
- Kyoto Encyclopedia of Genes and Genomes[http://www.genome.jp/kegg]
- Wolf YI, Rogozin IB, Kondrashov AS, Koonin EV: Genome alignment, evolution of prokaryotic genome organization, and prediction of gene function using genomic context. Genome Res 2001, 11(3):356–372. 10.1101/gr.GR-1619RView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.