Dissection of a Complex Disease Susceptibility Region Using a Bayesian Stochastic Search Approach to Fine Mapping

PLoS Genetics

Volumn 11, Issue 6, 2015, Pages

  Wallace, Chris ^a,b   Cutler, Antony J ^a   Pontikos, Nikolas ^a   Pekalski, Marcin L ^a   Burren, Oliver S ^a   Cooper, Jason D ^a,g   García, Arcadio Rubio ^a   Ferreira, Ricardo C ^a   Guo, Hui ^a,c   Walker, Neil M ^a   Smyth, Deborah J ^a   Rich, Stephen S ^d   Onengut Gumuscu, Suna ^d,e   Sawcer, Stephen J ^f   Ban, Maria ^f   Richardson, Sylvia ^b   Todd, John A ^a   Wicker, Linda S ^a  

^a UNIVERSITY OF CAMBRIDGE   (United Kingdom)

^b UNIVERSITY OF CAMBRIDGE   (United Kingdom)

^c UNIVERSITY OF MANCHESTER   (United Kingdom)

^d UNIVERSITY OF VIRGINIA   (United States)

^e UNIVERSITY OF VIRGINIA   (United States)

^f UNIVERSITY OF CAMBRIDGE   (United Kingdom)

^g UNIVERSITY OF CAMBRIDGE   (United Kingdom)

Author keywords

[No Author keywords available]

Indexed keywords

INTERLEUKIN 2 RECEPTOR ALPHA;

ARTICLE; BAYES THEOREM; BAYESIAN EVOLUTIONARY STOCHASTIC SEARCH ALGORITHM; CD4+ T LYMPHOCYTE; CONTROLLED STUDY; EVOLUTIONARY ALGORITHM; GENE EXPRESSION; GENE LINKAGE DISEQUILIBRIUM; HUMAN; HUMAN CELL; IL2RA GENE; INSULIN DEPENDENT DIABETES MELLITUS; MULTIPLE SCLEROSIS; PROMOTER REGION; PROTEIN EXPRESSION; SINGLE NUCLEOTIDE POLYMORPHISM; T LYMPHOCYTE SUBPOPULATION; ALGORITHM; CHROMOSOMAL MAPPING; GENETIC PREDISPOSITION; GENETICS; HAPLOTYPE; MARKOV CHAIN; PROCEDURES; STATISTICS AND NUMERICAL DATA;

ALGORITHMS; BAYES THEOREM; CHROMOSOME MAPPING; DIABETES MELLITUS, TYPE 1; GENETIC PREDISPOSITION TO DISEASE; HAPLOTYPES; HUMANS; INTERLEUKIN-2 RECEPTOR ALPHA SUBUNIT; LINKAGE DISEQUILIBRIUM; MULTIPLE SCLEROSIS; POLYMORPHISM, SINGLE NUCLEOTIDE; STOCHASTIC PROCESSES;

EID: 84937791662     PISSN: 15537390     EISSN: 15537404     Source Type: Journal    
DOI: 10.1371/journal.pgen.1005272     Document Type: Article

References (61)

1. Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, et al. (2012) The accessible chromatin landscape of the human genome. Nature 489: 75–82. doi: 10.1038/nature11232 22955617
2. McCarthy MI, Hirschhorn JN, (2008) Genome-wide association studies: potential next steps on a genetic journey. Hum Mol Genet 17: R156–R165. doi: 10.1093/hmg/ddn289 18852205
3. Miller AJ, (1984) Selection of subsets of regression variables. Journal of the Royal Statistical Society Series A (General) 147: pp. 389–425. doi: 10.2307/2981576
4. Maller JB, McVean G, Byrnes J, Vukcevic D, Wellcome Trust Case Control Consortium, et al. (2012) Bayesian refinement of association signals for 14 loci in 3 common diseases. Nat Genet 44: 1294–1301. doi: 10.1038/ng.2435 23104008
5. Pickrell JK, (2014) Joint analysis of functional genomic data and genome-wide association studies of 18 human traits. Am J Hum Genet 94: 559–573. doi: 10.1016/j.ajhg.2014.03.004 24702953
6. Giambartolomei C, Vukcevic D, Schadt EE, Franke L, Hingorani AD, et al. (2014) Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet 10: e1004383. doi: 10.1371/journal.pgen.1004383 24830394
7. International Multiple Sclerosis Genetics Consortium (IMSGC) (2013) Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis. Nat Genet 45: 1353–1360. doi: 10.1038/ng.2770 24076602
8. Trynka G, Hunt KA, Bockett NA, Romanos J, Mistry V, et al. (2011) Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease. Nat Genet 43: 1193–1201. doi: 10.1038/ng.998 22057235
9. Eyre S, Bowes J, Diogo D, Lee A, Barton A, et al. (2012) High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet 44: 1336–1340. doi: 10.1038/ng.2462 23143596
10. Hinks A, Cobb J, Marion MC, Prahalad S, Sudman M, et al. (2013) Dense genotyping of immune-related disease regions identifies 14 new susceptibility loci for juvenile idiopathic arthritis. Nat Genet 45: 664–669. doi: 10.1038/ng.2614 23603761
11. Servin B, Stephens M, (2007) Imputation-based analysis of association studies: candidate regions and quantitative traits. PLoS Genet 3: e114. doi: 10.1371/journal.pgen.0030114 17676998
12. Bottolo L, Richardson S, (2010) Evolutionary stochastic search for bayesian model exploration. Bayesian Analysis 5: 583–618. doi: 10.1214/10-BA523
13. Bottolo L, Chadeau-Hyam M, Hastie DI, Zeller T, Liquet B, et al. (2013) GUESS-ing polygenic associations with multiple phenotypes using a GPU-based evolutionary stochastic search algorithm. PLoS Genet 9: e1003657. doi: 10.1371/journal.pgen.1003657 23950726
14. Malek TR, Castro I, (2010) Interleukin-2 receptor signaling: at the interface between tolerance and immunity. Immunity 33: 153–165. doi: 10.1016/j.immuni.2010.08.004 20732639
15. Vella A, Cooper JD, Lowe CE, Walker N, Nutland S, et al. (2005) Localization of a type 1 diabetes locus in the IL2RA/CD25feng region by use of tag single-nucleotide polymorphisms. Am J Hum Genet 76: 773–779. doi: 10.1086/429843 15776395
16. International Multiple Sclerosis Genetics Consortium (IMSGC) (2007) Risk alleles for multiple sclerosis identified by a genomewide study. N Engl J Med 357: 851–862. doi: 10.1056/NEJMoa073493 17660530
17. Barton A, Thomson W, Ke X, Eyre S, Hinks A, et al. (2008) Rheumatoid arthritis susceptibility loci at chromosomes 10p15, 12q13 and 22q13. Nat Genet 40: 1156–1159. doi: 10.1038/ng.218 18794857
18. Lowe CE, Cooper JD, Brusko T, Walker NM, Smyth DJ, et al. (2007) Large-scale genetic fine mapping and genotype-phenotype associations implicate polymorphism in the IL2RA region in type 1 diabetes. Nat Genet 39: 1074–1082. doi: 10.1038/ng2102 17676041
19. Dendrou CA, Plagnol V, Fung E, Yang JHM, Downes K, et al. (2009) Cell-specific protein phenotypes for the autoimmune locus IL2RA using a genotype-selectable human bioresource. Nat Genet 41: 1011–1015. doi: 10.1038/ng.434 19701192
20. Orrù V, Steri M, Sole G, Sidore C, Virdis F, et al. (2013) Genetic variants regulating immune cell levels in health and disease. Cell 155: 242–256. doi: 10.1016/j.cell.2013.08.041 24074872
21. Ye CJ, Feng T, Kwon HK, Raj T, Wilson MT, et al. (2014) Intersection of population variation and autoimmunity genetics in human T cell activation. Science 345: 1254665. doi: 10.1126/science.1254665 25214635
22. Garg G, Tyler JR, Yang JHM, Cutler AJ, Downes K, et al. (2012) Type 1 diabetes-associated IL2RA variation lowers IL-2 signaling and contributes to diminished CD4+ CD25+ regulatory T cell function. J Immunol 188: 4644–4653. doi: 10.4049/jimmunol.1100272 22461703
23. Onengut-Gumuscu S, Chen WM, Burren O, Cooper NJ, Quinlan AR, et al. (in press) Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nat Genet.
24. Kass R, Raftery A, (1995) Bayes factors. Journal of the American Statistical Association 90: 773–795. doi: 10.1080/01621459.1995.10476572
25. Mousavi K, Zare H, Dell’orso S, Grontved L, Gutierrez-Cruz G, et al. (2013) eRNAs promote transcription by establishing chromatin accessibility at defined genomic loci. Mol Cell 51: 606–617. doi: 10.1016/j.molcel.2013.07.022 23993744
26. Maier LM, Lowe CE, Cooper J, Downes K, Anderson DE, et al. (2009) IL2RA genetic heterogeneity in multiple sclerosis and type 1 diabetes susceptibility and soluble interleukin-2 receptor production. PLoS Genet 5: e1000322. doi: 10.1371/journal.pgen.1000322 19119414
27. Cerosaletti K, Schneider A, Schwedhelm K, Frank I, Tatum M, et al. (2013) Multiple autoimmune-associated variants confer decreased IL-2R signaling in CD4+ CD25(hi) T cells of type 1 diabetic and multiple sclerosis patients. PLoS One 8: e83811. doi: 10.1371/journal.pone.0083811 24376757
28. Dendrou CA, Wicker LS, (2008) The IL-2/CD25 pathway determines susceptibility to T1D in humans and NOD mice. J Clin Immunol 28: 685–696. doi: 10.1007/s10875-008-9237-9 18780166
29. Yang Z, Fujii H, Mohan SV, Goronzy JJ, Weyand CM, (2013) Phosphofructokinase deficiency impairs ATP generation, autophagy, and redox balance in rheumatoid arthritis T cells. J Exp Med 210: 2119–2134. doi: 10.1084/jem.20130252 24043759
30. Arden C, Hampson LJ, Huang GC, Shaw JAM, Aldibbiat A, et al. (2008) A role for PFK-2/FBPase-2, as distinct from fructose 2,6-bisphosphate, in regulation of insulin secretion in pancreatic beta-cells. Biochem J 411: 41–51. doi: 10.1042/BJ20070962 18039179
31. Davison LJ, Wallace C, Cooper JD, Cope NF, Wilson NK, et al. (2012) Long-range dna looping and gene expression analyses identify DEXI as an autoimmune disease candidate gene. Hum Mol Genet 21: 322–333. doi: 10.1093/hmg/ddr468 21989056
32. Carbonetto P, Stephens M, (2013) Integrated enrichment analysis of variants and pathways in genome-wide association studies indicates central role for IL-2 signaling genes in type 1 diabetes, and cytokine signaling genes in Crohn’s disease. PLoS Genet 9: e1003770. doi: 10.1371/journal.pgen.1003770 24098138
33. Hoggart CJ, Whittaker JC, De Iorio M, Balding DJ, (2008) Simultaneous analysis of all SNPs in genome-wide and re-sequencing association studies. PLoS Genet 4: e1000130. doi: 10.1371/journal.pgen.1000130 18654633
34. Wu TT, Chen YF, Hastie T, Sobel E, Lange K, (2009) Genome-wide association analysis by lasso penalized logistic regression. Bioinformatics 25: 714–721. doi: 10.1093/bioinformatics/btp041 19176549
35. Ayers KL, Cordell HJ, (2010) SNP selection in genome-wide and candidate gene studies via penalized logistic regression. Genet Epidemiol 34: 879–891. doi: 10.1002/gepi.20543 21104890
36. Guan Y, Stephens M, (2011) Bayesian variable selection regression for genome-wide association studies and other large-scale problems. Ann Appl Stat 5: 1780–1815. doi: 10.1214/11-AOAS455
37. Leng C, Lin Y, Wahba G, (2006) A note on the lasso and related procedures in model selection. Stat Sin.
38. Zhao P, Yu B, (2006) On model selection consistency of lasso. J Mach Learn Res 7: 2541–2563.
39. Zou H, Hastie T, (2005) Regularization and variable selection via the elastic net. J R Statist Soc B 67: 301–320. doi: 10.1111/j.1467-9868.2005.00503.x
40. Yuan M, Lin Y, (2006) Model selection and estimation in regression with grouped variables. J R Stat Soc Series B Stat Methodol 68: 49–67. doi: 10.1111/j.1467-9868.2005.00532.x
41. Xiong HY, Alipanahi B, Lee LJ, Bretschneider H, Merico D, et al. (2015) RNA splicing. the human splicing code reveals new insights into the genetic determinants of disease. Science 347: 1254806.
42. Chen H, Wilkins LM, Aziz N, Cannings C, Wyllie DH, et al. (2006) Single nucleotide polymorphisms in the human interleukin-1b gene affect transcription according to haplotype context. Hum Mol Genet 15: 519–529. doi: 10.1093/hmg/ddi469 16399797
43. Lewontin RC, (1964) The interaction of selection and linkage. i. general considerations; heterotic models. Genetics 49: 49–67. 17248194
44. Chapman JM, Cooper JD, Todd JA, Clayton DG, (2003) Detecting disease associations due to linkage disequilibrium using haplotype tags: A class of tests and the determinants of statistical power. Hum Hered 56: 18–31. doi: 10.1159/000073729 14614235
45. Morris AP, (2006) A flexible bayesian framework for modeling haplotype association with disease, allowing for dominance effects of the underlying causative variants. Am J Hum Genet 79: 679–694. doi: 10.1086/508264 16960804
46. Hemani G, Shakhbazov K, Westra HJ, Esko T, Henders AK, et al. (2014) Detection and replication of epistasis influencing transcription in humans. Nature 508: 249–253. doi: 10.1038/nature13005 24572353
47. Wood AR, Tuke MA, Nalls MA, Hernandez DG, Bandinelli S, et al. (2014) Another explanation for apparent epistasis. Nature 514: E3–5. doi: 10.1038/nature13691 25279928
48. Downes K, Pekalski M, Angus KL, Hardy M, Nutland S, et al. (2010) Reduced expression of IFIH1 is protective for type 1 diabetes. PLoS ONE 5: e12646. doi: 10.1371/journal.pone.0012646 20844740
49. Bottolo L, Petretto E, Blankenberg S, Cambien F, Cook SA, et al. (2011) Bayesian detection of expression quantitative trait loci hot spots. Genetics 189: 1449–1459. doi: 10.1534/genetics.111.131425 21926303
50. Okada Y, Wu D, Trynka G, Raj T, Terao C, et al. (2014) Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature 506: 376–381. doi: 10.1038/nature12873 24390342
51. Howie BN, Donnelly P, Marchini J, (2009) A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet 5: e1000529. doi: 10.1371/journal.pgen.1000529 19543373
52. Liquet B, Bottolo L, Campanella G, Richardson S, Chadeau-Hyam M, (in press) R2GUESS: GPU-based R package for Bayesian variable selection regression of multivariate responses. J Stat Softw.
53. Tibshirani R, (1996) Optimal reinsertion:regression shrinkage and selection via the lasso. J R Statist Soc B 58: 267–288.
54. Friedman JH, Hastie T, Tibshirani R, (2010) Regularization paths for generalized linear models via coordinate descent. Journal of Statistical Software 33: 1–22. 20808728
55. Yang Y, Zou H, (2014) A fast unified algorithm for solving group-lasso penalize learning problems. Stat Comput: 1–13. doi: 10.1080/00949655.2014.975226
56. Pekalski ML, Ferreira RC, Coulson RMR, Cutler AJ, Guo H, et al. (2013) Postthymic expansion in human CD4 naive T cells defined by expression of functional high-affinity IL-2 receptors. J Immunol 190: 2554–2566. doi: 10.4049/jimmunol.1202914 23418630
57. Romanoski CE, Glass CK, Stunnenberg HG, Wilson L, Almouzni G, (2015) Epigenomics: Roadmap for regulation. Nature 518: 314–316. doi: 10.1038/518314a 25693562
58. Anders S, Pyl PT, Huber W (2014) HTSeq—a Python framework to work with high-throughput sequencing data. bioRxiv: 10.1101/002824.
59. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, et al. (2013) STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29: 15–21. doi: 10.1093/bioinformatics/bts635 23104886
60. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, et al. (2009) The sequence alignment/map format and SAMtools. Bioinformatics 25: 2078–2079. doi: 10.1093/bioinformatics/btp352 19505943
61. Love MI, Huber W, Anders S (2014) Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. bioRxiv: 10.1101/002832.