- Research article
- Open Access
- Published:

# Algebraic comparison of metabolic networks, phylogenetic inference, and metabolic innovation

*BMC Bioinformatics*
**volume 7**, Article number: 67 (2006)

## Abstract

### Background

Comparison of metabolic networks is typically performed based on the organisms' enzyme contents. This approach disregards functional replacements as well as orthologies that are misannotated. Direct comparison of the structure of metabolic networks can circumvent these problems.

### Results

Metabolic networks are naturally represented as directed hypergraphs in such a way that metabolites are nodes and enzyme-catalyzed reactions form (hyper)edges. The familiar operations from set algebra (union, intersection, and difference) form a natural basis for both the pairwise comparison of networks and identification of distinct metabolic features of a set of algorithms. We report here on an implementation of this approach and its application to the procaryotes.

### Conclusion

We demonstrate that metabolic networks contain valuable phylogenetic information by comparing phylogenies obtained from network comparisons with 16S RNA phylogenies. The algebraic approach to metabolic networks is suitable to study metabolic innovations in two sets of organisms, free living microbes and *Pyrococci*, as well as obligate intracellular pathogens.

## Background

The metabolic networks of a wide variety of organisms, in particular procaryotes, have been reconstructed by means of a combination of genomic annotations with biochemical and physiological data, see e.g. [1]. These networks are compiled in databases, in particular in the KEGG resource [2].

Large scale bacterial phylogenies that are based on single genes are notoriously plagued by gene transfer, gene duplication, gene deletion, and functional replacement of genes. The same holds for various approaches towards utilizing gene content for phylogenetic purposes, discussed e.g. by [3–8]. A recent article by [9] addressed this issue by considering the presence or absence of 64 individual subpathways that were identified based on the COG division [10] of the National Center for Biotechnology Information. A related approach, based on comparison of individual pathways was discussed by [11–13] and [14]. The pathways necessary for such approaches can be derived from a given metabolic network either "by hand" or using automated procedures such as metabolic flux analysis, see e.g. [15–19].

Instead of attempting to first reconstruct individual pathways, we take here a more global view by grounding our analysis in the direct comparison of the metabolic networks. While the application of generic graph distances or similarity measures (see e.g. [20]) is certainly appealing, they cannot be used in a straightforward manner for metabolic networks. The reason is that chemical reaction networks do not have a simple representation as graphs, at least not when metabolites are represented as nodes and reactions as edges. Instead, a metabolic network is naturally described by a directed hypergraph [21], or, equivalently, by a directed bipartite graph, in which metabolites and reactions (or, equivalently, the enzymes that catalyze the reactions), are represented by two different types of vertices. More global comparison of metabolic networks, in terms of various network indices and networks motifs, can be found in [22].

This contribution is organized as follows: In the next section we summarize an algebraic approach to comparison and manipulation of chemical reaction networks that is motivated by set theory. We briefly describe the C library that implements this approach. We then demonstrate that the symmetric difference of two metabolic networks can be used to derive a distance measure that is suitable for reconstructing phylogenetic relationship from metabolic network data. More interestingly, however, the same approach can be used directly to extract those subnetworks of the metabolism that are innovations in the particular subtree of the phylogeny. We illustrate our approach using pathogenic procaryotes as an example.

### The algebra of directed hypergraphs

A metabolic network is defined by its metabolites and the system of reactions that inter-converts them. We denote the set of metabolites by *X*. A chemical reaction *E* can be described as a pair of multisets (*E*^{-}, *E*^{+}), where *E*^{-} ⊆ *X* is the set of educts in the reaction and *E*^{+} ⊆ *X* is the set of reaction products. Slightly more generally, we can replace the multisets by ordinary sets and instead define the multiplicities of product and educt metabolites by means of the stoichiometric coefficients ${n}_{x,E}^{+}$ and ${n}_{x,E}^{-}$ of the products and educts, respectively. A metabolic network is thus a pair (*X*, $\mathcal{E}$) where $\mathcal{E}$ is a set of reactions. Such a structure is known as a directed hypergraph $\mathfrak{M}$(*X*, $\mathcal{E}$), see e.g. [21]. The *stoichiometric matrix* **S** of the network has the entries

