BACKGROUND
in recent years, large protein-protein interaction  datasets have allowed the description of relationships between proteins in complex biological systems  <cit> .
these data are commonly derived from yeast two hybrid , co-immunoprecipitation or tandem affinity purification  assays and have been obtained from a variety of unicellular and multicellular organisms  <cit> .
recent advances in high resolution mass spectrometry have further contributed to the rapid accumulation of ppi data  <cit> .
unfortunately, large scale experimental discovery of ppis remains difficult, expensive and beyond the means of many experimentalists.
currently, ppi data is almost exclusively limited to a few popular model organisms and amongst these coverage of the captured interactions is often biased to a particular domain and incomplete.
for many organisms, ppi data lags behind or is non-existent compared with genome sequence data.
in an attempt to address the relative paucity of data, a number of computational techniques have been proposed to predict and prioritise ppis  <cit> .
while the number of such methods is large, we focus here on methods that transfer functional information using cross-species orthology projection  <cit> .
in essence, the rationale is that for interacting proteins x and y in organism  we expect  that their orthologues x' and y' in organism  will also interact.
such conserved interaction pairs are called 'interologs'  <cit> .
the potential use of interolog mapping has been explored in several organisms including homo sapiens  <cit> , helicobacter pylori  <cit> , saccharomyces cerevisiae  <cit> , plasmodium falciparum  <cit>  and magnaporthe grisea  <cit> .
additionally a number of quantification methods have been developed to assess the confidence of predicted interologs  <cit> .
several web interfaces to interolog databases have been developed , but these are essentially ad hoc efforts.
they consider a small set of organisms for interolog prediction in restricted biological domains, thus hindering more widespread use.
underlying data sets are often frozen at the moment of publication or curated for a limited period of time , are dependent on other projects that are not updated  or are based on algorithms that are not state-of-the-art.
considering that both orthology projection methods and interaction data are continually updated, such static databases are destined to obsolescence.
to date the only project that provides the option to discover interologs with up to date data is openppi predictor  <cit> .
this represents a step forward, but relies on the user providing both the orthology relationships and known ppi data and only performs mappings between two species, without ranking or prioritising the putative ppi network obtained.
to address the lack of tools for performing multi-species interolog prediction 'on the fly' we created the perl module bio::homology::interologwalk.
the tool relies on bioperl  <cit> , the ensembl perl core and compara apis  <cit>  and the ebi-intact psi common query interface  enabled web service  <cit>  for its operation.
bio::homology::interologwalk is freely available under the gnu general public licence at the comprehensive perl archive network   <cit> .
bio::homology::interologwalk accepts as input a list of ensembl gene accession numbers from any of the vertebrate or metazoan genomes in ensembl and also for all species in the ensembl pan-taxonomic compara database.
the tool searches the ensembl core and compara databases and the psicquic-enabled ebi-intact ppi database to collect and analyse gene orthology and ppi data, together with ancillary information.
it then provides the option of filtering the putative interactions to retain those with strong experimental or phylogenetic support.
in addition the user can query the ppi database directly to collect all known interactions for the input gene list.
this allows the ready comparison of putative ppis from interolog projection to known ppis.
the software outputs plain text tab-separated files and can also output network representations of the ppi data and their attributes in a format compatible with the widely used biological network analysis tool cytoscape  <cit> .
we demonstrate the use of the software to investigate the potential of interolog projection on the genome of the fruit fly, drosophila melanogaster  <cit> .
the analysis  generates a novel putative ppi network that strengthens the connectivity of the known ppi network  proposes new interaction candidates.
we further calculate an interaction prioritisation index  for each of the ppis and use these to create a sub-network centred on a core of  <dig> dna replication proteins.
implementation
overview
a high-level schematic describing our implementation of the interolog walk concept is shown in figure  <dig>  the main purpose of bio::homology::interologwalk is to obtain a list of putative ppis given a set of user-selected gene identifiers in one genome of interest.
in order to be compatible with the module, the initial dataset must be a list of ensembl ids belonging to species in ensembl vertebrates, ensemblgenomes metazoa or ensembl pan-taxonomic compara databases.
concept.
schematics illustrating the principle behind interolog mapping.
proteins x and y are known to interact in a reference genome.
if they have orthologues x' and y' in the genome of interest, under certain conditions the existence of a putative x' - y' interaction can be assumed.
bio::homology::interologwalk implements this in a three-steps algorithm.
<dig>  get orthologues of genes of interest in reference genome.
queries the initial gene list against one or more ensembl databases and retrieves their orthologues.
options can be set to specify stringency of retrieved hits.
ancillary data fields are computed.
<dig>  get interactions in reference genome.
queries the orthology list built in  against psicquic-enabled ppi databases using rest.
this step will enrich the dataset built in  with the interactors of those orthologues, if any, plus ancillary data -- including parameters describing the nature and origin of the annotated interaction.
<dig>  get orthologues from reference genome back in genome of interest.
in this step the interactor list built in  is queried against one or more ensembl databases  to find orthologues back in the original genome of interest.
as in , a number of supplementary information fields are computed.
to carry out an interolog walk, bio::homology::interologwalk will first query the gene identifiers chosen by the user against the ensembl databases using the ensembl compara api  <cit> , retrieving a list of orthologous gene ids.
next, the algorithm will use the representational state transfer  interface  <cit>  to interrogate a psicquic-compliant ppi database with the list of orthologues returned by ensembl, to retrieve the list of known ppis involving them.
while there are already several interaction databases implementing the psicquic interface  <cit> , bio::homology::interologwalk currently relies on ebi intact  <cit>  as its source of experimental interactions.
having obtained a list of interactors for the orthologues of the initial gene set, in the last step of the main data mining procedure bio::homology::interologwalk will project the interactions retrieved  back to the original species of interest.
the final output is a list of putative interactors for the initial gene set and several fields of supporting data for the forward orthology map, the ppi data collection, and the backward orthology map.
these metadata fields can be analysed by a sub-module of the tool, to calculate a prioritisation index for the predicted ppis .
the procedure is organised as a pipeline of related data-processing activities.
the output of the basic pipeline can be further processed with the help of other methods in the module: it is possible to scan the results and compute counts, check for duplicate entries, isolate new gene ids  and generate cytoscape-compatible network representations of the data.
the package documentation gives a detailed account of post-processing options.
an additional stand-alone functionality of the module is the direct ppi pipeline: it is possible to use bio::homology::interologwalk to mine all the experimental ppis involving the initial gene list in the genome of interest .
this dataset is a 'snapshot' of the current experimental ppi network for the input dataset.
as such it is useful both by itself, because it tells what is currently know about the ppis involving the initial genes, and as a term of comparison for the putative ppi - because it can be used to evaluate the amount of overlap between the known and putative networks, as well as the novelty of the putative data.
additional file 1: 'simplified schematics of the interolog walk pipeline' illustrates both the putative and the experimental pipeline in detail.
one of the defining features of bio::homology::interologwalk is that the retrieval of both orthology data and protein interaction data happens on-the-fly.
the user inputs a list of gene ids plus a number of set-up parameters, and the data will be downloaded through web-service interfaces each time the program is run.
data sources
orthology predictions from ensembl compara
bio::homology::interologwalk uses the ensembl perl api http://www.ensembl.org/info/data/api.html to access the comparative biology data provided by the ensembl project through ensembl compara.
the orthology prediction method used by ensembl compara is described by vilella et al.
<cit>  and identifies several classes of homology association between genes .
ppi data from ebi-intact
bio::homology::interologwalk currently uses ebi intact  <cit>  as its source of experimental interactions.
as of june  <dig>  v.  <dig> .
<dig> of the intact database contains more than  <dig>  curated binary protein interaction evidences  <cit> .
bio::homology::interologwalk queries intact using the restful-based psicquic  implementation and retrieves data in psi-mi mitab <dig> tab-delimited format  <cit>  .
options for prioritisation of putative interactions
filtering
depending on the size of the input dataset and on the amount of information available through homology mapping, bio::homology::interologwalk can produce large numbers of putative interactions.
in such cases it might be beneficial to filter and prioritise these to generate a smaller set of putative interactions for further study.
the bio::homology::interologwalk module is composed of a number of functions that can be executed in sequence to create pipelines for retrieving interologs.
each of these functions offers options to filter the types of orthologues and interactions that are retrieved.
<dig>  spoke interactions: the user can choose whether to return any 'spoke' interactions when using interaction retrieval functions.
spoke interactions are binary interactions that are inferred from a complex of proteins that have been isolated together and as such the evidence for the interaction is indirect.
several of the most widely used protein-protein interaction data repositories including the two largest intact and biogrid explicitly draw the user's attention to the presence of spoke  interactions and provide the option of excluding them at an early stage.
<dig>  one-to-one orthology: for each of the orthology mapping functions the user can choose whether to restrict the mapping to explicit 1: <dig> relationships.
this is likely to significantly reduce the number of orthologues retrieved as the evolutionary distance between mapped species increases.
restricting mappings to direct orthologues increases the likelihood that the mapped proteins retain some common functionality.
conversely considering 1-to-many or many-to-many relationships that have arisen through duplication events risks connecting proteins and interactions whose functions have diverged  <cit> .
<dig>  experimental interactions: the user can specify whether to restrict the interactions retrieved to those that have been identified by experimental methods rather than by inference or prediction .
<dig>  physical interactions: the user can choose to retrieve only those interactions that test for direct physical association between proteins .
prioritisation
we have created an interaction prioritisation index  and a ppi conservation score  that can optionally be used to aid in the prioritisation of putative interactions.
the ipx summarises the contribution of several pieces of heterogeneous information that are collected during orthology projection and interaction retrieval.
it is not intended to be a quantitative measure of interaction reliability, but rather an integration of biological information such as orthology type, phylogenetic distance , percentage identity , interaction type and detection method .
this is similar to the work of huang et al.
<cit>  and yu et al.
<cit> .
yu et al.
used sequence similarity between the orthologous proteins to build a join similarity score, while huang et al.
proposed a scoring framework based on go functional annotation, domain information, tissue specificity and sub-cellular localisation to rank interolog-based human putative ppis obtained from six eukaryotes.
some of the indicators evaluated to create the ipx are:
•orthology type.
the kind of orthology relationship existing between an id in the genome of interest and its orthologue in the reference genome.
this feature indicates if there is a one-to-one mapping of orthologues, or if in-paralogy events in one or both sides mean we are considering a one-to-many, many-to-one or many-to-many orthologous mapping.
as explained in the filtering section, we particularly value putative ppis where both orthology relationships are of the one-to-one kind.
it has been shown  <cit>  that gene duplication is correlated with sub-functionalisation and neo-functionalisation.
when the two orthologous pairs in the interolog walk are of the one-to-one kind we set a boolean variable, Θ, to a non-negative value in the score.
we set Θ =  <dig> otherwise.
•expanded complex.
indicates whether the binary interaction has been extracted from a complex using the spoke expansion model.
a boolean non-negative term, ∑, is added to the score to reward each true binary interaction.
∑ =  <dig> for spoke-expanded binary interactions.
•opi.
overall percentage identity.
a numerical index representing the percentage identity of the conserved columns between the two orthology members' sequences.
given n total samples, we define a joint opi as the geometric mean of the two opis
•node to node distance.
a numerical indicator of the node-to-node distance in the consensus phylogenetic/species tree built by ensembl compara using genetrees  <cit>  .
we consider
where nnd <dig> is the node-to-node distance between the two orthologues in the forward projection, nnd <dig> is the node-to-node distance between the two orthologues in the backward orthology projection and we set
•interaction type & interaction detection method.
features based on psi-mi controlled vocabulary terms indicating, respectively, the type of interaction and the detection method used, within the hupo psi-mi hierarchy .
•ppi obtained with multiple methods & annotated in multiple organisms.
experimental ppis reconfirmed through the usage of further detection methods and/or observed in multiple reference genomes are acknowledged by this feature .
overall, the putative ppi interolog prioritisation index is
in this expression,
agglomerates the terms relative to the ppi in the reference organism: i is a feature scoring the interaction type and d is a feature scoring the interaction detection method.
mdm acknowledges those experimental ppis present in the database more than once, with different detection methods .
mtaxa is set to the number of reference genomes that possess an experimental interaction projecting back to the same putative ppi .
the four features are normalised to make sure their values are comparable.
normalisation constants are explained in additional file 2: 'definitions'.
the terms relative to the two orthology projections are combined in sort:
we set ωi,= ωo =  <dig>  optimisation of these two weights based on training data will allow to reward either the interaction component or the orthology component of the score to optimise performance on a case-by-case basis.
lastly, ∑ and Θ are boolean terms and we set ∑ =  <dig> whenever the putative ppi has been inferred from a binary ppi derived from a spoke-expanded complex , while Θ = n whenever the putative ppi has been inferred based exclusively on one-to-one orthology paths .
∑ and Θ are boolean flags not normalised in the ipx expression.
this is done to obtain a gross selection of putative ppi samples based on co-orthology/no co-orthology and spoke/no spoke information, prior to looking at other secondary metadata features.
the value n was chosen to be the smallest integer bigger than the maximum spread of the distribution of the normalised ipx features.
the ipx is composed of  <dig> features, f = , where  <dig> ≤ fi ≤  <dig>  ∀i ∈  <dig>  ...,  <dig> and so n =  <dig>
allowing Θ and ∑ to be one order of magnitude bigger than other ipx features means the ipx distribution will take a roughly three-modal shape, depending on the combinatorial values of ∑ and Θ, as follows:
<dig>  ∑ =  <dig>  Θ =  <dig>  - the experimental interaction is spoke-expanded and at least one of the two orthology projections is not one-to-one.
<dig>   Θ   - either the experimental interaction is spoke-expanded or at least one of the two orthology projections is not one-to-one.
<dig>  ∑ = n, Θ = n  - the experimental interaction is not expanded from a spoke-complex and the orthology projections are both one-to-one.
visual inspection of the modes in the ipx distribution can be used as strategy to filter out different sets of putative ppis, depending on the dataset considered and on the distribution of samples within the modes of the histogram.
the choice of n provides good visual separation of the modes in the ipx distribution to facilitate inspection.
we refer to the module code for further details.
the ppi conservation score  focuses on the potential for evolutionary conservation in the projected ppi by examining the density of the sub-network from which each experimental ppi is extracted.
it has been shown that the connectivity of well-conserved proteins in ppi networks is negatively correlated with their rate of evolution  <cit>  and, as a consequence, more connected proteins evolve at lower rate because they are subject to higher pressure to co-evolve with interacting proteins.
the ppi conservation score quantifies the degree of connectivity of the sub-network to which each experimental  ppi used for the interolog walk participates .
a binary interaction part of a very well-connected sub-network in the reference genome is more likely to have retained its functional characterisation after the projection to the organism of interest.
in our implementation, we follow the method suggested by huang and colleagues  <cit> , and define the ppi conservation score as
where γ =  <dig> · e/ and n and e are, respectively, the number of nodes and edges in the sub-network.
since the γ-connectedness measure is biased towards maximally connected small sub-networks, is relaxed by weighting it with the number of edges e.
schematics illustrating the ipx and pcs are shown in figure 2a-d and 2e respectively, and a detailed description of both can be found in additional file 2: 'definitions'.
there are many ways that an interolog could be prioritised.
we aim for bio::homology::interologwalk to be compatible with a diverse range of data and useful for many different kinds of users.
any prioritisation metric will be context-dependent and for this reason we offer a number of options to configure the process to suit the users requirements and the coverage and quality of the data available to them.
as such the generalised and customisable prioritisation scheme we provide here should provide the necessary exibility to allow application across a broad range of biological domains.
RESULTS
validation
retrieving known interactions through orthology walking
we tested the functionality of the bio::homology::interologwalk package by recovering known interactions using the orthologue walking principle .
to identify known interologs for the validation analyses, we obtained the complete genomes for five well-annotated species  from ensembl v.  <dig>  then, we extracted all the known experimental protein-protein associations for each of the five genomes gi  from ebi-intact.
we define  to be the set of the  experimental protein-protein interaction pairs in gi:
next, we selected five pairwise genome combinations gigj: mouse-human, human-yeast, human-fly, fly-yeast and yeast-worm.
for each gigj, we define the known positive evidence dataset  as the following subset of :
where ortho is the orthology operator.
is the set of all binary ppis in gi that match through orthology in gj .
the gene ids in the five ppi sets in  were used as input for the module.
to validate the ability of bio::homology::interologwalk to recover known interologs , we compared the degree of overlap between predicted nodes  and edges  and known positive nodes and edges, for each of the five sets .
for each venn diagram, the grey set represents the known positive set , while the white set corresponds to the algorithm's predictions.
bio::homology::interologwalk successfully retrieves 100% of the positive ppis in all cases considered.
in addition, putative ppi predictions are retrieved, which are candidates that can be prioritised for experimental evaluation.
interestingly, the known positive sets appear smaller than might be expected between closely related organisms like human and mouse.
this might be due to a combination of factors such as  the parameters for orthology classification used by ensembl are very stringent,  there are biases in experimental research across organisms   experimental ppi data will likely contain false positive interactions, which will not map through orthology.
it is also interesting to note that in the case of the yeast-worm pair  the number of novel ids and novel ppis retrieved is one order of magnitude smaller than in the other four cases.
this is consistent with the relatively limited amount of experimental ppi data available for c. elegans.
assessing the ipx using receiver operating curves
using the known positive datasets in kp from the previous section we next calculated roc curves to assess the performance of the ipx for each of the five species pairs.
for all five datasets, the area under the curve auc > <dig>  , demonstrating that there is a positive relationship between known positives and the ipx.
it is important to note, however, that this correlation varies depending on the dataset.
the reason for this is likely to be differences in the completeness of the known positive sets.
for all five datasets, the 'real' positive sets are unknown and the disparity between genome size and the number of known positives means that they are likely to represent a small proportion of the 'real' positive set.
as a consequence, the auc values are likely to underestimate the retrieval capability of the algorithm.
this also suggests that the ipx may not be optimised.
we anticipate that as coverage and categorisation of protein-protein interaction data becomes available it will be possible to optimise the ipx, improve these auc values and the utility of the metric.
the reason why a number of known positives have a low index lies in the nature of the ipx.
it is designed to reward functionally conserved interologs obtained from binary experimental ppis.
as stated above, the ipx penalises putative ppis that are from orthology projections where co-orthologues exist or from binary interactions that have been artificially extracted from protein complexes.
some known positives will fall into one or both of these two categories.
additional file 3: 'interolog prioritisation index histograms' shows ipx distributions for the five known positive datasets in kp.
additionally, we show the distribution of the known positives within the ipx histograms in additional file 4: 'distribution of positive samples within the ipx histograms'.
this chart shows, for each dataset, how many positive samples are in the low , average  and high  tiers of the ipx distribution.
for all but the mouse-human genome pair, most known positives fall in the second tier, and the mouse-human dataset is the only one to have most of its positives in the high tier.
we examined the relationship between the ipx and the loss of known positives for the five sample datasets .
the mouse-human dataset preserves 80% of the positives with an ipxthr =  <dig>  at the same threshold value, all of the other datasets fare significantly worse .
these results reflect the closer phylogenetic distance between mouse and human, in which less gene duplication will have occurred since divergence from their common ancestor in comparison to the other species pairs.
it is important to stress that the ipx is not a necessary and sufficient condition to assess the reliability of a putative ppi: a highly ranking interolog may not turn out to be an interesting candidate, but as the index is calculated using both experimental and phylogenetic measures, it would be logical to assess such interologs as candidates before those with lower values.
as is the case with any biological scoring metric, a feature with a low score may turn out to be a good candidate and should not be excluded solely on the basis of the score alone.
example — exploring the protein interactome of drosophila melanogaster
to demonstrate the use of bio::homology::interologwalk, we retrieved a list of all drosophila melanogaster genes  from ensembl-compara release  <dig>  for the interolog walk the reference genomes were not restricted to any specific species, but included all  <dig> available taxa even though we expected a proportion of the species in the sets to provide 'dead end' orthologues where no significant experimental ppi data existed.
for the interolog walk we filtered by retrieving only one-to-one orthologues in the orthology mapping phases and discarded  all interactions that were inferred from complexes  and  all interactions that were not experimental physical associations.
as a reference we also performed a direct mine of known interactions for the starting gene list with the same interactions filters.
the table shows a comparison of the ppi datasets and networks obtained for the drosophila melanogaster interactome dataset .
results obtained using the two available
bio::homology::interologwalk pipelines - putative and experimental - are shown.
in the putative pipeline, the data shown are relative to interactions obtained through interolog mapping.
in the experimental pipeline, the initial dataset has been queried against ebi intact to gather all known, experimental molecular associations available.
the field 'total interactions' indicates the total number of final entries of the form i =  obtained.
since i can be observed several times through different orthology paths, the field unique pp pairs shows the number of unique  pairs observed.
<dig>  net_ds_dmel_known  - the network consisting of all the experimentally-obtained physical associations involving genes in ds_dmel, according to ebi-intact;
<dig>  net_ds_dmel_putative  - the network consisting of putative interactions involving genes in ds_dmel according to bio::homology::interologwalk ;
<dig>  net_ds_dmel_union  - the network obtained computing the union of  and  where:
• each node is a node of net_ds_dmel_known, net_ds_dmel_putative, or both;
• each edge is either an edge of net_ds_dmel_known or an edge of net_ds_dmel_putative .
in order to explore the results of the interolog walk we analysed the networks using the network tool cytoscape  <cit> .
due to the size and complexity of the genome scale interaction networks we decided for the purposes of this illustration to focus on the  <dig> nodes in net_ds_dmel_known that were annotated with the term 'dna replication' in the gene ontology  <cit> .
in order to allow clear visualisation of the data, we further restricted this to a subset of  <dig> randomly selected genes.
we then retrieved all of their nearest neighbours in net_ds_dmel_known which produced five disconnected networks  the biggest of which features  <dig> dna replication genes .
to illustrate the utility of the interolog walk we performed the same procedure as above using net_ds_dmel_union.
in this case we retrieved a set of  <dig> 'dna replication' genes, a superset of the  <dig> found before, meaning that  <dig> additional dna replication genes are drawn in through the putative pipeline.
as before, we selected the sub-network of net_ds_dmel_union composed of the  <dig> core dna replication genes and their nearest neighbours.
the resulting sub-network, net_dnarep, composed of  <dig> nodes and  <dig> edges  has greatly increased connectivity .
indeed the main connected component in net_dnarep now comprises  <dig> genes and  <dig> interactions, and wires together  <dig> of the  <dig> core dna replication genes.
a broad analysis of net_dnarep reveals that the introduction of putative ppi data has allowed us to:
<dig>  collect evidence about new genes, not known to be related to dna replication before;
<dig>  increase the connectivity of the go-annotated dna replication proteins.
proteins that were known to be involved in dna replication now interact with proteins for which no evidence for dna replication involvement existed, meaning new potential candidates are drawn in to build a more complete picture of the domain.
using the ipx to refine the sub-network
given the relatively high number of interactions and participating genes in net_dnarep, we carried out a refinement of the interaction candidates obtained, using the ipx.
as described earlier, bio::homology::interologwalk can optionally calculate a prioritisation index for each of the putative ppis produced.
these can be employed to define a 'core' network for which there is strong biological and experimental evidence by removing nodes that are connected by putative ppis with low ipx indices.
in order to look at the composition of putative ppis in net_dnarep, we set ,  and mapped score values to edge thickness in net_dnarep .
figure 7d shows the sub-network backbone possessing the highest confidence according to the algorithm.
interestingly, the connectedness of the main component still remains at this threshold level.
four putative ppis survive the strict threshold: , ,  and .
all these pairs, apart from , are known to interact experimentally  and the putative interaction  reconfirms these results.
the  pair, on the other hand, has not been shown to interact in drosophila melanogaster and represents a candidate for lab validation.
CONCLUSIONS
in this paper, we present bio::homology::interologwalk, a perl module to retrieve, prioritise and visualise putative protein-protein interactions using interolog mapping.
unlike previous efforts, this perl library  automatically connects to orthology/ppi data web-services to generate up-to-date predictions 'on the fly'.
outputs its predictions in the form of simple text files, allowing to use its methods, or the data it produces, within the context of pipeline-based work flows of wider scope  optionally flags the predictions on the basis of related biological metadata through a prioritisation index, allowing the selection of a subset of candidates for in vivo validation.
we formally validate the efficacy of the tool and present roc curves to assess the association between the ipx and 'true-positive' interactions across several inter-species 'true-positive' sets.
we test the potential of the tool to retrieve putative ppis on the drosophila melanogaster genome and look more closely at one dna replication related sub-network identifying several novel nodes and interactions.
we conclude by using ipx thresholds to create a 'core' network for which there is strong biological and experimental support.
our interaction prioritisation index  is designed to encapsulate biologically relevant principles that relate directly to the assessments currently made manually by many researchers using interaction data.
we stress, however, that the ipx measure for an interaction is not fully explored here and that a full validation is not possible due to the current poor coverage of ppi data across species.
in our experience the ipx has proven to be a useful summary of biological metadata for ppis.
when performing an interolog walk we recommend the user first uses filtering options to restrict the interactions retrieved and then uses the ipx as a pragmatic aid to candidate prioritisation.
bio::homology::interologwalk is freely available for non-commercial purposes on the comprehensive perl archive network  and modifiable under the gnu gpl license.
the package includes full documentation and example scripts to simplify usage.
availability and requirements
project name : bio::homology::interologwalk
project home page : http://search.cpan.org/~ggallone/bio-homology-interologwalk/
programming language : perl
other requirements : bioperl, ensembl api.
the module also relies on a number of pre-requisite perl libraries.
see manual on website.
license : gnu gpl
list of abbreviations
api: application programming interface; cog: cluster of orthologs; cpan: comprehensive perl archive network; ebi: european bioinformatics institute; ipx: interolog prioritisation index; go: gene ontology; hupo-psi: human proteome organisation - proteomics standards initiative; mitab: molecular interactions tab delimited data exchange format; pcs: ppi conservation score; ppi: protein-protein interaction; psicquic: protemics standard initiative common query interface; psi-mi: proteomics standard initiative - molecular interactions; rest: representational state transfer; roc: receiver operating characteristic; tap: tandem affinity purification; treebest:  tree building guided by species tree.
competing interests
the authors declare no competing interests.
authors' contributions
gg designed and wrote the bio::homology::interologwalk module and performed all of the analyses.
apj, tis and jda supervised the work, and all authors contributed to the preparation of the final manuscript.
supplementary material
additional file 1
simplified schematics of the interolog walk pipelines.
flow diagram documenting the structure and data sources on which the bio::homology::interologwalk pipeline implementation is based.
click here for file
additional file 2
definitions.
supplementary text providing bio::homology::interologwalk implementation details, design decisions and mathematical background.
click here for file
additional file 3
interolog prioritisation index histograms.
ipx histograms for the five putative ppi datasets built from the positive datasets.
click here for file
additional file 4
distribution of positive samples within the ipx histograms.
distribution of known positive samples in the ipx histograms.
the chart shows, for each of the datasets in kp, the number of known positive samples in the low , average  and high  tiers of the ipx distribution.
click here for file
additional file 5
tpr, fpr and ipx threshold.
relationship between tpr, fpr and ipx threshold for the five putative ppi datasets obtained from the positive datasets through bio::homology::interologwalk.
click here for file
additional file 6
putative sub-network based on  <dig> core dna replication genes.
click here for file
acknowledgements and funding
funding was provided by the wellcome trust  under grant number  <dig>  by the erasysbio+ scheme  and by an mrc capacity building studentship in systems biology , number g68293x.
this work has made use of resources provided by ensembl http://www.ensembl.org and the ebi http://www.ebi.ac.uk.
gg wishes to thank lysimachos zografos for interesting discussions.

 in recent years, large protein-protein interaction datasets have allowed the description of relationships between proteins in complex biological systems. these data are commonly derived from yeast two hybrid, co-immunoprecipitation or tandem affinity purification assays and have been obtained from a variety of unicellular and multicellular organisms. recent advances in high resolution mass spectrometry have further contributed to the rapid accumulation of ppi data. unfortunately, large scale experimental discovery of ppis remains difficult, expensive and beyond the means of many experimentalists. currently, ppi data is almost exclusively limited to a few popular model organisms and amongst these coverage of the captured interactions is often biased to a particular domain and incomplete. for many organisms, ppi data lags behind or is non-existent compared with genome sequence data. in an attempt to address the relative paucity of data, a number of computational techniques have been proposed to predict and prioritise ppis. while the number of such methods is large, we focus here on methods that transfer functional information using cross-species orthology projection. in essence, the rationale is that for interacting proteins x and y in organism we expect that their orthologues x' and y' in organism will also interact. such conserved interaction pairs are called'interologs'. the potential use of interolog mapping has been explored in several organisms including homo sapiens, helicobacter pylori, saccharomyces cerevisiae, plasmodium falciparum and magnaporthe grisea. additionally a number of quantification methods have been developed to assess the confidence of predicted interologs. several web interfaces to interolog databases have been developed, but these are essentially ad hoc efforts. they consider a small set of organisms for interolog prediction in restricted biological domains, thus hindering more widespread use. underlying data sets are often frozen at the moment of publication or curated for a limited period of time, are dependent on other projects that are not updated or are based on algorithms that are not state-of-the-art. considering that both orthology projection methods and interaction data are continually updated, such static databases are destined to obsolescence. to date the only project that provides the option to discover interologs with up to date data is openppi predictor. this represents a step forward, but relies on the user providing both the orthology relationships and known ppi data and only performs mappings between two species, without ranking or prioritising the putative ppi network obtained. to address the lack of tools for performing multi-species interolog prediction'on the fly' we created the perl module bio::homology::interologwalk. the tool relies on bioperl, the ensembl perl core and compara apis and the ebi-intact psi common query interface enabled web service for its operation. bio::homology::interologwalk is freely available under the gnu general public licence at the comprehensive perl archive network. bio::homology::interologwalk accepts as input a list of ensembl gene accession numbers from any of the vertebrate or metazoan genomes in ensembl and also for all species in the ensembl pan-taxonomic compara database. the tool searches the ensembl core and compara databases and the psicquic-enabled ebi-intact ppi database to collect and analyse gene orthology and ppi data, together with ancillary information. it then provides the option of filtering the putative interactions to retain those with strong experimental or phylogenetic support. in addition the user can query the ppi database directly to collect all known interactions for the input gene list. this allows the ready comparison of putative ppis from interolog projection to known ppis. the software outputs plain text tab-separated files and can also output network representations of the ppi data and their attributes in a format compatible with the widely used biological network analysis tool cytoscape. we demonstrate the use of the software to investigate the potential of interolog projection on the genome of the fruit fly, drosophila melanogaster. the analysis generates a novel putative ppi network that strengthens the connectivity of the known ppi network proposes new interaction candidates. we further calculate an interaction prioritisation index for each of the ppis and use these to create a sub-network centred on a core of dna replication proteins. implementation overview a high-level schematic describing our implementation of the interolog walk concept is shown in figure the main purpose of bio::homology::interologwalk is to obtain a list of putative ppis given a set of user-selected gene identifiers in one genome of interest. in order to be compatible with the module, the initial dataset must be a list of ensembl ids belonging to species in ensembl vertebrates, ensemblgenomes metazoa or ensembl pan-taxonomic compara databases. concept. schematics illustrating the principle behind interolog mapping. proteins x and y are known to interact in a reference genome. if they have orthologues x' and y' in the genome of interest, under certain conditions the existence of a putative x' - y' interaction can be assumed. bio::homology::interologwalk implements this in a three-steps algorithm. get orthologues of genes of interest in reference genome. queries the initial gene list against one or more ensembl databases and retrieves their orthologues. options can be set to specify stringency of retrieved hits. ancillary data fields are computed. get interactions in reference genome. queries the orthology list built in against psicquic-enabled ppi databases using rest. this step will enrich the dataset built in with the interactors of those orthologues, if any, plus ancillary data -- including parameters describing the nature and origin of the annotated interaction. get orthologues from reference genome back in genome of interest. in this step the interactor list built in is queried against one or more ensembl databases to find orthologues back in the original genome of interest. as in, a number of supplementary information fields are computed. to carry out an interolog walk, bio::homology::interologwalk will first query the gene identifiers chosen by the user against the ensembl databases using the ensembl compara api, retrieving a list of orthologous gene ids. next, the algorithm will use the representational state transfer interface to interrogate a psicquic-compliant ppi database with the list of orthologues returned by ensembl, to retrieve the list of known ppis involving them. while there are already several interaction databases implementing the psicquic interface, bio::homology::interologwalk currently relies on ebi intact as its source of experimental interactions. having obtained a list of interactors for the orthologues of the initial gene set, in the last step of the main data mining procedure bio::homology::interologwalk will project the interactions retrieved back to the original species of interest. the final output is a list of putative interactors for the initial gene set and several fields of supporting data for the forward orthology map, the ppi data collection, and the backward orthology map. these metadata fields can be analysed by a sub-module of the tool, to calculate a prioritisation index for the predicted ppis. the procedure is organised as a pipeline of related data-processing activities. the output of the basic pipeline can be further processed with the help of other methods in the module: it is possible to scan the results and compute counts, check for duplicate entries, isolate new gene ids and generate cytoscape-compatible network representations of the data. the package documentation gives a detailed account of post-processing options. an additional stand-alone functionality of the module is the direct ppi pipeline: it is possible to use bio::homology::interologwalk to mine all the experimental ppis involving the initial gene list in the genome of interest. this dataset is a'snapshot' of the current experimental ppi network for the input dataset. as such it is useful both by itself, because it tells what is currently know about the ppis involving the initial genes, and as a term of comparison for the putative ppi - because it can be used to evaluate the amount of overlap between the known and putative networks, as well as the novelty of the putative data. additional file 1:'simplified schematics of the interolog walk pipeline' illustrates both the putative and the experimental pipeline in detail. one of the defining features of bio::homology::interologwalk is that the retrieval of both orthology data and protein interaction data happens on-the-fly. the user inputs a list of gene ids plus a number of set-up parameters, and the data will be downloaded through web-service interfaces each time the program is run. data sources orthology predictions from ensembl compara bio::homology::interologwalk uses the ensembl perl api http://www.ensembl.org/info/data/api.html to access the comparative biology data provided by the ensembl project through ensembl compara. the orthology prediction method used by ensembl compara is described by vilella et al. and identifies several classes of homology association between genes. ppi data from ebi-intact bio::homology::interologwalk currently uses ebi intact as its source of experimental interactions. as of june v.. of the intact database contains more than curated binary protein interaction evidences. bio::homology::interologwalk queries intact using the restful-based psicquic implementation and retrieves data in psi-mi mitab tab-delimited format. options for prioritisation of putative interactions filtering depending on the size of the input dataset and on the amount of information available through homology mapping, bio::homology::interologwalk can produce large numbers of putative interactions. in such cases it might be beneficial to filter and prioritise these to generate a smaller set of putative interactions for further study. the bio::homology::interologwalk module is composed of a number of functions that can be executed in sequence to create pipelines for retrieving interologs. each of these functions offers options to filter the types of orthologues and interactions that are retrieved. spoke interactions: the user can choose whether to return any'spoke' interactions when using interaction retrieval functions. spoke interactions are binary interactions that are inferred from a complex of proteins that have been isolated together and as such the evidence for the interaction is indirect. several of the most widely used protein-protein interaction data repositories including the two largest intact and biogrid explicitly draw the user's attention to the presence of spoke interactions and provide the option of excluding them at an early stage. one-to-one orthology: for each of the orthology mapping functions the user can choose whether to restrict the mapping to explicit 1: relationships. this is likely to significantly reduce the number of orthologues retrieved as the evolutionary distance between mapped species increases. restricting mappings to direct orthologues increases the likelihood that the mapped proteins retain some common functionality. conversely considering 1-to-many or many-to-many relationships that have arisen through duplication events risks connecting proteins and interactions whose functions have diverged. experimental interactions: the user can specify whether to restrict the interactions retrieved to those that have been identified by experimental methods rather than by inference or prediction. physical interactions: the user can choose to retrieve only those interactions that test for direct physical association between proteins. prioritisation we have created an interaction prioritisation index and a ppi conservation score that can optionally be used to aid in the prioritisation of putative interactions. the ipx summarises the contribution of several pieces of heterogeneous information that are collected during orthology projection and interaction retrieval. it is not intended to be a quantitative measure of interaction reliability, but rather an integration of biological information such as orthology type, phylogenetic distance, percentage identity, interaction type and detection method. this is similar to the work of huang et al. and yu et al.. yu et al. used sequence similarity between the orthologous proteins to build a join similarity score, while huang et al. proposed a scoring framework based on go functional annotation, domain information, tissue specificity and sub-cellular localisation to rank interolog-based human putative ppis obtained from six eukaryotes. some of the indicators evaluated to create the ipx are: •orthology type. the kind of orthology relationship existing between an id in the genome of interest and its orthologue in the reference genome. this feature indicates if there is a one-to-one mapping of orthologues, or if in-paralogy events in one or both sides mean we are considering a one-to-many, many-to-one or many-to-many orthologous mapping. as explained in the filtering section, we particularly value putative ppis where both orthology relationships are of the one-to-one kind. it has been shown that gene duplication is correlated with sub-functionalisation and neo-functionalisation. when the two orthologous pairs in the interolog walk are of the one-to-one kind we set a boolean variable, Θ, to a non-negative value in the score. we set Θ = otherwise. •expanded complex. indicates whether the binary interaction has been extracted from a complex using the spoke expansion model. a boolean non-negative term, ∑, is added to the score to reward each true binary interaction. ∑ = for spoke-expanded binary interactions. •opi. overall percentage identity. a numerical index representing the percentage identity of the conserved columns between the two orthology members' sequences. given n total samples, we define a joint opi as the geometric mean of the two opis •node to node distance. a numerical indicator of the node-to-node distance in the consensus phylogenetic/species tree built by ensembl compara using genetrees. we consider where nnd is the node-to-node distance between the two orthologues in the forward projection, nnd is the node-to-node distance between the two orthologues in the backward orthology projection and we set •interaction type & interaction detection method. features based on psi-mi controlled vocabulary terms indicating, respectively, the type of interaction and the detection method used, within the hupo psi-mi hierarchy. •ppi obtained with multiple methods & annotated in multiple organisms. experimental ppis reconfirmed through the usage of further detection methods and/or observed in multiple reference genomes are acknowledged by this feature. overall, the putative ppi interolog prioritisation index is in this expression, agglomerates the terms relative to the ppi in the reference organism: i is a feature scoring the interaction type and d is a feature scoring the interaction detection method. mdm acknowledges those experimental ppis present in the database more than once, with different detection methods. mtaxa is set to the number of reference genomes that possess an experimental interaction projecting back to the same putative ppi. the four features are normalised to make sure their values are comparable. normalisation constants are explained in additional file 2:'definitions'. the terms relative to the two orthology projections are combined in sort: we set ωi,= ωo = optimisation of these two weights based on training data will allow to reward either the interaction component or the orthology component of the score to optimise performance on a case-by-case basis. lastly, ∑ and Θ are boolean terms and we set ∑ = whenever the putative ppi has been inferred from a binary ppi derived from a spoke-expanded complex, while Θ = n whenever the putative ppi has been inferred based exclusively on one-to-one orthology paths. ∑ and Θ are boolean flags not normalised in the ipx expression. this is done to obtain a gross selection of putative ppi samples based on co-orthology/no co-orthology and spoke/no spoke information, prior to looking at other secondary metadata features. the value n was chosen to be the smallest integer bigger than the maximum spread of the distribution of the normalised ipx features. the ipx is composed of features, f =, where ≤ fi ≤ ∀i ∈..., and so n = allowing Θ and ∑ to be one order of magnitude bigger than other ipx features means the ipx distribution will take a roughly three-modal shape, depending on the combinatorial values of ∑ and Θ, as follows: ∑ = Θ = - the experimental interaction is spoke-expanded and at least one of the two orthology projections is not one-to-one. Θ - either the experimental interaction is spoke-expanded or at least one of the two orthology projections is not one-to-one. ∑ = n, Θ = n - the experimental interaction is not expanded from a spoke-complex and the orthology projections are both one-to-one. visual inspection of the modes in the ipx distribution can be used as strategy to filter out different sets of putative ppis, depending on the dataset considered and on the distribution of samples within the modes of the histogram. the choice of n provides good visual separation of the modes in the ipx distribution to facilitate inspection. we refer to the module code for further details. the ppi conservation score focuses on the potential for evolutionary conservation in the projected ppi by examining the density of the sub-network from which each experimental ppi is extracted. it has been shown that the connectivity of well-conserved proteins in ppi networks is negatively correlated with their rate of evolution and, as a consequence, more connected proteins evolve at lower rate because they are subject to higher pressure to co-evolve with interacting proteins. the ppi conservation score quantifies the degree of connectivity of the sub-network to which each experimental ppi used for the interolog walk participates. a binary interaction part of a very well-connected sub-network in the reference genome is more likely to have retained its functional characterisation after the projection to the organism of interest. in our implementation, we follow the method suggested by huang and colleagues, and define the ppi conservation score as where γ = · e/ and n and e are, respectively, the number of nodes and edges in the sub-network. since the γ-connectedness measure is biased towards maximally connected small sub-networks, is relaxed by weighting it with the number of edges e. schematics illustrating the ipx and pcs are shown in figure 2a-d and 2e respectively, and a detailed description of both can be found in additional file 2:'definitions'. there are many ways that an interolog could be prioritised. we aim for bio::homology::interologwalk to be compatible with a diverse range of data and useful for many different kinds of users. any prioritisation metric will be context-dependent and for this reason we offer a number of options to configure the process to suit the users requirements and the coverage and quality of the data available to them. as such the generalised and customisable prioritisation scheme we provide here should provide the necessary exibility to allow application across a broad range of biological domains.  validation retrieving known interactions through orthology walking we tested the functionality of the bio::homology::interologwalk package by recovering known interactions using the orthologue walking principle. to identify known interologs for the validation analyses, we obtained the complete genomes for five well-annotated species from ensembl v. then, we extracted all the known experimental protein-protein associations for each of the five genomes gi from ebi-intact. we define to be the set of the experimental protein-protein interaction pairs in gi: next, we selected five pairwise genome combinations gigj: mouse-human, human-yeast, human-fly, fly-yeast and yeast-worm. for each gigj, we define the known positive evidence dataset as the following subset of: where ortho is the orthology operator. is the set of all binary ppis in gi that match through orthology in gj. the gene ids in the five ppi sets in were used as input for the module. to validate the ability of bio::homology::interologwalk to recover known interologs, we compared the degree of overlap between predicted nodes and edges and known positive nodes and edges, for each of the five sets. for each venn diagram, the grey set represents the known positive set, while the white set corresponds to the algorithm's predictions. bio::homology::interologwalk successfully retrieves 100% of the positive ppis in all cases considered. in addition, putative ppi predictions are retrieved, which are candidates that can be prioritised for experimental evaluation. interestingly, the known positive sets appear smaller than might be expected between closely related organisms like human and mouse. this might be due to a combination of factors such as the parameters for orthology classification used by ensembl are very stringent, there are biases in experimental research across organisms experimental ppi data will likely contain false positive interactions, which will not map through orthology. it is also interesting to note that in the case of the yeast-worm pair the number of novel ids and novel ppis retrieved is one order of magnitude smaller than in the other four cases. this is consistent with the relatively limited amount of experimental ppi data available for c. elegans. assessing the ipx using receiver operating curves using the known positive datasets in kp from the previous section we next calculated roc curves to assess the performance of the ipx for each of the five species pairs. for all five datasets, the area under the curve auc >, demonstrating that there is a positive relationship between known positives and the ipx. it is important to note, however, that this correlation varies depending on the dataset. the reason for this is likely to be differences in the completeness of the known positive sets. for all five datasets, the'real' positive sets are unknown and the disparity between genome size and the number of known positives means that they are likely to represent a small proportion of the'real' positive set. as a consequence, the auc values are likely to underestimate the retrieval capability of the algorithm. this also suggests that the ipx may not be optimised. we anticipate that as coverage and categorisation of protein-protein interaction data becomes available it will be possible to optimise the ipx, improve these auc values and the utility of the metric. the reason why a number of known positives have a low index lies in the nature of the ipx. it is designed to reward functionally conserved interologs obtained from binary experimental ppis. as stated above, the ipx penalises putative ppis that are from orthology projections where co-orthologues exist or from binary interactions that have been artificially extracted from protein complexes. some known positives will fall into one or both of these two categories. additional file 3:'interolog prioritisation index histograms' shows ipx distributions for the five known positive datasets in kp. additionally, we show the distribution of the known positives within the ipx histograms in additional file 4:'distribution of positive samples within the ipx histograms'. this chart shows, for each dataset, how many positive samples are in the low, average and high tiers of the ipx distribution. for all but the mouse-human genome pair, most known positives fall in the second tier, and the mouse-human dataset is the only one to have most of its positives in the high tier. we examined the relationship between the ipx and the loss of known positives for the five sample datasets. the mouse-human dataset preserves 80% of the positives with an ipxthr = at the same threshold value, all of the other datasets fare significantly worse. these results reflect the closer phylogenetic distance between mouse and human, in which less gene duplication will have occurred since divergence from their common ancestor in comparison to the other species pairs. it is important to stress that the ipx is not a necessary and sufficient condition to assess the reliability of a putative ppi: a highly ranking interolog may not turn out to be an interesting candidate, but as the index is calculated using both experimental and phylogenetic measures, it would be logical to assess such interologs as candidates before those with lower values. as is the case with any biological scoring metric, a feature with a low score may turn out to be a good candidate and should not be excluded solely on the basis of the score alone. example — exploring the protein interactome of drosophila melanogaster to demonstrate the use of bio::homology::interologwalk, we retrieved a list of all drosophila melanogaster genes from ensembl-compara release for the interolog walk the reference genomes were not restricted to any specific species, but included all available taxa even though we expected a proportion of the species in the sets to provide'dead end' orthologues where no significant experimental ppi data existed. for the interolog walk we filtered by retrieving only one-to-one orthologues in the orthology mapping phases and discarded all interactions that were inferred from complexes and all interactions that were not experimental physical associations. as a reference we also performed a direct mine of known interactions for the starting gene list with the same interactions filters. the table shows a comparison of the ppi datasets and networks obtained for the drosophila melanogaster interactome dataset. results obtained using the two available bio::homology::interologwalk pipelines - putative and experimental - are shown. in the putative pipeline, the data shown are relative to interactions obtained through interolog mapping. in the experimental pipeline, the initial dataset has been queried against ebi intact to gather all known, experimental molecular associations available. the field'total interactions' indicates the total number of final entries of the form i = obtained. since i can be observed several times through different orthology paths, the field unique pp pairs shows the number of unique pairs observed. net_ds_dmel_known - the network consisting of all the experimentally-obtained physical associations involving genes in ds_dmel, according to ebi-intact; net_ds_dmel_putative - the network consisting of putative interactions involving genes in ds_dmel according to bio::homology::interologwalk; net_ds_dmel_union - the network obtained computing the union of and where: • each node is a node of net_ds_dmel_known, net_ds_dmel_putative, or both; • each edge is either an edge of net_ds_dmel_known or an edge of net_ds_dmel_putative. in order to explore the results of the interolog walk we analysed the networks using the network tool cytoscape. due to the size and complexity of the genome scale interaction networks we decided for the purposes of this illustration to focus on the nodes in net_ds_dmel_known that were annotated with the term'dna replication' in the gene ontology. in order to allow clear visualisation of the data, we further restricted this to a subset of randomly selected genes. we then retrieved all of their nearest neighbours in net_ds_dmel_known which produced five disconnected networks the biggest of which features dna replication genes. to illustrate the utility of the interolog walk we performed the same procedure as above using net_ds_dmel_union. in this case we retrieved a set of'dna replication' genes, a superset of the found before, meaning that additional dna replication genes are drawn in through the putative pipeline. as before, we selected the sub-network of net_ds_dmel_union composed of the core dna replication genes and their nearest neighbours. the resulting sub-network, net_dnarep, composed of nodes and edges has greatly increased connectivity. indeed the main connected component in net_dnarep now comprises genes and interactions, and wires together of the core dna replication genes. a broad analysis of net_dnarep reveals that the introduction of putative ppi data has allowed us to: collect evidence about new genes, not known to be related to dna replication before; increase the connectivity of the go-annotated dna replication proteins. proteins that were known to be involved in dna replication now interact with proteins for which no evidence for dna replication involvement existed, meaning new potential candidates are drawn in to build a more complete picture of the domain. using the ipx to refine the sub-network given the relatively high number of interactions and participating genes in net_dnarep, we carried out a refinement of the interaction candidates obtained, using the ipx. as described earlier, bio::homology::interologwalk can optionally calculate a prioritisation index for each of the putative ppis produced. these can be employed to define a'core' network for which there is strong biological and experimental evidence by removing nodes that are connected by putative ppis with low ipx indices. in order to look at the composition of putative ppis in net_dnarep, we set, and mapped score values to edge thickness in net_dnarep. figure 7d shows the sub-network backbone possessing the highest confidence according to the algorithm. interestingly, the connectedness of the main component still remains at this threshold level. four putative ppis survive the strict threshold:,, and. all these pairs, apart from, are known to interact experimentally and the putative interaction reconfirms these results. the pair, on the other hand, has not been shown to interact in drosophila melanogaster and represents a candidate for lab validation.  in this paper, we present bio::homology::interologwalk, a perl module to retrieve, prioritise and visualise putative protein-protein interactions using interolog mapping. unlike previous efforts, this perl library automatically connects to orthology/ppi data web-services to generate up-to-date predictions'on the fly'. outputs its predictions in the form of simple text files, allowing to use its methods, or the data it produces, within the context of pipeline-based work flows of wider scope optionally flags the predictions on the basis of related biological metadata through a prioritisation index, allowing the selection of a subset of candidates for in vivo validation. we formally validate the efficacy of the tool and present roc curves to assess the association between the ipx and'true-positive' interactions across several inter-species'true-positive' sets. we test the potential of the tool to retrieve putative ppis on the drosophila melanogaster genome and look more closely at one dna replication related sub-network identifying several novel nodes and interactions. we conclude by using ipx thresholds to create a'core' network for which there is strong biological and experimental support. our interaction prioritisation index is designed to encapsulate biologically relevant principles that relate directly to the assessments currently made manually by many researchers using interaction data. we stress, however, that the ipx measure for an interaction is not fully explored here and that a full validation is not possible due to the current poor coverage of ppi data across species. in our experience the ipx has proven to be a useful summary of biological metadata for ppis. when performing an interolog walk we recommend the user first uses filtering options to restrict the interactions retrieved and then uses the ipx as a pragmatic aid to candidate prioritisation. bio::homology::interologwalk is freely available for non-commercial purposes on the comprehensive perl archive network and modifiable under the gnu gpl license. the package includes full documentation and example scripts to simplify usage. availability and requirements project name: bio::homology::interologwalk project home page: http://search.cpan.org/~ggallone/bio-homology-interologwalk/ programming language: perl other requirements: bioperl, ensembl api. the module also relies on a number of pre-requisite perl libraries. see manual on website. license: gnu gpl list of abbreviations api: application programming interface; cog: cluster of orthologs; cpan: comprehensive perl archive network; ebi: european bioinformatics institute; ipx: interolog prioritisation index; go: gene ontology; hupo-psi: human proteome organisation - proteomics standards initiative; mitab: molecular interactions tab delimited data exchange format; pcs: ppi conservation score; ppi: protein-protein interaction; psicquic: protemics standard initiative common query interface; psi-mi: proteomics standard initiative - molecular interactions; rest: representational state transfer; roc: receiver operating characteristic; tap: tandem affinity purification; treebest: tree building guided by species tree. competing interests the authors declare no competing interests. authors' contributions gg designed and wrote the bio::homology::interologwalk module and performed all of the analyses. apj, tis and jda supervised the work, and all authors contributed to the preparation of the final manuscript. supplementary material additional file 1 simplified schematics of the interolog walk pipelines. flow diagram documenting the structure and data sources on which the bio::homology::interologwalk pipeline implementation is based. click here for file additional file 2 definitions. supplementary text providing bio::homology::interologwalk implementation details, design decisions and mathematical background. click here for file additional file 3 interolog prioritisation index histograms. ipx histograms for the five putative ppi datasets built from the positive datasets. click here for file additional file 4 distribution of positive samples within the ipx histograms. distribution of known positive samples in the ipx histograms. the chart shows, for each of the datasets in kp, the number of known positive samples in the low, average and high tiers of the ipx distribution. click here for file additional file 5 tpr, fpr and ipx threshold. relationship between tpr, fpr and ipx threshold for the five putative ppi datasets obtained from the positive datasets through bio::homology::interologwalk. click here for file additional file 6 putative sub-network based on core dna replication genes. click here for file acknowledgements and funding funding was provided by the wellcome trust under grant number by the erasysbio+ scheme and by an mrc capacity building studentship in systems biology, number g68293x. this work has made use of resources provided by ensembl http://www.ensembl.org and the ebi http://www.ebi.ac.uk. gg wishes to thank lysimachos zografos for interesting discussions.
