Module detection in complex networks using integer optimisation
- Gang Xu^{1},
- Laura Bennett^{2},
- Lazaros G Papageorgiou^{1} and
- Sophia Tsoka^{2}Email author
https://doi.org/10.1186/1748-7188-5-36
© Xu et al; licensee BioMed Central Ltd. 2010
Received: 27 April 2010
Accepted: 12 November 2010
Published: 12 November 2010
Abstract
Background
The detection of modules or community structure is widely used to reveal the underlying properties of complex networks in biology, as well as physical and social sciences. Since the adoption of modularity as a measure of network topological properties, several methodologies for the discovery of community structure based on modularity maximisation have been developed. However, satisfactory partitions of large graphs with modest computational resources are particularly challenging due to the NP-hard nature of the related optimisation problem. Furthermore, it has been suggested that optimising the modularity metric can reach a resolution limit whereby the algorithm fails to detect smaller communities than a specific size in large networks.
Results
We present a novel solution approach to identify community structure in large complex networks and address resolution limitations in module detection. The proposed algorithm employs modularity to express network community structure and it is based on mixed integer optimisation models. The solution procedure is extended through an iterative procedure to diminish effects that tend to agglomerate smaller modules (resolution limitations).
Conclusions
A comprehensive comparative analysis of methodologies for module detection based on modularity maximisation shows that our approach outperforms previously reported methods. Furthermore, in contrast to previous reports, we propose a strategy to handle resolution limitations in modularity maximisation. Overall, we illustrate ways to improve existing methodologies for community structure identification so as to increase its efficiency and applicability.
Background
Networks - i.e. groups of entities (nodes or vertices) pairs of which are linked through a form of common property (edges or links) - have formed an efficient representation framework for a variety of complex systems such as social groupings and internet connectivity [1]. The analysis of biological data in systems biology studies through the formalisms of network theory have received particular attention recently, due to the potential benefits that such methodologies can confer in mining the intricate relationships in metabolic networks [2–4], signaling pathways [5], gene regulatory networks [6] or other forms of protein interactions [7]. In general, the abstractions offered by graph theory representations (i) facilitate the analysis of network performance, (ii) provide a unifying framework for comparisons of features across different systems and (iii) assist the mathematical characterisation of system properties and dynamics.
Topological properties of networks are particularly important in revealing the organisational principles of nodes within the context of the entire system [8]. Community structures or modules are defined when a larger density of links exists within a specific part of the network than outside it [9]. Each of such modules can be regarded as a discrete entity whose function or properties are in some way separable from other modules. Modular structure underlies (i) the adaptability of a system to new conditions [10] and (ii) the robustness (or conversely the vulnerability) of the system to external attack or other form of change in topological features [11, 12]. The analysis of pairwise or even longer-range relationships in networks can reveal how preferential attachment of new nodes influences community structure [13, 14], giving rise to small-world or scale-free architectures [15].
In light of the above, the detection of modules and the analysis of community structure in networks has the potential to reveal the design principles of complex systems and provide important insights into how such systems are organised, how they evolve and how their components interact. For example, in biological networks, the analysis of enzyme connectedness may reveal participation in the same biological pathway, and module detection in protein interactions reflects protein function type or evolutionary properties [3, 7, 16]. Importantly, the characterisation of gene products with previously unknown functional properties through community detection has the potential to aid function assignment. Recent reviews have reported on the community detection problem in comprehensive manner [17, 18].
The two major avenues to detect community structure have been graph partitioning [19, 20] and hierarchical clustering methods [16, 20–23]. Major disadvantage in the case of graph partitioning is the absence of a termination criterion in the bisection process, while in hierarchical clustering there is no clear indication of where the tree should be split to yield the optimal partitioning. Such shortcomings result in either sub-optimal partitioning or unsatisfactory implementations for large networks. In addition to graph partitioning and hierarchical methods, which are particularly suitable for standard partitioning (each node belongs to a single community), other methodologies exist to detect overlapping communities as nodes may belong to several communities, for example, the clique percolation method, [24].
An important breakthrough in the community detection problem has taken the form of a quantitative measure to express the quality of community presence, namely modularity. Network modularity is defined as the fraction of all edges that lie within communities minus the expected value of the same quantity in a graph in which the vertices have the same degrees but edges are placed randomly [25–27]. Usually, in our experience, network modularity values of around 0.4-0.8 indicate strong community presence. Use of the modularity metric has transformed the community structure identification problem into an optimisation task where community structures can be determined by maximising the network modularity through various optimisation techniques [25].
As modularity optimisation is NP-hard [28, 29], efficient algorithms to find the maximum modularity values are unlikely to exist. Therefore, most approaches employ heuristics that aim at finding near-optimal solutions with modest computational cost. Crucial considerations in assessing the performance of modularity optimisation approaches are: (i) the scale and optimality handled by modularity optimisation methods and (ii) the resolution limit problem for small-size modules in large networks.
First, there seems to be a trade-off between network size and optimality achieved through modularity optimisation. Specifically, methods that guarantee global optimal solutions for modularity maximisation are able to operate only in small to medium-sized networks [30]. Divisive algorithms [25] were found to be prohibitively computationally expensive for large networks. On the other hand, methods that can be used on large networks, such as stochastic optimisation through simulated annealing [3, 31] and extremal optimisation [32], may yield sub-optimal solutions and so may suffer poor performance. In our own work, we have previously reported a rigorous mixed integer quadratic programming (MIQP) formulation to optimise the modularity metric with a set of linear constraints and mixed binary/continuous optimisation variables [30]. Due to the convexity properties of the model, global optimal solutions are achieved through the standard branch-and-bound procedure with commercial optimisation solvers, but use of this optimisation framework is limited to small-medium scale networks due to NP-hardness.
Second, doubts have been raised over the use of modularity optimisation for community detection recently, due to the observation that such procedures can reach a resolution limit [33]. This effect essentially implies that modules smaller than a specific scale are not detected, as the optimisation process combines smaller communities into larger ones in order to achieve better modularity. Some remedial procedures have been suggested through re-optimising each module [33, 34], tuning a resolution parameter [35], or implementing quantitative measures other than modularity [36].
Here, we aim to enhance the application of mathematical programming to community structure identification by: (i) developing an efficient methodology for module detection that is capable of handling large size networks and (ii) incorporating strategies for dealing systematically with the problem of a resolution limit in module detection through modularity optimisation approaches. Below, a two-stage solution approach for community identification using mathematical programming is described, the resolution limit in modularity optimisation is addressed via the introduction of an iterative procedure and the applicability of the proposed approaches is demonstrated through a number of network examples and comparisons with literature.
Methods
Stage 1: Initial network partition
where L_{ m }denotes the number of links in module m and D_{ m }is the degree of all nodes in module m. The modularity metric, Q, measures the difference between the fraction of links within communities and the expected fraction values when links are allocated randomly [25, 26]. The objective function employed here is the maximisation of the network modularity metric shown in equation (1).
where Y_{ nm }is a binary variable taking the value of 1 if node n is allocated to module m; 0 otherwise.
where CN_{ n }is the set of nodes e connected to node n.
Since global optimality of non-convex MINLP models cannot be guaranteed, different initial solutions are tested and the partition with the largest value of Q, is chosen as the best division from the set of candidate solutions. MINLP_Mod is performed for a given number of runs, N^{ max }, from random initial points and the node-module allocation with the maximum modularity value is stored and denoted by set I_{ m }. Using N^{ max }= 100 provides a good representation of solution space.
Stage 2: Iterative improvement of network partition
Having selected a node-module association with maximum modularity from the previous stage (i.e. I_{ m }), module allocation may be improved further through an iterative fixing and releasing scheme. The general idea is to solve a reduced MIQP formulation of modularity optimisation that was previously proposed in [28]. Most of the Y_{ nm }variables are fixed, which reduces the number of variables, thus resulting in a more tractable model. Sets of nodes are released in the sense that they are free to be re-allocated to a different module in subsequent executions of the MIQP model.
This stage initially adopts the node-module allocation obtained from stage 1 by fixing all the relevant Y_{ nm }binary variables in I_{ m }to the value of one. For the first module, the set of nodes in the module is denoted as RE_{ m }and the set of nodes to be released (or 'un-fixed') is denoted as Δ, where the size of Δ is N_R, a value chosen according to criteria described below. The reduced MIQP (OptMod) is solved (for details see [30]), with all nodes in Δ released and all other nodes fixed. I_{ m }is then updated with the solution from the reduced MIQP.
The above scheme is applied sequentially for remaining modules, which completes one round of the major improvement iteration, k, with network modularity value, Q_{ k }. The same strategy starts again, retaining the order of the modules, until no improvement of the modularity value is reported for two successive major iterations.
Comparing the single-level MIQP model, OptMod, to the reduced MIQP models as implemented here, the latter strategy involves fewer variables and constraints and can be terminated efficiently even in cases of larger size networks, as discussed in the Results and Discussion section. To justify why an iterative reduced MIQP is preferred over MINLPs, it should be mentioned that we achieved improved solutions by solving a series of reduced MIQP models, while no improvements have been observed when solving reduced MINLP models iteratively.
where Aver_{ m }denotes the average degree in module m without considering the inter-module links and U is a user-defined parameter. Here, we used a value of U = 200 which was shown to provide satisfactory results for all examples studied. As a result, the actual number of released nodes, N_R, will be the smaller value between MAX_{ m }^{ r }and the number of remaining nodes to be released in RE_{ m }(i.e. $N\_R=\mathrm{min}\{\left|R{E}_{m}\right|,MA{X}_{m}^{r}\}$).
In other words, if the number of nodes in module m is greater than $MA{X}_{m}^{r}$, the first N_R nodes, Δ, in module m will be released and the reduced MIQP solved. I_{ m }is updated and RE_{ m }becomes RE_{ m }|Δ. If the updated RE_{ m }is still greater than $MA{X}_{m}^{r}$, a further set of nodes of size N_R is released, otherwise all remaining nodes are released. The reduced MIQP is solved once again and I_{ m }and RE_{ m }updated accordingly. This is repeated until all nodes in the module have been released at one point and the procedure moves on to the next module. In order to determine the set of N_R nodes to be released in modules, we use a simple rule by first sorting nodes with non-decreasing indices and then assigning higher priority to nodes with smaller indices.
The above scheme is applied to all modules detected during Stage 1, with the sequential order of the modules maintained throughout the whole procedure. Future research can investigate the effect of changing the sequence of module reallocation, the appropriate selection of MAX_{ m }^{ r }and node prioritisation. Figure 1 illustrates the entire module detection strategy, iMod, encompassing Stages 1 and 2 of the mathematical programming algorithm reported above.
Procedure to address resolution limitations
Although the modularity metric has been widely accepted as a standard measure to quantify the community composition in networks and detect modules, resolution limit problems can hinder its application. Such effects entail the failure of modularity optimisation to detect modules smaller than a scale which depends on the size of the network and the degree of inter-connectedness of the modules, as the algorithm tends to merge small modules to achieve larger modularity values [33, 37, 38]. Methodologies that aim to overcome resolution limits can provide deeper insights into finer structures of modules in complex networks and a more accurate depiction of community structure on the basis of the modularity measure.
In this section, we report a solution procedure (ResMod) that allows smaller modules that may not be detected in the initial modularity optimisation to become apparent. First, the two-stage approach for module detection, iMod, is applied to the whole network to obtain a partition into several modules. In order to determine if these modules comprise smaller modules, each module is considered as a disjoint subnetwork, ignoring links with other modules, and iMod is then applied once to each subnetwork.
The partition of the subnetwork into smaller modules is accepted as part of the community structure of the original network if its modularity as a disconnected entity (i.e. only considering the links involved in the subnetwork) is greater than an enforced threshold. If the partition of the subnetwork yields a value less than this threshold, the new decomposition is not accepted and the subnetwork remains intact as a community of the original network.
Here, a threshold of value of 0.3 is adopted as a representative community structure indicator, in accordance to previous reports [25, 26, 33, 34]. This criterion is implemented to avoid over-partitioning that may hinder method applicability. We should note here that this implementation of a single and unvarying threshold to determine whether partitioning is required may not be enough to capture cases where random graphs (or partitions obtained by chance) have a modularity higher than 0.3. Ideally, this criterion should be complemented with an estimate of the statistical significance of the modularity achieved (see [4, 34]) to ensure that this value is above a fluctuation margin. However, in practice, even this coarse-grained approach to resolution limitation problems seems to work well in proposing finer community structures for common complex networks and it is a good first step into research for improving modularity maximisation methods.
Results and Discussion
Computational results comparing the performance of modularity optimisation methodologies across several network examples.
Networks | iMod | EB | EIG | C_{3}/C_{4} | EO | SA | QCUT | Greedy | ||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Name | N | L | Median Q | Best Q | M | Q | ||||||
Zachary | 34 | 78 | 0.420 | 0.420 | 4 | 0.401 | 0.419 | 0.417 | 0.419 | 0.420 | 0.419 | |
Dolphin | 62 | 159 | 0.529 | 0.529 | 5 | 0.520 | 0.518 | 0.519 | ||||
Les Miserables | 77 | 254 | 0.560 | 0.560 | 6 | 0.540 | 0.560 | 0.556 | ||||
P53 | 104 | 226 | 0.535 | 0.535 | 7 | 0.522 | 0.531 | |||||
Jazz | 198 | 2742 | 0.445 | 0.445 | 4 | 0.405 | 0.442 | 0.441 | 0.445 | 0.445 | 0.443 | |
E. coli | 418 | 519 | 0.780 | 0.781 | 19 | 0.766 | 0.752 | 0.776 | 0.779 | |||
S. cerevisiae | 688 | 1079 | 0.768 | 0.775 | 25 | 0.759 | 0.740 | 0.766 | 0.764 | |||
C. elegans | 453 | 2025 | 0.451 | 0.453 | 9 | 0.403 | 0.435 | 0.422 | 0.434 | 0.433 | 0.441 | |
1133 | 5451 | 0.575 | 0.580 | 9 | 0.532 | 0.572 | 0.567 | 0.574 | 0.576 | 0.543 |
A number of networks identified from the literature serve as test cases to showcase the efficiency of the computational methodology. Table 1 summarises all networks considered, their sizes and indicative results of the methodologies tested. Overall, nine examples were used with varying sizes, in terms of total number of nodes and links. These cases are inspired from social or biological relationships and represent well-studied cases in network analysis and related algorithm development.
Our methodologies for community detection and resolution limitations are compared against the most widely used approaches that employ modularity maximisation. Here, a brief account of such previously developed algorithms is given, together with reference to the original publications for more details on the main properties of each algorithm.
An algorithm based on edge-betweenness (EB) [25] that involves the iterative removal of edges with the highest betweenness score to split the network into communities has been one of the very first attempts to use modularity maximisation for module detection. The eigenvector approach (EIG) was later proposed by the same group, where network modularity was rewritten as eigenvectors of a modularity matrix and lead to a spectral algorithm for community detection [9]. Edge-betweenness has recently been extended through the use of edge weights defined by the edge-clustering coefficient (C_{3}/C_{4}) to improve module detections [51]. Popular optimisation methodologies have been proposed as efficient means to achieve modularity maximisation, namely extremal optimisation (EO) [32] and simulated annealing (SA) [33]. Recently, heuristic algorithms have been proposed, i.e. one that relies on spectral graph partitioning and local search (QCUT) [34] and a greedy method for iterative grouping of nodes into communities (Greedy) [52]. Both of these methods show good performance compared to previous approaches. For Greedy and QCUT, we used the relevant software to evaluate modules and estimate the resulting modularity. For all other methods, the reported results are taken from the relevant published papers.
Extensive comparisons of performance across all above methodologies show that iMod achieved network partitions with the highest modularity (Table 1). Consistently better performance was noted for iMod throughout all examples studied. It is important to mention that even small improvements in modularity can differentiate between good and exceptional methods, as has been noted previously [9].
For the example of the p53 network (Figure 2), modules were mapped onto KEGG pathways and pathway enrichment was calculated against the human genome through SubpathwayMiner [53]. We compared enriched pathways for the iMod and Greedy partitions and, even for small differences in community structures detected, more pathways were significantly enriched in the iMod partition. Even though clearly more work is need along these lines, this is an early indication that module detection through iMod may be more meaningful biologically.
Comparative analyses are hindered to some extent by missing values in Table 1, as different network examples were assessed through each of the reported methodologies. For instance, the p53 example has been implemented in three methods (iMod, QCUT and Greedy), the Dolphin and Les Miserables networks were considered by four methodologies (iMod, EB, QCUT and Greedy), the E. coli and S. cerevisiae by five (iMod, EIG, SA, QCUT and Greedy) and the remaining four networks have been tested by different combinations of seven methodologies out of eight community detection algorithms considered in total. Such missing values indicate an impediment in related comparison efforts and it is suggested that the definition of network examples as standards, where algorithm development and evaluation can be benchmarked, is needed in order to facilitate and improve comparative analyses [54]. However, it should also be noted that this is one of the most comprehensive comparisons of module detection methodologies employing modularity maximisation, to our knowledge.
Benchmarking was also extended to simulated networks to illustrate the efficiency of iMod. A large number of artificial networks with known community structure was generated, as described previously [25]. These synthetic networks comprise 128 nodes and are partitioned into four communities of 32 nodes with degree equal to 16. In addition, we considered the case where degree was set equal to 5, as this represented a more realistic estimate of the average node degree in real networks (see Table 1). The mixing parameter, μ, i.e. the fraction of all links in a particular module that end outside this module, was varied from 0.1 to 0.5. Increasing the mixing parameter makes the modules of the 'true' community structure less well defined and the communities less easily detected. Testing for a mixing parameter greater than 0.5 was not deemed necessary, as it would contradict the definition of community structure, where more intra-community links than inter-community links should exist.
Overall, iMod performed better for all examples tested. For node degree equal to 16, iMod and the Greedy method manage to retrieve the exact partition for all values of μ up to 0.35. Thereafter, iMod outperforms the Greedy method by continuing to extract the exact partition whereas the Greedy method's performance declines rapidly. In the case of degree equal to 5, iMod still achieves higher similarity to the known structure than the Greedy method for all values of μ.
Detection of resolution limitations
Improved community structures are not achieved solely through maximisation of modularity; further refinement by addressing resolution limits of modularity maximisation is critically important. Network modules obtained with the iMod algorithm were further partitioned as described in the ResMod procedure, ignoring all inter-module links, as outlined above.
Computational results for modularity optimisation and resolution limits in simulated network examples.
Networks | iMod | ResMod | |||||
---|---|---|---|---|---|---|---|
Name | N | L | Median Q | Best Q | M | Q_Reso | M |
Ring | 30 | 40 | 0.6750 | 0.6750 | 5 | 0.6500 | 10 |
Y-shape | 50 | 40 | 0.5426 | 0.5426 | 3 | 0.5416 | 4 |
The first example is a ring-shaped network composed of 10 identical complete graphs of three nodes each, represented by circles inter-connected by the minimal number of links (Figure 4a). This graph is an example of maximal modularity, since modularity converges to one as the number of complete graphs reaches infinity [33, 56]. Modularity maximisation using iMod initially suggests the existence of 5 modules, in accordance to other approaches [33]. Through implementation of ResMod to correct for resolution limits by optimising each of the five communities further without considering the inter-module links, the two smaller groups within each module become apparent and the total number of modules is correctly identified as ten.
The second synthetic example comprises four groups of nodes (y-shaped, Figure 4b). Each group, denoted by a circle, consists of completely connected graphs: the two leftmost groups comprise 20 nodes and the two on the right consist of 5 nodes each [33]. Methods that perform modularity maximisation tend to merge the two smallest groups to yield the highest possible modularity value at the cost of an inaccurate detection of underlying community structure. Partitioning the network through iMod and optimising each module through ResMod yields the accurate number of four modules and the network is partitioned correctly.
Computational results for module detection without correction for resolution problems (iMod) and after accounting for resolution (ResMod).
Network | iMod | ResMod | ||
---|---|---|---|---|
Name | Q | M | Q | M |
Zachary | 0.420 | 4 | 0.420 | 4 |
Dolphin | 0.529 | 5 | 0.504 | 7 |
Les Miserables | 0.560 | 6 | 0.560 | 6 |
P53 | 0.535 | 7 | 0.469 | 16 |
Jazz | 0.445 | 4 | 0.445 | 4 |
E. coli | 0.781 | 19 | 0.675 | 79 |
S. cerevisiae | 0.775 | 25 | 0.693 | 66 |
C. elegans | 0.453 | 9 | 0.366 | 44 |
0.580 | 9 | 0.432 | 72 |
Another computational methodology that accounts for resolution limitations is the use of simulated annealing (SA) for modularity maximisation, where simulated annealing is applied to each detected module to find out whether any sub-modules can be identified [33]. In comparison, the E. coli network was partitioned into 79 modules with a modularity of 0.675 with ResMod, compared to 76 modules with modularity of 0.661 in SA. For the yeast network, ResMod achieves 66 modules with modularity of 0.693, as opposed to 57 modules with a total community modularity of 0.677 in SA. In both cases, ResMod succeeded in further dividing a higher number of modules while still achieving partitions with better overall modularity scores. Furthermore, another possible advantage of the methodology presented here is an explicit account to avoid over-partitioning, through implementation of the threshold value. However, further work is planned in the future to address: (i) quality control measures in module discovery to assess whether the detected community structure reflects phenotypic properties well, and (ii) further development of measures to avoid over-partitioning when accounting for resolution problems.
Conclusions
Community structure identification through modularity maximisation is hindered by (i) the NP-hard properties of the related optimisation problem and (ii) the resolution limitations introduced through the modularity measure. We have previously reported the detection of community structure in small to medium networks through a mixed integer quadratic programming procedure that guarantees global optimal solutions for modularity maximisation [30]. Here, we extend this work to tackle large size networks through an iterative optimisation procedure that performs well as evidenced through comparative analyses. As a further improvement, we also report methodological details of identifying and addressing resolution limitations, thus retrieving a more accurate representation of community structure from data.
Despite significant advances in the area of module detection through modularity optimisation, it is important to mention some caveats. First, modularity may not be the most appropriate measure of topological network features, as it can introduce limitations in practical applications [33, 57]. Alternative measures have been proposed [36] and will be studied in future work in terms of their ability to enhance module detection. It is worth noting that solution procedures presented here are generic and can be implemented with any mathematical expression of community presence other than modularity.
Furthermore, the use of coarse-grained topological features as a means to represent a complex network may not always be sufficient in delineating the intricate relationships and phenotypic properties of the system at hand. For example, in biological networks modularity is a phenomenon linked to a varying contribution of evolutionary inheritance of features, genome organisation properties and functional attributes [8]. Enriched network abstractions (e.g. edge weight and directionality), development of more accurate fitness functions (for example to capture cooperation effects [58]), as well as methodologies incorporating dynamic features can all contribute to future advances.
Network theory and related computational approaches have significantly enhanced our ability to offer deep insights into the principles governing complex systems. Analysis of protein interactions has shed light into mechanisms of disease [59–61], the association of genetic to phenotypic properties [7, 62] and biological species [2]. In this respect, the role of an accurate computational procedure to reveal the relations between the structure and functions in complex systems is important. Methodologies that allow communities to be detected both optimally and unambiguously, such as the ones presented in this paper, can greatly assist in this direction.
Declarations
Acknowledgements
The authors thank Mark Newman (Zachary, dolphin, Les Miserables, Metabolic and Email datasets, http://www-personal.umich.edu/~mejn/netdata/) and Uri Alon (E. coli and S. cerevisiae datasets, http://www.weizmann.ac.il/mcb/UriAlon) for providing datasets and Prof Christos Ouzounis for comments. GX acknowledges financial support from ORSAS (Overseas Research Students Awards Scheme) and the Centre for Process Systems Engineering. LB thanks the School of Physical Sciences and Engineering and the Systems Biomedicine Graduate program for financial support.
Authors’ Affiliations
References
- Barabasi AL, Albert R: Emergence of scaling in random networks. Science. 1999, 286: 509-512. 10.1126/science.286.5439.509PubMedView ArticleGoogle Scholar
- Spirin V, Gelfand MS, Mironov AA, Mirny LA: A metabolic network in the evolutionary context: multiscale structure and modularity. Proc Natl Acad Sci USA. 2006, 103: 8774-8779. 10.1073/pnas.0510258103PubMedPubMed CentralView ArticleGoogle Scholar
- Guimera R, Amaral LAN: Functional cartography of complex metabolic networks. Nature. 2005, 433: 895-900. 10.1038/nature03288PubMedPubMed CentralView ArticleGoogle Scholar
- Sales-Pardo M, Guimera R, Moreira AA, Amaral LA: Extracting the hierarchical organization of complex systems. Proc Natl Acad Sci USA. 2007, 104: 15224-15229. 10.1073/pnas.0703740104PubMedPubMed CentralView ArticleGoogle Scholar
- Milo R, Itzkovitz S, Kashtan N, Levitt R, Shen-Orr S, Ayzenshtat I, Sheffer M, Alon U: Superfamilies of evolved and designed networks. Science. 2004, 303: 1538-1542. 10.1126/science.1089167PubMedView ArticleGoogle Scholar
- Grunwald S, Speer A, Ackermann J, Koch I: Petri net modelling of gene regulation of the Duchenne muscular dystrophy. Biosystems. 2008, 92: 189-205. 10.1016/j.biosystems.2008.02.005PubMedView ArticleGoogle 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.2136809100PubMedPubMed CentralView 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.1073374PubMedView ArticleGoogle Scholar
- Newman MEJ: Modularity and community structure in networks. Proc Natl Acad Sci USA. 2006, 103: 8577-8582. 10.1073/pnas.0601602103PubMedPubMed CentralView ArticleGoogle Scholar
- Tamames J, Moya A, Valencia A: Modular organization in the reductive evolution of protein-protein interaction networks. Genome Biol. 2007, 8: R94- 10.1186/gb-2007-8-5-r94PubMedPubMed CentralView ArticleGoogle Scholar
- Dartnell L, Simeonidis E, Hubank M, Tsoka S, Bogle ID, Papageorgiou LG: Robustness of the p53 network and biological hackers. FEBS Lett. 2005, 579: 3037-3042. 10.1016/j.febslet.2005.03.101PubMedView ArticleGoogle Scholar
- Albert R, Jeong H, Barabasi AL: Error and attack tolerance of complex networks. Nature. 2000, 406: 378-382. 10.1038/35019019PubMedView ArticleGoogle Scholar
- Koonin EV, Wolf YI, Karev GP: The structure of the protein universe and genome evolution. Nature. 2002, 420: 218-223. 10.1038/nature01256PubMedView ArticleGoogle Scholar
- Kunin V, Pereira-Leal JB, Ouzounis CA: Functional Evolution of the Yeast Protein Interaction Network. Molecular Biology and Evolution. 2004, 21: 1171-1176. 10.1093/molbev/msh085PubMedView ArticleGoogle Scholar
- Newman MEJ: The structure of scientific collaboration networks. Proc Natl Acad Sci USA. 2001, 98: 404-409. 10.1073/pnas.021544898PubMedPubMed CentralView ArticleGoogle Scholar
- Holme P, Huss M, Jeong H: Subnetwork hierarchies of biochemical pathways. Bioinformatics. 2003, 19: 532-538. 10.1093/bioinformatics/btg033PubMedView ArticleGoogle Scholar
- Fortunato S, Castellano C: Community structure in graphs. Encyclopedia of Complexity and System Science. 2008, arXiv:0712.2716,Google Scholar
- Fortunato S: Community detection in graphs. Physics Reports. 2010, 486: 75-174. 10.1016/j.physrep.2009.11.002View ArticleGoogle Scholar
- Garey MR, Johnson DS: Computers and Intractability: A Guide to the Theory of NP-Completeness. 1979, San Francisco: W.H.Freeman & Co Ltd,Google Scholar
- Boccaletti S, Latora V, Moreno Y, Chavez M, Hwang DU: Complex networks: Structure and dynamics. Physics Reports. 2006, 424: 175-308. 10.1016/j.physrep.2005.10.009View ArticleGoogle Scholar
- Rives AW, Galitski T: Modular organization of cellular networks. Proc Natl Acad Sci USA. 2003, 100: 1128-1133. 10.1073/pnas.0237338100PubMedPubMed CentralView ArticleGoogle Scholar
- Gustafsson M, Hornquist M, Lombardi A: Comparison and validation of community structures in complex networks. Physica A. 2006, 367: 559-576. 10.1016/j.physa.2005.12.017View ArticleGoogle Scholar
- Grafahrend-Belau E, Schreiber F, Heiner M, Sackmann A, Junker BH, Grunwald S, Speer A, Winder K, Koch I: Modularization of biochemical networks based on classification of Petri net t-invariants. BMC Bioinformatics. 2008, 9: 90- 10.1186/1471-2105-9-90PubMedPubMed CentralView ArticleGoogle Scholar
- Palla G, Derenyi I, Farkas I, Vicsek T: Uncovering the overlapping community structure of complex networks in nature and society. Nature. 2005, 435: 814-818. 10.1038/nature03607PubMedView ArticleGoogle Scholar
- Newman MEJ, Girvan M: Finding and evaluating community structure in networks. Phys Rev E. 2004, 69: 026113-10.1103/PhysRevE.69.026113. 10.1103/PhysRevE.69.026113View ArticleGoogle Scholar
- Newman MEJ: Fast algorithm for detecting community structure in networks. Phys Rev E. 2004, 69: 066133-10.1103/PhysRevE.69.066133. 10.1103/PhysRevE.69.066133View ArticleGoogle Scholar
- Girvan M, Newman MEJ: Community structure in social and biological networks. Proc Natl Acad Sci USA. 2002, 99: 7821-7826. 10.1073/pnas.122653799PubMedPubMed CentralView ArticleGoogle Scholar
- Brandes U, Delling D, Gaertler M, Gorke R, Hoefer M, Nikoloski Z: On Modularity Clustering. IEEE Transactions on Knowledge and Data Engineering. 2008, 20: 172-188. 10.1109/TKDE.2007.190689View ArticleGoogle Scholar
- Cafieri S, Hansen P, Liberti L: Edge ratio and community structure in networks. Phys Rev E. 2010, 81: 026105-10.1103/PhysRevE.81.026105. 10.1103/PhysRevE.81.026105View ArticleGoogle Scholar
- Xu G, Tsoka S, Papageorgiou LG: Finding community structures in complex networks using mixed integer optimisation. Eur Phys J B. 2007, 60: 231-239. 10.1140/epjb/e2007-00331-0View ArticleGoogle Scholar
- Medus A, Acuna G, Dorso CO: Detection of community structures in networks via global optimization. Physica A. 2005, 358: 593-604. 10.1016/j.physa.2005.04.022View ArticleGoogle Scholar
- Duch J, Arenas A: Community detection in complex networks using extremal optimization. Phys Rev E. 2005, 72: 027104-027200. 10.1103/PhysRevE.72.027104View ArticleGoogle Scholar
- Fortunato S, Barthelemy M: Resolution limit in community detection. Proc Natl Acad Sci USA. 2007, 104: 36-41. 10.1073/pnas.0605965104PubMedPubMed CentralView ArticleGoogle Scholar
- Ruan J, Zhang W: Identifying network communities with a high resolution. Phys Rev E. 2008, 77: 016104-10.1103/PhysRevE.77.016104. 10.1103/PhysRevE.77.016104View ArticleGoogle Scholar
- Arenas A, Fernandez A, Gomez S: Analysis of the structure of complex networks at different resolution levels. New Journal of Physics. 2008, 10: 053039-10.1088/1367-2630/10/5/053039. 10.1088/1367-2630/10/5/053039View ArticleGoogle Scholar
- Li Z, Zhang S, Wang RS, Zhang XS, Chen L: Quantitative function for community detection. Phys Rev E. 2008, 77: 036109-10.1103/PhysRevE.77.036109. 10.1103/PhysRevE.77.036109View ArticleGoogle Scholar
- Reichardt J, Bornholdt S: Detecting Fuzzy Community Structures in Complex Networks with a Potts Model. Phys Rev Lett. 2004, 93: 218701- 10.1103/PhysRevLett.93.218701PubMedView ArticleGoogle Scholar
- Kumpula JM, Saramaki J, Kaski K, Kertesz J: Limited resolution in complex network community detection with Potts model approach. Eur Phys J B. 2007, 56: 41-45. 10.1140/epjb/e2007-00088-4View ArticleGoogle Scholar
- Brooke A, Kendrick D, Meeraus A, Raman R: GAMS: A user's guide. 2003, Washington, DC: GAMS Development Corp,Google Scholar
- The GAMS SBB solver. http://www.gams.com/dd/docs/solvers/sbb.pdf
- , : ILOG CPLEX10.0 User's Manual. 2006,Google Scholar
- Zachary WW: An information flow model for conflict and fission in small groups. J Anthropol Res. 1977, 33: 452-473.Google Scholar
- Lusseau D: The emergent properties of a dolphin social network. Biology Letters. 2003, 270: S186-S188.Google Scholar
- Lusseau D, Schneider K, Boisseau OJ, Haase P, Slooten E, Dawson SM: The bottlenose dolphin community of Doubtful Sound features a large proportion of long-lasting associations. Can geographic isolation explain this unique trait?. Behavioral Ecology and Sociobiology. 2003, 54: 396-405. 10.1007/s00265-003-0651-yView ArticleGoogle Scholar
- Knuth DE: The Stanford graphbase: a platform for combinatorial computing. 1993, Reading, MA: Addison-Wesley,Google Scholar
- Gleiser PM, Danon L: Community Structure in Jazz. Advances in Complex Systems. 2003, 6: 565-574. 10.1142/S0219525903001067View ArticleGoogle Scholar
- Guimera R, Danon L, Diaz-Guilera A, Giralt F, Arenas A: Self-similar community structure in a network of human interactions. Phys Rev E. 2003, 68: 065103-10.1103/PhysRevE.68.065103. 10.1103/PhysRevE.68.065103View ArticleGoogle Scholar
- Shen-Orr SS, Milo R, Mangan S, Alon U: Network motifs in the transcriptional regulation network of Escherichia coli. Nat Genet. 2002, 31: 64-68. 10.1038/ng881PubMedView ArticleGoogle Scholar
- Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks. Science. 2002, 298: 824-827. 10.1126/science.298.5594.824PubMedView ArticleGoogle Scholar
- Jeong H, Tombor B, Albert R, Oltvai ZN, Barabasi AL: The large-scale organization of metabolic networks. Nature. 2000, 407: 651-654. 10.1038/35036627PubMedView ArticleGoogle Scholar
- Xiang J, Hu K, Tang Y: A class of improved algorithms for detecting communities in complex networks. Physica A. 2008, 387: 3327-3334. 10.1016/j.physa.2008.01.105View ArticleGoogle Scholar
- Blondel VD, Guillaume JL, Lambiotte R, Lefebvre E: Fast unfolding of communities in large networks. J Stat Mech. 2008, P10008-10.1088/1742-5468/2008/10/P10008.Google Scholar
- Li C, Li X, Miao Y, Wang Q, Jiang W, Xu C, Li J, Han J, Zhang F, Gong B, Xu L: SubpathwayMiner: a software package for flexible identification of pathways. Nucleic Acids Res. 2009, 37: e131- 10.1093/nar/gkp667PubMedPubMed CentralView ArticleGoogle Scholar
- Lancichinetti A, Fortunato S, Radicchi F: Benchmark graphs for testing community detection algorithms. Phys Rev E. 2008, 78: 046110-10.1103/PhysRevE.78.046110. 10.1103/PhysRevE.78.046110View ArticleGoogle Scholar
- Lancichinetti A, Fortunato S, Kertesz J: Detecting the overlapping and hierarchical community structure in complex networks. New Journal of Physics. 2009, 11: 033015-10.1088/1367-2630/11/3/033015. 10.1088/1367-2630/11/3/033015View ArticleGoogle Scholar
- Danon L, Diaz-Guilera A, Duch J, Arenas A: Comparing community structure identification. J Stat Mech. 2005, P09008-10.1088/1742-5468/2005/09/P09008.Google Scholar
- Rosvall M, Bergstrom CT: An information-theoretic framework for resolving community structure in complex networks. Proc Natl Acad Sci USA. 2007, 104: 7327-7331. 10.1073/pnas.0611034104PubMedPubMed CentralView ArticleGoogle Scholar
- Lozano S, Arenas A, Sanchez A: Mesoscopic structure conditions the emergence of cooperation on social networks. PLoS One. 2008, 3: e1892- 10.1371/journal.pone.0001892PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor IW, Linding R, Warde-Farley D, Liu Y, Pesquita C, Faria D, Bull S, Pawson T, Morris Q, Wrana JL: Dynamic modularity in protein interaction networks predicts breast cancer outcome. Nat Biotechnol. 2009, 27: 199-204. 10.1038/nbt.1522PubMedView ArticleGoogle Scholar
- Jiang X, Liu B, Jiang J, Zhao H, Fan M, Zhang J, Fan Z, Jiang T: Modularity in the genetic disease-phenotype network. FEBS Lett. 2008, 582: 2549-2554. 10.1016/j.febslet.2008.06.023PubMedView ArticleGoogle Scholar
- Lee DS, Park J, Kay KA, Christakis NA, Oltvai ZN, Barabasi AL: The implications of human metabolic network topology for disease comorbidity. Proc Natl Acad Sci USA. 2008, 105: 9880-9885. 10.1073/pnas.0802208105PubMedPubMed CentralView ArticleGoogle Scholar
- Slonim N, Elemento O, Tavazoie S: Ab initio genotype-phenotype association reveals intrinsic modularity in genetic networks. Mol Syst Biol. 2006, 2: 2006 0005,Google Scholar
Copyright
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.