Analysis of human mini-exome sequencing data from Genetic Analysis Workshop 17 using a Bayesian hierarchical mixture model

详细信息    查看全文

• 作者：Julio S Bueno Filho (1) (2)
  Gota Morota (1)
  Quoc Tran (3)
  Matthew J Maenner (4)
  Lina M Vera-Cala (4) (5)
  Corinne D Engelman (4)
  Kristin J Meyers (4)
  
• 刊名：BMC Proceedings
• 出版年：2011
• 出版时间：December 2011
• 年：2011
• 卷：5
• 期：9-supp
• 全文大小：221KB
• 参考文献：1. Manolio TA: Genomewide association studies and assessment of the risk of disease. / New Engl J Med 2010, 363:166鈥?76. CrossRef
  2. Gorlov IP, Gorlova OY, Sunyaev SR, Spitz MR, Amos CI: Shifting paradigm of association studies: Value of rare single-nucleotide polymorphisms. / Am J Hum Genet 2008, 82:100鈥?12. j.ajhg.2007.09.006">CrossRef
  3. Ng SB, Turner EH, Robertson PD, Flygare SD, Bigham AW, Lee C, Shaffer T, Wong M, Bhattacharjee A, Eichler EE, / et al.: Targeted capture and massively parallel sequencing of 12 human exomes. / Nature 2009, 461:272鈥?76. CrossRef
  4. Dering C, Pugh E, Ziegler A: Statistical analysis of rare sequence variants: an overview of collapsing methods. / Genet Epidemiol 2011,X(suppl X):X-X.
  5. Meuwissen TH, Hayes BJ, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. / Genetics 2001, 157:1819鈥?829.
  6. Yi N, George V, Allison DB: Stochastic search variable selection for identifying multiple quantitative trait loci. / Genetics 2003, 164:1129鈥?138.
  7. Meuwissen THE, Solberg TR, Shepherd R, Woolliams JA: A fast algorithm for BayesB type of prediction of genome-wide estimates of genetic value. / Genet Select Evol 2009, 41:2. CrossRef
  8. Raftery AE, Lewis SM: One long run with diagnostics: implementation strategies for Markov chain Monte Carlo. / Stat Sci 1992, 7:493鈥?97. CrossRef
  9. Raftery AE, Lewis SM: The number of iterations, convergence diagnostics, and generic Metropolis algorithms. In / Practical Markov Chain Monte Carlo. Edited by: WR Gilks, DJ Spiegelhalter, S Richardson. London, Chapman & Hall; 1995.
  10. R Development Core Team: R: a language and environment for statistical computing. [ject.org" class="a-plus-plus">http://www.R-project.org] Vienna, Austria, R Foundation for Statistical Computing; 2010.
  11. Meuwissen T, Goddard M: Accurate prediction of genetic values for complex traits by whole-genome resequencing. / Genetics 2010, 185:623鈥?31. CrossRef
  12. George EI, McCulloch RE: Approaches for Bayesian variable selection. / Stat Sinica 1997, 7:339鈥?73.
  13. Luedtke A, Powers S, Petersen A, Sitarik A, Bekmetjev A, Tintle NL: Evaluating methods for the analysis of rare variants in sequence data. / BMC Proc 2011,5(suppl 9):S119.
  
• 作者单位：Julio S Bueno Filho (1) (2)
  Gota Morota (1)
  Quoc Tran (3)
  Matthew J Maenner (4)
  Lina M Vera-Cala (4) (5)
  Corinne D Engelman (4)
  Kristin J Meyers (4)
  
  1. Department of Dairy Science, University of Wisconsin鈥揗adison, 444 Animal Science Building, 1675 Observatory Drive, Madison, WI, 53706-1284, USA
  2. Departamento de Ci锚ncias Exatas, Universidade Federal de Lavras, PO Box 3037, Lavras, MG, 37200-000, Brazil
  3. Department of Statistics, University of Wisconsin鈥揗adison, 1300 University Avenue, Madison, WI, 53706, USA
  4. Department of Population Health Sciences, University of Wisconsin鈥揗adison, 707 WARF Building, 610 North Walnut Street, Madison, WI, 53726, USA
  5. Departamento de Salud P煤blica Universidad Industrial de Santander, Carrera 32 #29-31 Piso 3, Bucaramanga, Santander, 680002, Colombia
  

文摘

Next-generation sequencing technologies are rapidly changing the field of genetic epidemiology and enabling exploration of the full allele frequency spectrum underlying complex diseases. Although sequencing technologies have shifted our focus toward rare genetic variants, statistical methods traditionally used in genetic association studies are inadequate for estimating effects of low minor allele frequency variants. Four our study we use the Genetic Analysis Workshop 17 data from 697 unrelated individuals (genotypes for 24,487 autosomal variants from 3,205 genes). We apply a Bayesian hierarchical mixture model to identify genes associated with a simulated binary phenotype using a transformed genotype design matrix weighted by allele frequencies. A Metropolis Hasting algorithm is used to jointly sample each indicator variable and additive genetic effect pair from its conditional posterior distribution, and remaining parameters are sampled by Gibbs sampling. This method identified 58 genes with a posterior probability greater than 0.8 for being associated with the phenotype. One of these 58 genes, PIK3C2B was correctly identified as being associated with affected status based on the simulation process. This project demonstrates the utility of Bayesian hierarchical mixture models using a transformed genotype matrix to detect genes containing rare and common variants associated with a binary phenotype.