- Research article
- Open Access
Environmental variability and modularity of bacterial metabolic networks
BMC Evolutionary Biology volume 7, Article number: 169 (2007)
Biological systems are often modular: they can be decomposed into nearly-independent structural units that perform specific functions. The evolutionary origin of modularity is a subject of much current interest. Recent theory suggests that modularity can be enhanced when the environment changes over time. However, this theory has not yet been tested using biological data.
To address this, we studied the relation between environmental variability and modularity in a natural and well-studied system, the metabolic networks of bacteria. We classified 117 bacterial species according to the degree of variability in their natural habitat. We find that metabolic networks of organisms in variable environments are significantly more modular than networks of organisms that evolved under more constant conditions.
This study supports the view that variability in the natural habitat of an organism promotes modularity in its metabolic network and perhaps in other biological systems.
Biological systems often show modularity, in the sense that they can be separated into nearly-independent structural subsystems, each of which performs a specific function [1–12]. The origin and preservation of modularity in biology is a subject of current interest [13–18].
One approach to study the origin of modularity employs computer simulations of biological evolution . Such simulations make random changes to a set of networks, and select those that best satisfy a given goal. Simulations towards a constant goal usually lead to non-modular networks. Even an initially modular network rapidly evolves connections that reduce modularity but increase fitness . This raises the question of which evolutionary forces can generate and maintain modular structure.
One possible force that can lead to modularity was suggested by evolutionary simulations where the goal changes over time [16, 18]. Modular networks spontaneously evolved when the goal changed with time in a way that preserves the same sub-goals but in different combinations . Under such 'modularly varying goals', the networks evolved and maintained modular structure, with a module for each sub-goal. When the goal changed, the connections between these modules were rapidly rewired to adapt to the new goal.
Here, we attempt to test these findings in a natural biological system, by asking whether there is a correlation between metabolic network's modularity and the variability of the environment in which it evolved. We concentrated on the metabolic networks of bacteria due to the availability and quality of the data. These networks can be systematically compared between species that live in environments that differ in their degree of variability.
Metabolic networks represent the set of biochemical metabolic reactions within a living cell . Metabolic networks of diverse bacterial species have been reconstructed based on their genomic sequence and additional biochemical data [21–23]. Previous structural comparisons of metabolic networks focused on differences between kingdoms or phyla (for example, between archea and bacteria, prokaryotes and eukaryotes) [23–28]. Here, we ask whether there exists a correlation between the degree of modularity of the metabolic network of an organism and the variability in its environment.
We analyzed the metabolic networks of 117 bacteria species living in broad range of habitats including oceans, salt lakes, thermal vents, soil and within hosts. The species were classified according to the degree of variability in their environment. We measured the modularity levels and additional related structural parameters of bacteria metabolic networks, finding that the level of variability in organism's environment correlates with the modularity of the networks: The more variable the environment, the more modular the metabolic network. Our study therefore supports the view that environmental variability promotes modularity in biological networks.
Classification of variability of bacterial environments
The natural environment of 117 bacterial species was classified based on the NCBI classification for bacterial lifestyle  (see methods). The classification includes six classes: Obligate bacteria[30, 31] that are obligately associated with a host, either intracellulary or extracellulary. An example is Buchnera that lives in symbiosis inside aphids and has little contact with the outside world. Specialized bacteria that live in specialized environments such as marine thermal vents. Aquatic bacteria, that live in fresh or seawater environment, and are not associated with hosts. Facultative bacteria, free living bacteria such as E. coli that often associate with a host. Multiple bacteria, that live in multiple different kinds of environments such as bacteria with a wide host range, and Terrestrial Bacteria, that live in the soil.
We ranked the variability of the different environments considering physical conditions (such as temperature, osmolarity, acidity, oxygen availability, etc) and heterogeneity of species . Obligate bacteria are thought to have the most constant environment as these bacteria live within a biochemically controlled and isolated environment , usually with only few other species or even no other species. Specialized and aquatic bacteria are adapted to a strict realm of ecological conditions, yet their habitat is less protected and exhibits higher species heterogeneity than the obligate class. Facultative bacteria that live both in hosts and in the outside world are exposed to a more variable environment than the former classes, but less variable than the multiple class that spans widely different habitats. Finally, the terrestrial class is often considered the most variable class , since the soil is highly heterogeneous and has a diverse ecology.
In addition to this qualitative ordering, we sought a quantitative measure that may reflect environmental variability. One such measure is the fraction of transcription factors out of the total number of genes in the organism. The reason for using this measure is that theoretical analysis based on cost-benefit analysis suggests that transcription factors are more strongly selected in variable environments than in constant ones (chapter 10 of Ref ). Hence, more regulators per gene are expected the more variable the environment. It was shown  that the benefit of a transcription factor is highest when it can appropriately regulate gene expression in response to environmental changes, offsetting the cost of production and maintenance of the transcription factor protein and associated sensory systems [3, 36]. Studies by Moran and colleagues indicate that under constant conditions, genes for transcription factors tend to be lost from the genome [37, 38]. Genes of organisms in nearly constant conditions are deleted, or are constitutively expressed and do not require transcription control, as was experimentally demonstrated for Buchnera  and Rickettsia .
In agreement with the above, we find that the fraction of transcription factors, as well as their total number, increases with the expected variability of the bacterial lifestyle classification (Fig 1a). Note also that the size of the metabolic networks (=number of metabolites) tends to increase with the expected variability of the environment (Fig 1b).
Modularity of metabolic networks correlates with variability in the environment
The metabolic networks of the different species were obtained from KEGG database. In the networks, each metabolite is a node and edges represent metabolic reactions (see Methods).
We used a standard measure of modularity [18, 41], to evaluate the modularity of each metabolic network. We find that modularity increases with the variability in the environment (Fig. 2). The lowest modularity is found for the obligate class, higher modularity is found for the specialized and aquatic classes, and the highest modularity is found for facultative, multiple and terrestrial classes (Correlation coefficient c = 0.63, p-value p < 10-4).
We further considered a control for the effects of network size on modularity, by comparing metabolic networks of the same size. This was achieved by constructing subnetworks of each network containing n = 60 nodes (metabolites), comparable to the smallest networks in the dataset. These sub-networks were constructed by contracting linear pathways and cycles, and removing dangling nodes, until a network of the required size was obtained. We find that the modularity of these equal-sized metabolic networks also significantly increased with environmental variability (c = 0.59, p < 10-4, additional file 2, section 1.3).
As a final control for the effect of network size, we computed the Pearson partial correlation  between modularity, environmental variability, and network size. We find that the correlation between modularity and variability is significant also when removing the effect of network size (c = 0.24, p = 0.02, see additional file 2 – section 1.4).
Each lifestyle class includes bacteria from different branches of the phylogenetic tree. We find low correlation between the similarity in the modularity of pairs of species and their distance on the phylogenetic tree (c = 0.1, p < 10-4, see additional file 2, section 3). Thus, modularity seems to be more correlated with environment than with evolutionary relatedness.
We find that networks from species in constant environments tend to be more tree-like than those in variable environments (Fig 3). This can be quantitatively seen by considering two topological measures, betweenness-centrality  and cyclic coefficients , both properly normalized to account for different network sizes. We find that the former increases and the latter decreases the more constant the environment, suggesting tree like structures (see additional file 2, sections 2.2–2.3). The tree-like structure seems to result from a lack of alternate metabolic paths in networks from constant environments.
Modules in organisms from variable environments are more functionally pure
We used the Newman-Girvan algorithm  to define structural modules in each network, and tested whether the identified modules correspond to well-defined metabolic functions . A structural module in a network was considered to be also a functional module if it was enriched for substrates that had a shared metabolic function such as: central metabolism, biosynthesis of amino acids, carbohydrates, lipids etc. (according to the KEGG pathway classification [44, 45]). For each network we evaluated two measures: functionality, defined as the fraction of structural modules that were significantly enriched for a metabolic function, and coverage – the fraction of biological functions that could be mapped to structural modules.
We find that the modules in networks of organisms from constant environments usually do not correspond to a defined metabolic function, but rather to mixtures of several biological categories. In contrast, the modules of networks from variable environment classes usually corresponded to a unique function and the majority of the metabolic functions could be assigned to at least one structural module (Fig. 4 and section 4 in additional file 2).
These results further indicate that the present approach for identifying modules yields modules that have biologically significant function .
This study indicates that variability in the environment correlates with enhanced modular organization of metabolic networks, while constant environment correlate with a less modular structure.
One interpretation of these findings can be made in the context of previous simulation studies of evolution in modularly varying environments . The metabolic goal that a bacterium faces can be considered as a combination of sub-goals. An example of a sub-goal is the biosynthesis of an amino acid such as histidine. If histidine is missing in the environment, the bacterium must synthesize it. If histidine is present, the bacterium can down-regulate the biosynthesis pathway and instead import this metabolite. When the environment changes over time, it introduces a different combination of such metabolic sub-goals.
Simulations suggested that varying the sub-goals leads to the evolution of networks with a modular structure, where each module corresponds to one of the sub-goals . Modular structure evolves despite the fact that it is less optimal than non-modular solutions . In contrast, evolution under a goal that is constant over time leads to non-modular networks, in which many nodes participate in several functions . The present findings may be interpreted within this context: Bacteria that live under varying environments typically evolve a functional module for each of the varying sub-goals. Bacteria under constant conditions tend to evolve towards a less modular design.
It is interesting to note that some metabolic goals are held relatively constant even when the environment changes. An example is energy metabolism, which is needed for growth by all of the bacterial species studied, in all environments. Analysis of the metabolic networks shows that the part of the network responsible for energy metabolism (central metabolism) is less modular than other parts of the network (such as biosynthesis of amino acids, nucleotides, vitamins etc) [44, 47–49]. More generally, the fraction of the metabolic network devoted to constant goals (such as central metabolism) seems to increases as the environment becomes more constant (Fig S8c in additional file 2).
An additional observation in the computer simulations  is that initially modular networks rapidly degrade into non-modular but more optimal structures when the goal becomes constant over time. Examples of such a degeneration of modularity can be seen by comparing the closely related species E. coli and Buchnera. E. coli lives in a variable environment, moving between its mammalian host and the external world. Buchnera lives in a more constant environment, as an endosymbiont of aphids. Buchnera is found to fuse together pathways that are separate in E. coli, and thus to achieve its metabolic goals with a smaller set of enzymes . One example occurs in the histidine and purine modules (Fig 5). Both modules convert the metabolic substrates PRPP to AICAR in two distinct parhways. E. coli seems to maintain these alternative pathways because under different environments (histidine/purine rich environments) only one of the pathways is utilized. For Buchnera, on the other hand, as an endosymbiont that supplies amino-acids to its host, histidine biosynthesis is a fixed goal and under no regulation. Here, the purine module can count on the histidine module for AICAR production. The two pathways were thus combined into a single module, in which many of the genes are used for both functions [38, 51]. It would be interesting to uncover other mechanisms that degenerate or enhance modularity by comparing networks of closely related species with different environments.
One limitation of the present study is the limited knowledge of metabolic networks for diverse species. The reconstructed networks, based on genomic data, were used to generate information about putative non-directed metabolic interactions. The present network representation ignore: i) directionality of reactions ii) reaction stochiometry iii) that only a fraction of the reactions are active under given environmental conditions (hence at best it offers only a static view on modularity). The above mentioned problems can be handled by more sophisticated network analyses . Such studies employ correlated reaction sets as mathematically defined modules in biochemical reaction networks. They constitute groups of reactions in a network that always appear together in functional states of that network and therefore represent a functional module of the reaction network. Previous work has shown that these sets can include non-obvious groups of reactions and differ from groupings of reactions based on structural analyses of network topology [53, 54]. One drawback of these latter approaches is that they require carefully annotated, genome-scale metabolic network, of which is only available for a handful number of species.
This study indicates that the modularity of metabolic networks correlates with the variability of the environment. Such a correlation supports the view that variability in the natural habitat promotes modularity. It would be important to test this more fully as data on metabolic and regulatory networks of diverse species becomes more complete.
We currently know more about the structure of metabolic networks than about the ecology of the organisms. It is a challenge to see how far one can go in what might be termed 'reverse ecology' : inferring from the structure of biological system information about the environment in which it evolved.
Lifestyles of bacterial species
The present classification of species can be found at additional file 1 and also at the corresponding author website. The bacterial environmental data is based on . The original classification was adjusted in the present study, to better reflect environmental variability. The main adjustment was the split of the "Host Associated" class into obligate bacteria and facultative bacteria. Further corrections of specific bacteria species were performed based on literature search. For example, vibrio species were changed from aquatic to facultative because they are pathogens of fish. The numbers of species in the classification were as following: 35 obligate bacteria, 5 specialized bacteria, 4 aquatic bacteria, 42 facultative bacteria, 28 multiple bacteria and 3 terrestrial bacteria.
Metabolic networks construction
Reconstructed metabolic networks of 117 bacteria were taken from the KEGG Database . Each network was represented as a substrate graph, where each node corresponds to a metabolite and an edge corresponds to a reaction. Highly connected metabolites (such as ATP, NADH, H20 etc) were removed – a crucial step for topological analysis . Similar results were found also for a bipartite graph representation including both reactions and metabolites as nodes (additional file 2, section 1.2). Analysis was preformed on the giant connected component of the networks, to avoid bias from small isolated components. Similar results were found also when the entire network was analyzed (additional file 2, section 1.1).
Quantitative measure of network modularity
To quantify the modularity of a network we used the normalized Q m measure of Kashtan et al. . This measure is based on the commonly used Newmann and Girvan modularity measure Q , defined as the fraction of edges that lie within modules rather than between modules relative to that expected by chance.
A module, therefore, is as a group of nodes with many interactions between them, and few interactions to the other nodes. The algorithm attempts to find the division of nodes into modules that maximizes Q. This measure sums over all modules in the network and hence scales with network size. To allow comparison of the modularity of networks with different size and connectivity, this parameter needs to be further normalized :
Q m = (Q real - Q rand )/(Q max - Q rand ) (1)
Where Q real is the Q value of the network, Q rand is the mean Q value of randomized networks, and Q max is the upper bound of Q for a given network ensemble. We obtained Q rand by averaging Q over hundred random networks that preserve the degree distribution of the real network . Q max was estimated as 1 - 1/M  where M is the number of modules in the real network (usually the lower bound for the number of modules in the random network). For the present substrate networks, M ranged from 4 (for small networks ~100 nodes) to 22 (for large networks ~400 nodes).
Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature. 1999, 402: C47-52. 10.1038/35011540.
Wagner G, Schlosser G: Modularity in Development and Evolution. 2004, University of Chicago Press
Alon U: An Introduction to Systems Biology: Design Principles of Biological circuits. 2006, Chapman & Hall/CRC, 1
Stelling J, Klamt S, Bettenbrock K, Schuster S, Gilles ED: Metabolic network structure determines key aspects of functionality and regulation. Nature. 2002, 4 (20): 190-193. 10.1038/nature01166.
Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297: 1551-1555. 10.1126/science.1073374.
Wolf DM, Arkin AP: Motifs, modules and games in bacteria. Curr Opin Microbiol. 2003, 6: 125-134. 10.1016/S1369-5274(03)00033-X.
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: 15428-15433. 10.1073/pnas.2136809100.
Ma HW, Zhao XM, Yuan YJ, Zeng AP: Decomposition of metabolic network into functional modules based on the global connectivity structure of reaction graph. Bioinformatics. 2004, 20: 1870-1876. 10.1093/bioinformatics/bth167.
Stelling J: Mathematical models in microbial systems biology. Curr Opin Microbiol. 2004, 7: 513-518. 10.1016/j.mib.2004.08.004.
Zhao J, Yu H, Luo JH, Cao ZW, Li YX: Hierarchical modularity of nested bow-ties in metabolic networks. BMC Bioinformatics. 2006, 7: 386-10.1186/1471-2105-7-386.
Kurata H, El-Samad H, Iwasaki R, Ohtake H, Doyle JC, Grigorova I, Gross CA, Khammash M: Module-based analysis of robustness tradeoffs in the heat shock response system. PLoS Comput Biol. 2006, 2: e59-10.1371/journal.pcbi.0020059.
Bhattacharyya RP, Remenyi A, Yeh BJ, Lim WA: Domains, Motifs, and Scaffolds: The Role of Modular Interactions in the Evolution and Wiring of Cell Signaling Circuits. Annu Rev Biochem. 2006
Wagner GP, Altenberg L: Complex Adaptations and the Evolution of Evolvability. Evolution. 1996, 50: 967-976. 10.2307/2410639.
Schlosser G, Thieffry D: Modularity in development and evolution. Bioessays. 2000, 22: 1043-1045. 10.1002/1521-1878(200011)22:11<1043::AID-BIES11>3.0.CO;2-C.
Lipson H, Pollack JB, Suh NP: Promoting Modularity In Evolutionary Design. Proceedings of DETC'01 2001 ASME Design Engineering Technical Conferences; Pittsburgh, Pennsylvania, USA. 2001
Lipson H, Pollack JB, Suh NP: On the origin of modular variation. Evolution Int J Org Evolution. 2002, 56: 1549-1556.
Variano EA, McCoy JH, Lipson H: Networks, dynamics, and modularity. Phys Rev Lett. 2004, 92: 188-701. 10.1103/PhysRevLett.92.188701.
Kashtan N, Alon U: Spontaneous evolution of modularity and network motifs. Proc Natl Acad Sci USA. 2005, 1 (02): 13773-13778. 10.1073/pnas.0503610102.
Lenski RE, Ofria C, Pennock RT, Adami C: The evolutionary origin of complex features. Nature. 2003, 423: 139-144. 10.1038/nature01568.
Fell DA, Wagner A: The small world of metabolism. Nat Biotechnol. 2000, 18: 1121-1122. 10.1038/81025.
Karp PD, Riley M, Paley SM, Pelligrini-Toole A: EcoCyc: an encyclopedia of Escherichia coli genes and metabolism. Nucleic Acids Res. 1996, 24: 32-39. 10.1093/nar/24.1.32.
Karp PD, Riley M, Paley SM, Pellegrini-Toole A: The MetaCyc Database. Nucleic Acids Res. 2002, 30: 59-61. 10.1093/nar/30.1.59.
Ma H, Zeng AP: Reconstruction of metabolic networks from genome data and analysis of their global structure for various organisms. Bioinformatics. 2003, 19: 270-277. 10.1093/bioinformatics/19.2.270.
Zhu D, Qin ZS: Structural comparison of metabolic networks in selected single cell organisms. BMC Bioinformatics. 2005, 6: 8-10.1186/1471-2105-6-8.
Ebenhoh O, Handorf T, Heinrich R: A cross species comparison of metabolic network functions. Genome Inform. 2005, 16: 203-213.
Zhang Y, Li S, Skogerbo G, Zhang Z, Zhu X, Sun S, Lu H, Shi B, Chen R: Phylophenetic properties of metabolic pathway topologies as revealed by global analysis. BMC Bioinformatics. 2006, 7: 252-10.1186/1471-2105-7-252.
Tun K, Dhar PK, Palumbo MC, Giuliani A: Metabolic pathways variability and sequence/networks comparisons. BMC Bioinformatics. 2006, 7: 24-10.1186/1471-2105-7-24.
Forst CV, Flamm C, Hofacker IL, Stadler PF: Algebraic comparison of metabolic networks, phylogenetic inference, and metabolic innovation. BMC Bioinformatics. 2006, 7: 67-10.1186/1471-2105-7-67.
Entrez Genome Project. [http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi]
Nilsson AI, Koskiniemi S, Eriksson S, Kugelberg E, Hinton JC, Andersson DI: Bacterial genome size reduction by experimental evolution. Proc Natl Acad Sci USA. 2005, 102: 12112-12116. 10.1073/pnas.0503654102.
Hosokawa T, Kikuchi Y, Nikoh N, Shimada M, Fukatsu T: Strict host-symbiont cospeciation and reductive genome evolution in insect gut bacteria. PLoS Biol. 2006, 4: e337-10.1371/journal.pbio.0040337.
Robertson BR, Button DK: Characterizing aquatic bacteria according to population, cell size, and apparent DNA content by flow cytometry. Cytometry. 1989, 10: 70-76. 10.1002/cyto.990100112.
Xu J: Microbial ecology in the age of genomics and metagenomics: concepts, tools, and recent advances. Mol Ecol. 2006, 15: 1713-1731. 10.1111/j.1365-294X.2006.02882.x.
Ochman H, Moran NA: Genes lost and genes found: evolution of bacterial pathogenesis and symbiosis. Science. 2001, 292: 1096-1099. 10.1126/science.1058543.
Bardgett RD: Causes and consequences of biological diversity in soil. Zoology (Jena). 2002, 105: 367-374.
Dekel E, Mangan S, Alon U: Environmental selection of the feed-forward loop circuit in gene-regulation networks. Phys Biol. 2005, 2: 81-88. 10.1088/1478-3975/2/2/001.
Moran NA: Tracing the evolution of gene loss in obligate bacterial symbionts. Curr Opin Microbiol. 2003, 6: 512-518. 10.1016/j.mib.2003.08.001.
Moran NA, Dunbar HE, Wilcox JL: Regulation of transcription in a reduced bacterial genome: nutrient-provisioning genes of the obligate symbiont Buchnera aphidicola. J Bacteriol. 2005, 187: 4229-4237. 10.1128/JB.187.12.4229-4237.2005.
Gil R, Sabater-Munoz B, Latorre A, Silva FJ, Moya A: Extreme genome reduction in Buchnera spp.: toward the minimal genome needed for symbiotic life. Proc Natl Acad Sci USA. 2002, 99: 4454-4458. 10.1073/pnas.062067299.
Andersson SG, Zomorodipour A, Andersson JO, Sicheritz-Ponten T, Alsmark UC, Podowski RM, Naslund AK, Eriksson AS, Winkler HH, Kurland CG: The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature. 1998, 396: 133-140. 10.1038/24094.
Newman ME: Fast algorithm for detecting community structure in networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2004, 69: 066-133.
Cramer D: A cautionary tale of two statistics: Partial correlation and standardixed partial regression. Journal of Psychology. 2003, 137 (5): 507-511.
Kim HJ, Kim JM: Cyclic topology in complex networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2005, 72: 036-109.
Guimera R, Nunes Amaral LA: Functional cartography of complex metabolic networks. Nature. 2005, 433: 895-900. 10.1038/nature03288.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.
Fischer E, Sauer U: Large-scale in vivo flux analysis shows rigidity and suboptimal performance of Bacillus subtilis metabolism. Nat Genet. 2005, 37 (63): 640-6.
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.
Dandekar T, Schuster S, Snel B, Huynen M, Bork P: Pathway alignment: application to the comparative analysis of glycolytic enzymes. Biochem J. 1999, 343 (Pt 1): 115-124. 10.1042/0264-6021:3430115.
Csete M, Doyle J: Bow ties, metabolism and disease. Trends Biotechnol. 2004, 22: 446-450. 10.1016/j.tibtech.2004.07.007.
Pal C, Papp B, Lercher MJ, Csermely P, Oliver SG, Hurst LD: Chance and necessity in the evolution of minimal metabolic networks. Nature. 2006, 440: 667-670. 10.1038/nature04568.
Zientz E, Dandekar T, Gross R: Metabolic interdependence of obligate intracellular bacteria and their insect hosts. Microbiol Mol Biol Rev. 2004, 68: 745-770. 10.1128/MMBR.68.4.745-770.2004.
Papin JA, Reed JL, Palsson BO: Hierarchical thinking in network biology: the unbiased modularization of biochemical networks. Trends Biochem Sci. 2004, 29: 641-647. 10.1016/j.tibs.2004.10.001.
Papin JA, Price ND, Wiback SJ, Fell DA, Palsson BO: Metabolic pathways in the post-genome era. Trends Biochem Sci. 2003, 28: 250-258. 10.1016/S0968-0004(03)00064-1.
Papin JA, Price ND, Palsson BO: Extreme pathway lengths and reaction participation in genome-scale metabolic networks. Genome Res. 2002, 12: 1889-1900. 10.1101/gr.327702.
Newman ME: Modularity and community structure in networks. Proc Natl Acad Sci USA. 2006, 103: 8577-8582. 10.1073/pnas.0601602103.
Guimera R, Sales-Pardo M, Amaral LA: Modularity from fluctuations in random graphs and complex networks. Phys Rev E. 2004, 70: 025-101. 10.1103/PhysRevE.70.025101.
Ziv E, Middendorf M, Wiggins CH: Information-theoretic approach to network modularity. Phys Rev E. 2005, 71: 046-117.
We thank S.G Oliver, T. Shlomi, E. Sharon and all the members of our lab for discussions. We thank the NIH and the Kahn Family Foundation for support.
The author(s) declares that there are no competing interests.
All authors designed the research and analyzed the results. MP performed the computational assays. All authors contributed to writing the manuscript.
Electronic supplementary material
Additional file 1: Microbial Environment database. Refined classification of 117 bacteria into 6 environmental classes. (XLS 25 KB)
Additional file 2: Supplementary online material (SOM). The SOM supplies additional networks analysis, control for network size effect and descriptions for reduced networks' construction. (PDF 339 KB)
About this article
Cite this article
Parter, M., Kashtan, N. & Alon, U. Environmental variability and modularity of bacterial metabolic networks. BMC Evol Biol 7, 169 (2007). https://doi.org/10.1186/1471-2148-7-169
- Metabolic Network
- Network Size
- Modular Structure
- Environmental Variability