Environmental variability and modularity of bacterial metabolic networks
© Parter et al; licensee BioMed Central Ltd. 2007
Received: 17 May 2007
Accepted: 23 September 2007
Published: 23 September 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 .
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 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.
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.
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).
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).
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.
- Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature. 1999, 402: C47-52. 10.1038/35011540.View ArticlePubMedGoogle Scholar
- Wagner G, Schlosser G: Modularity in Development and Evolution. 2004, University of Chicago PressGoogle Scholar
- Alon U: An Introduction to Systems Biology: Design Principles of Biological circuits. 2006, Chapman & Hall/CRC, 1Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Wolf DM, Arkin AP: Motifs, modules and games in bacteria. Curr Opin Microbiol. 2003, 6: 125-134. 10.1016/S1369-5274(03)00033-X.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Stelling J: Mathematical models in microbial systems biology. Curr Opin Microbiol. 2004, 7: 513-518. 10.1016/j.mib.2004.08.004.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 2006Google Scholar
- Wagner GP, Altenberg L: Complex Adaptations and the Evolution of Evolvability. Evolution. 1996, 50: 967-976. 10.2307/2410639.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Lipson H, Pollack JB, Suh NP: Promoting Modularity In Evolutionary Design. Proceedings of DETC'01 2001 ASME Design Engineering Technical Conferences; Pittsburgh, Pennsylvania, USA. 2001Google Scholar
- Lipson H, Pollack JB, Suh NP: On the origin of modular variation. Evolution Int J Org Evolution. 2002, 56: 1549-1556.View ArticleGoogle Scholar
- Variano EA, McCoy JH, Lipson H: Networks, dynamics, and modularity. Phys Rev Lett. 2004, 92: 188-701. 10.1103/PhysRevLett.92.188701.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Lenski RE, Ofria C, Pennock RT, Adami C: The evolutionary origin of complex features. Nature. 2003, 423: 139-144. 10.1038/nature01568.View ArticlePubMedGoogle Scholar
- Fell DA, Wagner A: The small world of metabolism. Nat Biotechnol. 2000, 18: 1121-1122. 10.1038/81025.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Ebenhoh O, Handorf T, Heinrich R: A cross species comparison of metabolic network functions. Genome Inform. 2005, 16: 203-213.PubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ochman H, Moran NA: Genes lost and genes found: evolution of bacterial pathogenesis and symbiosis. Science. 2001, 292: 1096-1099. 10.1126/science.1058543.View ArticlePubMedGoogle Scholar
- Bardgett RD: Causes and consequences of biological diversity in soil. Zoology (Jena). 2002, 105: 367-374.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Newman ME: Fast algorithm for detecting community structure in networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2004, 69: 066-133.Google Scholar
- Cramer D: A cautionary tale of two statistics: Partial correlation and standardixed partial regression. Journal of Psychology. 2003, 137 (5): 507-511.View ArticlePubMedGoogle Scholar
- Kim HJ, Kim JM: Cyclic topology in complex networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2005, 72: 036-109.Google Scholar
- Guimera R, Nunes Amaral LA: Functional cartography of complex metabolic networks. Nature. 2005, 433: 895-900. 10.1038/nature03288.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Csete M, Doyle J: Bow ties, metabolism and disease. Trends Biotechnol. 2004, 22: 446-450. 10.1016/j.tibtech.2004.07.007.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Newman ME: Modularity and community structure in networks. Proc Natl Acad Sci USA. 2006, 103: 8577-8582. 10.1073/pnas.0601602103.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ziv E, Middendorf M, Wiggins CH: Information-theoretic approach to network modularity. Phys Rev E. 2005, 71: 046-117.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.