Integrating Genome-Wide Genetic Variations and Monocyte Expression Data Reveals -Regulated Gene Modules in Humans

English version České info

One major expectation from the transcriptome in humans is to characterize the biological basis of associations identified by genome-wide association studies. So far, few cis expression quantitative trait loci (eQTLs) have been reliably related to disease susceptibility. Trans-regulating mechanisms may play a more prominent role in disease susceptibility. We analyzed 12,808 genes detected in at least 5% of circulating monocyte samples from a population-based sample of 1,490 European unrelated subjects. We applied a method of extraction of expression patterns—independent component analysis—to identify sets of co-regulated genes. These patterns were then related to 675,350 SNPs to identify major trans-acting regulators. We detected three genomic regions significantly associated with co-regulated gene modules. Association of these loci with multiple expression traits was replicated in Cardiogenics, an independent study in which expression profiles of monocytes were available in 758 subjects. The locus 12q13 (lead SNP rs11171739), previously identified as a type 1 diabetes locus, was associated with a pattern including two cis eQTLs, RPS26 and SUOX, and 5 trans eQTLs, one of which (MADCAM1) is a potential candidate for mediating T1D susceptibility. The locus 12q24 (lead SNP rs653178), which has demonstrated extensive disease pleiotropy, including type 1 diabetes, hypertension, and celiac disease, was associated to a pattern strongly correlating to blood pressure level. The strongest trans eQTL in this pattern was CRIP1, a known marker of cellular proliferation in cancer. The locus 12q15 (lead SNP rs11177644) was associated with a pattern driven by two cis eQTLs, LYZ and YEATS4, and including 34 trans eQTLs, several of them tumor-related genes. This study shows that a method exploiting the structure of co-expressions among genes can help identify genomic regions involved in trans regulation of sets of genes and can provide clues for understanding the mechanisms linking genome-wide association loci to disease.

Vyšlo v časopise: Integrating Genome-Wide Genetic Variations and Monocyte Expression Data Reveals -Regulated Gene Modules in Humans. PLoS Genet 7(12): e32767. doi:10.1371/journal.pgen.1002367
Kategorie: Research Article
prolekare.web.journal.doi_sk: https://doi.org/10.1371/journal.pgen.1002367

Souhrn

Zdroje

1. ManolioTACollinsFSCoxNJGoldsteinDBHindorffLA 2009 Finding the missing heritability of complex diseases. Nature 461 747 753

2. CooksonWLiangLAbecasisGMoffattMLathropM 2009 Mapping complex disease traits with global gene expression. Nat Rev Genet 10 184 194

3. GoringHHCurranJEJohnsonMPDyerTDCharlesworthJ 2007 Discovery of expression QTLs using large-scale transcriptional profiling in human lymphocytes. Nat Genet 39 1208 1216

4. DixonALLiangLMoffattMFChenWHeathS 2007 A genome-wide association study of global gene expression. Nat Genet 39 1202 1207

5. StrangerBENicaACForrestMSDimasABirdCP 2007 Population genomics of human gene expression. Nat Genet 39 1217 1224

6. EmilssonVThorleifssonGZhangBLeonardsonASZinkF 2008 Genetics of gene expression and its effect on disease. Nature 452 423 428

7. SchadtEEMolonyCChudinEHaoKYangX 2008 Mapping the genetic architecture of gene expression in human liver. PLoS Biol 6 e107 doi:10.1371/journal.pbio.0060107

8. IdaghdourYCzikaWShiannaKVLeeSHVisscherPM 2010 Geographical genomics of human leukocyte gene expression variation in southern Morocco. Nat Genet 42 62 67

9. KangHMYeCEskinE 2008 Accurate discovery of expression quantitative trait loci under confounding from spurious and genuine regulatory hotspots. Genetics 180 1909 1925

10. SchadtEEZhangBZhuJ 2009 Advances in systems biology are enhancing our understanding of disease and moving us closer to novel disease treatments. Genetica 136 259 269

11. YvertGBremRBWhittleJAkeyJMFossE 2003 Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors. Nat Genet 35 57 64

12. ZhuJZhangBSmithENDreesBBremRB 2008 Integrating large-scale functional genomic data to dissect the complexity of yeast regulatory networks. Nat Genet 40 854 861