the bio::homology::interologwalk module, sample scripts and full documentation are freely available from the comprehensive perl archive network under the gnu public license. our interolog prediction tool employs the ensembl perl api and psicquic enabled protein interaction data sources to generate up to date interologs'on-the-fly'. this represents a significant advance on previous methods for interolog prediction as it allows the use of the latest orthology and protein interaction data for all of the genomes in ensembl. we introduce bio::homology::interologwalk, a perl module to retrieve, prioritise and visualise putative protein-protein interactions through an orthology-walk method. in order to extend the utility of existing ppis, computational methods can be used that exploit functional conservation between orthologous proteins across taxa to predict putative ppis or'interologs'. the module uses orthology and experimental interaction data to generate putative ppis and optionally collates meta-data into an interaction prioritisation index that can be used to help prioritise interologs for further analysis. the module outputs simple text files, making it easy to customise the results by post-processing, allowing the putative ppi datasets to be easily integrated into existing analysis workflows. to date most interolog prediction efforts have been restricted to specific biological domains with fixed underlying data sources and there are no software tools available that provide a generalised framework for'on-the-fly' interolog prediction. we show the application of our interolog prediction method to the genomic interactome of the fruit fly, drosophila melanogaster.
 protein-protein interaction data are widely used to generate network models that aim to describe the relationships between proteins in biological systems. the fidelity and completeness of such networks is primarily limited by the paucity of protein interaction information and by the restriction of most of these data to just a few widely studied experimental organisms. in order to extend the utility of existing ppis, computational methods can be used that exploit functional conservation between orthologous proteins across taxa to predict putative ppis or'interologs'. to date most interolog prediction efforts have been restricted to specific biological domains with fixed underlying data sources and there are no software tools available that provide a generalised framework for'on-the-fly' interolog prediction.  we introduce bio::homology::interologwalk, a perl module to retrieve, prioritise and visualise putative protein-protein interactions through an orthology-walk method. the module uses orthology and experimental interaction data to generate putative ppis and optionally collates meta-data into an interaction prioritisation index that can be used to help prioritise interologs for further analysis. we show the application of our interolog prediction method to the genomic interactome of the fruit fly, drosophila melanogaster. we analyse the resulting interaction networks and show that the method proposes new interactome members and interactions that are candidates for future experimental investigation.  our interolog prediction tool employs the ensembl perl api and psicquic enabled protein interaction data sources to generate up to date interologs'on-the-fly'. this represents a significant advance on previous methods for interolog prediction as it allows the use of the latest orthology and protein interaction data for all of the genomes in ensembl. the module outputs simple text files, making it easy to customise the results by post-processing, allowing the putative ppi datasets to be easily integrated into existing analysis workflows. the bio::homology::interologwalk module, sample scripts and full documentation are freely available from the comprehensive perl archive network under the gnu public license. 