Highly interconnected genes in disease-specific networks are enriched for disease-associated polymorphisms
- Fredrik Barrenäs†1,
- Sreenivas Chavali†2,
- Alexessander Couto Alves3,
- Lachlan Coin3,
- Marjo-Riitta Jarvelin3, 4,
- Rebecka Jörnsten5,
- Michael A Langston6,
- Adaikalavan Ramasamy3, 7,
- Gary Rogers6,
- Hui Wang1 and
- Mikael Benson1, 8Email author
© Barrenäs et al.; licensee BioMed Central Ltd. 2012
Received: 27 February 2012
Accepted: 15 June 2012
Published: 15 June 2012
Complex diseases are associated with altered interactions between thousands of genes. We developed a novel method to identify and prioritize disease genes, which was generally applicable to complex diseases.
We identified modules of highly interconnected genes in disease-specific networks derived from integrating gene-expression and protein interaction data. We examined if those modules were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies. First, we analyzed publicly available gene expression microarray and genome-wide association study (GWAS) data from 13, highly diverse, complex diseases. In each disease, highly interconnected genes formed modules, which were significantly enriched for genes harboring disease-associated SNPs. To test if such modules could be used to find novel genes for functional studies, we repeated the analyses using our own gene expression microarray and GWAS data from seasonal allergic rhinitis. We identified a novel gene, FGF2, whose relevance was supported by functional studies using combined small interfering RNA-mediated knock-down and gene expression microarrays. The modules in the 13 complex diseases analyzed here tended to overlap and were enriched for pathways related to oncological, metabolic and inflammatory diseases. This suggested that this union of the modules would be associated with a general increase in susceptibility for complex diseases. Indeed, we found that this union was enriched with GWAS genes for 145 other complex diseases.
Modules of highly interconnected complex disease genes were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies.
Medical research often focuses on individual diseases and genes. However, complex diseases show considerable comorbidity and are associated with altered interactions between thousands of genes. This suggests a need to find generally applicable principles to study multiple diseases and genes. One solution may be to map differentially expressed, disease-associated genes on to the human protein-protein interaction (PPI) network. Gene expression microarray studies of several complex diseases have shown that differentially expressed genes tend to form modules of interacting and functionally related genes [1–5]. Those modules may help to identify genes harboring disease-associated SNPs . The identification is, however, complicated by the involvement of multiple modules in the same complex disease.
In this study, we hypothesized that modules containing the most interconnected complex disease-associated genes would be enriched for disease-associated SNPs (note that highly interconnected disease genes have many interactions with other disease genes, while hub genes have interactions with any other gene). This hypothesis was based on recent discoveries in network medicine. Firstly, the effects of disease-associated SNPs tend to propagate through the PPI network, affecting the local neighborhood around the SNP-harboring genes [7–9]. Secondly, genes harboring disease-associated SNPs tend to form modules in the PPI network [10–13]. Those studies are mainly based on rare hereditary diseases. A recent meta-analysis of genome-wide association studies (GWAS), however, suggests that genes harboring disease-associated SNPs in complex diseases are also highly interconnected . Taken together, previous studies showed that genes harboring disease-associated SNPs tend to form modules, and that the same is true for differentially expressed genes . In this study, we integrated these two findings and used modules formed by differentially expressed disease genes to find genes harboring disease-associated SNPs. For this, first we defined modules in disease-specific networks for 13 complex diseases and show that these highly interconnected genes in these modules are enriched for disease genes identified through GWAS. To test the general applicability of our findings, the selected diseases were highly diverse and included oncological, metabolic and inflammatory diseases. Using in-house generated gene-expression and GWAS data, we showed that such modules could be used for identifying novel genes for functional studies, using seasonal allergic rhinitis (SAR) as a disease model. Finally, we show that overlapping modules of the complex diseases are generally enriched for genes harboring disease-associated SNPs, especially pleiotropic genes, identified by GWAS of 145 complex diseases.
Disease-specific core susceptibility modules are enriched for disease-associated genes in 13 complex diseases
Overview of SuMs in 13 complex diseases
GWAS genes in PPI
GWAS genes in SuM
GWAS genes in core SuM
Genes in SuM
Chronic lymphocytic leukemia
Type 2 diabetes
Core SuMs could be used to find novel genes for functional studies
We then determined whether core SuMs could be used to find novel genes for functional studies. We repeated the analyses using our own gene expression microarray and GWAS data from patients with SAR. This is an ideal model of complex diseases because it is possible to mimic and analyze the disease process in allergen-challenged cells from patients  (Extended background in Additional file 2). Novel genes can be functionally examined by combining small interfering RNA (siRNA) and gene expression microarrays in Th2 polarized cells [17–19], (Extended background in Additional file 2).
We constructed a SuM and a core SuM for SAR using differentially expressed genes obtained by performing gene expression microarray analysis of allergen-challenged CD4+ cells in samples from 12 patients as previously described [19, 20] (Additional file 3). The SuM included 622 of the 2,822 differentially expressed genes in the PPI network, and 1,191 genes in total (Figure 1a in Additional file 4). Next, we tested if they could be replicated in the additional study material. The repeated analyses resulted in highly similar SuMs and core SuMs (P < 10-15 in both cases, determined by a χ2-test; Figure 1b in Additional file 4). We also compared gene expression microarray data from allergen-challenged CD4+ cells from patients with SAR and allergen-challenged CD4+ cells from healthy controls and found differences in disease-relevant pathways and genes (Extended results in Additional file 2).
To determine whether the genes in the SuM and core SuM were enriched for SAR-associated SNPs, we analyzed an independent GWAS of 4,772 individuals in the North Finland Birth Cohort (Figure 1e). We found that intragenic SNPs within the 119 genes in the core SuM were 3.4 times more likely to be disease-associated than expected by chance (P = 1 × 10-5). This led to the identification of two novel genes, FGF2 and MAPK8 (Additional file 5). These findings were appropriately supported by false discovery rate calculations (Extended results in Additional file 2; Additional file 6). These genes had not been previously associated with SAR . While MAPK8 has a known role in type 1 allergic inflammation [22–24], FGF2 is a novel gene. We tested the functional relevance of FGF2 by siRNA-mediated knock-down of this gene in Th2 polarized cells, followed by gene expression microarrays (Extended methods and Extended results in Additional file 2; Additional file 7). The knock-down resulted in altered expression of several pathways of potential relevance for type 1 allergic inflammation, as well as individual genes of known relevance for type 1 allergic inflammation, including MAFB and NFKB1 [25–30] (Additional file 8). In contrast to the core SuM, no enrichment of SNPs was found in the SuM.
The union of core SuMs of different diseases was enriched for pleiotropic genes
We have developed a novel method to define SuMs and core SuMs for complex diseases by combining gene-expression microarray and PPI data. To show the general applicability of the method, we analyzed diseases with highly divergent phenotypes, rather than focusing on a specific subset of phenotypically related diseases. We found that SuMs, and in particular core SuMs, were enriched for GWAS genes. By comparison, no enrichment of GWAS genes was found when analyzing all differentially expressed genes.
To test if core SuMs could be used to find novel genes for functional studies, we analyzed our own gene expression microarray and GWAS data from patients with SAR. This is an optimal model of complex diseases in that it has a clearly defined phenotype that occurs at a given time point each year and the external cause (pollen) and key cell-type (lymphocytes) are known (Extended background in Additional file 2). We identified a SuM and a core SuM, which were reproduced in an independent study. We tested for enrichment of GWAS genes by analyzing a large population-based GWAS. This led to the identification of two novel genes in the core SuM. One of those genes, FGF2, was novel in type 1 allergic inflammation. The relevance of this gene was supported by transcriptomal analysis following siRNA-mediated knock-down. The knock-down of FGF2 resulted in altered expression of pathways and genes of potential or known relevance for type 1 allergic inflammation. In contrast, we found no enrichment of GWAS genes among all differentially expressed genes. These findings support the hypothesis that core SuMs can be used to identify GWAS genes with moderate effect sizes as well as novel genes for functional studies.
Studies of complex diseases often focus only on specific diseases and the genes associated with them. However, complex diseases show considerable phenotypic and genotypic overlap. Moreover, the effects of individual complex disease genes are generally small, while their collective contribution may be large . This has led to increasing interest in studying groups of diseases and genes [31, 32, 34, 35]. In the final part of our study, we made a corresponding change of scale and considered whether core SuMs from the 13 studied diseases were associated with an increased susceptibility for complex diseases, in general. We found that core SuMs tended to overlap, suggesting that they were involved in shared pathogenic mechanisms in complex diseases. This was supported by the union of the core SuMs being enriched for pathways involved in complex diseases. Finally, we tested if the union of the core SuMs was generally enriched for GWAS genes from 145 complex diseases, which represented hundreds of thousands of patients. Indeed, that union was highly enriched for such GWAS genes. Interestingly, the enrichment was greater when considering only GWAS genes associated with more than one disease. Taken together, our findings showed that the core SuMs of highly interconnected disease genes were associated with increased susceptibility for complex diseases.
We propose that the pathways in the core SuMs increase our understanding of how shared pathogenic mechanisms contribute to complex diseases, and also help explain why many of those diseases show phenotypic overlap . From a therapeutic perspective, the core SuMs may be used to prioritize therapeutic novel genes. It should, however, be noted that even though the number of genes in the core SuMs was considerably smaller than all the differentially expressed genes, it may be difficult to find individual target genes. Instead, perhaps drugs targeting combinations of core SuM genes will be required. Another interesting possibility is that the same novel genes might be exploited to develop drugs that target more than one disease. From a diagnostic perspective, different combinations of SuM and core SuM proteins may be useful as diagnostic markers. We propose that such studies can be performed in both specific disease groups and diseases.
A limitation of this study is that it is mainly based on GWAS genes reported in a public database. This may result in missing GWAS genes in some core SuMs, due to too stringent cutoffs. In this study, we analyzed only intragenic SNPs, while intergenic SNPs may have regulatory and disease-causing roles. Another limitation is that the study is based on known physical and functional PPIs. Thus, our results may be confounded by knowledge bias. Moreover, some of the gene expression microarray studies used in this effort were performed in tissues that contained mixed cell populations. We anticipate that the increasing accuracy of PPIs as well as availability of gene expression microarray and GWAS data will lead to more accurate identification of SuMs.
SuMs and core SuMs may be used to find novel genes for functional studies, as well as to increase understanding of the specific and shared pathogenic mechanisms in complex diseases and how they relate to phenotypic manifestations.
Materials and methods
Definition of susceptibility modules of complex diseases
SuMs were defined by integrating PPI network data and differentially expressed genes for each disease. The PPI network was assembled from a large set of functional and physical PPIs obtained from STRING (version 8, using interactions with a confidence score ≥ 0.7). The SuMs were identified using a step-wise process. First, maximal cliques were extracted from the PPI network (A clique is a complete sub-network, that is, a sub-network with links connecting every pair of its nodes (Figure 1a) . Such a clique is maximal if it is not properly contained within another clique.) For this task we employed our custom clique extraction tools as previously applied [38, 39]. We noted that maximal cliques tend to be highly overlapping and we used all cliques down to a minimum size of 2 that were not part of other cliques. Each clique was tested for enrichment (P < 0.05) of differentially expressed genes in the disease using Fisher's exact test (Figure 1b). Differentially expressed genes between the patient and control samples were determined using a Student's t-test with a P < 0.05. Finally, overlapping enriched modules/cliques were mapped onto the PPI network so that each gene was represented only once, and overlapping enriched cliques could be identified as SuMs (Figure 1c). The SuMs included both differentially expressed genes and their neighbors in these cliques. In the SuMs, we identified highly interconnected genes and defined them as core SuMs (Figure 1d). Interconnectivity was measured by calculating the average shortest path length using the Network-Analyzer (v2.6.1) plug-in in Cytoscape (v2.6.0).
Enrichment of core SuMs for GWAS genes
GWAS genes for all complex diseases were obtained from 'a catalog of published Genome-wide association studies' . We tested for enrichment of GWAS genes using a permutation test. In each permutation, every GWAS gene was replaced with a random gene. The number of random genes present in the corresponding SuMs or core SuMs was noted. This process was repeated 100,000 times. The total number of random genes present in the disease SuMs represented the probability of finding GWAS genes in the SuMs by chance. The significance of the overrepresentation of GWAS-identified genes in the SuMs was measured as the ratio of random permutations that included as many or more SuM genes than the true GWAS-identified genes for all 13 diseases. Similarly, enrichment of GWAS genes in the union of core SuMs was also determined for 145 diseases. Pathway enrichment for the union of core SuM genes was determined using Ingenuity Pathway Analysis .
Gene-expression microarray analysis in seasonal allergic rhinitis
Peripheral blood mononuclear cells from patients with SAR and healthy controls were prepared and stimulated with grass pollen extract or diluent for seven days . For gene expression studies, T helper cells were enriched from the allergen-challenged peripheral blood mononuclear cells using anti-CD4-coated paramagnetic microbeads and a MACS (magnetic cell sorter) system according to the instructions of the manufacturer (Miltenyi Biotec GmbH, Bergisch Gladbach, Germany). cRNA was extracted from 200 ng total RNA using Ambion's Illumina RNA TotalPrep Amplification kit (Ambion, Inc., USA). In vitro transcription reaction and cRNA biotinylation were performed overnight (14 h). The RNA/cRNA concentrations where checked using Nanodrop ND-1000 before and after the amplifications. cRNA quality was controlled by BioRad's Experion electrophoresis station (Bio-Rad Laboratories, Inc., CA, USA). Transcriptional profiling in 12 patients was performed using Illumina's Sentrix® Human-6 Expression BeadChips (Illumina Inc., San Diego, CA, USA) according to the manufacturer's instructions. (The data can be obtained from the Gene Expression Omnibus under accession number GSE18574.) Probes with a detection score below 0.95 were discarded prior to differential expression analysis. Differentially expressed genes were determined using lmFit from the Bioconductor package Limma  (Additional file 3). Genes with a P-value < 0.05 after correction for multiple comparisons (false discovery rate) were determined to be differentially expressed. To validate results obtained from the data above, an additional set of patients (n = 3) was analyzed on an Affymetrix U133A platform as previously described .
GWAS analysis for seasonal allergic rhinitis
This cohort (NFBC1966) included individuals from the provinces of Oulu and Lapland . In 1997 (when participants were aged 31 years), 8,463 survivors were sent postal questionnaires and invited to clinical examination with a 71% response rate. DNA was collected for the majority of participants and a total of 4,772 individuals were successfully genotyped. All aspects of the study were reviewed and approved by the Ethics Committee of the University of Oulu and participants gave written informed consent. Genotyping was done using the Illumina HumanCNV370-Duo chip. The data were imputed to approximately 2.5 million SNPs using NCBI HapMap II CEU build 35 version 21 after pre-filtering SNPs (genotyping rate > 95%, P-value for HWE deviation > 10-4, minor allele frequency > 1%, imputation with a confidence call of R2 > 0.5) using IMPUTE .
The SAR phenotype was defined as a positive response to the 'have you ever had allergic rhinitis' question in the main questionnaire and a positive skin prick test for grass. Controls had a negative response to the question and a negative skin prick test. A positive skin prick test was defined as a mean wheal reaction to grass extract of at least 3 mm. Participants with a positive reaction to negative control (diluent of allergen extracts) or a negative reaction to positive control (10 mg/ml histamine dihydrochloride) were excluded. There were 456 patients that had SAR, while 2,569 individuals were controls.
where p(x, k, m, n)x is the P-value of the hypergeometric distribution when x out of the k significant SNPs are observed within the core SuM, when the core SuM holds m of the total n SNPs. The number of true positive SNP associations (TP) is given by TP = P - FP. The odds ratio of the enrichment was defined as the odds of a significant SNP within the core SuM (or any given subset of genes) as compared to the remaining genes. The linkage disequilibrium within the core SuM is assumed to be similar to the expected linkage disequilibrium in the rest of the genome; therefore, GWAS P-value distributions within and outside the core SuM should be comparable.
genome-wide association studies
seasonal allergic rhinitis
small interfering RNA
We thank Petter Holme for his insightful comments on this manuscript. This work was supported by the European Community's Seventh Framework Programme ([FP7/2007-2012] under grant agreement 223367 (MultiMod)), by the US National Institutes of Health under grants R01-AA-018776 and 3P20MD000516-07S1, and by the Swedish Research Council. SC is an EMBO Long Term Fellow.
- Ge H, Liu Z, Church GM, Vidal M: Correlation between transcriptome and interactome mapping data from Saccharomyces cerevisiae. Nat Genet. 2001, 29: 482-486. 10.1038/ng776.PubMedView ArticleGoogle Scholar
- Minguez P, Dopazo J: Assessing the biological significance of gene expression signatures and co-expression modules by studying their network properties. PLoS One. 2011, 6: e17474-10.1371/journal.pone.0017474.PubMedPubMed CentralView ArticleGoogle Scholar
- Narayanan M, Vetta A, Schadt EE, Zhu J: Simultaneous clustering of multiple gene expression and physical interaction datasets. PLoS Comput Biol. 2010, 6: e1000742-10.1371/journal.pcbi.1000742.PubMedPubMed CentralView ArticleGoogle Scholar
- Ramani AK, Li Z, Hart GT, Carlson MW, Boutz DR, Marcotte EM: A map of human protein interactions derived from co-expression of human mRNAs and their orthologs. Mol Syst Biol. 2008, 4: 180-PubMedPubMed CentralView ArticleGoogle Scholar
- Chen J, Liang H, Fernandez A: Protein structure protection commits gene expression patterns. Genome Biol. 2008, 9: R107-10.1186/gb-2008-9-7-r107.PubMedPubMed CentralView ArticleGoogle Scholar
- Pujana MA, Han JD, Starita LM, Stevens KN, Tewari M, Ahn JS, Rennert G, Moreno V, Kirchhoff T, Gold B, Assmann V, Elshamy WM, Rual JF, Levine D, Rozek LS, Gelman RS, Gunsalus KC, Greenberg RA, Sobhian B, Bertin N, Venkatesan K, Ayivi-Guedehoussou N, Sole X, Hernandez P, Lazaro C, Nathanson KL, Weber BL, Cusick ME, Hill DE, Offit K, et al: Network modeling links breast cancer susceptibility and centrosome dysfunction. Nat Genet. 2007, 39: 1338-1349. 10.1038/ng.2007.2.PubMedView ArticleGoogle Scholar
- Barabasi AL, Gulbahce N, Loscalzo J: Network medicine: a network-based approach to human disease. Nat Rev Genet. 2011, 12: 56-68. 10.1038/nrg2918.PubMedPubMed CentralView ArticleGoogle Scholar
- Goldstein DB: Common genetic variation and human traits. N Engl J Med. 2009, 360: 1696-1698. 10.1056/NEJMp0806284.PubMedView ArticleGoogle Scholar
- Schadt EE: Molecular networks as sensors and drivers of common human diseases. Nature. 2009, 461: 218-223. 10.1038/nature08454.PubMedView ArticleGoogle Scholar
- Lage K, Karlberg EO, Storling ZM, Olason PI, Pedersen AG, Rigina O, Hinsby AM, Tumer Z, Pociot F, Tommerup N, Moreau Y, Brunak S: A human phenome-interactome network of protein complexes implicated in genetic disorders. Nat Biotechnol. 2007, 25: 309-316. 10.1038/nbt1295.PubMedView ArticleGoogle Scholar
- Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL: The human disease network. Proc Natl Acad Sci USA. 2007, 104: 8685-8690. 10.1073/pnas.0701361104.PubMedPubMed CentralView ArticleGoogle Scholar
- Oti M, Snel B, Huynen MA, Brunner HG: Predicting disease genes using protein-protein interactions. J Med Genet. 2006, 43: 691-698. 10.1136/jmg.2006.041376.PubMedPubMed CentralView ArticleGoogle Scholar
- Feldman I, Rzhetsky A, Vitkup D: Network properties of genes harboring inherited disease mutations. Proc Natl Acad Sci USA. 2008, 105: 4323-4328. 10.1073/pnas.0701722105.PubMedPubMed CentralView ArticleGoogle Scholar
- Barrenas F, Chavali S, Holme P, Mobini R, Benson M: Network properties of complex human disease genes identified through genome-wide association studies. PLoS One. 2009, 4: e8090-10.1371/journal.pone.0008090.PubMedPubMed CentralView ArticleGoogle Scholar
- Chuang HY, Lee E, Liu YT, Lee D, Ideker T: Network-based classification of breast cancer metastasis. Mol Syst Biol. 2007, 3: 140-PubMedPubMed CentralView ArticleGoogle Scholar
- Woodfolk JA: T-cell responses to allergens. J Allergy Clin Immunol. 2007, 119: 280-294. 10.1016/j.jaci.2006.11.008. quiz 295-286PubMedView ArticleGoogle Scholar
- Bruhn S, Barrenas F, Mobini R, Andersson BA, Chavali S, Egan BS, Hovig E, Sandve GK, Langston MA, Rogers G, Wang H, Benson M: Increased expression of IRF4 and ETS1 in CD4(+) cells from patients with intermittent allergic rhinitis. Allergy. 2012, 67: 33-40. 10.1111/j.1398-9995.2011.02707.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Banfield G, Watanabe H, Scadding G, Jacobson MR, Till SJ, Hall DA, Robinson DS, Lloyd CM, Nouri-Aria KT, Durham SR: CC chemokine receptor 4 (CCR4) in human allergen-induced late nasal responses. Allergy. 2010, 65: 1126-1133.PubMedPubMed CentralGoogle Scholar
- Benson M, Langston MA, Adner M, Andersson B, Torinssson-Naluai A, Cardell LO: A network-based analysis of the late-phase reaction of the skin. J Allergy Clin Immunol. 2006, 118: 220-225. 10.1016/j.jaci.2006.03.006.PubMedView ArticleGoogle Scholar
- Benson M, Carlsson L, Guillot G, Jernas M, Langston MA, Rudemo M, Andersson B: A network-based analysis of allergen-challenged CD4+ T cells from patients with allergic rhinitis. Genes Immun. 2006, 7: 514-521. 10.1038/sj.gene.6364322.PubMedView ArticleGoogle Scholar
- Ramasamy A, Curjuric I, Coin LJ, Kumar A, McArdle WL, Imboden M, Leynaert B, Kogevinas M, Schmid-Grendelmeier P, Pekkanen J, Wjst M, Bircher AJ, Sovio U, Rochat T, Hartikainen AL, Balding DJ, Jarvelin MR, Probst-Hensch N, Strachan DP, Jarvis DL: A genome-wide meta-analysis of genetic variants associated with allergic rhinitis and grass sensitization and their interaction with birth order. J Allergy Clin Immunol. 2011, 128: 996-1005. 10.1016/j.jaci.2011.08.030.PubMedView ArticleGoogle Scholar
- Dong C, Yang DD, Tournier C, Whitmarsh AJ, Xu J, Davis RJ, Flavell RA: JNK is required for effector T-cell function but not for T-cell activation. Nature. 2000, 405: 91-94. 10.1038/35011091.PubMedView ArticleGoogle Scholar
- Dong C, Yang DD, Wysk M, Whitmarsh AJ, Davis RJ, Flavell RA: Defective T cell differentiation in the absence of Jnk1. Science. 1998, 282: 2092-2095.PubMedView ArticleGoogle Scholar
- Chu X, Ci X, He J, Wei M, Yang X, Cao Q, Li H, Guan S, Deng Y, Pang D, Deng X: A novel anti-inflammatory role for ginkgolide B in asthma via inhibition of the ERK/MAPK signaling pathway. Molecules. 2011, 16: 7634-7648. 10.3390/molecules16097634.PubMedView ArticleGoogle Scholar
- Hashizume H, Hamalainen H, Sun Q, Sucharczuk A, Lahesmaa R: Downregulation of mafB expression in T-helper cells during early differentiation in vitro. Scand J Immunol. 2003, 57: 28-34. 10.1046/j.1365-3083.2003.01181.x.PubMedView ArticleGoogle Scholar
- Chen Q, Dowhan DH, Liang D, Moore DD, Overbeek PA: CREB-binding protein/p300 co-activation of crystallin gene expression. J Biol Chem. 2002, 277: 24081-24089. 10.1074/jbc.M201821200.PubMedView ArticleGoogle Scholar
- Cornejo-Garcia JA, Fernandez TD, Torres MJ, Carballo M, Hernan I, Antunez C, Blanca M, Mayorga C: Differential cytokine and transcription factor expression in patients with allergic reactions to drugs. Allergy. 2007, 62: 1429-1438. 10.1111/j.1398-9995.2007.01542.x.PubMedView ArticleGoogle Scholar
- Erpenbeck VJ, Hagenberg A, Krentel H, Discher M, Braun A, Hohlfeld JM, Krug N: Regulation of GATA-3, c-maf and T-bet mRNA expression in bronchoalveolar lavage cells and bronchial biopsies after segmental allergen challenge. Int Arch Allergy Immunol. 2006, 139: 306-316. 10.1159/000091602.PubMedView ArticleGoogle Scholar
- Tenda Y, Yamashita M, Kimura MY, Hasegawa A, Shimizu C, Kitajima M, Onodera A, Suzuki A, Seki N, Nakayama T: Hyperresponsive TH2 cells with enhanced nuclear factor-kappa B activation induce atopic dermatitis-like skin lesions in Nishiki-nezumi Cinnamon/Nagoya mice. J Allergy Clin Immunol. 2006, 118: 725-733. 10.1016/j.jaci.2006.05.024.PubMedView ArticleGoogle Scholar
- El-Hashim AZ, Renno WM, Abduo HT, Jaffal SM, Akhtar S, Benter IF: Effect of inhibition of the ubiquitin-proteasome-system and IkappaB kinase on airway inflammation and hyperresponsiveness in a murine model of asthma. Int J Immunopathol Pharmacol. 2011, 24: 33-42.PubMedGoogle Scholar
- Ramos PS, Criswell LA, Moser KL, Comeau ME, Williams AH, Pajewski NM, Chung SA, Graham RR, Zidovetzki R, Kelly JA, Kaufman KM, Jacob CO, Vyse TJ, Tsao BP, Kimberly RP, Gaffney PM, Alarcon-Riquelme ME, Harley JB, Langefeld CD: A Comprehensive analysis of shared loci between systemic lupus erythematosus (SLE) and sixteen autoimmune diseases reveals limited genetic overlap. PLoS Genet. 2011, 7: e1002406-10.1371/journal.pgen.1002406.PubMedPubMed CentralView ArticleGoogle Scholar
- Sivakumaran S, Agakov F, Theodoratou E, Prendergast JG, Zgaga L, Manolio T, Rudan I, McKeigue P, Wilson JF, Campbell H: Abundant pleiotropy in human complex diseases and traits. Am J Hum Genet. 2011, 89: 607-618. 10.1016/j.ajhg.2011.10.004.PubMedPubMed CentralView ArticleGoogle Scholar
- Petretto E, Liu ET, Aitman TJ: A gene harvest revealing the archeology and complexity of human disease. Nat Genet. 2007, 39: 1299-1301. 10.1038/ng1107-1299.PubMedView ArticleGoogle Scholar
- Sirota M, Schaub MA, Batzoglou S, Robinson WH, Butte AJ: Autoimmune disease classification by inverse association with SNP alleles. PLoS Genet. 2009, 5: e1000792-10.1371/journal.pgen.1000792.PubMedPubMed CentralView ArticleGoogle Scholar
- Torkamani A, Dean B, Schork NJ, Thomas EA: Coexpression network analysis of neural tissue reveals perturbations in developmental processes in schizophrenia. Genome Res. 2010, 20: 403-412. 10.1101/gr.101956.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Chavali S, Barrenas F, Kanduri K, Benson M: Network properties of human disease genes with pleiotropic effects. BMC Syst Biol. 2010, 4: 78-10.1186/1752-0509-4-78.PubMedPubMed CentralView ArticleGoogle Scholar
- Bomze J, Budinich M, Pardalos P, Pelillo M: The maximum clique problem. Handbook of Combinatorial Optimization. Edited by: Du D-Z, Pardalos PM. 1999, Kluwer Academic Publishers, 4:Google Scholar
- Chesler EJ, Lu L, Shou S, Qu Y, Gu J, Wang J, Hsu HC, Mountz JD, Baldwin NE, Langston MA, Threadgill DW, Manly KF, Williams RW: Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function. Nat Genet. 2005, 37: 233-242. 10.1038/ng1518.PubMedView ArticleGoogle Scholar
- Voy BH, Scharff JA, Perkins AD, Saxton AM, Borate B, Chesler EJ, Branstetter LK, Langston MA: Extracting gene networks for low-dose radiation using graph theoretical algorithms. PLoS Comput Biol. 2006, 2: e89-10.1371/journal.pcbi.0020089.PubMedPubMed CentralView ArticleGoogle Scholar
- Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA: Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci USA. 2009, 106: 9362-9367. 10.1073/pnas.0903103106.PubMedPubMed CentralView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMedGoogle Scholar
- Rantakallio P, Leino S: Social and biological factors adding to the risks of infants. Katilolehti. 1969, 74: 520-522.PubMedGoogle Scholar
- Marchini J, Howie B, Myers S, McVean G, Donnelly P: A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007, 39: 906-913. 10.1038/ng2088.PubMedView ArticleGoogle Scholar
- Hubbard T, Andrews D, Caccamo M, Cameron G, Chen Y, Clamp M, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T, Down T, Durbin R, Fernandez-Suarez XM, Gilbert J, Hammond M, Herrero J, Hotz H, Howe K, Iyer V, Jekosch K, Kahari A, Kasprzyk A, Keefe D, Keenan S, Kokocinsci F, London D, Longden I, McVicker G, et al: Ensembl 2005. Nucleic Acids Res. 2005, 33: D447-453. 10.1093/nar/gki378.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.