13. AyrolesJFCarboneMAStoneEAJordanKWLymanRF 2009 Systems genetics of complex traits in Drosophila melanogaster. Nat Genet 41 299 307

14. MehrabianMAllayeeHStocktonJLumPYDrakeTA 2005 Integrating genotypic and expression data in a segregating mouse population to identify 5-lipoxygenase as a susceptibility gene for obesity and bone traits. Nat Genet 37 1224 1233

15. SchadtEELambJYangXZhuJEdwardsS 2005 An integrative genomics approach to infer causal associations between gene expression and disease. Nat Genet 37 710 717

16. GhazalpourADossSZhangBWangSPlaisierC 2006 Integrating genetic and network analysis to characterize genes related to mouse weight. PLoS Genet 2 e130 doi:10.1371/journal.pgen.0020130

17. BiswasSStoreyJDAkeyJM 2008 Mapping gene expression quantitative trait loci by singular value decomposition and independent component analysis. BMC Bioinformatics 9 244

18. HyvarinenAOjaE 2000 Independent component analysis: algorithms and applications. Neural Netw 13 411 430

19. LiebermeisterW 2002 Linear modes of gene expression determined by independent component analysis. Bioinformatics 18 51 60

20. ZellerTWildPSzymczakSRotivalMSchillertA 2010 Genetics and beyond—The transcriptome of human monocytes and disease susceptibility. PLoS ONE 5 e10693 doi:10.1371/journal.pone.0010693

21. LeeSIBatzoglouS 2003 Application of independent component analysis to microarrays. Genome Biol 4 R76

22. CarpentierASRivaATisseurPDidierGHenautA 2004 The operons, a criterion to compare the reliability of transcriptome analysis tools: ICA is more reliable than ANOVA, PLS and PCA. Comput Biol Chem 28 3 10

23. StrimmerK 2008 A unified approach to false discovery rate estimation. BMC Bioinformatics 9 303

24. LyonsPAKoukoulakiMHattonADoggettKWoffendinHB 2007 Microarray analysis of human leucocyte subsets: the advantages of positive selection and rapid purification. BMC Genomics 8 64

25. WatkinsNAGusnantoAde BonoBDeSMiranda-SaavedraD 2009 A HaemAtlas: characterizing gene expression in differentiated human blood cells. Blood 113 e1 9

26. SoranzoNSpectorTDManginoMKuhnelBRendonA 2009 A genome-wide meta-analysis identifies 22 loci associated with eight hematological parameters in the HaemGen consortium. Nat Genet 41 1182 1190

27. SatoSSanjoHTakedaKNinomiya-TsujiJYamamotoM 2005 Essential function for the kinase TAK1 in innate and adaptive immune responses. Nat Immunol 6 1087 1095

28. HeinigMPetrettoEWallaceCBottoloLRotivalM 2010 A conserved trans-acting regulatory locus underlies a proinflammatory gene expression network and susceptibility to autoimmune type 1 diabetes. Nature 467 460 464

29. The 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

30. ToddJAWalkerNMCooperJDSmythDJDownesK 2007 Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes. Nat Genet 39 857 864

31. PlagnolVSmythDJToddJAClaytonDG 2009 Statistical independence of the colocalized association signals for type 1 diabetes and RPS26 gene expression on chromosome 12q13. Biostatistics 10 327 334

32. HanninenATaylorCStreeterPRStarkLSSarteJM 1993 Vascular addressins are induced on islet vessels during insulitis in nonobese diabetic mice and are involved in lymphoid cell binding to islet endothelium. J Clin Invest 92 2509 2515

33. HanninenAJaakkolaIJalkanenS 1998 Mucosal addressin is required for the development of diabetes in nonobese diabetic mice. J Immunol 160 6018 6025

34. PhillipsJMHaskinsKCookeA 2005 MAdCAM-1 is needed for diabetes development mediated by the T cell clone, BDC-2.5. Immunology 116 525 531

35. XiongYChenXChenZWangXShiS RNA sequencing shows no dosage compensation of the active X-chromosome. Nat Genet 42 1043 1047

36. HuntKAZhernakovaATurnerGHeapGAFrankeL 2008 Newly identified genetic risk variants for celiac disease related to the immune response. Nat Genet 40 395 402

