- Open Access
Improving the prediction of yeast protein function using weighted protein-protein interactions
Theoretical Biology and Medical Modellingvolume 8, Article number: 11 (2011)
Bioinformatics can be used to predict protein function, leading to an understanding of cellular activities, and equally-weighted protein-protein interactions (PPI) are normally used to predict such protein functions. The present study provides a weighting strategy for PPI to improve the prediction of protein functions. The weights are dependent on the local and global network topologies and the number of experimental verification methods. The proposed methods were applied to the yeast proteome and integrated with the neighbour counting method to predict the functions of unknown proteins.
A new technique to weight interactions in the yeast proteome is presented. The weights are related to the network topology (local and global) and the number of identified methods, and the results revealed improvement in the sensitivity and specificity of prediction in terms of cellular role and cellular locations. This method (new weights) was compared with a method that utilises interactions with the same weight and it was shown to be superior.
A new method for weighting the interactions in protein-protein interaction networks is presented. Experimental results concerning yeast proteins demonstrated that weighting interactions integrated with the neighbor counting method improved the sensitivity and specificity of prediction in terms of two functional categories: cellular role and cell locations.
Determining protein functions is an important challenge in the post-genomic era and Automated Function Prediction is currently one of the most active research fields. Previously, researchers have attempted to determine protein functions using the structure of the protein and comparing it with similar proteins. Similarities between the protein and homologues from other organisms have been investigated to predict functions. However, the diversity of homologues meant that these time-consuming methods were inaccurate. Other techniques to predict protein functions including analyzing gene expression patterns [1, 2], phylogenetic profiles [3–5], protein sequences [6, 7] and protein domains [8, 9] have been utilised, but these technologies have high error rates, leading to the use of integrated multi-sources [10, 11]. The computational approach was designed to resolve the inaccuracy of protein prediction, using information gained from physical and genetic interaction maps to predict protein functions. Recently, researchers have introduced various techniques to determine the probability of protein function prediction using information extracted from PPI. Results from these trials have been promising, but they do not address effective problems including function correlation [12–14], network topology and strength of interaction.
Network topology represents an interaction between proteins and the mechanism of that interaction. Therefore, much information can be extracted from these networks with regards to the strength of the interaction and its contribution to new function prediction, i.e. weighted contribution. A PPI network can be described as a complex system of proteins linked by interactions, and the computational analysis of PPI networks begins with the representation of the PPI network structure [15, 16]. The simplest representation takes the form of a network graph consisting of nodes and edges . Proteins are represented as nodes and two proteins that interact physically are represented as adjacent nodes connected by an edge . On the basis of this graphical representation, various computational approaches including data mining, machine learning and statistical methods can be performed to reveal the PPI networks at different levels.
The computational analysis of PPI networks is challenging and faces major problems. The first problem concerns the unreliability of protein interactions derived from large-scale experiments, which have yielded numerous false positive results (Y2H). Secondly, a protein can have more than one function and could be considered in one or more functional groups, leading to overlapping function clusters. The third problem concerns the fact that proteins with different functions may interact. Therefore, a PPI has connections between proteins in different functional groups, leading to expansion of the topological complexity of the PPI networks. Neighbour counting is a method proposed by Schwikowski et al.  to infer the functions of an un-annotated protein from the PPI. This method locates the neighbour proteins and predicts their assigned functions and the frequency of these functions; the functions are arranged in descending order according to their frequencies. The first k functions are considered and assigned to the un-annotated protein. Some papers used this technique with k equalling three. This method makes use of information from the neighbours, but it has several drawbacks: (1) it considers the interactions to be of equal weights, which is not logical; (2) it does not consider the nature of the function and whether it is dominant; (3) it does not provide a confidence level for assigning a function to the protein. The problem of confidence levels was addressed in , where the authors used chi-square statistics to calculate significance levels on the basis of the probability that various functions are present. The chi-square method provides a deeper analysis than the neighbour counting method, but it is less sensitive and specific.
Deng et al.  considered various situations for the presence of a certain function for a protein of interest: (1) number of proteins having this function; (2) number of protein pairs (interacting) having the function; (3) number of protein pairs where one has the function and the other does not; (4) number of protein pairs without this function. A weighted sum of these numbers is calculated according to the random Markov field algorithm, which assigns different weights to interactions and overcomes the above problems by considering the entire interaction network . This method considers the frequency of proteins having the function of interest and the neighbours, with less weight being placed on neighbours that are further away. Therefore, it can be used to calculate the probability that an un-annotated protein has a function of interest, and the results are more accurate than those obtained by using neighbour counting or the chi-square method.
This paper presents a new method for predicting protein function based on estimating a weight for the strength of the interaction between proteins in the PPI. The similarity between protein interactions and the connected routers in a certain autonomous number of networks was explored. Applying the idea of a network linked list of protocols such as OSPF (Open Shortest Path First) can allow information concerning surrounding routers to be obtained, according to the principles of cost and level (hop count) [22, 23]. The suggested algorithm was compared with the equal weight interactions method to indicate differences in the accuracy of prediction.
The proposed approach was applied to infer the functions of un-annotated proteins in yeast and used weighting interactions rather than free weights (equal interactions). In YPD, proteins are assigned functions based on three criteria: "Biochemical function", "Subcellular location" and "Cellular role". The numbers of annotated and un-annotated proteins, based on the three functional categories, are presented in Table 1. The accuracy of the predictions was measured by the leave-one-out method. For each annotated protein with at least one annotated interaction partner, it was assumed to be un-annotated and functions were predicted using the weighted neighbour counting method. The predicted results were compared with the annotations of the protein. Repeating the leave-one-out experiment for all such proteins allowed the specificity (SP) and sensitivity (SN) to be defined . The corresponding values of overlapped proteins for "Biochemical function", "Subcellular location" and "Cellular role" were 1145, 1129 and 1407, respectively. In the first three Figures, the relationship between sensitivity and specificity was implemented for biochemical function, cell location and cellular role, respectively. In terms of the prediction method (neighbour counting method), a fixed number of the highest frequency functions can be compared. In the present study, although one data set is used, k (number of interactions) had a variety of values (from 2 to 5). Figures 1a-d demonstrate the specificity and sensitivity in terms of biochemical function when k equals 2, 3, 4 and 5. In terms of biochemical functions (Figure 1), the sensitivity of a proposed algorithm is higher when specificity values are low. However, for higher specificity the weightless technique (W0) has good sensitivity. Therefore, an established technique is sufficient for predicting biochemical function. As demonstrated in Figures 2 and 3, the sensitivity and specificity for all weights (new suggested techniques W1-W5) were higher than W0 for all values of k. It can be demonstrated that in the cell location function category, W2 (weight relating to IG1) is the best weight to use when the number of interactions for each protein is two. W3 (weights for IG2), W1 (weights for number of experimental method) and W5 (PCA for the basic three weights (W1, W2, W3)) were the best weights when the numbers of interactions for each protein were 3, 4 or 5, respectively. Furthermore, W2 was the best weight for the cellular role function category when the number of interactions was two, and W3 (weights of IG2) were the best weights for the cellular role function category when the numbers of interactions were 3, 4 or 5. There were overlaps between some weights on the indicated curves (overlap curves), but there was a small variation in terms of detecting these weights.
The majority of methods concerning the estimation of protein functions through protein-protein interactions (PPI) use the same weights for all interactions. Such methods do not consider the various situations for each interaction including the number of experimental methods used to identify the interactions, the number of leaves connected to the interaction (whether or not the protein is sticky) and the most common graphs for the studied species within the network. Therefore, this research introduces new weights for protein interactions to enhance protein function prediction. These weights are W1-W5, and W1 depends of the number of experimental methods that identify the interaction. W1 has high confidence (100%) when the number of experimental methods used is more than one. W2 depends on the number of leaves connected to the studied interactions, which indicates whether the protein is sticky or not. The high confidence of W2 is apparent when the IG1 value is less than three (the protein is not sticky). W3 relates to the value of IG2, which indicates the global topology of the network of the studied species; its value is highly confident when IG2 is less than zero. In addition, there are two estimated weights, W4 and W5. W4 is the average of the basic weights (W1, W2 and W3), and W5 is the PCA value for the same weights. Applying the suggested weights to yeast protein functions and integrating these weights with the neighbor counting method led to enhanced results in two function categories: cell location and cellular role. The sensitivity and specificity of every point on the curves of the two function categories were higher than those obtained using the weightless technique (free or equal weights (W0)). W3 was the best weight to use in the cellular role category when the numbers of interactions were 3, 4 or 5. The cell location function category did not have a common weight for all cases but in each case (number of interactions), there was a better weight compared with other methods.
This paper introduces a novel algorithm by comparing the proteins in protein-protein interaction networks to the connected routers in the same autonomous number of networking. The protein acts as a router, and the node and edge (interaction between two proteins) act as the connection between two routers (Figures 4 and 5), where routers have up to 100 interactions (29 interactions are the maximum in the yeast proteome). As presented in Figure 4a, a group of routers and their movable messages are indicated, and the connected routers are presented in Figure 4b. In Figure 5, the group of proteins are connected using different experimental methods. The routing system can be introduced by various types of connections (LAN, WAN, Serial) as different experimental methods of interactions in the protein system. Initially, the router will be unaware of neighbour routers on the link. Therefore, the linked state protocol will be applied to the routing system, where a link is an interface on a router and the protocols are the control system of all connected routers. The protocol includes information concerning the interface's IP address/mask, the type of network (ethernet (broadcast) or serial point-to-point link), the cost of that link and any neighbour routers on that link. In the protein system, a generic protocol is followed that identifies the protein by name (gene name, locus name, accession name etc...), ID (determined number for each protein), sequence (amino acids in given number and order) and functions (if known). The type of network will be elucidated; interaction between two proteins (protein pair) or dense interactions (cluster), and the weight of the interaction (our contribution). Furthermore, neighbours of the adjacent protein (known interactions in the network) are identified (Table 2). The protein interactions are calculated until the second level. The algorithm is performed following four steps: (1)- determining the level and degree for each adjacent protein, (2)- calculating the weight (cost) for each interaction (interaction with high cost/weight is strong), (3)- integrating these data to predict the function of the un-annotated proteins using the neighbourhood counting method, and (4)- calculating the sensitivity and specificity for the different weights.
There is a difference between the degree and the level of any node. The degree of a node (protein) is defined as the total number of connected nodes or proteins directly surrounding this node (protein A has degree equal to six) as shown in Figure 6; the level of a node is the layer of nodes related to the main one. The directed nodes have a level equal to one, and their neighbours are the second level as presented in Figure 6. The red nodes are the first level of protein A (black), the second level of proteins are the yellow coloured nodes (nodes connected to protein's A neighbours). The last (third) level is the group of proteins coloured green. In router networks, the hop count principle is performed to determine the router level. In this paper, the second level was assumed to be sufficient for extracting the most important information about the function of a protein. The concept of node level was applied to 2559 protein-protein interactions between 6416 proteins collected from the Munich Information center of Protein Sequences (MIPS, http://mips.gsf.de) for the yeast Saccharomyces cerevisiae . As demonstrated in Figure 7, proteins with ID numbers 1913, 3246 and 3517 had a level equal to one for the studied protein number 1, and the yellow nodes are second degree.
PPI weight calculation
Protein-protein interaction weights are introduced and each interaction has a specific weight. Three basic methods were considered in terms of calculating the weights of all the interactions and overcoming problems affecting the interaction network. The first method concerns the number of experimental methods. Protein-protein interactions are identified by high-throughput experimental methods such as Y2H [25–29], mass spectrometry of co-immunoprecipitated protein complexes (Co-IP) [30, 31], gene co-expression, TAP purification cross link, co-purification and biochemical methods. Challenging technical problems arise using the first two methods, which lead to spurious interactions due to self activation in Y2H and abundant contaminants with CO-IP. These problems lead to false positive interactions . Therefore, a quantitative method for evaluating the pathway through proteomics data is required. A number of experimental and computational approaches have been implemented for large-scale mapping of PPIs to realize the potential of protein networks for systems analysis. One method utilizes multiple independent sets of training positives to reduce the potential bias of using a single training set; this method uses association with publishing identifiers or foundation in two or more species, otherwise PPIs must have an expression correlation more than 0.6 . Another technique also obtains conserved patterns of protein interactions in multiple species . There are several methods for determining the reliability of interactions [35–38]. In this paper, the reliability or confidence is introduced by counting the number of experimental methods for each interaction; some interactions have been identified using many experimental methods and others identified by just one. In yeast proteins, approximately ten experimental methods can be used to identify protein-protein interactions (Edge between Protein (YBR0904) and Protein (YDR356W) can be identified by ten experimental methods where protein (AAC1) and protein (YHR005C-A) can be identified by one method). As demonstrated in Figure 8, approximately 750 interactions from 2559 proteins have been identified by more than one experimental method. More than half of all the interactions have been identified by just one method (~1800 interactions); researchers have high confidence (100%) concerning those interactions identified by more than one method and 50% confidence for the others (one method identification). Table 3 presents the yeast protein interactions, the number of experimental methods used to identify them and the identification value for each one. This method does not depend on clear points on computational algorithms, but reflects the strength of interaction from the laboratory viewpoint. Another approach for estimating the reliability of experimental methods concerns calculating the stability of every method.
The second method for calculating weights of interactions is the IG1 concept (Interaction Generality 1) [39–41]. A new method for assessing the reliability of protein-protein interactions (local topology) is obtained in biological experiments by calculating the number of proteins involved in a given interaction (number of protein leaves connecting to the two studied proteins incremented by one) as shown in Figure 9. IG1 assumes that complicated interaction networks are likely to be true positives. By implementing the IG1 on the collected data (yeast protein interactions), the range of IG1 was between one and 21 (Figure 10), meaning that some interactions have many leaves. According to the IG1 concept, increasing values leads to false positive interactions. In the suggested algorithm, it is assumed that interactions with an IG1 value less than four (as threshold) have high confidence (100%) and those with more than four have low confidence (Table 4). For example, the interaction between proteins YMR056C and YHRS01C has an IG1 value of three (weight = 100%) when the interaction between proteins YMR056C and YDR167W has an IG1 value of four (weight = 50%). However, the interaction between proteins YDL043C and YMR117C has an IG1 value of 21 (weight = 50%).
The third method for calculating the weight uses the IG2 concept (Interaction Generality 2), [42, 43]. This algorithm explores the five major sub-graphs of a network to obtain information concerning the global topology of the network. After collecting the five values for each interaction according to Figure 11, principal component analysis (PCA) has been implemented. The PCA concept for the previous major topologies of yeast protein networks was implemented and IG2 values ranged from -281 up to ~27 (Table 5). By determining the threshold (19) as the margin of reliability, it is assumed that IG2 values less than the threshold are more accurate than those above the threshold. Regarding the three previous methods for calculating weights, high confidence interactions can be collected compared with low confidence interactions (Figure 12). After collecting the weights from the three previous methods (number of experimental methods, IG1 and IG2), new weights strategies can be created using an average of the three values or PCA. Five different weights for each interaction were collected. As indicated in Table 6, interactions between proteins AAC1 and YHR005C-A had a W1 = 0.5, which means that only one method was used to identify it; W2 = 1, therefore it has more than three leaves in IG1 (IG1 < 4), W3 = 0.5 indicating that IG2 was more than 19, W4 is the average of the three weights which equalled 0.66 (1/3 Σ wi, i = 1..3), and W5 (PCA of the three weights with threshold equal zero) was 0.5, indicating that its value is more than zero. This example demonstrates a weak interaction (edge) between protein ID 1 (AAC1) and protein ID 1913 (YHR005C-A). Another example concerning high confidence is shown in the second row and concerns a protein interaction (edge) between ANC1 and SNF5, where the weights are 1, 1, 0.5, 0.83 and 1 for W1, W2, W3, W4 and W5, respectively. Relating to the main three measurements, many weights can be created by applying AND/OR processes. However, each weight can be multiplied by the coefficient relating to the importance of its role in determining the edge; 0.35, - 0.2 and - 0.4 for W1, W2 and W3, respectively. The integration was performed on the five weights explored (W1-W5). The neighbour counting method was applied six times, once for the basic weight (equal weights or traditional method) and once for each of the five estimated weights.
After collecting the levels for each protein and the five different weights, the neighbour counting method (frequencies of interaction partners having certain functions of interest) was implemented to predict the functions of unknown proteins. The new weights and weightless (edges with equal weights (W0)) algorithms were compared for proteins having up to five interactions. This demonstrated that for most selected new weights at a specific specificity (SP), the sensitivity (SN) was higher than obtained using W0. As demonstrated in Figures 1, 2 and 3, the sensitivity and the specificity of the three function categories of yeast (Biochemical, Cell location and Cellular role) are indicated respectively. Equations 1 and 2 present the basic laws of SN and SP.
The sensitivity (SN) and specificity (SP) can be defined as:
Where n i is the number of observed functions for protein P i
m i is the number of predicted functions for protein P i ,
and k i are the overlaps between them.
Zhao XM, Wang Y, Chen L, Aihara K: Gene function prediction using labeled and unlabeled data. BMC Bioinformatics. 2008, 9: 57-71. 10.1186/1471-2105-9-57.
Zhao H, Wu B: DNA-Protein Binding and gene expression patterns. Lecture Notes-Monograph Series. 2003, Statistics and Science: A Festschrift for Terry Speed, 40: 259-274.
Morin M: Phylogenetic Networks Simulation, Characterization, and Reconstruction. New Mexico. 2007
Sun J, Zhao Z: Construction of phylogenetic profiles based on the genetic distance of hundreds of genomes. Biochem Biophys Res Commun. 2007, 355 (3): 849-853. 10.1016/j.bbrc.2007.02.048.
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 (8): 4285-4288. 10.1073/pnas.96.8.4285.
Harrington E, Singh A: Quantitative assessment of protein function prediction from meta genomics shotgun sequences. PNAS. 2007, 104 (35): 13913-13918. 10.1073/pnas.0702636104.
Spriggs R, Murakami Y: Protein function annotation from sequence: prediction of residues interacting with RNA. Bioinformatics. 2009, 25 (12): 1492-1497. 10.1093/bioinformatics/btp257.
Friedberg I: Automated protein function prediction---the genomic challenge. Bioinformatics. 2006, 7 (3): 225-242.
Nariai N, Kolaczyk ED, Kasif S: Probabilistic Protein Function Prediction from Heterogeneous Genome-Wide Data. PLoS ONE. 2007, 2: e337-10.1371/journal.pone.0000337.
Gupta R GT: Comparative study of various genomic data sets for protein function prediction and enhancements using association analysis. the 7th SIAM International Conference on Data Mining; 28-30 April; Minneapolis. 2007
Liu Y, Kim I, Zhao H: Protein interaction predictions from diverse sources. Drug Discov Today. 2008, 13 (9-10): 409-416. 10.1016/j.drudis.2008.01.005.
Khaled S, Nahed S, Yasser K: Estimation of the correlation between protein sub-function categories based on overlapping proteins. 27th National Radio Science conference (NRSC). 2010, 323-331.
Khaled S, Nahed S, Yasser K: Exploring Protein Functions Correlation Based On Overlapping Proteins and Cluster Interactions. 1st middle east conference for biomedical engineering; 21-24 Feb; Sharqah. 2011
Khaled S, Nahed S, Yasser K: Determining The Relations Between Protein Sub-Function Categories Based On Overlapping Proteins. JCC. 2011,
Sharan R, Ulitsky I, Shamir R: Network-based prediction of protein function. Mol Syst Biol. 2007, 3-88.
Khaled S, Nahed S, Yasser K: Comparison between different methods for protein function prediction. 1st International Joint Conference (NRC), Cairo. 2009, 116-117.
Wagner A: How the global structure of protein interaction networks evolves. Proc Biol Sci. 2003, 270 (1514): 457-466. 10.1098/rspb.2002.2269.
Aidong Z: PROTEIN INTERACTION NETWORKS: Computational Analysis. 2009, New York, Press
Schwikowski B, Uetz P, Fields S: A network of protein-protein interactions in yeast. Nat Biotechnol. 2000, 18 (12): 1257-1261. 10.1038/82360.
Hishigaki H, Nakai K, Ono T, Tanigami A, Takagi T: Assessment of prediction accuracy of protein function from protein--protein interaction data. Yeast. 2001, 18 (6): 523-531. 10.1002/yea.706.
Deng M, Zhang K, Mehta S, Chen T, Sun F: Prediction of protein function using protein-protein interaction data. J Comput Biol. 2003, 10 (6): 947-960. 10.1089/106652703322756168.
Edde N: CCNA: Cisco certified network associate study guide. 2007, Canada press
Cisco "CCNA ICND2 Official Exam Certification Guide". 2007, USA press
The yeast protein interactions database. [http://mips.gsf.de]
Dobrin N, Randall W: New and Emerging Proteomic Technique. 2006, New Jersey Press
McDonald PN: Two-hybrid systems. Methods and protocols. Introduction. Methods Mol Biol. 2001, 177: v-viii.
Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, Lockshon D, Narayan V, Srinivasan M, Pochart P: A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. Nature. 2000, 403 (6770): 623-627. 10.1038/35001009.
Ito T, Chiba T, Ozawa R, Yoshida M, Hattori M, Sakaki Y: A comprehensive two-hybrid analysis to explore the yeast protein interactome. Proc Natl Acad Sci USA. 2001, 98 (8): 4569-4574. 10.1073/pnas.061034498.
Tong AH, Drees B, Nardelli G, Bader GD, Brannetti B, Castagnoli L, Evangelista M, Ferracuti S, Nelson B, Paoluzi S: A combined experimental and computational strategy to define protein interaction networks for peptide recognition modules. Science. 2002, 295 (5553): 321-324. 10.1126/science.1064987.
Gavin AC, Bosche M, Krause R, Grandi P, Marzioch M, Bauer A, Schultz J, Rick JM, Michon AM, Cruciat C: Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature. 2002, 415 (6868): 141-147. 10.1038/415141a.
Ho Y, Gruhler A, Heilbut A, Bader GD, Moore L, Adams SL, Millar A, Taylor P, Bennett K, Boutilier K: Systematic identification of protein complexes in Saccharomyces cerevisiae by mass spectrometry. Nature. 2002, 415 (6868): 180-183. 10.1038/415180a.
Mrowka R, Patzak A, Herzel H: Is there a bias in proteome research?. Genome Res. 2001, 11 (12): 1971-1973. 10.1101/gr.206701.
Yu J, Finley RL: Combining multiple positive training sets to generate confidence scores for protein-protein interactions. Bioinformatics. 2009, 25 (1): 105-111.
Sharan R, Suthram S, Kelley RM, Kuhn T, McCuine S, Uetz P, Sittler T, Karp RM, Ideker T: Conserved patterns of protein interaction in multiple species. Proc Natl Acad Sci USA. 2005, 102 (6): 1974-1979. 10.1073/pnas.0409522102.
Okada K: Flow model of the protein-protein interaction network for finding credible interactions. 2007, 317-326.
Bader JS, Chaudhuri A, Rothberg JM, Chant J: Gaining confidence in high-throughput protein interaction networks. Nat Biotechnol. 2004, 22 (1): 78-85. 10.1038/nbt924.
Aytuna AS, Gursoy A, Keskin O: Prediction of protein-protein interactions by combining structure and sequence conservation in protein interfaces. Bioinformatics. 2005, 21 (12): 2850-2855. 10.1093/bioinformatics/bti443.
Deng M, Sun F, Chen T: Assessment of the reliability of protein-protein interactions and protein function prediction. Pac Symp Biocomput. 2003, 140-151.
Saito R, Suzuki H, Hayashizaki Y: Interaction generality, a measurement to assess the reliability of a protein-protein interaction. Nucleic Acids Res. 2002, 30 (5): 1163-1168. 10.1093/nar/30.5.1163.
Nabieva E, Jim K, Agarwal A, Chazelle B, Singh M: Whole-proteome prediction of protein function via graph-theoretic analysis of interaction maps. Bioinformatics. 2005, 21 (Suppl 1): i302-310.
Chua HN, Sung WK, Wong L: Exploiting indirect neighbours and topological weight to predict protein function from protein-protein interactions. Bioinformatics. 2006, 22 (13): 1630-1623.
Saito R, Tomita M: Interaction Generality, a Measurement to Assess the Reliability of a Protein-Protein Interaction. Genome Informatics. 2002, 13: 324-325.
Saito R, Suzuki H, Hayashizaki Y: Construction of reliable protein-protein interaction networks with a new interaction generality measure. Bioinformatics. 2003, 19 (6): 756-763. 10.1093/bioinformatics/btg070.
The authors declare that they have no competing interests.
KS collected the biological data, participated data implementation and drafted the manuscript. NS participated in the statistical analysis. YK conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.
Khaled S Ahmed, Nahed H Saloma contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.