BACKGROUND
recent studies indicate that metabolic networks evolve at the local level through patchwork evolution and retrograde evolution  <cit> .
patchwork evolution, which is likely to be more important, occurs when an enzyme evolves from a broad spectrum enzyme to an enzyme with a highly specialized activity  <cit> .
retrograde evolution is a process where the depletion of a substrate from the environment leads to the evolution of an enzyme which can accept a new substrate and catalyze the production of the depleted substance  <cit> .
metabolic networks and other complex networks such as the film actor collaboration network, the world wide web, protein domain networks and protein-protein interaction networks are small-world networks with some properties which are consistent with scale-free networks  <cit> .
the small-worldness of the metabolic network of e. coli has recently been contested for an alternative network representation where carbon atomic traces in metabolic reactions were used  <cit> .
a small-world network is characterized by 1) short path lengths between any two nodes in the network and 2) a high clustering coefficient, which means that the neighbors of a certain node of the network are often connected to each other thereby forming clusters.
a scale-free network, in this context, has a power-law connectivity  distribution, i.e.
there are many nodes which have very low connectivities and a handful of nodes with much higher connectivities , see figure 1a.
scale-free networks are robust networks in the sense that they often remain intact when a large fraction of randomly chosen nodes is eliminated from the network  <cit> .
however, if a small fraction of the hubs of the network is eliminated the network is likely to become fragmented into several components.
it has been suggested that the scale-free character of biological networks has evolved through natural selection for the advantage of robustness and error-tolerance that the scale-free network topology confers to the organism  <cit> .
a study by gleiss et al showed that chemical reaction networks, which clearly have not been subjected to natural selection, also show scale-free characteristics thereby showing that scale-free networks can arise without natural selection and may be a general feature of chemical reaction networks  <cit> .
networks with scale-free properties have been shown to evolve when two simple rules are applied: 1) the network grows by the addition of new nodes.
2) preferential attachment: new nodes are more likely to become connected to well connected nodes in the network  <cit> .
while preferential attachment is often at the root of scale-freeness, a network with an power-law degree distribution might be produced through other mechanisms.
preferential attachment in the context of genetic networks may take place partly through gene duplication  <cit> .
in agreement with preferential attachment eisenberg and levanon  <cit>  showed that the proteins which have homologs in all  <dig> domains of life, which are likely to be of ancient origin, have higher connectivities in the protein-protein interaction network of s. cerevisiae.
in contrast, kunin et al  <cit>  recently showed that the most highly connected proteins date to after the evolution of primordial eukaryotes but before the radiation of eukaryotes to plants, metazoa and protista.
here, we investigate the evidence for preferential attachment and the role of horizontal gene transfer in the metabolic network evolution of e. coli.
RESULTS
connectivity and phylogenetic group
if preferential attachment is an important mechanism in the evolution of metabolic networks older enzymes should have a higher average connectivity  than younger enzymes.
in order to investigate this prediction we extracted the enzymes and the reactions in e. coli from the ecocyc  <cit>  and kegg  <cit>  databases.
the network representation of the metabolic network of e. coli was constructed using ecocyc, see methods.
the nodes in our graph represent the enzymes  catalyzing the reactions and the edges represent one or more compounds involved in the reactions.
there is an edge from enzyme el to enzyme e <dig> if el catalyzes a reaction where compound a is produced and then e <dig> uses a as substrate.
there can be at most one edge in each direction between the nodes in the graph.
the connectivity of a node is defined as the number of edges connecting the node to other nodes in the network.
enzymes  were collected from kegg orthology  <cit>  and were found in  <dig> different organisms .
among these the e. coli enzymes with representatives in  <dig> organisms  of reasonably well understood phylogenies were extracted for further studies, see figure  <dig>  the enzymes were divided into five age groups, see table  <dig>  the enzymes in group  <dig> are likely to be among the oldest since they have representatives in eukaryotes, archaea and bacteria.
group  <dig> contains the enzymes with representatives in eukaryotes and bacteria but not in archaea while group  <dig> contains the enzymes with representatives in archaea and bacteria but not in eukaryotes.
group  <dig> contains e. coli enzymes without representatives in eukaryotes or archaea but with representatives in bacteria other than γβ-proteobacteria and group  <dig> contains enzymes with representatives in γβ-proteobacteria only.
enzymes belonging to group 1– <dig> are probably ancient enzymes since they exist in at least two domains of life.
enzymes which are found in group  <dig> are only found in bacteria which are comparatively close relatives of e. coli, in the γβ-proteobacteria group, see figure  <dig>  which indicates that they are relatively recent additions to the metabolic repertory of e. coli.
it is possible that some of the enzymes in groups 1– <dig> could have evolved relatively recently and subsequently been horizontally transferred to the other domains of life but the evidence of horizontal gene transfer between organisms belonging to different domains of life is not abundant.
horizontal gene transfer between bacterial species is believed to be more common  <cit>  and therefore there may be enzymes in group  <dig> which have been transferred from γβ-proteobacteria.
gene loss, which is estimated to be three times more common than horizontal gene transfer in prokaryotes  <cit> , is probably a more important source of error in this study since genes which serve an important function in bacteria but not in archaea or eukaryotes may have been lost in the archaeal and eukaryotic lineages and as a consequence groups  <dig> and  <dig> may not exclusively contain relatively recently evolved enzymes.
the average connectivities for the e. coli enzymes with representatives in the five groups were calculated, see figure 3a.
the enzymes which are represented in all  <dig> domains of life  have the highest average connectivity together with the enzymes which occur in bacteria and eukaryotes  while the enzymes that occur only in γβ-proteobacteria  have the lowest connectivities.
the average connectivity was 40% higher for group  <dig> than for group  <dig> enzymes.
in order to estimate the significance of the results  <dig>  <dig> randomized networks were generated through shuffling the group numbers while preserving the network topology.
the z-score, see methods, of the enzymes that occur in eukaryotes is substantially higher than the z-score for the enzymes which only occur in archaea and bacteria, see figure 3b.
in a similar manner  <dig> other organisms were investigated.
we collected the metabolic reactions for these organisms from kegg.
it should be noted that these metabolic networks have not been investigated to the extent of the e. coli metabolic network.
therefore, the connectivities of the enzymes and the results are probably not of the same reliability.
we found that the larger prokaryotic genomes often have a particularly strong correlation between connectivity and domain presence while the correlation in smaller prokaryotes and eukaryotes is weaker, see figure  <dig>  many of the smaller prokaryotes are obligate intracellular parasites or symbionts.
the genomes of obligate parasites and symbionts have been metabolically reduced and many metabolic functions, such as the amino acid metabolism, are frequently provided by the host.
since enzymes involved in amino acid metabolism are often ancient proteins with high connectivities, see figure 5a and 5b, the absence of genes coding for these enzymes in the genomes of the obligate symbionts/parasites can account for the lack of correlation between connectivity and domain presence in smaller genomes.
in conclusion we found that e. coli enzymes which have representatives in all domains of life, and in eukaryotes but not archaea, have a higher average connectivity in the metabolic network of e. coli than the presumably younger enzymes which only have representatives in γβ-proteobacteria.
this finding lends support for one of the predictions of the mechanism of preferential attachment.
connectivity and horizontal gene transfer
it has been suggested that the scale-free properties of biological networks may arise, at least partially, as a result of preferential attachment of new nodes to highly connected nodes through gene duplication  <cit> .
preferential attachment by gene duplication may take place according to the following scenario; initially, the duplicated gene has exactly the same function and position in the network as the template gene.
since many genes are connected to the hub of the network, the duplicated gene is by chance likely to be connected to the hub of the network.
subsequently, the duplicate gene may evolve towards another functionality but it could retain some of its original function.
for instance, a multi-domain protein could loose one of its domains through deletion but retain the other domains and possibly part of its original functionality.
in such a scenario the older proteins are more likely to be highly connected than the younger proteins.
an alternative scenario is preferential attachment by horizontal gene transfer ; a new, or alternative, enzyme is introduced through hgt.
the new enzyme is more likely to be retained in the metabolic repertory if it confers a new or improved function at a central, rather than peripheral, position of the metabolism – such as if it is connected to a highly connected enzyme or if it is itself highly connected.
this is a consideration which may be particularly important in the metabolism of bacteria since some bacteria are prone to delete dispensable genes from their genomes  <cit> .
arguably, connectivity is a measure which indicates the centrality and importance of an enzyme in which case horizontally transferred genes should frequently be highly connected or be connected to highly connected enzymes.
according to this scenario, horizontally transferred enzymes would be preferentially attached to highly connected enzymes and/or be preferentially replacing highly connected enzymes.
although the extent of the evolutionary impact of hgt is still under debate  <cit> , it is generally accepted as an important evolutionary process in microbial species  <cit> .
roughly 18% of the protein coding genes in e. coli are likely to have been introduced into the e. coli genome by hgt since the species diverged from the salmonella lineage according to an analysis by lawrence and ochman where base composition and codon usage patterns were used to identify the horizontally transferred sequences  <cit> .
using this data set we found that  <dig> of the  <dig> e. coli enzymes used in this study are likely to be examples of hgt .
the average connectivity for the hgt enzymes in the metabolic network of e. coli is  <dig>  while the other enzymes  have an average connectivity of  <dig> .
the most striking difference between the connectivity distributions of the hgt enzymes and the non-hgt enzymes is that only 15% of the hgt enzymes have connectivities between 0– <dig> and 37% of the hgt enzymes have higher connectivities than  <dig> while 28% of the non-hgt enzymes have connectivities of 0– <dig> and 26% have connectivities higher than  <dig>  see figure  <dig>  these results indicate that horizontally transferred enzymes in the metabolic network of e. coli are either introduced into the repertory of the organism as a comparatively high connectivity enzyme or acquires many connections during the evolution of the network.
we classified the hgt enzymes into the five phylogenetic groups and determined the average connectivities for each group, see figure  <dig> and table  <dig>  only hgt enzymes belonging to group  <dig> and  <dig> have higher connectivities than the non-hgt enzymes belonging to the same group.
from figure 7a it is clear that the high average connectivities of group  <dig> and  <dig> enzymes, which was seen in figure 3a, is partly but not solely due to hgt enzymes.
it should be noted that the horizontal transfers of enzymes which are specific for βγ-proteobacteria only  may be underestimated since hgt events between closely related species are hard to detect  <cit> .
most horizontally transferred genes go through the process of amelioration, the adjustment of the transferred sequence to the base composition and codon usage of the resident genome.
therefore, most detectable hgts have taken place relatively recently in the history of e. coli  <cit> .
consequently, we can conclude that while it is true that the highly connected enzymes in the metabolic network of e. coli are often old in the sense that they are enzymes with representatives in eukaryotes, and which therefore probably originated in the last common ancestor of eukaryotes and bacteria, they are also overrepresented among the enzymes which have been introduced recently into the e. coli genome through hgt.
these findings suggest that horizontally transferred genes are introduced and retained preferentially at central positions of the metabolism of e. coli.
the connectivity of essential enzymes and isozymes
jeong et al  <cit>  showed that the highly connected proteins in the protein-protein interaction network of s. cerevisiae are more likely to be indispensable to the organism than less well connected proteins.
we wished to study if there was a similar correlation in the metabolic network of e. coli.
we used the essentiality classification from the study of gerdes et al  <cit>  of e. coli under aerobic growth in nutrition rich medium.
we calculated the mean connectivity for the essential and the dispensable enzymes respectively and found that the essential enzymes do not show a higher connectivity than expected .
it is possible that the relatively small size of the metabolic networks compared to the protein-protein interaction network is the reason a similar correlation could not be found in the metabolic network of e. coli.
the hubs are the most important nodes for the integrity of the network.
if a fraction of the hubs are removed the network is likely to become fragmented into smaller components.
since these enzymes are very important for the robustness of the network it might be suspected that the ec numbers with the highest connectivities could have more than one representatives in the genome, i.e.
that there are two or more isozymes representing these highly connected nodes.
isozymes in multicellular organisms are often active in different tissues while isozymes in single cellular organisms frequently have different substrate specificities or are activated in different environments .
we here designated a pair of enzymes as isozymes if they catalyze the same reaction but are coded for by different genes, which are not part of the same enzyme complex.
we used expasy  <cit> , sgd  <cit>  and ecocyc  <cit>  to determine which enzymes in the metabolic networks of e. coli and s. cerevisiae occur as isozymes.
we found  <dig> ec numbers that were associated with isozymes in e. coli and  <dig> ec numbers that were associated with isozymes in s. cerevisiae, see additional files.
the mean connectivities for the isozymes and the non-isozymes were determined and the result was compared to randomized networks.
we found that the isozymes do not have a noticeably higher mean connectivity than non-isozymes .
the result may indicate that isozymes are not necessarily crucial for the integrity of the metabolic network.
in accordance with our result it has recently been shown that the isozymes of s. cerevisiae are not overrepresented among essential enzymes  <cit> .
connectivity and function
kunin et al showed that the functional classes in the protein-protein network of s. cere-visiae display distinctly different connectivity levels  <cit> .
in a similar manner we investigated whether enzymes belonging to different functional groups are characterized by distinct connectivities.
we classified the enzymes into  <dig> functional classes according to ecocyc  <cit> ; lipid metabolism, nucleotide metabolism, amino acid metabolism, sugar metabolism, energy metabolism, other biosynthesis and other degradation and calculated the mean connectivities for the different functional classes, see figure 5b and table  <dig>  the mean connectivities for the enzymes involved in nucleotide, amino acid, other degradation and energy metabolism are higher than expected.
the amino acid metabolism and nucleotide metabolism enzymes are clearly over represented in  <dig> domains of life while enzymes involved in energy metabolism are slightly more common in  <dig> domains of life than expected by chance.
many of the pathways involved in energy metabolism, such as the citric acid cycle and glycolysis, are believed to be very old.
however, there are substantial variations in the energy metabolism between different species and domains of life  <cit> .
therefore, the observation that energy metabolism enzymes are not overrepresented in three domains of life is not surprising.
contrastingly, enzymes involved in lipid and sugar metabolism are on average half as well connected as the enzymes involved in nucleotide, amino acid and energy metabolism.
the group of enzymes involved in lipid metabolism is less than half the size of the second smallest functional group and due to its small size the z-score for this functional group is less reliable than for the other functional groups.
the sugar metabolism enzymes are clearly over represented among the enzymes that occur in bacteria only, see figure 5a, which was anticipated since there are many bacterial specific enzymes involved in sugar transportation  <cit> .
network growth through preferential attachment
according to the mechanism of preferential attachment new enzymes in the network should be preferentially attached to already well connected nodes.
we do not have access to the last universal common ancestor  that existed before the  <dig> domains of life evolved.
however, a rough representation of the metabolic network of that organism was created by extracting the enzymes that occur in all domains of life.
the connectivities of the enzymes in that derived network were determined.
we then calculated the number of enzyme nodes that have been added to an enzyme by subtracting the connectivity of the enzyme in the current e. coli network by the connectivity of the enzyme in the ancient network.
we found that the enzymes which have higher connectivities in the ancient network gain new connections at a higher rate than the enzymes with lower connectivities, see figure  <dig>  the correlation between connectivity in the ancient network and the connectivity increase appears to be linear  following the equation f =  <dig>  +  <dig> x, where x is the connectivity in the ancient network.
we can therefore conclude that the addition of new nodes to the metabolic network of e. coli occurs in a manner which is consistent with preferential attachment.
discussion
we have investigated two predictions generated from the mechanism of preferential attachment in the evolution of the metabolic network of e. coli.
first, if preferential attachment is of any significance in the evolution of the metabolic network of e. coli, the older enzymes in the network should have a higher average connectivity.
we have found that e. coli enzymes which are represented in three domains of life, and in eukaryotes but not archaea, have a higher average connectivity than expected by chance.
second, another prediction generated from the hypothesis of network evolution through preferential attachment is that highly connected nodes should gain new edges at a faster rate than nodes with low connectivities.
to investigate this prediction we extracted the enzymes with representatives in  <dig> domains of life and determined the network representing luca's metabolic network.
in accordance with the mechanism of preferential attachment we found a positive linear correlation between connectivity in the ancient network and number of connections gained through evolution.
further, we found that the e. coli enzymes which are believed to have undergone horizontal gene transfer  have a higher average connectivity than other enzymes .
this is especially true for the hgt enzymes with representatives in eukaryotes, which is the most highly connected group of e. coli enzymes.
this result suggests that the highly connected enzymes are often old in the sense that they are likely to have originated in luca and been part of the bacterial metabolic repertory for a long time.
however, these ancient enzymes are sometimes relatively recent additions to the metabolic network of e. coli.
it is possible that bacteria such as e. coli are adjusting their metabolic networks to a changing environment by replacing the relatively central enzymes, with high connectivities, for better adapted orthologs from other prokaryotic species.
CONCLUSIONS
it is well known that many novel functions in organisms are obtained through gene duplication, followed by subfunctionalization and neofunctionalization.
therefore, a possible biological explanation for the preferential attachment growth of metabolic networks, which we have now found some support for, could be that novel enzymes, which are created through gene duplication, maintain some compounds involved in the reaction catalyzed by the original enzyme throughout its future evolution.
as a supplementary explanation we propose that horizontally transferred enzymes are introduced preferentially at central positions of the metabolic network of e. coli.

 recent studies indicate that metabolic networks evolve at the local level through patchwork evolution and retrograde evolution. patchwork evolution, which is likely to be more important, occurs when an enzyme evolves from a broad spectrum enzyme to an enzyme with a highly specialized activity. retrograde evolution is a process where the depletion of a substrate from the environment leads to the evolution of an enzyme which can accept a new substrate and catalyze the production of the depleted substance. metabolic networks and other complex networks such as the film actor collaboration network, the world wide web, protein domain networks and protein-protein interaction networks are small-world networks with some properties which are consistent with scale-free networks. the small-worldness of the metabolic network of e. coli has recently been contested for an alternative network representation where carbon atomic traces in metabolic reactions were used. a small-world network is characterized by 1) short path lengths between any two nodes in the network and 2) a high clustering coefficient, which means that the neighbors of a certain node of the network are often connected to each other thereby forming clusters. a scale-free network, in this context, has a power-law connectivity distribution, i.e. there are many nodes which have very low connectivities and a handful of nodes with much higher connectivities, see figure 1a. scale-free networks are robust networks in the sense that they often remain intact when a large fraction of randomly chosen nodes is eliminated from the network. however, if a small fraction of the hubs of the network is eliminated the network is likely to become fragmented into several components. it has been suggested that the scale-free character of biological networks has evolved through natural selection for the advantage of robustness and error-tolerance that the scale-free network topology confers to the organism. a study by gleiss et al showed that chemical reaction networks, which clearly have not been subjected to natural selection, also show scale-free characteristics thereby showing that scale-free networks can arise without natural selection and may be a general feature of chemical reaction networks. networks with scale-free properties have been shown to evolve when two simple rules are applied: 1) the network grows by the addition of new nodes. 2) preferential attachment: new nodes are more likely to become connected to well connected nodes in the network. while preferential attachment is often at the root of scale-freeness, a network with an power-law degree distribution might be produced through other mechanisms. preferential attachment in the context of genetic networks may take place partly through gene duplication. in agreement with preferential attachment eisenberg and levanon showed that the proteins which have homologs in all domains of life, which are likely to be of ancient origin, have higher connectivities in the protein-protein interaction network of s. cerevisiae. in contrast, kunin et al recently showed that the most highly connected proteins date to after the evolution of primordial eukaryotes but before the radiation of eukaryotes to plants, metazoa and protista. here, we investigate the evidence for preferential attachment and the role of horizontal gene transfer in the metabolic network evolution of e. coli.  connectivity and phylogenetic group if preferential attachment is an important mechanism in the evolution of metabolic networks older enzymes should have a higher average connectivity than younger enzymes. in order to investigate this prediction we extracted the enzymes and the reactions in e. coli from the ecocyc and kegg databases. the network representation of the metabolic network of e. coli was constructed using ecocyc, see methods. the nodes in our graph represent the enzymes catalyzing the reactions and the edges represent one or more compounds involved in the reactions. there is an edge from enzyme el to enzyme e if el catalyzes a reaction where compound a is produced and then e uses a as substrate. there can be at most one edge in each direction between the nodes in the graph. the connectivity of a node is defined as the number of edges connecting the node to other nodes in the network. enzymes were collected from kegg orthology and were found in different organisms. among these the e. coli enzymes with representatives in organisms of reasonably well understood phylogenies were extracted for further studies, see figure the enzymes were divided into five age groups, see table the enzymes in group are likely to be among the oldest since they have representatives in eukaryotes, archaea and bacteria. group contains the enzymes with representatives in eukaryotes and bacteria but not in archaea while group contains the enzymes with representatives in archaea and bacteria but not in eukaryotes. group contains e. coli enzymes without representatives in eukaryotes or archaea but with representatives in bacteria other than γβ-proteobacteria and group contains enzymes with representatives in γβ-proteobacteria only. enzymes belonging to group 1– are probably ancient enzymes since they exist in at least two domains of life. enzymes which are found in group are only found in bacteria which are comparatively close relatives of e. coli, in the γβ-proteobacteria group, see figure which indicates that they are relatively recent additions to the metabolic repertory of e. coli. it is possible that some of the enzymes in groups 1– could have evolved relatively recently and subsequently been horizontally transferred to the other domains of life but the evidence of horizontal gene transfer between organisms belonging to different domains of life is not abundant. horizontal gene transfer between bacterial species is believed to be more common and therefore there may be enzymes in group which have been transferred from γβ-proteobacteria. gene loss, which is estimated to be three times more common than horizontal gene transfer in prokaryotes, is probably a more important source of error in this study since genes which serve an important function in bacteria but not in archaea or eukaryotes may have been lost in the archaeal and eukaryotic lineages and as a consequence groups and may not exclusively contain relatively recently evolved enzymes. the average connectivities for the e. coli enzymes with representatives in the five groups were calculated, see figure 3a. the enzymes which are represented in all domains of life have the highest average connectivity together with the enzymes which occur in bacteria and eukaryotes while the enzymes that occur only in γβ-proteobacteria have the lowest connectivities. the average connectivity was 40% higher for group than for group enzymes. in order to estimate the significance of the results randomized networks were generated through shuffling the group numbers while preserving the network topology. the z-score, see methods, of the enzymes that occur in eukaryotes is substantially higher than the z-score for the enzymes which only occur in archaea and bacteria, see figure 3b. in a similar manner other organisms were investigated. we collected the metabolic reactions for these organisms from kegg. it should be noted that these metabolic networks have not been investigated to the extent of the e. coli metabolic network. therefore, the connectivities of the enzymes and the results are probably not of the same reliability. we found that the larger prokaryotic genomes often have a particularly strong correlation between connectivity and domain presence while the correlation in smaller prokaryotes and eukaryotes is weaker, see figure many of the smaller prokaryotes are obligate intracellular parasites or symbionts. the genomes of obligate parasites and symbionts have been metabolically reduced and many metabolic functions, such as the amino acid metabolism, are frequently provided by the host. since enzymes involved in amino acid metabolism are often ancient proteins with high connectivities, see figure 5a and 5b, the absence of genes coding for these enzymes in the genomes of the obligate symbionts/parasites can account for the lack of correlation between connectivity and domain presence in smaller genomes. in conclusion we found that e. coli enzymes which have representatives in all domains of life, and in eukaryotes but not archaea, have a higher average connectivity in the metabolic network of e. coli than the presumably younger enzymes which only have representatives in γβ-proteobacteria. this finding lends support for one of the predictions of the mechanism of preferential attachment. connectivity and horizontal gene transfer it has been suggested that the scale-free properties of biological networks may arise, at least partially, as a result of preferential attachment of new nodes to highly connected nodes through gene duplication. preferential attachment by gene duplication may take place according to the following scenario; initially, the duplicated gene has exactly the same function and position in the network as the template gene. since many genes are connected to the hub of the network, the duplicated gene is by chance likely to be connected to the hub of the network. subsequently, the duplicate gene may evolve towards another functionality but it could retain some of its original function. for instance, a multi-domain protein could loose one of its domains through deletion but retain the other domains and possibly part of its original functionality. in such a scenario the older proteins are more likely to be highly connected than the younger proteins. an alternative scenario is preferential attachment by horizontal gene transfer; a new, or alternative, enzyme is introduced through hgt. the new enzyme is more likely to be retained in the metabolic repertory if it confers a new or improved function at a central, rather than peripheral, position of the metabolism – such as if it is connected to a highly connected enzyme or if it is itself highly connected. this is a consideration which may be particularly important in the metabolism of bacteria since some bacteria are prone to delete dispensable genes from their genomes. arguably, connectivity is a measure which indicates the centrality and importance of an enzyme in which case horizontally transferred genes should frequently be highly connected or be connected to highly connected enzymes. according to this scenario, horizontally transferred enzymes would be preferentially attached to highly connected enzymes and/or be preferentially replacing highly connected enzymes. although the extent of the evolutionary impact of hgt is still under debate, it is generally accepted as an important evolutionary process in microbial species. roughly 18% of the protein coding genes in e. coli are likely to have been introduced into the e. coli genome by hgt since the species diverged from the salmonella lineage according to an analysis by lawrence and ochman where base composition and codon usage patterns were used to identify the horizontally transferred sequences. using this data set we found that of the e. coli enzymes used in this study are likely to be examples of hgt. the average connectivity for the hgt enzymes in the metabolic network of e. coli is while the other enzymes have an average connectivity of. the most striking difference between the connectivity distributions of the hgt enzymes and the non-hgt enzymes is that only 15% of the hgt enzymes have connectivities between 0– and 37% of the hgt enzymes have higher connectivities than while 28% of the non-hgt enzymes have connectivities of 0– and 26% have connectivities higher than see figure these results indicate that horizontally transferred enzymes in the metabolic network of e. coli are either introduced into the repertory of the organism as a comparatively high connectivity enzyme or acquires many connections during the evolution of the network. we classified the hgt enzymes into the five phylogenetic groups and determined the average connectivities for each group, see figure and table only hgt enzymes belonging to group and have higher connectivities than the non-hgt enzymes belonging to the same group. from figure 7a it is clear that the high average connectivities of group and enzymes, which was seen in figure 3a, is partly but not solely due to hgt enzymes. it should be noted that the horizontal transfers of enzymes which are specific for βγ-proteobacteria only may be underestimated since hgt events between closely related species are hard to detect. most horizontally transferred genes go through the process of amelioration, the adjustment of the transferred sequence to the base composition and codon usage of the resident genome. therefore, most detectable hgts have taken place relatively recently in the history of e. coli. consequently, we can conclude that while it is true that the highly connected enzymes in the metabolic network of e. coli are often old in the sense that they are enzymes with representatives in eukaryotes, and which therefore probably originated in the last common ancestor of eukaryotes and bacteria, they are also overrepresented among the enzymes which have been introduced recently into the e. coli genome through hgt. these findings suggest that horizontally transferred genes are introduced and retained preferentially at central positions of the metabolism of e. coli. the connectivity of essential enzymes and isozymes jeong et al showed that the highly connected proteins in the protein-protein interaction network of s. cerevisiae are more likely to be indispensable to the organism than less well connected proteins. we wished to study if there was a similar correlation in the metabolic network of e. coli. we used the essentiality classification from the study of gerdes et al of e. coli under aerobic growth in nutrition rich medium. we calculated the mean connectivity for the essential and the dispensable enzymes respectively and found that the essential enzymes do not show a higher connectivity than expected. it is possible that the relatively small size of the metabolic networks compared to the protein-protein interaction network is the reason a similar correlation could not be found in the metabolic network of e. coli. the hubs are the most important nodes for the integrity of the network. if a fraction of the hubs are removed the network is likely to become fragmented into smaller components. since these enzymes are very important for the robustness of the network it might be suspected that the ec numbers with the highest connectivities could have more than one representatives in the genome, i.e. that there are two or more isozymes representing these highly connected nodes. isozymes in multicellular organisms are often active in different tissues while isozymes in single cellular organisms frequently have different substrate specificities or are activated in different environments. we here designated a pair of enzymes as isozymes if they catalyze the same reaction but are coded for by different genes, which are not part of the same enzyme complex. we used expasy, sgd and ecocyc to determine which enzymes in the metabolic networks of e. coli and s. cerevisiae occur as isozymes. we found ec numbers that were associated with isozymes in e. coli and ec numbers that were associated with isozymes in s. cerevisiae, see additional files. the mean connectivities for the isozymes and the non-isozymes were determined and the result was compared to randomized networks. we found that the isozymes do not have a noticeably higher mean connectivity than non-isozymes. the result may indicate that isozymes are not necessarily crucial for the integrity of the metabolic network. in accordance with our result it has recently been shown that the isozymes of s. cerevisiae are not overrepresented among essential enzymes. connectivity and function kunin et al showed that the functional classes in the protein-protein network of s. cere-visiae display distinctly different connectivity levels. in a similar manner we investigated whether enzymes belonging to different functional groups are characterized by distinct connectivities. we classified the enzymes into functional classes according to ecocyc; lipid metabolism, nucleotide metabolism, amino acid metabolism, sugar metabolism, energy metabolism, other biosynthesis and other degradation and calculated the mean connectivities for the different functional classes, see figure 5b and table the mean connectivities for the enzymes involved in nucleotide, amino acid, other degradation and energy metabolism are higher than expected. the amino acid metabolism and nucleotide metabolism enzymes are clearly over represented in domains of life while enzymes involved in energy metabolism are slightly more common in domains of life than expected by chance. many of the pathways involved in energy metabolism, such as the citric acid cycle and glycolysis, are believed to be very old. however, there are substantial variations in the energy metabolism between different species and domains of life. therefore, the observation that energy metabolism enzymes are not overrepresented in three domains of life is not surprising. contrastingly, enzymes involved in lipid and sugar metabolism are on average half as well connected as the enzymes involved in nucleotide, amino acid and energy metabolism. the group of enzymes involved in lipid metabolism is less than half the size of the second smallest functional group and due to its small size the z-score for this functional group is less reliable than for the other functional groups. the sugar metabolism enzymes are clearly over represented among the enzymes that occur in bacteria only, see figure 5a, which was anticipated since there are many bacterial specific enzymes involved in sugar transportation. network growth through preferential attachment according to the mechanism of preferential attachment new enzymes in the network should be preferentially attached to already well connected nodes. we do not have access to the last universal common ancestor that existed before the domains of life evolved. however, a rough representation of the metabolic network of that organism was created by extracting the enzymes that occur in all domains of life. the connectivities of the enzymes in that derived network were determined. we then calculated the number of enzyme nodes that have been added to an enzyme by subtracting the connectivity of the enzyme in the current e. coli network by the connectivity of the enzyme in the ancient network. we found that the enzymes which have higher connectivities in the ancient network gain new connections at a higher rate than the enzymes with lower connectivities, see figure the correlation between connectivity in the ancient network and the connectivity increase appears to be linear following the equation f = + x, where x is the connectivity in the ancient network. we can therefore conclude that the addition of new nodes to the metabolic network of e. coli occurs in a manner which is consistent with preferential attachment. discussion we have investigated two predictions generated from the mechanism of preferential attachment in the evolution of the metabolic network of e. coli. first, if preferential attachment is of any significance in the evolution of the metabolic network of e. coli, the older enzymes in the network should have a higher average connectivity. we have found that e. coli enzymes which are represented in three domains of life, and in eukaryotes but not archaea, have a higher average connectivity than expected by chance. second, another prediction generated from the hypothesis of network evolution through preferential attachment is that highly connected nodes should gain new edges at a faster rate than nodes with low connectivities. to investigate this prediction we extracted the enzymes with representatives in domains of life and determined the network representing luca's metabolic network. in accordance with the mechanism of preferential attachment we found a positive linear correlation between connectivity in the ancient network and number of connections gained through evolution. further, we found that the e. coli enzymes which are believed to have undergone horizontal gene transfer have a higher average connectivity than other enzymes. this is especially true for the hgt enzymes with representatives in eukaryotes, which is the most highly connected group of e. coli enzymes. this result suggests that the highly connected enzymes are often old in the sense that they are likely to have originated in luca and been part of the bacterial metabolic repertory for a long time. however, these ancient enzymes are sometimes relatively recent additions to the metabolic network of e. coli. it is possible that bacteria such as e. coli are adjusting their metabolic networks to a changing environment by replacing the relatively central enzymes, with high connectivities, for better adapted orthologs from other prokaryotic species.  it is well known that many novel functions in organisms are obtained through gene duplication, followed by subfunctionalization and neofunctionalization. therefore, a possible biological explanation for the preferential attachment growth of metabolic networks, which we have now found some support for, could be that novel enzymes, which are created through gene duplication, maintain some compounds involved in the reaction catalyzed by the original enzyme throughout its future evolution. as a supplementary explanation we propose that horizontally transferred enzymes are introduced preferentially at central positions of the metabolic network of e. coli. 