37. GudbjartssonDFBjornsdottirUSHalapiEHelgadottirASulemP 2009 Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction. Nat Genet 41 342 347

38. Newton-ChehCJohnsonTGatevaVTobinMDBochudM 2009 Genome-wide association study identifies eight loci associated with blood pressure. Nat Genet

39. GaneshSKZakaiNAvan RooijFJSoranzoNSmithAV 2009 Multiple loci influence erythrocyte phenotypes in the CHARGE Consortium. Nat Genet 41 1191 1198

40. LevyDEhretGBRiceKVerwoertGCLaunerLJ 2009 Genome-wide association study of blood pressure and hypertension. Nat Genet

41. VelazquezLChengAMFlemingHEFurlongerCVeselyS 2002 Cytokine signaling and hematopoietic homeostasis are disrupted in Lnk-deficient mice. J Exp Med 195 1599 1611

42. KadrmasJLBeckerleMC 2004 The LIM domain: from the cytoskeleton to the nucleus. Nat Rev Mol Cell Biol 5 920 931

43. Lanningham-FosterLGreenCLLangkamp-HenkenBDavisBANguyenKT 2002 Overexpression of CRIP in transgenic mice alters cytokine patterns and the immune response. Am J Physiol Endocrinol Metab 282 E1197 1203

44. HaoJSerohijosAWNewtonGTassoneGWangZ 2008 Identification and rational redesign of peptide ligands to CRIP1, a novel biomarker for cancers. PLoS Comput Biol 4 e1000138 doi:10.1371/journal.pcbi.1000138

45. SchulzeJMWangAYKoborMS 2009 YEATS domain proteins: a diverse family with many links to chromatin modification and transcription. Biochem Cell Biol 87 65 75

46. KellerMPChoiYWangPDavisDBRabagliaME 2008 A gene expression network model of type 2 diabetes links cell cycle regulation in islets with diabetes susceptibility. Genome Res 18 706 716

47. ChenYZhuJLumPYYangXPintoS 2008 Variations in DNA elucidate molecular networks that cause disease. Nature 452 429 435

48. PlaisierCLHorvathSHuertas-VazquezACruz-BautistaIHerreraMF 2009 A systems genetics approach implicates USF1, FADS3, and other causal candidate genes for familial combined hyperlipidemia. PLoS Genet 5 e1000642 doi:10.1371/journal.pgen.1000642

49. BrunetJPTamayoPGolubTRMesirovJP 2004 Metagenes and molecular pattern discovery using matrix factorization. Proc Natl Acad Sci U S A 101 4164 4169

50. GrieveICDickensNJPravenecMKrenVHubnerN 2008 Genome-wide co-expression analysis in multiple tissues. PLoS ONE 3 e4033 doi:10.1371/journal.pone.0004033

51. FriedmanNLinialMNachmanIPe'erD 2000 Using Bayesian networks to analyze expression data. J Comput Biol 7 601 620

52. ZhangXWYapYLWeiDChenFDanchinA 2005 Molecular diagnosis of human cancer type by gene expression profiles and independent component analysis. Eur J Hum Genet 13 1303 1311

53. TeschendorffAEJourneeMAbsilPASepulchreRCaldasC 2007 Elucidating the altered transcriptional programs in breast cancer using independent component analysis. PLoS Comput Biol 3 e161 doi:10.1371/journal.pcbi.0030161

54. LutterDUgocsaiPGrandlMOrsoETheisF 2008 Analyzing M-CSF dependent monocyte/macrophage differentiation: expression modes and meta-modes derived from an independent component analysis. BMC Bioinformatics 9 100

55. Shen-OrrSSTibshiraniRKhatriPBodianDLStaedtlerF 2010 Cell type-specific gene expression differences in complex tissues. Nat Methods 7 287 289

56. AbbasARWolslegelKSeshasayeeDModrusanZClarkHF 2009 Deconvolution of blood microarray data identifies cellular activation patterns in systemic lupus erythematosus. PLoS ONE 4 e6098 doi:10.1371/journal.pone.0006098

57. AulchenkoYSRipkeSIsaacsAvan DuijnCM 2007 GenABEL: an R library for genome-wide association analysis. Bioinformatics 23 1294 1296