Investigating the validity of current network analysis on static conglomerate networks by protein network stratification
© Zhang and Lu; licensee BioMed Central Ltd. 2010
Received: 2 February 2010
Accepted: 16 September 2010
Published: 16 September 2010
A molecular network perspective forms the foundation of systems biology. A common practice in analyzing protein-protein interaction (PPI) networks is to perform network analysis on a conglomerate network that is an assembly of all available binary interactions in a given organism from diverse data sources. Recent studies on network dynamics suggested that this approach might have ignored the dynamic nature of context-dependent molecular systems.
In this study, we employed a network stratification strategy to investigate the validity of the current network analysis on conglomerate PPI networks. Using the genome-scale tissue- and condition-specific proteomics data in Arabidopsis thaliana, we present here the first systematic investigation into this question. We stratified a conglomerate A. thaliana PPI network into three levels of context-dependent subnetworks. We then focused on three types of most commonly conducted network analyses, i.e., topological, functional and modular analyses, and compared the results from these network analyses on the conglomerate network and five stratified context-dependent subnetworks corresponding to specific tissues.
We found that the results based on the conglomerate PPI network are often significantly different from those of context-dependent subnetworks corresponding to specific tissues or conditions. This conclusion depends neither on relatively arbitrary cutoffs (such as those defining network hubs or bottlenecks), nor on specific network clustering algorithms for module extraction, nor on the possible high false positive rates of binary interactions in PPI networks. We also found that our conclusions are likely to be valid in human PPI networks. Furthermore, network stratification may help resolve many controversies in current research of systems biology.
In the contemporary systems biology, the cell itself can be viewed as a complex network of interacting proteins, nucleic acids, and other biomolecules [1–6]. The network representation has been widely applied to describe various molecular systems including protein interaction maps, metabolites and reactions, transcriptional regulation maps, signal transduction pathways and functional association networks [3, 7–14]. Network approaches have been proven useful in predicting protein functions, guiding large-scale experiments, facilitating drug discovery and design, and expediting novel biomarker identification [15–21].
Due to the relative scarcity of the protein-protein interaction (PPI) data, a common practice is to assemble all available binary PPIs of a certain organism into a combined static conglomerate network in order to gain a systems-level view on this organism [8, 18, 22, 23]. Similar approaches are also applied to build other types of molecular networks, such as gene regulatory networks and metabolic networks [24, 25]. These combined conglomerate networks are often an integration of diverse datasets generated from high-throughput and small-scale experiments, predictive computational methods, and expert curations [6, 18, 26, 27]. Network analyses are then performed on these conglomerate networks. Many structural characteristics of molecular networks have been revealed [3, 24, 28–30].
However, since static conglomerate networks correspond to the collective of all available data components under various temporal and spatial conditions for certain organisms, studies using such networks often ignore the dynamic nature of molecular systems, and thus lose context-dependent information. Indeed, recent pioneering efforts in revealing network dynamics have suggested that the structures of static conglomerate networks might differ significantly from those of context-dependent networks under specific conditions. For example, Luscombe et al. examined the yeast transcriptional networks and discovered that topological characteristics, network motifs and crucial transcription factors as hub nodes (nodes that have high degree values), respectively, are different under five conditions (i.e., cell cycle, sporulation, diauxic shifts, DNA damage, stress response) . Using the same condition-specific networks, Zhang et al. further discovered that regulatory patterns for transcription factor hubs changed substantially under different contexts . Han et al. reported the status of hubs in yeast PPI network changes under different temporal conditions, thereby, can be classified into "party hubs" and "date hubs" . Overlaying time-specific gene expression data onto yeast cell cycle protein interactions, de Lichtenberg et al. reported dynamic variations of protein complexes during the yeast cell cycle based on a time-dependent interaction network . Recent studies have also shown that house-keeping genes and tissue-specific genes have different topological properties in the human PPI network [34, 35].
A recently published genome-scale tissue- and condition-specific proteomics data of Arabidopsis thaliana, as well as the availability of several large-scale A. thaliana PPI datasets has provided us with an excellent opportunity to carry out this investigation. In this study, we employed an approach called network stratification to extract context-dependent subnetworks from a large-scale conglomerate network by integrating genome-scale context-dependent data (see Methods). Particularly, conglomerate networks can be stratified in at least three dimensions: time (disparate time points or periods), space (diverse organisms, tissues, or subcellular localizations), and condition (cell cycle, stress response, DNA damage, etc.). A stratified context-dependent subnetwork, for example a PPI network containing genes that only exist in the mouse kidney or a human PPI network consisting of differentially expressed genes of primary breast cancer samples, should still encode systems-level information. Compared with unstratified conglomerate networks, stratified context-dependent subnetworks correspond to relatively homogenous contexts, and therefore network stratification may help reveal the dynamic nature of biological systems and underlying biological processes.
Using this network stratification strategy, we investigated three types of most commonly conducted network analyses: topological, functional and modular analyses . Topological analysis consists of 1) various network topology characteristics such as node degree and its distribution, clustering coefficient, eccentricity, and node betweenness; 2) compatibility, or agreement of nodes and interactions between network pairs; and 3) roles of important nodes such as network hubs and bottlenecks (nodes with high betweenness centrality values). Functional analysis includes the identification of specifically enriched functions in different networks as well as differences in functional enrichment among networks. Modular analysis considers the agreement among modules or clusters extracted from different networks and functional enrichment of members in modules. The goal of this study is to address the questions of whether the results based on conglomerate networks and context-dependent subnetworks are different, and to what extent the current practice can be used. In this study, the results are mainly based on an unstratified network as well as five stratified tissue-specific subnetworks.
The novelties of our study are: First, although researchers have started to study dynamic or tissue-specific networks, the above questions have not yet been systematically addressed. Here these questions will be addressed by performing systematic comparative analyses between the unstratified network and stratified subnetworks. Second, we use proteomics data instead of mRNA expression data for integration with protein interactions, while previous works on dynamic or tissue-specific networks all used gene expression data to build subnetworks. Although there is a general trend for protein concentration to vary with mRNA expression levels, observations suggest individual gene expression intensities may not be well correlated with corresponding protein abundance values [38, 39]. Therefore, it is more direct and precise to stratify protein networks with proteomics data than using gene expression data.
Network Assembly and Subnetwork Construction
We assembled a conglomerate PPI network of A. thaliana by combining predicted, verified and curated PPI datasets from TAIR and AtPID databases [40, 41] (see Methods). The assembled conglomerate A. thaliana PPI network consists of 13,136 proteins and 42,131 interactions. We then constructed 18 A. thaliana PPI networks with respect to specific or combinatorial temporal and spatial conditions, according to the total and 17 context-dependent protein lists of the plant described in  (see Methods). These 18 networks form a 4-level hierarchy [Supp. Figure 1 in Additional file 1]. The level-1 network is the part of the assembled conglomerate network that contains all proteins included in the proteomics data. It corresponds to the collective temporal dynamics and tissue specificity of both A. thaliana organs and cell cultures; henceforth referred to as the unstratified total network, or the total network. Level-2 organs and cell culture networks are coarse-stratified subnetworks representing A. thaliana PPIs in organs and cell culture, respectively. Level-3 tissue- and cell culture condition-specific networks are moderate-stratified subnetworks corresponding to comparatively more specific contexts. And level-4 sub-tissue-specific networks are fine-stratified subnetworks related to the most specific tissue-related contexts in this study.
Like many large-scale molecular networks, the 18 unstratified and stratified networks all display scale-free characteristics, i.e., power-law in node degree distribution [Table S1 in Additional file 2]. However, stratified subnetworks and the unstratified total network generally differ from one another in most network statistics including average node degree, eccentricity, and node betweenness [Table S2 in Additional file 2]. Generally, significant differences of average node degree values exist among 18 networks. Pairwised Welch's t-tests show that 112 out of 153 t-scores between any pair of 18 networks have p-values < 0.01, many of which are well below 1e-5 [Table S2 in Additional file 2]. In addition, average degrees of subnetworks are significantly different from those of randomly stratified subnetworks by node sampling that are scale-free and comparable in total nodes and interactions to the actual subnetworks (p-value < 0.01). The maximum node degree of different networks has a substantial difference intuitively, due to the fact that the unstratified total network is a union of stratified subnetworks. For instance, the protein with the maximum degree value in the unstratified total PPI network (SUMO1, or AT4G26840) has 146 interacting partners, while the highest node degree in fine-stratified cotyledons and moderate-stratified seeds PPI networks are only 41 (PBD1, or AT3G22630) and 54 (IMPA-6, or AT1G02690), respectively. Similarly, network statistics such as eccentricity and node betweenness among networks all show measurable differences with statistical significance (p-value < 0.01) [Table S2 in Additional file 2]. Average clustering coefficient values are relatively similar between the total network and stratified subnetworks, all being significantly higher than random expectations by node sampling (p-value < 0.01), which indicates the universal existence of modularity in all these networks. Overall, the results suggest that the majority of network statistics calculated from conglomerate networks are not transferable when evaluating context-dependent subnetworks.
Topology Analysis: Hubs and Bottlenecks May Change Status during Stratification
Identification of important nodes, for example, network hubs (i.e., nodes with high degree) and bottlenecks (i.e., nodes with high betweenness centrality) , are commonly conducted in molecular network topology analysis. If stratified PPI subnetworks differ significantly in network topology from the unstratified total network, we should expect status changes as hubs or bottlenecks for the same proteins in different networks during stratification. Such status changes may reveal dynamics of networks under specific temporal, spatial or environmental conditions. Here we focus on these two types of important nodes: network hubs and bottlenecks.
The numbers/percentages of protein hub status change, for both "hub to non-hub" and "non-hub to hub", are generally significantly different from those of randomized networks of similar sizes (p-value < 0.05) [Table S4A in Additional file 2]. A smaller number of hub status change compared with randomization indicates that in PPI networks, the roles of hubs and non-hubs are more likely to be preserved than those in randomized networks. This is possibly due to the evolved organizational structure in PPI networks that facilitates the preservation of hub proteins among organisms or tissues, and such structure does not exist in random networks . Results showed consistency when using different percentage cutoffs (5% or 10%) for hub definition [Table S4A and S5 in Additional file 2].
Scenario 1: Hubs in the total network remain hubs in stratified subnetworks
Conserved in all five level-3 stratified tissue-specific subnetworks, SAC52 (AT1G14320, suppressor of acaulis 52) is a big hub with 25 interacting partners in the unstratified total network, which remains a hub in each of the tissue-specific subnetworks [Figure 4A]. The biological process function of SAC52, translation, is enriched in each of the five tissues [36, 42], which explains the central role of this gene in any specific tissue. Similarly, many proteins with tissue-independently enriched functions (such as translation, glycolysis), for example GAPCP-1 (AT1G79530) and RPS8A (AT5G20290), fall in this category. Interact with many partners under various contexts, these pan-tissue hubs may include the previously defined "party hubs" that have simultaneous interactions, as well as a portion of "date hubs" that are still "date hubs" in subnetworks .
Scenario 2: Hubs in the total network become non-hubs in subnetworks
This scenario includes three sub-categories: tissue-preferred hubs, tissue-specific hubs and tissue non-hubs. Tissue-preferred hubs are those that remain hubs in at least one tissue, while losing their hub status in other tissues. For example, AP19 (protein binding/protein transporter; AT2G17380) is a hub in the total network (degree 11) [Figure 4B]. Although AP19 is still a hub in the stratified roots subnetwork (degree 10), it becomes a non-hub in stratified leaves, flowers, siliques, and seeds subnetworks (with degree of 6, 8, 6, 5, respectively). The main biological process functions of AP19 are intracellular protein transport (GO:0006886) and protein transport (GO:0015031), the latter being significantly enriched in the root of the plant [36, 42]. This may explain why it becomes less important in other tissues of the plant. Proteins in this category usually have functions enriched in certain specific tissues, for example chloroplast- and cold response-related protein CPN60B (AT1G55490; hubs in the total network and stratified leaves subnetwork), thereby becoming less important and losing their hub status in other tissues.
Tissue-specific hubs are the hubs in the total network that only exist in one specific tissue and remain hubs in the tissue. VIP4 (AT5G61150; Vernalization independence 4) is a hub in the unstratified total network [Figure 4C]. When stratifying the total network, VIP4 disappears in any other tissue-specific subnetworks, yet remains a hub in stratified flowers subnetwork. The functions of VIP4 include negative regulation of flower development (GO:0009910) and vernalization response (GO:0010048) [40, 43], which may explain its particular role in flowers of the plant. Although VIP4 only exists in flowers, because it is highly connected, it becomes a hub in the total network.
Tissue non-hubs lose their hub status in all tissue-specific subnetworks after stratifying the total network. Such proteins may interact with different partners and perform different functions under different contexts. For instance, RPS15 (AT1G04270; Cytosolic ribosomal protein S15) is a hub in the total network (degree 10), which becomes a non-hub in any tissue-specific subnetwork [Figure 4D]. Though not studied in this work, we also expect hubs with transient interactions (i.e., dynamic interactions that are time-dependent) to fall in this category. Because a conglomerate network combines all context-dependent interactions, the corresponding protein may appear a hub in the conglomerate network.
In this scenario, the hubs identified based on a conglomerate network are incorrect in at least one context-dependent network.
Scenario 3: Non-hubs in the total network are hubs in subnetworks
Appearing as a non-hub in the total network (degree 9, top 25% but not top 20%), AK-HSDH_I (AT1G31230; Aspartate kinase-homoserine dehydrogenase I) is actually a hub in the stratified leaves subnetwork (degree 8, top 20%). This non-hub to hub change, albeit not substantial in terms of degree percentage ranking, may have biological significance because the majority of its interacting partners preserve in the smaller leaves subnetwork [Figure 4E]. In fact, AK-HSDH_I have functions chloroplast (GO:0009507) and chloroplast stroma (GO:0009570) of GO biological processes, which are both leaf-related functions [44, 45]. This type of hubs in a tissue-specific subnetwork is largely equivalent to the notion of "local" hubs in .
In order to quantitatively measure the variation of node degree, or node "hubbiness" between the unstratified total network and stratified subnetworks, we further classified nodes into five exclusive classes according to their degree values. In each network, using the top 5, 10, 20, and 50 percent of degree values as cutoffs, nodes/proteins are assigned into one of the five exclusive classes of different "hubbiness". After assigning each protein in all 18 different networks of four levels into one of the five classes, a substantial portion of proteins, 38.8% (2,562 out of 6,606) change in degree class among all networks, i.e., "hubbiness" of these nodes changes in one of the stratified subnetworks compared with that of the total network. If considering a change of two or more in degree class for the same nodes in two different networks as a "leap" change, using the same classification, 5.4% (354 out of 6,606) proteins have leap changes in their degree classes [Table S6A in Additional file 2]. We applied different cutoffs and different numbers of classes for classification, and the results showed consistency [Table S7 in Additional file 2].
While hubs as important nodes denote high local connectivity in networks , bottlenecks are nodes with high betweenness centrality that may represent connectors between clusters, or communities in networks. Among the 18 networks in the 4-level hierarchy, the overlap between network hubs and bottlenecks ranges from 0.39 to 0.55 measured by the Jaccard Index. The rather high overlap between hubs and bottlenecks is consistent with a previous study that degree and betweenness are highly correlated in PPI networks . We adopted the same approach as that of hub status change measurement to measure bottleneck status change between the total network and stratified tissue-specific subnetworks. The results are similar to those of hub status change except that in contrast to the minor "non-hub to hub" changes, a higher number/percentage (9% on average, which is calculated by dividing the number of non-bottlenecks that change into bottlenecks in a subnetwork by the total number of bottlenecks in the same subnetwork) of "non-bottleneck to bottleneck" changes is observable [Figure 3, Table S4B and Table S5 in Additional file 2], which indicates that modules in the total network as densely connected subnetworks are relatively sensitive to the change of network topology. These results suggest that the majority of bottlenecks and non-bottlenecks in conglomerate networks remain unchanged during stratification, though bottleneck status change appears in a higher degree than hubs. Therefore, the majority of bottlenecks and non-bottlenecks identified based on conglomerate networks are still valid in context-dependent networks.
The Total Network and Stratified Subnetworks Differ in Functional Enrichment
Because a protein may have multiple GO annotations but only perform some of its annotated functions in a specific subnetwork, it would be more precise to consider interacting pairs in a network that share the same GO terms to determine their true functions in the network. Therefore, we selected interacting protein pairs that share the same biological process functions and evaluated the distribution of such pairs in each network into the same pool of biological process annotations. Fisher's exact tests identified a larger portion of highly enriched functions with generally more significant p-values than the results from single-node/protein function analysis in each of the stratified subnetworks compared with interacting protein pairs in the total network [Supp. Figure 2 in Additional file 1]. While photosynthesis in leaves, protein transport in roots, and response to heat in seeds are still highly enriched, protein folding, protein catabolic process, as well as ubiquitin-dependent protein catabolic process become universally and significantly enriched in all five tissues of the plant. The enriched functions in each of the tissue-specific subnetworks are consistent with previous studies .
We further investigated general differences in protein functional enrichment among the total network and stratified subnetworks, using GO slim high-level annotations by TAIR and TIGR [42, 47] [Table S8 in Additional file 2]. If proteins in the total network differ in functional enrichment from those of level-2 coarse-stratified organs and cell culture subnetworks and the five level-3 stratified tissue-specific subnetworks, we should expect a difference in among-functional category distributions between GO functional annotations from each of the stratified subnetworks and those from the total network. χ2 tests showed such among-biological processes GO term distributions are significantly different (p-value < 1e-5) between protein annotation counts from the unstratified total network and each of the stratified subnetworks [Figure 6]. Using interacting pair annotation counts, χ2 tests identified similar yet more statistically significant results for comparing the distributions of protein interacting pair annotations among networks [Supp. Figure 2 in Additional file 1]. In addition, a different among-network distribution between proteins of a specific GO term and the total number of annotated proteins is observable. For instance, protein annotation counts in all 15 GO slim terms of Molecular Function category form a distribution among 18 networks. Annotation counts in GO slim term nucleotide binding of Molecular Function form another distribution among networks, which is significantly different from the distribution of all Molecular Function protein annotations, with a p-value of 2.0e-8 by a χ2 test. Similar tests show that the majority of GO terms in biological process, molecular function, and cellular component all have those significantly different distributions [Figure 6]. Using interacting pair annotation counts, instead of single protein annotation counts, an even larger portion of GO terms shows significantly different distributions [Supp. Figure 2 in Additional file 1].
These results suggest that over-represented protein functions in the total network are often different from those in subnetworks; therefore, the enriched functions of a conglomerate network do not reflect the true functions of context-dependent networks.
The Total Network and Stratified Subnetworks Differ in Modular Structures
Modules, or communities in molecular networks, are expected to correspond to functional units [1, 48–53], which are commonly analyzed. To extract modules from networks, we used three different clustering algorithms to partition each of the 18 networks in the 4-level hierarchy: a simulated annealing-based algorithm to optimize a defined modularity score ; a density-based network clustering algorithm , which identifies clique-like components (densely connected subnetworks) in a network as modules; an edge betweenness-based partitioning method  partitioning a network into modular structures by iteratively removing interactions of the highest betweenness scores.
The Total Network and Stratified Subnetworks Differ in Modular Functional Enrichment
We evaluated the GO function distribution of single proteins as well as protein interacting pairs in modules extracted from each network, based on the same biological process annotation terms as those used for functional enrichment analysis in stratified subnetworks. Using modules extracted by the simulated-annealing based method as an example, Fisher's exact tests show that, in general, proteins in the majority of modules clustered from the total network and the five level-3 stratified tissue-specific subnetworks are enriched in at least one GO biological process term. The percentages of modules enriched for the total network, and roots, leaves, flowers, siliques and seeds subnetworks are 83.6% (51 out of 61 modules), 90.9% (30 out of 33), 100% (33 out of 33), 85.1% (40 out of 47), 86.4% (38 out of 44), and 96.6% (28 out of 29), respectively. Considering the heterogeneity of components in modules extracted from the total network, i.e., proteins and interactions in the total network are combined from those under different contexts, it is reasonable that the percentage of functionally enriched modules is the lowest for the total network. Twenty-six biological process functional terms are universally enriched by proteins in modules from all six networks (the total network and five tissue-specific subnetworks), such as glycolysis, translation, transport, and protein folding; whereas, modules extracted from each of the five tissue-specific subnetworks have a number of exclusively enriched GO terms, representing tissue-specific functions, for example proton transport in roots, photosynthetic electron transport chain in leaves, negative regulation of flower development in flowers, phospholipid biosynthetic process in siliques, and anaerobic respiration in seeds [Table S10A in Additional file 2 and 3]. In addition, a number of modules from the total network are observed to have enriched GO terms that are not over-represented in any module from any tissue-specific subnetwork [Table S10A in Additional file 2 and 3]. This suggests that such modules and corresponding enriched functions may be artifacts due to the assembly of conglomerate networks by combining PPIs from various contexts, as illustrated in [Figure 1]. The results are consistent when using protein interacting pair annotations for modular functional enrichment [Table S10B in Additional file 2 and 3]. Similar results were obtained when evaluating modular functional enrichment by using GO slim high-level terms [Table S8 in Additional file 2] of all molecular function, biological process and cellular component categories [Table S11 in Additional file 2].
In summary, with respect to modular structures and modular functional enrichment, results of modular analysis based on conglomerate networks are mostly not transferable to context-dependent networks.
Effects of False Interactions
It is widely suspected that the PPIs determined from high-throughput experiments may contain a large number of false interactions, or false positives [54, 55]. In a recent publication, Huang et al. estimated the false discovery rates to be up to 20% for large-scale PPI screens . In another study, Wuchty also used 20% as the false interaction rate in his simulation . Therefore, we expect a similar percentage of false PPIs to exist in our assembled conglomerate network. To address the issue of the possible influence of these false interactions in the assembled conglomerate network on the result of our analyses, we randomly replaced 10 to 20% of PPIs from the total network and performed network stratification and the same comparative analyses on networks with replaced PPIs. The results are generally consistent with those from network analysis without replaced PPIs. For example, the percentages of network hub status change with 10% replaced interactions are highly similar to those of original networks, for both "hubs to non-hubs" and "non-hubs to hubs" changes (see Additional file 4). When replaced interactions were increased to 20%, a slight increase about 3% in "non-hubs to hubs" change and a slight decrease about 5% in "hubs to non-hubs" change are observable [Supp. Figure 3 in Additional file 1]. When replaced interactions were further increased to 30 or 40%, although characteristics of modified networks become deviated from the original ones, for example the clustering coefficient drops significantly from ~0.3 to 0.06 when 40% interactions were replaced in the modified networks, the analysis results are still generally consistent while being different from random expectations (Additional file 4). These results indicate the robustness of our analyses against false PPIs.
On the other hand, currently identified PPIs may be only a portion of the entire interactome with a large number of PPIs yet to be discovered. The effects of these unidentified PPIs, or false negatives, are difficult to assess directly by simulation. In this study, in order to limit the effects of false negatives, we have incorporated as many PPIs as available into the conglomerate PPI network (including validated, predicted, and curated interactions) to achieve high data coverage. We expect that the network structure and properties of the conglomerate network provide good approximations for those of the entire interactome. Therefore, the observations made in this study will likely be valid for the entire interactome. However, this can only be verified when more true PPIs become available in the future.
Stratification Using Human PPI Networks
We also investigated whether the conclusions obtained from A. thaliana PPI networks can be applicable in other organisms. Therefore, we also conducted network stratification and analyses on human PPI networks. Due to the incompleteness and heterogeneity of context-specific proteomics data in human, we limited our preliminary study in human to stratifying conglomerate human PPI networks into two subnetworks of human brain and kidney. We applied network stratification using human brain and kidney proteomics data, and comparatively analyzed various aspects among stratified brain- and kidney-specific subnetworks and static conglomerate human PPI networks. Similar approaches and methods were adopted as those applied in A. thaliana network analyses, and results showed that most of our conclusions from A. thaliana network analyses are also valid in human PPI networks. For example, network hub and bottleneck status changes during stratification in human networks are similar to those in A. thaliana networks [Supp. Figure 4 in Additional file 1]. Generally, stratified brain and kidney subnetworks are different from conglomerate networks in network statistical characteristics, topological properties, including hub/bottleneck status changes, and modular structures. Detailed data can be found in Additional file 5.
The main reason we postulate that conglomerate networks do not reflect the true structural and functional properties of context-dependent networks is that context-dependent and meaningful biological information may be lost or misinterpreted in conglomerate networks. As implied by our results, when analyzing molecular networks, for example PPI networks, for studies focusing on specific conditions or tissues, using conglomerate networks for analysis may create errors. Ideally, rather than analyzing a conglomerate network, using stratified context-specific networks, meaningful information can be truly retrieved. On the other hand, over-stratification, i.e., network stratification based on a limited amount of highly context-specific datasets, should be avoided, because the incompleteness of context-specific information may lead to additional errors in the system-level view of the corresponding dynamics.
Summary: the validity of current network analysis on static conglomerate networks
How different are the results?
Network statistics except average clustering coefficient are generally significantly different between the total network and stratified subnetworks
To some degree
Functional Enrichment in Network
Proteins in different networks are enriched in context-specific functions
Hubs in Conglomerate Networks
To some degree
16% hubs in the total network change into non-hubs in subnetworks on average
Non-hubs in Conglomerate Networks
0.8% hubs in subnetworks are non-hubs in the total network on average, which correspond to a trivial portion (< 0.1%) of non-hubs in the total network
Bottlenecks in Conglomerate Networks
To some degree
13% bottlenecks in the total network change into non-bottlenecks in subnetworks on average
Non-bottlenecks in Conglomerate Networks
9% bottlenecks in subnetworks are non-bottlenecks in the total network on average, which correspond to a trivial portion of non-bottlenecks in the total network
Modules from the total network and those from subnetworks have low Modular Compatibility
Functional Enrichment in Modules
Modules from different networks differ in over-represented context-specific functions
Network stratification may help resolve many controversies in current network biology. For instance, as implied from our results, unreasonably big hubs (nodes that have a huge number of interacting partners in a network) in conglomerate networks  may be a result of combining interactions under different conditions, under which these hubs may have only a moderately large number of interactions. Another example is that, contrary to the common belief that network modules should correspond to functional units, structural modules extracted from yeast PPI networks showed lack of functional enrichment . This controversy may be resolved if the study were performed on condition-specific stratified subnetworks, because based on our results, modules extracted from subnetworks have higher functional enrichment than those from the static network.
Recently, there are debates over whether network hubs fall into two distinct categories of different topological characteristics, specifically named "party" or "date" hubs, as defined by having simultaneous or asynchronous interactions, respectively [32, 59, 61–63]. From the stratification perspective, our results suggest that the hub status of individual proteins may change among networks corresponding to diverse contexts, specific or combinatorial, and therefore different types of hubs may exist, only when taking into account context-dependent subnetworks but not static conglomerate networks alone even with reliable interactions.
Recent studies on revisiting false positive rates in large-scale experimental datasets [54, 56] suggest false positive rates in large-scale data may not be as high as initially assessed . This may be explained from the stratification perspective. A large fraction of interactions in a conglomerate network may no longer exist when they are validated under a specific experimental condition under which such interactions do not take place. These interactions could be taken as false positives; however, when considering other contexts, they may well be true interactions.
A recent study has shown that curated PPI data might not be suitable for topological analysis due to their potential sociological biases . We have included these curated interactions because these interactions are commonly incorporated in conglomerate networks [18, 35]. In this study, although a small portion of curated PPI data are included in the conglomerate network (1,177 vs. 42,131, or 0.028%), we expect their influence to be limited due to the small portion.
Finally, we admit that the most rigorous way to draw conclusions on the differences among networks is to compare a conglomerate network versus PPI networks determined by using the same experimental approach under each specific condition. However, such data are currently unavailable and we do not expect such data to become available soon due to the cost and expense to carry out a large-scale screen of PPIs. Our network stratification approach has provided a close approximation of the subnetworks under specific conditions and thus helps reveal useful insights on the structural and functional differences between a conglomerate network and the context-dependent networks.
Using PPI networks and proteomics data in A. thaliana, we have carried out a proof-of-principle analysis on the structural and functional differences between a conglomerate network and the context-specific networks. We stratified a conglomerate A. thaliana PPI network into context-dependent networks with genome-scale tissue- and condition-specific proteomics data, and systematically analyzed topological, functional and modular differences among the conglomerate and context-dependent networks. We found that the results based on the conglomerate PPI network are often significantly different from the context-dependent subnetworks corresponding to specific tissues or conditions with respect to topological statistics, functional enrichment, and modular components. This conclusion is not particularly dependent on relatively arbitrary cutoffs (such as those defining network hubs or bottlenecks), nor on different module extraction algorithms, nor on the possible high false positive rates in the PPI networks. The consistency in our results suggests the robustness of our analyses.
Due to the limited data availability, we have also performed preliminary analysis on human PPI networks. We found that most of our conclusions are likely true in human as well. We speculate that the observed differences are also likely to exist in other molecular networks, such as gene regulatory networks. In directed networks, we also expect significant changes in network motifs between conglomerate and context-dependent networks. Our results may have implications in various other types of networks such as social networks, collaboration networks, and World Wide Web. Each of these networks could possibly be stratified into subnetworks in the dimension of time, location, or condition while context-dependent information is supplied.
Data Sources for Assembling the Conglomerate A. thaliana PPI Network
To assemble a conglomerate PPI network of A. thaliana, we combined three PPI datasets: predicted PPIs from TAIR containing 3,617 proteins and 19,979 interactions , curated PPIs from TAIR consisting of 770 proteins and 1,177 interactions , and predicted PPIs from AtPID containing 11,706 proteins and 24,418 interactions . The assembled conglomerate network built has a total of 13,136 proteins and 42,131 interactions. We selected the part of the assembled conglomerate network that contains all proteins included in the proteomics data as the total network with 6,606 proteins and 22,165 interactions.
Context-Dependent A. thaliana Proteomics Data
The large-scale A. thaliana proteomic dataset from Baerenfaller et al. includes 18 protein lists with respect to specific or combinatorial contexts, which form a 4-level hierarchy . The top level protein list is the total protein list containing 13,029 proteins identified in any context, which corresponds to the collective of diverse temporal-spatial contexts. The two second level protein lists consist of a list of 10,902 proteins in any specific tissue/organ of A. thaliana, and a list of 8,698 proteins in any of the three cell culture conditions (dark, light, and light small). The third level protein lists include either tissue-specific proteins in the root, leaf, flower, silique and seed of the plant (6,125, 4,853, 9,075, 5,779 and 3,789 proteins, respectively), or condition-specific proteins that express under one of the three cell culture conditions (6,547, 6,474 and 4,472 proteins for dark, light and light small cell culture conditions, respectively). The fourth level protein lists contain sub-tissue-specific proteins, for example, 5,159 proteins that express in the root on the tenth day, or 5,215 proteins identified in the open flowers.
Network Stratification Based on the Conglomerate A. thaliana PPI Network and Context-Dependent Proteomics Data
Generally, network stratification can be performed by first overlaying context-dependent data onto a conglomerate network, followed by selecting and combining interactions that exist in the same contexts to build context-dependent subnetworks. Here we constructed 17 A. thaliana PPI subnetworks corresponding to specific or combinatorial contexts [Figure 2], by overlaying large-scale proteomic lists described in Baerenfaller et al.  onto the assembled conglomerate PPI network, and selecting interacting proteins that both exist under certain conditions. Selected interactions under the same conditions are then combined to form subnetworks corresponding to specific contexts. The coverage of proteins in each of the networks built ranges from 43.5% (1,596 out of 3,665 proteins for the fine-stratified cotyledons subnetwork) to 53.5% (2,390 out of 4,466 proteins for the fine-stratified 23-day roots subnetwork) compared with source protein lists. For the total network, 50.7% (6,606 out of 13,029) of proteins from the source protein list are covered in the network. In this study, network stratification is performed in combined dimensions of time (proteins from different tissues are sampled at different time), location (different tissues/organs of the plant) and condition (different cell culture conditions).
Data Sources for Human PPI Network Stratification
For the stratification of human PPI networks, we used two candidate conglomerate networks: all non-redundant PPIs from the seventh version of the HPRD database , including 9,305 proteins and 35,021 interactions, as well as a pooled dataset of 11,203 proteins and 57,235 interactions assembled in Chuang et al. , which combines PPIs from yeast two-hybrid experiments [16, 65], computational predictions via orthology and co-citation , and curation of the literature [66–68]. Human proteomics data we used are human brain and kidney proteomic datasets from HUPO projects [69, 70]. The brain proteomics data include 1,832 proteins with gene symbols that express in human brain, of which 803 proteins and 1,616 interactions, and 937 proteins and 2,674 interactions can be mapped onto the HPRD PPI network and the pooled PPI network, respectively. The kidney protein list consists of 2,950 non-redundant proteins with gene symbols, of which 1,417 proteins and 3,117 interactions, and 1,593 proteins and 6,498 interactions were present and selected from the two conglomerate networks, respectively.
Identification of Network Hubs and Bottlenecks
In general, network hubs and bottlenecks are nodes with relatively high degree and betweenness values, respectively. Node degree is the number of edges connecting to the node; node betweenness is measured by the number of shortest paths between any pair of nodes in the network that go through the node. In this study, we quantitatively defined nodes with the top 20% of degree values as network hubs, and thus the remaining 80% of nodes are considered non-hubs. Results showed consistency in topology analysis when changing the cutoff percentage to 5 or 10% for network hub definition. For network bottlenecks, we also defined the top 5, 10 or 20% of nodes with highest betweenness values as bottlenecks, and the analysis results showed consistency. We used network tool tYNA to calculate node degree and betweenness values in each network (http://tyna.gersteinlab.org/tyna/; ).
Node Degree/Betweenness Class Assignment
To measure node degree or betweenness value change among networks, we used top 5, 10, 20, 50 percent degree or betweenness values as cutoffs, and put each protein from each of the networks into one of the five classes, representing nodes of similar degree or betweenness values, respectively. To address the issue of rather arbitrariness of such classification, we also used top 10, 20, 50 as percentage cutoffs for four classes, or top 5, 15, 25, 60 as percentage cutoffs for five classes, for node degree and betweenness values, respectively. The analysis results showed consistency, which means node degree/betweenness variations for same proteins during network stratification generally exist.
Compatibility Measurement between Two Subnetworks
where Aν, Bν are node sets, and A e , B e are edge sets of two networks, respectively. This index quantitatively measures the overlap of nodes or edges between two networks.
Module Identification through Network Clustering
where N is the number of modules, L is the total number of edges in the network, l s is the number of edges within module s, and k s is the sum of the degrees of the nodes in modules [49, 72]. Modularity by definition seeks to maximize intra-connections and minimize inter-connections of modules in the network. This algorithm uses simulated annealing strategy to optimize the Modularity score. 3) An edge betweenness centrality-based algorithm that partitions the network iteratively to extract modules . Based on the fact that edges of high betweenness centrality are often connectors between communities, this algorithm iteratively removes an edge with the highest betweenness centrality in the network. When a cutoff is selected, for example, 40 percent of edge removal, the rest of the network, usually partitioned as subnetworks, contains the resulting modules.
Compatibility Measurement between Two Sets of Modules Extracted from Networks
where M and N are numbers of modules extracted from two networks, i and j denote module index, T1 i (or T2 j ) is the number of nodes in each module from one of the two networks, P1 i (or P2 j ) is the number of nodes in the most compatible module of T1 i (or T2 j ) from the other network, and C1 i (or C2 j ) is the number of overlapping nodes between T1 i and P1 i (or between T2 j and P2 j ). Given a module A from a network, we restrict that the most compatible module B for A from another set of modules should have at least 50% of nodes in module A. In this study, a small number of modules are perfectly conserved while stratified, i.e., two modules extracted from two networks (or from one network with different clustering parameters) are identical. In addition, a large number of modules have compatible modules.
Randomization in Network Compatibility Analysis
When measuring network compatibility, we also calculated Jaccard Indices between randomized network pairs for comparison purpose. Randomized networks were built by randomly selecting from the total network the same number of nodes or edges as those of real networks. 500 randomized networks were built for statistical tests and p-value calculations.
Randomization in Network Statistics and Hub/Bottleneck Status Change Analysis
In order to evaluate whether the difference in network statistics, for example average degree or clustering coefficient, between stratified subnetworks and the total network, and whether "hub to non-hub", "non-hub to hub", or "bottleneck to non-bottleneck", "non-bottleneck to bottleneck" changes depend on the change of network sizes, we built randomized networks, assessed the number and percentage of hub/bottleneck status changes from the total network, and calculated p-values by Welch's t-tests and Student's t-tests for statistics and hub status changes, respectively. Twenty randomized networks were constructed by randomly selecting a similar number of proteins as in each of the tissue-specific protein lists, overlaying these random proteins onto the total network, and selecting interacting protein pairs that both exist. (In order to approximate both the node and edge numbers of the original networks, the numbers of nodes and edges in the random networks are often not identical to those in the original networks). The total degree of each of the randomized network must be within 10% of the original network. These randomized networks are scale-free, and comparable in total nodes and interactions to the original tissue-specific subnetworks.
Replacing Interactions in Networks
In order to address the issue about the possible influence of false interactions in the assembled conglomerate A. thaliana PPI network. We randomly removed 10 to 40% of interactions in the unstratified total network, and replaced them with the same number of interactions between random proteins that do not exist in the unstratified network, so that the resulting networks contain 10 to 40% of replaced interactions while maintaining the same size as the total network. This approach is similar to the method used in a former study .
List of abbreviations
We thank Katja Bäerenfaller for sharing A. thaliana tissue- and condition-specific proteomics data with detailed description. We thank Raj Bhatnagar for reading the manuscript critically. We thank Han-Yu Chuang and Trey Ideker for sharing their preprocessed data of an assembled human PPI network. We thank three anonymous referees for their valuable suggestions. MZ was supported in part by ACS grant IRG-92-026-12 and CCHMC Trustee grant awarded to LJL.
- Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature 1999, 402: C47–52. 10.1038/35011540View ArticlePubMedGoogle Scholar
- Eisenberg D, Marcotte EM, Xenarios I, Yeates TO: Protein function in the post-genomic era. Nature 2000, 405: 823–826. 10.1038/35015694View ArticlePubMedGoogle Scholar
- Barabasi AL, Oltvai ZN: Network biology: understanding the cell's functional organization. Nat Rev Genet 2004, 5: 101–113. 10.1038/nrg1272View ArticlePubMedGoogle Scholar
- Ideker T, Sharan R: Protein networks in disease. Genome Res 2008, 18: 644–652. 10.1101/gr.071852.107View ArticlePubMedPubMed CentralGoogle Scholar
- Braun P, Rietman E, Vidal M: Networking metabolites and diseases. Proc Natl Acad Sci USA 2008, 105: 9849–9850. 10.1073/pnas.0805644105View ArticlePubMedPubMed CentralGoogle Scholar
- Xia Y, Yu H, Jansen R, Seringhaus M, Baxter S, Greenbaum D, Zhao H, Gerstein M: Analyzing cellular biochemistry in terms of molecular networks. Annu Rev Biochem 2004, 73: 1051–1087. 10.1146/annurev.biochem.73.011303.073950View ArticlePubMedGoogle Scholar
- Gavin AC, Bosche M, Krause R, Grandi P, Marzioch M, Bauer A, Schultz J, Rick JM, Michon AM, Cruciat CM, et al.: Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature 2002, 415: 141–147. 10.1038/415141aView ArticlePubMedGoogle Scholar
- Li S, Armstrong CM, Bertin N, Ge H, Milstein S, Boxem M, Vidalain PO, Han JD, Chesneau A, Hao T, et al.: A map of the interactome network of the metazoan C. elegans. Science 2004, 303: 540–543. 10.1126/science.1091403View ArticlePubMedPubMed CentralGoogle Scholar
- Lemke N, Heredia F, Barcellos CK, Dos Reis AN, Mombach JC: Essentiality and damage in metabolic networks. Bioinformatics 2004, 20: 115–119. 10.1093/bioinformatics/btg386View ArticlePubMedGoogle Scholar
- Ma HW, Zeng AP: The connectivity structure, giant strong component and centrality of metabolic networks. Bioinformatics 2003, 19: 1423–1430. 10.1093/bioinformatics/btg177View ArticlePubMedGoogle Scholar
- Arita M: The metabolic world of Escherichia coli is not small. Proc Natl Acad Sci USA 2004, 101: 1543–1547. 10.1073/pnas.0306458101View ArticlePubMedPubMed CentralGoogle Scholar
- Balazsi G, Barabasi AL, Oltvai ZN: Topological units of environmental signal processing in the transcriptional regulatory network of Escherichia coli. Proc Natl Acad Sci USA 2005, 102: 7841–7846. 10.1073/pnas.0500365102View ArticlePubMedPubMed CentralGoogle Scholar
- Ma'ayan A, Jenkins SL, Neves S, Hasseldine A, Grace E, Dubin-Thaler B, Eungdamrong NJ, Weng G, Ram PT, Rice JJ, et al.: Formation of regulatory patterns during signal propagation in a Mammalian cellular network. Science 2005, 309: 1078–1083. 10.1126/science.1108876View ArticlePubMedPubMed CentralGoogle Scholar
- Tong AH, Lesage G, Bader GD, Ding H, Xu H, Xin X, Young J, Berriz GF, Brost RL, Chang M, et al.: Global mapping of the yeast genetic interaction network. Science 2004, 303: 808–813. 10.1126/science.1091317View ArticlePubMedGoogle Scholar
- Sharan R, Ulitsky I, Shamir R: Network-based prediction of protein function. Mol Syst Biol 2007, 3: 88. 10.1038/msb4100129View ArticlePubMedPubMed CentralGoogle Scholar
- Rual JF, Venkatesan K, Hao T, Hirozane-Kishikawa T, Dricot A, Li N, Berriz GF, Gibbons FD, Dreze M, Ayivi-Guedehoussou N, et al.: Towards a proteome-scale map of the human protein-protein interaction network. Nature 2005, 437: 1173–1178. 10.1038/nature04209View ArticlePubMedGoogle Scholar
- Qiu P, Wang ZJ, Liu KJ, Hu ZZ, Wu CH: Dependence network modeling for biomarker identification. Bioinformatics 2007, 23: 198–206. 10.1093/bioinformatics/btl553View ArticlePubMedGoogle Scholar
- Chuang HY, Lee E, Liu YT, Lee D, Ideker T: Network-based classification of breast cancer metastasis. Mol Syst Biol 2007, 3: 140. 10.1038/msb4100180View ArticlePubMedPubMed CentralGoogle Scholar
- Yang K, Ma W, Liang H, Ouyang Q, Tang C, Lai L: Dynamic simulations on the arachidonic acid metabolic network. PLoS Comput Biol 2007, 3: e55. 10.1371/journal.pcbi.0030055View ArticlePubMedPubMed CentralGoogle Scholar
- Hopkins AL: Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol 2008, 4: 682–690. 10.1038/nchembio.118View ArticlePubMedGoogle Scholar
- Ma H, Goryanin I: Human metabolic network reconstruction and its impact on drug discovery and development. Drug Discov Today 2008, 13: 402–408. 10.1016/j.drudis.2008.02.002View ArticlePubMedGoogle Scholar
- Kelley R, Ideker T: Systematic interpretation of genetic interactions using protein networks. Nat Biotechnol 2005, 23: 561–566. 10.1038/nbt1096View ArticlePubMedPubMed CentralGoogle Scholar
- Jeong H, Mason SP, Barabasi AL, Oltvai ZN: Lethality and centrality in protein networks. Nature 2001, 411: 41–42. 10.1038/35075138View ArticlePubMedGoogle Scholar
- Guelzim N, Bottani S, Bourgine P, Kepes F: Topological and causal structure of the yeast transcriptional regulatory network. Nat Genet 2002, 31: 60–63. 10.1038/ng873View ArticlePubMedGoogle Scholar
- Herrgard MJ, Swainston N, Dobson P, Dunn WB, Arga KY, Arvas M, Bluthgen N, Borger S, Costenoble R, Heinemann M, et al.: A consensus yeast metabolic network reconstruction obtained from a community approach to systems biology. Nat Biotechnol 2008, 26: 1155–1160. 10.1038/nbt1492View ArticlePubMedPubMed CentralGoogle Scholar
- Luscombe NM, Babu MM, Yu H, Snyder M, Teichmann SA, Gerstein M: Genomic analysis of regulatory network dynamics reveals large topological changes. Nature 2004, 431: 308–312. 10.1038/nature02782View ArticlePubMedGoogle Scholar
- Prill RJ, Iglesias PA, Levchenko A: Dynamic properties of network motifs contribute to biological network organization. PLoS Biol 2005, 3: e343. 10.1371/journal.pbio.0030343View ArticlePubMedPubMed CentralGoogle 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.1089167View ArticlePubMedGoogle Scholar
- Teichmann SA, Babu MM: Gene regulatory network growth by duplication. Nat Genet 2004, 36: 492–496. 10.1038/ng1340View ArticlePubMedGoogle Scholar
- Yu H, Kim PM, Sprecher E, Trifonov V, Gerstein M: The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLoS Comput Biol 2007, 3: e59. 10.1371/journal.pcbi.0030059View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Z, Liu C, Skogerbo G, Zhu X, Lu H, Chen L, Shi B, Zhang Y, Wang J, Wu T, Chen R: Dynamic changes in subgraph preference profiles of crucial transcription factors. PLoS Comput Biol 2006, 2: e47. 10.1371/journal.pcbi.0020047View ArticlePubMedPubMed CentralGoogle Scholar
- Han JD, Bertin N, Hao T, Goldberg DS, Berriz GF, Zhang LV, Dupuy D, Walhout AJ, Cusick ME, Roth FP, Vidal M: Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature 2004, 430: 88–93. 10.1038/nature02555View ArticlePubMedGoogle Scholar
- de Lichtenberg U, Jensen LJ, Brunak S, Bork P: Dynamic complex formation during the yeast cell cycle. Science 2005, 307: 724–727. 10.1126/science.1105103View ArticlePubMedGoogle Scholar
- Lin WH, Liu WC, Hwang MJ: Topological and organizational properties of the products of house-keeping and tissue-specific genes in protein-protein interaction networks. BMC Syst Biol 2009, 3: 32. 10.1186/1752-0509-3-32View ArticlePubMedPubMed CentralGoogle Scholar
- Bossi A, Lehner B: Tissue specificity and the human protein interaction network. Mol Syst Biol 2009, 5: 260. 10.1038/msb.2009.17View ArticlePubMedPubMed CentralGoogle Scholar
- Baerenfaller K, Grossmann J, Grobei MA, Hull R, Hirsch-Hoffmann M, Yalovsky S, Zimmermann P, Grossniklaus U, Gruissem W, Baginsky S: Genome-scale proteomics reveals Arabidopsis thaliana gene models and proteome dynamics. Science 2008, 320: 938–941. 10.1126/science.1157956View ArticlePubMedGoogle Scholar
- Zhang M, Deng J, Fang C, Zhang X, Lu LJ: Biomolecular network analysis and applications. In Knowledge-Based Bioinformatics: From analysis to interpretation. Edited by: Alterovitz G, Ramoni M. John Wiley & Sons; 2010:253–288.Google Scholar
- Greenbaum D, Luscombe NM, Jansen R, Qian J, Gerstein M: Interrelating different types of genomic data, from proteome to secretome: 'oming in on function. Genome Res 2001, 11: 1463–1468. 10.1101/gr.207401View ArticlePubMedGoogle Scholar
- Kislinger T, Cox B, Kannan A, Chung C, Hu P, Ignatchenko A, Scott MS, Gramolini AO, Morris Q, Hallett MT, et al.: Global survey of organ and organelle protein expression in mouse: combined proteomic and transcriptomic profiling. Cell 2006, 125: 173–186. 10.1016/j.cell.2006.01.044View ArticlePubMedGoogle Scholar
- Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, et al.: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucleic Acids Res 2008, 36: D1009–1014. 10.1093/nar/gkm965View ArticlePubMedPubMed CentralGoogle Scholar
- Cui J, Li P, Li G, Xu F, Zhao C, Li Y, Yang Z, Wang G, Yu Q, Shi T: AtPID: Arabidopsis thaliana protein interactome database--an integrative platform for plant systems biology. Nucleic Acids Res 2008, 36: D999–1008. 10.1093/nar/gkm844View ArticlePubMedPubMed CentralGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25: 25–29. 10.1038/75556View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang H, van Nocker S: The vernalization independence 4 gene encodes a novel regulator of flowering locus C. Plant J 2002, 31: 663–673. 10.1046/j.1365-313X.2002.01380.xView ArticlePubMedGoogle Scholar
- Zybailov B, Rutschow H, Friso G, Rudella A, Emanuelsson O, Sun Q, van Wijk KJ: Sorting signals, N-terminal modifications and abundance of the chloroplast proteome. PLoS One 2008, 3: e1994. 10.1371/journal.pone.0001994View ArticlePubMedPubMed CentralGoogle Scholar
- Rutschow H, Ytterberg AJ, Friso G, Nilsson R, van Wijk KJ: Quantitative proteomics of a chloroplast SRP54 sorting mutant and its genetic interactions with CLPC1 in Arabidopsis. Plant Physiol 2008, 148: 156–175. 10.1104/pp.108.124545View ArticlePubMedPubMed CentralGoogle Scholar
- Missiuro PV, Liu K, Zou L, Ross BC, Zhao G, Liu JS, Ge H: Information flow analysis of interactome networks. PLoS Comput Biol 2009, 5: e1000350. 10.1371/journal.pcbi.1000350View ArticlePubMedPubMed CentralGoogle Scholar
- Berardini TZ, Mundodi S, Reiser L, Huala E, Garcia-Hernandez M, Zhang P, Mueller LA, Yoon J, Doyle A, Lander G, et al.: Functional annotation of the Arabidopsis genome using controlled vocabularies. Plant Physiol 2004, 135: 745–755. 10.1104/pp.104.040071View ArticlePubMedPubMed CentralGoogle 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/nature03607View ArticlePubMedGoogle Scholar
- Guimera R, Nunes Amaral LA: Functional cartography of complex metabolic networks. Nature 2005, 433: 895–900. 10.1038/nature03288View ArticlePubMedPubMed CentralGoogle Scholar
- Girvan M, Newman ME: Community structure in social and biological networks. Proc Natl Acad Sci USA 2002, 99: 7821–7826. 10.1073/pnas.122653799View ArticlePubMedPubMed CentralGoogle Scholar
- Samanta MP, Liang S: Predicting protein functions from redundancies in large-scale protein interaction networks. Proc Natl Acad Sci USA 2003, 100: 12579–12583. 10.1073/pnas.2132527100View ArticlePubMedPubMed CentralGoogle Scholar
- Cho YR, Hwang W, Ramanathan M, Zhang A: Semantic integration to identify overlapping functional modules in protein interaction networks. BMC Bioinformatics 2007, 8: 265. 10.1186/1471-2105-8-265View ArticlePubMedPubMed CentralGoogle Scholar
- Yu H, Paccanaro A, Trifonov V, Gerstein M: Predicting interactions in protein networks by completing defective cliques. Bioinformatics 2006, 22: 823–829. 10.1093/bioinformatics/btl014View ArticlePubMedGoogle Scholar
- Yu H, Braun P, Yildirim MA, Lemmens I, Venkatesan K, Sahalie J, Hirozane-Kishikawa T, Gebreab F, Li N, Simonis N, et al.: High-quality binary protein interaction map of the yeast interactome network. Science 2008, 322: 104–110. 10.1126/science.1158684View ArticlePubMedPubMed CentralGoogle Scholar
- von Mering C, Krause R, Snel B, Cornell M, Oliver SG, Fields S, Bork P: Comparative assessment of large-scale data sets of protein-protein interactions. Nature 2002, 417: 399–403. 10.1038/nature750View ArticlePubMedGoogle Scholar
- Huang H, Bader JS: Precision and recall estimates for two-hybrid screens. Bioinformatics 2009, 25: 372–378. 10.1093/bioinformatics/btn640View ArticlePubMedPubMed CentralGoogle Scholar
- Wuchty S: Evolution and topology in the yeast protein interaction network. Genome Res 2004, 14: 1310–1314. 10.1101/gr.2300204View ArticlePubMedPubMed CentralGoogle Scholar
- Calvano SE, Xiao W, Richards DR, Felciano RM, Baker HV, Cho RJ, Chen RO, Brownstein BH, Cobb JP, Tschoeke SK, et al.: A network-based analysis of systemic inflammation in humans. Nature 2005, 437: 1032–1037. 10.1038/nature03985View ArticlePubMedGoogle 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.1522View ArticlePubMedGoogle Scholar
- Wang Z, Zhang J: In search of the biological significance of modular structures in protein networks. PLoS Comput Biol 2007, 3: e107. 10.1371/journal.pcbi.0030107View ArticlePubMedPubMed CentralGoogle Scholar
- Batada NN, Reguly T, Breitkreutz A, Boucher L, Breitkreutz BJ, Hurst LD, Tyers M: Stratus not altocumulus: a new view of the yeast protein interaction network. PLoS Biol 2006, 4: e317. 10.1371/journal.pbio.0040317View ArticlePubMedPubMed CentralGoogle Scholar
- Bertin N, Simonis N, Dupuy D, Cusick ME, Han JD, Fraser HB, Roth FP, Vidal M: Confirmation of organized modularity in the yeast interactome. PLoS Biol 2007, 5: e153. 10.1371/journal.pbio.0050153View ArticlePubMedPubMed CentralGoogle Scholar
- Batada NN, Reguly T, Breitkreutz A, Boucher L, Breitkreutz BJ, Hurst LD, Tyers M: Still stratus not altocumulus: further evidence against the date/party hub distinction. PLoS Biol 2007, 5: e154. 10.1371/journal.pbio.0050154View ArticlePubMedPubMed CentralGoogle Scholar
- Ramani AK, Bunescu RC, Mooney RJ, Marcotte EM: Consolidating the set of known human protein-protein interactions in preparation for large-scale mapping of the human interactome. Genome Biol 2005, 6: R40. 10.1186/gb-2005-6-5-r40View ArticlePubMedPubMed CentralGoogle Scholar
- Stelzl U, Worm U, Lalowski M, Haenig C, Brembeck FH, Goehler H, Stroedicke M, Zenkner M, Schoenherr A, Koeppen S, et al.: A human protein-protein interaction network: a resource for annotating the proteome. Cell 2005, 122: 957–968. 10.1016/j.cell.2005.08.029View ArticlePubMedGoogle Scholar
- Peri S, Navarro JD, Amanchy R, Kristiansen TZ, Jonnalagadda CK, Surendranath V, Niranjan V, Muthusamy B, Gandhi TK, Gronborg M, et al.: Development of human protein reference database as an initial platform for approaching systems biology in humans. Genome Res 2003, 13: 2363–2371. 10.1101/gr.1680803View ArticlePubMedPubMed CentralGoogle Scholar
- Alfarano C, Andrade CE, Anthony K, Bahroos N, Bajec M, Bantoft K, Betel D, Bobechko B, Boutilier K, Burgess E, et al.: The Biomolecular Interaction Network Database and related tools 2005 update. Nucleic Acids Res 2005, 33: D418–424. 10.1093/nar/gki051View ArticlePubMedPubMed CentralGoogle Scholar
- Joshi-Tope G, Gillespie M, Vastrik I, D'Eustachio P, Schmidt E, de Bono B, Jassal B, Gopinath GR, Wu GR, Matthews L, et al.: Reactome: a knowledgebase of biological pathways. Nucleic Acids Res 2005, 33: D428–432. 10.1093/nar/gki072View ArticlePubMedPubMed CentralGoogle Scholar
- Hamacher M, Apweiler R, Arnold G, Becker A, Bluggel M, Carrette O, Colvis C, Dunn MJ, Frohlich T, Fountoulakis M, et al.: HUPO Brain Proteome Project: summary of the pilot phase and introduction of a comprehensive data reprocessing strategy. Proteomics 2006, 6: 4890–4898. 10.1002/pmic.200600295View ArticlePubMedGoogle Scholar
- Miyamoto M, Yoshida Y, Taguchi I, Nagasaka Y, Tasaki M, Zhang Y, Xu B, Nameta M, Sezaki H, Cuellar LM, et al.: In-depth proteomic profiling of the normal human kidney glomerulus using two-dimensional protein prefractionation in combination with liquid chromatography-tandem mass spectrometry. J Proteome Res 2007, 6: 3680–3690. 10.1021/pr070203nView ArticlePubMedGoogle Scholar
- Yip KY, Yu H, Kim PM, Schultz M, Gerstein M: The tYNA platform for comparative interactomics: a web tool for managing, comparing and mining multiple networks. Bioinformatics 2006, 22: 2968–2970. 10.1093/bioinformatics/btl488View ArticlePubMedGoogle Scholar
- Newman ME, Girvan M: Finding and evaluating community structure in networks. Phys Rev E Stat Nonlin Soft Matter Phys 2004, 69: 026113.View ArticlePubMedGoogle 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.