interestingly, the enzymes which have been proposed as candidates for horizontal gene transfer have a higher average connectivity than the other enzymes. in addition, we found that enzymes which are candidates for horizontal gene transfer have a higher average connectivity than other enzymes. this indicates that while new enzymes are attached preferentially to highly connected enzymes, these highly connected enzymes have sometimes been introduced into the e. coli genome by horizontal gene transfer. furthermore, it was found that new edges are added to the highly connected enzymes at a faster rate than to enzymes with low connectivities which is consistent with preferential attachment. a possible biological explanation for preferential attachment growth of metabolic networks is that novel enzymes created through gene duplication maintain some of the compounds involved in the original reaction, throughout its future evolution. e. coli enzymes which have representatives in eukaryotes have a higher average connectivity while enzymes which are represented only in the prokaryotes, and especially the enzymes only present in βγ-proteobacteria, have lower connectivities than expected by chance. here, we have found indications of preferential attachment in the metabolic network of e. coli. we speculate that e. coli has adjusted its metabolic network to a changing environment by replacing the relatively central enzymes for better adapted orthologs from other prokaryotic species. many biological networks show some characteristics of scale-free networks. scale-free networks can evolve through preferential attachment where new nodes are preferentially attached to well connected nodes. here we have investigated preferential attachment in the context of metabolic networks. in networks which have evolved through preferential attachment older nodes should have a higher average connectivity than younger nodes. the connectivities of the enzymes in the metabolic network of escherichia coli were determined and representatives for these enzymes were located in eukaryotes, archaea and bacteria.
 many biological networks show some characteristics of scale-free networks. scale-free networks can evolve through preferential attachment where new nodes are preferentially attached to well connected nodes. in networks which have evolved through preferential attachment older nodes should have a higher average connectivity than younger nodes. here we have investigated preferential attachment in the context of metabolic networks.  the connectivities of the enzymes in the metabolic network of escherichia coli were determined and representatives for these enzymes were located in eukaryotes, archaea and bacteria. e. coli enzymes which have representatives in eukaryotes have a higher average connectivity while enzymes which are represented only in the prokaryotes, and especially the enzymes only present in βγ-proteobacteria, have lower connectivities than expected by chance. interestingly, the enzymes which have been proposed as candidates for horizontal gene transfer have a higher average connectivity than the other enzymes. furthermore, it was found that new edges are added to the highly connected enzymes at a faster rate than to enzymes with low connectivities which is consistent with preferential attachment.  here, we have found indications of preferential attachment in the metabolic network of e. coli. a possible biological explanation for preferential attachment growth of metabolic networks is that novel enzymes created through gene duplication maintain some of the compounds involved in the original reaction, throughout its future evolution. in addition, we found that enzymes which are candidates for horizontal gene transfer have a higher average connectivity than other enzymes. this indicates that while new enzymes are attached preferentially to highly connected enzymes, these highly connected enzymes have sometimes been introduced into the e. coli genome by horizontal gene transfer. we speculate that e. coli has adjusted its metabolic network to a changing environment by replacing the relatively central enzymes for better adapted orthologs from other prokaryotic species. 