Abstract
High-throughput detection of protein interactions has had a major impact in our understanding of the intricate molecular machinery underlying the living cell, and has permitted the construction of very large protein interactomes. The protein networks that are currently available are incomplete and a significant percentage of their interactions are false positives. Fortunately, the structural properties observed in good quality social or technological networks are also present in biological systems. This has encouraged the development of tools, to improve the reliability of protein networks and predict new interactions based merely on the topological characteristics of their components. Since diseases are rarely caused by the malfunction of a single protein, having a more complete and reliable interactome is crucial in order to identify groups of inter-related proteins involved in disease etiology. These system components can then be targeted with minimal collateral damage. In this article, an important number of network mining tools is reviewed, together with resources from which reliable protein interactomes can be constructed. In addition to the review, a few representative examples of how molecular and clinical data can be integrated to deepen our understanding of pathogenesis are discussed.
1. Introduction
The existence of living cells is not possible without organized and coordinated communication between proteins. Failure of the control mechanisms that underlie these delicate relationships can lead to disease or even death (Lesk, ). This highlights that the study of the complex network of interactions between proteins is crucial to improve our understanding of the intricate mechanisms that make life possible (Lesk, ; Loscalzo and Barabasi, ). To ease the analysis of this involved biological machine, it is commonly represented as a network of nodes, linked to each other if there is evidence of their physical or functional relationship.
Today we have access to vast Protein Interaction Networks (PINs) from different organisms, due to high-throughput experimental techniques that are often an improved variation of yeast-two-hybrid screenings, or of co-immunoprecipitation followed by mass spectrometry (Vidal et al., 2011). Nevertheless, these networks are incomplete and contain a significant number of false positive interactions (Kuchaiev et al., ). However, it is fortunate that their structural properties are not different from those observed in good quality social or technological networks (Albert and Barabási, ; Liu et al., ; Cannistraci et al., ) (Figure 1A). These topological similarities have prompted the development of tools, based on node-connectivity properties, aimed at improving the reliability and completeness of complex networks (Cannistraci et al., ).
Figure 1
The reliability indices and predictions resulting from the application of these methods can be integrated with other sources of high-quality protein interactions (PIs). With these, one can construct reliable PINs that, together with clinical and genetic data, represent the fundamental pieces of information used in the emerging field of network medicine (Barabási et al.,
2. Topological reliability assessment and prediction of protein interactions
The observable network topologies of biological systems are not complete and contain spurious interactions. In addition, the mechanisms that lead two proteins to interact are not fully understood yet. As a consequence, traditional machine learning algorithms cannot be easily applied to PINs. Not only is the definition of features to discriminate between interacting and non-interacting proteins a challenging task, but also the construction of positive and negative sets of interactions to train these algorithms. For example, two unlinked proteins in the observable network cannot be considered as part of the negative set: it could very well be that they are disconnected due to experimental constraints that prevented scientists from observing their interaction. Alternatively it could be that, two linked proteins represent a false positive that is part of the dataset because one of the interactors is, for example, a sticky protein (Saito et al.,
In this context, the assignment of likelihood scores to connected and disconnected pairs of proteins, on the mere basis of the observable network topology, is a convenient means to improve the degree of confidence and completeness of PINs (Cannistraci et al.,
2.1. Neighborhood-based techniques
In 2001, Newman found that the relative probability of collaboration between scientists increases with their number of common acquaintances (Newman,
Newman's findings triggered the development of a myriad of neighborhood-based approaches (Lü and Zhou,
In 2013, Cannistraci and colleagues introduced a paradigm shift in topological link prediction, by noting that the presence of a tightly connected set of CNs increases the probability of interaction between non-adjacent nodes (Cannistraci et al.,
Although the above mentioned techniques can be applied to PINs, they were formulated for networks in general and do not consider any particular biological assumption. The pioneers of PI reliability assessment and prediction are Saito and colleagues. In 2002, after observing that the partners of sticky proteins and self-activators do not interact with anything else in PINs, they proposed the Interaction Generality index (IG1), which assigns low reliability scores to protein pairs whose neighbors have very few partners (Saito et al.,
Another two indices put forward in the context of protein interactomes are the Interaction Reliability by Alternative Paths index (IRAP) and its successor IRAP* (Chen et al.,
All the afore-mentioned techniques represent, in general, an efficient and accurate way to identify protein pairs that are good candidates for interaction (see the formulation of some of them and their application to a toy example in Figure 1C). However, they all strongly depend on topological information to work properly. As a consequence, they perform poorly when applied to very sparse networks, like the PINs of non-model or poorly annotated organisms (You et al., 2010).
2.2. Maximum likelihood techniques
Maximum likelihood approaches, introduced mainly for link prediction, rely on the underlying community structure of complex networks. In the Hierarchical Random Graph (Clauset et al.,
2.3. Network embedding techniques
Data analysts are regularly faced with the problem of finding meaningful low-dimensional representations of high-dimensional data. Algorithms such as Multidimensional Scaling or Principal Component Analysis embed data to low dimensions by preserving inter-sample distances or covariances but, if the dataset under study contains non-linear structure, they fail to provide useful mappings (Tenenbaum et al., 2000). To solve this issue, non-linear dimensionality reduction algorithms, such as Isometric Feature Mapping (ISOMAP), are commonly employed. Under the hypothesis that the biological features that lead to a PI are complex and non-linear, one could assume that PINs are shaped over a manifold embedded in a high-dimensional space, where interacting proteins are geometrically close to each other and disconnected pairs are far apart (Kuchaiev et al.,
Nataša Pržulj and her colleagues are pioneers in the modeling of PINs with geometric graphs. Their computational experiments show close matches between important topological properties of PINs and geometric random graphs (Przulj et al.,
Around the same time period, a group of physicists and network scientists were independently developing a framework to model complex networks, resting on the assumption that a hidden metric space underlies them and shapes their topology (Boguñá et al.,
Finally, in the non-centered Minimum Curvilinear Embedding (ncMCE), a technique that has been successfully applied in different fields (Alanis-Lobato et al.,
2.4. General framework for measuring the effectiveness of these techniques
In order to benchmark the accuracy of a link prediction technique, the following framework is commonly employed:
Remove L randomly selected PIs from the observable network topology.
Assign confidence scores to disconnected protein pairs in the pruned network with a topological technique and sort them decreasingly (best candidate interactions positioned at the top of this list).
Take L protein pairs from the top of the sorted list and compute the proportion present in the set of interactions removed in 1. This is a measure of the technique's precision.
Repeat steps 1–3 t times, removing different sets of randomly selected PIs.
Repeat steps 1–4 removing 2L, 3L, etc. interactions, up to the point where the network loses connectivity. This allows for the construction of a sparsification curve (SC), whose points are the mean precisions of the technique applied at each sparsification level.
This evaluation depicts the ability of a topological approach to predict accurately under the presence of less and less network information. Nonetheless, it has an intrinsic problem because, as discussed above, some of the candidate interactions with high confidence scores may not be part of the randomly removed set of PIs. However, they may represent good candidates that current technologies cannot measure. Moreover, members of the removed set of links may be false positives that good link predictors are correctly discarding by giving them low scores. Subsequently, researchers have opted for using Gene Ontology (GO) similarities (Yu et al., 2010) to discriminate between good and bad candidate PIs. This is based on the guilt-by-association principle (Oliver,
Figure 1E presents the minimum area under the SC for most of the topological techniques described in this section, when they are applied to four yeast networks for the link prediction task (Cannistraci et al.,
3. Resources for high confidence protein interactions
Proteins with a high likelihood to interact can considerably reduce the universe of possible pairs to test in the lab and guide wet-lab validations. These interactions can then be integrated with available repositories of high-quality PIs that attach confidence scores to each reported interaction (see Figure 1F). One of such resources is the Search Tool for the Retrieval of Interacting Genes (STRING), which provides a combined score that indicates higher confidence when more than one source of evidence supports an interaction (Szklarczyk et al., 2011). STRING evidence sources include computational associations (neighborhood-based, co-occurrence, co-expression, text mining), high-throughput experiments, other databases, and interactions identified in other organisms. The current version of STRING (available at http://string-db.org) provides an interactive network viewer and access to interactions between almost 10 million proteins, from more than 2000 organisms (Szklarczyk et al., 2015).
The Human Integrated Protein-Protein Interaction rEference (HIPPIE) retrieves interactions from major expert-curated databases and calculates a score for each PI, reflecting its combined experimental evidence. This score is a function of the number of studies supporting the interaction, the quality of the experimental techniques used to measure it and the number of organisms in which it is present (Schaefer et al.,
Another worth-mentioning resource is INstruct (http://instruct.yulab.org/). It collects interactions from eight major expert-curated databases and filters out low-quality PIs, to keep only those supported by domain-domain interactions obtained from co-crystal structures (Wang et al., 2012; Meyer et al.,
It is important to stress that when querying interactions from these resources, high-confidence should be preferred over size. In a recent study, Rolland and colleagues assembled PIs from 7 public databases and found that interactions supported by multiple sources can be validated at rates that are significantly higher than those of PIs supported by a single method (Rolland et al.,
4. Protein interaction networks in health and disease
It is possible that the first work that advocated for a systems-based approach to disease is the one by Goh et al. (
Figure 2

(A) A bipartite network of diseases and their associated genes or symptoms can be mapped to the disease or gene/symptom space by linking nodes of one type that are connected with the same nodes of the other. The weight of the edges in the resulting projection indicates the number of such common nodes. (B) The application of a community detection algorithm to the Autoimmune Disease Network, mapped to the gene space, reveals groups of genes associated with similar disorders and high levels of co-morbidity (adapted from Alanis-Lobato et al.,
A similar analysis, focused on the gene projection, was performed considering only autoimmune diseases (Alanis-Lobato et al.,
Given a set of proteins associated with a patient's phenotype, Lage and co-workers are able to rank disease-causing proteins as the top candidates with the help of a phenotype similarity score. This also allows them to identify previously unknown disease-causing complexes (Lage et al.,
In 2014, Zhou and colleagues extracted disease and symptom terms from the Medical Subject Headings (MeSH) in PubMed and linked diseases with symptoms via bibliographic records (Figure 2A middle). Instead of simply mapping this network to the disease space, they describe each disease with a vector of symptoms, with entries quantifying the strength of association between each symptom and the disease. Later, they compute a pairwise cosine-similarity matrix between these vectors and only the most significant similarities are considered to construct a network of weighted links between diseases (Zhou et al., 2014). Analysis of the resulting network shows that disease pairs with high symptom similarity are more likely to share associated genes and PIs. This symptom-based disease network is also organized in highly interconnected communities of similar diseases, which shows that similar symptoms imply similar disorders.
The recent work of Menche and colleagues is quite relevant, as it shows that, despite its incompleteness and biases, the current human PIN can be mined and integrated with disease data to uncover pathobiological relationships between disorders and better understand their etiology. After compiling a network of roughly 140k interactions between more than 13k human proteins, nodes are labeled with their associated diseases with the help of OMIM and a set of 299 disorders defined by MeSH. Although the disease module hypothesis predicts that proteins associated with the same trait should be highly interconnected (Barabási et al.,
5. Conclusion
Viewing the relationships between cell compartments and their constituting molecules as a complex circuitry of tightly interconnected components is widespread in systems biology (Vidal et al., 2011). This has led to breakthroughs that the study of the individual system components would not have made possible (Takahashi and Yamanaka, 2006; Levine and Oren,
Since the identification of biological features to distinguish between interacting and non-interacting proteins is very difficult, mining the topological characteristics of PINs is useful in the reliability assessment and prediction of PIs (Cannistraci et al.,
As the quantity and quality of molecular datasets increase, network science offers a new means to analysing interacting gene products at a systems level (Loscalzo and Barabasi,
Conflict of interest statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Statements
Acknowledgments
The author would like to thank Antonia Forshaw for proofreading the article, as well as Prof. Timothy Ravasi and Dr. Carlo Cannistraci for encouraging the preparation of this review.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AdamicL. A.AdarE. (2003). Friends and neighbors on the Web. Soc. Netw.25, 211–230. 10.1016/S0378-8733(03)00009-1
2
Alanis-LobatoG.CannistraciC. V.ErikssonA.ManicaA.RavasiT. (2015). Highlighting nonlinear patterns in population genetics datasets. Sci. Rep.5:8140. 10.1038/srep08140
3
Alanis-LobatoG.CannistraciC. V.RavasiT. (2013). Exploitation of genetic interaction network topology for the prediction of epistatic behavior.Genomics102, 202–208. 10.1016/j.ygeno.2013.07.010
4
Alanis-LobatoG.CannistraciC. V.RavasiT. (2014). Exploring the genetics underlying autoimmune diseases with network analysis and link prediction, in Middle East Conference on Biomedical Engineering (MECBME) (Doha: IEEE), 167–170.
5
AlbertR.BarabásiA.-L. (2002). Statistical mechanics of complex networks. Rev. Mod. Phys.74, 47–97. 10.1103/RevModPhys.74.47
6
BarabásiA.-L.GulbahceN.LoscalzoJ. (2011). Network medicine: a network-based approach to human disease.Nat. Rev. Genet.12, 56–68. 10.1038/nrg2918
7
BoguñáM.KrioukovD.ClaffyK. C. (2009). Navigability of complex networks. Nat. Phys.5, 74–80. 10.1038/nphys1130
8
CannistraciC. V.Alanis-LobatoG.RavasiT. (2013a). From link-prediction in brain connectomes and protein interactomes to the local-community-paradigm in complex networks. Sci. Rep.3, 1–13. 10.1038/srep01613
9
CannistraciC. V.Alanis-LobatoG.RavasiT. (2013b). Minimum curvilinearity to enhance topological prediction of protein interactions by network embedding.Bioinformatics29, i199–i209. 10.1093/bioinformatics/btt208
10
ChenJ.ChuaH. N.HsuW.LeeM.-L.NgS.-K.SaitoR.et al. (2006a). Increasing confidence of protein-protein interactomes.Genome Informat.17, 284–297.
11
ChenJ.HsuW.LeeM. L.NgS.-K. (2006b). Increasing confidence of protein interactomes using network topological metrics.Bioinformatics22, 1998–2004. 10.1093/bioinformatics/btl335
12
ChuaH. N.SungW.-K.WongL. (2006). Exploiting indirect neighbours and topological weight to predict protein function from protein-protein interactions.Bioinformatics22, 1623–1630. 10.1093/bioinformatics/btl145
13
ClausetA.MooreC.NewmanM. E. (2008). Hierarchical structure and the prediction of missing links in networks.Nature453, 98–101. 10.1038/nature06830
14
DiceL. R. (1945). Measures of the amount of ecologic association between species. Ecology26, 297–302.
15
GohK.-I.CusickM. E.ValleD.ChildsB.VidalM.BarabásiA.-L. (2007). The human disease network. Proc. Natl. Acad. Sci. U.S.A.104, 8685–8690. 10.1073/pnas.0701361104
16
GuimeràR.Sales-PardoM. (2009). Missing and spurious interactions and the reconstruction of complex networks. Proc. Natl. Acad. Sci. U.S.A.106, 1–6. 10.1073/pnas.0709640104
17
HighamD. J.RasajskiM.PrzuljN. (2008). Fitting a geometric graph to a protein-protein interaction network.Bioinformatics24, 1093–1099. 10.1093/bioinformatics/btn079
18
JaccardP. (1912). The distribution of flora in the alpine zone. New Phytol.11, 37–50.
19
KatzL. (1953). A new status index derived from sociometric analysis. Psychometrika18, 39–43.
20
KrioukovD.KitsakM.SinkovitsR. S.RideoutD.MeyerD.BoguñáM. (2012). Network cosmology.Sci. Rep.2, 1–6. 10.1038/srep00793
21
KrioukovD.PapadopoulosF.KitsakM.VahdatA.BoguñáM. (2010). Hyperbolic geometry of complex networks. Phys. Rev. E82:036106. 10.1103/physreve.82.036106
22
KuchaievO.RasajskiM.HighamD. J.PrzuljN. (2009). Geometric de-noising of protein-protein interaction networks.PLoS Comput. Biol.5:e1000454. 10.1371/journal.pcbi.1000454
23
LageK.KarlbergE. O.StørlingZ. M.ÓlasonP. I.PedersenA. G.RiginaO.et al. (2007). A human phenome-interactome network of protein complexes implicated in genetic disorders. Nat. Biotechnol.25, 309–316. 10.1038/nbt1295
24
LeskA. M. (2007). Systems biology, in Introduction to Genomics, Chapter 7, 1st Edn. (New York, NY: Oxford University Press), 359–405.
25
LevineA. J.OrenM. (2009). The first 30 years of p53: growing ever more complex. Nat. Rev. Cancer9, 749–758. 10.1038/nrc2723
26
LiuG.WongL.ChuaH. N. (2009). Complex discovery from weighted PPI networks.Bioinformatics25, 1891–1897. 10.1093/bioinformatics/btp311
27
LiuY.-Y.SlotineJ.-J.BarabásiA.-L. (2011). Controllability of complex networks.Nature473, 167–173. 10.1038/nature10011
28
LoscalzoJ.BarabasiA.-L. (2011). Systems biology and the future of medicine. Wiley Interdiscip. Rev. Syst. Biol. Med.3, 619–627. 10.1002/wsbm.144
29
LüL.JinC.-H.ZhouT. (2009). Similarity index based on local paths for link prediction of complex networks. Phys. Rev. E80:046122. 10.1103/PhysRevE.80.046122
30
LüL.ZhouT. (2011). Link prediction in complex networks: a survey. Phys. A390, 1150–1170. 10.1016/j.physa.2010.11.027
31
MencheJ.SharmaA.KitsakM.GhiassianS. D.VidalM.LoscalzoJ.et al. (2015). Uncovering disease-disease relationships through the incomplete interactome. Science347, 1257601–1257601. 10.1126/science.1257601
32
MeyerM. J.DasJ.WangX.YuH. (2013). INstruct: a database of high-quality 3D structurally resolved protein interactome networks. Bioinformatics29, 1577–1579. 10.1093/bioinformatics/btt181
33
NewmanM. (2001). Clustering and preferential attachment in growing networks. Phys. Rev. E64, 1–4. 10.1103/PhysRevE.64.025102
34
OliverS. (2000). Guilt-by-association goes global.Nature403, 601–603. 10.1038/35001165
35
PapadopoulosF.KitsakM.SerranoM. A.BoguñáM.KrioukovD. (2012). Popularity versus similarity in growing networks.Nature489, 537–540. 10.1038/nature11459
36
PrzuljN.CorneilD. G.JurisicaI. (2004). Modeling interactome: scale-free or geometric?Bioinformatics20, 3508–3515. 10.1093/bioinformatics/bth436
37
RavasiT.SuzukiH.CannistraciC. V.KatayamaS.BajicV. B.TanK.et al. (2010). An atlas of combinatorial transcriptional regulation in mouse and man. Cell140, 744–752. 10.1016/j.cell.2010.01.044
38
ResnikP. (1999). Semantic similarity in a taxonomy: an information-based measure and its application to problems of ambiguity in natural language. J. Artif. Intell. Res.11, 95–130.
39
RollandT.TasM.SahniN.YiS.LemmensI.FontanilloC.et al. (2014). A proteome-scale map of the human interactome network. Cell159, 1212–1226. 10.1016/j.cell.2014.10.050
40
SaitoR.SuzukiH.HayashizakiY. (2002). Interaction generality, a measurement to assess the reliability of a protein-protein interaction.Nucleic Acids Res.30, 1163–1168. 10.1093/nar/30.5.1163
41
SaitoR.SuzukiH.HayashizakiY. (2003). Construction of reliable protein-protein interaction networks with a new interaction generality measure. Bioinformatics19, 756–763. 10.1093/bioinformatics/btg070
42
SchaeferM. H.FontaineJ. F.VinayagamA.PorrasP.WankerE. E.Andrade-NavarroM. A. (2012). HIPPIE: integrating protein interaction networks with experiment based quality scores. PLoS ONE7:e31826. 10.1371/journal.pone.0031826
43
SchaeferM. H.LopesT. J. S.MahN.ShoemakerJ. E.MatsuokaY.FontaineJ.-F.et al. (2013). Adding protein context to the human protein-protein interaction network to reveal meaningful interactions. PLoS Comput. Biol.9:e1002860. 10.1371/journal.pcbi.1002860
44
SzklarczykD.FranceschiniA.KuhnM.SimonovicM.RothA.MinguezP.et al. (2011). The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored.Nucleic Acids Res.39, D561–D568. 10.1093/nar/gkq973
45
SzklarczykD.FranceschiniA.WyderS.ForslundK.HellerD.Huerta-CepasJ.et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res.43, D447–D452. 10.1093/nar/gku1003
46
TakahashiK.YamanakaS. (2006). Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors. Cell126, 663–676. 10.1016/j.cell.2006.07.024
47
TenenbaumJ. B.de SilvaV.LangfordJ. C. (2000). A global geometric framework for nonlinear dimensionality reduction.Science290, 2319–2323. 10.1126/science.290.5500.2319
48
VidalM.CusickM. E.BarabásiA.-L. (2011). Interactome networks and human disease.Cell144, 986–998. 10.1016/j.cell.2011.02.016
49
WangJ.DuZ.PayattakoolR.YuP.ChenC.-F. (2007). A new method to measure the semantic similarity of GO terms.Bioinformatics23, 1274–1281. 10.1093/bioinformatics/btm087
50
WangX.WeiX.ThijssenB.DasJ.LipkinS. M.YuH. (2012). Three-dimensional reconstruction of protein networks provides insight into human genetic disease.Nat. Biotechnol.30, 159–164. 10.1038/nbt.2106
51
WuX.JiangR.ZhangM. Q.LiS. (2008). Network-based global inference of human disease genes. Mol. Syst. Biol.4:189. 10.1038/msb.2008.27
52
YangH.NepuszT.PaccanaroA. (2012). Improving GO semantic similarity measures by exploring the ontology beneath the terms and modelling uncertainty.Bioinformatics28, 1383–1389. 10.1093/bioinformatics/bts129
53
YouZ.-H.LeiY.-K.GuiJ.HuangD.-S.ZhouX. (2010). Using manifold embedding for assessing and predicting protein interactions from high-throughput experimental data.Bioinformatics26, 2744–2751. 10.1093/bioinformatics/btq510
54
YuG.LiF.QinY.BoX.WuY.WangS. (2010). GOSemSim: an R package for measuring semantic similarity among GO terms and gene products.Bioinformatics26, 976–978. 10.1093/bioinformatics/btq064
55
ZhouT.LüL.ZhangY.-C. (2009). Predicting missing links via local information. Eur. Phys. J. B71, 623–630. 10.1140/epjb/e2009-00335-8
56
ZhouX.MencheJ.BarabásiA.-L.SharmaA. (2014). Human symptoms-disease network. Nat. Commun.5:4212. 10.1038/ncomms5212. Available online at: http://www.nature.com/ncomms/2014/140626/ncomms5212/full/ncomms5212.html
57
ZhuY.ZhangX.-F.DaiD.-Q.WuM.-Y. (2013). Identifying spurious interactions and predicting missing interactions in the protein-protein interaction networks via a generative network model. IEEE/ACM Trans. Comput. Biol. Bioinformatics10, 219–225. 10.1109/TCBB.2012.164
Summary
Keywords
interactome, proteome, network, reliability, prediction, medicine, disease, pathogenesis
Citation
Alanis-Lobato G (2015) Mining protein interactomes to improve their reliability and support the advancement of network medicine. Front. Genet. 6:296. doi: 10.3389/fgene.2015.00296
Received
29 May 2015
Accepted
07 September 2015
Published
23 September 2015
Volume
6 - 2015
Edited by
Spyros Petrakis, Aristotle University of Thessaloniki, Greece
Reviewed by
Anshu Bhardwaj, Council of Scientific and Industrial Research, India; Maria I. Klapa, Foundation for Research & Technology Hellas, Greece
Updates

Check for updates
Copyright
© 2015 Alanis-Lobato.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gregorio Alanis-Lobato, Institute of Molecular Biology, Ackermannweg 4, 55128 Mainz, Germany g.alanis-lobato@imb-mainz.de
This article was submitted to Systems Biology, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.