Skip to main content

Scaling probabilistic models of genetic variation to millions of humans

Author(s): Gopalan, Prem; Hao, Wei; Blei, David M; Storey, John D

To refer to this page use:
Abstract: A major goal of population genetics is to quantitatively understand variation of genetic polymorphisms among individuals. The aggregated number of genotyped humans is currently on the order millions of individuals, and existing methods do not scale to data of this size. To solve this problem we developed TeraStructure, an algorithm to fit Bayesian models of genetic variation in structured human populations on tera-sample-sized data sets (1012 observed genotypes, e.g., 1M individuals at 1M SNPs). TeraStructure is a scalable approach to Bayesian inference in which subsamples of markers are used to update an estimate of the latent population structure between samples. We demonstrate that TeraStructure performs as well as existing methods on current globally sampled data, and we show using simulations that TeraStructure continues to be accurate and is the only method that can scale to tera-sample-sizes.
Publication Date: Dec-2016
Electronic Publication Date: 7-Nov-2016
Citation: Gopalan, Prem, Hao, Wei, Blei, David M, Storey, John D. (2016). Scaling probabilistic models of genetic variation to millions of humans. Nature Genetics, 48 (12), 1587 - 1590. doi:10.1038/ng.3710
DOI: doi:10.1038/ng.3710
ISSN: 1061-4036
EISSN: 1546-1718
Pages: 1587 - 1590
Type of Material: Journal Article
Journal/Proceeding Title: Nature Genetics
Version: Author's manuscript

Items in OAR@Princeton are protected by copyright, with all rights reserved, unless otherwise indicated.