BioMed Research International

Volume 2013 (2013), Article ID 523415, 8 pages

http://dx.doi.org/10.1155/2013/523415

Research Article

Identification of Age-Related Macular Degeneration Related Genes by Applying Shortest Path Algorithm in Protein-Protein Interaction Network

Jian Zhang,¹ Min Jiang,² Fei Yuan,² Kai-Yan Feng,³ Yu-Dong Cai,⁴ Xun Xu,¹ and Lei Chen⁵

¹Department of Ophthalmology, Shanghai First People’s Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai 200080, China

²State Key Laboratory of Medical Genomics, Institute of Health Sciences, Shanghai Jiaotong University School of Medicine and Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200025, China

³Beijing Genomics Institute, Shenzhen Beishan Industrial Zone, Shenzhen 518083, China

⁴Institute of Systems Biology, Shanghai University, Shanghai 200444, China

⁵College of Information Engineering, Shanghai Maritime University, Shanghai 201306, China

Received 22 October 2013; Accepted 27 November 2013

Academic Editor: Tao Huang

Copyright © 2013 Jian Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This study attempted to find novel age-related macular degeneration (AMD) related genes based on 36 known AMD genes. The well-known shortest path algorithm, Dijkstra’s algorithm, was applied to find the shortest path connecting each pair of known AMD related genes in protein-protein interaction (PPI) network. The genes occurring in any shortest path were considered as candidate AMD related genes. As a result, 125 novel AMD genes were predicted. The further analysis based on betweenness and permutation test indicates that there are 10 genes involved in the formation or development of AMD and may be the actual AMD related genes with high probability. We hope that this contribution would promote the study of age-related macular degeneration and discovery of novel effective treatments.

1. Introduction

Macular degeneration generally refers to age-related macular degeneration (AMD or ARMD). It is a major cause of blindness and visual impairment in older individuals (>50 years) in Western countries [1, 2]. In China, the disease incidence is 6.04%–11.19% among age 60–69 and rises with the increase of aging population [3]. Its main symptom is central visual damage. Known as one of the hardest eye diseases to treat, the mechanism underlying AMD has not been clear yet.

However, many risk factors have been identified related to AMD occurrence, such as old age, smoking [4], hypertension, oxidative stress [5], and high-energy visible light [6, 7]. Family history and gene mutations are genetic elements for AMD occurrence. The lifetime risk of developing late-stage macular degeneration is 50% for people who have a relative with macular degeneration, versus 12% for people who do not have relatives with macular degeneration [8, 9].

As a complex disease, identification of disease-related genes is prerequisite and persistent. Through previous genetic analysis and experimental validation, numerous genomic loci and a large number of candidate genes have been shown to involve in AMD [10, 11]. However, they just account for part of the AMD pathogenesis. The majority of AMD cases cannot be explained through these known mechanisms. Recently, through meta-analysis of genome-wide association studies (GWAS) for advanced AMD, it was estimated that currently identified loci account for approximately 55% of the heritability of advanced AMD [12]. So, to identify novel genes involved in the formation and development of AMD is necessary and will improve our knowledge about additional pathways and pathological mechanisms of the disease, as well as methods for more effective treatments.

Since the number of human genes is very huge, it is impossible to screen them to discover novel AMD related genes by experiment alone. Computational methods give an alternative way to help scientists pick out genes that are AMD related genes with higher probability than others. On the other hand, computational methods have been successfully used to tackling various problems in many biological areas, such as protein attributions prediction [13–17], drug design [18–22], and analysis of complicated biological network [23]. In this study, we proposed a novel computational method to identify novel AMD related genes. Based on the current known AMD related genes, retrieved from Retina International and some previous published documents, some novel candidate AMD related genes were discovered by applying Dijkstra’s algorithm [24] in the protein-protein interaction network. Further GO and KEGG pathway analysis indicates that some candidate genes are involved in the formation and development of AMD. Thus, they may be the actual AMD related genes with high possibility, which should be paid more attention by scientists. It is hopeful that our contribution would help to uncover the mechanism of this disease and discover novel effective treatments.

2. Materials and Methods

2.1. Known AMD Related Genes

