Abstract
Background
Alternative representations of biochemical networks emphasise different aspects of the data and contribute to the understanding of complex biological systems. In this study we present a variety of automated methods for visualisation of a protein-protein interaction network, using the basic helix-loop-helix (bHLH) family of transcription factors as an example.Results
Network representations that arrange nodes (proteins) according to either continuous or discrete information are investigated, revealing the existence of protein sub-families and the retention of interactions following gene duplication events. Methods of network visualisation in conjunction with a phylogenetic tree are presented, highlighting the evolutionary relationships between proteins, and clarifying the context of network hubs and interaction clusters. Finally, an optimisation technique is used to create a three-dimensional layout of the phylogenetic tree upon which the protein-protein interactions may be projected.Conclusion
We show that by incorporating secondary genomic, functional or phylogenetic information into network visualisation, it is possible to move beyond simple layout algorithms based on network topology towards more biologically meaningful representations. These new visualisations can give structure to complex networks and will greatly help in interpreting their evolutionary origins and functional implications. Three open source software packages (InterView, TVi and OptiMage) implementing our methods are available.Free full text
![Logo of bmcbioi](/http/europepmc.org/corehtml/pmc/pmcgifs/logo-bmcbioi.png)
An exploration of alternative visualisations of the basic helix-loop-helix protein interaction network
Associated Data
Abstract
Background
Alternative representations of biochemical networks emphasise different aspects of the data and contribute to the understanding of complex biological systems. In this study we present a variety of automated methods for visualisation of a protein-protein interaction network, using the basic helix-loop-helix (bHLH) family of transcription factors as an example.
Results
Network representations that arrange nodes (proteins) according to either continuous or discrete information are investigated, revealing the existence of protein sub-families and the retention of interactions following gene duplication events. Methods of network visualisation in conjunction with a phylogenetic tree are presented, highlighting the evolutionary relationships between proteins, and clarifying the context of network hubs and interaction clusters. Finally, an optimisation technique is used to create a three-dimensional layout of the phylogenetic tree upon which the protein-protein interactions may be projected.
Conclusion
We show that by incorporating secondary genomic, functional or phylogenetic information into network visualisation, it is possible to move beyond simple layout algorithms based on network topology towards more biologically meaningful representations. These new visualisations can give structure to complex networks and will greatly help in interpreting their evolutionary origins and functional implications. Three open source software packages (InterView, TVi and OptiMage) implementing our methods are available.
Background
"Graphics reveal data. Indeed graphics can be more precise and revealing than conventional statistical computations." Edward R. Tufte, The Visual Display of Quantitative Information [1].
The effective visual representation of complex data is an integral but perhaps undervalued part of a bioinformatician's job [2]. For an increasing number of researchers, this largely concerns the representation of networks, defined as sets of nodes (also called vertices) with corresponding sets of connections (undirected edges or directed arcs) between nodes. Methodologies that make the depiction of biological networks more accessible to biologists need to be developed in order to make these complex data sets as meaningful, and useful, as possible.
Biological networks
Biological networks come in many shapes and sizes. Signalling networks, food webs, metabolic pathways and gene regulation networks are examples of network data sets that are models of biological systems, often encapsulating knowledge representing many decades of experimental work [3]. Other types of network are derived from computations on genomic data, via literature mining or from the results of high-throughput experiments, and are therefore only indirectly related to the underlying biological system [4]. This latter class would include gene co-expression and co-mention networks, most protein-protein interaction (PPI) data sets and networks constructed using phylogenetic profiles [5] or gene fusion data [6].
Network layout algorithms
Network layout algorithms automatically produce visual representations of the linked nodes of a network. The aim of these algorithms is to provide easily interpretable layouts [7]. There are many aesthetic goals for these algorithms, including minimizing the number of edge crossings, minimizing the total area of the graph and maximizing symmetry [8]. Authoritative accounts of the wide variety of network layout algorithms are given in [9].
Both generic and bioinformatics-specific software are currently used for the visualisation of biological networks (Table (Table1).1). Such applications are often limited in the number of nodes and interactions that can be displayed clearly at the same time. They are capable of showing the topology of a network, but are usually devoid of meaningful biological context. Our aim here is to present some automated methods for visualisation of a protein-protein interaction network that incorporate biological information.
Table 1
Examples of network visualisation programs used in bioinformatics
Software | Description |
Pajek [37] | A visualisation and analysis application for large graphs, used primarily for social network analysis. |
GraphViz [38] | Implements a number of common graph layout algorithms. |
Otter [39] | A tool developed for visualisation of internet data. |
H3Viewer [40] | Provides layout and interactive navigation of graphs in three-dimensional hyperbolic space. |
Biolayout(Java) [41] | Higher-level biological networks (e.g. metabolic pathways and regulatory networks) may be visualised. Basic network statistics are reported. |
Osprey [42] | Builds data-rich graphical representations that are colour-coded for gene function and experimental interaction data. Web interfaces are used to retrieve up-to-date interaction data. |
Cytoscape [43] | Integrates PPI networks with microarray and other gene expression data. Allows analysis of such networks by filtering subsets of nodes or interactions. |
MAGE [44] | Used primarily in molecular modelling but has also been used to visualise networks. Produces three-dimensional images that can be zoomed and rotated in real time. |
VisANT [45] | Developed to provide interactive visual mining of biological interaction data sets. |
Java applet [46] | Displays protein interactions organised by network distance and biological function. |
Integration of biological information into network visualisation
The application of visualisation technology to network data can provide important insights into a system's structure and function [2]. In particular, integrating protein-protein interaction networks with supplementary information about the biological relationships between proteins makes it possible to display the network in a more meaningful way [10]. This extra information could be in the form of a phylogenetic tree, genomic location, known functional relationships, cellular compartments etc.
The bHLH gene family
In order to explore the alternative methods by which network data may be organised meaningfully, using biological information, a data set was needed that was rich in protein interactions and additional information, such as phylogeny. Our data set of choice was the bHLH transcription factor family that was previously studied by our group [11].
The basic helix-loop-helix (bHLH) proteins are a complex multi-gene family of transcription factors with a wide role in the developmental processes of an organism, including neurogenesis, myogenesis, and sex determination [12]. The characteristic bHLH domain is approximately 60 amino acids long and has a DNA binding region followed by two α-helices, separated by a variable-length loop. This HLH domain promotes dimerization, allowing the formation of homodimers (a complex of two identical protein molecules) or heterodimers (a complex of two different proteins) between different group members [11]. bHLH proteins are found in eukaryotic lineages but not in prokaryotes. It is assumed that the animal bHLH group expanded by gene duplications at the origin of animal multi-cellularity [12]. Based on previous studies, the mammalian bHLH proteins have been classified into five sub-families according to both phylogenetic relationships and PPI network topology [11,13,14].
In our previous analysis of the evolution of the bHLH transcription factor family [11], it became apparent that the visualisation of the protein interaction network topology alone provides limited biological insight. The integration of phylogenetic data resulted in the network becoming more ordered and biologically more meaningful. During this work, it became clear that an automated visualisation tool was needed. Here we explore several alternative types of evolutionary information – pairwise sequence diversity, discrete phylogenetic groupings, the inferred evolutionary tree itself – and discuss the biological insight that can be implied from the different types of representation.
Results
Spring-embedded layout
Figure Figure11 shows a typical layout produced by our spring-embedded viewer, using only the network topology as input. Nodes are coloured according to their bHLH sub-family (or group). This representation of the network would be similar to the output of many of the tools for network visualisation listed in Table Table1.1. Although this view emphasises the topological features of the PPI network, the nodes are not organised according to any extra biological information and in general the different colours are randomly scattered across the page.
![An external file that holds a picture, illustration, etc.
Object name is 1471-2105-8-289-1.jpg An external file that holds a picture, illustration, etc.
Object name is 1471-2105-8-289-1.jpg](/http/europepmc.org/articles/PMC1963338/bin/1471-2105-8-289-1.jpg)
Spring-embedded layout for the bHLH PPI network. Nodes represent proteins and are coloured according to sub-family (or group): E2A/A group (red), MAX/B group (green), ARNT/C group (blue), ID/G group (yellow) and HES/E group (purple). Edges show reported physical interactions between proteins taken from [11].
Spring-embedded layout incorporating evolutionary distances
When we adapt our viewer to distribute nodes according to the evolutionary distances between protein sequences, the extra constraints make it more difficult for the spring-embedding algorithm to reach a globally optimal solution without manual assistance. However, repeated visualisations using random starting positions of the nodes gave qualitatively similar clustering of the bHLH sub-families within the PPI network, confirming that the method is reliable and reproducible. By dragging nodes around the screen it is possible to explore alternative network arrangements to determine if more stable layouts may be reached and to test how well such layouts agree with secondary information such as the sub-family groups. Figure Figure22 shows a typical output for the bHLH network. As expected the ARNT and HES groups (blue and purple) form distinct clusters, whilst E2A and ID (red and yellow) appear to be closely related, and the paraphyletic group MAX (shown in green) is more dispersed.
Cluster-based layout
In many circumstances, a known classification of the nodes in a network is available that can provide important biological context to the network visualisation. In these cases, it may be appropriate to partition the nodes into discrete clusters and present these as separate groups. We implemented such a view by a further adaptation of our spring-embedded network viewer, creating a "container" for each cluster within which nodes belonging to that cluster are constrained to lie. The network layout works in the same way as the original spring-embedded viewer, so that the nodes and their containers automatically arrange themselves on the screen to produce an easily interpreted view that can be manipulated by the user. Figure Figure33 shows the output of this layout program, where the bHLH PPI network has been clustered according to the identified protein groups. Each grey "container" circle has an area proportional to the number of nodes it contains.
Phylogenetic interaction matrix
In contrast to the traditional nodes-and-edges view, a PPI network may also be represented as an interaction matrix, where the proteins are ordered in a list to reflect a phylogenetic tree and each cell represents a protein-protein interaction. This type of visualisation has been used in several recent publications [15-18]. The ATV tree viewer program [19] was modified to produce an interactive view of a phylogenetic matrix (Figure (Figure4).4). Coloured cells represent interactions between proteins belonging to the same sub-family and grey cells the interactions between proteins in different sub-families. To help explore the relationship between the PPI data and the phylogeny, the tree may be rearranged by re-rooting at a selected node or by swapping the order of the branches.
![An external file that holds a picture, illustration, etc.
Object name is 1471-2105-8-289-4.jpg An external file that holds a picture, illustration, etc.
Object name is 1471-2105-8-289-4.jpg](/http/europepmc.org/articles/PMC1963338/bin/1471-2105-8-289-4.jpg)
Phylogenetic interaction matrix for the bHLH PPI network. The inferred bHLH phylogenetic tree from [11] is shown on the left (branch lengths are not to scale). Intra-sub-family interactions are shown in the matrix to the right as filled cells in the same sub-family colour designations as figure 1; interactions between sub-families are in grey.
Tree and Arcs
To investigate another method for visualising protein interactions in relation to phylogeny, the ATV program [19] was further modified to display the protein-protein interactions as arcs against a phylogenetic tree. Figure Figure55 shows such a view of the bHLH PPI network. Grey arcs represent interactions between proteins in different sub-families; coloured arcs represent interactions within the same sub-family. The phylogenetic tree may be rearranged in the same way as for the matrix view, both to help explore the PPI network topology in relation to the tree and to minimise the number of arc crossings in the view.
Tree layout in three-dimensions
Using a force-directed optimisation method, we were able to produce tree layouts in three-dimensions upon which the protein interactions could be projected. This method is not guaranteed to find the globally optimal solution, so different random starting positions for nodes converge to different results. However, all runs produced final tree layouts that were qualitatively very similar, showing a clear separation of the bHLH sub-families in evolutionary space as shown in Figure Figure6.6. It is easy to explore the network structure in this view by rotating and zooming the layout (Additional File 1) using the KiNG three-dimensional visualisation software [20].
![An external file that holds a picture, illustration, etc.
Object name is 1471-2105-8-289-6.jpg An external file that holds a picture, illustration, etc.
Object name is 1471-2105-8-289-6.jpg](/http/europepmc.org/articles/PMC1963338/bin/1471-2105-8-289-6.jpg)
Force-directed tree layout in three-dimensions. Adding an extra dimension to the tree layout separates the sub-families according to evolutionary distance and makes their interactions easier to explore. The colours have the same designations as figure 1. For clarity the tree layout (A) and interactions (B) are shown separately. The three-dimensional Kinemage file is available as an additional file.
Discussion
The methods used in this study have been chosen to illustrate how secondary information may be applied to organise a PPI network into a biologically meaningful visualisation. In principle, this extra information could be in the form of phylogenetic profiles, genomic location, functional similarity, cellular compartment, gene co-expression or any other discrete or continuous property of proteins (or pairs of proteins). Since in this case study we are primarily interested in using visualisation methods to investigate the evolution of interactions in the bHLH gene sub-family where proteins are paralogous, we have concentrated on the use of phylogenetic information, namely an evolutionary distance matrix and a phylogenetic tree.
Evolutionary distances
Using evolutionary distance data as constraints between protein pairs within the PPI network (Figure (Figure2)2) consistently produced a meaningful clustering of the proteins with respect to their sub-families. Our evolutionary distance visualisation of the bHLH network organises the nodes into well-defined groups that correspond to the protein sub-families. This can be taken as evidence that this method is successful in producing node arrangements that are meaningful representations of the evolutionary relationships between the proteins.
Ramani and Marcotte [17] used a similar method in three-dimensional space to organise PPI networks according to distance matrix data. These networks were used to substantiate that interacting protein partners exhibit correlated evolution and therefore have similar phylogenetic trees, since proteins that cluster in distinct regions in space mirror the adjacent placement of orthologues in a phylogenetic tree.
In the bHLH interaction network, the clusters formed by the different sub-families have been shown to have distinct functionalities [11], and it is in general true that a set of proteins organised according to sequence similarity can directly reveal protein function [10]. Protein interaction networks are already being used to uncover biological roles or functional classifications for unknown proteins, for example in the popular "guilt by association" method [21]. This gives rise to the notion of using network visualisations organised by evolutionary distances as exploratory tools for the prediction of function for unannotated proteins.
Cluster-based layout
Grouping proteins according to secondary information provides an additional level of abstraction by which the PPI network may be organised (Figure (Figure3).3). Partitioning the network in such a way clarifies the visualisation of interactions both within and between groups, and may help to resolve ambiguous group assignments. Discrete clustering may also be useful in cases where the biological property of interest is not subject to continuous variation, for example Gene Ontology functional category or DNA binding motif.
This view of the bHLH network emphasises the connectivity between the different groups, showing for example that the ARNT cluster (in blue) connects only with the HES group (purple), whereas the ID proteins (yellow) connect to both the HES and E2A (red) groups. The hubs within each group can also be clearly identified, and in general the layout is clearer, with fewer edge crossings and less clutter than the original spring-embedded layout (Figure (Figure11).
Discretely partitioning the network allows us to apply many different types of secondary information to its layout. This becomes particularly important whenever biological networks are considered to have a modular structure [22]. By separating modules into different groups, the degree of modularity can be visualised objectively by the relative numbers of inter- and intra-module edges. Each module is clearly shown as a separate unit, operating in the context of the whole system.
Phylogenetic trees
Drawing protein interactions against the phylogenetic tree maps function to phylogeny, permitting explicit interaction relationships within and between sub-families to be explored. This may also provide a visual insight into the evolutionary processes responsible for the PPI network [11]. Our automated implementation juxtaposes the interaction matrix against a phylogenetic tree (Figure (Figure4).4). There are a number of features of the PPI matrix that make it an ideal tool for investigating network characteristics. For example, self-interactions – shown as filled cells on the diagonal of the matrix – are particularly clear using this representation. These interactions can be difficult to distinguish using other visualisation methods. Clusters of PPIs connecting proteins from the same sub-family are arranged near the diagonal and are coloured in the appropriate sub-family colour. Interactions that are between proteins in different sub-families are shown in grey. The relationship between the protein family phylogeny and the interaction matrix may also be explored by manipulating the tree layout. The tree may be re-rooted at a specified node, or branches swapped, thus preserving the tree topology whilst changing the order of the interaction matrix.
This phylogeny and matrix visualisation reveals hidden clusters of interactions, not apparent in other representations of the network. A further advantage of this method is that larger numbers of proteins and interactions can be mapped and visualised than in other representations, where high densities of interactions can obscure the detail of the network topology. Nodes with special characteristics can reveal themselves much more clearly in this representation than in a busy classical graph layout. Proteins acting as "hubs" are indicated by the presence of many filled cells within the same row or column, as seen for example with E2A and ARNT. A block of filled cells highlights interactions that have been retained by a group of closely-related proteins, for example the interactions of the "linker" proteins ID1 to ID4. Hongchao et al. [15] also successfully combined this technique with a clustering algorithm to visualise the PPI network in yeast. The interaction matrix representation of a PPI network is a good starting point for the exploration of very large graphs. Such visualisation is analogous to an aerial satellite map, providing a good overview of the complete network of interactions, which can then be "zoomed in" to explore the detail of specific areas.
In an alternative visualisation of the protein interactions in relation to phylogeny, the protein-protein interactions were drawn as arcs against the same phylogenetic tree (Figure (Figure5).5). As in the matrix view, the interactions may be plotted against the entire tree or against a subtree in order to show more detail. An important first insight to the PPI network may be gained with this method, revealing the context of potentially important proteins such as hubs or linkers within the network. Although the tree-and-arcs view tends to be more cluttered than the matrix view, with arcs crossing and partially obscuring each other, it has a much greater visual impact and emphasises the number of connections to each individual node. This particular visualisation has been instrumental in revealing hidden order in the topology of another protein dimerisation network, that of the bZIP transcription factors: in particular it revealed a link between redox control of DNA binding and the architecture of the network [23]. Phylogenetic tree and network visualisation have also been combined in the TreeDyn package [24].
Moving to a three-dimensional representation of the phylogenetic tree (Figure (Figure6)6) produces a clear separation of the proteins according to evolutionary distance. Although this view is probably the most information dense of all those considered in this study, use of an interactive viewer such as KiNG allows the network to be explored by zooming or rotating the view. In addition, the various components (tree, protein nodes, interactions and labels) may be displayed or hidden independently to emphasise different aspects of the data. Exploring other approaches to visualisation in three-dimensions [25,26] has the potential to be very useful for the layout of biological networks.
Insights into the evolution of the bHLH transcription factors
Several models of evolution in PPI networks have been based upon ideas of gene duplication and domain reshuffling [27-30]. Other models have assumed that an existing ancestral network is duplicated when all genes coding for interacting proteins are duplicated simultaneously [31]. This may occur during whole-genome duplication, or tandem duplications, where a fraction of the interactions between duplications may become lost.
Our automated visualisations support an interesting mechanism of evolution of the bHLH network proposed by Amoutzias et al. [11], namely that the network has evolved its sub-family structure by single domain rearrangements and then duplication of these, rather than generation of new sub-families by large-scale duplication events such as whole genome duplication. Note, this does not preclude a role for large-scale gene duplication in increasing the total number of paralogues in the network, subsequent to the duplication of the precursor bHLH sub-family members. Each sub-family is characterised by a distinct domain arrangement, and most form well-defined phylogenetic groups. The exception to this is the MAX group (green), which is paraphyletic. The five sub-families of this network are distinguishable in all of the visualisations, though the MAX group is clearly more dispersed than the other clusters in the evolutionary distance representations (Figures (Figures22 and and6)6) and forms two distinct clades in the phylogenetic tree (Figures (Figures44 and and5).5). This evolutionary relationship indicates that the MAX (or B group) most probably corresponds to the ancestral group as defined by its shared domain architecture.
The alternative representations of the network also make the evolutionary origins of its topological structures much clearer than the basic spring-embedded view shown in Figure Figure1.1. For example, repeated duplication of the ancestral MAX-binding MYC and MAD proteins has led to MAX gaining a hub-like character. Figures Figures22 and and66 show the MYC and MAD proteins as distinct groupings, all attached to MAX. The same information is shown in Figure Figure44 by two localised groups of filled green cells in the MAX row, and in Figure Figure55 as two sets of parallel arcs all linking to MAX. Several other examples of the retention of protein interactions following gene duplication can be seen in these figures.
The cluster-based view emphasises the connectivity of the five sub-families (Figure (Figure3).3). The HES group of repressors (purple) acts as a set of "bridge" proteins between the E2A (red) and ARNT (blue) sub-families. The ID repressors (yellow) mainly interact with the E2A group, but also bind HES1. Finally, the MAX group (green) is almost completely independent of the rest of the network: only USF1 and USF2 bind a protein from another sub-family (E2A).
Conclusion
Although many different programs are available for the visualisation of networks in bioinformatics, a major disadvantage of these existing tools is their lack of biological context. Producing layouts based solely on network topology gives only the most basic impression of a network's functional implications, and may actually obscure important relationships between the biological entities represented.
In this study we have demonstrated that the application of alternative network visualisation techniques can reveal different aspects of what are usually very complex data sets. The incorporation of secondary information in the form of distance metrics, trees or discrete groupings of nodes can provide insights into evolutionary processes and may help to define modules within hierarchically-structured systems. Using complementary visualisations as exploratory tools will assist in the analysis of network data sets of all sizes and types, giving us the means to put bioinformatic networks into their proper biological perspective.
Methods
bHLH sequence and interaction data
The mammalian bHLH multiple sequence alignments and PPIs were taken from our earlier work [11]. The bHLH interaction data were collated from the published literature. Multiple sequence alignments were created using CLUSTAL W [32] and gap-stripped to remove columns consisting of more than 50% gaps. For the purposes of investigating the use of phylogenetic data as constraints for network visualisation, pairwise evolutionary distances were calculated from the multiple alignment using the Jones-Taylor-Thornton substitution model for amino acid replacements per site [33]. The resulting distance matrix was used to infer a phylogenetic tree with the program BIONJ [34].
Spring-embedded layout
As an example of a classic spring-embedded network layout, a protein network viewer was developed based on the TouchGraph Java library [35]. TouchGraph attempts to optimise the network layout by minimising the lengths of edges whilst making nodes repel each other. Crucially, the user is able to manipulate the resulting layout by clicking and dragging nodes around the screen, following which the network will "relax" to its preferred shape. Network visualisations produced using TouchGraph may therefore be described as "semi-automated", in the sense that the user is able to assist the layout algorithm in producing an aesthetically pleasing result that remains consistent with the constraints imposed by spring embedding.
Spring-embedded layout incorporating evolutionary distances
The evolutionary distance matrix calculated from the bHLH sequence alignment was used as the basis for a weighted all-against-all network where each edge is weighted by the evolutionary distance between the two proteins that it connects. This network was used as a secondary input to the TouchGraph-based layout program such that each edge in the all-against-all network had an equilibrium length proportional to its weight, but only the edges present in the original PPI network were made visible. The spring-embedded layout algorithm will therefore attempt to optimise the distance between each pair of proteins so that it reflects as closely as possible their estimated evolutionary distance.
Cluster-based layout
A further method of visualisation was investigated for PPI networks that have been clustered into discrete groups using biological information. The protein viewer using the TouchGraph library was used as the basis of another spring-embedded network viewer that constrains nodes to lie within circular "containers", with area proportional to the number of nodes in the corresponding cluster. For the bHLH network being studied, proteins were clustered according to their sub-family as defined in [11].
Phylogenetic interaction matrix
Another method of visualisation that provides a direct comparison between the PPI network and the phylogenetic data is to draw the protein-protein interactions against the phylogenetic tree itself. The ATV software for the layout and manipulation of trees [19] was extended to plot an interaction matrix of the bHLH PPI network, where proteins are listed in the order corresponding to a layout of the neighbor-joining phylogenetic tree. Each filled cell within the matrix represents an interaction between two proteins. Different colours are used to differentiate interactions within or between protein sub-families.
Tree and Arcs
As an alternative to the phylogenetic matrix view, the ATV software was further adapted to allow protein-protein interactions to be directly visualised on a phylogenetic tree, represented as arcs going from source to target protein. As in the matrix view, arcs were coloured to differentiate inter- and intra-sub-family interactions.
Tree layout in three-dimensions
The visualisation of evolutionary distances in the PPI network was extended to a third dimension by using the optimization procedure described below. This finds a distribution of nodes in three-dimensional space such that their pairwise inter-node distances are as close as possible to a set of defined constraints. The constraint between each pair of nodes was calculated as the sum of branch lengths between the corresponding nodes in the tree.
The optimization procedure for inter-node distances is based on a force-directed approach, but is designed to be inherently 'noisy', including a progressive damping factor so that the layout gradually converges towards an optimal configuration whilst avoiding spurious local minima.
For a network with n nodes and a set of (positive) pairwise constraints cij, the method can be summarised as follows:
1. Start with a damping factor f = 0.5 and a random distribution of the nodes in three-dimensional space.
2. Repeat n times:
Set t = 0.
Choose a node x at random.
Calculate distances dxj from x to every other node.
Calculate the relative errors
For the node y with the greatest error exy:
Add exy2 to t.
Calculate the vector
Move x by the damped vector f
3. Calculate
4. If q is lower than any previously found value, store the current layout as the best.
5. If no improvement in q has been seen for 10 iterations, reduce f by 0.5%.
6. If no improvement in q has been seen for 100 iterations or the maximum number of iterations has been reached, report the best layout and stop.
7. Go to 2.
The resulting three-dimensional layout was output as a kinemage file [36] and visualised using the KiNG software [20].
Availability and Requirements
Project name: Network Visualisation.
Project home page: http://www.manchester.ac.uk/bioinformatics/resources. The methods described in this paper are implemented in the following open source software packages:
InterView
This software (based on TouchGraph [35]) includes implementations of the two-dimensional spring-embedded, distance-constrained spring-embedded and cluster-based layouts. Nodes and/or edges may be coloured to represent different functions or other relevant information and both undirected and directed networks are supported. Layouts produced with InterView may be exported as PDF images. InterView may be downloaded from http://www.manchester.ac.uk/bioinformatics/interview.
TVi
This software (based on ATV [19]) includes implementations of the phylogenetic tree/matrix and tree/arcs layouts. Nodes may be coloured according to defined groupings and the tree layout may be manipulated by swapping adjacent branches. Layouts produced with TVi may be exported as PDF images. TVi may be downloaded from http://www.manchester.ac.uk/bioinformatics/tvi.
OptiMage
This software produces visualisations of networks in three-dimensional space, allowing nodes to be distributed in one of three ways: according to network topology, using a phylogenetic tree (with or without associated branch lengths), or using distance matrix data. Directed or undirected interaction networks may be used, and subsets of nodes and/or edges may be coloured according to defined groupings. Output is presented in kinemage format, making it easy to share via the Web. OptiMage is available both as a command-line program and via a web server at http://www.manchester.ac.uk/bioinformatics/optimage.
Operating systems: Platform independent; Programming language: Java; Licence: GNU GPL (OptiMage), Apache-style open source (InterView/TVi).
Authors' contributions
BJH and JWP wrote the interaction matrix and tree-and-arcs visualisation software (TVi). JWP wrote the TouchGraph spring-embedded viewer (InterView, based on a prototype written by BJH) and 3D optimization software (OptiMage, based on a prototype written by SCL). GDA provided the data for the bHLH family of proteins. DLR conceived and supervised the project with the help of SCL and GDA. JWP and BJH drafted the manuscript. All authors read and approved the final manuscript.
Supplementary Material
bHLH.kin. Kinemage file containing the interactive three-dimensional view corresponding to Figure Figure6.6. This file may be viewed using the freely available KiNG software [20].
Acknowledgements
Thanks to Amelie Veron for helpful comments. JWP is funded by a BBSRC research grant (BB/C515412/1) to DLR. GDA and BJH were supported by an EPSRC platform grant (GR/R80810/01) for part of this work. BJH wishes to thank DLR, GDA and Heather Vincent for the supervision of the MSc thesis upon which this work was based.
References
- Tufte ER. The Visual Display of Quantitative Information. Cheshire, Connecticut , Graphics Press; 1983. [Google Scholar]
- Tao Y, Liu Y, Friedman C, Lussier YA. Information visualization techniques in bioinformatics during the postgenomic era. Drug Discovery Today: BioSilico. 2004;2:237–245. 10.1016/S1741-8364(04)02423-0. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Kitano H. Computational systems biology. Nature. 2002;420:206–210. 10.1038/nature01254. [Abstract] [CrossRef] [Google Scholar]
- Forst CV. Network genomics - a novel approach for the analysis of biological systems in the post-genomic era. Mol Biol Rep. 2002;29:265–280. 10.1023/A:1020437311167. [Abstract] [CrossRef] [Google 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.4285. [Europe PMC free article] [Abstract] [CrossRef] [Google 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/47056. [Abstract] [CrossRef] [Google Scholar]
- Battista GD, Eades P, Tamassia R, Tollis IG. Annotated bibliography on graph drawing. Computational Geometry: Theory and Applications. 1994;4:235–282. [Google Scholar]
- Purchase C, Carrington DA, Allder JA. Empirical evaluation of aesthetics-based graph layout. Empirical Software Engineering. 2002;7:233–255. 10.1023/A:1016344215610. [CrossRef] [Google Scholar]
- Tollis IG, Eades P, Tamassia R. Graph Drawing: Algorithms for the Visualization of Graphs. Prentice Hall; 1999. [Google Scholar]
- Adai A, Date S, Wieland S, Marcotte E. LGL: creating a map of protein function with an algorithm for visualizing very large biological networks. . J Mol Biol. 2004;340:179–190. 10.1016/j.jmb.2004.04.047. [Abstract] [CrossRef] [Google Scholar]
- Amoutzias GD, Robertson DL, Oliver SG, Bornberg-Bauer E. Convergent evolution of gene networks by single-gene duplications in higher eukaryotes. EMBO Rep. 2004;5:274–279. 10.1038/sj.embor.7400096. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Ledent V, Paquet O, Vervoort M. Phylogenetic analysis of the human basic helix-loop-helix proteins. Genome Biology. 2002;3:RESEARCH0030. 10.1186/gb-2002-3-6-research0030. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Atchley WR, Fitch WM. A natural classification of the basic helix-loop-helix class of transcription factors. Proc Natl Acad Sci USA. 1997;94:5172–5176. 10.1073/pnas.94.10.5172. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Ledent V, Vervoort M. The basic helix-loop-helix protein family: comparative genomics and phylogenetic analysis. Genome Res. 2001;11:754–770. 10.1101/gr.177001. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Hongchao L, Xiaopeng Z, Haifeng L, Geir S, Jingfen Z, Yong Z, Lun C, Yi Z, Shiwei S, Jingyi X, Dongbo B, Runsheng C. The interactome as a tree—an attempt to visualize the protein–protein interaction network in yeast. Nucleic Acids Res. 2004;32:4804–4811. 10.1093/nar/gkh814. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Newman J, Keating AE. Comprehensive identification of human bZIP interactions with coiled-coil arrays. Science. 2003;300:2097–2101. 10.1126/science.1084648. [Abstract] [CrossRef] [Google Scholar]
- Ramani A, Marcotte M. Exploiting the Co-evolution of interacting proteins to discover interaction specificity. J Mol Biol. 2003;327:273–284. 10.1016/S0022-2836(03)00114-1. [Abstract] [CrossRef] [Google Scholar]
- de Folter S, Immink RG, Kieffer M, Parenicova L, Henz SR, Weigel D, Busscher M, Kooiker M, Colombo L, Kater MM, Davies B, Angenent GC. Comprehensive interaction map of the Arabidopsis MADS Box transcription factors. Plant Cell. 2005;5:1424–1433. 10.1105/tpc.105.031831. [Abstract] [CrossRef] [Google Scholar]
- Zmasek CM, Eddy SR. ATV: display and manipulation of annotated phylogenetic trees. Bioinformatics. 2001;17:383–384. 10.1093/bioinformatics/17.4.383. [Abstract] [CrossRef] [Google Scholar]
- KiNG: Kinemage Next Generation http://kinemage.biochem.duke.edu/software/king.php
- Bader GD, Heilbut A, Andrews B, M T, T H, Boone C. Functional genomics and proteomics: charting a multidimensional map of the yeast cell. Trends in Cell Biology. 2003;13:344–356. 10.1016/S0962-8924(03)00127-2. [Abstract] [CrossRef] [Google Scholar]
- Hartwell LH, Hopfield JJ, Leibler S, Murray AW. From molecular to modular cell biology. Nature. 1999;402:C47–C52. 10.1038/35011540. [Abstract] [CrossRef] [Google Scholar]
- Amoutzias GD, Bornberg-Bauer E, Oliver SG, Robertson DL. Reduction/oxidation-phosphorylation control of DNA binding in the bZIP dimerization network. BMC Genomics. 2006;7:107. 10.1186/1471-2164-7-107. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Chevenet F, Brun C, Banuls AL, Jacq B, Christen R. TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics. 2006;7:439. 10.1186/1471-2105-7-439. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Hughes T, Hyun Y, Liberles DA. Visualising very large phylogenetic trees in three dimensional hyperbolic space. BMC Bioinformatics. 2004;5:48. 10.1186/1471-2105-5-48. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Ruths DA, Chen ES, Ellis L. Arbor 3D: an interactive environment for examining phylogenetic and taxonomic trees in multiple dimensions. Bioinformatics. 2000;16:1003–1009. 10.1093/bioinformatics/16.11.1003. [Abstract] [CrossRef] [Google Scholar]
- Wagner A. How the global structure of protein interaction networks evolves. Proc Biol Sci. 2003;270:457–466. 10.1098/rspb.2002.2269. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Amoutzias GD, Robertson DL, Bornberg-Bauer E. The evolution of protein interaction networks in regulatory proteins. Comparative and Functional Genomics. 2004;5:79–84. 10.1002/cfg.365. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Pastor-Satorras R, Smith E, Sole RV. Evolving protein interaction networks through gene duplication. J Theor Biol. 2003;222:199–210. 10.1016/S0022-5193(03)00028-6. [Abstract] [CrossRef] [Google Scholar]
- Hallinan J. Gene duplication and hierarchical modularity in intracellular interaction networks. Biosystems. 2004;74:51–62. 10.1016/j.biosystems.2004.02.004. [Abstract] [CrossRef] [Google Scholar]
- Evlampiev K, Isambert H. Evolution of protein interaction networks by whole genome duplication and domain shuffling. arXiv:q-bioMN/0606036. 2006.
- Thompson J, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–4680. 10.1093/nar/22.22.4673. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992;8:275–282. [Abstract] [Google Scholar]
- Gascuel O. BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997;14:685–695. [Abstract] [Google Scholar]
- TouchGraph: Dynamic Graph Layout http://www.touchgraph.com
- Richardson DC, Richardson JS. The kinemage: a tool for scientific communication. Protein Sci. 1992;1:3–9. [Europe PMC free article] [Abstract] [Google Scholar]
- Batagelj V, Mrvar A. Pajek - program for large network analysis. Connections. 1998. pp. 47–57.
- Gansner ER, North SC. An open graph visualization system and its applications to software engineering. Software - Practice and Experience. 2000;21:1129–1164. [Google Scholar]
- Otter: A general-purpose network visualization tool http://www.caida.org/tools/visualization/otter/
- Munzner T. Exploring large graphs in 3D hyperbolic space. IEEE Computer Graphics and Applications. 1998;18:18–23. 10.1109/38.689657. [CrossRef] [Google 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-00009. [Abstract] [CrossRef] [Google Scholar]
- Breitkreutz BJ, Stark C, Tyers M. Osprey: a network visualization system. Genome Biology. 2003;4:R22. 10.1186/gb-2003-4-3-r22. [Europe PMC free article] [Abstract] [CrossRef] [Google 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.1239303. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Freeman LC, Webster CM, Kirke DM. Exploring social structure using dynamic three-dimensional colour images. Social Networks. 1998;20:109–118. 10.1016/S0378-8733(97)00016-6. [CrossRef] [Google Scholar]
- Hu Z, Mellor J, Wu J, Yamada T, Holloway D, Delisi C. VisANT: data-integrating visual framework for biological networks and modules. Nucleic Acids Res. 2005;33:W352–W357. 10.1093/nar/gki431. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
- Mrowka R. A Java applet for visualizing protein-protein interaction. Bioinformatics. 2001;17:669–671. 10.1093/bioinformatics/17.7.669. [Abstract] [CrossRef] [Google Scholar]
Articles from BMC Bioinformatics are provided here courtesy of BMC
Full text links
Read article at publisher's site: https://doi.org/10.1186/1471-2105-8-289
Read article for free, from open access legal sources, via Unpaywall:
https://bmcbioinformatics.biomedcentral.com/counter/pdf/10.1186/1471-2105-8-289
Citations & impact
Impact metrics
Citations of article over time
Smart citations by scite.ai
Explore citation contexts and check if this article has been
supported or disputed.
https://scite.ai/reports/10.1186/1471-2105-8-289
Article citations
JNets: exploring networks by integrating annotation.
BMC Bioinformatics, 10:95, 26 Mar 2009
Cited by: 8 articles | PMID: 19323810 | PMCID: PMC2674432
Review Free full text in Europe PMC
Cataloguing the HIV type 1 human protein interaction network.
AIDS Res Hum Retroviruses, 24(12):1497-1502, 01 Dec 2008
Cited by: 82 articles | PMID: 19025396 | PMCID: PMC2655106
Review Free full text in Europe PMC
Origins of the Xylella fastidiosa prophage-like regions and their impact in genome differentiation.
PLoS One, 3(12):e4059, 31 Dec 2008
Cited by: 27 articles | PMID: 19116666 | PMCID: PMC2605562
Reconstruction of ancestral protein interaction networks for the bZIP transcription factors.
Proc Natl Acad Sci U S A, 104(51):20449-20453, 12 Dec 2007
Cited by: 26 articles | PMID: 18077348 | PMCID: PMC2154451
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
BioStudies: supplemental material and supporting data
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
2.5D visualisation of overlapping biological networks.
J Integr Bioinform, 5(1), 10 Nov 2008
Cited by: 2 articles | PMID: 20134057
Visualisation and navigation methods for typed protein-protein interaction networks.
Appl Bioinformatics, 2(3 suppl):S19-24, 01 Jan 2003
Cited by: 4 articles | PMID: 15130812
Exploration of biological network centralities with CentiBiN.
BMC Bioinformatics, 7:219, 21 Apr 2006
Cited by: 81 articles | PMID: 16630347 | PMCID: PMC1524990
Review Free full text in Europe PMC
Discerning molecular interactions: A comprehensive review on biomolecular interaction databases and network analysis tools.
Gene, 642:84-94, 10 Nov 2017
Cited by: 59 articles | PMID: 29129810
Review
Biomolecular interaction network database.
Brief Bioinform, 6(2):194-198, 01 Jun 2005
Cited by: 36 articles | PMID: 15975228
Review
Funding
Funders who supported this work.
Biotechnology and Biological Sciences Research Council (1)
Analysis and enhanced representation of protein interaction networks: towards understanding the evolution of complex intracellular systems
Professor David Robertson, The University of Manchester
Grant ID: BB/C515412/1