{"paper_id":"e593cb62-126b-44b8-997e-1b2bc83699bc","body_text":"During the past decades, significant progress has been made in our understanding of human diseases [ 1 ]. However, the genetic architectures of complex diseases are still largely unclear. Many common diseases tend to be related to each other, and it is speculated that they may share common genetic origin. Thus, studying the correlations of human diseases has the potentials of better understanding the genotype to phenotype mapping [ 2 ,  3 ] and better predicting disease association genes [ 4 ,  5 ,  6 ,  7 ,  8 ]. Moreover, learning which diseases are correlated can help use existing drugs to treat multiple similar diseases [ 9 ,  10 ,  11 ,  12 ,  13 ].\nMeanwhile, network science is a rising field where entities and their complex relationships are studied on a global scale [ 14 ,  15 ,  16 ], and has seen increasing applications to perform advanced analysis on biomedical data [ 17 ,  18 ,  19 ,  20 ,  21 ,  22 ,  23 ,  24 ]. There are various cellular components in the human body that interact with each other within the same cell or across different cells [ 15 ]. A network called the  human interactome  can be constructed according to the interactions of those different cellular components. Each component can be represented as a vertex in the network and interactions among them can be captured as links (or edges) connecting pairs of the cellular components. Those cellular components can be proteins or metabolites, and the network refers to protein-protein interaction (PPI) network [ 25 ,  26 ,  27 ] or metabolic network [ 28 ,  29 ,  30 ].\nSome studies aimed at identifying the correlations among diseases through network analysis [ 15 ,  31 ,  32 ]. Goh  et al . [ 33 ] constructed a human disease network (HDN) by connecting pairs of diseases when they share common association genes. Of 1,284 diseases in the HDN, 867 have at least one link to other diseases, and 516 form a giant component, suggesting that the genetic origins of most diseases, to some extent, are shared with other diseases. Moreover, the HDN naturally and visibly clustered according to major disease classes such as cancer cluster and neurological disease cluster. Zhou  et al . [ 34 ] extracted over twenty million bibliographic records from PubMed [ 35 ] in order to obtain 147,978 connections between 322 symptoms and 4,219 diseases. A human symptoms-disease network (HSDN) was then constructed and was able to show the symptom similarity between all pairs of diseases (7,488,851 links) in the network. The weight of links represented the similarity of symptoms between two diseases. They showed that the correlations among diseases were significantly related to the genetic associations that each pair of diseases had in common as well as the interactions between their related proteins. Lee  et al . [ 36 ] built a disease metabolism network in order to study disease comorbidity for better disease prediction and prevention. Two diseases are connected if enzymes associated with them catalyze adjacent metabolic reactions. Their results show that diseases with higher degrees, i.e., connecting with many other diseases, have a higher rate of prevalence and mortality.\nMeasuring the centrality of vertices helps identify important vertices in the network in terms of connecting to all other vertices. Centrality measures have been used frequently to analyze biological networks over the past decades [ 37 ,  38 ,  39 ]. The most common centrality measures include degree (the total number of neighbors), closeness (the total distance to all other vertices), and betweenness (the fraction of locating on the shortest paths of all pairs of vertices) [ 40 ]. Despite wide applications in biological networks, these centrality measures are rather general and may not be able to capture all the properties of vertices in the context of biological networks. Furthermore, closeness and betweenness have high computational complexity due to the fact that pair-wise shortest paths in a network need to be enumerated in order to compute the centralities. Therefore, carefully tailored and more efficient centrality measures are needed for specific network of interest, in this study, the human disease network.\nKöhler  et al . [ 41 ] proposed a vertex importance measure for disease genes in the context of PPI networks. They used a random walk strategy to assess the distance between vertices in the network, and reported improved performance comparing with conventional distance-based centrality measures. Wu  et al . [ 42 ] integrated PPI networks with gene expression data in order to rank disease genes associated with various cancers. They showed that their method was able to find replicable high-rank genes using different datasets. Martinez  et al . [ 43 ] proposed a generic vertex prioritization method using the idea of propagating information across data networks and measuring the correlation between the propagated values for a query and a target set of entities. The authors tested their method by ranking disease genes associated with Alzheimer’s disease, diabetes mellitus type 2 and breast cancer. They reported some new high-rank association genes that could bring new insights into the diseases.\nIn this article, we propose a new method for the construction of a  weighted human disease network (WHDN) and a new centrality measure to identify the most important diseases. First we use a large database of disease-gene associations to build a weighted bipartite disease-gene network, and then construct a weighted disease network where link weights capture the strengths of the pairwise disease correlations. After the backbone extraction of the WHDN, we design a centrality measure for the context of the WHDN that considers not only the degree of a vertex but also the importance of its incident edges. Then we compare our new centrality measure with degree, closeness and betweenness by evaluating the network efficiency decline rate with the removal of top-ranked vertices by each centrality measurement. Finally, we present the top 30 diseases ranked by our centrality measure in our WHDN and discuss their biological implications.\n\nGiven the multiple-step pipeline structure of this study, we show the result of each step after the description of the corresponding method. The source code of our analysis and network files are accessible through the Github link:  https://github.com/MIBlab-MUN/vertex-centrality-DILW .\nThe data used in this project describe disease-gene associations (DGAs) from multiple curated databases including UNIPROT [ 44 ], CTD (human subset) [ 45 ], PsyGeNET [ 46 ], Orphanet [ 47 ], and HPO [ 48 ]. The disease-gene association data are collected by DisGeNet group, available on DisGeNET v4.0 [ 49 ]. The current version of the data set contains 130,821 DGAs, between 13,075 diseases and 8,949 genes. Each DGA is assigned with a score  a i j , for disease  i  and gene  j , within the range of [0, 1] based on its level of evidence, the number and the type of database sources supporting the DGA, and the number of publications verifying the association between the gene and the disease [ 49 ]. We first clean up the data in order to ensure that all diseases and genes in the dataset are unique and that there is no replication of disease-gene associations. Next, since we would like to consider the correlation among all diseases, we keep diseases and syndromes in the dataset for our analysis and remove injuries or poisonings, anatomical abnormalities, acquired abnormalities, mental or behavioral dysfunctions, signs or symptoms, findings, congenital abnormalities, neoplastic processes, and pathologic functions. We use DisGeNet web-based application [ 49 ] for this filtering.\nThe best representation for depicting the associations among genes and diseases is a bipartite graph, which is called the  disease-gene association network  in this research. The bipartite graph contains two different sets of vertices. One set includes diseases and the other one contains genes. By definition, no edge is allowed to connect a pair of vertices in the same set of vertices in a bipartite graph. That is, there can be no link either between a pair of diseases or a pair of genes. There is an edge between a gene and a disease if there is an association between them. Their link weight is assigned as the score  a i j , for disease  i  and gene  j , computed in the DGA database described in the previous section. A sample subgraph of the bipartite network is shown in  Fig 1 .\nThe bipartite network has two sets of vertices, i.e., genes and diseases, represented by rectangle and gray ellipses respectively. An edge connects a disease and a gene if there is a known association between them. The weight of an edge reflects the strength of the DGA  a i j  between disease  i  and gene  j .\nFig 2  depicts the degree distributions of diseases and genes in the bipartite disease-gene association network. For the set of diseases, the maximum degree is 564, of the disease  epilepsy , and the average degree is 5.43. In  Fig 2a ), the degree distribution of the diseases is right-skewed and heavy-tailed, indicated by the straight linear fit on a log-log scale. For the set of genes, the maximum degree is 111, of the gene  LMNA , and the average degree is 5.81.\nDegree distribution of a) diseases and b) genes in the bipartite disease-gene association network. The distributions are shown on a log-log scale.\nThe bipartite network is comprised of multiple connected components with a single giant component.  Fig 3  shows its distribution of the size of connected components. The giant component has 10,212 vertices consisting of 5,278 diseases and 4,934 genes. Apart from the giant component, all other connected components are small with a size varying from two to nine, and most of them are only single pairs of one disease and one gene.  Fig 3  shows that there is a considerable number of components with two vertices, i.e., 844 isolated disease-gene pairs. Since we are interested in investigating the large-scale genetic correlations of human diseases, we focus on the giant component of the disease-gene bipartite network in the downstream analyses.\nThe network has a single giant component with 10,212 vertices, and the majority of other connected components are of size two, i.e., consisting of only one disease and one gene.\nWe construct the WHDN using the giant connected component of the bipartite disease-gene network. We use  D  and  G  to denote sets of 5,278 diseases and 4,934 genes respectively in the giant connected component. In the WHDN, an edge links two diseases  i  and  j  if they have at least one association gene in common, and the weight of the edge,  w ij , is computed based on the number of shared association genes, as well as the strengths of those associations.\nSuch a weight definition is inspired by Newman’s study on scientific collaboration networks [ 14 ], where vertices are scientists and two scientists are connected by an unweighted edge if they have coauthored one or more scientific papers together. To define the strength of the tie between two connected scientists, two factors are considered. First, two scientists whose names appear on a paper together with many other coauthors know one another less well on average than two who are the sole authors of a paper. Thus, the collaborative ties are weighted inversely according to the number of coauthors of a paper. Second, authors who have written many papers together will know one another better on average than those who have written few papers together. Thus, all coauthored papers are added up to account for the tie strength of two scientists.\nHere, similarly, first we consider that the correlation of two diseases through a gene is stronger when they are the sole associated diseases with this gene than when there are many other diseases associated with the same gene. Second, the correlation of two diseases is considered stronger when they share more genes through stronger associations than less genes or weaker associations. Thus, we extend Newman’s method to weighted graph and define the weight of edge  w ij  between two diseases  i  and  j  as\n w i j = ∑ g ∈ G δ i g δ j g ( a i g + a j g ) s g , (1) \nwhere  δ i g  is one if disease  i  and gene  g  have a DGA, and zero otherwise.  a i g  is the score of their DGA assessed by DisGeNET as discussed in the previous section, and  s g  is the strength of gene  g  as a vertex in the bipartite disease-gene network, defined as the sum of the scores of the DGAs between gene  g  and its directly linked diseases,\n s g = ∑ i ∈ D a i g . (2) \nSuch a weight definition indicates that the correlation strength of two diseases is weighted inversely according to the strengths of the genes they share, and is proportional to the total number of genes they share and the strengths of their DGAs.\nFor example, in  Fig 1 , the weight between diseases  contact dermatitis  (CD) and  white sponge nevus 1  (WSN1) is calculated as follows,\n w C D , W S N 1 = ∑ g ∈ G δ C D g δ W S N 1 g ( a C D g + a W S N 1 g ) s g = a C D K R T 4 + a W S N 1 K R T 4 s K R T 4 = 0 . 2 + 0 . 48 0 . 881 = 0 . 7718 .\nNote that the weight of two diseases can be greater than one when they share multiple genes. For example the weight between diseases WSN1 and  hereditary mucosal Leukokeratosis  (HML) is calculated as follows,\n w W S N 1 , H M L = ∑ g ∈ G δ W S N 1 g δ H M L g ( a W S N 1 g + a H M L g ) s g = a W S N 1 K R T 4 + a H M L K R T 4 s K R T 4 + a W S N 1 K R T 13 + a H M L K R T 13 s K R T 13 = 0 . 48 + 0 . 201 0 . 881 + 0 . 2 + 0 . 2008 0 . 6008 = 0 . 7729 + 0 . 6671 = 1 . 44 .\nSince the WHDN is constructed using vertices from the giant component of the bipartite disease-gene association network, it only has a single connected component with all 5,278 vertices in the disease set  D . Two vertices have an edge connecting them if the represented two diseases have at least one shared gene, and the edge weight is assessed as described above. The WHDN has 11,2324 edges and an average vertex degree of 42.56. That is, a disease correlates with on average 42.56 other diseases with varying strengths.  Fig 4  depicts the distribution of all the edge weights in the WHDN. As we can see that a large number of edge weights are of small values and may not be particularly interesting for the subsequent analysis. Those weak edges not only add computational overhead to the network analysis, but also render the network difficult to interpret. Therefore, next we perform an edge reduction and only extract the most meaningful structure of the network.\nThe weight of an edge quantifies the shared genetic background of two connected diseases. There are 112,324 edges in the graph with weights ranging from 0.0152 to 22.4506.\nThe most straightforward strategy for network reduction may be to use a global weight threshold and remove all links that have weights lower than the threshold. However, such a global thresholding strategy is somewhat arbitrary and may overlook the network information present below the cutoff scale. Here, to preserve the multi-scale backbone of the WHDN while removing less relevant and meaningful edges we use a multi-scale filtering method proposed by Serrano  et al . [ 50 ]. Such a multi-scale backbone exaction algorithm has been used to reduce the network size while preserving the meaningful structure of biological networks in multiple studies [ 34 ,  51 ,  52 ,  53 ].\nFirst, the weight of edge linking vertex  i  with its neighbor  j  can be normalized as\n N i j = w i j s i , (3) \nwhere  s i  is the vertex strength, i.e., the sum of weights incident to vertex  i , similar to  Eq (2)  and defined as\n s i = ∑ j ∈ Γ i w i j , (4) \nwhere Γ i  is the set of vertex  i ’s neighbors. Therefore, there are two different normalized values for a link  e ij  using the strengths of its two end vertices  s i  and  s j  as the denominator.\nSecond, a null model is used to assess the expectation if the weights of links connecting to a particular vertex were distributed randomly. That is, the normalized weight  N ij  that corresponds to the link connecting to a certain vertex of degree  k  is produced by a random assignment from an uniform distribution. Thus the probability density function for the variable taking a particular value  x  is\n p ( x ) d x = ( k − 1 ) ( 1 − x ) k − 2 d x . (5)\nThen, to identify whether the probability,  β ij , of link weight  N ij  is compatible with the null model with a threshold  β  is given as\n β i j = 1 − ( k − 1 ) ∫ 0 N i j ( 1 − x ) k − 2 d x < β . (6)\nAll links with computed  β ij  lower than a given threshold  β  are preserved in the network. Note that each edge has two different values  β ij  and  β ji . For solving this problem, OR and AND rules can be used. Under the first rule, if either  β ij  and  β ji  is lower than  β , the link will be preserved. In the second case, an edge is preserved if both  β ij  and  β ji  are lower than  β . Darabos  et al . [ 51 ] empirically found that the AND rule preserve the network features better than using the OR rule in the context of human phenotype networks. In this article, the AND rule is adopted to reduce the size of the network by removing the links which are less relevant.\nTo find the best cutoff for  β , we calculate clustering coefficient, percentage of remaining vertices and links, and total weight of the networks as a function of  β  in the range [0, 1].  Fig 5  shows the results of network metrics as a function of  β  cutoffs. We choose a  β  cutoff when the clustering coefficient and the remaining vertices and weights are maximally preserved while as many links are removed as possible. Accordingly, the cutoff  β  = 0.501 can be determined, shown as the vertical dashed line in the figure.\nCC represents clustering coefficient, %Vertices is the percentage of remaining vertices, %Weights is the percentage of weights left after removing links, and %Links is the percentage of remaining links.\nAfter the backbone extraction, the WHDN has 4,898 vertices and 38,275 edges. Those vertices are no longer connected in a single component.  Fig 6  shows the size distribution of its connected components. There is a giant component with 4,810 vertices and its degree distribution is shown in  Fig 7 . Again the degree distribution is heavy tailed and resembles a power-law relationship. The vertex  epilepsy  has the highest degree of 576. This giant component will be the focus for our next step analysis, i.e., measuring vertex importance in order to find the most central diseases in terms of correlating with other diseases.\nThe network has a single giant component with 4,810 vertices.\nThe distribution is shown on a log-log scale.\nAlthough various vertex centrality measures have been proposed in the literature [ 37 ,  38 ,  40 ,  41 ,  54 ], the quantification of the importance of a vertex in a network is often context-specific. For some networks, measuring degree may suffice since a vertex can be considered important when its number of neighbors is the sole criterion. For some networks, e.g., information communication networks, a vertex may be considered more important if its distances to all other vertices are short, then closeness centrality serves this purpose well. For our WHDN, a disease is considered important if it correlates with many other diseases (degree) as well as if the correlations are themselves very important (edge importance).\nWe propose a vertex importance measure for WHDN by extending a centrality measure for unweighted networks proposed by Liu  et al . [ 54 ]. This measure assesses the centrality of a vertex based on both its degree and the importance of its incident links (DIL centrality). For its extension on weighted graphs, we name it the DIL-W centrality.\nFirst, in the context of unweighted graph, the importance of a link  e ij  that connects vertex  v i  and  v j  can be calculated as follows:\n I e i j = U e i j λ e i j , (7) \nwhere  U e i j = ( k i − t − 1 ) ( k j − t − 1 )  and  λ e i j = t 2 + 1 . Following the convention,  k i  and  k j  are the degrees of vertex  v i  and  v j , respectively, and  t  is the number of triangles with one edge being  e ij .\nSubsequently, the contribution that vertex  v i  makes to the importance of  e ij  is computed as\n C v i v j = I e i j · k i − 1 k i + k j − 2 , (8) \nwhere  j  ∈ Γ i , and Γ i  is the neighborhood of vertex  i .\nThen, the DIL centrality of vertex  v i  is calculated by combining both its degree and the importance of its incident links,\n DIL v i = k i + ∑ v j ∈ Γ i C v i v j . (9)\nFor weighted networks, we modify the computation of  U  in  Eq (7)  as\n U e i j = ( s i − t i ) · ( s j − t j ) , (10) \nwhere  s i  is the strength of vertex  v i , calculated as in  Eq (4) , and  t i  is the weight sum of links incident to vertex  v i  that form triangles with  e ij . This follows the intuition that first an edge is considered more important when its two end vertices have higher strengths. Second, the importance of an edge is reduced when it has alternative two-hop paths connecting the same set of end vertices. Therefore, we subtract  t i  from  s i  in  Eq (10) .\nWe define λ for weighted graphs as\n λ e i j = t i + t j 2 + 1 . (11)\nFinally, the importance of a vertex can be measured by\n DIL-W v i = s i + ∑ v j ∈ Γ i C v i v j , (12) \nwhere  C v i v j  is defined as\n C v i v j = I e i j · s i s i + s j . (13) \nNote that, if we remove the second component in the definition of DIL-W, the centrality measure simply becomes vertex strength, i.e., weighted degree.\nIn the weighted graph given in  Fig 8 , vertex  a  has a higher strength but a lower degree than vertex  b . We compute their DIL-W centralities and investigate which one is more central when both factors are considered.\nFirst we have their strength values  s a  = 0.9 + 0.3 + 0.5 + 0.6 = 2.3, and  s b  = 0.2 + 0.11 + 0.2 + 0.7 + 0.5 = 1.71. Their neighborhoods are Γ a  = { b ,  c ,  d ,  g } and Γ b  = { a ,  c ,  e ,  f ,  g }. For vertex  a ,\n ∑ v j ∈ Γ a C a v j = C a b + C a c + C a d + C a g , \nwhere\n C a b = I e a b · s a s a + s b , \nand\n I e a b = U e a b λ e a b = ( s a − t a ) · ( s b − t b ) t a + t b 2 + 1 . \nWe have\n t a = w a c + w a g = 0 . 3 + 0 . 6 = 0 . 9 , \nand\n t b = w b c + w b g = 0 . 2 + 0 . 7 = 0 . 9 . \nSo\n C a b = ( s a − t a ) · ( s b − t b ) t a + t b 2 + 1 · s a s a + s b = ( 2 . 3 − 0 . 9 ) · ( 1 . 71 − 0 . 9 ) 0 . 9 + 0 . 9 2 + 1 · 2 . 3 2 . 3 + 1 . 71 = 0 . 3423 \nWe can also have\n C a c = 0 . 3285 , C a d = 1 . 4878 , and C a g = 0 . 4312 . \nThen\n DIL-W a = s a + ∑ v j ∈ Γ a C a v j = 2 . 3 + ( 0 . 3423 + 0 . 3285 + 1 . 4878 + 0 . 4312 ) = 4 . 8898 . \nSimilarly, we can compute the DIL-W centrality of vertex  b  DIL-W b  = 2.8916. Therefore, based on both the degree and importance of incident edges, vertex  a  is considered more important than vertex  b .\nWe apply the DIL-W centrality measurement to the giant component of the backbone of WHDN, the distribution is shown in  Fig 9 . The DIL-W scores have a high dynamic range, from 0.0610 to 80688.1129. The majority of the vertices have low scores and a few number of vertices can have scores that are greater by orders of magnitude.\nWe compare our DIL-W measurement with three most commonly used centralities, i.e., degree, closeness, and betweenness, when applied to the giant component of the backbone of WHDN. For weighted graphs, degree centrality is calculated as vertex strength given by  Eq (4) . Closeness and betweenness are shortest-path-based centralities. Shortest path computation can be extended for weighted graph as follows,\n d i j w = m i n ( 1 w i h + . . . + 1 w h j ) . (14) \nHere  d i j w  denotes the weighted distance between vertex  i  and  j , and  w ih  is the weight of the edge linking vertex  i  and  h , where  h  is the intermediate vertex between vertices  i  and  j . Since in our WHDN edge weight suggests strength, the distance between two vertices is the minimum sum of the inverse of edge weight along the path connecting them. Once the weighted distance is defined, closeness and betweenness can be calculated by their original definitions.\nFig 10  shows the correlation of DIL-W scores with a) degree, b) closeness, and c) betweenness centralities. As we can see, there is a positive correlation between DIL-W measure and all other three vertex centrality measures. The Spearman’s rank correlation coefficient is 0.672 comparing DIL-W with closeness, is 0.71 comparing DIL-W with betweenness, and is 0.947 comparing DIL-W with degree.\nCorrelation of DIL-W scores with a) degree centrality, b) closeness centrality, and c) betweenness centrality in the WHDN.\nTo evaluate our new vertex importance quantification method, DIL-W, we measure the network efficiency before and after we remove the most important vertices in the WHDN. In the context of the WHDN, the network efficiency indicates the extend to which the original connectivity of the network is maintained. We calculate the decline rate of network efficiency after removing  m  top-rank vertices. The network efficiency [ 55 ] is computed based on the connectivity of a network. A higher connectivity suggests a higher network efficiency. The network efficiency is defined by\n η = 1 n ( n − 1 ) ∑ v i ≠ v j ∈ V 1 d i j , (15) \nwhere  n  is the total number of vertices in the network,  V  is the vertex set, and  d ij  is the weighted distance between vertex  v i  and  v j . Thus, the decline rate of the network efficiency is calculated as\n μ = 1 − η η 0 , (16) \nwhere  η 0  is the efficiency of the original network, and  η  is the network efficiency after some vertices are removed.\nWhen a more importance vertex is removed, we expect to see a greater decline rate of the network efficiency. Thus we can use  μ  as an indicator for the actual impact of removing a vertex in the network.  Fig 11  shows the decline rate of the network efficiency when we remove each of the top 40 vertices ranked by a) degree (DC), b) closeness (CC), c) betweenness (BC), and d) DIL-W. Further removal of top ranked vertices could be investigated but was not included in the current study given the high computational demand. As shown in the figure, we do not observe a monotonic relationship across all four centrality methods. However, the correlation analysis shows that our method, DIL-W, has a slighter stronger negative correlation between the decline rate and the rank of the removed vertex than the other three. The Spearman’s rank correlation coefficient,  ρ , for degree, closeness, and betweenness is −0.1801, −0.0017, and −0.0679, respectively. In comparison, DIL-W has a negative correlation coefficient −0.2698.\nDecline rate of network efficiency after removing a single vertex ranked by a) degree centrality (DC), b) closeness centrality (CC), c) betweenness centrality (BC), and d) DIL-W.\nWe also consider removing all  m  top-rank vertices at once and see how this accumulative removal affects the efficiency of the network.  Fig 12  shows the decline rate of the network efficiency after removing all top  m  vertices ranked by different centrality measures. The graph shows that the proposed method, DIL-W, has the highest decline rate of network efficiency for 57.5% of the data points, while betweenness, closeness, and degree have 27.5%, 10%, and 5%, respectively. This suggests that DIL-W is able to select a set of more important vertices comparing with the other three centrality measures. As seen in  Fig 12 , the four methods are very comparable until the top 11 diseases are removed from the network. Then DIL-W has a significant higher network efficiency decline rate than the rest. Betweenness centrality catches up around point 30 and becomes very comparable afterwards.\nSince one main contribution of our study is to add edge weights to the HDN, we collect another set of results by computing vertex centralities without the consideration of edge weights. That is, the network structure remains the same but edges now do not carry weights, then the weighted DC, CC, and BC simply become their original definitions for un-weighted graphs, and DIL-W is replaced by the original DIL. The comparison is depicted in  S1 Fig , which shows that excluding edge weights results in very similar vertex rankings by various centrality measures and essentially no significant difference in evaluation.\nTable 1  shows the top 30 diseases ranked by our DIL-W method, their degrees, and their neighbors that have the strongest correlations (i.e., edge weights). References that support the known comorbidity of the disease pairs are also given.\nThe table shows the diseases, along with their rankings, their degrees in the WHDN, their direct neighbors with the strongest edge weight, and literature references that have discussed the correlations of the disease pairs.\nIn addition, we compare the top 30 diseases ranked by different centrality measures (see  Fig 13 ). The figure shows the top 30 diseases ranked by our proposed DIL-W (x-axis), as well as their rankings by other three centrality measures. If a disease is not among the top 30 ranks by a centrality measure, the data point will be shown as a zero on the x-axis. We see that 18 out of 30 top diseases ranked by DIL-W are not picked up by at least one other centrality measures. This comparison result further justifies the utility of our proposed centrality measure on finding central vertices (diseases) undetected by other conventional centrality measures.\nThe diseases on the x-axis are ordered based on their ranks by DIL-W. A data point landing on the x-axis indicates that the corresponding disease was not among the top 30 ranked by a centrality measure.\n\nIn this article, we use a network-based analysis to identify important human diseases that share genetic background with many other diseases through strong associations. We collect a large number of known disease-gene associations (DGAs) using DisGeNET in order to construct a bipartite disease-gene network. Subsequently, a weighted human disease network (WHDN) is built by connecting pairs of diseases that share associated genes and the edge weights reflect the number of genes they share as well as the strength of the DGAs. Then we develop a new vertex centrality measure for the WHDN, degree and importance of link centrality (DIL-W), which considers both the degree of a vertex and the importance of its incident edges in weighted graphs. Our network-based analysis methods are shown to be able to identify more important diseases comparing to degree, closeness and betweenness centralities. The identified disease-disease correlations include previous knowledge supported by published literature as well as less known and novel correlations that can be valuable for further studies.\nThe contribution of our study is two fold, the construction of the WHDN and the importance measurement of a vertex considering both its degree and edge(s). First, comparing to the HDN (an un-weighted graph) proposed by Goh et al. [ 33 ], the mechanism of including vertices and edges is the same, but we add the consideration of the confidence and strength of disease-disease correlations and add weights to edges of the HDN. Such a WHDN allows us to prune the network using a vertex disparity filter [ 50 ], which considerably reduces the complexity of the network by removing less-significant edges (from 112,324 to 38,275 before and after the back-bone extraction), while preserving most of the vertices (from 5,278 to 4,898, respectively).\nSecond, we further extend a new vertex centrality measure DIL-W for the WHDN, which quantifies the importance of a vertex by considering its degree and the aggregative importance of its attached edge(s), with the inspiration that a disease should be considered important if it is correlated with many other diseases (i.e., its degree) and these correlations are themselves strong and significant (i.e., edge importance).\nDIL-W only uses local information of a vertex for its importance assessment, and its computational complexity is  O ( | V | k ¯ 2 ) , where | V | is the total number of vertices and  k ¯  is the average degree of vertices in a network. Thus, DIL-W can be efficient to compute for large and sparse networks.\nUpon application to the WHDN, DIL-W is shown to outperform three commonly used centrality measures, degree, closeness and betweenness, and has identified top diseases including  epilepsy ,  anemia , and  obesity .  Table 1  shows the degree in the WHDN and the most correlated disease of those 30 top-rank diseases. We are also able to find previous publications that verify almost all the correlations of those pairs of diseases, shown as references in the table. Besides some very well-known correlations such as  heart failure — obesity  and  diabetes — obesity , the table also reports some less known but interesting correlations. For instance, Savin [ 58 ] showed that  atypical retinitis pigmentosa  is correlated with  obesity . Moreover, the correlation between  anemia  and  pediatric failure to thrive  had not been reported in the literature until recently Dimmock  et al . [ 57 ] suggested  anemia  as one of the novel causes of  failure to thrive  in children. Zimmerman [ 61 ] studied the cause of different types of  cirrhosis  resulting from different drug-induced injuries. This supports our finding on the correlation between  cirrhosis  and  chemical and drug induced liver injury .\nThe disease-gene associations come from DisGeNet [ 49 ] only. While this is a valuable resource, it is merely one of the many databases that have disease gene information (including Jensen Lab’s DISEASES [ 80 ] and DiseaseConnect [ 81 ] databases), all of which have their own disease association scoring convention. The alternative databases will be explored in our future study.\nAnother future direction we would like to explore is to implement our proposed centrality measure DIL-W for other networks and to test its utility. Centrality measures essentially tell us how important a vertex is in the context of a network structure, and this “importance” can take different meanings in various types of networks. For instance, in Internet, vertices are physical routers, servers, and computers that are responsible for information transportation, therefore, vertex importance should reflect how much a vertex controls and its remove influences the traffic flow. We expect DIL-W to find useful venues for  weighted  networks that consider vertices as important when they are connected with many others through strong relationships.\nOur understanding of human diseases is still largely unclear and the disease-gene associations are far from being complete. Future studies could explore the utilization of multiple types of data and more powerful computational tools to better cluster and categorize human diseases and to predict new genes and other factors that can explain diseases.\n\n(EPS)\nClick here for additional data file.","source_license":"CC-BY-4.0","license_restricted":false}