The AMD related genes were collected from the following resources: Sixteen genes are found in Retina International at the website http://www.retina-international.org/files/sci-news/remacdy.htm (recent update from March 24, 2010), with “age-related macular dystrophy” selected as key words. The genes for the complement system proteins factor H (CFH), factor B (CFB), and factor 3 (C3) are strongly associated with a person’s risk for developing AMD. HTRA1 (encoding a secreted serine protease), which has implications for the disease, was identified in 2006 [25, 26]. Mutational analysis of ABCR in juvenile macular degeneration (STGD) families revealed a total of 19 different mutations including homozygous mutations in two families with consanguineous parentage. These data indicate that ABCR is the causal gene of STGD [27]. Deletion of the complement factor H-related genes CFHR3 and CFHR1 protects against AMD [28, 29]. The AMD gene consortium identified 19 loci with associations reaching < 5 × 10⁻⁸, including seven new loci after genome-wide association studies (GWAS) and meta-analysis. The 23 nearby genes were selected [30]. Finally, we obtained 36 known AMD related genes after integration of all resources, which are listed in Table 1.

Table 1: AMD related genes.

Index
Gene name
Index
Gene name

1
ABCR
2
ADAMTS9

3
APOE
4
ARMS2

5
B3GALTL
6
C2

7
C3
8
CETP

9
CFB
10
CFH

11
CFH1
12
CFHR1

13
CFHR3
14
CFI

15
COL10A1
16
COL8A1

17
DDR1
18
ERCC6

19
FBLN5
20
FBLN6

21
FILIP1L
22
HMCN1

23
HTRA1
24
IER3

25
LIPC
26
LOC387715

27
QRX
28
RAD51B

29
RAXL1
30
SLC16A8

31
TGFBR1
32
TIMP3

33
TLR3
34
TLR4

35
TNFRSF10A
36
VEGFA

2.2. Protein-Protein Interaction (PPI) Network