${S}_{xE}={n}_{x,E}^{+}-{n}_{x,E}^{-}\left(1\right)$

where *x* ∈ *X* is a metabolite, and *E* is a reaction. For completeness, we remark that the set *E*^{c} = *E*^{+} ∩ *E*^{-} are the catalysts of the reaction *E*. Furthermore, a reaction is autocatalytic if ${n}_{x,E}^{+}$ - ${n}_{x,E}^{-}$ ≠ 0 for some *x* ∈ *E*^{c}. By abuse of notation we write *E* = *E*^{+} ∪ *E*^{-} for the set of metabolites involved in the reaction *E*. Furthermore, we write supp$\mathcal{E}$ = ∪{*E*|*E* ∈ $\mathcal{E}$} for the set metabolites that actually take part in the reactions. We call a network $\mathfrak{M}$(*X*, $\mathcal{E}$) *clean* if *X* = supp$\mathcal{E}$ and define the *clean up operator* as $\mathfrak{M}$ = (supp$\mathcal{E}$, $\mathcal{E}$). Furthermore, for a given set $\mathcal{E}$ of reactions and set *A* metabolites we define

[*A*] = {*E* ∈ $\mathcal{E}$ | (*E*^{+} ∪ *E*^{-}) ⊆ *A*} (2)

The restriction of a network $\mathfrak{M}$(*X*, *E*) to a set *A* of metabolites is defined as the clean network

[*A*] = (*A*, $\mathcal{E}$[*A*]). (3)

For short we write $\mathfrak{M}$ [$\mathcal{E}$] = $\mathfrak{M}$ [supp$\mathcal{E}$] for the restriction with respect to a set of reactions. The number of reactions in a network $\mathfrak{M}$ will be denoted by ||$\mathfrak{M}$||.

In order to compare networks in a systematic way, we need to be able to determine the differences and the commonalities of two networks. Inspired by the usual algebra of sets, which is based on the operations "union", "intersection", and "difference", we introduce analogous mathematical constructions for chemical reaction networks. Fig. 1 illustrates the basic operations of this "network algebra" which are formally defined below.

In the following let

' (*X'*, ε') and '' (*X''*, ε'') be two networks. Of course we have ' = ''iff *X'* = *X''* and ε' = ε''. The empty network will be denoted by ∅.

#### Union

