In this study we investigated the advantage of including network information in prioritizing disease genes of type 1 diabetes (T1D). First, a naïve Bayesian network (NBN) model was developed to integrate information from multiple data sources and to define a T1D-involvement probability score (PS) for each individual gene. The algorithm was validated using known functional candidate genes as a benchmark. Genes with higher PS were found to be more likely to appear in T1D-related publications. Next a new network activity metric was proposed to evaluate the T1D relevance of protein-protein interaction (PPI) subnetworks. The metric considered the contribution both from individual genes and from network topological characteristics. The predictions were confirmed by several independent datasets, including a genome wide association study (GWAS), and two large-scale human gene expression studies. We found that novel candidate genes in the T1D subnetworks showed more significant associations with T1D than genes predicted using PS alone. Interestingly, most novel candidates were not encoded within the human leukocyte antigen (HLA) region, and their expression levels showed correlation with disease only in cohorts with low-risk HLA genotypes. The results suggested the importance of mapping disease gene networks in dissecting the genetics of complex diseases, and offered a general approach to network-based disease gene prioritization from multiple data sources.
Get full access to this article
View all access options for this article.
References
1.
AboumradE., MadecA.M., ThivoletC.2007. The CXCR4/CXCL12 (SDF-1) signalling pathway protects non-obese diabetic mouse from autoimmune diabetes. Clin Exp Immunol, 148:432–439.
2.
AertsS., LambrechtsD., MaityS.et al.2006. Gene prioritization through genomic data fusion. Nat Biotechnol, 24:537–544.
3.
AraguesR., SanderC., OlivaB.2008. Predicting cancer involvement of genes from heterogeneous data. BMC Bioinformatics, 9:172.
4.
BaranziniS.E., GalweyN.W., WangJ.et al.2009. Pathway and network-based analysis of genome-wide association studies in multiple sclerosis. Hum Mol Genet, 18:2078–2090.
5.
BarrettJ.C., ClaytonD.G., ConcannonP.et al.2009. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet, 41:703–707.
6.
BotsteinD., RischN.2003. Discovering genotypes underlying human phenotypes: past successes for mendelian disease, future approaches for complex disease. Nat Genet, 33,Suppl:228–237.
7.
CarterS.L., BrechbuhlerC.M., GriffinM., BondA.T.2004. Gene co-expression network topology provides a framework for molecular characterization of cellular state. Bioinformatics (Oxf, England), 20:2242–2250.
8.
ChuangH.Y., LeeE., LiuY.T., LeeD., IdekerT.2007. Network-based classification of breast cancer metastasis. Mol Systems Biol, 3:140.
9.
ColonnaL., CatalanoG., ChewC.et al.2010. Therapeutic targeting of Syk in autoimmune diabetes. J Immunol, 185:1532–1543.
10.
DaoP., ColakR., SalariR.et al.2010. Inferring cancer subnetwork markers using density-constrained biclustering. Bioinformatics, 26:i625–i631.
11.
DomingosP., PazzaniM.1997. On the optimality of the simple Bayesian classifier under zero-one loss. Machine Learning, 29:103–130.
12.
EavesI.A., WickerL.S., GhandourG.et al.2002. Combining mouse congenic strains and microarray gene expression analyses to study a complex trait: the NOD model of type 1 diabetes. Genome Res, 12:232–243.
13.
EnglishS.B., ButteA.J.2007. Evaluation and integration of 49 genome-wide experiments and the prediction of previously unknown obesity-related genes. Bioinformatics (Oxf, England), 23:2910–2917.
14.
FrankeL., BakelH., FokkensL., De JongE.D., Egmont-PetersenM., WijmengaC.2006. Reconstruction of a functional human gene network, with an application for prioritizing positional candidate genes. Am J Hum Genet, 78:1011–1025.
GaoS., WangX.2009. Predicting type 1 diabetes candidate genes using human protein-protein interaction networks. J Comput Sci Syst Biol, 2:133–146.
17.
GaoS., WangX.2011. Quantitative utilization of prior biological knowledge in the Bayesian network modeling of gene expression data. BMC Bioinformatics, 12:359.
GaultonK.J., WillerC.J., LiY.et al.2008. Comprehensive association study of type 2 diabetes and related quantitative traits with 222 candidate genes. Diabetes, 57:3136–3144.
20.
GeorgeR.A., LiuJ.Y., FengL.L., Bryson-RichardsonR.J., FatkinD., WoutersM.A.2006. Analysis of protein sequence and interaction data for candidate disease gene prediction. Nucleic Acids Res, 34:e130.
21.
GohK.I., CusickM.E., ValleD., ChildsB., VidalM., BarabasiA.L.2007. The human disease network. Proc Natl Acad Sci USA, 104:8685–8690.
22.
GustafsonT.A., HeW., CraparoA., SchaubC.D., O'NeillT.J.1995. Phosphotyrosine-dependent interaction of SHC and insulin receptor substrate 1 with the NPEY motif of the insulin receptor via a novel non-SH2 domain. Molec Cellular Biol, 15:2500–2508.
23.
HillN.J., StotlandA., SolomonM., SecrestP., GetzoffE., SarvetnickN.2007. Resistance of the target islet tissue to autoimmune destruction contributes to genetic susceptibility in type 1 diabetes. Biol Direct, 2:5.
24.
HungJ.H., WhitfieldT.W., YangT.H., HuZ., WengZ., DelisiC.2010. Identification of functional modules that correlate with phenotypic difference: the influence of network topology. Genome Biol, 11:R23.
25.
HunterS., ApweilerR., AttwoodT.K.et al.2009. InterPro: the integrative protein signature database. Nucleic Acids Res, 37:D211–D215.
26.
IdekerT., OzierO., SchwikowskiB., SiegelA.F.2002. Discovering regulatory and signalling circuits in molecular interaction networks. Bioinformatics (Oxf, England), 18,Suppl 1:S233–S240.
27.
Jimenez-SanchezG., ChildsB., ValleD.2001. Human disease genes. Nature, 409:853–855.
28.
JonesC.E., BrownA.L., BaumannU.2007. Estimating the annotation error rate of curated GO database sequence annotations. BMC Bioinformatics, 8:170.
29.
KaizerE.C., GlaserC.L., ChaussabelD., BanchereauJ., PascualV., WhiteP.C.2007. Gene expression in peripheral blood mononuclear cells from children with diabetes. J Clin Endocrinol Metab, 92:3705–3711.
30.
KannM.G.2007. Protein interactions and disease: computational approaches to uncover the etiology of diseases. Briefings Bioinformatics, 8:333–346.
31.
LageK., KarlbergE.O., StorlingZ.M.et al.2007. A human phenome-interactome network of protein complexes implicated in genetic disorders. Nat Biotech, 25:309–316.
32.
LiQ., XuB., MichieS.A., RubinsK.H., SchreriberR.D., McDevittH.O.2008. Interferon-α initiates type 1 diabetes in nonobese diabetic mice. Proc Natl Assn Sci USA, 105:12439–12444.
33.
LiuC.C., ChenW.S., LinC.C.et al.2006. Topology-based cancer classification and related pathway mining using microarray data. Nucleic Acids Res, 34:4069–4080.
34.
LiuwantaraD., ElliotM., SmithM.W.et al.2006. Nuclear factor-kappaB regulates beta-cell death: a critical role for A20 in beta-cell protection. Diabetes, 55:2491–2501.
35.
LozaM.J., McCallC.E., LiL., IsaacsW.B., XuJ., ChangB.L.2007. Assembly of inflammation-related genes for pathway-focused genetic analysis. PloS One, 2:e1035.
MaH., SchadtE.E., KaplanL.M., ZhaoH.2011. COSINE: COndition-SpecIfic sub-NEtwork identification using a global optimization method. Bioinformatics (Oxf, England), 27:1290–1298.
38.
MassaM.S., ChiognaM., RomualdiC.2010. Gene set analysis exploiting the topology of a pathway. BMC Systems Biol, 4:121.
39.
MathisD., VenceL., BenoistC.2001. Beta-cell death during progression to diabetes. Nature, 414:792–798.
40.
McCarthyM.I., AbecasisG.R., CardonL.R.et al.2008. Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet, 9:356–369.
41.
MishraG.R., SureshM., KumaranK.et al.2006. Human protein reference database—2006 update. Nucleic Acids Res, 34:D411–D414.
42.
MorahanG., MehtaM., JamesI.et al.2011. Tests for genetic interactions in type 1 diabetes: linkage and stratification analyses of 4,422 affected sib-pairs. Diabetes, 60:1030–1040.
43.
NitschD., TrancheventL.C., ThienpontB.et al.2009. Network analysis of differential expression for the identification of disease-causing genes. PloS One, 4:e5526.
44.
OtiM., BrunnerH.G.2007. The modular nature of genetic diseases. Clin Genet, 71:1–11.
45.
PflegerC., MortensenH.B., HansenL.et al.2008. Association of IL-1ra and adiponectin with C-peptide and remission in patients with type 1 diabetes. Diabetes, 57:929–937.
46.
PociotF., AkolkarB., ConcannonP.et al.2010. Genetics of type 1 diabetes: what's next?Diabetes, 59:1561–1571.
47.
RichS.S., AkolkarB., ConcannonP.et al.2009. Overview of the Type I Diabetes Genetics Consortium. Genes Immun, 10,Suppl 1:S1–S4.
48.
Ruefli-BrasseA.A., LeeW.P., HurstS., DixitV.M.2004. Rip2 participates in Bcl10 signaling and T-cell receptor-mediated NF-kappaB activation. J Biological Chem, 279:1570–1574.
49.
SharpC.D., HuangM., GlaweJ.et al.2008. Stromal cell-derived factor-1/CXCL12 stimulates chemorepulsion of NOD/LtJ T-cell adhesion to islet microvascular endothelium. Diabetes, 57:102–112.
50.
SmithN.G., Eyre-WalkerA.2003. Human disease genes: patterns and predictions. Gene, 318:169–175.
51.
SmootM., OnoK., IdekerT., MaereS.2011. PiNGO: a Cytoscape plugin to find candidate genes in biological networks. Bioinformatics (Oxf, England), 27:1030–1031.
52.
SuJ., YoonB.J., DoughertyE.R.2010. Identification of diagnostic subnetwork markers for cancer in human protein-protein interaction network. BMC Bioinformatics, 11,Suppl 6:S8.
53.
ThomasR., GohlkeJ., StopperG., ParhamF., PortierC.2009. Choosing the right path: enhancement of biologically relevant sets of genes or proteins using pathway structure. Genome Biol, 10:R44.
54.
TiffinN., AdieE., TurnerF.et al.2006. Computational disease gene identification: a concert of methods prioritizes type 2 diabetes and obesity candidate genes. Nucleic Acids Res, 34:3067–3081.
55.
VeyrierasJ.-B., KudaravalliS., KimS.Y.et al.2008. High-resolution mapping of expression-QTLs yields insight into human gene regulation. PLoS Genet, 4:e1000214.
56.
WangK., LiM., BucanM.2007. Pathway-based approaches for analysis of genomewide association studies. Am J Hum Genet, 81:1278–1283.
57.
WangX., JiaS., GeoffreyR., AlemzadehR., GhoshS., HessnerM.J.2008. Identification of a molecular signature in human type 1 diabetes mellitus using serum and functional genomics. J Immunol, 180:1929–1937.
58.
WangY., XiaY.2008. Condition specific subnetwork identification using an optimization model. The Second International Symposium on Optimization and Systems Biology.
59.
Wellcome Trust Case-Control Consortium. 2007. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature, 447:661–678.
60.
XuJ., LiY.2006. Discovering disease-genes by topological features in human protein-protein interaction network. Bioinformatics (Oxf, England), 22:2800–2805.
61.
ZeitlmannL., SirimP., KremmerE., KolanusW.2001. Cloning of ACP33 as a novel intracellular ligand of CD4. J Biol Chem, 276:9123–9132.
62.
ZhongH., BeaulaurierJ., LumP.Y.et al.2010. Liver and adipose expression associated SNPs are enriched for association to type 2 diabetes. PLoS Genet, 6:e1000932.
63.
ZhuM., ZhaoS.2007. Candidate gene identification approach: progress and challenges. Int J Biol Sci, 3:420–427.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