The PPI network was constructed according to the information retrieved from the well-known database STRING (Search Tool for the Retrieval of Interacting Genes/Proteins, http://string.embl.de/) (version 9.0) [31], a large database containing direct (physical) and indirect (functional) interactions. In the database, each interaction consists of two proteins and a score quantifying the likelihood that an interaction may occur. For convenience, let denote the score of the interaction between two proteins and . The constructed network took proteins as its nodes, and the edge between any two nodes existed if and only if the score of the interaction between the corresponding proteins was greater than zero. To reflect the difference of interactions, each edge with endpoints and in the network was labeled an edge weight defined as follows:

where and were corresponding proteins of nodes and , respectively. Some previous studies have shown that proteins with an edge connecting them in the network, known as interactive proteins, are more likely to share common features than those without edge connecting them [14–16, 32]. Accordingly, it can be deduced from this fact that proteins with small distance in the network may share similar features.

2.3. Calculation of Betweenness Using Dijkstra’s Algorithm

As described in Section 2.1, 36 genes were collected, which are related to AMD and must have some common features related to AMD. According to Section 2.2, it can be deduced that the proteins in the shortest path connecting any pair of actual AMD related genes may all share some common features related to AMD. To obtain the shortest path, the Dijkstra’s algorithm, a well-known shortest path algorithm conceived by Dijkstra in 1956, was employed [24].

For the given node in a network, its betweenness is a value that is related to the number of the shortest paths containing it as an inner node among all shortest paths connecting all pair of nodes [33]. Since betweenness accounts for direct and indirect influences of proteins at distant network [34], it has been used for studying various natural and man-made networks [33, 35–38]. Here, we used the concept of betweenness to identify novel AMD related genes from known 36 genes. Compared to the original concept of betweenness, we only considered the shortest paths, finding by Dijkstra’s algorithm in PPI network, which connected two known AMD related genes instead of all shortest paths. The betweenness of a node in this study is defined as the number of shortest paths containing as an inner node among these shortest paths.

2.4. Identification of Related Genes Based on Betweenness and Permutation Test

The betweenness of a node/gene in this study means the number of shortest paths that connect all pairs of known AMD related genes/nodes and contain the node/gene as an inner node. Therefore, it is possible that genes/nodes with high betweenness may share more features related to AMD than those with low betweenness. As for proteins with betweenness equal to 0, we considered the likelihood that they are the novel AMD related genes to be zero. Accordingly, we picked out proteins with betweenness greater than 0 and they are termed as shortest path genes.

The betweenness of some nodes/genes may be influenced by the essential structure of the network. For example, the cut-vertex of the network may always receive high betweenness regardless of the distribution of known nodes/genes. To avoid this situation, a permutation test was conducted to further screen the shortest path genes. We randomly selected the same number of nodes/genes as the actual number of AMD related genes from PPI network 500 times and recalculated the shortest paths between these randomly selected genes. The permutation FDR of the shortest path genes was defined as

where was the number of shortest paths among actual AMD related genes across shortest path gene , was the number of shortest paths among randomly selected genes across shortest path gene , and was the count of times when was greater than . According to Jiang et al.’s work [38], smaller permutation FDR of one shortest path gene indicates that it is the actual AMD related gene with high possibility.

2.5. KEGG and GO Enrichment Analysis

Functional annotation tool of DAVID [39] was used for KEGG pathway and GO enrichment analysis. The enrichment value was corrected to control family-wide false discovery rate under certain rate (e.g., ≤0.05) with Benjamin multiple testing correction method [40]. All the genes in the human genome were selected as background during the enrichment analysis. There are 13 items in the output of DAVID; their meanings are listed as below:

(1)category: DAVID category, that is, KEGG or GO;(2)term: gene set name;(3)count: the number of genes associated with this gene set;(4)percentage: calculated by “gene associated with this gene set”/“total number of query genes;”(5) value: modified Fisher exact value;(6)genes: the list of genes from your query set that are annotated to this gene set;(7)list total: the number of genes in your query list mapped to any gene set in this ontology;(8)pop hits: the number of genes annotated to this gene set on the background list;(9)pop total: the number of genes on the background list mapped to any gene set in this ontology;(10)fold enrichment: the ratio of the proportions on query genes and the background information which are associated with the gene set;(11)Bonferroni: Bonferroni adjusted value;(12)Benjamini: Benjamini adjusted value;(13)FDR: FDR adjusted value.

3. Results and Discussion

3.1. Candidate Genes Filtered by Betweenness and Permutation Test

For the 36 known AMD related genes, we searched the shortest path connecting any pair of them. After counting the betweenness of inner nodes in these paths, we obtained 168 shortest path genes with betweenness greater than 0, where 4 are known AMD related genes. These 168 genes are available in Supplementary Material I available online at http://dx.doi.org/10.1155/2013/523415. To further screen these genes, the permutation test was conducted and permutation FDRs (refer to (2)) of these genes were calculated, which are also listed in Supplementary Material I. It can be seen that 125 genes were with permutation FDRs less than 0.01, which are considered to be greatly related to AMD.

3.2. Results of Functional Annotation Tool

Functional annotation tool of DAVID was used to analyze the 125 genes based on KEGG pathways and GO terms. The GO and KEGG analysis results can be found in Supplementary Materials II and III, respectively. The following sections give the detailed discussion based on these results.

3.3. GO Enrichment Analysis

It is observed from Supplementary Material II that 865 GO terms were enriched by the 125 genes. The top 1% of GO terms (totally 9) sorted by “Count” item (refer to Figure 1 for the “Count” item of these GO terms) were investigated and discussed as below.

Figure 1: The top 9 GO terms shared by 125 genes. The -axis lists GO’s ID and name, while the -axis represents the number of genes shared the GO term among the 125 genes.

Among these 9 GO terms, three of them are cellular component (CC) GO terms and six biological process (BP) GO terms. The three CC terms included: GO:0005886: plasma membrane (“count” = 60), GO:0044459: plasma membrane part (“count” = 45), and GO:0005829: cytosol (“count” = 40). Early AMD usually takes place accompanied by accumulation of drusen, a kind of cellular debris, between the retina and the choroid [41]. It was speculated that exosomes produced by cellular lysosome formed the drusen [42]. They carried many extra- or intracellular proteins, lipids, and cellular components wrapped by or included plasma membrane, and gathered other extracellular factors to deposit as a accumulative core [43]. These indicated the role of cytosol parts and plasma membrane in the pathogenesis of AMD in the early stage.

The six BP terms were GO:0007166: cell surface receptor linked signal transduction (“count” = 49), GO:0007242: intracellular signaling cascade (“count” = 42), GO:0010604: positive regulation of macromolecule metabolic process (“count” = 37), GO:0042981: regulation of apoptosis (“count” = 36), GO:043067: regulation of programmed cell death (“count” = 36), and GO:0010941: regulation of cell death (“count” = 36). The first two may represent all the processes of signal transduction. The macromolecule metabolic process was relatively critical in the development of AMD. As mentioned above, metabolism of proteins and other biological macromolecules cannot go without lysosomes. And if that does not work, the exosomes carrying these undigested molecules would induce the formation of drusen or further cause increased local immune responses and promote the AMD. The next three were all related with cell death. It was reported that necrosis can be responsible for the cell loss in the retinal degeneration [44]. Immune responses can protect cells from death and kill them when they were impaired. Advanced AMD in mice model involved death of photoreceptor cells [45]. These genes may not be the causative factors in AMD, but what they have in common would give us a suggestion about how the AMD happens or what problems may occur in the development.

3.4. KEGG Pathway Enrichment Analysis

It can be observed from Supplementary Material III that there were 56 pathways which were significantly enriched in by the 125 genes potentially being related with AMD pathogenesis. The KEGG pathways with “count” item at least 14 were investigated. We found that 5 of them are related to AMD. Figure 2 shows the “count” item of these pathways. In details, T-cell receptor signaling pathway (hsa04660) and neurotrophin signaling pathway (hsa04722) ranked in the top five, covered 22 genes and 21 genes (refer to Figure 2), respectively. Seventeen genes (“percentage” = 13.6%) enriched in the Chemokine signaling pathway (hsa04062), 15 genes (“percentage” = 12%) in the B cell receptor signaling pathway (hsa04662), and 14 genes (“percentage” = 11.2%) in the Toll-like receptor signaling pathway (hsa04620). Although the exact mechanism of AMD is unknown, decades of studies have highlighted the critical role of immune processes in the development, progression, and treatment of AMD [46]. Advanced AMD is of great relevance for the retinal neovascularization and atrophic lesions involving macrophages, lymphocytes, microglia, and mast cells, as well as fibroblasts [47, 48].

Figure 2: The 5 KEGG pathways shared by 125 genes. The -axis lists pathway’s ID and name, while the -axis represents the number of genes shared the pathway among the 125 genes.

Undeniably, all immune responses have dual, opposing roles, normally preventing the host while promoting disease under disorder conditions. The Chemokine signaling pathway is initiated by cytokine-cytokine receptor interactions. It was found that animal models with CC-chemokine ligand 2 (CCL2), CC-chemokine receptor 2 (CCR2), and/or CX3C-chemokine receptor 1 (CX3CR1) deficiency had AMD-like features, including thickening Bruch membrane, deposition of drusen, retinal pigmented epithelium (RPE) hypertrophy, and retinal thinning [49, 50]. T cell receptor signaling pathway and Toll-like receptor signaling pathway are both activated to release immunologic factors, such as IL-1β, IL-2, and IL-18. RPE cells can release IL-1β in responding to kinds of chemokines which may be dangerous signals and were abundant in AMD [51–53]. While IL- 1β and Toll-like receptor 2 (TLR2) ligands can promote neovascularization, activating TLR3 would decrease choroidal neovascularization and promote RPE degeneration [54–56].

B cells do not show a direct role in mediating angiogenesis or tissue damage in AMD pathogenesis, but it cannot be excluded in which there is some relationship between B cell immune and AMD. For example, the Ccl2−/− Cx3cr1−/− mouse model was found having retina-specific autoantibodies [49, 50]. The neurotrophin signaling pathway involves many nerve growth factors (NGFs) and neurotrophins (NT3, NT4). Microenvironment balance of the eye would guarantee the health of eye organisms. It was successful to relieve patients with AMD with NGF eye drop, which suggested the lack of NGF may be a cause of the AMD symptoms and may involve in the disease progress [57]. All the above indicated the relationship of pathways with AMD and may expand the avenues to explore new mechanisms in AMD occurrence and development.

3.5. Analysis of 10 Significant AMD Related Genes

There were 12 genes (listed in Table 2) whose betweenness was larger than 100 among the 125 shortest path genes, where 1 gene is the known AMD related gene. Among the rest 11 genes, we found that 10 of them are related to the formation or development of AMD.

Table 2: Shortest path genes with betweenness greater than 100.

Ensemble ID of genes
Gene name
Betweenness

ENSP00000245907a
C3
336

ENSP00000313419
CD19
319

ENSP00000356024
CR2
319

ENSP00000264033
CBL
295

ENSP00000275493
EGFR
195

ENSP00000344456
CTNNB1
172

ENSP00000344818
UBC
158

ENSP00000269305
TP53
114

ENSP00000415941
C4B
111

ENSP00000284981
APP
104

ENSP00000326366
PSEN1
104

ENSP00000356016
CR1
101

It is the known AMD related gene.

Genetic studies had implicated immune responses especially the complement system in AMD pathogenesis and development [43, 58, 59]. The genes for the complement system proteins factor H (CFH), factor B (CFB), and factor 3 (C3) are strongly associated with AMD susceptibility. For instance, the Y402H variant in CFH gene is significantly associated with AMD in Asian populations [60]. C3 can protect the aging retina by complement system [61]. There were two complement component receptors (CR1 and CR2) genes in our result. As c3b/c4b or c3d receptors, the signaling was definitely disturbed if they were dysfunctional. CD19 was a B-lymphocyte antigen. It was reported that the RPE can inhibit the B-cell activation [62], maybe through the B-cell surface antigens, such as cd19.

PSEN1, also known as PS1, was an Alzheimer’s disease (AD) gene. The transgenic mice with this gene can express AMD features under some conditions [63], which suggested some internal relationship between the gene and AMD. APP was also an AD gene and would be involved in AMD too. The UBC (ubiquitin C) belonged to the ubiquitin-proteasome system; one terminal hydrolase of which was closely related with photoreceptor cell apoptosis [64]. CBL (also named E3 ubiquitin protein ligase), EGFR, and P53 were both tumor related genes, so their ordinary condition was all necessary for cells growth or organ development, especially the photoreceptor cells’ development [65–67]. The Wnt/β-Catenin Signaling had been involved in retina development in many studies, and its activities can protect photoreceptors from damage and take part in the retina cell regeneration [68, 69]. Although these genes did not show directly relationship with AMD, their roles indicated some connection with the disease and that would be a cue for further researches to find out more AMD related factors.

4. Conclusion

Identification of disease genes is one of the most important and challenging problems in biomedicine and genomics. Since AMD has a high risk factor which results in a loss of vision for elder adults, it is eager to find novel AMD related genes, thereby understanding its mechanism well and discovering effective treatments. This study proposed a computation method to identify AMD related genes by applying the shortest path algorithm in the PPI network. Through the GO and KEGG pathway analysis of the genes identified by our method, they are significantly enriched for GO terms and KEGG pathways related to AMD, which implies they have direct or indirect relationship with the formation or development of AMD. It is hopeful that this contribution would promote the study of this disease.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This paper is supported by the National Science and Technology Major Projects of the Twelfth Five-Year Plan (no. 2011ZX09302-007-02), National Basic Research Program of China (no. 2011CB510102 and no. 2011CB510101), National Natural Science Foundation of China (no. 81273424, no. 81170862, no. 31371335, no. 61202021, and no. 61373028), Innovation Program of Shanghai Municipal Education Commission (no. 12ZZ087 and no. 12YZ120), Shanghai Educational Development Foundation (no. 12CG55), and grant of “The First-class Discipline of Universities in Shanghai”.

References

K. M. Gehrs, D. H. Anderson, L. V. Johnson, and G. S. Hageman, “Age-related macular degeneration—emerging pathogenetic and therapeutic concepts,” Annals of Medicine, vol. 38, no. 7, pp. 450–471, 2006.

R. Klein, C.-F. Chou, B. E. K. Klein, X. Zhang, S. M. Meuer, and J. B. Saaddine, “Prevalence of age-related macular degeneration in the US population,” Archives of Ophthalmology, vol. 129, no. 1, pp. 75–80, 2011.

X. Tong, Macular Degeneration, Baidu Baike, 2013.

C. M. Cheung and T. Y. Wong, “Treatment of age-related macular degeneration,” The Lancet, vol. 382, no. 9900, pp. 1230–1232, 2013.

R. Sarangarajan and S. P. Apte, “Melanin aggregation and polymerization: possible implications in age-related macular degeneration,” Ophthalmic Research, vol. 37, no. 3, pp. 136–141, 2005.

C. Glazer-Hockstein and J. L. Dunaief, “Could blue light-blocking lenses decrease the risk of age-related macular degeneration?” Retina, vol. 26, no. 1, pp. 1–4, 2006.

T. H. Margrain, M. Boulton, J. Marshall, and D. H. Sliney, “Do blue light filters confer protection against age-related macular degeneration?” Progress in Retinal and Eye Research, vol. 23, no. 5, pp. 523–531, 2004.

M. B. Gorin, “Genetic insights into age-related macular degeneration: controversies addressing risk, causality, and therapeutics,” Molecular Aspects of Medicine, vol. 33, pp. 467–486, 2012.

U. Chakravarthy, T. Y. Wong, A. Fletcher et al., “Clinical risk factors for age-related macular degeneration: a systematic review and meta-analysis,” BMC Ophthalmology, vol. 10, no. 1, article 31, 2010.

A. Y. C. Ting, T. K. M. Lee, and I. M. MacDonald, “Genetics of age-related macular degeneration,” Current Opinion in Ophthalmology, vol. 20, no. 5, pp. 369–376, 2009.

K. Antoniak, W. Bienias, and J. Z. Nowak, “Age-related macular degeneration—A complex genetic disease,” Klinika Oczna, vol. 110, no. 4–6, pp. 211–218, 2008.

Y. Yu, T. R. Bhangale, J. Fagerness et al., “Common variants near FRK/COL10A1 and VEGFA are associated with advanced age-related macular degeneration,” Human Molecular Genetics, vol. 20, no. 18, pp. 3699–3709, 2011.

B. Niu, Y. D. Cai, W. C. Lu, G.-Z. Li, and K.-C. Chou, “Predicting protein structural class with AdaBoost Learner,” Protein and Peptide Letters, vol. 13, no. 5, pp. 489–492, 2006.

P. Gao, Q. P. Wang, L. Chen, and T. Huang, “Prediction of human genes regulatory functions based on proteinprotein interaction network,” Protein and Peptide Letters, vol. 19, pp. 910–916, 2012.

L. Hu, T. Huang, X. Shi, W.-C. Lu, Y. D. Cai, and K.-C. Chou, “Predicting functions of proteins in mouse based on weighted protein-protein interaction network and protein hybrid properties,” PLoS ONE, vol. 6, no. 1, Article ID e14556, 2011.

L. Hu, T. Huang, X.-J. Liu, and Y. D. Cai, “Predicting protein phenotypes based on protein-protein interaction network,” PLoS ONE, vol. 6, no. 3, Article ID e17668, 2011.

J. Cheng, J. Li, Z. Wang, J. Eickholt, and X. Deng, “The MULTICOM toolbox for protein structure prediction,” BMC Bioinformatics, vol. 13, article 65, 2012.

L. Chen, W. M. Zeng, Y.D. Cai, K.-Y. Feng, and K.-C. Chou, “Predicting anatomical therapeutic chemical (ATC) classification of drugs by integrating chemical-chemical interactions and similarities,” PLoS ONE, vol. 7, no. 4, Article ID e35254, 2012.

Y. Yamanishi, M. Araki, A. Gutteridge, W. Honda, and M. Kanehisa, “Prediction of drug-target interaction networks from the integration of chemical and genomic spaces,” Bioinformatics, vol. 24, no. 13, pp. i232–i240, 2008.

Y. Yamanishi, M. Kotera, M. Kanehisa, and S. Goto, “Drug-target interaction prediction from chemical, genomic and pharmacological data in an integrated framework,” Bioinformatics, vol. 26, no. 12, pp. i246–i254, 2010.

L. Chen, T. Huang, J. Zhang, et al., “Predicting drugs side effects based on chemical-chemical interactions and protein-chemical interactions,” BioMed Research International, vol. 2013, Article ID 485034, 8 pages, 2013.

L. Chen, J. Lu, X. Luo, and K.-Y. Feng, “Prediction of drug target groups based on chemical-chemical similarities and chemical-chemical/protein connections,” Biochimica et Biophysica Acta, 2013.

L. Chen, W. M. Zeng, Y. D. Cai, and T. Huang, “Prediction of metabolic pathway using graph property, chemical functional group and chemical structural set,” Current Bioinformatics, vol. 8, pp. 200–207, 2013.

T. H. Gormen, C. E. Leiserson, R. L. Rivest, and C. Stein, Eds., Introduction to Algorithms, MIT Press, Cambridge, Mass, USA, 1990.

Z. Yang, N. J. Camp, H. Sun et al., “A variant of the HTRA1 gene increases susceptibility to age-related macular degeneration,” Science, vol. 314, no. 5801, pp. 992–993, 2006.

A. DeWan, M. Liu, S. Hartman et al., “HTRA1 promoter polymorphism in wet age-related macular degeneration,” Science, vol. 314, no. 5801, pp. 989–992, 2006.

T. P. Dryja, C. E. Briggs, E. L. Berson, P. J. Rosenfeld, and M. Abitbol, “ABCR gene and age-related macular degeneration,” Science, vol. 279, pp. 1107–1107, 1998.

A. E. Hughes, N. Orr, H. Esfandiary, M. Diaz-Torres, T. Goodship, and U. Chakravarthy, “A common CFH haplotype, with deletion of CFHR1 and CFHR3, is associated with lower risk of age-related macular degeneration,” Nature Genetics, vol. 38, no. 10, pp. 1173–1177, 2006.

L. G. Fritsche, N. Lauer, A. Hartmann et al., “An imbalance of human complement regulatory proteins CFHR1, CFHR3 and factor H influences risk for age-related macular degeneration (AMD),” Human Molecular Genetics, vol. 19, no. 23, pp. 4694–4704, 2010.

L. G. Fritsche, W. Chen, M. Schu, et al., “Seven new loci associated with age-related macular degeneration,” Nature Genetics, vol. 45, pp. 433–439, 2013.

L. J. Jensen, M. Kuhn, M. Stark et al., “STRING 8—a global view on proteins and their functional interactions in 630 organisms,” Nucleic Acids Research, vol. 37, no. 1, pp. D412–D416, 2009.

Y. F. Gao, L. Chen, Y. D. Cai, et al., “Predicting metabolic pathways of small molecules and enzymes based on interaction information of chemicals and proteins,” PLoS ONE, vol. 7, Article ID e45944, 2012.

J. Davis and M. Goadrich, “The relationship between precision-recall and ROC curves,” in Proceedings of the 23rd International Conference on Machine Learning (ICML '06), pp. 233–240, New York, NY, USA, June 2006.

J. B. M. Craven, Markov Networks for Detecting Overlapping Elements in Sequence Data, The MIT Press, Cambridge, Mass, USA, 2005.

R. Bunescu, R. Ge, R. J. Kate et al., “Comparative experiments on learning information extractors for proteins and their interactions,” Artificial Intelligence in Medicine, vol. 33, no. 2, pp. 139–155, 2005.

D. E. Johnson and G. H. I. Wolfgang, “Predicting human safety: screening and computational approaches,” Drug Discovery Today, vol. 5, no. 10, pp. 445–454, 2000.

B. Q. Li, B. Niu, L. Chen, et al., “Identifying chemicals with potential therapy of HIV based on protein-protein and protein-chemical interaction network,” PLoS ONE, vol. 8, Article ID e65207, 2013.

M. Jiang, Y. Chen, Y. Zhang, et al., “Identification of hepatocellular carcinoma related genes with k-th shortest paths in a protein-protein interaction network,” Molecular BioSystems, vol. 9, pp. 2720–2728, 2013.

D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.

Y. Benjamini and D. Yekutieli, “The control of the false discovery rate in multiple testing under dependency,” Annals of Statistics, vol. 29, no. 4, pp. 1165–1188, 2001.

G. S. Hageman, P. J. Luthert, N. H. Victor Chong, L. V. Johnson, D. H. Anderson, and R. F. Mullins, “An integrated hypothesis that considers drusen as biomarkers of immune-mediated processes at the RPE-Bruch's membrane interface in aging and age-related macular degeneration,” Progress in Retinal and Eye Research, vol. 20, no. 6, pp. 705–732, 2001.

A. L. Wang, T. J. Lukas, M. Yuan, N. Du, M. O. Tso, and A. H. Neufeld, “Autophagy and exomoses in the aged retinal pigment epithelium: possible relevance to drusen formation and age-related macular degeneration,” PLoS ONE, vol. 4, no. 1, Article ID e4160, 2009.

D. H. Anderson, R. F. Mullins, G. S. Hageman, and L. V. Johnson, “A role for local inflammation in the formation of drusen in the aging eye,” The American Journal of Ophthalmology, vol. 134, no. 3, pp. 411–431, 2002.

Y. Murakami, H. Matsumoto, M. Roh, et al., “Programmed necrosis, not apoptosis, is a key mediator of cell loss and DAMP-mediated inflammation in dsRNA-induced retinal degeneration,” Cell Death and Differentiation, 2013.

M. Chen, J. V. Forrester, and H. Xu, “Dysregulation in retinal para-inflammation and age-related retinal degeneration in CCL2 or CCR2 deficient mice,” PLoS ONE, vol. 6, no. 8, Article ID e22818, 2011.

J. Ambati, J. P. Atkinson, and B. D. Gelfand, “Immunology of age-related macular degeneration,” Nature Reviews Immunology, vol. 13, pp. 438–451, 2013.

S. M. Whitcup, A. Sodhi, J. P. Atkinson, et al., “The role of the immune response in age-related macular degeneration,” International Journal of Inflammation, vol. 2013, Article ID 348092, 10 pages, 2013.

P. L. Penfold, M. C. Madigan, M. C. Gillies, and J. M. Provis, “Immunological and aetiological aspects of macular degeneration,” Progress in Retinal and Eye Research, vol. 20, no. 3, pp. 385–414, 2001.

J. Ambati, A. Anand, S. Fernandez et al., “An animal model of age-related macular degeneration in senescent Ccl-2- or Ccr-2-deficient mice,” Nature Medicine, vol. 9, no. 11, pp. 1390–1397, 2003.

R. J. Ross, M. Zhou, D. Shen et al., “Immunological protein expression profile in Ccl2/Cx3cr1 deficient mice with lesions similar to age-related macular degeneration,” Experimental Eye Research, vol. 86, no. 4, pp. 675–683, 2008.

W. A. Tseng, T. Thein, K. Kinnunen, et al., “NLRP3 inflammasome activation in retinal pigment epithelial cells by lysosomal destabilization: implications for age-related macular degeneration,” Investigative Ophthalmology and Visual Science, vol. 54, pp. 110–120, 2013.

A. Kauppinen, H. Niskanen, T. Suuronen, et al., “Oxidative stress activates NLRP3 inflammasomes in ARPE-19 cells—implications for age-related macular degeneration (AMD),” Immunology Letters, vol. 147, pp. 29–33, 2012.

R. T. Liu, J. Gao, S. Cao, et al., “Inflammatory mediators induced by amyloid-beta in the retina and RPE in vivo: implications for inflammasome activation in age-related macular degeneration,” Investigative Ophthalmology and Visual Science, vol. 54, pp. 2225–2237, 2013.

M. E. Kleinman, K. Yamada, A. Takeda et al., “Sequence-and target-independent angiogenesis suppression by siRNA via TLR3,” Nature, vol. 452, no. 7187, pp. 591–597, 2008.

M. E. Kleinman, H. Kaneko, W. G. Cho et al., “Short-interfering RNAs induce retinal degeneration via TLR3 and IRF3,” Molecular Therapy, vol. 20, no. 1, pp. 101–108, 2012.

S. Lavalette, W. Raoul, M. Houssier et al., “Interleukin-1β inhibition prevents choroidal neovascularization and does not exacerbate photoreceptor degeneration,” The American Journal of Pathology, vol. 178, no. 5, pp. 2416–2423, 2011.

A. Lambiase, M. Coassin, P. Tirassa, F. Mantelli, and L. Aloe, “Nerve growth factor eye drops improve visual acuity and electrofunctional activity in Age-related macular degeneration: a case report,” Annali dell'Istituto Superiore di Sanita, vol. 45, no. 4, pp. 439–442, 2009.

C. Baudouin, G. A. Peyman, D. Fredj-Reygrobellet et al., “Immunohistological study of subretinal membranes in age-related macular degeneration,” Japanese Journal of Ophthalmology, vol. 36, no. 4, pp. 443–451, 1992.

D. H. Anderson, M. J. Radeke, N. B. Gallo et al., “The pivotal role of the complement system in aging and age-related macular degeneration: hypothesis re-visited,” Progress in Retinal and Eye Research, vol. 29, no. 2, pp. 95–112, 2010.

N. Kondo, H. Bessho, S. Honda, and A. Negi, “Complement factor H Y402H variant and risk of age-related macular degeneration in Asians: a systematic review and meta-analysis,” Ophthalmology, vol. 118, no. 2, pp. 339–344, 2011.

J. Hoh Kam, E. Lenassi, T. H. Malik, M. C. Pickering, and G. Jeffery, “Complement component c3 plays a critical role in protecting the aging retina in a murine model of age-related macular degeneration,” The American Journal of Pathology, vol. 183, pp. 480–492, 2013.

S. Sugita, S. Horie, Y. Yamada, and M. Mochizuki, “Inhibition of B-Cell activation by retinal pigment epithelium,” Investigative Ophthalmology and Visual Science, vol. 51, no. 11, pp. 5783–5788, 2010.

Z. Dong, J. Li, Y. Leng et al., “Cyclic intensive light exposure induces retinal lesions similar to age-related macular degeneration in APPswe/PS1 bigenic mice,” BMC Neuroscience, vol. 13, article 34, 2012.

Y. Sano, A. Furuta, R. Setsuie et al., “Photoreceptor cell apoptosis in the retinal degeneration of Uchl3-deficient mice,” The American Journal of Pathology, vol. 169, no. 1, pp. 132–141, 2006.

J. Wan, R. Ramachandran, and D. Goldman, “HB-EGF is necessary and sufficient for Müller glia dedifferentiation and retina regeneration,” Developmental Cell, vol. 22, no. 2, pp. 334–347, 2012.

L. Vuong, D. E. Brobst, I. Ivanovic, D. M. Sherry, and M. R. Al-Ubaidi, “p53 selectively regulates developmental apoptosis of rod photoreceptors,” PLoS ONE, vol. 8, Article ID e67381, 2013.

H. Meisner, A. Daga, J. Buxton et al., “Interactions of Drosophila Cbl with epidermal growth factor receptors and role of Cbl in R7 photoreceptor cell development,” Molecular and Cellular Biology, vol. 17, no. 4, pp. 2217–2225, 1997.

B. M. Braunger, A. Ohlmann, M. Koch, et al., “Constitutive overexpression of Norrin activates Wnt/beta-catenin and endothelin-2 signaling to protect photoreceptors from light damage,” Neurobiology of Disease, vol. 50, pp. 1–12, 2013.

D. Sanges, N. Romo, G. Simonte, U. di Vicino, A. D. Tahoces, et al., “Wnt/beta-catenin signaling triggers neuron reprogramming and regeneration in the mouse retina,” Cell Reports, vol. 4, pp. 271–286, 2013.