The union $\mathfrak{M}$ = ${\mathfrak{M}}^{\prime}$ ∪ ${\mathfrak{M}}^{\u2033}$ is defined as the network (*X'* ∪ *X"*, ${\mathcal{E}}^{\prime}$ ∪ ${\mathcal{E}}^{\u2033}$). Note that $\mathfrak{M}$ is clean if both ${\mathfrak{M}}^{\prime}$ and ${\mathfrak{M}}^{\u2033}$ are clean.

#### Intersection

The intersection $\mathfrak{M}$ = ${\mathfrak{M}}^{\prime}$ ∩ ${\mathfrak{M}}^{\u2033}$ is defined as the clean network

= (*X'* ∩ *X"*, ${\mathcal{E}}^{\prime}$ ∩ ${\mathcal{E}}^{\u2033}$) (4)

Note that (${\mathcal{E}}^{\prime}$ ∩ ${\mathcal{E}}^{\u2033}$) [*X'* ∩ *X"*] = ${\mathcal{E}}^{\prime}$ ∩ ${\mathcal{E}}^{\u2033}$.

#### Difference

The *difference* $\mathfrak{M}$ = ${\mathfrak{M}}^{\prime}$\${\mathfrak{M}}^{\u2033}$ is defined as the clean network

= (supp(${\mathcal{E}}^{\prime}$\${\mathcal{E}}^{\u2033}$), ${\mathcal{E}}^{\prime}$\${\mathcal{E}}^{\u2033}$) (5)

The difference network contains all reactions occurring in ${\mathfrak{M}}^{\prime}$ but not ${\mathfrak{M}}^{\u2033}$, and all metabolites occurring in the remaining reactions.

The *strict difference* $\mathfrak{M}$ = ${\mathfrak{M}}^{\prime}$\\${\mathfrak{M}}^{\u2033}$ is the clean network

= (*X'*\*X"*, (${\mathcal{E}}^{\prime}$\${\mathcal{E}}^{\u2033}$), [*X'*\*X"*]) (6)

The new network contains only those metabolites occurring in ${\mathfrak{M}}^{\prime}$ but not ${\mathfrak{M}}^{\u2033}$, and only those reactions from ${\mathfrak{M}}^{\prime}$ that can be performed with the remaining metabolites. Thus, we have ||${\mathfrak{M}}^{\prime}$\\${\mathfrak{M}}^{\u2033}$|| ≤ ||${\mathfrak{M}}^{\prime}$\${\mathfrak{M}}^{\u2033}$||.

#### Symmetric difference

The *symmetric difference* $\mathfrak{M}$ = ${\mathfrak{M}}^{\prime}$ △ ${\mathfrak{M}}^{\u2033}$ is defined as the clean network $\mathfrak{M}$ = (${\mathfrak{M}}^{\prime}$ ∪ ${\mathfrak{M}}^{\u2033}$)\(${\mathfrak{M}}^{\prime}$ ∩ ${\mathfrak{M}}^{\u2033}$).

#### Strict symmetric difference

The *strict symmetric difference* $\mathfrak{M}$ = ${\mathfrak{M}}^{\prime}$ ◇ ${\mathfrak{M}}^{\u2033}$ is $\mathfrak{M}$ = (${\mathfrak{M}}^{\prime}$ ∪ ${\mathfrak{M}}^{\u2033}$)\\(${\mathfrak{M}}^{\prime}$ ∩ ${\mathfrak{M}}^{\u2033}$).

The Vienna Reaction Network Library Vienna-RNL implements these basic set-theoretic operations on chemical reaction networks. It is available under the GNU Public License from [32] and as additional file: 1. The library is written in platform independent ANSI C and provides basic data structures for chemical reactions and their networks, IO routines for reading and writing and various formats, as well as set operations such as the union, intersection, or difference of two chemical reaction networks. It is intended for the use in conjunction with the user's own C programs or PERL scripts.

The library contains IO-Routines for reading SBML [23], an XML based dialect for the standardized representation of systems biology models, writing of SBML is currently being implemented. The capability of reading and writing SBML will make the functionality of the Vienna Reaction Network Library accessible to about 80 other software systems which support SBML [33].

### Phylogenies from networks

Datasets were retrieved from the KEGG database on metabolic networks [2], which holds genomic and network data of about 20 Archaea, 200 Bacteria and 20 Eucarya, where in particular the data of many Eukaryotes are incomplete. In a preparatory step we decomposed the individual KEGG-pathways into their chemical reactions and combined these to a complete network for each organism.

The simplest approach to inferring phylogenetic relationships from metabolic networks is to use a distance measure *d* on the set of reaction networks. We use here

$d({\mathfrak{M}}^{\prime},{\mathfrak{M}}^{\u2033})=\frac{\Vert {\mathfrak{M}}^{\prime}\u25b3{\mathfrak{M}}^{\u2033}\Vert}{\Vert {\mathfrak{M}}^{\prime}\Vert +\Vert {\mathfrak{M}}^{\u2033}\Vert -\Vert {\mathfrak{M}}^{\prime}\cap {\mathfrak{M}}^{\u2033}\Vert}=\frac{\Vert {\mathfrak{M}}^{\prime}\u25b3{\mathfrak{M}}^{\u2033}\Vert}{\Vert {\mathfrak{M}}^{\prime}\cup {\mathfrak{M}}^{\u2033}\Vert}\left(7\right)$

Alternatively, the strict symmetric difference ${\mathfrak{M}}^{\prime}$ ◇ ${\mathfrak{M}}^{\u2033}$could be used to define a difference measure. Furthermore, other normalizations of the difference measure could be used. By calculating tree distances between the 16S rRNA phylogenies and the network phylogenies, using the treedist program from the phylip package [24], we have observed, however, that equ.(7) performs best with respect to reproducing trusted 16S RNA phylogenies. Tree distances between the 16S rRNA phylogeny and the network phylogeny calculated with the symmetric difference are 2 and 0.023, using symmetric difference and branch score distance measure, respectively. Compared with tree distances between the 16S rRNA phylogeny and the network phylogeny utilizing the strict symmetric difference yields a much larger tree distance of 10 and 0.093, respectively.

Distance-based network phylogenies are computed using the Fitch algorithm [25] implemented in the phylip package as well as using the splits-decomposition algorithm from the SplitsTree package [26].

An example comprising a selection of bacterial and archaeal metabolic networks is shown in Fig. 2, see Table 1 for the list of species used. The phylogeny inferred from the metabolic networks conforms almost perfectly with the maximum parsimony tree computed from the 16S rRNA sequences of the same organisms. The rRNA sequences were aligned using clustalx. The minor discrepancies are due to poorly resolved nodes as can be seen in the split-decomposition network below. The congruence of rRNA trees and network-based phylogenies demonstrates that the structure of metabolic networks encodes phylogenetic information. Metabolic networks thus constitute a source of phylogenetic signal that is completely independent of sequence information.

The use of distance measures reduces the available information on the network structure already in the first step. We therefore complement distance based phylogenetic analysis with parsimony methods using reaction content: For a given set of organisms, we first calculated the union $\mathfrak{U}$ of their metabolic networks. For each organism we then constructed a reaction profile reflecting presence or absence of each reaction in $\mathfrak{U}$ in the metabolic network of the respective organism. This approach thus is reduced to reconstructing phylogenies from character-tables that represent the presence/absence of particular reactions in the reaction network. It should be noted that this is similar, but not quite the same, as using the presence or absence of orthologous enzymes (see e.g. [3–8]). The main difference is that the network based approach tolerates functional replacements that may occur, e.g., through horizontal gene transfer [9].

### Metabolic innovations

The algebraic approach to metabolic network evolution can also be used in a straightforward way to trace the history of metabolic innovations. To this end, consider a (trusted) unrooted phylogenetic tree *T* in which each leaf of *T* is labeled with the metabolic network $\mathfrak{M}$_{
k
}of the corresponding taxon *k*. Each edge *e* of *T* defines a split, i.e., a bipartition *σ*_{
e
}= {*U*_{
e
}, _{
e
}} of the set of taxa. Here we regard splits as directed. Note that mathematically we can define innovations at each split in both directions. One of the two subsets *U* or , however, contains the ancestral state, hence only one direction makes biological sense: this is the one where the ancestral state (root of the tree) is located in the sub-set . This knowledge has to be provided externally.

Consider an (arbitrary) directed split *σ* = (*U*, ) on the given set of taxa, i.e., a pair of sets of taxa (*U*, ) such that *U* θ, θ, and *U* ∩ = θ. We define the differential metabolic network

The network $\mathcal{D}$(*σ*) describes the *metabolic innovations* in *U* relative to the "background" $\overline{U}$.

As discussed in the previous section, network phylogenies are rather sensitive with respect to life-style and environmental constraints. The organisms whose metabolic networks have been utilized to compute the 16s rRNA tree shown in Fig. 2 are able to freely live in the environment with a reasonably large capacity for adaptation.

As a first example we analyzed the unique metabolic network from the *Pyrococcus* genus. Figure 3 shows the network phylogeny from Fig. 2 with the *Pyrococcus spp*. clade highlighted. The resulting differential network indicates reactions present in *Pyrococcus spp*. but absent in all other organisms of the phylogeny (Figure 3). For example, reaction *R* 01087 is catalyzed by Maleate cis-trans-isomerase which is utilized in maleate assimilating and high-temperature bacteria. A second sub-network involving both ADP-forming acetate and propanoate CoA ligases is potentially used in the organisms to convert between acetate and propanoate and their corresponding CoA forms.

As a second example we analyzed our set of reference organism (Figure 2) with obligatory intracellular pathogens. Figure 4 shows the phylogeny with the selected pathogens (dashed oval). Interestingly, *Mollicutes*, such as *Mycoplasmae* and *Spirochaetes*, such as *Treponema* are grouped together. They all possess a minimal gene-set, and thus a highly optimized and host-dependent metabolic network. Surprisingly, this set of organisms has specific reactions that are absent in the remaining organisms of the phylogeny. Figure 5 shows the corresponding differential network which consists of five sub-network. The two largest networks involve sugar-conversions and parts of glycolysis. Smaller networks correspond to formylation of tetrahydrofolate as well as cholin and carnitine pathways.

## Discussion

The Vienna Reaction Network Library introduced above treats chemical reaction networks, and metabolic networks in particular, as directed hypergraphs. A framework borrowed from set algebra provides natural definitions of unions, intersections, and differences that can be used to compare the metabolic networks of different organisms. We have demonstrated that metabolic networks convey phylogenetic information and can indeed be used to infer phylogenetic relationships of free-living organisms in a way that is similar to gene-content based approaches. In contrast to the latter, however, metabolic network based phylogenies are less sensitive to the effects of functional replacement, e.g., through horizontal gene transfer.

Differences of metabolic networks among subtrees of a trusted phylogeny, or more generally, along any split of interest in a set of organisms can be computed directly, making it easy to study metabolic innovations in particular clades. A first application of our network phylogeny analysis involved three members of the *Pyrococcus spp*. clade. The metabolic reactions resulting from the split between the *Pyrococci* and the remaining organisms involve the maleate cis-trans-isomerase reaction, ADP-forming acetate and propanoate CoA ligase reactions as well as beta-D-Glucose:NAD(P)^{+} 1-oxoreductase.

Our second example considers a class of intra-cellular pathogens that includes *Mycoplasmae, Ureaplasmae*, and *Spirochete*. Their restricted repertoire of metabolic reactions reflect the specialized life-style. Many metabolic pathways are not required in such a rich environment and have been lost in the course of evolution. On the other hand, constructing a network phylogeny including these microbes, we observe metabolic reactions assembling an unconnected network that is present in this set of intracellular pathogens and absent in the remaining organisms. Such reactions include phosphorilization and conversions of sugars and derivates, deaminating lyase reactions, and reactions involving carnitine, choline and tetrahydrofolate.

At present, metabolic network data are compiled by a multitude of methods, and at least in part are constructed by genomic similarity with other organisms. Strictly speaking, therefore, we cannot view metabolic network data such as those complied in the KEGG database as independent from genomic data. With the recent advances of experimental techniques in metabolomics (see e.g. [27–29]), however, the situation is rapidly improving.

## Conclusion

Our comparative approach to metabolic network analysis, which focuses on individual reactions rather than on aggregate features such as pathways, simplifies the identification of metabolic innovations and, in particular, facilitates the recognition of organisms as potential biological threat agents based on their metabolic repertoire. Furthermore, the ability to easily identify differences in metabolic capacity between pathogens should be useful for devising a refined classification of pathogenicity based on metabolic capabilities.

In this contribution we have restricted ourselves to unweighted networks. Distance measures between networks, however, could be refined by attaching weights to both vertices and (hyper-)edges without requiring significant algorithmic changes. These could reflect, e.g., how essential a reaction or a metabolite is for each organism. With the increasing amount and accuracy of available data it might also be feasible to devise a stochastic model of the evolution of metabolic networks, which could then be turned into a scoring scheme for a generalized version of (local) graph alignment along the lines of [30].

## References

- 1.
Becker SA, Palsson BO:

**Genome-scale reconstruction of the metabolic network in***Staphylococcus aureus N315:***an initial draft to the two-dimensional annotation.***BMC Microbiology*2005,**5:**8. [epub] 10.1186/1471-2180-5-8 - 2.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M:

**The KEGG resource for deciphering the genome.***Nucl Acids Res*2004,**32:**D277-D280. 10.1093/nar/gkh063 - 3.
Fitz-Gibbon ST, House CH:

**Whole genome-based phylogenetic analysis of free-living microorganisms.***Nucleic Acids Res*1999,**27:**4218–4222. 10.1093/nar/27.21.4218 - 4.
Ma HW, Zeng AP:

**Phylogenetic comparison of metabolic capacities of organisms at genome level.***Mol Phylogenet Evol*2004,**31:**204–213. 10.1016/j.ympev.2003.08.011 - 5.
Snel B, Bork P, Huynen MA:

**Genome phylogeny based on gene content.***Nat Genet*1999,**21:**108–110. 10.1038/5052 - 6.
Snel B, Bork P, Huynen MA:

**Genomes in flux: the evolution of archaeal and proteobacterial gene content.***Genome Res*2002,**12:**17–25. 10.1101/gr.176501 - 7.
Wolf YI, Rogozin IB, Grishin NV, Tatusov RL, Koonin EV:

**Genome trees constructed using five different approaches suggest new major bacterial clades.***BMC Evol Biol*2001,**1:**8. [epub] 10.1186/1471-2148-1-8 - 8.
Yang S, Doolittle RF, Bourne PE:

**Phylogeny determined by protein domain content.***Proc Natl Acad Sci USA*2005,**102:**373–378. 10.1073/pnas.0408810102 - 9.
Hong SH, Kim TY, Lee SY:

**Phylogenetic analysis based on genome-scale metabolic pathway reaction content.***Appl Microbiol Biotechnol*2004,**65:**203–210. 10.1007/s00253-004-1641-3 - 10.
Tatusov RL, Koonin EV, Lipman DJ:

**A Genomic Perspective on Protein Families.***Science*1997,**278:**631–637. 10.1126/science.278.5338.631 - 11.
Dandekar T, Schuster S, Snel B, Huynen MA, Bork P:

**Pathway Alignment: application to the comparative analysis of glycolytic enzymes.***Biochem J*1999,**343:**115–124. 10.1042/0264-6021:3430115 - 12.
Forst CV, Schulten K:

**Evolution of Metabolism: A New Method for the Comparison of Metabolic Pathways Using Genomic Information.***J Comp Biol*1999,**6:**343–360. 10.1089/106652799318319 - 13.
Forst CV, Schulten K:

**Phylogenetic analysis of metabolic pathways.***J Mol Evol*2001,**52:**471–489. - 14.
Heymans M, Singh AK:

**Deriving phylogenetic trees from the similarity analysis of metabolic pathways.***Bioinformatics*2003,**19**(Suppl 1):i138-i346. 10.1093/bioinformatics/btg1018 - 15.
Schilling CH, Palsson B∅:

**The underlying pathway structure of biochemical reaction networks.***Proc Natl Acad Sci USA*1998,**95:**4193–4198. 10.1073/pnas.95.8.4193 - 16.
Schilling CH, Letscher D, Palsson B∅:

**Theory for the systematic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective.***J Theor Biol*2000,**203:**229–248. 10.1006/jtbi.2000.1073 - 17.
Schuster S, Fell DA, Dandekar T:

**A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks.***Nature Biotechnol*2000,**18:**326–332. 10.1038/73786 - 18.
Xiong M, Zhao J, Xiong H:

**Network-based regulatory pathways analysis.***Bioinformatics*2004,**20:**2056–2066. 10.1093/bioinformatics/bth201 - 19.
Gagneur J, Klamt S:

**Computation of elementary modes: a unifying framework and the new binary approach.***BMC Bioinformatics*2004,**5:**175. [epub] 10.1186/1471-2105-5-175 - 20.
Bunke H, Shearer K:

**On a relation between graph edit distance and maximum common subgraph.***Pattern Rec Let*1997,**18:**689–694. 10.1016/S0167-8655(97)00060-3 - 21.
Zeigarnik AV:

**On Hypercycles and Hypercircuits in Hypergraphs.**In*Discrete Mathematical Chemistry, DIMACS series in discrete mathematics and theoretical computer science*.*Volume 51*. Edited by: Hansen P, Fowler PW, Zheng M. Providence, RI: American Mathematical Society; 2000:377–383. - 22.
Zhu D, Qin ZS:

**Structural comparison of metabolic networks in selected single cell organisms.***BMC Bioinformatics*2005,**6:**8. [epub] [Doi:10.1186/1471–2105–6-8] 10.1186/1471-2105-6-8 - 23.
Hucka M, Finney A, Bornstein BJ, Keating SM, Shapiro BE, Matthews J, Kovitz BL, J SM, Funahashi A, Doyle JC, Kitano H:

**Evolving a lingua franca and associated software infrastructure for computational systems biology: the Systems Biology Markup Language (SBML) project.***Syst Biol*2004,**1:**41–52. [Doi: 10.1049/sb:20045008] 10.1049/sb:20045008 - 24.
Felsenstein J:

**Inferring phylogenies from protein sequences by parsimony, distance, and likelihood methods.***Methods Enzymol*1996,**266:**418–427. - 25.
Fitch W, Margoliash E:

**Construction of phylogenetic trees.***Science*1967,**155:**279–284. - 26.
Huson DH:

**SplitsTree: analyzing and visualizing evolutionary data.***Bioinformatics*1998,**14:**68–73. 10.1093/bioinformatics/14.1.68 - 27.
Brown SC, Kruppa G, Dasseux JL:

**Metabolomics applications of FT-ICR mass spectrometry.***Mass Spectrom Rev*2005,**24:**223–231. 10.1002/mas.20011 - 28.
Griffin JL:

**Metabolic profiles to define the genome: can we hear the phenotypes?***Philos Trans R Soc Lond B Biol Sci*2004,**359:**857–871. 10.1098/rstb.2003.1411 - 29.
Sumner LW, Mendes P, Dixon RA:

**Plant metabolomics: large-scale phytochemistry in the functional genomics era.***Phytochemistry*2003,**62:**817–836. 10.1016/S0031-9422(02)00708-2 - 30.
Berg J, Lässig M:

**Local graph alignment and motif search in biological networks.***Proc Natl Acad Sci USA*2004,**101:**14689–14694. 10.1073/pnas.0305199101 - 31.
Huynen MA, Dandekar T, Bork P:

**Variation and evolution of the citric-acid cycle: a genomic perspective.***Trends Microbiol*1999,**7:**281–291. 10.1016/S0966-842X(99)01539-5 - 32.
- 33.

## Acknowledgements

This work was supported in part by the Austrian *Fonds zur Förderung der Wissenschaftlichen Forschung*, Project No. P15893, by the German *DFG* Bioinformatics Initiative BIZ-6/1-2, by the Austrian *Gen-AU bioinformatics integration network* sponsored by BM-BWK and BM-WA, and by the Laboratory Directed Research and Development program of the Los Alamos National Laboratory, Project No. 20040184ER.

## Author information

### Affiliations

### Corresponding author

## Additional information

### Authors' contributions

The idea on a metabolic network algebra has been conceived by ILH and PFS during a visit at the Santa Fe Institute. The implementation was substantially improved by CF. CVF developed the distance measurements and applications to microbial metabolic networks. All authors collaborated closely in writing the manuscript.

## Electronic supplementary material

### 12859_2005_806_MOESM1_ESM.gz

Additional file 1: The Vienna Reaction Network Library Vienna-RNL implements these basic set-theoretic operations on chemical reaction networks. It is available under the GNU Public License from [32] (GZ 232 KB)

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

**Open Access**
This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License (
https://creativecommons.org/licenses/by/2.0
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Forst, C.V., Flamm, C., Hofacker, I.L. *et al.* Algebraic comparison of metabolic networks, phylogenetic inference, and metabolic innovation.
*BMC Bioinformatics* **7, **67 (2006). https://doi.org/10.1186/1471-2105-7-67

Received:

Accepted:

Published:

### Keywords

- Metabolic Network
- Reaction Network
- Symmetric Difference
- Tree Distance
- Functional Replacement