An efficient hierarchical generalized linear mixed model for pathway analysis of genome-wide association studies.
BACKGROUND
In genome-wide association studies (GWAS) of complex diseases, genetic variants having real but weak associations often fail to be detected at the stringent genome-wide significance level. Pathway analysis, which tests disease association with combined association signals from a group of variants in the same pathway, has become increasingly popular. However, because of the complexities in genetic data and the large sample sizes in typical GWAS, pathway analysis remains to be challenging. We propose a new statistical model for pathway analysis of GWAS. This model includes a fixed effects component that models mean disease association for a group of genes, and a random effects component that models how each gene's association with disease varies about the gene group mean, thus belongs to the class of mixed effects models.


RESULTS
The proposed model is computationally efficient and uses only summary statistics. In addition, it corrects for the presence of overlapping genes and linkage disequilibrium (LD). Via simulated and real GWAS data, we showed our model improved power over currently available pathway analysis methods while preserving type I error rate. Furthermore, using the WTCCC Type 1 Diabetes (T1D) dataset, we demonstrated mixed model analysis identified meaningful biological processes that agreed well with previous reports on T1D. Therefore, the proposed methodology provides an efficient statistical modeling framework for systems analysis of GWAS.


BACKGROUND
The software code for mixed models analysis is freely available at http://biostat.mc.vanderbilt.edu/LilyWang.
Department of Biostatistics, Vanderbilt University, Nashville, TN 37232, USA. lily.wang@vanderbilt.edu
Lily
L
Wang
Peilin
P
Jia
Russell D
RD
Wolfinger
Xi
X
Chen
Britney L
BL
Grayson
Thomas M
TM
Aune
Zhongming
Z
Zhao
D003198:Computer Simulation; D055106:Genome-Wide Association Study / Q000379:methods*; D005838:Genotype; D006801:Humans; D016014:Linear Models; D015810:Linkage Disequilibrium; D020641:Polymorphism, Single Nucleotide; D012984:Software
Nat Genet. 2009 Jun;41(6):703-7
19430480
Genet Epidemiol. 2010 Nov;34(7):716-24
20842628
PLoS Genet. 2008 Jul;4(7):e1000115
18852846
PLoS One. 2009 Nov 30;4(11):e8068
19956648
Nature. 2007 Jun 7;447(7145):661-78
17554300
Bioinformatics. 2008 Jan 1;24(1):146-8
18024970
Diabetes. 1986 Oct;35(10):1119-23
3530842
Genet Epidemiol. 2010 Apr;34(3):222-231
20013942
Am J Hum Genet. 2010 Jan;86(1):6-22
20074509
Schizophr Res. 2010 Sep;122(1-3):38-42
20659789
Proc Natl Acad Sci U S A. 2009 Jun 9;106(23):9362-7
19474294
Am J Hum Genet. 2007 Aug;81(2):397-404
17668389
Nat Genet. 2000 May;25(1):25-9
10802651
Lancet. 1974 Oct 12;2(7885):864-6
4137711
Genomics. 2008 Nov;92(5):265-72
18722519
Am J Hum Genet. 2007 Dec;81(6):1278-83
17966091
Am J Hum Genet. 2010 Feb 12;86(2):113-25
20060087
Genet Epidemiol. 2008 Nov;32(7):658-68
18481796
Am J Hum Genet. 2010 Jun 11;86(6):860-71
20560206
Math Biosci. 2002 Mar;176(1):35-51
11867082
J Comput Biol. 2001;8(6):625-37
11747616
Nat Genet. 2003 Jul;34(3):267-73
12808457
PLoS Genet. 2008 Aug 29;4(8):e1000167
18769715
Proc Natl Acad Sci U S A. 2005 Oct 25;102(43):15545-50
16199517
Bioinformatics. 2009 Jan 15;25(2):237-42
19029127
Nucleic Acids Res. 1999 Jan 1;27(1):29-34
9847135
Lancet. 1978 Apr 29;1(8070):908-10
76848
Am J Hum Genet. 2009 Jul;85(1):13-24
19539887
Nat Rev Genet. 2008 May;9(5):356-69
18398418
Bioinformatics. 2007 Oct 1;23(19):2581-8
17785348
Nat Genet. 2010 Apr;42(4):355-60
20208535
J Clin Invest. 2008 May;118(5):1590-605
18451988
Am J Hum Genet. 2010 Jun 11;86(6):929-42
20560208
Genet Epidemiol. 2009 Dec;33(8):700-9
19333968
Diabetes. 1984 Jun;33(6):552-5
6609855
J Am Stat Assoc. 2010 Sep 1;105(491):1042-1055
21052523
Am J Hum Genet. 2003 Feb;72(2):351-63
12525994
Nat Genet. 2007 Sep;39(9):1045-51
17728769
Nat Genet. 2006 Feb;38(2):203-8
16380716