The Symbiosis Interactome: a computational approach reveals novel components, functional interactions and modules in Sinorhizobium meliloti
- Ignacio Rodriguez-Llorente1,
- Miguel A Caviedes1,
- Mohammed Dary1,
- Antonio J Palomares^1,
- Francisco M Cánovas2 and
- José M Peregrín-Alvarez2, 3Email author
© Rodriguez-Llorente et al; licensee BioMed Central Ltd. 2009
Received: 29 January 2009
Accepted: 16 June 2009
Published: 16 June 2009
Rhizobium-Legume symbiosis is an attractive biological process that has been studied for decades because of its importance in agriculture. However, this system has undergone extensive study and although many of the major factors underpinning the process have been discovered using traditional methods, much remains to be discovered.
Here we present an analysis of the 'Symbiosis Interactome' using novel computational methods in order to address the complex dynamic interactions between proteins involved in the symbiosis of the model bacteria Sinorhizobium meliloti with its plant hosts. Our study constitutes the first large-scale analysis attempting to reconstruct this complex biological process, and to identify novel proteins involved in establishing symbiosis. We identified 263 novel proteins potentially associated with the Symbiosis Interactome. The topology of the Symbiosis Interactome was used to guide experimental techniques attempting to validate novel proteins involved in different stages of symbiosis. The contribution of a set of novel proteins was tested analyzing the symbiotic properties of several S. meliloti mutants. We found mutants with altered symbiotic phenotypes suggesting novel proteins that provide key complementary roles for symbiosis.
Our 'systems-based model' represents a novel framework for studying host-microbe interactions, provides a theoretical basis for further experimental validations, and can also be applied to the study of other complex processes such as diseases.
Plant-microbe interactions play an important role in agriculture and a lot of effort has been dedicated to analyse these interactions in detail. One of these interactions is the Rhizobium-Legume symbiosis, a process that allows the growth of the plant in the absence of externally supplied nitrogen. This is a well studied agronomically important process that is also used as a model to study general genetic aspects of plant-microbe interactions [1, 2]. Rhizobial bacteria and legumes have evolved complex signal exchange mechanisms in which a lot of genes are involved . To probe this complexity further we chose to study the model rhizobial symbiont genome Sinorhizobium meliloti . S. meliloti is a model bacterium that can engage in a symbiotic interaction by infecting the roots of members of the genera Medicago and Melilotus, being the S. meliloti-Medicago truncatula interaction the model system for indeterminate type nodules .
The sequencing of hundreds of complete genomes from diverse species is having a tremendous impact on our understanding of biology by enabling the identification of all proteins and the analysis of their function. Despite the vast body of literature about the Rhizobium-legume interaction there have been no systematic large-scale attempts to identify its components and function using a systems biology perspective, and most studies have been restricted to the analysis of individual proteins. However, biological functions results from the interactions of proteins so that understanding the network of biological linkages utilizing functional genomics information is becoming a hot topic in current research projects [6–11]. The main advantage of creating these networks lies in the ability to understand biological processes from a system level perspective. This would ideally require the application of computational and experimental techniques to combine experimental observations of protein-protein interactions (PPIs) and computational predictions derived from different data sources. To date a variety of methods have been developed to derive large scale networks of PPIs for a variety of organisms. These range from experimental methods such as yeast two-hybrid screens, or tandem affinity purification coupled with mass spectrometry [6, 8, 9, 12], to computational methods such as genome context methods [13, 14]. The integration of these types of data helps to provide a complete overview of gene networks of high value for characterizing many biological processes, and ultimately, for understanding the basis of host-microbe interactions including diseases [15–17]. However, experimental information is sometimes missing and deriving gene networks from different computational approaches is not an easy task. Computational predictions such as those obtained by applying genome context methods usually measure functional interactions between proteins. The assumption is that proteins are most likely to interact if: a) their proteins are either present or absent together across multiple genomes (the Phylogenetic Profile method) ; b) a gene fusion event occurred in other species (the Gene Fusion or Rosetta Stone method) [19, 20]; c) the genes are in physical proximity (the Gene Cluster method) ; or d) the genes are conserved in physical proximity and in phylogenetically distant genomes (the Gene Neighbor method) . These methods have the advantage over experimental methods and other computational methods based on protein conservation such as Interologs  or literature mining , that they are not biased towards well studied or conserved proteins or interactions . Therefore, genome context methods are able to highlight organism-specific features since they just rely on genome structure. The outputs derived by these methods can be computationally integrated in order to reconstruct network models of the relations between genes [13, 14]. Data integration for inferring protein associations is advantageous for two main reasons. First, combining data from diverse studies and methods generates data sets of higher quality, and second, integration effectively captures different aspects of organism's biology [25–27]. Further exploiting the topological properties of these networks, clustering algorithms have subsequently allowed proteins to be organized into discrete interconnected units known as functional modules representing either protein complexes or biochemical pathways [28, 29]. In addition, integration of additional functional and comparative genomics data sets are further providing insights into how these modules and their components are co-ordinated and how they may have evolved [9, 30].
Due to the scarcity of large-scale experimental assays aiming to study this important microorganism-host interaction, we chose to apply a systems-based computational approach to evaluate and organize our current knowledge about this complex biological process further. Here, we first reconstruct an extensive and accurate functional network in S. meliloti by integrating the functional associations present in the two well known databases PROLINKS  and STRING  (see methods). These databases host functional linkage predictions obtained mainly by the four different computational genome context approaches described above. Second, we present an analysis of the 'Symbiosis Interactome' (a detailed functional interaction network of the proteins involved in the S. meliloti-Legume symbiosis) by first mapping proteins known to be involved in symbiosis on top of the S. meliloti network, and secondly, by extending this resulting network by means of a novel method, referred here as 'phenotypic profiling', which is further extended by incorporating data from the computational prediction of functional modules. This computational approach potentially revealed the complex interplay of functional interactions between proteins involved in S. meliloti-Medicago symbiosis providing a way to expand the current understanding of symbiosis by enabling hypothesis generation based on our predicted network. Finally, since one of the major advantages of constructing PPI networks is the ability to predict functions for proteins based on their association with well known proteins, we identified and tested the functions of candidate proteins and demonstrate that novel Symbiosis Interactome proteins can still be discovered despite the many decades of effort dedicated to study this important and complex biological process.
The S. meliloti network
The S. meliloti network demonstrated to have properties of scale-free network [see Additional file 1] like other biological networks, the Internet and social networks . Most of the proteins had few interacting partners, where a subset of 'hubs' form a far greater number of connections. Scale-free networks are predicted to be robust against random node removal but vulnerable to hub removal, a property that might be preserved across evolution . Furthermore, the average clustering coefficient (ACC) of the intersection network and its diameter or average shortest path length (L) (see methods) suggests properties of a small-word network (L ~ Lrandom, ACC >> ACCrandom) typical of intracellular network in which the nodes are connected when they are involved in the same biological processes .
Prediction of functional modules
While defining accurate PPI networks is important, the ultimate goal of interactome analyses is to identify the functional modules in these networks, that is, proteins with related functions that tend to be clustered into highly interconnected subnetworks [10, 33, 34], and to validate them. To assess if our network could also be clustered into such subnetworks, we first tested the capacity of the S. meliloti network to form groups of highly interconnected proteins, as indicated by its Average Clustering Coefficient (ACC) (see methods). Indeed, the ACC of the S. meliloti network is much higher (ACC = 0.41) than other large-scale E. coli (ACC = 0.15  and ACC = 0.08 ) and H. pylori  (ACC = 0.02) experimental, and random networks (ACC = 0.0002) suggesting the organization of the S. meliloti network in functional modules.
The Symbiosis Interactome network
We first undertook an exhaustive literature-search analysis to identify and compile a list of bacterial proteins whose role in the symbiosis Rhizobium-Legume has been widely studied (Additional file 2 and methods). These proteins were classified as "classical-known" proteins in different categories according to the stage of symbiosis they are involved in.
Prediction of functional annotation and stage of symbiosis
A major goal for many functional genomics and proteomics projects is the generation of accurate functional information for every gene and its product. Although tremendous progress has been made through the application of such systematic studies, we found that within the S. meliloti proteome 3,376 (54%) proteins were not assigned to a functional category according to COGs, 290 (5%) have been assigned category S (function unknown), and a further 307 (5%) proteins have only been assigned into category 'R' ('general function prediction'). There has been recent progress in the development of novel methods of functional inference based on network connectivity . The availability of our S. meliloti functional network thus provides a valuable resource for future studies aimed at predicting the functions of these high number of functionally 'orphan' proteins. In order to test the ability of our functional network to accurately infer reliable functional annotations and the stage of symbiosis where components of the Symbiosis Interactome may participate, we investigated a basic network-based approach based on functional category membership within predicted functional modules. To provide estimates of the accuracy of functional modules on inferring reliable functional annotations, we applied a cross-validation procedure to predict functional annotations (see methods). We were able to identify correct annotations for 87%–100% of the proteins contained in modules depending on the stringency of COGs category assignments (see methods and Fig. 3c). The accuracy of this type of functional module predictions has been found to be superior to other methods based merely on direct interacting partners [Peregrín-Alvarez JM, Xiong X, Su C, Parkinson J: The modular organization of protein interactions in Escherichia coli, submitted]. These findings highlight both the quality of the network and the predicted functional modules for hypothesis generation and future experimental validation.
Based on these results, module 266, for example, includes three proteins Q92QS6 (Smc01792), Q92QS4 (SMc01794) and Q92VP9 (Smb21071) [see Additional file 2]. The first two proteins are involved in M (cell wall/membrane/envelope biogenesis) while the third one has no COGs category assignment. We therefore predict the latter is potentially involved in this biological process. Furthermore, interestingly, we correctly identify the stage of symbiosis for 92%–100% of the proteins contained in modules depending on stringency (see methods and Fig. 3c). Again based on these promising results, module 208, for example, includes two nodulation proteins: nodP2 (Smb21223) and nodQ2 (Smb21224); and the novel protein Q92VH5 (SMb21225) [see Additional file 2], therefore, being tempting to speculate the participation of the latter in nodulation.
The conservation and evolution of the Symbiosis Interactome network
To investigate the conserved nature and evolution of our predicted Symbiosis Interactome network, the classical-known and novel Symbiosis Interactome components were classified into different node ages according to their phylogenetic distribution (see methods). A total of 313 (~ 68%) proteins were classified as old nodes (with broad phenotypic profiles (i.e with homologs in 7 or 8 phenotypic categories) suggesting an old evolutionary origin for symbiosis [8, 46]. Furthermore, of the 92 classical-known proteins previously identified as components of the Symbiosis Interactome 62 (~ 67%) had homologs with distantly related genomes, indicating that these highly conserved proteins were a valid system from which to derive a model of symbiosis. In addition, highly conserved genes tend to involve essential genes [8, 9]. Since most of the genes known to be involved in symbiosis are highly conserved [see Additional file 2] this suggests that these genes could be essential for organism's survival or at least determinant for symbiosis. Indeed, many of the novel genes predicted by our approach are missing from a S. meliloti mutant collection recently published  (data not shown) suggesting an essential role for many of these novel genes. It has also been shown that nodes with high network connectivity tend to be essential nodes [8, 9, 15]. Since most of the 'classical-known' and other novel Symbiosis Interactome proteins have multiple interacting partners (315 (~ 68%) and 341 (~ 74%) proteins using the Symbiosis Interactome and the complete intersection S. meliloti network, respectively, interact with more than one protein in the network) (see methods), this suggests that these proteins may indeed have a key role in this important biological process. It follows from these findings that the number of interactions of the Symbiosis Interactome proteins are positively correlated with its conservation [see Additional file 1] supporting a model of evolution of the Symbiosis Interactome from core components by adding additional ones over time .
M. sativa plants were inoculated with S. meliloti strains mutated at these genes, using S. meliloti 1021 as control wild-type strain (see methods). We could not observe any difference in nodulation phenotypes between plants inoculated with the strain mutated in Q92TC2 and the 1021 control strain (Fig. 5b). On the other hand, differences in nodulation were observed when plants were inoculated with the other mutants. A 20–30% decrease in nodule number (depending on the experiment these are maximum and minimum values) was observed in plants inoculated with the strain mutated in etfB1, and a 20–25% decrease in nodule number in plants inoculated with the mutant in Q92P53. These differences have been shown as biologically significant in other symbiosis studies [48–50]. In addition, it is important to notice that a high percentage of small nodules (white and probably non-fixing nodules) was also observed in plants inoculated with etfB1 mutant. Surprisingly, plants inoculated with the strain mutated in msbA1 showed a 20–25% increase in nodule number when compared with control strain (Fig. 5b). In summary, these results clearly suggest that still there could be a number of non-described proteins involved in the Rhizobium-Legume interaction.
Further functional predictions
Based on our experimental results and the interactions of the novel targeted proteins, etfB1 acts in a module involved in energy production and conversion, and we predict it to be potentially involved in nitrogen fixation [see Additional file 2]; in fact, the high percentage of small non-fixing nodules induced by the strain mutated in this gene is consistent with this role. msbA1 functions in a module together with ndvA and exsA genes and is potentially involved in glucan synthesis; and Q92P53 is functioning within a module involved in lipid transport and metabolism in coordination with nod genes, and may be potentially involved in the regulation of the first stages of nodule formation. These novel findings only represents hypothesis and still have to be analysed in more detail to shed more light on their precise biological role and mechanistic details but, nonetheless, the predictions highlighted here represent a tempting guide for further experimental validation.
The building of our final 'Symbiosis Interactome network' complemented our initial classical-known list in many different ways. First, we extended the initial set from 92 to 163 known components (92 from the intersection and 71 from the union network). Second, we identified 263 potential novel Symbiosis Interactome components, representing ideal targets for further experimental validation [see Additional file 2]. Third, the incorporation of functional modules in the network provides additional information concerning the structure and functional organization of the Symbiosis Interactome. Interestingly, functional modules tend to be formed by proteins involved in the same stage of symbiosis [see Additional file 2] suggesting that distinct symbiosis-stages are organized and coordinated as distinct functional modules. Therefore, the incorporation of modules apart from providing another structural dimension to the Symbiosis Interactome also allows the prediction of both protein function and the symbiosis-stage a novel component may participate (Fig. 3c). This highlight both the quality of the network and the functional modules we predicted as guide for direct experimental validation. The final 'Symbiosis Interactome network', therefore, hosts the organization of the Symbiosis Interactome into functional interactions and modules, and constitutes the first attempt toward the representation of this complex biological process (Fig. 4).
Novel predicted components include many conserved proteins of unknown functions and others participating in a variety of cellular processes (Fig. 4). Novel proteins may represent false negatives components not identified by current experimental techniques perhaps because they are highly specialized components or maybe recruited to the Symbiosis Interactome under specific conditions that have escaped from detection and are therefore absent from our 'classical-known' preliminary data. Our experimental results yielded a preliminary notable success (3 positive cases out of 4 proteins tested experimentally) for predicting novel S. meliloti-M. sativa symbiotic components by using our computational approach. The results also provide tempting clues in regard to the predictive potential of our approach for hypothesis generation and guiding future experimental validation. For example, the two module-network scenarios presented here suggest high accuracy at predicting novel components and functional modules. Furthermore, high scored interactions based on our probability scores are experimentally validated as opposed to low quality interactions for which we could not find any direct experimental evidence, at least not for the gene Q92TC2 tested here. For this particular protein and the remaining 259 non-tested novel proteins, it is difficult to determine how many of them could be really involved in this complex biological process. It has been described that mutations in some bacterial nodulation genes do not have any influence in the symbiotic properties of the bacteria. For example, S. meliloti cells mutated in fixT gene are not affected in nodulation with M. sativa host plants . The expression of this fixation gene is regulated by FixH protein, which is essential for nodulation (mutations in fixH gives a Fix- phenotype, that is, non-fixing nodules). It has been suggested that some nodulation proteins could have a role in symbiosis when the expression of essential proteins is blocked. In the same manner, there are proteins that could be essential for nodulation in special situations, such as biotic and abiotic stress. In addition, there are proteins that could be involved in the symbiotic competitiveness of the rhizobial strain. Finally, another alternative explanation is that the potential involvement of the gene Q92TC2 in symbiosis might be compensated by other genes performing similar functions. Indeed, a gene family analysis by using sequence similarity clustering through the MCL algorithm  (see methods) revealed an intriguing gene family expansion in this particular case (31 genes in this family), whereas in the other 3 mutated genes we do not observe such drastic family expansions (with 1 (singleton family), 3, and 14 gene family members, for the genes Q92P53, etfB1, and msbA1, respectively). This interesting result suggests that other members of this large gene family might rescue its potential role in symbiosis through the establishment of backup circuits, such as occurs in other well studied model organisms . There is evidence of direct backup compensation between gene duplicates with overlapping functions where one gene can cover for the loss of its paralogue, and sometimes these compensations occur only for certain functions under given conditions . In all these situations, the single mutation of these genes in conventional laboratory conditions would not be the best experiment to assess their role in symbiosis. We believe this novel finding supports the model of network robustness through gene duplication , and it also has very interesting implications regarding the selection of the right candidate genes and experimental method in future validation studies.
While the functional network presented here provides valuable clues about the components of the bacterial Symbiosis Interactome, the main limitation of our study is the lack of experimental information on PPIs which made us to consider as input only computationally derived functional genomics data. Integration of computational approaches with recently published  and future experimental interaction data would likely improve the quality of our network and the prediction of novel components. This can be done by using Bayesian or probabilistic models shown to result in accurate confidence scoring systems [[26, 27, 55], Peregrín-Alvarez JM, Xiong X, Su C, Parkinson J: The modular organization of protein interactions in Escherichia coli, submitted]. Furthermore, although we believe we have been very flexible by allowing interactions between proteins with potential phenotypic profiles and not directly interacting with the giant-central network component, our Symbiosis Interactome network can still serve as a platform to add other interactions and components potentially involved in symbiosis. For example, we can choose other proteins with other interesting phenotypic profiles to extend our network such as those profiles showing homologs in other symbionts and/or pathogenic species since these bacteria often use the same core molecular mechanism to maintain their associations with hosts . Future analyses will also include further network extensions based on recently characterized symbiosis components [57–59], inclusion of other interesting phenotypic profiles (see above), a larger-scale experimental validation of the novel components predicted to be involved in symbiosis, and further analyses of the components and pathways involved in host-microbe, and host (i.e. plant) interactions. Finally, through an iterative process, novel Symbiosis Interactome components once experimentally confirmed, can be then added to the known set, potentially increasing the list of novel components and finally revealing the complete picture of the Symbiosis Interactome network.
The essential contribution of symbiosis to understand host-microbe interactions underscores the importance of further studying the structure and organization of the Symbiosis Interactome. Here we presented a novel 'systems-based model' that provided for the very first time new insights into the functional organization of the S. meliloti Symbiosis Interactome and the necessary framework on which to build, in an iterative manner, to further our understanding of symbiosis. We have identified 263 potential novel symbiosis components, and have demonstrated experimentally the participation of novel proteins involved in this important process. These novel proteins might not be essential for symbiosis but still determinant for the microbe-plant interaction since most of the essential components for this process have been described through decades of effort. Understanding the biology of this important model organism is essential not only for having a network view of how this biological process functions at a molecular level but also for the development of anti-microbial drugs since many of the proteins and modules involved in bacterial-symbiosis may be conserved, and thus, performing similar functions, in other microbial pathogens . Furthermore, we can use our network as a template to derive other Symbiosis Interactome networks for other bacteria-related species which is particularly important given the difficulty and cost of obtaining high throughput screens. Those maps should provide an useful starting point for predicting functional interactions and modules, and the function of unknown proteins. It remain to be seen which of these interactions and components do indeed occur and what is the specific role they play in each of these organisms. We believe that this model adds a new view and dimension to our understanding of host-microbe interactions, and can be extended to study other complex biological processes such as those involved in diseases.
An initial list of proteins known to be involved in the Rhizobium-Legume symbiosis was obtained and manually curated using PubMed, Google, journal-specific searches, and literature reviews and citations. We have called this list the 'classical-known' set.
We used S. meliloti genome context data from the PROLINKS  and STRING  databases. While both databases use the same genome context methods to derive functional linkages they both differ in the statistical procedures and scoring systems they use to provide high quality interactions. We reasoned that the overlap between both databases (intersection) represents interactions more likely to be true positives, and that the union of both databases represents a dataset with higher coverage (see below). We used all medium-to-high confidence functional linkages provided by the STRING database. From PROLINKS database we used those functional linkages in S. meliloti over 0.6 confidence. This cut-off provided a true positive rate similar to the one obtained by using the medium-to-high confidence data from the STRING database. The genome context data obtained from these two databases were combined into two single non-redundant datasets: one based on the overlapping between these two databases (the intersection dataset), and another one based on the union of the databases (the union dataset).
The confidence scores associated to each functional linkage provided by the original STRING and PROLINKS databases were re-scored according to the following criteria: STRING provides unified scores representing the confidence of a given functional linkage. The bigger the score, the more reliable the interaction. We reasoned that those interactions present in both databases are the most reliable ones, and we tested it by calculating ROC curves (see below). STRING scores were transformed into a scoring scale 0 – 0.5, the closer to 0.5, the bigger the confidence of the interaction. PROLINKS provides independent confidence scores for each applied independent genome context method. The scores were combined into an unified score by summing all confidence scores for a particular functional linkage and transforming the resulting number to a 0 – 0.5 scale. This procedure resulted in a 0 – 1 confidence score for those functional linkages present in both databases (the intersection data set) and a 0 – 0.5 confidence score for those interactions present in only one of the databases.
The validity of our re-scoring approach and the integrated networks was tested by calculating Receiving Operating Curves (ROC) and the Area Under the Curve (AUC) of the intersection, union, PROLINKS and STRING data sets as a measure of accuracy.
To be able to calculate accurate ROC curves and AUCs it is crucial to complement a positive gold standard set with a negative one. Because a reference set of known interactions is not available for S. meliloti, here we consider as positive set those functional linkages belonging to the same COGs functional category [37, 60]. The construction of a negative set is rather problematic because it is impossible to be sure that two proteins do not interact. However, by using those pairs of proteins that are present in different COGs functional categories and do not colocalize in the same cell compartment it is possible to make a list of protein pairs that are unlikely to interact, thus representing a good approximation to a negative set. COGs annotations were mapped to functional linkages and the periplasmic location of all the proteins was predicted (see below).
The periplasmic location of all S. meliloti proteins was predicted using the SIGCLEAVAGE software . The proteins were considered as periplasmic if they contained at least one predicted signal sequence within 50 residues from the N terminus. The proteins that did not contain any signal sequence throughout the entire sequence were considered cytoplasmic, and the remaining proteins were not classified.
The protein IDs of the functional linkage data from the intersection and union networks were converted to gene names using UNIPROT database  and these were used to map the list of classical-known proteins in S. meliloti.
For each S. meliloti sequence, a BLASTP  search was performed against 200 complete genome datasets. Both S. meliloti and other complete genomes were downloaded from the COGENT database  [see Additional file 2]. Homologs for each S. meliloti protein were determined based on a raw bit score threshold of 50, and were used to generate phenotypic profiles as follow: the complete genomes were manually curated and assigned to the following 8 phenotype categories [see Additional file 2] using PubMed, Google, and other web-specific searches: C, root colonizing bacteria; Fn, nitrogen-fixing bacteria in symbiosis with plants; Fl, free living nitrogen-fixing bacteria; P, pathogen; Pp, plant pathogen; S, soil cohabitant; Sy, symbiont/commensal; and O, other organisms. The only restriction for categorizing was that the genome in question have to be classified into one category only, the one with the most relevant phenotype for the study of symbiosis [see Additional file 2]. For example, if a genome could be classified as C and S, we considered only the category C because all C are also category S; or if a genome could be classified as P and S, P was considered more important for our analysis and thus classified as P only.
We then built phylogenetic profiles  for each S. meliloti protein and mapped the phenotypic data on top of the phylogenetic profiles yielding what we term 'phenotypic profiles'. For example, a protein with a phenotypic profile "FnFl" stands for a protein with homologs in plant nitrogen fixing bacteria and free-living nitrogen fixing bacteria only, thus representing a protein that may be potentially involved in symbiosis.
Unless otherwise noted network analyses were performed using Perl scripts developed in house. The degree (k) of a node (protein) in an interaction network is defined by the number of interactions of the node with other nodes in the network. For a node of degree k, its clustering coefficient (CC) is defined as 2N/k(k-1), where N is the number of interactions between the node's k neighbors and k(k-1)/2 is the number of possible interactions between its neighbors. A CC of 1 means that all the neighbors of a node are fully interconnected. The shortest path length between two nodes in the network is the number of edges in a shortest path connecting them. The shortest path length is infinity if there are no paths between two nodes. Network diameters were obtained using Pajek , and cluster coefficients and shortest path lengths were obtained using tYNA .
To act as controls, random networks were created by randomly selecting equal numbers of proteins (compared with the comparator network) from the S. meliloti network and randomly connecting them with equal numbers of interactions.
where A and B represents different networks, SAB the similarity (i.e. the frequency of common interactions) of A versus B, and SBA the similarity of B versus A.
Detection of functional modules
We identified highly connected functional modules operating within the intersection S. meliloti network by using the Markov Cluster (MCL) algorithm . MCL was applied to our S. meliloti network by testing several inflation operators, and settling on values that provided the highest clusters size, and the best overlap (semantic similarity)  of the computed clusters with the functional categories of the highly curated database COGs .
To compute the significant of finding specific COGs modules, a p-value for each module was calculated based on the distribution of 10,000 random module sets of the same size (assuming a normal distribution) and our module predictions, therefore, representing the probability of seeing such modules at chance. COGs categories with general function prediction, unknown or unassigned were not considered in this analysis. Only modules with at least three components with COGs assignments were statistically computed.
Prediction of functional annotation and stage of symbiosis
Predictions of functional annotation and the stage of symbiosis were performed using enrichment of COGs terms in functional modules (see above). Module prediction for a protein employed the predicted functional modules and derived COGs/symbiosis-stage annotations for the target proteins based on the highest percentage of common COGs/symbiosis-stage terms among the different components of the functional module. Correct COGs/symbiosis-stage assignments additionally required at least 20% of the interaction module components to have the same COGs/symbiosis-stage category. Two measures of stringency were employed: high stringency predictions required the majority of interaction module components to be assigned to the same COGs/symbiosis-stage category; low stringency predictions only required any of the interaction module components to possess the same COGs/symbiosis-stage category (albeit with the additional proviso that at least 20% of the module partners were so annotated). To measure the accuracy of module predictions we used a leave-one-out (LOO) cross-validation procedure, i.e. only proteins which itself and one of its module components possessed an annotation were used in cross-validation. The LOO method randomly selects a protein and compares its known annotation with that predicted by the functional module method.
Gene family analyses
S. meliloti mutants
S. meliloti mutants were obtained from a Mini-Tn5 transposon library constructed in the Lehrstuhl für Genetik (Bielefeld University, Germany) . Based on four network scenarios (see results) we selected the following S. meliloti mutants for experimental validation of our approach: 2011mTn5 STM.3.02.D12_transposon(etfB1), 2011mTn5 STM.4.10.F09_transposon(Q92P53), 2011mTn5 STM.3.08.C10_transposon(Q92TC2), and 2011mTn5 STM.1.06.E11_transposon (msbA1).
Seeds of alfalfa (Medicago sativa L. ecotype. Aragon) were surface sterilised on 70% ethanol for 10 minutes, exhaustively washed on distilled water and placed in water-agar plates for 36 hours at 22°C in the dark. 0.5–1 cm root pre-germinated seedlings were carefully transferred to squared plates containing a slope of BNM-agar medium . Seedlings were inoculated with 100 μl of an overnight culture of S. meliloti mutants or the strain 1021 as control. The lower part of the plate was covered with black paper in order to avoid the roots getting exposed to light. Plates were placed on an Ibercex G-28 plant growth cabinet at 22°C with 16 hours photoperiod. Plants were taken out of the plates at 28 days post-inoculation (dpi) for nodule analyses (counting, size, color, etc). Three independent experiments with 50 plants per experiment were done (150 plants in total). General aspects of plants were also analysed.
This work was supported by operating funds from the Canadian Institutes of Health Research (CIHR). P-A.JM also acknowledges the Ramon y Cajal Program from the Ministry of Science and Technology, Spain. R-L.I, C.MA and D.M also acknowledge Dr. Eloisa Pajuelo and the Spanish Ministry of Education (Plant Biotechnology Program) for financial support. We would like to thank Dr. Anke Becker (Lehrstuhl für Genetik, Bielefeld University, Germany) for providing S. meliloti mutants.
- Graham PH, Vance CP: Legumes: importance and constraints to greater use. Plant Physiol. 2003, 131: 872-877. 10.1104/pp.017004PubMed CentralView ArticlePubMedGoogle Scholar
- Stacey G, Libault M, Brechenmacher L, Wan J, May GD: Genetics and functional genomics of legume nodulation. Curr Opin Plant Biol. 2006, 9: 110-121. 10.1016/j.pbi.2006.01.005View ArticlePubMedGoogle Scholar
- Jones KM, Kobayashi H, Davies BW, Taga ME, Walker GC: How rhizobial symbionts invade plants: the Sinorhizobium-Medicago model. Nat Rev Microbiol. 2007, 5: 619-33. 10.1038/nrmicro1705PubMed CentralView ArticlePubMedGoogle Scholar
- Galibert F, Finan TM, Long SR, Puhler A, Abola P, Ampe F, Barloy-Hubler F, Barnett MJ, Becker A, Boistard P, Bothe G, Boutry M, Bowser L, Buhrmester J, Cadieu E, Capela D, Chain P, Cowie A, Davis RW, Dreano S, Federspiel NA, Fisher RF, Gloux S, Godrie T, Goffeau A, Golding B, Gouzy J, Gurjal M, Hernandez-Lucas I, Hong A, et al.: The composite genome of the legume symbiont Sinorhizobium meliloti. Science. 2001, 293: 668-72. 10.1126/science.1060966View ArticlePubMedGoogle Scholar
- Cook DR: Medicago truncatula – a model in the making!. Curr Opin Plant Biol. 1999, 2: 301-304. 10.1016/S1369-5266(99)80053-3View ArticlePubMedGoogle Scholar
- Rain JC, Selig L, De Reuse H, Battaglia V, Reverdy C, Simon S, Lenzen G, Petel F, Wojcik J, Schächter V, Chemama Y, Labigne A, Legrain P: The protein-protein interaction map of Helicobacter pylori. Nature. 2001, 409: 211-215. 10.1038/35051615View ArticlePubMedGoogle Scholar
- Arifuzzaman M, Maeda M, Itoh A, Nishikata K, Takita C, Saito R, Ara T, Nakahigashi K, Huang HC, Hirai A, Tsuzuki K, Nakamura S, Altaf-Ul-Amin M, Oshima T, Baba T, Yamamoto N, Kawamura T, Ioka-Nakamichi T, Kitagawa M, Tomita M, Kanaya S, Wada C, Mori H: Large-scale identification of protein-protein interaction of Escherichia coli K-12. Genome Res. 2006, 16: 686-691. 10.1101/gr.4527806PubMed CentralView ArticlePubMedGoogle Scholar
- Butland G, Peregrín-Alvarez JM, Li J, Yang W, Yang X, Canadien V, Starostine A, Richards D, Beattie B, Krogan N, Davey M, Parkinson J, Greenblatt J, Emili A: Interaction network containing conserved and essential protein complexes in Escherichia coli. Nature. 2005, 433: 531-537. 10.1038/nature03239View ArticlePubMedGoogle Scholar
- Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, Li J, Pu S, Datta N, Tikuisis AP, Punna T, Peregrín-Alvarez JM, Shales M, Zhang X, Davey M, Robinson MD, Paccanaro A, Bray JE, Sheung A, Beattie B, Richards DP, Canadien V, Lalev A, Mena F, Wong P, Starostine A, Canete MM, Vlasblom J, Wu S, Orsi C, et al.: Global landscape of protein complexes in the yeast Saccharomyces cerevisiae. Nature. 2006, 440: 637-643. 10.1038/nature04670View ArticlePubMedGoogle Scholar
- Li S, Armstrong CM, Bertin N, Ge H, Milstein S, Boxem M, Vidalain PO, Han JD, Chesneau A, Hao T, Goldberg DS, Li N, Martinez M, Rual JF, Lamesch P, Xu L, Tewari M, Wong SL, Zhang LV, Berriz GF, Jacotot L, Vaglio P, Reboul J, Hirozane-Kishikawa T, Li Q, Gabel HW, Elewa A, Baumgartner B, Rose DJ, Yu H, et al.: A map of the interactome network of the metazoan C. elegans. Science. 2004, 303: 540-543. 10.1126/science.1091403PubMed CentralView ArticlePubMedGoogle Scholar
- Yellaboina S, Goyal K, Mande SC: Inferring genome-wide functional linkages in E. coli by combining improved genome context methods: comparison with high-throughput experimental data. Genome Res. 2007, 17: 527-535. 10.1101/gr.5900607PubMed CentralView ArticlePubMedGoogle Scholar
- Ho Y, Gruhler A, Heilbut A, Bader GD, Moore L, Adams SL, Millar A, Taylor P, Bennett K, Boutilier K, Yang L, Wolting C, Donaldson I, Schandorff S, Shewnarane J, Vo M, Taggart J, Goudreault M, Muskat B, Alfarano C, Dewar D, Lin Z, Michalickova K, Willems AR, Sassi H, Nielsen PA, Rasmussen KJ, Andersen JR, Johansen LE, Hansen LH, et al.: Systematic identification of protein complexes in Saccharomyces cerevisiae by mas spectrometry. Nature. 2002, 415: 180-183. 10.1038/415180aView ArticlePubMedGoogle Scholar
- Bowers PM, Pellegrini M, Thompson MJ, Fierro J, Yeates TO, Eisenberg D: Prolinks: a database of protein functional linkages derived from coevolution. Genome Biol. 2004, 5: R35- 10.1186/gb-2004-5-5-r35PubMed CentralView ArticlePubMedGoogle Scholar
- von Mering C, Jensen LJ, Kuhn M, Chaffron S, Doerks T, Krüger B, Snel B, Bork P: STRING 7 – recent developments in the integration and prediction of protein interactions. Nucleic Acids Res. 2007, 35: D358-362. 10.1093/nar/gkl825PubMed CentralView ArticlePubMedGoogle 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
- Ideker T, Thorsson V, Ranish JA, Christmas R, Buhler J, Eng JK, Bumgarner R, Goodlett DR, Aebersold R, Hood L: Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. Science. 2001, 292: 929-934. 10.1126/science.292.5518.929View ArticlePubMedGoogle Scholar
- Strong M, Mallick P, Pellegrini M, Thompson MJ, Eisenberg D: Inference of protein function and protein linkages in Mycobacterium tuberculosis based on prokaryotic genome organization: a combined computational approach. Genome Biol. 2003, 4: R59- 10.1186/gb-2003-4-9-r59PubMed CentralView ArticlePubMedGoogle Scholar
- Pellegrini M, Marcotte EM, Thompson MJ, Eisenberg D, Yeates TO: Assigning protein functions by comparative genome analysis: protein phylogenetic profiles. Proc Natl Acad Sci USA. 1999, 96: 4285-4288. 10.1073/pnas.96.8.4285PubMed CentralView ArticlePubMedGoogle Scholar
- Marcotte EM, Pellegrini M, Ng HL, Rice DW, Yeates TO, Eisenberg D: Detecting protein function and protein-protein interactions from genome sequences. Science. 1999, 285: 751-753. 10.1126/science.285.5428.751View ArticlePubMedGoogle Scholar
- Enright AJ, Iliopoulos I, Kyrpides NC, Ouzounis CA: Protein interaction maps for complete genomes based on gene fusion events. Nature. 1999, 402: 86-90. 10.1038/47056View ArticlePubMedGoogle Scholar
- Overbeek R, Fonstein M, D'Souza M, Pusch GD, Maltsev N: The use of gene clusters to infer functional coupling. Proc Natl Acad Sci USA. 1999, 96: 2896-2901. 10.1073/pnas.96.6.2896PubMed CentralView ArticlePubMedGoogle Scholar
- Yu H, Luscombe NM, Lu HX, Zhu X, Xia Y, Han JD, Bertin N, Chung S, Vidal M, Gerstein M: Annotation transfer between genomes: protein-protein interologs and protein-DNA regulogs. Genome Res. 2004, 14: 1107-1118. 10.1101/gr.1774904PubMed CentralView ArticlePubMedGoogle Scholar
- Hoffmann R, Valencia A: Implementing the iHOP concept for navigation of biomedical literature. Bioinformatics. 2005, 21 (Suppl 2): ii252-258. 10.1093/bioinformatics/bti1142View ArticlePubMedGoogle Scholar
- Hakes L, Pinney JW, Robertson DL, Lovell SC: Protein-protein interaction networks and biology-what's the connection?. Nat Biotech. 2008, 26: 69-72. 10.1038/nbt0108-69.View ArticleGoogle 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
- Jansen R, Yu H, Greenbaum D, Kluger Y, Krogan NJ, Chung S, Emili A, Snyder M, Greenblatt JF, Gerstein M: A Bayesian networks approach for predicting protein-protein interactions from genomic data. Science. 2003, 302: 449-453. 10.1126/science.1087361View ArticlePubMedGoogle Scholar
- Lee I, Date SV, Adai AT, Marcotte EM: A Probabilistic Functional Network of Yeast Genes. Science. 2004, 306: 1555-1558. 10.1126/science.1099511View ArticlePubMedGoogle Scholar
- Gavin AC, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, Rau C, Jensen LJ, Bastuck S, Dümpelfeld B, Edelmann A, Heurtier MA, Hoffman V, Hoefert C, Klein K, Hudak M, Michon AM, Schelder M, Schirle M, Remor M, Rudi T, Hooper S, Bauer A, Bouwmeester T, Casari G, Drewes G, Neubauer G, Rick JM, Kuster B, Bork P, et al.: Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006, 440: 631-636. 10.1038/nature04532View ArticlePubMedGoogle Scholar
- Pu S, Vlasblom J, Emili A, Greenblatt J, Wodak SJ: Identifying functional modules in the physical interactome of Saccharomyces cerevisiae. Proteomics. 2007, 7: 944-960. 10.1002/pmic.200600636View 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
- Salwinski L, Miller CS, Smith AJ, Pettit FK, Bowie JU, Eisenberg D: The Database of Interacting Proteins: 2004 update. Nucleic Acids Res. 2004, 32: D449-451. 10.1093/nar/gkh086PubMed CentralView ArticlePubMedGoogle Scholar
- Wuchty S, Oltvai ZN, Barabási AL: Evolutionary conservation of motif constituents in the yeast protein interaction network. Nat Genet. 2003, 35: 176-179. 10.1038/ng1242View ArticlePubMedGoogle Scholar
- Giot L, Bader JS, Brouwer C, Chaudhuri A, Kuang B, Li Y, Hao YL, Ooi CE, Godwin B, Vitols E, Vijayadamodar G, Pochart P, Machineni H, Welsh M, Kong Y, Zerhusen B, Malcolm R, Varrone Z, Collis A, Minto M, Burgess S, McDaniel L, Stimpson E, Spriggs F, Williams J, Neurath K, Ioime N, Agee M, Voss E, Furtak K, et al.: A protein interaction map of Drosophila melanogaster. Science. 2003, 302: 1727-1736. 10.1126/science.1090289View ArticlePubMedGoogle Scholar
- Schwikowski B, Uetz P, Fields S: A network of protein-protein interactions in yeast. Nature Biotech. 2000, 18: 1257-1261. 10.1038/82360.View ArticleGoogle Scholar
- van Dongen S: Graph clustering by flow simulation. PhD thesis. 2000, University of Utrecht, The NetherlandsGoogle Scholar
- Kelley BP, Sharan R, Karp RM, Sittler T, Root DE, Stockwell BR, Ideker T: Conserved pathways within bacteria and yeast as revealed by global protein network alignment. Proc Natl Acad Sci USA. 2003, 100: 11394-11399. 10.1073/pnas.1534710100PubMed CentralView ArticlePubMedGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4: 41- 10.1186/1471-2105-4-41PubMed CentralView ArticlePubMedGoogle Scholar
- Sarachu M, Colet M: wEMBOSS: a web interface for EMBOSS. Bioinformatics. 2005, 21: 540-541. 10.1093/bioinformatics/bti031View ArticlePubMedGoogle Scholar
- Djordjevic MA, Chen HC, Natera S, Van Noorden G, Menzel C, Taylor S, Renard C, Geiger O, Weiller GF, : A global analysis of protein expression profiles in Sinorhizobium meliloti : discovery of new genes for nodule occupancy and stress adaptation. Mol Plant Microbe Interact. 2003, 16: 508-524. 10.1094/MPMI.2003.16.6.508View ArticlePubMedGoogle Scholar
- Djordjevic MA: Sinorhizobium meliloti metabolism in the root nodule: a proteomic perspective. Proteomics. 2004, 4: 1859-1872. 10.1002/pmic.200300802View ArticlePubMedGoogle Scholar
- Mauchline TH, Fowler JE, East AK, Sartor AL, Zaheer R, Hosie AH, Poole PS, Finan TM: Mapping the Sinorhizobium meliloti 1021 solute-binding protein-dependent transportome. Proc Natl Acad Sci USA. 2006, 103: 17933-17938. 10.1073/pnas.0606673103PubMed CentralView ArticlePubMedGoogle Scholar
- Capela D, Filipe C, Bobik C, Batut J, Bruand C: Sinorhizobium meliloti differentiation during symbiosis with alfalfa: a transcriptomic dissection. Mol Plant Microbe Interact. 2006, 19: 363-372. 10.1094/MPMI-19-0363View ArticlePubMedGoogle Scholar
- Krol E, Becker A: Global transcriptional analysis of the phosphate starvation response in Sinorhizobium meliloti strains 1021 and 2011. Mol Genet Genomics. 2004, 272: 1-17. 10.1007/s00438-004-1030-8View ArticlePubMedGoogle Scholar
- Kuehn MJ, Kesty NC: Bacterial outer membrane vesicles and the host-pathogen interaction. Genes Dev. 2005, 19: 2645-55. 10.1101/gad.1299905View ArticlePubMedGoogle Scholar
- Sharan R, Ulitsky I, Shamir R: Network-based prediction of protein function. Mol Syst Biol. 2007, 3: 88- 10.1038/msb4100129PubMed CentralView ArticlePubMedGoogle Scholar
- Pereira-Leal JB, Audit B, Peregrín-Alvarez JM, Ouzounis CA: An exponential core in the heart of the yeast protein interaction network. Mol Biol Evol. 2005, 22: 421-425. 10.1093/molbev/msi024View ArticlePubMedGoogle Scholar
- Pobigaylo N, Wetter D, Szymczak S, Schiller U, Kurtz S, Meyer F, Nattkemper TW, Becker A: Construction of a large signature-tagged mini-Tn5 transposon library and its application to mutagenesis of Sinorhizobium meliloti. Appl Environ Microbiol. 2006, 72: 4329-4337. 10.1128/AEM.03072-05PubMed CentralView ArticlePubMedGoogle Scholar
- Lorio JC, Kim WS, Krishnan HB: NopB, a soybean cultivar-specificity protein from Sinorhizobium fredii USDA 257, is a type III secreted protein. Mol Plant Microbe Interact. 2004, 17: 1259-1268. 10.1094/MPMI.2004.17.11.1259View ArticlePubMedGoogle Scholar
- Townsend GE, Forsberg LS, Keating DH: Mesorhizobium loti produces nodPQ-dependent sulfated cell surface polysaccharides. J Bacteriol. 2006, 188: 8560-8572. 10.1128/JB.01035-06PubMed CentralView ArticlePubMedGoogle Scholar
- Piñero S, Rivera J, Romero D, Cevallos MA, Martínez A, Bolívar F, Gosset G: Tyrosinase from Rhizobium etli is involved in nodulation efficiency and symbiosis-associated stress resistance. J Mol Microbiol Biotechnol. 2007, 13: 35-44. 10.1159/000103595View ArticlePubMedGoogle Scholar
- Foussard M, Garnerone AM, Ni F, Soupène E, Boistard P, Batut J: Negative autoregulation of the Rhizobium meliloti fixK gene is indirect and requires a newly identified regulator, FixT. Mol Microbiol. 1997, 25: 27-37. 10.1046/j.1365-2958.1997.4501814.xView ArticlePubMedGoogle Scholar
- Stein A, Aloy P: A molecular interpretation of genetic interactions in yeast. FEBS Lett. 2008, 582: 1245-1250. 10.1016/j.febslet.2008.02.020View ArticlePubMedGoogle Scholar
- Gu X: Evolution of duplicate genes versus genetic robustness against null mutations. Trends Genet. 2003, 19: 354-356. 10.1016/S0168-9525(03)00139-2View ArticlePubMedGoogle Scholar
- Shimoda Y, Shinpo S, Kohara M, Nakamura Y, Tabata S, Sato S: A Large Scale Analysis of Protein-Protein Interactions in the Nitrogen-fixing Bacterium Mesorhizobium loti. DNA Res. 2008, 15: 13-23. 10.1093/dnares/dsm028PubMed CentralView ArticlePubMedGoogle Scholar
- Myers CL, Robson D, Wible A, Hibbs MA, Chiriac C, Theesfeld CL, Dolinski K, Troyanskaya OG: Discovery of biological networks from diverse functional genomic data. Genome Biol. 2005, 6: R114- 10.1186/gb-2005-6-13-r114PubMed CentralView ArticlePubMedGoogle Scholar
- Dale C, Moran NA: Molecular interactions between bacterial symbionts and their hosts. Cell. 2006, 126: 453-465. 10.1016/j.cell.2006.07.014View ArticlePubMedGoogle Scholar
- Pobigaylo N, Szymczak S, Nattkemper TW, Becker A: Identification of genes relevant to symbiosis and competitiveness in Sinorhizobium meliloti using signature-tagged mutants. Mol Plant Microbe Interact. 2008, 21: 219-31. 10.1094/MPMI-21-2-0219View ArticlePubMedGoogle Scholar
- Jones KM, Sharopova N, Lohar DP, Zhang JQ, VandenBosch KA, Walker GC: Differential response of the plant Medicago truncatula to its symbiont Sinorhizobium meliloti or an exopolysaccharide-deficient mutant. Proc Natl Acad Sci USA. 2008, 105: 704-9. 10.1073/pnas.0709338105PubMed CentralView ArticlePubMedGoogle Scholar
- van Noorden GE, Kerim T, Goffard N, Wiblin R, Pellerone FI, Rolfe BG, Mathesius U: Overlap of proteome changes in Medicago truncatula in response to auxin and Sinorhizobium meliloti. Plant Physiol. 2007, 144: 1115-31. 10.1104/pp.107.099978PubMed CentralView ArticlePubMedGoogle Scholar
- Jansen R, Gerstein M: Analyzing protein function on a genomic scale: the importance of gold-standard positives and negatives for network prediction. Curr Opin Microbiol. 2004, 7: 535-545. 10.1016/j.mib.2004.08.012View ArticlePubMedGoogle Scholar
- UniProt Consortium: The universal protein resource (UniProt). Nucleic Acids Res. 2008, D190-D195. 36 Database
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.View ArticlePubMedGoogle Scholar
- Janssen P, Enright AJ, Audit B, Cases I, Goldovsky L, Harte N, Kunin V, Ouzounis CA: COmplete GENome Tracking (COGENT): a flexible data environment for computational genomics. Bioinformatics. 2003, 19: 1451-1452. 10.1093/bioinformatics/btg161View ArticlePubMedGoogle Scholar
- Batagelj V, Mrvar A: Pajek – Program for large network analysis. Connections. 1998, 21: 47-57.Google 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
- Goldovsky L, Cases I, Enright AJ, Ouzounis CA: BioLayout(Java): versatile network visualisation of structural and functional relationships. Appl Bioinformatics. 2005, 4: 71-74. 10.2165/00822942-200504010-00009View ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303PubMed CentralView ArticlePubMedGoogle Scholar
- Lord PW, Stevens RD, Brass A, Goble CA: Investigating semantic similarity measures across the Gene Ontology: the relationship between sequence and annotation. Bioinformatics. 2003, 19: 1275-1283. 10.1093/bioinformatics/btg153View ArticlePubMedGoogle Scholar
- Ehrhardt DW, Atkinson EM, Long SR: Depolarization of alfalfa root hair membrane potential by Rhizobium meliloti Nod factors. Science. 1992, 256: 998-1000. 10.1126/science.10744524View ArticlePubMedGoogle Scholar